Theoretical and crystallographic data investigations of noncovalent S···O interactions

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

The presence of chalcogen bonds in native proteins was investigated on a non-redundant and high-resolution (≤ 1 Angstrom) set of protein crystal structures deposited in the Protein Data Bank. It was observed that about one half of the sulfur atoms of methionines and disulfide bridges from chalcogen bonds with nucleophiles (oxygen and sulfur atoms, and aromatic rings). This suggests that chalcogen bonds are a non-bonding interaction important for protein stability. Quite numerous chalcogen bonds involve water molecules. Interestingly, in the case of disulfide bridges, chalcogen bonds have a marked tendency to occur along the S-S bond extension rather than along the C-S bond extension. Additionally, it has been observed that closer residues have a higher probability of being connected by a chalcogen bonds, while the secondary structure of the two residues connected by a chalcogen bond do not correlate with its formation.Communicated by Ramaswamy H. Sarma.

Chalcogen Bonding Catalysis with Phosphonium Chalcogenide (PCH)

ConspectusThe exploration of new catalysis concepts and strategies to drive chemical reactions is of vital importance for the sustainable development of organic synthesis. Recently, chalcogen bonding catalysis has emerged as a new concept for organic synthesis and has been demonstrated to be an important synthetic tool capable of addressing elusive reactivity and selectivity issues. This Account describes our progress in the research field of chalcogen bonding catalysis, including (1) the discovery of phosphonium chalcogenide (PCH) as highly efficient chalcogen bonding catalyst; (2) the development of "chalcogen-chalcogen bonding catalysis" and "chalcogen···π bonding catalysis" modes; (3) the demonstration that chalcogen bonding catalysis with PCH can activate hydrocarbons to achieve cyclization and coupling reactions of alkenes; (4) the discovery of unusual results that chalcogen bonding catalysis with PCH can solve elusive reactivity and selectivity issues that are inaccessible by classic catalysis approaches; and (5) the elucidation of chalcogen bonding mechanisms.With PCH catalysts, we systematically studied their chalcogen bonding properties, the relationship between structure and catalysis, and their application in facilitating a diverse array of reactions. Enabled by chalcogen-chalcogen bonding catalysis, an efficient assembly reaction of three molecules of β-ketoaldehyde and one indole derivative in a single operation was realized, delivering heterocycles with a newly constructed seven-membered ring. In addition, a Se···O bonding catalysis approach achieved an efficient synthesis of calix[4]pyrroles. We developed a "dual chalcogen bonding catalysis" strategy to solve reactivity and selectivity issues in the Rauhut-Currier-type reactions and related cascade cyclizations, thus shifting conventionally covalent Lewis base catalysis to a cooperative Se···O bonding catalysis approach. This strategy enables the cyanosilylation of ketones to take place in the presence of a ppm-level amount of PCH catalyst loading. Furthermore, we established chalcogen···π bonding catalysis for catalytic transformation of alkenes. In the research field of supramolecular catalysis, the activation of hydrocarbons such as alkenes by weak interactions is a highly interesting unresolved topic. We showed that the Se···π bonding catalysis approach could efficiently activate alkenes to achieve both coupling and cyclization reactions. Chalcogen···π bonding catalysis with PCH catalysts is particularly highlighted by the capability of facilitating strong Lewis-acid inaccessible transformations, such as the controlled cross coupling of triple alkenes. Overall, this Account presents a panoramic view of our research on chalcogen bonding catalysis with PCH catalysts. The works described in this Account provide a significant platform to solve synthetic problems.

Interplay between hydrogen and chalcogen bonds in cysteine

Protein structures are stabilized by several types of chemical interactions between amino acids, which can compete with each other. This is the case of chalcogen and hydrogen bonds formed by the thiol group of cysteine, which can form three hydrogen bonds with one hydrogen acceptor and two hydrogen donors and a chalcogen bond with a nucleophile along the extension of the CS bond. A survey of the Protein Data Bank shows that hydrogen bonds are about 40-50 more common than chalcogen bonds, suggesting that they are stronger and, consequently, prevail, though not always. It is also observed that frequently a thiol group that forms a chalcogen bond is also involved, as a hydrogen donor, in a hydrogen bond.

Chalcogen Bonds: How to Characterize Them in Solution?

Chalcogen bonds (ChBs) occur between molecules containing Lewis acidic chalcogen substituents and Lewis bases. Recently, ChB emerged as a pivotal interaction in solution-based applications such as anion recognition, anion transport and catalysis. However, before moving to applications, the involvement of ChB must be established in solution. In this Concept article, we provide a brief review of the currently available experimental investigations of ChB in solution.

Various Sorts of Chalcogen Bonds Formed by an Aromatic System

