Binding interactions and in silico ADME prediction of isoconessimine derivatives as potent acetylcholinesterase inhibitors

In pursuit of new acetylcholinesterase (AChE) inhibitors for treating Alzheimer's disease (AD), a series of ten previously synthesized isoconessimine compounds (7a-7j) was in silico investigated for their binding interactions with AChE and pharmacokinetics based on absorption, distribution, metabolism, and excretion (ADME) properties using molecular docking, ONIOM (Our own N-layered Integrated molecular Orbital and molecular Mechanics) method and SwissADME tools. Docking experiments showed that all compounds bind within the active site gorge of AChE (PDB entry 1C2B), posing its aryloxy-substitutional ethyl group to catalytic site and conessine skeleton to peripheral anionic site. ONIOM interaction energy was used as an ONIOM score to improve docking score, and it ranked 7b as the most potent AChE inhibitor, in agreement with previous experiment. Residues, ASP74, TRP86, GLY122, GLU202, TRP286, GLU292, SER293, ILE294, TYR337, TYR341, and HIS447 were identified as important for the binding of the AChE-isoconessimine complex. The SwissADME investigation suggested that four compounds (7a, 7c, 7d and 7f) agree with the rules of drug-likeness. The steric and electronic effects on the aryloxy-substitutional ethyl group as important factors in the AChE inhibition were also discussed, which brings a better understanding of Alzheimer's disease drug development.

Inhibitory potential of furanocoumarins against cyclin dependent kinase 4 using integrated docking, molecular dynamics and ONIOM methods

Cyclin Dependent Kinase 4 (CDK4) is vital in the process of cell-cycle and serves as a G1 phase checkpoint in cell division. Selective antagonists of CDK4 which are in use as clinical chemotherapeutics cause various side-effects in patients. Furanocoumarins induce anti-cancerous effects in a range of human tumours. Therefore, targeting these compounds against CDK4 is anticipated to enhance therapeutic effectiveness. This work intended to explore the CDK4 inhibitory potential of 50 furanocoumarin molecules, using a comprehensive approach that integrates the processes of docking, drug-likeness, pharmacokinetic analysis, molecular dynamics simulations and ONIOM (Our own N-layered Integrated molecular Orbital and Molecular mechanics) methods. The top five best docked compounds obtained from docking studies were screened for subsequent analysis. The molecules displayed good pharmacokinetic properties and no toxicity. Epoxybergamottin, dihydroxybergamottin and notopterol were found to inhabit the ATP-binding zone of CDK4 with substantial stability and negative binding free energy forming hydrogen bonds with key catalytic residues of the protein. Notopterol exhibiting the highest binding energy was subjected to ONIOM calculations wherein the hydrogen bonding interactions were retained with significant negative interaction energy. Hence, through these series of computerised methods, notopterol was screened as a potent CDK4 inhibitor and can act as a starting point in successive processes of drug design.Communicated by Ramaswamy H. Sarma.

Highly robust quantum mechanics and umbrella sampling studies on inclusion complexes of curcumin and β-cyclodextrin

The poor aqueous solubility of curcumin (CUR) obstructs its wide utilization in nutraceuticals, cosmetics, and pharmaceutical companies. This study is dedicated to investigate the stability of CUR inside the hydrophobic pocket of β-cyclodextrin (β-CD), hydroxypropyl-β-CD (HP-β-CD), and 2,6-Di-O-methyl-β-CD (DM-β-CD). Initially, molecular mechanics (MM) calculations and subsequently quantum mechanical (QM) calculations were performed on inclusion complexes to strengthen the MM results. We performed microsecond timescale MD simulations to observe the structural dynamics of CUR inside the cavity of CDs. We elucidated the most stable binding orientations of CUR inside the cavity of CDs based on binding free energy obtained from the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) and umbrella sampling simulations. Furthermore, the two-layered ONIOM (B3LYP/6-311+G(2d,p):PM7) method with CPCM implicit water model was used to derive the complete energetics and thermodynamics of inclusion complexes at 1:1 stoichiometry. Each inclusion reaction was exothermic and spontaneous. The chemical reactivity and kinetic stability of inclusion complexes were described by HOMO-LUMO molecular orbital energies. In conclusion, our studies revealed that HP-β-CD had the highest binding affinity for CUR with the most negative complexation energy (-6520.69 kJ/mol) and Gibb's free energy change (-6448.20 kJ/mol). The atomic-level investigation of noncovalent interactions between CUR and the hydroxypropyl groups in HP-β-CD/CUR complex may be helpful to drive new derivatives of HP-β-CD with better host capacity. The computational strategy adopted here might serve as a benchmark for increasing the solubility of numerous clinically significant molecules.

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.

Dispersion-corrected DFT calculations and umbrella sampling simulations to investigate stability of Chrysin-cyclodextrin inclusion complexes

The study of inclusion complexes of Chrysin (ChR) with three forms of cyclodextrins (CDs) α-, β-, and γ-CD was accomplished to examine the stability of ChR inside the central cavities of CDs. The aim of study was to identify the most suitable form of CD to improve the hydro-solubility of poorly soluble ChR bioactive molecule. Microsecond timescale molecular dynamics (MD) simulations were performed on four inclusion complexes (α-CD/ChR, β-CD/ChR, and two conformations of γ-CD/ChR) to examine the dynamics of ChR inside the cavity of CDs. The first conformation of γ-CD/ChR inclusion complex (γ-CD1/ChR) was identified to possess the highest affinity between host and guest molecule on the basis of binding energy calculated by employing Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) and umbrella sampling simulations. To further strengthen the claims of classical and biased MD studies, Our own N-layered Integrated molecular Orbital and Molecular mechanics (ONIOM) (wB97XD/6-311+g(d,p):pm7) calculations were performed on the selected inclusion complexes. The ONIOM based complexation energy reaffirmed that ChR had highest affinity for the γ-CD1 host molecule. Further, the non-covalent interaction analysis was conducted using Multiwfn software on QM-optimized inclusion complexes with wB97XD/6-311+G(d,p) model chemistry, revealing non-covalent interactions between ChR and CDs. This atomic level information helped us to gain better insights into critical atoms of ChR and CD that participated in intermolecular interactions and identify γ-CD as a suitable host molecule for improving the hydro-solubulity of ChR. The structural insights would help to derive new derivatives of γ-CD with better host capacity.

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.

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.

