Five Hypotheses on the Origins of Temperature Dependence of 77Se NMR Shifts in Diselenides.

Notable thermal shifts in diselenides have been documented in 77Se NMR for more than 50 years, but no satisfactory explanation has been found. Here, five hypotheses are considered as possible explanations for the large temperature dependence of the 77Se chemical shifts of diaryl and dialkyl diselenides compared to monoselenides and selenols. Density functional theory calculations are provided to bolster hypotheses and better understand the effects of barrier height and dipole energies. It is proposed that the temperature dependence of diselenide 77Se NMR chemical shifts is due to rotation around the Se-Se bond and sampling of twisted conformers at higher temperatures. The molecular twisting is solvent dependent; here, DMSO-d6 and toluene-d8 were evaluated. No correlation was established between para-substituents on diaryl diselenides and the magnitude of the change in the 77Se NMR shift (Δδ) with temperature.

Antioxidant Chimeric Molecules: Are Chemical Motifs Additive? The Case of a Selenium-Based Ligand.

We set up an in silico experiment and designed a chimeric compound integrating molecular features from different efficient ROS (Reactive Oxygen Species) scavengers, with the purpose of investigating potential relationships between molecular structure and antioxidant activity. Furthermore, a selenium centre was inserted due to its known capacity to reduce hydroperoxides, acting as a molecular mimic of glutathione peroxidase; finally, since this organoselenide is a precursor of a N-heterocyclic carbene ligand, its Au(I) carbene complex was designed and examined. A validated protocol based on DFT (Density Functional Theory) was employed to investigate the radical scavenging activity of available sites on the organoselenide precursor ((SMD)-M06-2X/6-311+G(d,p)//M06-2X/6-31G(d)), as well as on the organometallic complex ((SMD)-M06-2X/SDD (Au), 6-311+G(d,p)//ZORA-BLYP-D3(BJ)/TZ2P), considering HAT (Hydrogen Atom Transfer) and RAF (Radical Adduct Formation) regarding five different radicals. The results of this case study suggest that the antioxidant potential of chemical motifs should not be considered as an additive property when designing a chimeric compound, but rather that the relevance of a molecular topology is derived from a chemical motif combined with an opportune chemical space of the molecule. Thus, the direct contributions of single functional groups which are generally thought of as antioxidants per se do not guarantee the efficient radical scavenging potential of a molecular species.

Indole-5,6-quinones display hallmark properties of eumelanin.

Melanins are ubiquitous biopolymers produced from phenols and catechols by oxidation. They provide photoprotection, pigmentation and redox activity to most life forms, and inspire synthetic materials with desirable optical, electronic and mechanical properties. The chemical structures of melanins remain elusive, however, creating uncertainty about their roles, and preventing the design of synthetic mimics with tailored properties. Indole-5,6-quinone (IQ) has been implicated as a biosynthetic intermediate and structural subunit of mammalian eumelanin pigments, but its instability has prevented its isolation and unambiguous characterization. Here we use steric shielding to stabilize IQ and show that 'blocked' derivatives exhibit eumelanin's characteristic ultrafast nonradiative decay and its ability to absorb light from the ultraviolet to the near-infrared. These new compounds are also redox-active and a source of paramagnetism, emulating eumelanin's unique electronic properties, which include persistent radicals. Blocked IQs are atomistically precise and tailorable molecules that can offer a bottom-up understanding of emergent properties in eumelanin and have the potential to advance the rational design of melanin-inspired materials.

Diphenyl Diselenide and SARS-CoV-2: in silico Exploration of the Mechanisms of Inhibition of Main Protease (Mpro) and Papain-like Protease (PLpro).

The SARS-CoV-2 pandemic has prompted global efforts to develop therapeutics. The main protease of SARS-CoV-2 (Mpro) and the papain-like protease (PLpro) are essential for viral replication and are key targets for therapeutic development. In this work, we investigate the mechanisms of SARS-CoV-2 inhibition by diphenyl diselenide (PhSe)2 which is an archetypal model of diselenides and a renowned potential therapeutic agent. The in vitro inhibitory concentration of (PhSe)2 against SARS-CoV-2 in Vero E6 cells falls in the low micromolar range. Molecular dynamics (MD) simulations and density functional theory (DFT) calculations [level of theory: SMD-B3LYP-D3(BJ)/6-311G(d,p), cc-pVTZ] are used to inspect non-covalent inhibition modes of both proteases via π-stacking and the mechanism of covalent (PhSe)2 + Mpro product formation involving the catalytic residue C145, respectively. The in vitro CC50 (24.61 µM) and EC50 (2.39 µM) data indicate that (PhSe)2 is a good inhibitor of the SARS-CoV-2 virus replication in a cell culture model. The in silico findings indicate potential mechanisms of proteases' inhibition by (PhSe)2; in particular, the results of the covalent inhibition here discussed for Mpro, whose thermodynamics is approximatively isoergonic, prompt further investigation in the design of antiviral organodiselenides.

Radical Scavenging Potential of Ginkgolides and Bilobalide: Insight from Molecular Modeling.

The reactive oxygen species (ROS) scavenging capacities of ginkgolides and bilobalide, which are the peculiar constituents of the extract of Ginkgo biloba, are investigated in silico (level of theory: (SMD)-M06-2X/6-311+G(d,p)//M06-2X/6-31G(d)). Unlike other popular antioxidant natural substances, the carbon backbones of these compounds are entirely aliphatic and exclusively single C-C bonds are present. The selectivity for alkoxyl radicals via hydrogen-atom transfer (HAT) is assessed; importantly, the scavenging of peroxyl radicals is also possible from a peculiar site, here labeled C10 both for ginkgolides and bilobalide. The energetics are described in detail, and the analysis discloses that the studied compounds are powerful scavengers, with thermodynamic and kinetic properties similar to those of Trolox and melatonin, and that, in addition, they display selectivity for peroxyl radicals. These are all chemical-reactivity features contributing to the therapeutic action of the extract of G. biloba.

In silico analysis of the antidepressant fluoxetine and similar drugs as inhibitors of the human protein acid sphingomyelinase: a related SARS-CoV-2 inhibition pathway.

Acid Sphingomyelinase (ASM) is a human phosphodiesterase that catalyzes the metabolism of sphingomyelin (SM) to ceramide and phosphocholine. ASM is involved in the plasma membrane cell repair and is associated with the lysosomal inner lipid membrane by nonbonding interactions. The disruption of those interaction would result in ASM release into the lysosomal lumen and consequent degradation of its structure. Furthermore, SARS-CoV-2 infection has been linked with ASM activation and with a ceramide domain formation in the outer leaflet of the plasma membrane that is thought to be crucial for the viral particles recognition by the host cells. In this study, we have explored in silico the behavior of fluoxetine and related drugs as potential inhibitors of ASM. Theoretically, these drugs would be able to overpass lysosomal membrane and reach the interactions that sustain ASM structure, breaking them and inhibiting the ASM. The analyses of docking data indicated that fluoxetine allocated mainly in the N-terminal saposin domain via nonbonding interactions, mostly of hydrophobic nature. Similar results were obtained for venlafaxine, citalopram, atomoxetine, nisoxetine and fluoxetine's main metabolite norfluoxetine. In conclusion, it was observed that the saposin allocation may be a good indicative of the drugs inhibition mechanism, once this domain is responsible for the binding of ASM to lysosomal membrane and some of those drugs have previously been reported to inhibit the phosphodiesterase by releasing its structure in the lysosomal lumen. Our MD data also provides some insight about natural ligand C18 sphingomyelin conformations on saposin.Communicated by Ramaswamy H. Sarma.

In silico studies of Mpro and PLpro from SARS-CoV-2 and a new class of cephalosporin drugs containing 1,2,4-thiadiazole.

The SARS-CoV-2 proteases Mpro and PLpro are important targets for the development of antivirals against COVID-19. The functional group 1,2,4-thiadiazole has been indicated to inhibit cysteinyl proteases, such as papain and cathepsins. Of note, the 1,2,4-thiadiazole moiety is found in a new class of cephalosporin FDA-approved antibiotics: ceftaroline fosamil, ceftobiprole, and ceftobiprole medocaril. Here we investigated the interaction of these new antibiotics and their main metabolites with the SARS-CoV-2 proteases by molecular docking, molecular dynamics (MD), and density functional theory (DFT) calculations. Our results indicated the PLpro enzyme as a better in silico target for the new antibacterial cephalosporins. The results with ceftaroline fosamil and the dephosphorylate metabolite compounds should be tested as potential inhibitor of PLpro, Mpro, and SARS-CoV-2 replication in vitro. In addition, the data here reported can help in the design of new potential drugs against COVID-19 by exploiting the S atom reactivity in the 1,2,4-thiadiazole moiety. Supplementary Information: The online version contains supplementary material available at 10.1007/s11224-022-02036-5.

Parameter free evaluation of SN2 reaction rates for halide substitution in halomethane.

We estimate the kinetic constants of a series of archetypal SN2 reactions, i.e., the nucleophilic substitutions of halides in halomethane. A parameter free, multiscale approach recently developed [Campeggio et al., Phys. Chem. Chem. Phys., 2020, 22, 3455] is employed. The protocol relies on quantum mechanical calculations for the description of the energy profile along the intrinsic reaction coordinate, which is then mapped onto a reaction coordinate conveniently built for the reactive process. A Kramers-Klein equation is used to describe the stochastic time evolution of the reaction coordinate and its velocity; friction is parameterized using a hydrodynamic model and Kramers theory is used to derive the rate constant of the reaction. The method is here applied to six SN2 reactions in water at 295.15 K, which differ in the nucleophile and the leaving group. The computed reaction rates are in good agreement with the experimental data and correlate well with the trends observed for the activation energies.

ROS-Scavenging Selenofluoxetine Derivatives Inhibit In Vivo Serotonin Reuptake.

While the neurochemistry that underpins the behavioral phenotypes of depression is the subject of many studies, oxidative stress caused by the inflammation comorbid with depression has not adequately been addressed. In this study, we described novel antidepressant-antioxidant agents consisting of selenium-modified fluoxetine derivatives to simultaneously target serotonin reuptake (antidepressant action) and oxidative stress. Excitingly, we show that one of these agents (1-F) carries the ability to inhibit serotonin reuptake in vivo in mice. We therefore present a frontier dual strategy that paves the way for the future of antidepressant therapies.

Drug Design: Where We Are and Future Prospects.

Medicinal chemistry is facing new challenges in approaching precision medicine. Several powerful new tools or improvements of already used tools are now available to medicinal chemists to help in the process of drug discovery, from a hit molecule to a clinically used drug. Among the new tools, the possibility of considering folding intermediates or the catalytic process of a protein as a target for discovering new hits has emerged. In addition, machine learning is a new valuable approach helping medicinal chemists to discover new hits. Other abilities, ranging from the better understanding of the time evolution of biochemical processes to the comprehension of the biological meaning of the data originated from genetic analyses, are on their way to progress further in the drug discovery field toward improved patient care. In this sense, the new approaches to the delivery of drugs targeted to the central nervous system, together with the advancements in understanding the metabolic pathways for a growing number of drugs and relating them to the genetic characteristics of patients, constitute important progress in the field.

Selenoxide Elimination Triggers Enamine Hydrolysis to Primary and Secondary Amines: A Combined Experimental and Theoretical Investigation.

We discuss a novel selenium-based reaction mechanism consisting in a selenoxide elimination-triggered enamine hydrolysis. This one-pot model reaction was studied for a set of substrates. Under oxidative conditions, we observed and characterized the formation of primary and secondary amines as elimination products of such compounds, paving the way for a novel strategy to selectively release bioactive molecules. The underlying mechanism was investigated using NMR, mass spectrometry and density functional theory (DFT).

Methylmercury Can Facilitate the Formation of Dehydroalanine in Selenoenzymes: Insight from DFT Molecular Modeling.

Experimental studies have indicated that electrophilic mercury forms (e.g., methylmercury, MeHg+) can accelerate the breakage of selenocysteine in vitro. Particularly, in 2009, Khan et al. (Environ. Toxicol. Chem. 2009, 28, 1567-1577) proposed a mechanism for the degradation of a free methylmercury selenocysteinate complex that was theoretically supported by Asaduzzaman et al. (Inorg. Chem. 2010, 50, 2366-2372). However, little is known about the fate of methylmercury selenocysteinate complexes embedded in an enzyme, especially in conditions of oxidative stress in which methylmercury target enzymes operate. Here, an accurate computational study on molecular models (level of theory: COSMO-ZORA-BLYP-D3(BJ)/TZ2P) was carried out to investigate the formation of dehydroalanine (Dha) in selenoenzymes, which irreversibly impairs their function. Methylselenocysteine as well as methylcysteine and methyltellurocysteine were included to gain insight on the peculiar behavior of selenium. Dha forms in a two-step process, i.e., the oxidation of the chalcogen nucleus followed by a syn-elimination leading to the alkene and the chalcogenic acid. The effect of an excess of hydrogen peroxide, which may lead to the formation of chalcogenones before the elimination, and of MeHg+, a severe toxicant targeting selenoproteins, which leads to the formation of methylmercury selenocysteinate, are also studied with the aim of assessing whether these pathological conditions facilitate the formation of Dha. Indeed, elimination occurs after chalcogen oxidation and MeHg+ facilitates the process. These results indicate a possible mechanism of toxicity of MeHg+ in selenoproteins.

Effect of Methylmercury Binding on the Peroxide-Reducing Potential of Cysteine and Selenocysteine.

Methylmercury (CH3Hg+) binding to catalytically fundamental cysteine and selenocysteine of peroxide-reducing enzymes has long been postulated as the origin of its toxicological activity. Only very recently, CH3Hg+ binding to the selenocysteine of thioredoxin reductase has been directly observed [Pickering, I. J. Inorg. Chem., 2020, 59, 2711-2718], but the precise influence of the toxicant on the peroxide-reducing potential of such a residue has never been investigated. In this work, we employ state-of-the-art density functional theory calculations to study the reactivity of molecular models of the free and toxified enzymes. Trends in activation energies are discussed with attention to the biological consequences and are rationalized within the chemically intuitive framework provided by the activation strain model. With respect to the free, protonated amino acids, CH3Hg+ binding promotes oxidation of the S or Se nucleus, suggesting that chalcogenoxide formation might occur in the toxified enzyme, even if the actual rate of peroxide reduction is almost certainly lowered as suggested by comparison with fully deprotonated amino acids models.

Radical scavenging activity of natural antioxidants and drugs: Development of a combined machine learning and quantum chemistry protocol.

Many natural substances and drugs are radical scavengers that prevent the oxidative damage to fundamental cell components. This process may occur via different mechanisms, among which, one of the most important, is hydrogen atom transfer. The feasibility of this process can be assessed in silico using quantum mechanics to compute ΔGHAT ○. This approach is accurate, but time consuming. The use of machine learning (ML) allows us to reduce tremendously the computational cost of the assessment of the scavenging properties of a potential antioxidant, almost without affecting the quality of the results. However, in many ML implementations, the description of the relevant features of a molecule in a machine-friendly language is still the most challenging aspect. In this work, we present a newly developed machine-readable molecular representation aimed at the application of automatized ML algorithms. In particular, we show an application on the calculation of ΔGHAT ○.

Antioxidant Potential of Psychotropic Drugs: From Clinical Evidence to In Vitro and In Vivo Assessment and toward a New Challenge for in Silico Molecular Design.

Due to high oxygen consumption, the brain is particularly vulnerable to oxidative stress, which is considered an important element in the etiopathogenesis of several mental disorders, including schizophrenia, depression and dependencies. Despite the fact that it is not established yet whether oxidative stress is a cause or a consequence of clinic manifestations, the intake of antioxidant supplements in combination with the psychotropic therapy constitutes a valuable solution in patients' treatment. Anyway, some drugs possess antioxidant capacity themselves and this aspect is discussed in this review, focusing on antipsychotics and antidepressants. In the context of a collection of clinical observations, in vitro and in vivo results are critically reported, often highlighting controversial aspects. Finally, a new challenge is discussed, i.e., the possibility of assessing in silico the antioxidant potential of these drugs, exploiting computational chemistry methodologies and machine learning. Despite the physiological environment being incredibly complex and the detection of meaningful oxidative stress biomarkers being all but an easy task, a rigorous and systematic analysis of the structural and reactivity properties of antioxidant drugs seems to be a promising route to better interpret therapeutic outcomes and provide elements for the rational design of novel drugs.

Chalcogen-mercury bond formation and disruption in model Rabenstein's reactions: A computational analysis.

Methylmercury is a highly toxic compound and human exposure is mainly related to consumption of polluted fish and seafood. The inactivation of thiol-based enzymes, promoted by the strong affinity binding of electrophilic mercuric ions to thiol and selenol groups of proteins, is likely an important factor explaining its toxicity. A key role is played by the chemistry and reactivity of the mercury-chalcogens bond, particularly HgS and HgSe, which is the focus of this computational work (level of theory: (COSMO)-ZORA-BLYP-D3(BJ)/TZ2P). We analyze nine ligand-exchange model reactions (the so-called Rabenstein's reactions) involving an entering ligand (methylchalcogenolate) and a substrate (methylchalcogenolatemethylmercury). Trends in reaction and activation energies are discussed and a change in mechanism is reported for all cases when going from gas phase to water, that is, from a single-well potential energy surface (PES) to a canonical SN 2-like mechanism. The reasons accounting for the biochemically challenging and desired displacement of methylmercury from a seleno/thiol protein can be found already in these model reactions, as can be seen from the similarities of the ligand exchange reactions in solution in thermodynamics and kinetics.

Multiscale modeling of reaction rates: application to archetypal SN2 nucleophilic substitutions.

We propose an approach to the evaluation of kinetic rates of elementary chemical reactions within Kramers' theory based on the definition of the reaction coordinate as a linear combination of natural, pseudo Z-matrix, internal coordinates of the system. The element of novelty is the possibility to evaluate the friction along the reaction coordinate, within a hydrodynamic framework developed recently [J. Campeggio et al., J. Comput. Chem. 2019, 40, 679-705]. This, in turn, allows to keep into account barrier recrossing, i.e. the transmission coefficient that is employed in correcting transition state theory evaluations. To test the capabilities and the flaws of the approach we use as case studies two archetypal SN2 reactions. First, we consider to the standard substitution of chloride ion to bromomethane. The rate constant at 295.15 K is evaluated to k/c⊖ = 2.7 × 10-6 s-1 (with c⊖ = 1 M), which compares well to the experimental value of 3.3 × 10-6 s-1 [R. H. Bathgate and E. A. Melwyn-Hughes, J. Chem. Soc 1959, 2642-2648]. Then, the method is applied to the SN2 reaction of methylthiolate to dimethyl disulfide in water. In biology, such an interconversion of thiols and disulfides is an important metabolic topic still not entirely rationalized. The predicted rate constant is k/c⊖ = 7.7 × 103 s-1. No experimental data is available for such a reaction, but it is in accord with the fact that the alkyl thiolates to dialkyl disulfides substitutions in water have been found to be fast reactions [S. M. Bachrach, J. M. Hayes, T. Dao and J. L. Mynar, Theor. Chem. Acc. 2002, 107, 266-271].

Fluoxetine scaffold to design tandem molecular antioxidants and green catalysts.

Fluoxetine finds application in the treatment of depression and mood disorders. This selective serotonin-reuptake inhibitor (SSRI) also contrasts oxidative stress by direct ROS scavenging, modulation of the endogenous antioxidant defense system, and/or enhancement of the serotonin antioxidant capacity. We synthesised some fluoxetine analogues incorporating a selenium nucleus, thus expanding its antioxidant potential by enabling a hydroperoxides-inactivating, glutathione peroxidase (GPx)-like activity. Radical scavenging and peroxidatic activity were combined in a water-soluble, drug-like, tandem antioxidant molecule. Selenofluoxetine derivatives were reacted with H2O2 in water, and the mechanistic details of the reaction were unravelled combining nuclear magnetic resonance (NMR), electrospray ionisation-mass spectrometry (ESI-MS) and quantum chemistry calculations. The observed oxidation-elimination process led to the formation of seleninic acid and cinnamylamine in a trans-selective manner. This mechanism is likely to be extended to other substrates for the preparation of unsaturated cinnamylamines.

Human elastin-like polypeptides as a versatile platform for exploitation of ultrasensitive bilirubin detection by UnaG.

A new, bifunctional recombinant protein was expressed as the fusion product of human elastin-like polypeptide (HELP) and the bilirubin-binding protein UnaG. The engineered product displays both the HELP-specific property of forming a functional hydrogel matrix and the UnaG-specific capacity of emitting green fluorescence upon ligand binding. The new fusion protein has been proven to be effective at detecting bilirubin in complex environments with high background noise. A cell culture model of the stress response, consisting of bilirubin released in the cell culture medium, was set up to assess the bilirubin-sensing properties of the functional matrix obtained by cross-linking the HELP moiety. Our engineered protein allowed us to monitor cell induction by the release of bilirubin in the culture medium on a nanomolar scale. This study shows that elastin-like protein fusion represents a versatile platform for the development of novel and commercially viable analytical and biosensing devices.

The 125Te Chemical Shift of Diphenyl Ditelluride: Chasing Conformers over a Flat Energy Surface.

The interest in diphenyl ditelluride (Ph2Te2) is related to its strict analogy to diphenyl diselenide (Ph2Se2), whose capacity to reduce organic peroxides is largely exploited in catalysis and green chemistry. Since the latter is also a promising candidate as an antioxidant drug and mimic of the ubiquitous enzyme glutathione peroxidase (GPx), the use of organotellurides in medicinal chemistry is gaining importance, despite the fact that tellurium has no recognized biological role and its toxicity must be cautiously pondered. Both Ph2Se2 and Ph2Te2 exhibit significant conformational freedom due to the softness of the inter-chalcogen and carbon⁻chalcogen bonds, preventing the existence of a unique structure in solution. Therefore, the accurate calculation of the NMR chemical shifts of these flexible molecules is not trivial. In this study, a detailed structural analysis of Ph2Te2 is carried out using a computational approach combining classical molecular dynamics and relativistic density functional theory methods. The goal is to establish how structural changes affect the electronic structure of diphenyl ditelluride, particularly the 125Te chemical shift.