The chalcogen Y atom in the aromatic ring of thiophene and its derivatives YC₄H₄ (Y = S, Se, Te) can engage in a number of different interactions with another such unit within the homodimer. Quantum calculations show that the two rings can be oriented perpendicular to one another in a T-shaped dimer in which the Y atom accepts electron density from the π-system of the other unit in a Y···π chalcogen bond (ChB). This geometry best takes advantage of attractions between the electrostatic potentials surrounding the two monomers. There are two other geometries in which the two Y atoms engage in a ChB with one another. However, instead of a simple interaction between a σ-hole on one Y and the lone pair of its neighbor, the interaction is better described as a pair of symmetrically equivalent Y···Y interactions, in which charge is transferred in both directions simultaneously, thereby effectively doubling the strength of the bond. These geometries differ from what might be expected based simply on the juxtaposition of the electrostatic potentials of the two monomers.

Participation of S and Se in hydrogen and chalcogen bonds

The heavier chalcogen atoms S, Se, and Te can each participate in a range of different noncovalent interactions. They can serve as both proton donor and acceptor in H-bonds. Each atom can also act as electron acceptor in a chalcogen bond.

What is the preferred geometry of sulfur–disulfide interactions?

Combined crystallographic and quantum chemical studies showed that in most cases, in crystal structures, interactions between sulphur atoms and disulphide bonds are bifurcated.

Theoretical Studies of IR and NMR Spectral Changes Induced by Sigma-Hole Hydrogen, Halogen, Chalcogen, Pnicogen, and Tetrel Bonds in a Model Protein Environment

Various types of σ-hole bond complexes were formed with FX, HFY, H2FZ, and H3FT (X = Cl, Br, I; Y = S, Se, Te; Z = P, As, Sb; T = Si, Ge, Sn) as Lewis acid. In order to examine their interactions with a protein, N-methylacetamide (NMA), a model of the peptide linkage was used as the base. These noncovalent bonds were compared by computational means with H-bonds formed by NMA with XH molecules (X = F, Cl, Br, I). In all cases, the A-F bond, which lies opposite the base and is responsible for the σ-hole on the A atom (A refers to the bridging atom), elongates and its stretching frequency undergoes a shift to the red with a band intensification, much as what occurs for the X-H bond in a H-bond (HB). Unlike the NMR shielding decrease seen in the bridging proton of a H-bond, the shielding of the bridging A atom is increased. The spectroscopic changes within NMA are similar for H-bonds and the other noncovalent bonds. The C=O bond of the amide is lengthened and its stretching frequency red-shifted and intensified. The amide II band shifts to higher frequency and undergoes a small band weakening. The NMR shielding of the O atom directly involved in the bond rises, whereas the C and N atoms both undergo a shielding decrease. The frequency shifts of the amide I and II bands of the base as well as the shielding changes of the three pertinent NMA atoms correlate well with the strength of the noncovalent bond.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

Improved Description of Sulfur Charge Anisotropy in OPLS Force Fields: Model Development and Parameterization

The atomic point-charge model used in most molecular mechanics force fields does not represent well the electronic anisotropy that is featured in many directional noncovalent interactions. Sulfur participates in several types of such interactions with its lone pairs and σ-holes. The current study develops a new model, via the addition of off-atom charged sites, for a variety of sulfur compounds in the OPLS-AA and OPLS/CM5 force fields to address the lack of charge anisotropy. Parameter optimization is carried out to reproduce liquid-state properties, torsional and noncovalent energetics from reliable quantum mechanical calculations, and electrostatic potentials. Significant improvements are obtained for computed free energies of hydration, reducing the mean unsigned errors from ca. 1.4 to 0.4-0.7 kcal/mol. Enhanced accuracy in directionality and energetics is also obtained for molecular complexes with sulfur-containing hydrogen and halogen bonds. Moreover, the new model reproduces the unusual conformational preferences of sulfur-containing compounds with 1,4-intramolecular chalcogen bonds. Transferability of the new force field parameters to cysteine and methionine is verified via molecular dynamic simulations of blocked dipeptides. The study demonstrates the effectiveness of using off-atom charge sites to address electronic anisotropy, and provides a parametrization methodology that can be applied to other systems.

S···O and S···N Sulfur Bonding Interactions in Protein-Ligand Complexes: Empirical Considerations and Scoring Function

Sulfur bonding interactions between organosulfur compounds and proteins were examined using crystal structures deposited to-date in the PDB. The data was analyzed as a function of sulfur-σ-hole-bonding (i.e., sulfur bonds) to main chain Lewis bases, viz. oxygen and nitrogen atoms of the backbone amide linkages. The analyses also included an examination of sulfur bonding to side chain Lewis bases (O, N, and S) and to the "non-classical" Lewis bases present in electron-rich aromatic amino acids as-well-as to donor-acceptor bond angle distributions. The interactions analyzed included those restricted to the sum of van der Waals radii of the respective atoms or to a distance of 4 Å. The surveyed data revealed that sulfur bonding tendencies (C-S-C bond angles) were impacted not only by steric effects but perhaps also by enthalpic features present in both the donor and acceptor participants. This knowledge is not only of fundamental interest but is also important in terms of materials and drug-design involving moieties incorporating the sulfur atom. Additionally, a new empirical scoring function was developed to address the anisotropy of sulfur in protein-ligand interactions. This newly developed scoring function is incorporated into AutoDock Vina molecular docking program and is valuable for modeling and drug design.