Structural basis for inhibition of a GH116 β-glucosidase and its missense mutants by GBA2 inhibitors: Crystallographic and quantum chemical study

The crystal structure of the Thermoanaerobacterium xylanolyticum in glycoside hydrolase family 116 (TxGH116) β-glucosidase provides a structural model for human GBA2 glucosylceramidase, an enzyme defective in hereditary spastic paraplegia and a potential therapeutic target for treating Gaucher disease. To assess the therapeutic potential of known inhibitors, the X-ray structure of TxGH116 in complex with isofagomine (IFG) was determined at 2.0 Å resolution and showed the IFG bound in a relaxed chair conformation. The binding of IFG and 7 other iminosugar inhibitors to wild-type and mutant enzymes (Asp508His and Arg786His) mimicking GBA2 pathogenic variants was then evaluated computationally by two-layered ONIOM calculations (at the B3LYP:PM7 level). Calculations showed that six charged residues, Glu441, Asp452, His507, Asp593, Glu777, and Arg786 influence inhibitor binding most. His507, Glu777 and Arg786, form strong hydrogen bonds with the inhibitors (∼1.4-1.6 Å). Thus, the missense mutation of one of these residues in Arg786His has a greater effect on the interaction energies for all inhibitors compared to Asp508His. In line with the experimental data for the inhibitors that have been tested, the favorable interaction energy between the inhibitors and the TxGH116 protein followed the trend: isofagomine > 1-deoxynojirimycin > glucoimidazole > N-butyl-deoxynojirimycin ≈ N-nonyl-deoxynojirimycin > conduritol B epoxide ≈ azepane 1 > azepane 2. The obtained structural and energetic properties and comparison to the GBA2 model can lead to understanding of structural requirement for inhibitor binding in GH116 to aid the design of high potency GBA2 inhibitors.

Quantum chemical studies on the binding domain of SARS-CoV-2 S-protein: human ACE2 interface complex

A two-layer ONIOM(B3LYP/6-31G*:PM7) method is used to model the binding of several drug/drug-like molecules (L) at the SARS-CoV-2 S-protein: human ACE2 protein interface cavity. The selected molecules include a set of thirty-five ligands from the study of Smith and Smith which showed a high docking score in the range of -7.0 to -7.7 kcal/mol and another set of seven repurposing drugs, viz. favipiravir, remdesivir, EIDD, galidesivir, triazavirin, ruxolitinib, and baricitinib. The ONIOM model of the cavity (M) showed a highly polarized electron distribution along its top-to-bottom direction while Ls with lengths in the range 1.0 - 1.5 nm fitted well inside the cavity in a head-to-tail fashion to yield ML complexes. The ligands showed a large variation in the ONIOM-level binding energy (Eb), in the range -2.7 to -85.4 kcal/mol. The Eb of ML complexes better than -40.0 kcal/mol is observed for myricetin, fidarestat, protirelin, m-digallic acid, glucogallin, benserazide hydrochlorideseradie, remdesivir, tazobactum, sapropterin, nitrofurantoin, quinonoid, pyruvic acid calcium isoniazid, and aspartame, and among them the highest Eb -85.4 kcal/mol is observed for myricetin. A hydroxy substitution is suggested for the phenyl ring of aspartame to improve its binding behavior at the cavity, and the resulting ligand 43 showed the best Eb -84.5 kcal/mol. The ONIOM-level study is found to be effective for the interpretation of the noncovalent interactions resulting from residues such as arginine, histidine, tyrosine, lysine, carboxylate, and amide moieties in the active site and suggests rational design strategies for COVID-19 drug development.Communicated by Ramaswamy H. Sarma.

Binding profile of a mixed-ligand silver(I) complex with DNA and Topoisomerase I

A new mixed-ligand Ag(I) complex, [Ag(daf)(phen)]NO₃ (daf = 4,5-diazafluoren-9-one and dian = N-(4,5-diazafluoren-9-ylidene)aniline), was synthesized. The elemental analysis, FTIR, ¹HNMR, UV-Vis spectroscopy, cyclic voltammetry, and DFT (Density Functional Theory) geometry optimization method were applied in order to predict the Ag(I) complex structure which concluded to a distorted tetrahedral N4 coordination around the Ag(I) center. A detailed in silico analysis of the bioaffinity of the complex to DNA and human DNA-Topoisomerase I was conducted using molecular docking simulations and ONIOM (Our own N-layered Integrated molecular Orbital and molecular Mechanics) techniques. In this overall scenario, the results suggest the dominance of π-π stacking interactions of the heteroaromatic ligands in the intercalating pocket of DNA and the active site of the enzyme and the rational correlation between being a good intercalator and a potent Topoisomerase I inhibitor. In vitro DNA-binding experiments by spectrophotometric, spectrofluorometric, Voltammetric, and viscometric techniques at physiological pH also confirmed the computational results. The complex inhibited MCF-7 cell growth in a dose-dependent manner while being nontoxic on HUVEC normal cells.

DFT and ONIOM Simulation of 1,3-Butadiene Polymerization Catalyzed by Neodymium-Based Ziegler-Natta System

