Insights into the Structural Modification of Selenium-Doped Derivatives with Narrowband Emissions: A Theory Study.

The research on boron/nitrogen (B/N)-based multiresonance thermally activated delayed fluorescence (MR-TADF) emitters has been a prominent topic due to their narrowband emission and high luminous efficiency. However, devices derived from the common types of narrowband TADF materials often experience an efficiency roll-off, which could be ascribed to their relatively slow triplet-singlet exciton interconversion. Since inserting the heavy Se atom into the B/N scheme has been a proven strategy to address the abovementioned issues, herein, extensive density functional theory (DFT) and time-dependent DFT (TD-DFT) simulations have been employed to explore the effects of the structural modification on a series of structurally modified selenium-doped derivatives. Furthermore, the two-layered ONIOM (QM/MM) model has been employed to study the pressure effects on the crystal structure and photophysical properties of the pristine CzBSe. The theoretical results found that the introduced tert-butyl units in Cz-BSeN could result in a shorter charge transfer distance and smaller reorganization energy than the parent CzBSe. In contrast to directly incorporating the o-carborane (Cb) unit to CzBSe, incorporating the bridged phenyl units is important in order to achieve narrowband emissions and high luminous efficiency. The lowest three triplet excited states of CzBSe, Cz-BSeN and PhCb-BSeN all contribute to their triplet-singlet exciton conversions, resulting in a high utilization of triplet excitons. The pressure has an evident influence on the photophysical properties of the aggregated CzBSe and is favored for obtaining narrowband emissions. Our work is promised to provide a feasible strategy for designing selenium-doped derivatives with narrowband emissions and rapid triplet-singlet exciton interconversions.

A mechanistic insight for the biosynthesis of N,N-dimethyltryptamine: An ONIOM theoretical approach.

Psychoactive natural products are potent serotonergic agonists capable of modulating brain functions such as memory and cognition. These substances have shown therapeutic potential for treating various mental disorders. The fact that N,N-dimethyltryptamine (DMT) is produced endogenously in several plants and animals, including humans, makes it particularly attractive. As an amino acid-derived alkaloid, the DMT biosynthetic pathway is part of the L-tryptophan biochemical cascade and can be divided into the decarboxylation by an aromatic L-amino acid decarboxylase (AADC) for tryptamine formation and the subsequent double-methylation by the indolethylamine-N-methyltransferase (INMT) through the cofactor S-adenosyl-L-methionine (SAM), a methyl donor. Unlike the decarboxylation mechanism of L-tryptophan, the molecular details of the double methylation of tryptamine have not been elucidated. Therefore, we propose an in silico model using molecular dynamics (MD), non-covalent interaction index (NCI) and density functional theory (DFT) calculations with the ONIOM QM:MM B3LYP/6-31+G(d,p):MM/UFF level of theory. Based on the obtained energetic data, the potential energy surface (PES) indicates an SN2 mechanism profile, with the second methylation energy barrier being the rate-limiting step with δG‡=60kJ∙mol-1 larger than the previous methylation, following the NCI analysis showing more repulsive interactions for the second transition state. In addition, the hybridization information of each reaction step provides geometric details about the double-methylation.

Calculating 13 C NMR chemical shifts of large molecules using the eXtended ONIOM method at high accuracy with a low cost.

Fragmentation-based methods for nuclear magnetic resonance (NMR) chemical shift calculations have become more and more popular in first-principles calculations of large molecules. However, there are many options for a fragmentation-based method to select, such as theoretical methods, fragmentation schemes, the number of levels of theory, etc. It is important to study the optimal combination of the options to achieve a good balance between accuracy and efficiency. Here we investigate different combinations of options used by a fragmentation-based method, the eXtended ONIOM (XO) method, for 13 C chemical shift calculations on a set of organic and biological molecules. We found that: (1) introducing Hartree-Fock exchange into density functional theory (DFT) could reduce the calculation error due to fragmentation in contrast to pure DFT functionals, while a hybrid functional, xOPBE, is generally recommended; (2) fragmentation schemes generated from the molecular tailoring approach (MTA) with small level parameter n, for example, n = 2 and the degree-based fragmentation method (DBFM) with n = 1, are sufficient to achieve satisfactory accuracy; (3) the two-level XO (XO2) NMR calculation is superior to the calculation with only one level of theory, as the second level (i.e., low level) of theory provides a way to well describe the long-range effect. These findings are beneficial to practical applications of fragmentation-based methods for NMR chemical shift calculations of large molecules.

