Unveiling the quantum secrets of triel metal triangles: a tale of stability, aromaticity, and relativistic effects.

Low lying electronic states of Al3-, Ga3-, In3-, and Tl3- have been characterized using high level multiconfigurational quasi degenerate perturbation theory on the multiconfigurational self-consistent field. Among these species, the singlet states emerge as the predominant energy minima, displaying remarkable stability. However, within the Tl3- series, our investigation leads to the identification of the high-spin , as the most stable spin state, a result corroborated by previous experimental detection via photoelectron spectroscopy. Similarly, we have also identified the singlet state of In3- as the signal detected previously experimentally. By applying Mandado's rules and an array of aromaticity indicators, it is conclusively demonstrated that both the singlet and quintet states exhibit multiple-fold aromaticity, while the triplets exhibit conflicting aromaticity. Furthermore, this investigation highlights the significant impact of relativistic effects, as they enhance the stability of the state relative to its singlet counterpart. These findings shed new light on the electronic structures and properties of these ions, offering valuable insights into their chemical behavior and potential applications.

Correction: The stability of biradicaloid versus closed-shell [E(µ-XR)]2 (E = P, As; X = N, P, As) rings. Does aromaticity play a role?

Correction for 'The stability of biradicaloid versus closed-shell [E(µ-XR)]2 (E = P, As; X = N, P, As) rings. Does aromaticity play a role?' by Rafael Grande-Aztatzi et al., Phys. Chem. Chem. Phys., 2016, 18, 11879-11884.

Correction: The aromaticity of dicupra[10]annulenes.

Correction for 'The aromaticity of dicupra[10]annulenes' by Rafael Grande-Aztatzi et al., Phys. Chem. Chem. Phys., 2017, 19, 9669-9675.

The aromaticity of dicupra[10]annulenes.

An extensive theoretical investigation of the electronic structure of a tested fair model dicupra[10]annulene compound, based on the analysis of atom-pair delocalization indices, Bader's molecular graph, the inspection of the canonical molecular orbitals, the z components of their Nuclear Independent Chemical Shifts, NICS(0)zz, and the normalized Giambiagi multicenter delocalization indices, concludes that the perimeter aromaticity of the dicupra[10]annulene ring is consistent with both 10 and 14 π-electron Hückel aromatic 10-membered rings. In either case, the 10-membered ring encloses two 6 π-electron aromatic inner rings, hinged at the Cu-Cu bond. This work demonstrates that the aromaticity of dicupra[10]annulenes closely resembles that of naphthalene. Hence, they are best regarded as metalla-polyacenes, which could make the building blocks of extended structures such as metalated nanotubes.

Spectroscopic and Theoretical Study of Cu(I) Binding to His111 in the Human Prion Protein Fragment 106-115.

The ability of the cellular prion protein (PrP(C)) to bind copper in vivo points to a physiological role for PrP(C) in copper transport. Six copper binding sites have been identified in the nonstructured N-terminal region of human PrP(C). Among these sites, the His111 site is unique in that it contains a MKHM motif that would confer interesting Cu(I) and Cu(II) binding properties. We have evaluated Cu(I) coordination to the PrP(106-115) fragment of the human PrP protein, using NMR and X-ray absorption spectroscopies and electronic structure calculations. We find that Met109 and Met112 play an important role in anchoring this metal ion. Cu(I) coordination to His111 is pH-dependent: at pH >8, 2N1O1S species are formed with one Met ligand; in the range of pH 5-8, both methionine (Met) residues bind to Cu(I), forming a 1N1O2S species, where N is from His111 and O is from a backbone carbonyl or a water molecule; at pH <5, only the two Met residues remain coordinated. Thus, even upon drastic changes in the chemical environment, such as those occurring during endocytosis of PrP(C) (decreased pH and a reducing potential), the two Met residues in the MKHM motif enable PrP(C) to maintain the bound Cu(I) ions, consistent with a copper transport function for this protein. We also find that the physiologically relevant Cu(I)-1N1O2S species activates dioxygen via an inner-sphere mechanism, likely involving the formation of a copper(II) superoxide complex. In this process, the Met residues are partially oxidized to sulfoxide; this ability to scavenge superoxide may play a role in the proposed antioxidant properties of PrP(C). This study provides further insight into the Cu(I) coordination properties of His111 in human PrP(C) and the molecular mechanism of oxygen activation by this site.