Using modern methods of quantum chemistry, a theoretical substantiation of the high cis-stereospecificity of 1,3-butadiene polymerization catalyzed by the neodymium-based Ziegler-Natta system was carried out. For DFT and ONIOM simulation, the most cis-stereospecific active site of the catalytic system was used. By analyzing the total energy, as well as the enthalpy and Gibbs free energy of the simulated catalytically active centers, it was found that the coordination of 1,3-butadiene in the trans-form was more favorable than in the cis-form by 11 kJ/mol. However, as a result of π-allylic insertion mechanism modeling, it was found that the activation energy of cis-1,3-butadiene insertion into the π-allylic neodymium-carbon bond of the terminal group on the reactive growing chain was 10-15 kJ/mol lower than the activation energy of trans-1,3-butadiene insertion. The activation energies did not change when both trans-1,4-butadiene and cis-1,4-butadiene were used for modeling. That is, 1,4-cis-regulation was due not to the primary coordination of 1,3-butadiene in its cis-configuration, but to its lower energy of attachment to the active site. The obtained results allowed us to clarify the mechanism of the high cis-stereospecificity of 1,3-butadiene polymerization by the neodymium-based Ziegler-Natta system.

A multiscale ONIOM study of the buckminsterfullerene (C60) Diels-Alder reaction: from model design to reaction path analysis

The hybrid ONIOM (Our own N-layered Integrated molecular Orbital and molecular Mechanics) formalism is employed to investigate the Diels-Alder reaction of the buckminsterfullerene C₆₀. Our computations suggest that the ONIOM2(M06-2X/6-31G(d): SVWN/STO-3G) model, enclosing both the diene and the pyracyclene fragment of C60 in the higher-layer, provides a reasonable trade-off between accuracy and computational cost as it comes to predicting reaction energetics. Moreover, the frontier molecular orbital (FMO) theory and activation strain model (ASM) are jointly relied on to rationalize the effect of -OH and -CN substituents on the activation barrier of this reaction. Finally, reaction paths are scrutinized to get insight into the various forces underpinning the process of cycloadduct formation.

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.

A new step in kinetic proofreading due to misacylated-tRNA during ribosomal peptide bond formation

The translational accuracy in protein synthesis is contributed to by several mechanisms in the ribosome, generally called kinetic proofreading. This process in the ribosome inhibits the non-cognate codon-anticodon interaction. However, it is not sufficient for fidelity of protein synthesis since a wrong amino acid can easily be added to the growing polypeptide chain if a tRNA while cognate to the mRNA, carries a non-cognate amino acid. Therefore, additional to the kinetic proofreading, there must be some hitherto unknown characteristic in misacylated-tRNAs to stop the process of protein synthesis if such misacylated-tRNA is accommodated in the ribosomal A-site. In order to understand this characteristic, we have performed computational quantum chemistry analysis on five different tRNA molecules, each one attached to five different amino acids with one being cognate to the tRNA and the other four non-cognate. This study shows the importance of aminoacyl-tRNA binding energy in ensuring fidelity of protein synthesis.

An insight into the interaction between α-ketoamide- based inhibitor and coronavirus main protease: A detailed in silico study

The search for therapeutic drugs that can neutralize the effects of COVID-2019 (SARS-CoV-2) infection is the main focus of current research. The coronavirus main protease (M_pro) is an attractive target for anti-coronavirus drug design. Further, α-ketoamide is proved to be very effective as a reversible covalent-inhibitor against cysteine proteases. Herein, we report on the non-covalent to the covalent adduct formation mechanism of α-ketoamide-based inhibitor with the enzyme active site amino acids by QM/SQM model (QM = quantum mechanical, SQM = semi-empirical QM). To uncover the mechanism, we focused on two approaches: a concerted and a stepwise fashion. The concerted pathway proceeds via deprotonation of the thiol of cysteine (here, Cys₁₄₅ SγH) and simultaneous reversible nucleophilic attack of sulfur onto the α-ketoamide warhead. In this work, we propose three plausible concerted pathways. On the contrary, in a traditional two-stage pathway, the first step is proton transfer from Cys₁₄₅ SγH to His₄₁ Nδ forming an ion pair, and consecutively, in the second step, the thiolate ion attacks the α-keto group to form a thiohemiketal. In this reaction, we find that the stability of the tetrahedral intermediate oxyanion/hydroxyl group plays an important role. Moreover, as the α-keto group has two faces Si or Re for the nucleophilic attack, we considered both possibilities of attack leading to S- and R-thiohemiketal. We computed the structural, electronic, and energetic parameters of all stationary points including transition states via ONIOM and pure DFT method. Additionally, to characterize covalent, weak noncovalent interaction (NCI) and hydrogen-bonds, we applied NCI-reduced density gradient (NCI-RDG) methods along with Bader's Quantum Theory of Atoms-in-Molecules (QTAIM) and natural bonding orbital (NBO) analysis.

Mechanistic study of L-6-hydroxynicotine oxidase by DFT and ONIOM methods

L-6-Hydroxynicotine oxidase (LHNO) is a member of monoamine oxidase (MAO) family and catalyzes conversion of (S)-6-hydroxynicotine to 6-hydroxypseudooxynicotine during bacterial degradation of nicotine. Recent studies indicated that the enzyme catalyzes oxidation of carbon-nitrogen bond instead of previously proposed carbon-carbon bond. Based on kinetics and mutagenesis studies, Asn166, Tyr311, and Lys287 as well as an active site water molecule have roles in the catalysis of the enzyme. A number of studies including experimental and computational methods support hydride transfer mechanism in MAO family as a common mechanism in which a hydride ion transfer from amine substrate to flavin cofactor is the rate-limiting step. In this study, we formulated computational models to study the hydride transfer mechanism using crystal structure of enzyme-substrate complex. The calculations involved ONIOM and DFT methods, and we evaluated the geometry and energetics of the hydride transfer process while probing the roles of active site residues. Based on the calculations involving hydride, radical, and polar mechanisms, it was concluded that hydride transfer mechanism is the only viable mechanism for LHNO.