Computational Insights on the Hydride and Proton Transfer Mechanisms of D-Arginine Dehydrogenase.

D-Arginine dehydrogenase from Pseudomonas aeruginosa (PaDADH) is an amine oxidase which catalyzes the conversion of D-arginine into iminoarginine. It contains a non-covalent FAD cofactor that is involved in the oxidation mechanism. Based on substrate, solvent, and multiple kinetic isotope effects studies, a stepwise hydride transfer mechanism is proposed. It was shown that D-arginine binds to the active site of enzyme as α-amino group protonated, and it is deprotonated before a hydride ion is transferred from its α-C to FAD. Based on a mutagenesis study, it was concluded that a water molecule is the most likely catalytic base responsible from the deprotonation of α-amino group. In this study, we formulated computational models based on ONIOM method to elucidate the oxidation mechanism of D-arginine into iminoarginine using the crystal structure of enzyme complexed with iminoarginine. The calculations showed that Arg222, Arg305, Tyr249, Glu87, His 48, and two active site water molecules play key roles in binding and catalysis. Model systems showed that the deprotonation step occurs prior to hydride transfer step, and active site water molecule(s) may have participated in the deprotonation process.

Virtual Screening and Quantum Chemistry Analysis for SARS-CoV-2 RNA-Dependent RNA Polymerase Using the ChEMBL Database: Reproduction of the Remdesivir-RTP and Favipiravir-RTP Binding Modes Obtained from Cryo-EM Experiments with High Binding Affinity.

The novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was identified as the pathogenic cause of coronavirus disease 2019 (COVID-19). The RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 is a potential target for the treatment of COVID-19. An RdRp complex:dsRNA structure suitable for docking simulations was prepared using a cryo-electron microscopy (cryo-EM) structure (PDB ID: 7AAP; resolution, 2.60 Å) that was reported recently. Structural refinement was performed using energy calculations. Structure-based virtual screening was performed using the ChEMBL database. Through 1,838,257 screenings, 249 drugs (37 approved, 93 clinical, and 119 preclinical drugs) were predicted to exhibit a high binding affinity for the RdRp complex:dsRNA. Nine nucleoside triphosphate analogs with anti-viral activity were included among these hit drugs, and among them, remdesivir-ribonucleoside triphosphate and favipiravir-ribonucleoside triphosphate adopted a similar docking mode as that observed in the cryo-EM structure. Additional docking simulations for the predicted compounds with high binding affinity for the RdRp complex:dsRNA suggested that 184 bioactive compounds could be anti-SARS-CoV-2 drug candidates. The hit bioactive compounds mainly consisted of a typical noncovalent major groove binder for dsRNA. Three-layer ONIOM (MP2/6-31G:AM1:AMBER) geometry optimization calculations and frequency analyses (MP2/6-31G:AMBER) were performed to estimate the binding free energy of a representative bioactive compound obtained from the docking simulation, and the fragment molecular orbital calculation at the MP2/6-31G level of theory was subsequently performed for analyzing the detailed interactions. The procedure used in this study represents a possible strategy for discovering anti-SARS-CoV-2 drugs from drug libraries that could significantly shorten the clinical development period for drug repositioning.

The Catalytic Mechanism of the Retaining Glycosyltransferase Mannosylglycerate Synthase.

