On the antibacterial photodynamic inactivation mechanism of Emodin and Dermocybin natural photosensitizers: A theoretical investigation.

A DFT and TDDFT study has been carried out on monomeric anthraquinones Emodin and Dermocybin (Em, Derm) recently proposed as natural antibacterial photosensitizers able to act also against gram-negative microbes. The computational study has been performed considering the relative amount of neutral and ionic forms of each compound in water, with the variation of pH. The occurrence of both Type I and Type II photoreactions has been explored computing the absorption properties of each species, the spin-orbit coupling constants (SOC), the vertical ionization potentials and the vertical electron affinities. The most plausible deactivation channels leading to the population of excited triplet states have been proposed. Our data indicate Emodin as more active than Dermocybin in antimicrobial photodynamic therapy throughout the Type II mechanism. Our data support a dual TypeI/II activity of the monomeric anthraquinones Emodin and Dermccybin in water, in all the considered protonation states.

Viral RNA Replication Suppression of SARS-CoV-2: Atomistic Insights into Inhibition Mechanisms of RdRp Machinery by ddhCTP.

The nonstructural protein 12, known as RNA-dependent RNA polymerase (RdRp), is essential for both replication and repair of the viral genome. The RdRp of SARS-CoV-2 has been used as a promising candidate for drug development since the inception of the COVID-19 spread. In this work, we performed an in silico investigation on the insertion of the naturally modified pyrimidine nucleobase ddhCTP into the SARS-CoV-2 RdRp active site, in a comparative analysis with the natural one (CTP). The modification in ddhCTP involves the removal of the 3'-hydroxyl group that prevents the addition of subsequent nucleotides into the nascent strand, acting as an RNA chain terminator inhibitor. Quantum mechanical investigations helped to shed light on the mechanistic source of RdRp activity on the selected nucleobases, and comprehensive all-atom simulations provided insights about the structural rearrangements occurring in the active-site region when inorganic pyrophosphate (PPi) is formed. Subsequently, the intricate pathways for the release of PPi, the catalytic product of RdRp, were investigated using Umbrella Sampling simulations. The results are in line with the available experimental data and contribute to a more comprehensive point of view on such an important viral enzyme.

On the Role of Temperature in the Depolymerization of PET by FAST-PETase: An Atomistic Point of View on Possible Active Site Pre-Organization and Substrate-Destabilization Effects.

Enzyme FAST-PETase, recently obtained by a machine learning approach, can depolymerize poly(ethylene terephthalate) (PET), a synthetic resin employed in plastics and in clothing fibers. Therefore it represents a promising solution for the recycling of PET-based materials. In this study, a model of PET was adopted to describe the substrate, and all-atoms classical molecular dynamics (MD) simulations on apo- and substrate-bound FAST-PETase were carried out at 30 and 50 °C to provide atomistic details on the binding step of the catalytic cycle. Comparative analysis shed light on the interactions occurring between the FAST-PETase and 4PET at 50 °C, the optimal working conditions of the enzyme. Pre-organization of the enzyme active and binding sites has been highlighted, while MD simulations of FAST-PETase:4PET pointed out the occurrence of solvent-inaccessible conformations of the substrate promoted by the enzyme. Indeed, neither of these conformations was observed during MD simulations of the substrate alone in solution performed at 30, 50 and 150 °C. The analysis led us to propose that, at 50 °C, the FAST-PETase is pre-organized to bind the PET and that the interactions occurring in the binding site can promote a more reactive conformation of PET substrate, thus enhancing the catalytic activity of the enzyme.

Degradation of Penicillinic Antibiotics and ß-Lactamase Enzymatic Catalysis in a Biomimetic Zn-Based Metal-Organic Framework.

ß-Lactam antibiotics are one of the most commonly prescribed drugs to treat bacterial infections. However, their use has been somehow limited given the emergence of bacteria with resistance mechanisms, such as ß-lactamases, which inactivate them by degrading their four-membered ß-lactam rings. So, a total knowledge of the mechanisms governing the catalytic activity of ß-lactamases is required. Here, we report a novel Zn-based metal-organic framework (MOF, 1), possessing functional channels capable to accommodate and interact with antibiotics, which catalyze the selective hydrolysis of the penicillinic antibiotics amoxicillin and ceftriaxone. In particular, MOF 1 degrades, very efficiently, the four-membered ß-lactam ring of amoxicillin, acting as a ß-lactamase mimic, and expands the very limited number of MOFs capable to mimic catalytic enzymatic processes. Combined single-crystal X-ray diffraction (SCXRD) studies and density functional (DFT) calculations offer unique snapshots on the host-guest interactions established between amoxicillin and the functional channels of 1. This allows to propose a degradation mechanism based on the activation of a water molecule, promoted by a Zn-bridging hydroxyl group, concertedly to the nucleophilic attack to the carbonyl moiety and the cleaving of C-N bond of the lactam ring.