819 molecular knot: a theoretical analysis of the electronic structure using an ONIOM approach

The present work analyzes the electronic and molecular properties of the 8₁₉ ([Fe(II)₄]Cℓ) and metal-free knot ligand complexes obtained from X-ray crystal structure of molecular 8₁₉ knot complex [Fe(II)₄(PF₆)₇]Cℓ. The studies were theoretically investigated by means of DFT, TD-DFT, and ONIOM approaches. Basis sets functions from all-electron calculations for bromine, iodine, and iron atoms were adapted to be used along with relativistic effective core potential, while H, C, N, O, and Cℓ atoms were described by Pople basis sets. The diffusion effect of halogen into the 8₁₉ cavity, UV-Vis, and Electronic Circular Dichroism spectra were also analyzed. All calculations were performed using solvent effect through the SCRF/SMD model and dispersion effects by Grimme methodology. The value of mean separation distance between Cℓ and iron atom (7.218 Å) is in good agreement with X-ray experimental result (7.258 Å). Circular dichroism spectrum of metal-free 8₁₉ knot ligand was calculated and the maximum absorption in 262 nm, Δ𝜖 obtained was 67 L mol^- 1 cm^- 1. These results are qualitatively similar to those obtained experimentally, 295 nm and 80 L mol^- 1 cm^- 1, respectively. In this study, we report the electronic and molecular properties of the 8₁₉ ([Fe(II)₄]Cl and metal-free knot ligand complexes and compare with the results obtained from X-ray crystallographic data of 8₁₉ knot complex [Fe(II)₄(PF₆)₇]Cl. The 8₁₉ knot were investigated by means of DFT, TD-DFT, and ONIOM approaches. Basis sets functions from all-electron for Br, I, and Fe atoms were adapted to be used along with relativistic effective core potential, while H, C, N, O, and Cl atoms were described by Pople basis sets. The objective was to understand the stability of the 8₁₉ knot as a function of the substitution of the central halogen atom (Cl), and the signal in the circular dichroism spectra. From the equilibrium geometries, we have obtained good results for values of the bond distance, bond angle, and dihedral angle along the molecular structure when these variables are compared with the results obtained from X-ray data. The diffusion effect of halogen into the 8₁₉ cavity, UV-Vis, and Electronic Circular Dichroism spectra was also analyzed. Circular dichroism spectrum of metal-free 8₁₉ knot ligand was calculated, and the maximum absorption is in good agreement with the experimental value. The ONIOM methodology combined with the relativistic effective core potential and the atomic basis sets provide good results for systems with a complex topology, such as knots.

Assessment of the orientation and conformation of pigments in protein binding sites from infrared difference spectra

Time resolved FTIR difference spectroscopy (DS) has been used to study photosystem I (PSI) with the disubstituted 1,4-naphthoquinones acequinocyl (AcQ) and lapachol (Lpc) incorporated into the A₁ binding site. AcQ is a 2-acetoxy-3-dodecyl-1,4-naphthoquinone, Lpc is a 2-hydroxy-3-(3-methyl-2-butenyl)-1,4-naphthoquinone. To assess whether the experimental spectra are specific to different orientations of the quinone and their substitutions ONIOM-type QM/MM vibrational frequency calculations were undertaken for various orientations of the pigments and side-chain conformations in the A₁ binding site. Comparison of calculated and experimental spectra for the reduced species (semiquinone anion) suggests that the orientation for the naphthoquinone ring in the binding site and specific side-chain conformations can be identified based on the spectra. In native PSI phylloquinone (PhQ) in the A₁ binding site binds with its phytyl chain ortho to the hydrogen bonded carbonyl group. This is not found to be the case for the hydrocarbon tail of AcQ, which is meta to the H-bonded carbonyl group. In contrast, Lpc in PSI binds with its hydrocarbon tail also ortho to the H-bonded carbonyl group. Furthermore, comparison of calculated and experimental spectra indicates which conformations the acetoxy group of AcQ and the hydroxy group of Lpc adopt in the A₁ binding site.

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.

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.

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 (Ala₂₅ ). 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 Ala₂₅ . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅ . 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 Ala₂₅ . 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.

Two-Photon Absorption Cross-Sections in Fluorescent Proteins Containing Non-canonical Chromophores Using Polarizable QM/MM

Multi-photon absorption properties, particularly two-photon absorption (2PA), of fluorescent proteins (FPs) have made them attractive tools in deep-tissue clinical imaging. Although the diversity of photophysical properties for FPs is wide, there are some caveats predominant among the existing FP variants that need to be overcome, such as low quantum yields and small 2PA cross-sections. From a computational perspective, Salem et al. (2016) suggested the inclusion of non-canonical amino acids in the chromophore of the red fluorescent protein DsRed, through the replacement of the tyrosine amino acid. The 2PA properties of these new non-canonical chromophores (nCCs) were determined in vacuum, i.e., without taking into account the protein environment. However, in the computation of response properties, such as 2PA cross-sections, the environment plays an important role. To account for environment and protein-chromophore coupling effects, quantum mechanical/molecular mechanical (QM/MM) schemes can be useful. In this work, the polarizable embedding (PE) model is employed along with time-dependent density functional theory to describe the 2PA properties of a selected set of chromophores made from non-canonical amino acids as they are embedded in the DsRed protein matrix. The objective is to provide insights to determine whether or not the nCCs could be developed and, thus, generate a new class of FPs. Results from this investigation show that within the DsRed environment, the nCC 2PA cross-sections are diminished relative to their values in vacuum. However, further studies toward understanding the 2PA limit of these nCCs using different protein environments are needed.

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 pK_a 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.

Towards chemically accurate QM/MM simulations on GPUs

Computational chemistry simulations are extensively used to model natural phenomena. To maintain performance similar to molecular mechanics, but achieve comparable accuracy to quantum mechanical calculations, many researchers are using hybrid QM/MM methods. In this article we evaluate our GPU-accelerated ONIOM implementation by measurements on the crambin and HIV integrase proteins with different size QM model systems. We demonstrate that by using a larger QM region, a better energy accuracy can be achieved at the expense of simulation time. This trade-off is important to consider for the researcher running QM/MM calculations. Furthermore, we show that the ONIOM energy monotonically approaches the pure quantum mechanical energy of the whole system. The experiments are made feasible by utilizing the cutting-edge BrianQC quantum chemistry module for Hartree-Fock level SCF and our GPU-accelerated MMFF94 force field implementation for molecular mechanics calculations.

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.

TriplatinNC and Biomolecules: Building Models Based on Non-covalent Interactions

The class of polynuclear platinum(II) compounds have demonstrated a great interest because their high activity against cancer cells. Among these new compounds, the TriplatinNC also called AH78, demonstrated surprising antitumor activity, in some cases equivalent to cisplatin. It is well-known that complex charge +8 favors interaction with DNA and other biomolecules non-covalently, through the hydrogen bonds with phosphate and sulfate groups present in these structures. The hydrogen atoms of the amine interact with the oxygen atoms of the phosphate and sulfate groups present in the DNA strand and heparan sulfate, respectively. These interactions can cause significant twists in double helix and inhibit the activity of these biomolecules. The present investigation is an attempt to provide a benchmark theoretical study about TriplatinNC. We have described the non-covalent interactions through small reliable mimetic models. The non-covalent interactions were also evaluated on larger models containing DNA fractions with six nitrogenous base pairs (CGCGAA) and fractions of the disaccharide that makes the HS evaluated by the hybrid QM/MM ONIOM methodology.

On the encapsulation of Olsalazine by β-cyclodextrin: A DFT-based computational and spectroscopic investigations

In this work, the supramolecular host-guest interaction of the prodrug Olsalazine (OLZ) and β-Cyclodextrin (β-CD) was examined experimentally and computationally. Experimentally, employing the UV-Vis spectroscopic method in aqueous media at various pH's, results obtained using the Benesi-Hilderbrand approach demonstrated that OLZ can form supramolecular inclusion complex with β-CD with stoichiometric ratio of 1:1. Furthermore, these results revealed that the formation of OLZ: β-CD complexes exhibited insignificant pH dependency in the range 5-8 with an average binding constant (K_b) of approximately 1×10³M^-1. Computationally, geometry optimization of 1:1 OLZ: β-CD complexes was performed employing the ONIOM (DFT((ωB97XB)/6-31+G(d)),SQM(PM3)) approach. Obtained results demonstrated that OLZ: β-CD complex is stabilized by the formation of intermolecular hydrogen bonds with an average length of approximately 1.8Å. Additionally, the stability of OLZ: β-CD complex was demonstrated employing ADMP molecular dynamic simulations over a timeframe of 500fs. The molecularity of the supramolecular host-guest interaction between OLZ and β-CD is presented and interpreted in the essence of TD-DFT and molecular orbitals analyses.

Revisiting the mechanism for the polar hydrochlorination of alkenes

Various kinetic studies of the addition of hydrogen halides to alkenes were carried out in the 1930s, 1940s, and 1970s. Since then, there have been theoretical analyses of several aspects of alkene reactivity and regioselectivity during hydrohalogenation. A few works have studied the influence of the hydrogen halide when it is acting as a catalyst (as well as a reactant) under third-order kinetics. However, there has not been any theoretical investigation of the possibility that there is a carbocation intermediate in the contact ion pair. Therefore, we revisited the mechanistic aspects of this reaction at different levels of theory and using IRC (in a gaseous medium) and ONIOM (with an explicit solvent model). We found that (i) there is only one transition state (one energy barrier) for all possible reaction mechanisms in an apolar medium, (ii) there is no carbocation intermediate when a hydrohalogenation reaction is performed in apolar and polar protic media, (iii) hydrochlorination in apolar and polar protic media occurs through an asynchronous concerted mechanism (not through the stepwise mechanism described in undergraduate/graduate literature), and (iv) there are three possible competing mechanisms (with second-, third-, and fourth-order kinetics) in an apolar medium; the mechanism with fourth-order kinetics has the smallest energy barrier, while that with second-order kinetics (a mechanism in which the hydrogen halide does not also act as a catalyst) has the highest energy barrier. Graphical abstract Hydrochlorination of 2,3-dimethylbutene through an asynchronous concerted mechanism.

Recognition of repulsive and attractive regions of selected serotonin receptor binding site using FMO-EDA approach

The complexes of selected long-chain arylpiperazines with homology models of 5-HT_1A, 5-HT_2A, and 5-HT₇ receptors were investigated using quantum mechanical methods. The molecular geometries of the ligand-receptor complexes were firstly optimized with the Our own N-layered Integrated molecular Orbital and molecular Mechanics (ONIOM) method. Next, the fragment molecular orbitals method with an energy decomposition analysis scheme (FMO-EDA) was employed to estimate the interaction energies in binding sites. The results clearly showed that orthosteric binding sites of studied serotonin receptors have both attractive and repulsive regions. In the case of 5-HT_1A and 5-HT_2A two repulsive areas, located in the lower part of the binding pocket, and one large area of attraction engaging many residues at the top of all helices were identified. Additionally, for the 5-HT₇ receptor, the third area of destabilization located at the extracellular end of the helix 6 was found.

Inhibition mechanism of L,D-transpeptidase 5 in presence of the β-lactams using ONIOM method

Tuberculosis (TB) is one of the world's deadliest diseases resulting from infection by the bacterium, Mycobacterium tuberculosis (M.tb). The L,D-transpeptidase enzymes catalyze the synthesis of 3 → 3 transpeptide linkages which are predominant in the peptidoglycan of the M.tb cell wall. Carbapenems is class of β-lactams that inactivate L,D-transpeptidases by acylation, although differences in antibiotic side chains modulate drug binding and acylation rates. Herein, we used a two-layered our Own N-layer integrated Molecular Mechanics ONIOM method to investigate the catalytic mechanism of L,D-transpeptidase 5 (Ldt_Mt5) by β-lactam derivatives. Ldt_Mt5 complexes with six β-lactams, ZINC03788344 (1), ZINC02462884 (2), ZINC03791246 (3), ZINC03808351 (4), ZINC03784242 (5) and ZINC02475683 (6) were simulated. The QM region (high-level) comprises the β-lactam, one water molecule and the Cys360 catalytic residue, while the rest of the Ldt_Mt5 residues were treated with AMBER force field. The activation energies (ΔG^#) were calculated with B3LYP, M06-2X and ωB97X density functionals with 6-311++G(2d, 2p) basis set. The ΔG^# for the acylation of Ldt_Mt5 by the selected β-lactams were obtained as 13.67, 20.90, 22.88, 24.29, 27.86 and 28.26 kcal mol^-1respectively. Several of the compounds showed an improved ΔG^# when compared to the previously calculated energies for imipenem and meropenem for the acylation step for Ldt_Mt5. This model provides further validation of the catalytic inhibition mechanism of LDTs with atomistic detail.

Can Adsorption on Graphene be Used for Isotopic Enrichment? A DFT Perspective

We have explored the theoretical applicability of adsorption on graphene for the isotopic enrichment of aromatic compounds. Our results indicate that for nonpolar molecules, like benzene, the model compound used in these studies shows a reasonable isotopic fractionation that is obtained only for the deuterated species. For heavier elements, isotopic enrichment might be possible with more polar compounds, e.g., nitro- or chloro-substituted aromatics. For benzene, it is also not possible to use isotopic fractionation to differentiate between different orientations of the adsorbed molecule over the graphene surface. Our results also allowed for the identification of theory levels and computational procedures that can be used for the reliable prediction of the isotope effects on adsorption on graphene. In particular, the use of partial Hessian is an attractive approach that yields acceptable values at an enormous increase of speed.

Do Better Quality Embedding Potentials Accelerate the Convergence of QM/MM Models? The Case of Solvated Acid Clusters

This study examines whether the use of more accurate embedding potentials improves the convergence of quantum mechanics/molecular mechanics (QM/MM) models with respect to the size of the QM region. In conjunction with density functional theory calculations using the ωB97X-D functional, various embedding potentials including the TIP3P water model, the effective fragment potential (EFP), and semi-empirical methods (PM6, PM7, and DFTB) were used to simulate the deprotonation energies of solvated acid clusters. The calculations were performed on solvated neutral (HA) and cationic (HB⁺) acids clusters containing 160 and 480 water molecules using configurations sampled from molecular dynamics simulations. Consistently, the ωB97X-D/EFP model performed the best when using a minimal QM region size. The performance for the other potentials appears to be highly sensitive to the charge character of the acid/base pair. Neutral acids display the expected trend that semi-empirical methods generally perform better than TIP3P; however, an opposite trend was observed for the cationic acids. Additionally, electronic embedding provided an improvement over mechanical embedding for the cationic systems, but not the neutral acids. For the best performing ωB97X-D/EFP model, a QM region containing about 6% of the total number of solvent molecules is needed to approach within 10 kJ mol^-1 of the pure QM result if the QM region was chosen based on the distance from the reaction centre.

Insights Into Receptor-Based Anesthetic Pharmacophores and Anesthetic-Protein Interactions

General anesthetics are thought to allosterically bind and potentiate the inhibitory currents of the GABA_A receptor through drug-specific binding sites. The physiologically relevant isoform of the GABA_A receptor is a transmembrane ligand-gated ion channel consisting of five subunits (γ-α-β-α-β linkage) symmetrically arranged around a central chloride-conducting pore. Although the exact molecular structure of this heteropentameric GABA_A receptor remains unknown, molecular modeling has allowed significant advancements in understanding anesthetic binding and action. Using the open-channel conformations of the homologous glycine and glutamate-gated chloride receptors as templates, a homology model of the GABA_A receptor was constructed using the Discovery Studio computational chemistry software suite. Consensus structural alignment of the homology templates allowed for the construction of a three-dimensional heteropentameric GABA_A receptor model with (γ₂-β₃-α₁-β₃-α₁) subunit linkage. An anesthetic binding site was identified within the transmembrane α/β intersubunit space by the convergence of three residues shown to be essential for anesthetic activity in previous studies with mutant mice (β₃-N265, β₃-M286, α₁-L232). Propofol derivatives docked into this binding site showed log-linear correlation with experimentally derived GABA_A receptor potentiation (EC₅₀) values, suggesting this binding site may be important for receptor activation. The receptor-based pharmacophore was analyzed with surface maps displaying the predominant anesthetic-protein interactions, revealing an amphiphilic binding cavity incorporating the three residues involved in anesthetic modulation. Quantum mechanics calculations of the bonding patterns found in complementary high-resolution receptor systems further elucidated the complex nature of anesthetic-protein interactions.

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.

Theoretical insight into the solvent effect of H2O and formamide on the cooperativity effect in HMX complex

The cooperativity effects of the H-bonding interactions in HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane)∙∙∙HMX∙∙∙FA (formamide), HMX∙∙∙HMX∙∙∙H₂O and HMX∙∙∙HMX∙∙∙HMX complexes involving the chair and chair-chair HMX are investigated by using the ONIOM2 (CAM-B3LYP/6-31++G(d,p):PM3) and ONIOM2 (M06-2X/6-31++G(d,p):PM3) methods. The solvent effect of FA or H₂O on the cooperativity effect in HMX∙∙∙HMX∙∙∙HMX are evaluated by the integral equation formalism polarized continuum model. The results show that the cooperativity and anti-cooperativity effects are not notable in all the systems. Although the effect of solvation on the binding energy of ternary system HMX∙∙∙HMX∙∙∙HMX is not large, that on the cooperativity of H-bonds is notable, which leads to the mutually strengthened H-bonding interaction in solution. This is perhaps the reason for the formation of different conformation of HMX in different solvent. Surface electrostatic potential and reduced density gradient are used to reveal the nature of the solvent effect on cooperativity effect in HMX∙∙∙HMX∙∙∙HMX. Graphical abstract RDG isosurface and electrostatic potential surface of HMX∙∙∙HMX∙∙∙HMX.