To protect their intracellular proteins, extremophile microorganisms synthesize molecules called compatible solutes. These molecules are the result of the attachment of a small negatively charged molecule to a sugar molecule. It has been found that these molecules, not only protect the microorganism against osmotic stress but also against other extreme conditions. They can also confer protection against extreme conditions to isolated enzymes from different organisms making them an exciting prospect for potential biotechnological applications. One of the most widespread compatible solute in hyperthermophile organisms is the molecule 2-O-α-D-mannosyl-D-glycerate (MG). In addition to confer protection to proteins against extreme conditions, MG was found to prevent Alzheimer's ß-amyloid aggregation and reduce α-synuclein fibril formation in Parkinson's disease. In this work we studied, using computational methods, the catalytic mechanism of the synthesis of MG by the enzyme mannosylglycerate synthase (MGS) from the thermophilic bacteria Rhodothermus marinus.

Theoretical Study of the Mechanism of Ribosomal Peptide Bond Formation Using the ONIOM Method.

Peptide bond formation in living cells occurs at the peptidyl transferase center (PTC) of the large ribosomal subunit and involves the transfer of the peptidyl group from peptidyl-tRNA to aminoacyl-tRNA. Despite numerous kinetic and theoretical studies, many details of this reaction -such as whether it proceeds via a stepwise or concerted mechanism- remain unclear. In this study, we calculated the geometry and energy of the transition states and intermediates in peptide bond formation in the PTC environment using the ONIOM (our own n-layered integrated molecular orbital and molecular mechanics) method. The calculations indicated that the energy of the transition states of stepwise mechanisms are lower than those of concerted mechanisms and suggested that the reaction involves a neutral tetrahedral intermediate that is stabilized through the hydrogen-bonding network in the PTC environment. The results will lead to a better understanding of the mechanism of peptidyl transfer reaction, and resolve fundamental questions of the steps and molecular intermediates involved in peptide bond formation in the ribosome.

Towards Quantum-Chemical Modeling of the Activity of Anesthetic Compounds.

The modeling of the activity of anesthetics is a real challenge because of their unique electronic and structural characteristics. Microscopic approaches relevant to the typical features of these systems have been developed based on the advancements in the theory of intermolecular interactions. By stressing the quantum chemical point of view, here, we review the advances in the field highlighting differences and similarities among the chemicals within this group. The binding of the anesthetics to their partners has been analyzed by Symmetry-Adapted Perturbation Theory to provide insight into the nature of the interaction and the modeling of the adducts/complexes allows us to rationalize their anesthetic properties. A new approach in the frame of microtubule concept and the importance of lipid rafts and channels in membranes is also discussed.

Prediction of pKa s of Late Transition-Metal Hydrides via a QM/QM Approach.

Three implicit solvation models, the conductor-like polarizable continuum model (C-PCM), the conductor-like screening model (COSMO), and universal implicit solvent model (SMD), combined with a hybrid two layer QM/QM approach (ONIOM), were utilized to calculate the pKa values, using a direct thermodynamic scheme, of a set of Group 10 transition metal (TM) hydrides in acetonitrile. To obtain the optimal combination of quantum methods for ONIOM calculations with implicit solvation models, the influence of factors, such as the choice of density functional and basis set, the atomic radii used to build a cavity in the solvent, and the size of the model system in an ONIOM scheme, was examined. Additionally, the impact of Grimme's empirical dispersion correction and exact exchange was also investigated. The results were calibrated by experimental data. This investigation provides insight about effective models for the prediction of thermodynamic properties of TM-containing complexes with bulky ligands. © 2019 Wiley Periodicals, Inc.

Elongation method with intermediate mechanical and electrostatic embedding for geometry optimizations of polymers.