The stability of biradicaloid versus closed-shell [E(µ-XR)]2 (E = P, As; X = N, P, As) rings. Does aromaticity play a role?

High-level multiconfigurational self-consistent field calculations, supplemented with multiconfigurational quasi-degenerate perturbation theory ab initio calculations with the aug-cc-pVTZ basis set, demonstrate that the [E(µ-XH)]2 (E = P, As; X = N, P, As) compounds possess one planar and one butterfly-like isomer. The calculations predict that for X = N, planar isomers, which bear substantial biradicaloid character, are more stable than their butterfly-like counterpart isomers, which feature closed-shell electronic structures. This has been ascribed to the fact that the increased bond angle strain at E-N-E is not compensated by the E-E σ (deformed) bond formation in the butterfly-like isomers, yielding the planar structures, which hold wider E-N-E bond angles, as the most stable isomers. As N is substituted by heavier atoms, either P or As, the E-P(As)-E bond angle strain is released and, additionally, as the formed E-E σ-bond is less deformed, the butterfly isomer becomes the most stable isomer. Subsequent evaluation of the normalized Giambiagi multicenter electron delocalization indices revealed no sign of electron delocalization in the four-membered rings and consequently, it is concluded that aromaticity does not play any role in the stabilization of the planar isomers.

Phosphorylation promotes Al(iii) binding to proteins: GEGEGSGG as a case study.

Aluminum, the third most abundant element in the Earth's crust and one of the key industrial components of our everyday life, has been associated with several neurodegenerative diseases due to its ability to promote neurofilament tangles and ß-amyloid peptide aggregation. However, the experimental characterization of aluminum speciation in vivo is a difficult task. In the present study, we develop a theoretical protocol that combines molecular dynamics simulations, clustering of structures, and density functional theory for the characterization of the binding of aluminum to the synthetic neurofilament analogue octapeptide GEGEGSGG and its phosphorylated variant. Our protocol is tested with respect to previous NMR experimental data, which allows for a full interpretation of the experimental information available and its relationship with key thermodynamic quantities. Our results demonstrate the importance of phosphorylation in the ability of a peptide to bind to aluminum. Thus, phosphorylation: (i) changes the binding pattern of aluminum to GEGEGSGG, shifting the preferential binding site from the C-terminal to S6(P); (ii) increases the binding affinity by a factor of around 15 kcal mol(-1) in free energy; and (iii) may cause significant changes in the secondary structure and stiffness of the polypeptide chain, specially in the case of bidentate binding modes. Our results shed light on the possibility of aluminum to induce aggregation of ß-amyloid proteins and neurofilament tangles.

Structural and optical properties of the naked and passivated Al5Au5 bimetallic nanoclusters.

The structural and optical properties of both the naked and passivated bimetallic Al5Au5 nanoclusters have been analyzed based on data obtained from ab initio density functional theory and quantum molecular dynamics simulations. It has been found that the Al5Au5 nanocluster possesses a hollow shaped minimum energy structure with segregated Al and Au layered domains, the former representing the electrophilic domain and the latter the nucleophilic domain. In particular, it has been shown that alkali metal cations attach in the nucleophilic domain and hop from one Au site to the next one in the picoseconds time scale, while anions are bound tightly to the Al atoms of the electrophilic domain. Simulating annealing studies are very suggestive of the proneness of the nanocluster towards coalescence into large cluster units, when the cluster is left unprotected by appropriate ligands. Further passivation studies with NaF salt suggest, nonetheless, the possibility of the isolation of the Al5Au5 cluster in molten salts or ionic liquids.

Planar pentacoordinate carbons in CBe5(4-) derivatives.

The potential energy surfaces of a series of clusters with formula CBe5Lin(n-4) (n = 1 to 5) have been systematically explored. Our computations show that the lithium cations preserve the CBe5(4-) pentagon, such that the global minimum structure for these series of clusters has a planar pentacoordinate carbon (ppC) atom. The systems are primarily connected via a network of multicenter σ-bonds, in which the C atom acts as σ-acceptor and this acceptance of charge is balanced by the donation of the 2pz electrons to the π-cloud. The induced magnetic field analysis suggests that the clusters with formula CBe5Lin(n-4) (n = 1 to 5) are fully delocalized. The fact that these ppC-containing clusters are the lowest-energy forms on the corresponding potential energy surfaces raises expectations that these species can be prepared experimentally in the gas phase.