Multistructural microiteration technique for geometry optimization and reaction path calculation in large systems

We propose a multistructural microiteration (MSM) method for geometry optimization and reaction path calculation in large systems. MSM is a simple extension of the geometrical microiteration technique. In conventional microiteration, the structure of the non-reaction-center (surrounding) part is optimized by fixing atoms in the reaction-center part before displacements of the reaction-center atoms. In this method, the surrounding part is described as the weighted sum of multiple surrounding structures that are independently optimized. Then, geometric displacements of the reaction-center atoms are performed in the mean field generated by the weighted sum of the surrounding parts. MSM was combined with the QM/MM-ONIOM method and applied to chemical reactions in aqueous solution or enzyme. In all three cases, MSM gave lower reaction energy profiles than the QM/MM-ONIOM-microiteration method over the entire reaction paths with comparable computational costs. © 2017 Wiley Periodicals, Inc.

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.

Modifications of the chromophore of Spinach aptamer based on QM:MM calculations

Spinach aptamer was developed as an RNA analog of the green fluorescent protein. The aptamer interacts with its ligand and modifies its electronic spectrum so that it fluoresces brightly at the wavelength of 501 nm. Song et al. investigated modifications of the ligand in their experimental study and found a molecule emitting at 523 nm upon creating a complex with the Spinach aptamer. The crystal structure of the aptamer in complex with its original ligand has been published, which enabled us to study the system computationally. In this article, we suggest several new modifications of the ligand that shift the emission maximum of the complex to even longer wavelengths. Our results are based on combined quantum mechanical/molecular mechanical calculations with DFT method used for geometry optimization and TD-DFT for calculations of absorption and emission energies.

