Chalcogen bonds formed by protein sulfur atoms in proteins. A survey of high-resolution structures deposited in the protein data bank

Anions as Lewis Acids in Noncovalent Bonds

The ability of an anion to serve as electron-accepting Lewis acid in a noncovalent bond is assessed via DFT calculations. NH3 is taken as the common base, and is paired with a host of ACln - anions, with central atom A=Ca, Sr, Mg, Te, Sb, Hg, Zn, Ag, Ga, Ti, Sn, I, and B. Each anion reacts through its σ or π-hole although the electrostatic potential of this hole is quite negative in most cases. Despite the contact between this negative hole and the negative region of the approaching nucleophile, the electrostatic component of the interaction energy of each bond is highly favorable, and accounts for more than half of the total attractive energy. The double negative charge of dianions precludes a stable complex with NH3.

Structure and Substrate Specificity of S-Methyl Thiourocanate Hydratase

Nicotinamide adenine dinucleotide (NAD+) is a common cofactor in enzyme-catalyzed reactions that involve hydride transfers. In contrast, urocanase and urocanase-like enzymes use NAD+ for covalent electrophilic catalysis. Deciphering avenues by which this unusual catalytic strategy has diversified by evolution may point to approaches for the design of novel enzymes. In this report, we describe the S-methyl thiourocanate hydratase (S-Me-TUC) from Variovorax sp. RA8 as a novel member of this small family of NAD+-dependent hydratases. This enzyme catalyzes the 1,4-addition of water to S-methyl thiourocanate as the second step in the catabolism of S-methyl ergothioneine. The crystal structure of this enzyme in complex with the cofactor and a product analogue identifies critical sequence motifs that explain the narrow and nonoverlapping substrate scopes of S-methyl thiourocanate-, urocanate-, thiourocanate-, and Nτ-methyl urocanate-specific hydratases. The discovery of a S-methyl ergothioneine catabolic pathway also suggests that S-methylation or alkylation may be a significant activity in the biology of ergothioneine.

Computational Insight into the Nature and Strength of the π-Hole Type Chalcogen∙∙∙Chalcogen Interactions in the XO2∙∙∙CH3YCH3 Complexes (X = S, Se, Te; Y = O, S, Se, Te)

In recent years, the non-covalent interactions between chalcogen centers have aroused substantial research interest because of their potential applications in organocatalysis, materials science, drug design, biological systems, crystal engineering, and molecular recognition. However, studies on π-hole-type chalcogen∙∙∙chalcogen interactions are scarcely reported in the literature. Herein, the π-hole-type intermolecular chalcogen∙∙∙chalcogen interactions in the model complexes formed between XO2 (X = S, Se, Te) and CH3YCH3 (Y = O, S, Se, Te) were systematically studied by using quantum chemical computations. The model complexes are stabilized via one primary X∙∙∙Y chalcogen bond (ChB) and the secondary C-H∙∙∙O hydrogen bonds. The binding energies of the studied complexes are in the range of -21.6~-60.4 kJ/mol. The X∙∙∙Y distances are significantly smaller than the sum of the van der Waals radii of the corresponding two atoms. The X∙∙∙Y ChBs in all the studied complexes except for the SO2∙∙∙CH3OCH3 complex are strong in strength and display a partial covalent character revealed by conducting the quantum theory of atoms in molecules (QTAIM), a non-covalent interaction plot (NCIplot), and natural bond orbital (NBO) analyses. The symmetry-adapted perturbation theory (SAPT) analysis discloses that the X∙∙∙Y ChBs are primarily dominated by the electrostatic component.

On the Importance of Halogen and Chalcogen Bonds in the Solid State of Nucleic Acids: A Combined Crystallographic and Theoretical Perspective

In this work, intra- and intermolecular halogen and chalcogen bonds (HlgBs and ChBs, respectively) present in the solid state of nucleic acids (NAs) have been studied at the RI-MP2/def2-TZVP level of theory. To achieve this, a Protein Data Bank (PDB) survey was carried out, revealing a series of structures in which Br/I or S/Se/Te atoms belonging to nucleobases or pentose rings were involved in noncovalent interactions (NCIs) with electron-rich species. The energetics and directionality of these NCIs were rationalized through a computational study, which included the use of Molecular Electrostatic Potential (MEP) surfaces, the Quantum Theory of Atoms in Molecules (QTAIM), and Non Covalent Interaction plot (NCIplot) and Natural Bonding Orbital (NBO) techniques.

Automated identification of chalcogen bonds in AlphaFold protein structure database files: is it possible?

Protein structure prediction and structural biology have entered a new era with an artificial intelligence-based approach encoded in the AlphaFold2 and the analogous RoseTTAfold methods. More than 200 million structures have been predicted by AlphaFold2 from their primary sequences and the models as well as the approach itself have naturally been examined from different points of view by experimentalists and bioinformaticians. Here, we assessed the degree to which these computational models can provide information on subtle structural details with potential implications for diverse applications in protein engineering and chemical biology and focused the attention on chalcogen bonds formed by disulphide bridges. We found that only 43% of the chalcogen bonds observed in the experimental structures are present in the computational models, suggesting that the accuracy of the computational models is, in the majority of the cases, insufficient to allow the detection of chalcogen bonds, according to the usual stereochemical criteria. High-resolution experimentally derived structures are therefore still necessary when the structure must be investigated in depth based on fine structural aspects.

Competition Between the Two σ-Holes in the Formation of a Chalcogen Bond

A chalcogen atom Y contains two separate σ-holes when in a R₁ YR₂ molecular bonding pattern. Quantum chemical calculations consider competition between these two σ-holes to engage in a chalcogen bond (ChB) with a NH₃ base. R groups considered include F, Br, I, and tert-butyl (tBu). Also examined is the situation where the Y lies within a chalcogenazole ring, where its neighbors are C and N. Both electron-withdrawing substituents R₁ and R₂ act cooperatively to deepen the two σ-holes, but the deeper of the two holes consistently lies opposite to the more electron-withdrawing group, and is also favored to form a stronger ChB. The formation of two simultaneous ChBs in a triad requires the Y atom to act as double electron acceptor, and so anti-cooperativity weakens each bond relative to the simple dyad. This effect is such that some of the shallower σ-holes are unable to form a ChB at all when a base occupies the other site.