The elongation method with intermediate mechanical and electrostatic embedding (ELG-IMEE) is proposed. The electrostatic embedding uses atomic charges generated by a charge sensitivity analysis (CSA) method and parameterized for three different population analyses, namely, the Merz-Singh-Kollman scheme, the charge model 5, and the atomic polar tensor. The obtained CSA models were tested on two model systems. Test calculations show that the electrostatic embedding provides several times of decrease in the difference of energies of testing and reference calculations in comparison with the conventional elongation approach (ELG). The mechanical embedding is implemented in a combination of the conventional elongation method and the ONIOM approach. Moreover, it was demonstrated that the geometry optimization with the ELG-IMEE reduces the errors in the optimized structures by about one order in root-mean-square deviation, when compared to ELG.

fromage: A library for the study of molecular crystal excited states at the aggregate scale.

The study of photoexcitations in molecular aggregates faces the twofold problem of the increased computational cost associated with excited states and the complexity of the interactions among the constituent monomers. A mechanistic investigation of these processes requires the analysis of the intermolecular interactions, the effect of the environment, and 3D arrangements or crystal packing on the excited states. A considerable number of techniques have been tailored to navigate these obstacles; however, they are usually restricted to in-house codes and thus require a disproportionate effort to adopt by researchers approaching the field. Herein, we present the FRamewOrk for Molecular AGgregate Excitations (fromage), which implements a collection of such techniques in a Python library complemented with ready-to-use scripts. The program structure is presented and the principal features available to the user are described: geometrical analysis, exciton characterization, and a variety of ONIOM schemes. Each is illustrated by examples of diverse organic molecules in condensed phase settings. The program is available at https://github.com/Crespo-Otero-group/fromage.

Theoretical analysis of the kinetic isotope effect on carboxylation in RubisCO.

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO) fixes atmospheric carbon dioxide into bioavailable sugar molecules. It is also well known that a kinetic isotope effect (KIE; CO2 carbon atoms) accompanies the carboxylation process. To describe the reaction and the KIE α, two different types of molecular dynamics (MD) simulations (ab initio MD and classical MD) have been performed with an Own N-layered Integrated molecular Orbitals and molecular Mechanics (ONIOM)-hybrid model. A channel structure for CO2 transport has been observed during the MD simulation in RubisCO, and assuming the reaction path from the inlet to the product through the coordinate complex with Mg2+ , simulations have been performed on several molecular configuration models fixing several distances between CO2 and ribulose-1,5-bisphosphate along the channel. Free energy analysis and diffusion coefficient analysis have been evaluated for different phases of the process. It is confirmed that the isotopic fractionation effect for CO2 containing either 13 C or 12 C would appear through the transiting path in the channel structure identified in RubisCO. The estimated isotope fractionation constant was quite close to the experimental value.

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide.

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala25 ). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala25 . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala25 . This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala25 . Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.

Mechanistic insights on the reduction of glutathione disulfide by protein disulfide isomerase.

We explore the enzymatic mechanism of the reduction of glutathione disulfide (GSSG) by the reduced a domain of human protein disulfide isomerase (hPDI) with atomistic resolution. We use classical molecular dynamics and hybrid quantum mechanics/molecular mechanics calculations at the mPW1N/6-311+G(2d,2p):FF99SB//mPW1N/6-31G(d):FF99SB level. The reaction proceeds in two stages: (i) a thiol-disulfide exchange through nucleophilic attack of the Cys53-thiolate to the GSSG-disulfide followed by the deprotonation of Cys56-thiol by Glu47-carboxylate and (ii) a second thiol-disulfide exchange between the Cys56-thiolate and the mixed disulfide intermediate formed in the first step. The Gibbs activation energy for the first stage was 18.7 kcal·mol-1, and for the second stage, it was 7.2 kcal·mol-1, in excellent agreement with the experimental barrier (17.6 kcal·mol-1). Our results also suggest that the catalysis by protein disulfide isomerase (PDI) and thiol-disulfide exchange is mostly enthalpy-driven (entropy changes below 2 kcal·mol-1 at all stages of the reaction). Hydrogen bonds formed between the backbone of His55 and Cys56 and the Cys56-thiol result in an increase in the Gibbs energy barrier of the first thiol-disulfide exchange. The solvent plays a key role in stabilizing the leaving glutathione thiolate formed. This role is not exclusively electrostatic, because an explicit inclusion of several water molecules at the density-functional theory level is a requisite to form the mixed disulfide intermediate. In the intramolecular oxidation of PDI, a transition state is only observed if hydrogen bond donors are nearby the mixed disulfide intermediate, which emphasizes that the thermochemistry of thiol-disulfide exchange in PDI is influenced by the presence of hydrogen bond donors.