Degradation of Penicillinic Antibiotics and ß-Lactamase Enzymatic Catalysis in a Biomimetic Zn-Based Metal-Organic Framework.

Invited for the cover of this issue are Jesús Ferrando-Soria, Donatella Armentano, Antonio Leyva-Pérez, Emilio Pardo and co-workers at University of Valencia, Technical University of Valencia and University of Calabria. The image depicts the crystal structure of a novel ZnII biological metal-organic framework that mimics ß-lactamase enzymes. Read the full text of the article at 10.1002/chem.202301325.

Periodicity of the Affinity of Lanmodulin for Trivalent Lanthanides and Actinides: Structural and Electronic Insights from Quantum Chemical Calculations.

Lanmodulin (LanM) is the first identified macrochelator that has naturally evolved to sequester ions of rare earth elements (REEs) such as Y and all lanthanides, reversibly. This natural protein showed a 106 times better affinity for lanthanide cations than for Ca, which is a naturally abundant and biologically relevant element. Recent experiments have shown that its metal ion binding activity can be further extended to some actinides, like Np, Pu, and Am. For this reason, it was thought that LanM could be adopted for the separation of REE ions and actinides, thus increasing the interest in its potential use for industry-oriented applications. In this work, a systematic study of the affinity of LanM for lanthanides and actinides has been carried out, taking into account all trivalent ions belonging to the 4f (from La to Lu) and 5f (from Ac to Lr) series, starting from their chemistry in solution. On the basis of a recently published nuclear magnetic resonance structure, a model of the LanM-binding site was built and a detailed structural and electronic description of initial aquo- and LanM-metal ion complexes was provided. The obtained binding energies are in agreement with the available experimental data. A possible reason that could explain the origin of the affinity of LanM for these metal ions is also discussed.

Incorporation of N7-Platinated Guanines into Thermus Aquaticus (Taq) DNA Polymerase: Atomistic Insights from Molecular Dynamics Simulations.

In this work, we elucidated some key aspects of the mechanism of action of the cisplatin anticancer drug, cis-[Pt(NH3)2Cl2], involving direct interactions with free nucleotides. A comprehensive in silico molecular modeling analysis was conducted to compare the interactions of Thermus aquaticus (Taq) DNA polymerase with three distinct N7-platinated deoxyguanosine triphosphates: [Pt(dien)(N7-dGTP)] (1), cis-[Pt(NH3)2Cl(N7-dGTP)] (2), and cis-[Pt(NH3)2(H2O)(N7-dGTP)] (3) {dien = diethylenetriamine; dGTP = 5'-(2'-deoxy)-guanosine-triphosphate}, using canonical dGTP as a reference, in the presence of DNA. The goal was to elucidate the binding site interactions between Taq DNA polymerase and the tested nucleotide derivatives, providing valuable atomistic insights. Unbiased molecular dynamics simulations (200 ns for each complex) with explicit water molecules were performed on the four ternary complexes, yielding significant findings that contribute to a better understanding of experimental results. The molecular modeling highlighted the crucial role of a specific α-helix (O-helix) within the fingers subdomain, which facilitates the proper geometry for functional contacts between the incoming nucleotide and the DNA template needed for incorporation into the polymerase. The analysis revealed that complex 1 exhibits a much lower affinity for Taq DNA polymerase than complexes 2-3. The affinities of cisplatin metabolites 2-3 for Taq DNA polymerase were found to be quite similar to those of natural dGTP, resulting in a lower incorporation rate for complex 1 compared to complexes 2-3. These findings could have significant implications for the cisplatin mechanism of action, as the high intracellular availability of free nucleobases might promote the competitive incorporation of platinated nucleotides over direct cisplatin attachment to DNA. The study's insights into the incorporation of platinated nucleotides into the Taq DNA polymerase active site suggest that the role of platinated nucleotides in the cisplatin mechanism of action may have been previously underestimated.