Two Binding Geometries for Risperidone in Dopamine D3 Receptors: Insights on the Fast-Off Mechanism through Docking, Quantum Biochemistry, and Molecular Dynamics Simulations

Risperidone is an atypical antipsychotic used in the treatment of schizophrenia and of symptoms of irritability associated with autism spectrum disorder (ASD). Its main action mechanism is the blockade of D2-like receptors acting over positive and negative symptoms of schizophrenia with small risk of extrapyramidal symptoms (EPS) at doses corresponding to low/moderate D2 occupancy. Such a decrease in the side effect incidence can be associated with its fast unbinding from D2 receptors in the nigrostriatal region allowing the recovery of dopamine signaling pathways. We performed docking essays using risperidone and the D3 receptor crystallographic data and results suggested two possible distinct orientations for risperidone at the binding pocket. Orientation 1 is more close to the opening of the binding site and has the 6-fluoro-1,2 benzoxazole fragment toward the bottom of the D3 receptor cleft, while orientation 2 is deeper inside the binding pocket with the same fragment toward to the receptor surface. In order to unveil the implications of these two binding orientations, classical molecular dynamics and quantum biochemistry computations within the density functional theory formalism and the molecular fractionation with conjugate caps framework were performed. Quantum mechanics/molecular mechanics suggests that orientation 2 (considering the contribution of Glu90) is slightly more energetically stable than orientation 1 with the main contribution coming from residue Asp110. The residue Glu90, positioned at the opening of the binding site, is closer to orientation 1 than 2, suggesting that it may have a key role in stability through attractive interaction with risperidone. Therefore, although orientations 1 and 2 are both likely to occur, we suggest that the occurrence of the first may contribute to the reduction of side effects in patients taking risperidone due to the reduction of dopamine receptor occupancy in the nigrostriatal region through a mechanism of fast dissociation. The atypical effect may be obtained simply by either delaying D3R full blockage by spatial hindrance of orientation 1 at the binding site or through an effective blockade followed by orientation 1 fast dissociation. While the molecular interpretation suggested in this work shed some light on the potential molecular mechanisms accounting for the reduced extrapyramidal symptoms observed during risperidone treatment, further studies are necessary in order to evaluate the implications of both orientations during the receptor activation/inhibition. Altogether these data highlight important hot spots in the dopamine receptor binding site bringing relevant information for the development of novel/derivative agents with atypical profile.

