1
|
He W, Yang H, Li Y, Cui Y, Wei L, Xu T, Li Y, Zhang M. Identifying the toxic mechanisms of emerging electronic contaminations liquid crystal monomers and the construction of a priority control list for graded control. THE SCIENCE OF THE TOTAL ENVIRONMENT 2024; 951:175398. [PMID: 39128516 DOI: 10.1016/j.scitotenv.2024.175398] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2024] [Revised: 08/05/2024] [Accepted: 08/06/2024] [Indexed: 08/13/2024]
Abstract
Liquid crystal monomers (LCMs) are identified as emerging organic contaminations with largely unexplored health impacts. To elucidate their toxic mechanisms, support the establishment of environmental discharge and management standards, and promote effective LCMs control, this study constructs a database covering 20,545 potential targets of 1431 LCMs, highlighting 9 key toxic target proteins that disrupt the nervous system and metabolic functions. GO and KEGG pathway analysis suggests LCMs severely affect nervous system, linked to neurodegenerative diseases and mental health disorders, with toxicity variations driven by electronegativity and structural complexity of LCM terminal groups. To achieve tiered control of LCMs, construct toxicity risk control lists for 9 key toxic target proteins, suitable for the graded control of LCMs, management recommendations are provided based on toxicity levels. These lists were validated for reliability and offer reliable toxicity predictions for LCMs. SHAP analysis points to electronic properties, molecular shape, and structural characteristics of LCMs as primary health impact factors. As the first study integrating machine learning with computational toxicology to outline LCMs health impacts, it aims to enhance public understanding of LCM toxicity risks and support the development of environmental standards, effective management of LCM production and emissions, and reduction of public exposure risks.
Collapse
Affiliation(s)
- Wei He
- MOE Key Laboratory of Resources Environmental Systems Optimization, North China Electric Power University, Beijing, China
| | - Hao Yang
- MOE Key Laboratory of Resources Environmental Systems Optimization, North China Electric Power University, Beijing, China
| | - Yunxiang Li
- MOE Key Laboratory of Resources Environmental Systems Optimization, North China Electric Power University, Beijing, China
| | - Yuhan Cui
- MOE Key Laboratory of Resources Environmental Systems Optimization, North China Electric Power University, Beijing, China
| | - Luanxiao Wei
- MOE Key Laboratory of Resources Environmental Systems Optimization, North China Electric Power University, Beijing, China
| | - Tingzhi Xu
- MOE Key Laboratory of Resources Environmental Systems Optimization, North China Electric Power University, Beijing, China.
| | - Yu Li
- MOE Key Laboratory of Resources Environmental Systems Optimization, North China Electric Power University, Beijing, China
| | - Meng Zhang
- College of Environmental Sciences and Engineering, State Environmental Protection Key Laboratory of All Material Fluxes in River Ecosystems, Peking University, Beijing 100871, China; The Key Laboratory of Water and Sediment Sciences, Ministry of Education, International Joint Laboratory for Regional Pollution Control, Ministry of Education, Beijing 100871, China.
| |
Collapse
|
2
|
Boča R, Imrich R, Štofko J, Vranovičová B, Rajnák C. Molecular properties of linear amino acids in water. Amino Acids 2024; 56:5. [PMID: 38300332 PMCID: PMC10834582 DOI: 10.1007/s00726-023-03365-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 12/13/2023] [Indexed: 02/02/2024]
Abstract
Four linear amino acids of increased separation of the carboxyl and amino groups, namely glycine (aminoacetic acid), β-alanine (3-aminopropanoic acid), GABA (4-aminobutanoic acid) and DAVA (5-aminopentanoic acid), have been studied by quantum chemical ab initio and DFT methods including the solvent effect in order to get electronic structure and molecular descriptors, such as ionisation energy, electron affinity, molecular electronegativity, chemical hardness, electrophilicity index, dipole moment, quadrupole moment and dipole polarizability. Thermodynamic functions (zero-point energy, inner energy, enthalpy, entropy, and the Gibbs energy) were evaluated after the complete vibrational analysis at the true energy minimum provided by the full geometry optimization. Reaction Gibbs energy allows evaluating the absolute redox potentials on reduction and/or oxidation. The non-local non-additive molecular descriptors were compared along the series showing which of them behave as extensive, varying in match with the molar mass and/or separation of the carboxyl and amino groups. Amino acidic forms and zwitterionic forms of the substances were studied in parallel in order to compare their relative stability and redox properties. In total, 24 species were investigated by B3LYP/def2-TZVPD method (M1) including neutral molecules, molecular cations and molecular anions. For comparison, MP2/def2-TZVPD method (M2) with full geometry optimization and vibrational analysis in water has been applied for 12 species; analogously, for 24 substances, DLPNO-CCSD(T)/aug-cc-pVTZ method (M3) has been applied in the geometry obtained by MP2 and/or B3LYP. It was found that the absolute oxidation potential correlates with the adiabatic ionisation energy; the absolute reduction potential correlates with the adiabatic electron affinity and the electrophilicity index. In order to validate the used methodology with experimental vertical ionisation energies and vibrational spectrum obtained in gas phase, calculations were done also in vacuo.
Collapse
Affiliation(s)
- Roman Boča
- Faculty of Health Sciences, University of SS Cyril and Methodius, 91701, Trnava, Slovakia.
| | - Richard Imrich
- Faculty of Health Sciences, University of SS Cyril and Methodius, 91701, Trnava, Slovakia
| | - Juraj Štofko
- Faculty of Health Sciences, University of SS Cyril and Methodius, 91701, Trnava, Slovakia
| | - Beáta Vranovičová
- Faculty of Natural Sciences, University of SS Cyril and Methodius, 91701, Trnava, Slovakia
| | - Cyril Rajnák
- Faculty of Natural Sciences, University of SS Cyril and Methodius, 91701, Trnava, Slovakia
| |
Collapse
|
3
|
Isert C, Atz K, Riniker S, Schneider G. Exploring protein-ligand binding affinity prediction with electron density-based geometric deep learning. RSC Adv 2024; 14:4492-4502. [PMID: 38312732 PMCID: PMC10835705 DOI: 10.1039/d3ra08650j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 01/19/2024] [Indexed: 02/06/2024] Open
Abstract
Rational structure-based drug design relies on accurate predictions of protein-ligand binding affinity from structural molecular information. Although deep learning-based methods for predicting binding affinity have shown promise in computational drug design, certain approaches have faced criticism for their potential to inadequately capture the fundamental physical interactions between ligands and their macromolecular targets or for being susceptible to dataset biases. Herein, we propose to include bond-critical points based on the electron density of a protein-ligand complex as a fundamental physical representation of protein-ligand interactions. Employing a geometric deep learning model, we explore the usefulness of these bond-critical points to predict absolute binding affinities of protein-ligand complexes, benchmark model performance against existing methods, and provide a critical analysis of this new approach. The models achieved root-mean-squared errors of 1.4-1.8 log units on the PDBbind dataset, and 1.0-1.7 log units on the PDE10A dataset, not indicating significant advantages over benchmark methods, and thus rendering the utility of electron density for deep learning models context-dependent. The relationship between intermolecular electron density and corresponding binding affinity was analyzed, and Pearson correlation coefficients r > 0.7 were obtained for several macromolecular targets.
Collapse
Affiliation(s)
- Clemens Isert
- ETH Zurich, Department of Chemistry and Applied Biosciences Vladimir-Prelog-Weg 4 8093 Zurich Switzerland +41 44 633 73 27
| | - Kenneth Atz
- ETH Zurich, Department of Chemistry and Applied Biosciences Vladimir-Prelog-Weg 4 8093 Zurich Switzerland +41 44 633 73 27
| | - Sereina Riniker
- ETH Zurich, Department of Chemistry and Applied Biosciences Vladimir-Prelog-Weg 4 8093 Zurich Switzerland +41 44 633 73 27
| | - Gisbert Schneider
- ETH Zurich, Department of Chemistry and Applied Biosciences Vladimir-Prelog-Weg 4 8093 Zurich Switzerland +41 44 633 73 27
| |
Collapse
|
4
|
Mei Y, Chen X, Wei R, Chang XY, Tao L, Liu LL. An Isolable Radical Anion Featuring a 2-Center-3-Electron π-Bond without a Clearly Defined σ-Bond. Angew Chem Int Ed Engl 2023; 62:e202315555. [PMID: 37942957 DOI: 10.1002/anie.202315555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 11/08/2023] [Accepted: 11/09/2023] [Indexed: 11/10/2023]
Abstract
Featuring an extra electron in the π* antibonding orbital, species with a 2-center-3-electron (2c3e) π bond without an underlying σ bond are scarcely known. Herein, we report the synthesis, isolation and characterization of a radical anion salt [K(18-C-6)]+ {[(HCNDipp)2 Si]2 P2 }⋅- (i.e. [K(18-C-6)]+ 3⋅- ) (18-C-6=18-crown-6, Dipp=2,6-diisopropylphenyl), in which 3⋅- features a perfectly planar Si2 P2 four-membered ring. This species represents the first example of a Si- and P-containing analog of a bicyclo[1.1.0]butane radical anion. The unusual bonding motif of 3⋅- was thoroughly investigated via X-ray diffraction crystallography, electron paramagnetic resonance spectroscopy (EPR), and calculations by density functional theory (DFT), which collectively unveiled the existence of a 2c3e π bond between the bridgehead P atoms and no clearly defined supporting P-P σ bond.
Collapse
Affiliation(s)
- Yanbo Mei
- Department of Chemistry and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen, 518055, China
- Department of Chemistry and Dongguan Key Laboratory for Data Science and Intelligent Medicine, Great Bay University, Dongguan, 523000, China
| | - Xiaodan Chen
- Guangdong Provincial Key Laboratory of Functional Supramolecular Coordination Materials and Applications, College of Chemistry and Materials Science, Jinan University, Guangzhou, 510632, China
| | - Rui Wei
- Department of Chemistry and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen, 518055, China
| | - Xiao-Yong Chang
- Department of Chemistry and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen, 518055, China
| | - Lizhi Tao
- Department of Chemistry and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen, 518055, China
| | - Liu Leo Liu
- Department of Chemistry and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen, 518055, China
- State Key Laboratory of Elemento-Organic Chemistry, Nankai University, Tianjin, 300071, China
| |
Collapse
|
5
|
Lan X, Zhang X, Mei Y, Hu C, Liu LL. Utilizing bis(imino)dihydroacridanide pincer ligands in p-block chemistry: synthesis and catalysis of an antimony monocation salt. Dalton Trans 2023; 52:15660-15664. [PMID: 37859530 DOI: 10.1039/d3dt03310d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2023]
Abstract
We present the synthesis and characterization of an Sb(III) monocation salt stabilized by a bulky bis(imino)dihydroacridanide pincer ligand. The Lewis acidity of the Sb cation is quantified using the Guttmann-Beckett method and confirmed by its reaction with 4-dimethylaminopyridine, which forms a Lewis acid-base adduct. This Sb cation exhibits catalytic activity in the cyanosilylation of arylketones. The electronic structure of the Sb cation as well as the mechanism of the catalytic transformation are explored by density functional theory computations.
Collapse
Affiliation(s)
- Xiaofang Lan
- Department of Chemistry, Guangdong Provincial Key Laboratory of Catalysis and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen 518055, China.
| | - Xin Zhang
- Department of Chemistry, Guangdong Provincial Key Laboratory of Catalysis and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen 518055, China.
| | - Yanbo Mei
- Department of Chemistry, Guangdong Provincial Key Laboratory of Catalysis and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen 518055, China.
| | - Chaopeng Hu
- Department of Chemistry, Guangdong Provincial Key Laboratory of Catalysis and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen 518055, China.
| | - Liu Leo Liu
- Department of Chemistry, Guangdong Provincial Key Laboratory of Catalysis and Research Center for Chemical Biology and Omics Analysis, College of Science, Southern University of Science and Technology, Shenzhen 518055, China.
| |
Collapse
|
6
|
Lefrancois-Gagnon KM, Mawhinney RC. Toward Universal Substituent Constants: Relating QTAIM Functional Group Descriptors to Substituent Effect Proxies. J Chem Inf Model 2023; 63:6068-6080. [PMID: 37729015 DOI: 10.1021/acs.jcim.3c00987] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/22/2023]
Abstract
Substituents modulate reactions, but their effects are commonly described by using proxies to their functional group properties. Substituent descriptors from the quantum theory of atoms in molecules, which are true functional group properties, are related here to these proxies, which have historically had chemically relevant meaning. Due to the large number of descriptors, multivariate analysis is used to intuit their significance. Multiple linear regression, principal component, and partial least squares regression analyses highlight that these substituent descriptors contain similar information to the proxies while being intrinsic, predictable substituent properties. Sources of error limiting quantitative reproduction of the proxy data include transferability, experimental accuracy, and solvation issues.
Collapse
Affiliation(s)
| | - Robert C Mawhinney
- Department of Chemistry, Lakehead University, 955 Oliver Road, Thunder Bay, Ontario P7B 5E1, Canada
| |
Collapse
|
7
|
Petrova VV, Domnin AV, Porozov YB, Kuliaev PO, Solovev YV. Implementation of machine learning protocols to predict the hydrolysis reaction properties of organophosphorus substrates using descriptors of electron density topology. J Comput Chem 2023. [PMID: 37772443 DOI: 10.1002/jcc.27227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2023] [Revised: 07/25/2023] [Accepted: 08/28/2023] [Indexed: 09/30/2023]
Abstract
Prediction of catalytic reaction efficiency is one of the most intriguing and challenging applications of machine learning (ML) algorithms in chemistry. In this study, we demonstrated a strategy for utilizing ML protocols applied to Quantum Theory of Atoms In Molecules (QTAIM) parameters to predict the ability of the A17 L47K catalytic antibody to covalently capture organophosphate pesticides. We found that the novel "composite" DFT functional B97-3c could be effectively employed for fast and accurate initial geometry optimization, aligning well with the input dataset creation. QTAIM descriptors proved to be well-established in describing the examined dataset using density-based and hierarchical clustering algorithms. The obtained clusters exhibited correlations with the chemical classes of the input compounds. The precise physical interpretation of the QTAIM properties simplifies the explanation of feature impact for both supervised and unsupervised ML protocols. It also enables acceleration in the search for entries with desired properties within large databases. Furthermore, our findings indicated that Ridge Regression with Laplacian kernel and CatBoost Regressor algorithms demonstrated suitable performance in handling small datasets with non-trivial dependencies. They were able to predict the actual reaction barrier values with a high level of accuracy. Additionally, the CatBoost Classifier proved reliable in discriminating between "active" and "inactive" compounds.
Collapse
Affiliation(s)
- Vlada V Petrova
- M.M. Shemyakin and Yu.A, Ovchinnikov Institute of Bioorganic Chemistry of the Russian Academy of Sciences, Moscow, Russia
- Quantum Chemistry Department, Institute of Chemistry, St. Petersburg State University, Saint Petersburg, Russia
| | - Anton V Domnin
- Quantum Chemistry Department, Institute of Chemistry, St. Petersburg State University, Saint Petersburg, Russia
| | - Yuri B Porozov
- St. Petersburg School of Physics, Mathematics, and Computer Science, HSE University, Saint Petersburg, Russia
- The Center of Bio- and Chemoinformatics, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| | - Pavel O Kuliaev
- Independent Researcher from Saint Petersburg, Saint Petersburg, Russia
| | - Yaroslav V Solovev
- M.M. Shemyakin and Yu.A, Ovchinnikov Institute of Bioorganic Chemistry of the Russian Academy of Sciences, Moscow, Russia
| |
Collapse
|
8
|
Kostal J. Making the Case for Quantum Mechanics in Predictive Toxicology─Nearly 100 Years Too Late? Chem Res Toxicol 2023; 36:1444-1450. [PMID: 37676849 DOI: 10.1021/acs.chemrestox.3c00171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/09/2023]
Abstract
The use of quantum mechanics (QM) has long been the norm to study covalent-binding phenomena in chemistry and biochemistry. The pharmaceutical industry leverages QM models explicitly in covalent drug discovery and implicitly to characterize short-range interactions in noncovalent binding. Predictive toxicology has resisted widespread adoption of QM, including in the pharmaceutical industry, despite its obvious relevance to the metabolic processes in the upstream of adverse outcome pathways and advances in both QM methods and computational resources, which support fit-for-purpose applications in reasonable timeframes. Here, we make the case for embracing QM as an indispensable part of a toxicologist's toolkit. We argue that QM provides the necessary orthogonality to alert-based expert systems and traditional QSARs, consistent with calls for animal-free integrated testing strategies for safety assessments of commercial chemicals. We outline existing roadblocks to this transition, including the need to train model developers in QM and the shift toward service-based toxicity models that utilize high-performance computing clusters. Lastly, we describe recent examples of successful implementations of QM in hazard assessments and propose how in silico toxicology can be further advanced by integrating QM with artificial intelligence.
Collapse
Affiliation(s)
- Jakub Kostal
- Designing Out Toxicity (DOT) Consulting LLC, 2121 Eisenhower Avenue, Alexandria, Virginia 22314, United States
- The George Washington University, 800 22nd Street NW, Washington, DC, 20052, United States
| |
Collapse
|
9
|
Bhowmik P, Baezid HM, Arabi II. Assessment of antidiabetic activity of three Phenylspirodrimanes from fungus Stachybotrys chartarum MUT 3308 by ADME, QSAR, molecular docking and molecular dynamics simulation studies against protein tyrosine phosphatase 1B (PTP1B). J Biomol Struct Dyn 2023:1-15. [PMID: 37698508 DOI: 10.1080/07391102.2023.2256410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2023] [Accepted: 08/24/2023] [Indexed: 09/13/2023]
Abstract
Phenylspirodrimanes (PSD) are the sesquiterpene quinone type meroterpenoids found in nature. PSDs are found to exhibit inhibitory activity against immunocomplex diseases, and tyrosine kinase receptors. Three of the different PSDs C1, C2, and C3 that have been reported to be isolated from the sponge-associated fungus Stachybotrys chartarum MUT 3308 are selected for studying their possible inhibitory effect against type 2 diabetes mellitus. Mechanistically, blocking protein tyrosine phosphatase 1B (PTP1B) helps to reduce the insulin resistance induction caused by the high expression of PTP1B. The QSAR, ADME, toxicity (T) study was carried out to predict the pharmacokinetic properties and the biological activities of the PSDs. PASS prediction web tool was used to find and select the target proteins 1NNY, and 2HNP. According to the molecular docking simulations, C1 and C2 showed better binding affinity of -8.5 kcal/mol, and -8.1 kcal/mol respectively against 1NNY compared to the control ligand. RMSD, RMSF, Rg, and SASA analysis revealed that both C1, and C2 showed better stability, minor conformational changes, and minor fluctuation upon binding to PTP1B. Protein contact analysis was carried out to validate the residues that are in contact with the ligands according to molecular docking studies. Overall, C1, and C2 could be proposed as novel natural hits to be developed and small modifications of these PSDs could result in inducing the binding affinity significantly, although experimental validation is required for further evaluation of the work.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Prasenjit Bhowmik
- Department of Chemistry-BMC, Biochemistry, Disciplinary Domain of Science and Technology, Uppsala University, Uppsala, Sweden
- Department of Chemistry, Faculty of Science, University of Chittagong, Chittagong, Bangladesh
- Department of Textile Engineering, Green University of Bangladesh, Narayanganj, Bangladesh
| | - Hossain Mohammad Baezid
- Department of Genetic Engineering and Biotechnology, Faculty of Biological Sciences, University of Chittagong, Chittagong, Bangladesh
| | - Ishmam Ibnul Arabi
- Department of Chemistry, Faculty of Science, University of Chittagong, Chittagong, Bangladesh
- Department of Textile Engineering, Green University of Bangladesh, Narayanganj, Bangladesh
| |
Collapse
|
10
|
Hoelm M, Adamczyk J, Wzgarda-Raj K, Palusiak M. Effect of a Substituent on the Properties of Salicylaldehyde Hydrazone Derivatives. J Org Chem 2023; 88:2132-2139. [PMID: 36735741 PMCID: PMC9942203 DOI: 10.1021/acs.joc.2c02547] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The present study investigates the effect of the substitution of salicylaldehyde hydrazones at two selected positions, i.e., the para-position with regard to the proton-donating and proton-accepting centers forming the hydrogen bridge. A detailed analysis of structural data obtained by theoretical approaches and X-ray experiments, together with original resonance Hammett's constants, indicates that the strength of the intramolecular hydrogen bonding present in salicylaldehyde hydrazones can be selectively modulated by substitution of the parent molecular system with the chemical group of known π-electron-donating or -accepting properties. Our findings provide an insight into planning synthesis pathways for salicylaldehyde hydrazone species and predicting their result with regard to their H-bonding and related physical and chemical properties.
Collapse
Affiliation(s)
- Marta Hoelm
- Department
of Physical Chemistry, Faculty of Chemistry, University of Lodz, Pomorska 163/165, Lodz90-236, Poland
| | - Justyna Adamczyk
- Department
of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, Lodz91-403, Poland
| | - Kinga Wzgarda-Raj
- Department
of Physical Chemistry, Faculty of Chemistry, University of Lodz, Pomorska 163/165, Lodz90-236, Poland
| | - Marcin Palusiak
- Department
of Physical Chemistry, Faculty of Chemistry, University of Lodz, Pomorska 163/165, Lodz90-236, Poland,
| |
Collapse
|
11
|
Reddy Manne M, Panicker RR, Ramakrishnan K, Hareendran HMK, Kumar Pal S, Kumar S, Pallepogu R, Desikan R, Sivaramakrishna A. Synthesis and Biological Evaluation of a Series of Quinoline‐Based Quinazolinones and Carbamic Anhydride Derivatives. ChemistrySelect 2023. [DOI: 10.1002/slct.202204508] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Madhava Reddy Manne
- Department of Chemistry School of Advanced Sciences Vellore Institute of Technology (VIT) Vellore 632 014 Tamil Nadu India
| | - Rakesh R Panicker
- Department of Chemistry School of Advanced Sciences Vellore Institute of Technology (VIT) Vellore 632 014 Tamil Nadu India
| | - Kumar Ramakrishnan
- Department of Chemistry School of Advanced Sciences Vellore Institute of Technology (VIT) Vellore 632 014 Tamil Nadu India
| | - Hima M. K. Hareendran
- Department of Chemistry School of Advanced Sciences Vellore Institute of Technology (VIT) Vellore 632 014 Tamil Nadu India
| | - Sudhir Kumar Pal
- Center for Bio-separation Technology Vellore Institute of Technology Vellore 632014 Tamil Nadu India
| | - Sanjit Kumar
- Center for Bio-separation Technology Vellore Institute of Technology Vellore 632014 Tamil Nadu India
| | - Raghavaiah Pallepogu
- Department of Chemistry Central University of Karnataka Kadaganchi Kalaburagi – 585 367 Karnataka India
| | - Rajagopal Desikan
- Department of Chemistry School of Advanced Sciences Vellore Institute of Technology (VIT) Vellore 632 014 Tamil Nadu India
| | - Akella Sivaramakrishna
- Department of Chemistry School of Advanced Sciences Vellore Institute of Technology (VIT) Vellore 632 014 Tamil Nadu India
| |
Collapse
|
12
|
How does electron exchange correlation influences reactivity of metallo-β-lactamase L1 against cephalosporin antibiotics. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111774] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
13
|
Iramain MA, Hidalgo JR, Sundius T, Brandán SA. A combined study on structures and vibrational spectra of the antiviral rimantadine using SQMFF and DFT calculations. Heliyon 2022; 8:e10102. [PMID: 36016527 PMCID: PMC9396557 DOI: 10.1016/j.heliyon.2022.e10102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2022] [Revised: 05/21/2022] [Accepted: 07/25/2022] [Indexed: 01/18/2023] Open
Abstract
In this research, a combined study on structures and vibrational spectra of antiviral rimantadine have been performed using hybrid B3LYP/6-311++G∗∗ calculations and the scaled quantum force field (SQMFF) procedure. Harmonic force fields and scaled force constants of Free Base (FB), Cationic (CA) and Hydrochloride (HCl) species derived from the antiviral rimantadine have been calculated in gas phase and in aqueous solution using normal internal coordinates and scaling factors. Good correlations were acquired comparing the theoretical IR, Raman, 1H- 13C-NMR and UV spectra of three species with the analogous experimental ones, suggesting probably, the presence of all them in both phases. The main force constants of three species have evidenced lower values than the corresponding to antiviral amantadine. The ionic character of N1-H33⋯Cl36 bond of HCl species in aqueous solution evidence positive Mulliken charge on N1 atom indicating that this species is as CA one. Rimantadine presents higher solvation energies in water than other antiviral species, such as chloroquin, niclosamide, cidofovir and brincidofovir. The FB and HCl species of rimantadine are slightly less reactive than the corresponding to amantadine while the opposite is observed for the CA species. The predicted ECD spectra for the FB and CA species show positive Cotton effect different from the negative observed for the HCl one. These different behaviours of three species of rimantadine could probably explain the differences observed in the intensities of bands predicted in the electronic spectra of these species.
Collapse
Affiliation(s)
- Maximiliano A. Iramain
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica, Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471 (4000), San Miguel de Tucumán, Tucumán, Argentina
| | - José Ruiz Hidalgo
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica, Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471 (4000), San Miguel de Tucumán, Tucumán, Argentina
| | - Tom Sundius
- Department of Physics, University of Helsinki, Finland
| | - Silvia Antonia Brandán
- Cátedra de Química General, Instituto de Química Inorgánica, Facultad de Bioquímica, Química y Farmacia, Universidad Nacional de Tucumán, Ayacucho 471 (4000), San Miguel de Tucumán, Tucumán, Argentina
| |
Collapse
|
14
|
Jara-Cortés J, Matta CF, Hernández-Trujillo J. A fast approximate extension of the interacting quantum atoms energy decomposition to excited states. J Comput Chem 2022; 43:1068-1078. [PMID: 35470908 DOI: 10.1002/jcc.26863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Revised: 02/25/2022] [Accepted: 03/23/2022] [Indexed: 11/11/2022]
Abstract
An approach is developed for the fast calculation of the interacting quantum atoms energy decomposition (IQA) from the information contained in the first order reduced density matrix only. The proposed methodology utilizes an approximate exchange-correlation density from Density Matrix Functional Theory without the need to evaluate the correlation-exchange contribution directly. Instead, weight factors are estimated to decompose the exact Vxc into atomic and pairwise contributions. In this way, the sum of the IQA contributions recovers the energy obtained from the electronic structure calculation. This method can, hence, be applied to obtain atomic contributions in excited states on the same footing as in their ground states using any method that delivers the reduced first-order density matrix. In this way, one can locate chromophores from first principles quantum chemical calculations. Test calculations on the ground and excited states of a set of small molecules indicate that the scaled atomic contributions reproduce vertical electronic transition energies calculated exactly. This approach may be useful to extend the applicability of the IQA approach in the study of large photochemical systems especially when the calculations of the second order reduced density matrices is prohibitive or not possible.
Collapse
Affiliation(s)
- Jesús Jara-Cortés
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic, Mexico
| | - Chérif F Matta
- Department of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia, Canada
| | - Jesús Hernández-Trujillo
- Departamento de Física y Química Teórica, Facultad de Química, UNAM. Circuito Escolar, Ciudad Universitaria, Mexico City, Mexico
| |
Collapse
|
15
|
Structural Insights and Supramolecular Description of Gliclazide and its Impurity F. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113707] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Lefrancois-Gagnon KM, Mawhinney RC. Toward universal substituent constants: Transferability of group descriptors from the quantum theory of atoms in molecules. J Comput Chem 2022; 43:265-278. [PMID: 34842294 DOI: 10.1002/jcc.26787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Revised: 09/29/2021] [Accepted: 11/06/2021] [Indexed: 11/06/2022]
Abstract
Traditionally, substituents are described not by their intrinsic properties, but by their effect elsewhere in a molecule. However, the quantum theory of atoms in molecules (QTAIM) provides a route to intrinsic substituent descriptors. Ideally, these descriptors would exhibit minimal change as the local environment changes, and hence be considered transferable. Whether this is true is an open question. Here, we evaluated the transferability of QTAIM functional group descriptors for 117 functional groups in a series of 17 substrates to determine whether descriptors obtained using hydrogen as substrate are transferable. The functional group volume has a strong, consistent, linear relationship throughout. All other hydrogen-based group descriptors exhibit a relatively strong linear relationship with those in carbon-based substrates and a reasonable linear relationship with those in non-carbon-based substrates. Outliers are readily interpretable in terms of substrate induced functional group geometry changes. As expected, directional descriptors lying along the substituent-substrate axis are the least conserved.
Collapse
|
17
|
Gallegos M, Guevara-Vela JM, Pendás ÁM. NNAIMQ: A neural network model for predicting QTAIM charges. J Chem Phys 2022; 156:014112. [PMID: 34998318 DOI: 10.1063/5.0076896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Atomic charges provide crucial information about the electronic structure of a molecular system. Among the different definitions of these descriptors, the one proposed by the Quantum Theory of Atoms in Molecules (QTAIM) is particularly attractive given its invariance against orbital transformations although the computational cost associated with their calculation limits its applicability. Given that Machine Learning (ML) techniques have been shown to accelerate orders of magnitude the computation of a number of quantum mechanical observables, in this work, we take advantage of ML knowledge to develop an intuitive and fast neural network model (NNAIMQ) for the computation of QTAIM charges for C, H, O, and N atoms with high accuracy. Our model has been trained and tested using data from quantum chemical calculations in more than 45 000 molecular environments of the near-equilibrium CHON chemical space. The reliability and performance of NNAIMQ have been analyzed in a variety of scenarios, from equilibrium geometries to molecular dynamics simulations. Altogether, NNAIMQ yields remarkably small prediction errors, well below the 0.03 electron limit in the general case, while accelerating the calculation of QTAIM charges by several orders of magnitude.
Collapse
Affiliation(s)
- Miguel Gallegos
- Depto. Química Física y Analítica, Universidad de Oviedo, 33006 Oviedo, Spain
| | - José Manuel Guevara-Vela
- Institute of Chemistry, National Autonomous University of Mexico, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán, Mexico City C.P. 04510, Mexico
| | - Ángel Martín Pendás
- Depto. Química Física y Analítica, Universidad de Oviedo, 33006 Oviedo, Spain
| |
Collapse
|
18
|
Momen R, Azizi A, Morales-Bayuelo A, Pazhoohesh M, Ji X. New insights of QTAIM and stress tensor to finding non-competitive/competitive torquoselectivity of cyclobutene. J Chem Phys 2021; 155:204305. [PMID: 34852485 DOI: 10.1063/5.0068694] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This study aims to investigate the phenomenon of torquoselectivity through three thermal cyclobutene ring-opening reactions (N1-N3). This research focuses on the nature of the chemical bond, electronic reorganization, predicting non-competitive or competitive reactions, and torquoselectivity preference within Quantum Theory of Atoms in Molecules (QTAIM) and stress tensor frameworks. Various theoretical analyses for these reactions, such as metallicity ξ(rb), ellipticity ε, total local energy density H(rb), stress tensor polarizability ℙσ, stress tensor eigenvalue λ3σ, and bond-path length, display differently for non-competitive and competitive reactions as well as for the conrotatory preferences either it is the transition state outward conrotatory (TSout) or transition state inward conrotatory (TSin) directions by presenting degeneracy or non-degeneracy in their results. The ellipticity profile provides the motion of the bond critical point locations due to the different substituents of cyclobutene. In agreement with experimental results, examinations demonstrated that N1 is a competitive reaction and N2-N3 are non-competitive reactions with TSout and TSin preference directions, respectively. The concordant results of QTAIM and stress tensor scalar and vectors with experimental results provide a better understanding of reaction mechanisms.
Collapse
Affiliation(s)
- Roya Momen
- State Key Laboratory of Powder Metallurgy, College of Chemistry and Chemical Engineering, Central South University Changsha, Changsha 410083, China
| | - Alireza Azizi
- College of Material Science and Engineering, Central South University, Changsha, Hunan 410083, China
| | | | - Mehdi Pazhoohesh
- School of Engineering and Sustainable Development, De Montfort University, The Gateway, Leicester LE1 9BH, United Kingdom
| | - Xiaobo Ji
- State Key Laboratory of Powder Metallurgy, College of Chemistry and Chemical Engineering, Central South University Changsha, Changsha 410083, China
| |
Collapse
|
19
|
Takács G, Sándor M, Szalai Z, Kiss R, Balogh GT. Analysis of the uncharted, druglike property space by self-organizing maps. Mol Divers 2021; 26:2427-2441. [PMID: 34709525 PMCID: PMC9532340 DOI: 10.1007/s11030-021-10343-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Accepted: 10/15/2021] [Indexed: 12/29/2022]
Abstract
Physicochemical properties are fundamental to predict the pharmacokinetic and pharmacodynamic behavior of drug candidates. Easily calculated descriptors such as molecular weight and logP have been found to correlate with the success rate of clinical trials. These properties have been previously shown to highlight a sweet-spot in the chemical space associated with favorable pharmacokinetics, which is superior against other regions during hit identification and optimization. In this study, we applied self-organizing maps (SOMs) trained on sixteen calculated properties of a subset of known drugs for the analysis of commercially available compound databases, as well as public biological and chemical databases frequently used for drug discovery. Interestingly, several regions of the property space have been identified that are highly overrepresented by commercially available chemical libraries, while we found almost completely unoccupied regions of the maps (commercially neglected chemical space resembling the properties of known drugs). Moreover, these underrepresented portions of the chemical space are compatible with most rigorous property filters applied by the pharma industry in medicinal chemistry optimization programs. Our results suggest that SOMs may be directly utilized in the strategy of library design for drug discovery to sample previously unexplored parts of the chemical space to aim at yet-undruggable targets.
Collapse
Affiliation(s)
- Gergely Takács
- Department of Chemical and Environmental Process Engineering, Budapest University of Technology and Economics, Műegyetem rakpart 3, Budapest, 1111, Hungary
- Mcule.com Kft, Bartók Béla út 105-113, Budapest, 1115, Hungary
| | - Márk Sándor
- Mcule.com Kft, Bartók Béla út 105-113, Budapest, 1115, Hungary
| | - Zoltán Szalai
- Mcule.com Kft, Bartók Béla út 105-113, Budapest, 1115, Hungary
| | - Róbert Kiss
- Mcule.com Kft, Bartók Béla út 105-113, Budapest, 1115, Hungary.
| | - György T Balogh
- Department of Chemical and Environmental Process Engineering, Budapest University of Technology and Economics, Műegyetem rakpart 3, Budapest, 1111, Hungary.
- Department of Pharmacodynamics and Biopharmacy, University of Szeged, Szeged, 6720, Hungary.
| |
Collapse
|
20
|
Vargas S, Hennefarth MR, Liu Z, Alexandrova AN. Machine Learning to Predict Diels-Alder Reaction Barriers from the Reactant State Electron Density. J Chem Theory Comput 2021; 17:6203-6213. [PMID: 34478623 DOI: 10.1021/acs.jctc.1c00623] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Reaction barriers are key to our understanding of chemical reactivity and catalysis. Certain reactions are so seminal in chemistry that countless variants, with or without catalysts, have been studied, and their barriers have been computed or measured experimentally. This wealth of data represents a perfect opportunity to leverage machine learning models, which could quickly predict barriers without explicit calculations or measurement. Here, we show that the topological descriptors of the quantum mechanical charge density in the reactant state constitute a set that is both rigorous and continuous and can be used effectively for the prediction of reaction barrier energies to a high degree of accuracy. We demonstrate this on the Diels-Alder reaction, highly important in biology and medicinal chemistry, and as such, studied extensively. This reaction exhibits a range of barriers as large as 270 kJ/mol. While we trained our single-objective supervised (labeled) regression algorithms on simpler Diels-Alder reactions in solution, they predict reaction barriers also in significantly more complicated contexts, such a Diels-Alder reaction catalyzed by an artificial enzyme and its evolved variants, in agreement with experimental changes in kcat. We expect this tool to apply broadly to a variety of reactions in solution or in the presence of a catalyst, for screening and circumventing heavily involved computations or experiments.
Collapse
Affiliation(s)
- Santiago Vargas
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1569, United States
| | - Matthew R Hennefarth
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1569, United States
| | - Zhihao Liu
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1569, United States
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1569, United States.,California NanoSystems Institute, University of California, Los Angeles, 570 Westwood Plaza, Los Angeles, California 90095-1569, United States
| |
Collapse
|
21
|
Bartashevich EV, Mukhitdinova SE, Tsirelson VG. Bond orders and electron delocalization indices for S–N, S–C and S–S bonds in 1,2,3-dithiazole systems. MENDELEEV COMMUNICATIONS 2021. [DOI: 10.1016/j.mencom.2021.09.029] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
|
22
|
Rajput A, Kumar M. Anti-Ebola: an initiative to predict Ebola virus inhibitors through machine learning. Mol Divers 2021; 26:1635-1644. [PMID: 34357513 PMCID: PMC8343361 DOI: 10.1007/s11030-021-10291-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2021] [Accepted: 07/28/2021] [Indexed: 01/17/2023]
Abstract
Ebola virus is a deadly pathogen responsible for a frequent series of outbreaks since 1976. Despite various efforts from researchers worldwide, its mortality and fatality are quite high. For antiviral drug discovery, the computational efforts are considered highly useful. Therefore, we have developed an 'anti-Ebola' web server, through quantitative structure-activity relationship information of available molecules with experimental anti-Ebola activities. Three hundred and five unique anti-Ebola compounds with their respective IC50 values were extracted from the 'DrugRepV' database. Later, the compounds were used to extract the molecular descriptors, which were subjected to regression-based model development. The robust machine learning techniques, namely support vector machine, random forest and artificial neural network, were employed using tenfold cross-validation. After a randomization approach, the best predictive model showed Pearson's correlation coefficient ranges from 0.83 to 0.98 on training/testing (T274) dataset. The robustness of the developed models was cross-evaluated using William's plot. The highly robust computational models are integrated into the web server. The 'anti-Ebola' web server is freely available at https://bioinfo.imtech.res.in/manojk/antiebola . We anticipate this will serve the scientific community for developing effective inhibitors against the Ebola virus.
Collapse
Affiliation(s)
- Akanksha Rajput
- Virology Unit and Bioinformatics Centre, Institute of Microbial Technology, Council of Scientific and Industrial Research (CSIR), Sector 39A, Chandigarh, 160036, India
| | - Manoj Kumar
- Virology Unit and Bioinformatics Centre, Institute of Microbial Technology, Council of Scientific and Industrial Research (CSIR), Sector 39A, Chandigarh, 160036, India. .,Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, 201002, India.
| |
Collapse
|
23
|
Verma P, Srivastava A, Srivastava K, Tandon P, Shimpi MR. Molecular Structure, Spectral Investigations, Hydrogen Bonding Interactions and Reactivity-Property Relationship of Caffeine-Citric Acid Cocrystal by Experimental and DFT Approach. Front Chem 2021; 9:708538. [PMID: 34381761 PMCID: PMC8350318 DOI: 10.3389/fchem.2021.708538] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2021] [Accepted: 06/21/2021] [Indexed: 11/27/2022] Open
Abstract
The pharmaceutical cocrystal of caffeine-citric acid (CAF-CA, Form II) has been studied to explore the presence of hydrogen bonding interactions and structure-reactivity-property relationship between the two constituents CAF and Citric acid. The cocrystal was prepared by slurry crystallization. Powder X-ray diffraction (PXRD) analysis was done to characterize CAF-CA cocrystal. Also, differential scanning calorimetry (DSC) confirmed the existence of CAF-CA cocrystal. The vibrational spectroscopic (FT-IR and FT-Raman) signatures and quantum chemical approach have been used as a strategy to get insights into structural and spectral features of CAF-CA cocrystal. There was a good correlation among the experimental and theoretical results of dimer of cocrystal, as this model is capable of covering all nearest possible interactions present in the crystal structure of cocrystal. The spectroscopic results confirmed that (O33-H34) mode forms an intramolecular (C25 = O28∙∙∙H34-O33), while (O26-H27) (O39-H40) and (O43-H44) groups form intermolecular hydrogen bonding (O26-H27∙∙∙N24-C22, O39-H40∙∙∙O52 = C51 and O43-H44∙∙∙O86 = C83) in cocrystal due to red shifting and increment in bond length. The quantum theory of atoms in molecules (QTAIM) analysis revealed (O88-H89∙∙∙O41) as strongest intermolecular hydrogen bonding interaction with interaction energy -12.4247 kcal mol-1 in CAF-CA cocrystal. The natural bond orbital analysis of the second-order theory of the Fock matrix highlighted the presence of strong interactions (N∙∙∙H and O∙∙∙H) in cocrystal. The HOMO-LUMO energy gap value shows that the CAF-CA cocrystal is more reactive, less stable and softer than CAF active pharmaceutical ingredients. The electrophilic and nucleophilic reactivities of atomic sites involved in intermolecular hydrogen bond interactions in cocrystal have been demonstrated by mapping electron density isosurfaces over electrostatic potential i.e. plotting molecular electrostatic potential (MESP) map. The molar refractivity value of cocrystal lies within the set range by Lipinski and hence it may be used as orally active form. The results show that the physicochemical properties of CAF-CA cocrystal are enhanced in comparison to CAF (API).
Collapse
Affiliation(s)
- Priya Verma
- Department of Physics, University of Lucknow, Lucknow, India
| | | | | | - Poonam Tandon
- Department of Physics, University of Lucknow, Lucknow, India
| | - Manishkumar R. Shimpi
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, Sweden
- Chemistry of Interfaces, Luleå University of Technology, Luleå, Sweden
| |
Collapse
|
24
|
Nandini Asha R, Ravindran Durai Nayagam B, Bhuvanesh N. Synthesis, molecular docking, and in silico ADMET studies of 4-benzyl-1-(2,4,6-trimethyl-benzyl)-piperidine: Potential Inhibitor of SARS-CoV2. Bioorg Chem 2021; 112:104967. [PMID: 33975232 PMCID: PMC8096530 DOI: 10.1016/j.bioorg.2021.104967] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2021] [Revised: 04/21/2021] [Accepted: 05/01/2021] [Indexed: 11/26/2022]
Abstract
Nowadays, over 200 countries face a wellbeing emergency because of epidemiological disease COVID-19 caused by the SARS-CoV-2 virus. It will cause a very high effect on the world's economy and the worldwide health sector. The present work is an investigation of the newly synthesized 4-benzyl-1-(2,4,6-trimethyl-benzyl)-piperidine (M1BZP) molecule's inhibitory potential against important protein targets of SARS-CoV-2 using computational approaches. M1BZP crystallizes in monoclinic type with P1211 space group. For the title compound M1BZP, spectroscopic characterization like 1H NMR, 13C NMR, FTIR, were carried out. The geometry of the compound had been optimized by the DFT method and its results were compared with the X-ray diffraction data. The calculated energies for the Highest Occupied Molecular Orbital (HOMO) and the Lowest Unoccupied Molecular Orbital (LUMO) showed the stability and reactivity of the title compound. Intermolecular interactions in the crystal network were determined using Hirshfeld surface analyses. The molecular electrostatic potential (MEP) picture was drawn using the same level of theory to visualize the chemical reactivity and charge distribution on the molecule. Molecular docking study performed for the synthesized compound revealed an efficient interaction with the COVID-19 protease and resulted in good activities. We hope the present study would help workers in the field to develop potential vaccines and therapeutics against the novel coronavirus. Virtual ADME studies were carried out as well and a relationship between biological, electronic, and physicochemical qualifications of the target compound was determined. Toxicity prediction by computational technique for the title compound was also carried out.
Collapse
Affiliation(s)
- R Nandini Asha
- Department of Chemistry and Research Centre, Pope's College (Autonomous), Sawyerpuram-628251, Affiliated to Manonmaniam Sundaranar University, Tirunelveli 627012, Tamil Nadu, India.
| | - B Ravindran Durai Nayagam
- Department of Chemistry and Research Centre, Pope's College (Autonomous), Sawyerpuram-628251, Affiliated to Manonmaniam Sundaranar University, Tirunelveli 627012, Tamil Nadu, India.
| | - Nattamai Bhuvanesh
- Department of Chemistry, Texas A&M University, College Station, TX 77842, USA.
| |
Collapse
|
25
|
A double bond with weak σ- and strong π-interactions is still a double bond. Nat Commun 2021; 12:4037. [PMID: 34188025 PMCID: PMC8241832 DOI: 10.1038/s41467-021-24238-x] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2020] [Accepted: 06/04/2021] [Indexed: 01/05/2023] Open
|
26
|
Shahamirian M, Azami SM. Strong intramolecular hydrogen bonding in confined amino acids. J Mol Graph Model 2021; 106:107913. [PMID: 33892298 DOI: 10.1016/j.jmgm.2021.107913] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 03/26/2021] [Accepted: 03/26/2021] [Indexed: 12/07/2022]
Abstract
Intramolecular hydrogen bonding is evaluated in three different amino acids encapsulated in C60 fullerene in the context of electron density analysis. While conventional intramolecular hydrogen bonding in isolated amino acids are dominated by electrostatic character, it is shown that strong intramolecular hydrogen bonding can be formed in confined amino acids so that in two cases covalent intramolecular hydrogen bonding is appeared in the confined species. Also, results show that zwitterionic amino acids are stable in confined state, where no implicit or explicit solvation is applied. Covalent character for intramolecular hydrogen bonding in amino acids have not yet been reported.
Collapse
Affiliation(s)
- M Shahamirian
- Department of Chemistry, Faculty of Science, Sarvestan Branch, Islamic Azad University, Sarvestan, 73451-173, Iran.
| | - S M Azami
- Department of Chemistry, Yasouj University, Yasouj, 75918-74934, Iran
| |
Collapse
|
27
|
Firme CL. Local potential energy density model (LPE): Applications and limitations to quantify intra/intermolecular interactions. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113143] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
28
|
Beck ME, Riplinger C, Neese F, Bistoni G. Unraveling individual host-guest interactions in molecular recognition from first principles quantum mechanics: Insights into the nature of nicotinic acetylcholine receptor agonist binding. J Comput Chem 2021; 42:293-302. [PMID: 33232540 DOI: 10.1002/jcc.26454] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2020] [Revised: 10/30/2020] [Accepted: 11/02/2020] [Indexed: 11/08/2022]
Abstract
Drug binding to a protein target is governed by a complex pattern of noncovalent interactions between the ligand and the residues in the protein's binding pocket. Here we introduce a generally applicable, parameter-free, computational method that allows for the identification, quantification, and analysis of the key ligand-residue interactions responsible for molecular recognition. Our strategy relies on Local Energy Decomposition analysis at the "gold-standard" coupled cluster DLPNO-CCSD(T) level. In the study case shown in this paper, nicotine and imidacloprid binding to the nicotinic acetylcholine receptor, our approach provides new insights into how individual amino acids in the active site determine sensitivity and selectivity of the ligands, extending and refining classical pharmacophore hypotheses. By inference, the method is applicable to any kind of host/guest interactions with potential applications in industrial biocatalysis and protein engineering.
Collapse
Affiliation(s)
- Michael Edmund Beck
- Department Computational Life Science, Bayer AG, Division Cropscience, Monheim am Rhein, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| |
Collapse
|
29
|
Kleemiss F, Wieduwilt EK, Hupf E, Shi MW, Stewart SG, Jayatilaka D, Turner MJ, Sugimoto K, Nishibori E, Schirmeister T, Schmidt TC, Engels B, Grabowsky S. Similarities and Differences between Crystal and Enzyme Environmental Effects on the Electron Density of Drug Molecules. Chemistry 2021; 27:3407-3419. [PMID: 33090581 PMCID: PMC7898524 DOI: 10.1002/chem.202003978] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Indexed: 01/28/2023]
Abstract
The crystal interaction density is generally assumed to be a suitable measure of the polarization of a low-molecular weight ligand inside an enzyme, but this approximation has seldomly been tested and has never been quantified before. In this study, we compare the crystal interaction density and the interaction electrostatic potential for a model compound of loxistatin acid (E64c) with those inside cathepsin B, in solution, and in vacuum. We apply QM/MM calculations and experimental quantum crystallography to show that the crystal interaction density is indeed very similar to the enzyme interaction density. Less than 0.1 e are shifted between these two environments in total. However, this difference has non-negligible consequences for derived properties.
Collapse
Affiliation(s)
- Florian Kleemiss
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Department of Chemistry and BiochemistryUniversity of BernFreiestrasse 3, 3012 BernSwitzerland
| | - Erna K. Wieduwilt
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019Université de Lorraine & CNRSBoulevard Arago, 57078 MetzFrance
| | - Emanuel Hupf
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
| | - Ming W. Shi
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Scott G. Stewart
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Dylan Jayatilaka
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Michael J. Turner
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Kunihisa Sugimoto
- Japan Synchrotron Radiation Research InstituteSPring-81-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5198Japan
- Institute for Integrated Cell-Materials SciencesKyoto UniversityYoshida-Ushinomiya-cho, Sakyo-ku, Kyoto 606-8501Japan
| | - Eiji Nishibori
- Division of Physics, Faculty of Pure and Applied Sciences, Tsukuba Research Center for Energy Materials ScienceUniversity of TsukubaTsukubaJapan
| | - Tanja Schirmeister
- Institute of Pharmaceutical and Biomedical SciencesJohannes-Gutenberg University MainzStaudingerweg 5, 55128 MainzGermany
| | - Thomas C. Schmidt
- Institute for Physical and Theoretical ChemistryJulius-Maximilians-University WürzburgEmil-Fischer-Str. 42, 97074 WürzburgGermany
| | - Bernd Engels
- Institute for Physical and Theoretical ChemistryJulius-Maximilians-University WürzburgEmil-Fischer-Str. 42, 97074 WürzburgGermany
| | - Simon Grabowsky
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Department of Chemistry and BiochemistryUniversity of BernFreiestrasse 3, 3012 BernSwitzerland
| |
Collapse
|
30
|
Wang H, Qin Z, Yan A. Classification models and SAR analysis on CysLT1 receptor antagonists using machine learning algorithms. Mol Divers 2021; 25:1597-1616. [PMID: 33534023 DOI: 10.1007/s11030-020-10165-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2020] [Accepted: 11/27/2020] [Indexed: 12/21/2022]
Abstract
Cysteinyl leukotrienes 1 (CysLT1) receptor is a promising drug target for rhinitis or other allergic diseases. In our study, we built classification models to predict bioactivities of CysLT1 receptor antagonists. We built a dataset with 503 CysLT1 receptor antagonists which were divided into two groups: highly active molecules (IC50 < 1000 nM) and weakly active molecules (IC50 ≥ 1000 nM). The molecules were characterized by several descriptors including CORINA descriptors, MACCS fingerprints, Morgan fingerprint and molecular SMILES. For CORINA descriptors and two types of fingerprints, we used the random forests (RF) and deep neural networks (DNN) to build models. For molecular SMILES, we used recurrent neural networks (RNN) with the self-attention to build models. The accuracies of test sets for all models reached 85%, and the accuracy of the best model (Model 2C) was 93%. In addition, we made structure-activity relationship (SAR) analyses on CysLT1 receptor antagonists, which were based on the output from the random forest models and RNN model. It was found that highly active antagonists usually contained the common substructures such as tetrazoles, indoles and quinolines. These substructures may improve the bioactivity of the CysLT1 receptor antagonists.
Collapse
Affiliation(s)
- Hongzhao Wang
- State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, University of Chemical Technology, Beijing, People's Republic of China
| | - Zijian Qin
- State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, University of Chemical Technology, Beijing, People's Republic of China
| | - Aixia Yan
- State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, University of Chemical Technology, Beijing, People's Republic of China.
| |
Collapse
|
31
|
Galindo-Murillo R, Winkler L, García-Ramos JC, Ruiz-Azuara L, Cortés-Guzmán F, Cheatham TE. Ancillary Ligand in Ternary Cu II Complexes Guides Binding Selectivity toward Minor-Groove DNA. J Phys Chem B 2020; 124:11648-11658. [PMID: 33320672 PMCID: PMC7770820 DOI: 10.1021/acs.jpcb.0c09296] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Copper-containing
compounds known as Casiopeínas are
biologically active molecules which show promising antineoplastic
effects against several cancer types. Two possible hypotheses regarding
the mode of action of the Casiopeínas have emerged from
the experimental evidence: the generation of reactive oxygen species
or the ability of the compounds to bind and interact with nucleic
acids. Using robust molecular dynamics simulations, we investigate
the interaction of four different Casiopeínas with the DNA duplex d(GCACGAACGAACGAACGC). The studied copper complexes contain either 4–7- or 5–6-substituted
dimethyl phenanthroline as the primary ligand and either glycinate
or acetylacetonate as the secondary ligand. For statistical significance
and to reduce bias in the simulations, four molecules of each copper
compound were manually placed at a distance of 10 Å away from
the DNA and 20 independent molecular dynamics simulations were performed,
each reaching at least 30 μs. This time scale allows us to reproduce
expected DNA terminal base-pair fraying and also to observe intercalation/base-pair
eversion events generated by the compounds interacting with DNA. The
results reveal that the secondary ligand is the guide toward the mode
of binding between the copper complex and DNA in which glycinate prefers
minor-groove binding and acetylacetonate produces base-pair eversion
and intercalation. The CuII complexes containing glycinate
interact within the DNA minor groove which are stabilized principally
by the hydrogen bonds formed between the amino group of the aminoacidate
moiety, whereas the compounds with the acetylacetonate do not present
a stable network of hydrogen bonds and the ligand interactions enhance
DNA breathing dynamics that result in base-pair eversion.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, Utah 84112, United States
| | - Lauren Winkler
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, Utah 84112, United States
| | - Juan Carlos García-Ramos
- Escuela de Ciencias de la Salud, Universidad Autónoma de Baja California, Blvd. Zertuche y Blvd. Los Lagos, Fraccionamiento Valle Dorado, Ensenada, Baja California 22890, Mexico
| | - Lena Ruiz-Azuara
- Departamento de Química Inorgánica y Nuclear. Facultad de Química. Universidad Nacional Autónoma de México. Avenida Universidad 3000, Ciudad Universitaria, Ciudad de México 04510, Mexico
| | - Fernando Cortés-Guzmán
- Departamento de Fisicoquímica. Instituto de Química. Universidad Nacional Autónoma de México. Avenida Universidad 3000, Ciudad Universitaria, Ciudad de México 04510, Mexico
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, Utah 84112, United States
| |
Collapse
|
32
|
Brovarets' OO, Muradova A, Hovorun DM. A Quantum-Mechanical Looking Behind the Scene of the Classic G·C Nucleobase Pairs Tautomerization. Front Chem 2020; 8:574454. [PMID: 33330362 PMCID: PMC7732530 DOI: 10.3389/fchem.2020.574454] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Accepted: 09/02/2020] [Indexed: 11/13/2022] Open
Abstract
For the first time, at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of theory, a comprehensive quantum-mechanical investigation of the physico-chemical mechanism of the tautomeric wobblization of the four biologically-important G·C nucleobase pairs by the participation of the monomers in rare, in particular mutagenic, tautomeric forms (marked with an asterisk) was provided. These novel tautomeric transformations (wobblization or shifting of the bases within the pair) are intrinsically inherent properties of the G·C nucleobase pairs. In this study, we have obtained intriguing results, lying far beyond the existing representations. Thus, it was shown that Löwdin's G*·C*(WC) base pair does not tautomerize according to the wobblization mechanism. Tautomeric wobblization of the G*·C*(rWC) (relative Gibbs free energy ΔG = 0.00/relative electronic energy ΔE = 0.00 kcal·mol-1) ("r"-means the configuration of the base pair in reverse position; "WC"-the classic Watson-Crick configuration) and G*t·C*(H) (ΔG = -0.19/ΔE = 0.29 kcal·mol-1) ("H"-Hoogsteen configuration;"t" denotes the O6H hydroxyl group in the trans position) base pairs are preceded by the stages of the base pairs tautomerization by the single proton transfer (SPT). It was established that the G*t·C*(rH) (ΔG = 2.21/ΔE = 2.81 kcal·mol-1) base pair can be wobbled through two different pathways via the traditional one-stage mechanism through the TSs, which are tight G+·C- ion pairs, stabilized by the participation of only two intermolecular H-bonds. It was found out that the G·C base pair is most likely incorporated into the DNA/RNA double helix with parallel strands in the G*·C*(rWC), G·C*(rwwc), and G*·C(rwwc) ("w"-wobble configuration of the pair) tautomeric forms, which are in rapid tautomeric equilibrium with each other. It was proven that the G*·C*(rWC) nucleobase pair is also in rapid tautomeric equilibrium with the eight tautomeric forms of the so-called Levitt base pair. It was revealed that a few cases of tautomerization via the DPT of the nucleobase pairs by the participation of the C8H group of the guanine had occurred. The biological role of the obtained results was also made apparent.
Collapse
Affiliation(s)
- Ol'ha O. Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Alona Muradova
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| |
Collapse
|
33
|
Zhang G, Chen D, Wang S, Chen H, Wei N, Chen G. Molecular Insight into the Discrepancy of Antitubercular Activity between 8‐Nitro and 8‐Cyano Benzothiazinones. ChemistrySelect 2020. [DOI: 10.1002/slct.202004156] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Gang Zhang
- State Key Laboratory of Bioactive Substances and Function of Natural Medicine Institute of Materia Medica, Peking Union Medical College and Chinese Academy of Medical Sciences Beijing 100050 China
| | - Daoxing Chen
- School of Pharmaceutical Sciences Wenzhou Medical University, Wenzhou Zhejiang 325035 China
| | - Sufang Wang
- The Environment Analysis and Testing Laboratory Chinese Research Academy of Environmental Sciences Beijing 100012 China
| | - Hualong Chen
- Beijing Institute of Brain Disorders, Laboratory of Brain Disorders, Ministry of Science and Technology Planning and Research Department, Beijing Science Center Collaborative Innovation Center for Brain Disorders, Capital Medical University Beijing 100029 China
| | - Ning Wei
- Division of Hematology-Oncology Cancer Therapeutics Program, UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh Department of Medicine, University of Pittsburgh School of Medicine Pittsburgh PA 15261 USA
| | - Guirong Chen
- Liaoning University of Traditional Chinese Medicine, Shenyang Liaoning 110847 China
| |
Collapse
|
34
|
Sadjadi S, Matta CF, Hamilton IP. Bonding and metastability for Group 12 dications. J Comput Chem 2020; 42:40-49. [PMID: 33063900 DOI: 10.1002/jcc.26431] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2020] [Revised: 09/25/2020] [Accepted: 09/28/2020] [Indexed: 11/06/2022]
Abstract
Electronic structure and bonding properties of the Group 12 dications M2 2+ (M = Zn, Cd, Hg) are investigated and electron density-derived quantities are used to characterize the metastability of these species. Of particular interest are the complementary descriptions afforded by the Laplacian of the electron density ∇2 ρ(r) and the one-electron Bohm quantum potential (Q = ∇ 2 ρ r / 2 ρ r ) along the bond path. Further, properties derived from the pair density including the localization-delocalization matrices (LDMs) and the interacting quantum atoms (IQA) energies are analyzed within the framework of the quantum theory of atoms in molecules (QTAIM). From the crossing points of the singlet (ground) and triplet (excited) potential energy curves, the barriers for dissociation (BFD) are estimated to be 25.2 kcal/mol (1.09 eV) for Zn2 2+ , 22.8 kcal/mol (0.99 eV) for Cd2 2+ , and 26.4 kcal/mol (1.14 eV) for Hg2 2+ . For comparison and benchmarking purposes, the case of N2 2+ is considered as a texbook example of metastability. At the equilibrium geometries, LDMs, which are used here as an electronic fingerprinting tool, discriminate and group together Group 12 M2 2+ from its isoelectronic Group 11 M2 . While "classical" bonding indices are inconclusive in establishing regions of metastability in the bonding, it is shown that the one-electron Bohm quantum potential is promising in this regard.
Collapse
Affiliation(s)
- SeyedAbdolreza Sadjadi
- Department of Physics, Faculty of Science, Laboratory for Space Research, The University of Hong Kong, Hong Kong SAR, China
| | - Chérif F Matta
- Department of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia, Canada.,Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Ian P Hamilton
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario, Canada
| |
Collapse
|
35
|
Characteristic Substitution Shift Model (CSSM) in the CCFDF (Charge-Charge Flux-Dipole Flux) model for the dipole moment derivatives of molecules X2CY (X = F, cl; Y = O, S). Chem Phys 2020. [DOI: 10.1016/j.chemphys.2020.110830] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
36
|
Abstract
Background: A principle of complementarity is a well-established concept in chemistry and biology. This concept is based on the overlap of electron clouds of the molecules in question. Materials & methods: In this article, one such approach (an in-house developed quantum free-orbital AlteQ method) was used to evaluate the complementarity of 51 CDK-ligand complexes. Results: A significant universally applicable correlation (adjusted R2 = 0.9749; p < 2.2 × 10-16) relating the product of ligand and enzyme electron densities to the product of distances between the contacting atomic centers and the type of atoms involved in the interaction was found. Conclusion: The terms calculated in this article can provide a good basis for prognosis of bioactivity and scientifically based molecular docking.
Collapse
|
37
|
Sowlati-Hashjin S, Karttunen M, Matta CF. Manipulation of Diatomic Molecules with Oriented External Electric Fields: Linear Correlations in Atomic Properties Lead to Nonlinear Molecular Responses. J Phys Chem A 2020; 124:4720-4731. [PMID: 32337997 DOI: 10.1021/acs.jpca.0c02569] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Oriented external electric fields (OEEFs) have been shown to have great potential in being able to provide unprecedented control of chemical reactions, catalysis, and selectivity with applications ranging from H2 storage to molecular machines. We report a theoretical study of the atomic origins of molecular changes because of OEEFs since understanding the characteristics of OEEF-induced couplings between atomic and molecular properties is an important step toward comprehensive understanding of the effects of strong external fields on the molecular structure, stability, and reactivity. We focus on the atomic and molecular (bond) properties of a set of homo- (H2, N2, O2, F2, and Cl2) and heterodiatomic (HF, HCl, CO, and NO) molecules under intense external electric fields in the context of quantum theory of atoms in molecules (QTAIM). It is shown that the atomic properties (atomic charges, energies, and localization indices) correlate linearly with the field strengths, but molecular properties (bond length, electron density at the bond critical point, and electron delocalization index) exhibit nonlinear responses to the imposed fields. In particular, the changes in the electron density distribution alter the shapes and locations of the zero-flux surfaces, atomic volumes, atomic electron population, and localization/delocalization indices. The topography and topology of the molecular electrostatic potential undergo dramatic changes. External fields also perturb the covalent-polar-ionic characteristic of the studied chemical bonds, hallmarking the impact of electric fields on the stability and reactivity of chemical compounds. The findings are well-rationalized within the framework of the QTAIM and form a coherent conceptual understanding of these effects in prototypical diatomic molecules.
Collapse
Affiliation(s)
- Shahin Sowlati-Hashjin
- Department of Chemistry, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,The Centre of Advanced Materials and Biomaterials Research, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,Department of Chemistry, Saint Mary's University, Halifax, Nova Scotia B3H 3C3, Canada.,Department of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia B3M 2J6, Canada
| | - Mikko Karttunen
- Department of Chemistry, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,The Centre of Advanced Materials and Biomaterials Research, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,Department of Applied Mathematics, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada
| | - Chérif F Matta
- Department of Chemistry, Saint Mary's University, Halifax, Nova Scotia B3H 3C3, Canada.,Department of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia B3M 2J6, Canada.,Department of Chemistry, Dalhousie University, Halifax, Nova Scotia B3H,4J3, Canada
| |
Collapse
|
38
|
Drug design by machine-trained elastic networks: predicting Ser/Thr-protein kinase inhibitors' activities. Mol Divers 2020; 25:899-909. [PMID: 32222890 DOI: 10.1007/s11030-020-10074-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2020] [Accepted: 03/11/2020] [Indexed: 12/23/2022]
Abstract
An elastic network model (ENM) represents a molecule as a matrix of pairwise atomic interactions. Rich in coded information, ENMs are hereby proposed as a novel tool for the prediction of the activity of series of molecules, with widely different chemical structures, but a common biological activity. The new approach is developed and tested using a set of 183 inhibitors of serine/threonine-protein kinase enzyme (Plk3) which is an enzyme implicated in the regulation of cell cycle and tumorigenesis. The elastic network (EN) predictive model is found to exhibit high accuracy and speed compared to descriptor-based machine-trained modeling. EN modeling appears to be a highly promising new tool for the high demands of industrial applications such as drug and material design.
Collapse
|
39
|
Ramírez-Palma LG, García-Jacas CR, García-Ramos JC, Almada-Monter R, Galindo-Murillo R, Cortés-Guzmán F. Pharmacophoric sites of anticancer metal complexes located using quantum topological atomic descriptors. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2019.127480] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
|
40
|
Martins FA, Freitas MP. Bioconformational modulation of a thymidine kinase enzyme ligand through F⋯HO intramolecular hydrogen bond. J Mol Graph Model 2020; 97:107545. [PMID: 32062583 DOI: 10.1016/j.jmgm.2020.107545] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2019] [Revised: 01/20/2020] [Accepted: 01/21/2020] [Indexed: 10/25/2022]
Abstract
While the induced-fit of a ligand towards an enzyme is pivotally dictated by intermolecular hydrogen bonds between the small molecule and amino acid residues in the binding site, the role of intramolecular hydrogen bond as contributing interaction for a bioactive conformation is not well understood. This work reports a theoretical conformational analysis of a thymidine kinase enzyme ligand (NMF) that is prone to experience an F⋯HO intramolecular hydrogen bond, inside and outside the biological binding site. This interaction stabilizes the most favorable conformations of NMF in the gas phase and, although it can be disrupted in a biological environment due to intermolecular hydrogen bonds in some cases, these interactions are competitive in other systems. Therefore, an intramolecular hydrogen bond can affect the conformational likeliness most related to the bioactivity. Moreover, isolated conformations governed by this interaction cannot be unequivocally used to generate molecular descriptors in 3D-QSAR (Quantitative Structure-Activity Relationships), since the bioactive conformation may not be determined only by intramolecular interactions.
Collapse
Affiliation(s)
- Francisco A Martins
- Department of Chemistry, Federal University of Lavras, 37200-900, Lavras, MG, Brazil
| | - Matheus P Freitas
- Department of Chemistry, Federal University of Lavras, 37200-900, Lavras, MG, Brazil.
| |
Collapse
|
41
|
Verma PL, Gejji SP. Electronic structure, vibrational spectra and 1H NMR chemical shifts of the ion pair composites within imidazolium functionalized geminal dicationic ionic liquids from density functional theory. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2019.127112] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
42
|
Bartashevich EV, Matveychuk YV, Mukhitdinova SE, Sobalev SA, Khrenova MG, Tsirelson VG. The common trends for the halogen, chalcogen, and pnictogen bonds via sorting principles and local bonding properties. Theor Chem Acc 2020. [DOI: 10.1007/s00214-019-2534-y] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
43
|
Ramírez-Palma DI, García-Jacas CR, Carpio-Martínez P, Cortés-Guzmán F. Predicting reactive sites with quantum chemical topology: carbonyl additions in multicomponent reactions. Phys Chem Chem Phys 2020; 22:9283-9289. [DOI: 10.1039/d0cp00300j] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The reactivity of an atom within a molecule depends mostly on the way the electron density polarizes reflected in the quadrupole moment of the reactive atom.
Collapse
Affiliation(s)
| | - Cesar R. García-Jacas
- Cátedras CONACYT – Departamento de Ciencias de la Computación
- Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE)
- Ensenada
- Mexico
| | | | | |
Collapse
|
44
|
Brovarets’ OO, Hovorun DM. A new era of the prototropic tautomerism of the quercetin molecule: A QM/QTAIM computational advances. J Biomol Struct Dyn 2019; 38:4774-4800. [PMID: 31711364 DOI: 10.1080/07391102.2019.1691660] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Ol’ha O. Brovarets’
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| |
Collapse
|
45
|
Protein polarization effects in the thermodynamic computation of vibrational Stark shifts. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2522-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
46
|
Brovarets’ OO, Hovorun DM. Intramolecular tautomerization of the quercetin molecule due to the proton transfer: QM computational study. PLoS One 2019; 14:e0224762. [PMID: 31751372 PMCID: PMC6874073 DOI: 10.1371/journal.pone.0224762] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2019] [Accepted: 10/20/2019] [Indexed: 12/21/2022] Open
Abstract
Quercetin molecule (3, 3', 4', 5, 7-pentahydroxyflavone, C15H10O7) is an important flavonoid compound of natural origin, consisting of two aromatic A and B rings linked through the C ring with endocyclic oxygen atom and five hydroxyl groups attached to the 3, 3', 4', 5 and 7 positions. This molecule is found in many foods and plants, and is known to have a wide range of therapeutic properties, like an anti-oxidant, anti-toxic, anti-inflammatory etc. In this study for the first time we have revealed and investigated the pathways of the tautomeric transformations for the most stable conformers of the isolated quercetin molecule (Brovarets' & Hovorun, 2019) via the intramolecular proton transfer. Energetic, structural, dynamical and polar characteristics of these transitions, in particular relative Gibbs free and electronic energies, characteristics of the intramolecular specific interactions-H-bonds and attractive van der Waals contacts, have been analysed in details. It was demonstrated that the most probable process among all investigated is the proton transfer from the O3H hydroxyl group of the C ring to the C2' carbon atom of the C2'H group of the B ring along the intramolecular O3H…C2' H-bond with the further formation of the C2'H2 group. It was established that the proton transfer from the hydroxyl groups to the carbon atoms of the neighboring CH groups is assisted at the transition states by the strong intramolecular HCH…O H-bond (~28.5 kcal∙mol-1). The least probable path of the proton transfer-from the C8H group to the endocyclic O1 oxygen atom-causes the decyclization of the C ring in some cases. It is shortly discussed the biological importance of the obtained results.
Collapse
Affiliation(s)
- Ol’ha O. Brovarets’
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| |
Collapse
|
47
|
The kernel energy method applied to quantum theory of atoms in molecules – energies of interacting quantum atoms. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.136650] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
48
|
Electronic structure and spectral characteristics of alkyl substituted imidazolium based dication-X2 (X = Br, BF4, PF6 and CF3SO3) complexes from theory. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111548] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
49
|
Brovarets' OO, Oliynyk TA, Hovorun DM. Novel Tautomerisation Mechanisms of the Biologically Important Conformers of the Reverse Löwdin, Hoogsteen, and Reverse Hoogsteen G *·C * DNA Base Pairs via Proton Transfer: A Quantum-Mechanical Survey. Front Chem 2019; 7:597. [PMID: 31620420 PMCID: PMC6759773 DOI: 10.3389/fchem.2019.00597] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Accepted: 08/12/2019] [Indexed: 11/22/2022] Open
Abstract
For the first time, in this study with the use of QM/QTAIM methods we have exhaustively investigated the tautomerization of the biologically-important conformers of the G*·C* DNA base pair-reverse Löwdin G*·C*(rWC), Hoogsteen G*'·C*(H), and reverse Hoogsteen G*'·C*(rH) DNA base pairs-via the single (SPT) or double (DPT) proton transfer along the neighboring intermolecular H-bonds. These tautomeric reactions finally lead to the formation of the novel G· C O 2 * (rWC), G N 2 * · C(rWC), G*'N2·C(rWC), G N 7 * · C(H), and G*'N7·C(rH) DNA base mispairs. Gibbs free energies of activation for these reactions are within the range 3.64-31.65 kcal·mol-1 in vacuum under normal conditions. All TSs are planar structures (Cs symmetry) with a single exception-the essentially non-planar transition state TSG*·C*(rWC)↔G+·C-(rWC) (C1 symmetry). Analysis of the kinetic parameters of the considered tautomerization reactions indicates that in reality only the reverse Hoogsteen G*'·C*(rH) base pair undergoes tautomerization. However, the population of its tautomerised state G*'N7·C(rH) amounts to an insignificant value-2.3·10-17. So, the G*·C*(rWC), G*'·C*(H), and G*'·C*(rH) base pairs possess a permanent tautomeric status, which does not depend on proton mobility along the neighboring H-bonds. The investigated tautomerization processes were analyzed in details by applying the author's unique methodology-sweeps of the main physical and chemical parameters along the intrinsic reaction coordinate (IRC). In general, the obtained data demonstrate the tautomeric mobility and diversity of the G*·C* DNA base pair.
Collapse
Affiliation(s)
- Ol'ha O. Brovarets'
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Pharmacology, Bohomolets National Medical University, Kyiv, Ukraine
| | - Timothy A. Oliynyk
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
- Department of Pathophysiology, Bohomolets National Medical University, Kyiv, Ukraine
| |
Collapse
|
50
|
Brovarets’ OO, Hovorun DM. Conformational diversity of the quercetin molecule: a quantum-chemical view. J Biomol Struct Dyn 2019; 38:2817-2836. [DOI: 10.1080/07391102.2019.1656671] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Ol’ha O. Brovarets’
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
| | - Dmytro M. Hovorun
- Department of Molecular and Quantum Biophysics, Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kyiv, Ukraine
- Department of Molecular Biotechnology and Bioinformatics, Institute of High Technologies, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine
| |
Collapse
|