Sulphur- and Selenium-for-Oxygen Replacement as a Strategy to Obtain Dual Type I/Type II Photosensitizers for Photodynamic Therapy.

The effect on the photophysical properties of sulfur- and selenium-for-oxygen replacement in the skeleton of the oxo-4-dimethylaminonaphthalimide molecule (DMNP) has been explored at the density functional (DFT) level of theory. Structural parameters, excitation energies, singlet-triplet energy gaps (ΔES-T), and spin-orbit coupling constants (SOC) have been computed. The determined SOCs indicate an enhanced probability of intersystem crossing (ISC) in both the thio- and seleno-derivatives (SDMNP and SeDMNP, respectively) and, consequently, an enhancement of the singlet oxygen quantum yields. Inspection of Type I reactions reveals that the electron transfer mechanisms leading to the generation of superoxide is feasible for all the compounds, suggesting a dual Type I/Type II activity.

On the Recognition of Natural Substrate CTP and Endogenous Inhibitor ddhCTP of SARS-CoV-2 RNA-Dependent RNA Polymerase: A Molecular Dynamics Study.

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 outbreak that is affecting the entire planet. As the pandemic is still spreading worldwide, with multiple mutations of the virus, it is of interest and of help to employ computational methods for identifying potential inhibitors of the enzymes responsible for viral replication. Attractive antiviral nucleotide analogue RNA-dependent RNA polymerase (RdRp) chain terminator inhibitors are investigated with this purpose. This study, based on molecular dynamics (MD) simulations, addresses the important aspects of the incorporation of an endogenously synthesized nucleoside triphosphate, ddhCTP, in comparison with the natural nucleobase cytidine triphosphate (CTP) in RdRp. The ddhCTP species is the product of the viperin antiviral protein as part of the innate immune response. The absence of the ribose 3'-OH in ddhCTP could have important implications in its inhibitory mechanism of RdRp. We built an in silico model of the RNA strand embedded in RdRp using experimental methods, starting from the cryo-electron microscopy structure and exploiting the information obtained by spectrometry on the RNA sequence. We determined that the model was stable during the MD simulation time. The obtained results provide deeper insights into the incorporation of nucleoside triphosphates, whose molecular mechanism by the RdRp active site still remains elusive.

Antioxidant and copper-chelating power of new molecules suggested as multiple target agents against Alzheimer's disease. A theoretical comparative study.

In this study, the scavenging activity against OOH radicals and the copper-chelating ability of two new synthesized molecules (named L1 and L2) that can act as multiple target agents against Alzheimer's disease have been investigated at the density functional theory level. The pKa and molar fractions at physiological pH have been predicted. The main antioxidant reaction mechanisms in lipid-like and water environments have been considered and the relative rate constants determined. The copper-chelating ability of the two compounds has also been explored at different coordination sites and computing the complexation kinetic constants. Results show the L1 compound is a more effective radical scavenging and copper-chelating agent than L2.

Chalcogen Effects in the Photophysical Properties of Dimethylamino-1,8-naphthalimide Dyes Revealed by DFT Investigation.

Thionation of carbonyl groups of known dyes is a rapidly emerging strategy to propose an advance toward heavy-atom-free photosensitizers to be used in photodynamic therapy (PDT). The sulfur-for-oxygen replacement has recently proved to enhance the singlet oxygen quantum yield of some existing fluorophores and to shift the absorption band at longer wavelengths. Drawing inspiration from this challenging evidence, the effect of both sulfur- and selenium-for-oxygen replacement in the skeleton of the oxo-4-dimethylamino-1,8-naphthalimide molecule (DMN) has been analyzed by means of a DFT study. The thio- and seleno-derivatives (SDMN and SeDMN, respectively) have been shown to offer the possibility to access a multitude of ISC (intersystem crossing) pathways involved in the triplet deactivation mechanisms with a consequent enhancement of the singlet oxygen production, also arising from the change of orbital type involved in the radiationless 1nπ* → 3ππ* transitions. Moreover, the change in nature from a 1ππ* to a 1nπ* observed in the SeDMN has been revealed to be crucial to reach more clinically useful regions of the spectrum suggesting that the selenium-for-oxygen replacement can be proposed as a strategy to achieve more suitable PDT agents while proposing an advance toward heavy-atom-free PSs.