Structural evolution of small gold clusters doped by one and two boron atoms.

The potential energy surfaces (PES) of a series of gold-boron clusters with formula Aun B (n = 1-8) and Aum B2 (m = 1-7) have been explored using a modified stochastic search algorithm. Despite the complexity of the PES of these clusters, there are well-defined growth patterns. The bonding of these clusters is analyzed using the adaptive natural density partitioning and the natural bonding orbital analyses. Reactivity is studied in terms of the molecular electrostatic potential.

Spectroscopic and electronic structure studies of copper(II) binding to His111 in the human prion protein fragment 106-115: evaluating the role of protons and methionine residues.

The prion protein (PrP(C)) is implicated in the spongiform encephalopathies in mammals, and it is known to bind Cu(II) at the N-terminal region. The region around His111 has been proposed to be key for the conversion of normal PrP(C) to its infectious isoform PrP(Sc). The principal aim of this study is to understand the role of protons and methionine residues 109 and 112 in the coordination of Cu(II) to the peptide fragment 106-115 of human PrP, using different spectroscopic techniques (UV-vis absorption, circular dichroism, and electron paramagnetic resonance) in combination with detailed electronic structure calculations. Our study has identified a proton equilibrium with a pK(a) of 7.5 associated with the Cu(II)-PrP(106-115) complex, which is ascribed to the deprotonation of the Met109 amide group, and it converts the site from a 3NO to a 4N equatorial coordination mode. These findings have important implications as they imply that the coordination environment of this Cu binding site at physiological pH is a mixture of two species. This study also establishes that Met109 and Met112 do not participate as equatorial ligands for Cu, and that Met112 is not an essential ligand, while Met109 plays a more important role as a weak axial ligand, particularly for the 3NO coordination mode. A role for Met109 as a highly conserved residue that is important to regulate the protonation state and redox activity of this Cu binding site, which in turn would be important for the aggregation and amyloidogenic properties of the protein, is proposed.

Theoretical characterization of Al(III) binding to KSPVPKSPVEEKG: Insights into the propensity of aluminum to interact with key sequences for neurofilament formation.

Classical molecular dynamic simulations and density functional theory are used to unveil the interaction of aluminum with various phosphorylated derivatives of the fragment KSPVPKSPVEEKG (NF13), a major multiphosphorylation domain of human neurofilament medium (NFM). Our calculations reveal the rich coordination chemistry of the resultant structures with a clear tendency of aluminum to form multidentate structures, acting as a bridging agent between different sidechains and altering the local secondary structure around the binding site. Our evaluation of binding energies allows us to determine that phosphorylation has an increase in the affinity of these peptides towards aluminum, although the interaction is not as strong as well-known chelators of aluminum in biological systems. Finally, the presence of hydroxides in the first solvation layer has a clear damping effect on the binding affinities. Our results help in elucidating the potential structures than can be formed between this exogenous neurotoxic metal and key sequences for the formation of neurofilament tangles, which are behind of some of the most important degenerative diseases.

Does phosphorylation increase the binding affinity of aluminum? A computational study on the aluminum interaction with serine and O-phosphoserine.

Several toxic effects arise from aluminum's presence in living systems, one of these effects is to alter the natural role of enzymes and non-enzyme proteins. Aluminum promotes the hyperphosphorylation of normal proteins. In order to assess the aluminum-binding abilities of phosphorylated proteins and peptides, the interaction of aluminum at different pH with serine and phosphoserine is studied by a Density Functional Theory study, combined with polarizable continuum models to account for bulk solvent effects, and the electronic structure of selected complexes are analyzed by Quantum Theory of "Atoms in Molecules". Our results confirm the high ability of aluminum to bind polypeptides as the pH lowers. Moreover, the phosphorylation of the building blocks increases the affinity for aluminum, in particular at physiological pH. Finally, aluminum shows a tendency to be chelated forming different size rings.

Oxidation of Acid, Base, and Amide Side-Chain Amino Acid Derivatives via Hydroxyl Radical.