Highly Dynamic Anion-Quadrupole Networks in Proteins

The dynamics of anion-quadrupole (or anion-π) interactions formed between negatively charged (Asp/Glu) and aromatic (Phe) side chains are for the first time computationally characterized in RmlC (Protein Data Bank entry 1EP0 ), a homodimeric epimerase. Empirical force field-based molecular dynamics simulations predict anion-quadrupole pairs and triplets (anion-anion-π and anion-π-π) are formed by the protein during the simulated trajectory, which suggests that the anion-quadrupole interactions may provide a significant contribution to the overall stability of the protein, with an average of -1.6 kcal/mol per pair. Some anion-π interactions are predicted to form during the trajectory, extending the number of anion-quadrupole interactions beyond those predicted from crystal structure analysis. At the same time, some anion-π pairs observed in the crystal structure exhibit marginal stability. Overall, most anion-π interactions alternate between an "on" state, with significantly stabilizing energies, and an "off" state, with marginal or null stabilizing energies. The way proteins possibly compensate for transient loss of anion-quadrupole interactions is characterized in the RmlC aspartate 84-phenylalanine 112 anion-quadrupole pair observed in the crystal structure. A double-mutant cycle analysis of the thermal stability suggests a possible loss of anion-π interactions compensated by variations of hydration of the residues and formation of compensating electrostatic interactions. These results suggest that near-planar anion-quadrupole pairs can exist, sometimes transiently, which may play a role in maintaining the structural stability and function of the protein, in an otherwise very dynamic interplay of a nonbonded interaction network as well as solvent effects.

The Halogen Bond

The halogen bond occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a halogen atom in a molecular entity and a nucleophilic region in another, or the same, molecular entity. In this fairly extensive review, after a brief history of the interaction, we will provide the reader with a snapshot of where the research on the halogen bond is now, and, perhaps, where it is going. The specific advantages brought up by a design based on the use of the halogen bond will be demonstrated in quite different fields spanning from material sciences to biomolecular recognition and drug design.

Dual functions of Lewis acid and base of Se in F2C=Se and their interplay in F 2CSe•••NH 3•••HX

High-level quantum chemical calculations of the ternary systems F2CSe∙∙∙NH3∙∙∙HX (X=BeH, BH2, OH, CN, OCH3, Cl, and F) and the corresponding binary systems have been carried out in view of geometries, vibrational frequencies, interaction energies, orbital interactions, and electron densities. The molecular electrostatic potentials of F2CSe demonstrate that the Se atom could play a dual role of Lewis acid and base to form a chalcogen bond with NH3 and a hydrogen bond or a covalent interaction with HX, respectively. The chalcogen bond can compete with the hydrogen bond for the complexes involving F2CSe, but the covalent interaction is far stronger than the chalcogen bond. In the ternary complexes, both types of interactions are strengthened by each other, characterized by a shorter binding distance, a larger electron density, and a stronger orbital interaction. The covalent interaction has a greater enhancing effect on the chalcogen bond than the hydrogen bond does, resulting in a prominent shortening of ~0.23 Å distance for the Se∙∙∙N distance in F2CSe∙∙∙NH3∙∙∙BH3. The enhancement of both interactions in the ternary complexes has been understood with the electrostatic potentials and orbital interactions. Graphical Abstract The dual functions of Lewis acid and base of Se in F2CSe are enhanced each other in the ternary complexes.

Hypervalent nonbonded interactions of a divalent sulfur atom. Implications in protein architecture and the functions

In organic molecules a divalent sulfur atom sometimes adopts weak coordination to a proximate heteroatom (X). Such hypervalent nonbonded S···X interactions can control the molecular structure and chemical reactivity of organic molecules, as well as their assembly and packing in the solid state. In the last decade, similar hypervalent interactions have been demonstrated by statistical database analysis to be present in protein structures. In this review, weak interactions between a divalent sulfur atom and an oxygen or nitrogen atom in proteins are highlighted with several examples. S···O interactions in proteins showed obviously different structural features from those in organic molecules (i.e., π(o) → σ(s)* versus n(o) → σ(s)* directionality). The difference was ascribed to the HOMO of the amide group, which expands in the vertical direction (π(o)) rather than in the plane (n(o)). S···X interactions in four model proteins, phospholipase A₂ (PLA₂), ribonuclease A (RNase A), insulin, and lysozyme, have also been analyzed. The results suggested that S···X interactions would be important factors that control not only the three-dimensional structure of proteins but also their functions to some extent. Thus, S···X interactions will be useful tools for protein engineering and the ligand design.