Vibrational analysis on the revised potential energy curve of the low-barrier hydrogen bond in photoactive yellow protein

Photoactive yellow protein (PYP) has a characteristic hydrogen bond (H bond) between p-coumaric acid chromophore and Glu46, whose OH bond length has been observed to be 1.21 Å by the neutron diffraction technique [Proc. Natl. Acad. Sci. 106, 440-4]. Although it has been expected that such a drastic elongation of the OH bond could be caused by the quantum effect of the hydrogen nucleus, previous theoretical computations including the nuclear quantum effect have so far underestimated the bond length by more than 0.07 Å. To elucidate the origin of the difference, we performed a vibrational analysis of the H bond on potential energy curve with O…O distance of 2.47 Å on the equilibrium structure, and that with O…O distance of 2.56 Å on the experimental crystal structure. While the vibrationally averaged OH bond length for equilibrium structure was underestimated, the corresponding value for crystal structure was in reasonable agreement with the corresponding experimental values. The elongation of the O…O distance by the quantum mechanical or thermal fluctuation would be indispensable for the formation of a low-barrier hydrogen bond in PYP.

Theoretical study of the mechanism of protein arginine deiminase 4 (PAD4) inhibition by F-amidine

Protein arginine deiminase 4 (PAD4) catalyzes the hydrolysis of a peptidylarginine residue to form a citrulline residue and ammonia during posttranslational modification. This process plays a pivotal role in rheumatoid arthritis (RA) and gene regulation. F-amidine belongs to a series of haloacetamidine compounds that are the most potent PAD4 inhibitors described to date. F-amidine acts as a mechanism-based inhibitor of PAD4, inactivating PAD4 by the covalent modification of the active site Cys645. In this manuscript, the fundamental mechanism of PAD4 inhibition by F-amidine is investigated using a QM/MM approach. Our simulations show that in the PAD4-F-amidine reactant complex, the active site Cys645 exists as a thiolate and His471 is protonated. This is consistent with the reverse protonation mechanism wherein the active site nucleophile, Cys645, in PAD4 exists as a thiolate in the active form of the enzyme. Inhibition of PAD4 by F-amidine is initiated by the nucleophilic addition of Sγ to the Cζ of F-amidine, leading to the formation of a tetrahedral intermediate. His471 serves as a proton donor, helping F to leave the fluoroacetamidine moiety of F-amidine; meanwhile, Sγ forms a three-membered ring with Cζ and Cη of F-amidine. Subsequently, the three-membered sulfonium ring collapses and rearranges to the final thioether product. His471 acts as a proton donor in the transition state and facilitates the inhibition reaction of PAD4.