Hydroxyl radical (•OH) is known to be highly reactive. Herein, we analyze the oxidation of acid (Asp and Glu), base (Arg and Lys), and amide (Asn and Gln) containing amino acid derivatives by the consecutive attack of two •OH. In this work, we study the reaction pathway by means of density functional theory. The oxidation mechanism is divided into two steps: (1) the first •OH can abstract a H atom or an electron, leading to a radical amino acid derivative, which is the intermediate of the reaction and (2) the second •OH can abstract another H atom or add itself to the formed radical, rendering the final oxidized products. The studied second attack of •OH is applicable to situations where high concentration of •OH is found, e.g., in vitro. Carbonyls are the best known oxidation products for these reactions. This work includes solvent dielectric and confirmation's effects of the reaction, showing that both are negligible. Overall, the most favored intermediates of the oxidation process at the side chain correspond to the secondary radicals stabilized by hyperconjugation. Intermediates show to be more stable in those cases where the spin density of the unpaired electron is lowered. Alcohols formed at the side chains are the most favored products, followed by the double-bond-containing ones. Interestingly, Arg and Lys side-chain scission leads to the most favored carbonyl-containing oxidation products, in line with experimental results.

Inducing magnetic communication in caged dinuclear Co(II) systems.

The synthesis, structural, electronic and magnetic characterization of five dinuclear Co(II) azacryptand compounds (1-5) bridged through different ions are reported. The magnetic exchange interactions, 2J values, obtained from theoretical computations show that the variation of the intermetallic angles and distances lead to antiferromagnetic behaviours. Magneto-structural correlations show a trend, where the angles Co(II)-bridge-Co(II) closer to 180° favour an increase in the superexchange pathway leading to higher AF interaction values.

Insights into the oxygen-based ligand of the low pH component of the Cu(2+)-amyloid-ß complex.

In spite of significant experimental effort dedicated to the study of Cu(2+) binding to the amyloid beta (Aß) peptide, involved in Alzheimer's disease, the nature of the oxygen-based ligand in the low pH component of the Cu(2+)-Aß(1-16) complex is still under debate. This study reports density-functional-theory-based calculations that explore the potential energy surface of Cu(2+) complexes including N and O ligands at the N-terminus of the Aß peptide, with a focus on evaluating the role of Asp1 carboxylate in copper coordination. Model conformers including 3, 6, and 17 amino acids have been used to systematically study several aspects of the Cu(2+)-coordination such as the Asp1 side chain conformation, local peptide backbone geometry, electrostatic and/or hydrogen bond interactions, and number and availability of Cu(2+) ligands. Our results show that the Asp1 peptide carbonyl binds to Cu(2+) only if the coordination number is less than four. In contrast, if four ligands are available, the most stable structures include the Asp1 carboxylate in equatorial position instead of the Asp1 carbonyl group. The two lowest energy Cu(2+)-Aß(1-17) models involve Asp1 COO(-), the N-terminus, and His6 and His14 as equatorial ligands, with either a carbonyl or a water molecule in the axial position. These models are in good agreement with experimental data reported for component I of the Cu(2+)-Aß(1-16) complex, including EXAFS- and X-ray-derived Cu(2+)-ligand distances, Cu(2+) EPR parameters, and (14)N and (13)C superhyperfine couplings. Our results suggest that at low pH, Cu(2+)-Aß species with Asp1 carboxylate equatorial coordination coexist with species coordinating the Asp1 carbonyl. Understanding the bonding mechanism in these species is relevant to gain a deeper insight on the molecular processes involving copper-amyloid-ß complexes, such as aggregation and redox activity.

Structural models for Cu(II) bound to the fragment 92-96 of the human prion protein.

The prion protein (PrP(C)) binds Cu(II) in its N-terminal region, and it is associated to a group of neurodegenerative diseases termed transmissible spongiform encephalopaties (TSEs). The isoform PrP(Sc), derived from the normal PrP(C), is the pathogenic agent of TSEs. Using spectroscopic techniques (UV-vis absorption, circular dichroism, and electron paramagnetic resonance) and electronic structure calculations, we obtained a structural description for the different pH-dependent binding modes of Cu(II) to the PrP(92-96) fragment. We have also evaluated the possibility of water molecule ligation to the His96-bound copper ion. Geometry-optimized structural models that reproduce the spectroscopic features of these complexes are presented. Two Cu(II) binding modes are relevant at physiological pH: 4N and 3NO equatorial coordination modes; these are best described by models with no participation of water molecules in the coordination sphere of the metal ion. In contrast, the 2N2O and N3O coordination modes that are formed at lower pH involve the coordination of an axial water molecule. This study underscores the importance of including explicit water molecules when modeling copper binding sites in PrP(C).