Computational Study Reveals the Role of Water Molecules in the Inhibition Mechanism of LAT1 by 1,2,3-Dithiazoles.

The L-type amino acid transporter LAT1, involved in many biological processes including the overexpression of some tumors, is considered a potential pharmacological target. The 1,2,3-Dithiazole scaffold was predicted to inhibit LAT1 by the formation of an intermolecular disulfide bond with the thiolate group of cysteine(s). As a result of the identification of these irreversible covalent inhibitors, we decided to deeply investigate the recognition stage and the covalent interaction, characterizing the chemical structures of the selected ligands. With the aim to provide new insights into the access of the ligands to the binding pocket and to reveal the residues involved in the inhibition, we performed docking, molecular dynamics simulations, and density functional theory-based investigation of three 1,2,3-dithiazoles against LAT1. Our computational analysis further highlighted the crucial role played by water molecules in the inhibition mechanism. The results here presented are consistent with experimental observations and provide insights that can be helpful for the rational design of new-to-come LAT1's inhibitors.

Theoretical investigation on bisarylselanylbenzo-2,1,3-selenadiazoles as potential photosensitizers in photodynamic therapy.

Density functional theory and time-dependent (TDDFT) calculations were carried out for recently reported bisarylselanylbenzo-2,1,3-selenadiazoles derivatives capable of producing singlet oxygen (1O2) under UV-Vis irradiation. Conformational behaviors, excitation energies, singlet-triplet energy gaps, and spin-orbit coupling constants were evaluated. The conformational analysis evidences that two different conformers have to be taken into consideration to completely describe the photophysical properties of this class of molecules. TDDFT results show that these compounds, though possessing absorption wavelengths that fall in the violet region, are characterized by singlet-triplet energy gaps greater than the energy required to excite the molecular oxygen, thus being able to produce the cytotoxic species, spin-orbit coupling constants large enough to ensure efficient singlet-triplet intersystem spin crossing, and even the highly reactive superoxide anion O2 •(-) by autoionization and subsequent electron transfer to molecular oxygen in its ground state.

The Se-S Bond Formation in the Covalent Inhibition Mechanism of SARS-CoV-2 Main Protease by Ebselen-like Inhibitors: A Computational Study.

The inhibition mechanism of the main protease (Mpro) of SARS-CoV-2 by ebselen (EBS) and its analog with a hydroxyl group at position 2 of the benzisoselenazol-3(2H)-one ring (EBS-OH) was studied by using a density functional level of theory. Preliminary molecular dynamics simulations on the apo form of Mpro were performed taking into account both the hydrogen donor and acceptor natures of the Nδ and Nε of His41, a member of the catalytic dyad. The potential energy surfaces for the formation of the Se-S covalent bond mediated by EBS and EBS-OH on Mpro are discussed in detail. The EBS-OH shows a distinctive behavior with respect to EBS in the formation of the noncovalent complex. Due to the presence of canonical H-bonds and noncanonical ones involving less electronegative atoms, such as sulfur and selenium, the influence on the energy barriers and reaction energy of the Minnesota hybrid meta-GGA functionals M06, M06-2X and M08HX, and the more recent range-separated hybrid functional wB97X were also considered. The knowledge of the inhibition mechanism of Mpro by the small protease inhibitors EBS or EBS-OH can enlarge the possibilities for designing more potent and selective inhibitor-based drugs to be used in combination with other antiviral therapies.

The Contribution of Density Functional Theory to the Atomistic Knowledge of Electrochromic Processes.

In this review, we provide a brief overview of the contribution that computational studies can offer to the elucidation of the electronic mechanisms responsible for the electrochromism phenomenon, through the use of the density functional theory (DFT) and its time-dependent formulation (TDDFT). Although computational studies on electrochromic systems are not as numerous as those for other physico-chemical processes, we will show their reliability and ability to predict structures, excitation energies, and redox potentials. The results confirm that these methods not only help in the interpretation of experimental data but can also be used for the rational design of molecules with interesting electrochromic properties to be initiated for synthesis and experimental characterization.

A Review on Coordination Properties of Al(III) and Fe(III) toward Natural Antioxidant Molecules: Experimental and Theoretical Insights.