NMR and Computational Studies as Analytical and High-Resolution Structural Tool for Complex Hydroperoxides and Diastereomeric Endo-Hydroperoxides of Fatty Acids in Solution-Exemplified by Methyl Linolenate.

A combination of selective 1D Total Correlation Spectroscopy (TOCSY) and 1H-13C Heteronuclear Multiple Bond Correlation (HMBC) NMR techniques has been employed for the identification of methyl linolenate primary oxidation products without the need for laborious isolation of the individual compounds. Complex hydroperoxides and diastereomeric endo-hydroperoxides were identified and quantified. Strongly deshielded C-O-O-H 1H-NMR resonances of diastereomeric endo-hydroperoxides in the region of 8.8 to 9.6 ppm were shown to be due to intramolecular hydrogen bonding interactions of the hydroperoxide proton with an oxygen atom of the five-member endo-peroxide ring. These strongly deshielded resonances were utilized as a new method to derive, for the first time, three-dimensional structures with an assignment of pairs of diastereomers in solution with the combined use of 1H-NMR chemical shifts, Density Functional Theory (DFT), and Our N-layered Integrated molecular Orbital and molecular Mechanics (ONIOM) calculations.

An ONIOM investigation of the effect of conformation on bond dissociation energies in peptides.

In the present study, we use the ONIOM strategy of Morokuma and coworkers to examine the various CH bond dissociation energies (BDEs) of a small peptide (2ONW) and compare these with values obtained for its component individual amino acid residues. To evaluate suitable methods for ONIOM-based geometry optimizations, we test an "internal consistency" approach against full B3-LYP//B3-LYP results, and find B3-LYP/6-31G(d):AM1 to be appropriate. We find that the BDEs at the α-carbon in 2ONW are generally larger than the corresponding values for the individual residues on their own. This is attributed to the constraints of the peptide backbone leading to conformations that are not ideal for captodative stabilization of the resulting α-radicals. At the more flexible ß- and γ-positions, the differences between the BDEs in 2ONW and the individual residues are smaller. Overall, the α-BDEs are smaller than the ß- and γ-BDEs in most cases. Thus, to rationalize the inertness of peptide backbones with respect to α-hydrogen abstraction that is frequently found experimentally, it is necessary to consider alternative protection mechanisms such as the polar effect. © 2018 Wiley Periodicals, Inc.

Hydrolysis of chemically distinct sites of human serum albumin by polyoxometalate: A hybrid QM/MM (ONIOM) study.

In this study, mechanisms of hydrolysis of all four chemically diverse cleavage sites of human serum albumin (HSA) by [Zr(OH)(PW11 O39 )]4- (ZrK) have been investigated using the hybrid two-layer QM/MM (ONIOM) method. These reactions have been proposed to occur through the following two mechanisms: internal attack (IA) and water assisted (WA). In both mechanisms, the cleavage of the peptide bond in the Cys392-Glu393 site of HSA is predicted to occur in the rate-limiting step of the mechanism. With the barrier of 27.5 kcal/mol for the hydrolysis of this site, the IA mechanism is found to be energetically more favorable than the WA mechanism (barrier = 31.6 kcal/mol). The energetics for the IA mechanism are in line with the experimentally measured values for the cleavage of a wide range of dipeptides. These calculations also suggest an energetic preference (Cys392-Glu393, Ala257-Asp258, Lys313-Asp314, and Arg114-Leu115) for the hydrolysis of all four sites of HSA. © 2018 Wiley Periodicals, Inc.