Comparison of calculated and experimental isotope edited FTIR difference spectra for purple bacterial photosynthetic reaction centers with different quinones incorporated into the QA binding site

Previously we have shown that ONIOM type (QM/MM) calculations can be used to simulate isotope edited FTIR difference spectra for neutral ubiquinone in the QA binding site in Rhodobacter sphaeroides photosynthetic reaction centers. Here we considerably extend upon this previous work by calculating isotope edited FTIR difference spectra for reaction centers with a variety of unlabeled and (18)O labeled foreign quinones incorporated into the QA binding site. Isotope edited spectra were calculated for reaction centers with 2,3-dimethoxy-5,6-dimethyl-1,4-benzoquinone (MQ0), 2,3,5,6-tetramethyl-1, 4-benzoquinone (duroquinone, DQ), and 2,3-dimethyl-l,4-naphthoquinone (DMNQ) incorporated, and compared to corresponding experimental spectra. The calculated and experimental spectra agree well, further demonstrating the utility and applicability of our ONIOM approach for calculating the vibrational properties of pigments in protein binding sites. The normal modes that contribute to the bands in the calculated spectra, their composition, frequency, and intensity, and how these quantities are modified upon (18)O labeling, are presented. This computed information leads to a new and more detailed understanding/interpretation of the experimental FTIR difference spectra. Hydrogen bonding to the carbonyl groups of the incorporated quinones is shown to be relatively weak. It is also shown that there is some asymmetry in hydrogen bonding, accounting for 10-13 cm(-1) separation in the frequencies of the carbonyl vibrational modes of the incorporated quinones. The extent of asymmetry in H-bonding could only be established by considering the spectra for various types of quinones incorporated into the QA binding site. The quinones listed above are "tail-less." Spectra were also calculated for reaction centers with corresponding "tail" containing quinones incorporated, and it is found that replacement of the quinone methyl group by a phytyl or prenyl chain does not alter ONIOM calculated spectra.

Binding affinity of substituted ureido-benzenesulfonamide ligands to the carbonic anhydrase receptor: a theoretical study of enzyme inhibition

The binding properties of a series of benzenesulfonamide inhibitors (4-substituted-ureido-benzenesulfonamides, UBSAs) of human carbonic anhydrase II (hCA II) enzyme with active site residues have been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) model. To account for the important docking interactions between the UBSAs ligand and hCA II enzyme, a molecular docking program AutoDock Vina is used. The molecular docking results obtained by AutoDock Vina revealed that the docked conformer has root mean square deviation value less than 1.50 Å compared to X-ray crystal structures. The inhibitory activity of UBSA ligands against hCA II is found to be in good agreement with the experimental results. The thermodynamic parameters for inhibitor binding show that hydrogen bonding, hydrophilic, and hydrophobic interactions play a major role in explaining the diverse inhibitory range of these derivatives. Additionally, natural bond orbital analysis is performed to characterize the ligand-metal charge transfer stability. The insights gained from this study have great potential to design new hCA-II inhibitor, 4-[3-(1-p-Tolyl-4-trifluoromethyl-1H-pyrazol-3-yl)-ureido]-benzenesulfonamide, which belongs to the family of UBSA inhibitors and shows similar type of inhibitor potency with hCA II. This work also reveals that a QM/MM model and molecular docking method are computationally feasible and accurate for studying substrate-protein inhibition.

Computational study on the conformation and vibration frequencies of β-sheet of ε-polylysine in vacuum

Two oligomers, each containing 3 l-lysine residues, were used as model molecules for the simulation of the β-sheet conformation of ɛ-polylysine (ɛ-PLL) chains. Their C terminals were capped with ethylamine and N terminals were capped with α-l-aminobutanoic acid, respectively. The calculations were carried out with the hybrid two-level ONOIM (B3LYP/6-31G:PM3) computational chemistry method. The optimized conformation was obtained and IR frequencies were compared with experimental data. The result indicated that the two chains were winded around each other to form a distinct cyclohepta structure through bifurcated hydrogen bonds. The groups of amide and α-amidocyanogen coming from one chain and the carbonyl group from the other chain were involved in the cyclohepta structure. The bond angle of the bifurcated hydrogen bonds was 66.6°. The frequency analysis at ONIOM [B3LYP/6-31G (d):PM3] level showed the IR absorbances of the main groups, such as the amide and amidocyanogen groups, were in accordance with the experimental data.