This review focuses on the ability of some natural antioxidant molecules (i.e., hydroxycinnamic acids, coumarin-3-carboxylic acid, quercetin, luteolin and curcumin) to form Al(III)- and Fe(III)-complexes with the aim of evaluating the coordination properties from a combined experimental and theoretical point of view. Despite the contributions of previous studies on the chemical properties and biological activity of these metal complexes involving such natural antioxidants, further detailed relationships between the structure and properties are still required. In this context, the investigation on the coordination properties of Al(III) and Fe(III) toward these natural antioxidant molecules might deserve high interest to design water soluble molecule-based metal carriers that can improve the metal's intake and/or its removal in living organisms.

Iodine substituted phosphorus corrole complexes as possible photosensitizers in photodynamic therapy: Insights from theory.

The search for new dyes to be used as photosensitizers in photodynamic therapy (PDT) is a field of great interest from both experimental and theoretical viewpoints. In this study, the main photophysical properties (excitation energies, singlet-triplet energy gap, and spin orbit coupling matrix elements) of some unsubstituted and iodine substituted phosphorus corrole complexes have been determined by using density functional theory and its time-dependent formulation. Results show that these compounds can be proposed as photosensitizers in PDT. The heavy atom effects have been rationalized on the basis of El-Sayed rules.

How Lanthanide Ions Affect the Addition-Elimination Step of Methanol Dehydrogenases.

The recently discovered methanol dehydrogenase, XoxF, is a widespread enzyme used by methylotrophic bacteria to oxidize methanol for carbon and energy, and requires lanthanide ions for its activity. This enzyme represents an essential component of methanol utilization by both methanol- and methane-utilizing bacteria. The present investigation looks on the electronic, energetic and geometrical behavior of the methanol dehydrogenase from Methylacidiphilum fumariolicum SolV, which is strictly dependent on early lanthanide metals with +3 oxidation states, by examining enzyme-substrate complexes of all the lanthanides. We focus on the catalytic reaction mechanism of two methanol dehydrogenases having as cofactor europium and ytterbium belonging to the mid- and later- series of lanthanides, in comparison with the methanol dehydrogenase containing the cerium, one early lanthanide. Our results provide evidence for the influence of the lanthanide contraction effect in all the elementary steps of the catalytic reaction mechanism. This indication may prove useful for developing new catalytic machineries of enzymes that adopt new-to-nature transformations.

On the Catalytic Activity of the Engineered Coiled-Coil Heptamer Mimicking the Hydrolase Enzymes: Insights from a Computational Study.

Recently major advances were gained on the designed proteins aimed to generate biomolecular mimics of proteases. Although such enzyme-like catalysts must still suffer refinements for improving the catalytic activity, at the moment, they represent a good example of artificial enzymes to be tested in different fields. Herein, a de novo designed homo-heptameric peptide assembly (CC-Hept) where the esterase activity towards p-nitro-phenylacetate was obtained for introduction of the catalytic triad (Cys-His-Glu) into the hydrophobic matrix, is the object of the present combined molecular dynamics and quantum mechanics/molecular mechanics investigation. Constant pH Molecular Dynamics simulations on the apoform of CC-Hept suggested that the Cys residues are present in the protonated form. Molecular dynamics (MD) simulations of the enzyme-substrate complex evidenced the attitude of the enzyme-like system to retain water molecules, necessary in the hydrolytic reaction, in correspondence of the active site, represented by the Cys-His-Glu triad on each of the seven chains, without significant structural perturbations. A detailed reaction mechanism of esterase activity of CC-Hept-Cys-His-Glu was investigated on the basis of the quantum mechanics/molecular mechanics calculations employing a large quantum mechanical (QM) region of the active site. The proposed mechanism is consistent with available esterases kinetics and structural data. The roles of the active site residues were also evaluated. The deacylation phase emerged as the rate-determining step, in agreement with esterase activity of other natural proteases.

The Generation of the Oxidant Agent of a Mononuclear Nonheme Fe(II) Biomimetic Complex by Oxidative Decarboxylation. A DFT Investigation.

The oxidative decarboxylation of the iron(II) α-hydroxy acid (mandelic acid) complex model, biomimetic of Rieske dioxygenase, has been investigated at the density functional level. The explored mechanism sheds light on the role of the α-hydroxyl group on the dioxygen activation. The potential energy surfaces have been explored in different electronic spin states. The rate-determining step of the process is the proton transfer. The oxidative decarboxylation preferentially takes place on the quintet state.