molUP: A VMD plugin to handle QM and ONIOM calculations using the gaussian software.

The notable advances obtained by computational (bio)chemistry provided its widespread use in many areas of science, in particular, in the study of reaction mechanisms. These studies involve a huge number of complex calculations, which are often carried out using the Gaussian suite of programs. The preparation of input files and the analysis of the output files are not easy tasks and often involve laborious and complex steps. Taking this into account, we developed molUP: a VMD plugin that offers a complete set of tools that enhance the preparation of QM and ONIOM (QM/MM, QM/QM, and QM/QM/MM) calculations. The starting structures for these calculations can be imported from different chemical formats. A set of tools is available to help the user to examine or modify any geometry parameter. This includes the definition of layers in ONIOM calculations, choosing fixed atoms during geometry optimizations, the recalculation or adjustment of the atomic charges, performing SCANs or IRC calculations, etc. molUP also extracts the geometries from the output files as well as the energies of each of them. All of these tasks are performed in an interactive GUI that is extremely helpful for the user. MolUP was developed to be easy to handle by inexperienced users, but simultaneously to be a fast and flexible graphical interface to allow the advanced users to take full advantage of this plugin. The program is available, free of charges, for macOS, Linux, and Windows at the PortoBioComp page https://www.fc.up.pt/PortoBioComp/database/doku.php?id=molup. © 2018 Wiley Periodicals, Inc.

DNA and HSA interaction of Vanadium (IV), Copper (II), and Zinc (II) complexes derived from an asymmetric bidentate Schiff-base ligand: multi spectroscopic, viscosity measurements, molecular docking, and ONIOM studies.

The interaction of three complexes [Zn(II), Cu(II), and V(IV)] derived from an asymmetric bidentate Schiff-base ligand with DNA and HSA was studied using fluorescence quenching, UV-Vis spectroscopy, viscosity measurements, and computational methods [molecular docking and our Own N-layered Integrated molecular Orbital and molecular Mechanics (ONIOM)]. The obtained results revealed that the DNA and HSA affinities for binding of the synthesized compounds follow as V(IV) > Zn(II) > Cu(II) and Zn(II) > V(IV) > Cu(II), respectively. The distance between these compounds and HSA was obtained based on the Förster's theory of non-radiative energy transfer. Furthermore, computational molecular docking was carried out to investigate the DNA- and HSA-binding pose of the compounds. Molecular docking calculations showed that H-bond, hydrophobic, and π-cation interactions have dominant role in stability of the compound-HSA complexes. ONIOM method was utilized to investigate the HSA binding of the compounds more precisely in which molecular-mechanics method (UFF) and semi-empirical method (PM6) were selected for the low layer and the high layer, respectively. The results show that the structural parameters of the compounds changed along with binding, indicating the strong interaction between the compounds with HSA and DNA. Viscosity measurements as well as computational docking data suggest that all metal complexes interact with DNA, presumably by groove-binding mechanism.

QM/MM Investigation for Protonation States in a Bilin Reductase PcyA-Biliverdin IXα Complex.

Herein we report quantum mechanical/molecular mechanical (QM/MM) studies to investigate the most probable protonation states of active site amino acids and bound substrate based on a recently reported neutron diffraction structure of phycocyanobilin:ferredoxin oxidoreductase (PcyA) by Unno et al. This structure was considered to be bound in its initial state of biliverdin IXα (BV), which has the C-pyrrole ring in the deprotonated state. The protonation state of BV suggested by neutron and spectroscopic studies is a stable, two-electron reduced complex with a bound hydronium ion. Several ambiguities in the neutron structure were observed which prompted a further theoretical analysis of the structure. This structural investigation provides new understanding of the PcyA and BV protonation states not previously reported in the literature. Our calculations suggest that the hydronium ion (H3 O+ ) is energetically unfavorable, preferentially protonating the neighboring His88 residue and that the C-ring of BV is not protonated.