1
|
Gómez-Tatay L, Hernández-Andreu JM. Xenobiology for the Biocontainment of Synthetic Organisms: Opportunities and Challenges. Life (Basel) 2024; 14:996. [PMID: 39202738 PMCID: PMC11355180 DOI: 10.3390/life14080996] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2024] [Revised: 08/07/2024] [Accepted: 08/09/2024] [Indexed: 09/03/2024] Open
Abstract
Since the development of recombinant DNA technologies, the need to establish biosafety and biosecurity measures to control genetically modified organisms has been clear. Auxotrophies, or conditional suicide switches, have been used as firewalls to avoid horizontal or vertical gene transfer, but their efficacy has important limitations. The use of xenobiological systems has been proposed as the ultimate biosafety tool to circumvent biosafety problems in genetically modified organisms. Xenobiology is a subfield of Synthetic Biology that aims to construct orthogonal biological systems based on alternative biochemistries. Establishing true orthogonality in cell-based or cell-free systems promises to improve and assure that we can progress in synthetic biology safely. Although a wide array of strategies for orthogonal genetic systems have been tested, the construction of a host harboring fully orthogonal genetic system, with all parts operating in an orchestrated, integrated, and controlled manner, still poses an extraordinary challenge for researchers. In this study, we have performed a thorough review of the current literature to present the main advances in the use of xenobiology as a strategy for biocontainment, expanding on the opportunities and challenges of this field of research.
Collapse
Affiliation(s)
- Lucía Gómez-Tatay
- Institute of Life Sciences, Universidad Católica de Valencia San Vicente Mártir, 46001 Valencia, Spain;
| | - José Miguel Hernández-Andreu
- Grupo de Investigación en Medicina Molecular y Mitocondrial, Universidad Católica de Valencia San Vicente Mártir, 46001 Valencia, Spain
| |
Collapse
|
2
|
Szczygiel M, Derewenda U, Scheiner S, Minor W, Derewenda ZS. A structural role for tryptophan in proteins, and the ubiquitous Trp C δ1-H...O=C (backbone) hydrogen bond. Acta Crystallogr D Struct Biol 2024; 80:551-562. [PMID: 38941144 PMCID: PMC11220837 DOI: 10.1107/s2059798324005515] [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: 03/12/2024] [Accepted: 06/09/2024] [Indexed: 06/29/2024] Open
Abstract
Tryptophan is the most prominent amino acid found in proteins, with multiple functional roles. Its side chain is made up of the hydrophobic indole moiety, with two groups that act as donors in hydrogen bonds: the Nϵ-H group, which is a potent donor in canonical hydrogen bonds, and a polarized Cδ1-H group, which is capable of forming weaker, noncanonical hydrogen bonds. Due to adjacent electron-withdrawing moieties, C-H...O hydrogen bonds are ubiquitous in macromolecules, albeit contingent on the polarization of the donor C-H group. Consequently, Cα-H groups (adjacent to the carbonyl and amino groups of flanking peptide bonds), as well as the Cϵ1-H and Cδ2-H groups of histidines (adjacent to imidazole N atoms), are known to serve as donors in hydrogen bonds, for example stabilizing parallel and antiparallel β-sheets. However, the nature and the functional role of interactions involving the Cδ1-H group of the indole ring of tryptophan are not well characterized. Here, data mining of high-resolution (r ≤ 1.5 Å) crystal structures from the Protein Data Bank was performed and ubiquitous close contacts between the Cδ1-H groups of tryptophan and a range of electronegative acceptors were identified, specifically main-chain carbonyl O atoms immediately upstream and downstream in the polypeptide chain. The stereochemical analysis shows that most of the interactions bear all of the hallmarks of proper hydrogen bonds. At the same time, their cohesive nature is confirmed by quantum-chemical calculations, which reveal interaction energies of 1.5-3.0 kcal mol-1, depending on the specific stereochemistry.
Collapse
Affiliation(s)
- Michal Szczygiel
- Department of Molecular Physiology and Biological PhysicsUniversity of Virginia1340 Jefferson Park AvenueCharlottesvilleVA22908-0736USA
| | - Urszula Derewenda
- Department of Molecular Physiology and Biological PhysicsUniversity of Virginia1340 Jefferson Park AvenueCharlottesvilleVA22908-0736USA
| | - Steve Scheiner
- Department of Chemistry and BiochemistryUtah State UniversityLoganUtahUSA
| | - Wladek Minor
- Department of Molecular Physiology and Biological PhysicsUniversity of Virginia1340 Jefferson Park AvenueCharlottesvilleVA22908-0736USA
| | - Zygmunt S. Derewenda
- Department of Molecular Physiology and Biological PhysicsUniversity of Virginia1340 Jefferson Park AvenueCharlottesvilleVA22908-0736USA
| |
Collapse
|
3
|
Peluzo BMTC, Moura RT, Kraka E. Extraction of uranyl from spent nuclear fuel wastewater via complexation-a local vibrational mode study. J Mol Model 2024; 30:216. [PMID: 38888814 DOI: 10.1007/s00894-024-06000-4] [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: 03/14/2024] [Accepted: 05/30/2024] [Indexed: 06/20/2024]
Abstract
CONTEXT The efficient extraction of uranyl from spent nuclear fuel wastewater for subsequent reprocessing and reuse is an essential effort toward minimization of long-lived radioactive waste. N-substituted amides and Schiff base ligands are propitious candidates, where extraction occurs via complexation with the uranyl moiety. In this study, we extensively probed chemical bonding in various uranyl complexes, utilizing the local vibrational modes theory alongside QTAIM and NBO analyses. We focused on (i) the assessment of the equatorial O-U and N-U bonding, including the question of chelation, and (ii) how the strength of the axial U = O bonds of the uranyl moiety changes upon complexation. Our results reveal that the strength of the equatorial uranium-ligand interactions correlates with their covalent character and with charge donation from O and N lone pairs into the vacant uranium orbitals. We also found an inverse relationship between the covalent character of the equatorial ligand bonds and the strength of the axial uranium-oxygen bond. In summary, our study provides valuable data for a strategic modulation of N-substituted amide and Schiff base ligands towards the maximization of uranyl extraction. METHOD Quantum chemistry calculations were performed under the PBE0 level of theory, paired with the relativistic NESCau Hamiltonian, currently implemented in Cologne2020 (interfaced with Gaussian16). Wave functions were expanded in the cc-pwCVTZ-X2C basis set for uranium and Dunning's cc-pVTZ for the remaining atoms. For the bonding properties, we utilized the package LModeA in the local modes analyses, AIMALL in the QTAIM calculations, and NBO 7.0 for the NBO analyses.
Collapse
Affiliation(s)
- Bárbara M T C Peluzo
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX, 75275-0314, USA
| | - Renaldo T Moura
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX, 75275-0314, USA
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraíba, Areia, 58397-000, Paraíba, Brazil
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX, 75275-0314, USA.
| |
Collapse
|
4
|
Venkataramanan NS, Suvitha A, Sahara R. Unveiling the Intermolecular Interactions between Drug 5-Fluorouracil and Watson-Crick/Hoogsteen Base Pairs: A Computational Analysis. ACS OMEGA 2024; 9:24831-24844. [PMID: 38882136 PMCID: PMC11170692 DOI: 10.1021/acsomega.4c01545] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/18/2024] [Revised: 05/09/2024] [Accepted: 05/14/2024] [Indexed: 06/18/2024]
Abstract
The adsorption of 5-fluorouracil (5FU) on Watson-Crick (WC) base pairs and Hoogsteen (HT) base pairs has been studied using the dispersion-corrected density functional theory (DFT). The adsorption, binding energy, and thermochemistry for the drug 5FU on the WC and HT base pairs were determined. The most stable geometries were near planar geometry, and 5FU has a higher preference for WC than HT base pairs. The adsorption energies of 5FU on nucleobase pairs are consistently higher than pristine nucleobase pairs, indicating that nucleobase pair cleavage is less likely during the adsorption of the 5FU drug. The enthalpy change for the formation of 5FU-DNA base pairs is higher than that for the formation of 5FU-nucleobases and is enthalpy-driven. The E gap of AT base pairs is higher, suggesting that their chemical reactivity toward further reaction would be less than that of GC base pairs. The electron density difference (EDD) analysis shows a significant decrease in electron density in aromatic regions on the purine bases (adenine/guanine) compared to the pyrimidine bases. The MESP diagram of the stable 5FU-nucleobase pair complexes shows a directional interaction, with the positive regions in a molecule interacting with the negative region of other molecules. The atoms in molecule analysis show that the ρ(r) values of C=O···H-N are higher than those of N···H/N-H···O. The N···H intermolecular bonds between the base pair/drug and nucleobases are weak, closed shell interactions and are electrostatic in nature. The noncovalent interaction analysis shows that several new spikes are engendered along with an increase in their strength, which indicates that the H-bonding interactions are stronger and play a dominant role in stabilizing the complexes. Energy decomposition analysis shows that the drug-nucleobase pair complex has a marginal increase in the electrostatic contributions compared to nucleobase pair complexes.
Collapse
Affiliation(s)
| | | | - Ryoji Sahara
- Research Center for Structural Materials, National Institute for Materials Science (NIMS), 1-2-1 Sengen, Tsukuba 305-0047, Japan
| |
Collapse
|
5
|
Mak CH. Hydration Waters Make Up for the Missing Third Hydrogen Bond in the A·T Base Pair. ACS PHYSICAL CHEMISTRY AU 2024; 4:180-190. [PMID: 38560756 PMCID: PMC10979491 DOI: 10.1021/acsphyschemau.3c00058] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 01/05/2024] [Accepted: 01/08/2024] [Indexed: 04/04/2024]
Abstract
Base pairing complementarity is central to DNA function. G·C and A·T pair specificity is thought to originate from the different number of hydrogen bonds the pairs make. Quantifying how many hydrogen bonds exist can be difficult because water molecules in the surrounding can make up for or disrupt direct hydrogen bonds, and the hydration structures around A·T and G·C pairs on duplex DNA are distinct. Large-scale computer simulations have been used here to create a detailed map for the hydration structure on A·T and G·C base pairs in water. The contributions of specific hydration waters to the free energy of each of the hydrogen bonds in the A·T and G·C pairs were computed. Using the equilibrium fractions of hydrated versus unhydrated states from the hydration profiles, the impact of specific bound waters on each hydrogen bond can be uniquely quantified using a thermodynamic construction. The findings suggest that hydration water in the minor groove of an A·T pair can provide up to about 2 kcal/mol of free energy advantage, effectively making up for the missing third hydrogen bond in the A·T pair compared to G·C, rendering the intrinsic thermodynamic stability of the A·T pair almost synonymous with G·C.
Collapse
Affiliation(s)
- Chi H. Mak
- Departments of Chemistry
and Quantitative and Computational Biology, and Center of Applied
Mathematical Sciences, University of Southern
California, Los Angeles, California 90089, United States
| |
Collapse
|
6
|
Dangat Y, Freindorf M, Kraka E. Mechanistic Insights into S-Depalmitolyse Activity of Cln5 Protein Linked to Neurodegeneration and Batten Disease: A QM/MM Study. J Am Chem Soc 2024; 146:145-158. [PMID: 38055807 DOI: 10.1021/jacs.3c06397] [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: 12/08/2023]
Abstract
Ceroid lipofuscinosis neuronal protein 5 (Cln5) is encoded by the CLN5 gene. The genetic variants of this gene are associated with the CLN5 form of Batten disease. Recently, the first crystal structure of Cln5 was reported. Cln5 shows cysteine palmitoyl thioesterase S-depalmitoylation activity, which was explored via fluorescent emission spectroscopy utilizing the fluorescent probe DDP-5. In this work, the mechanism of the reaction between Cln5 and DDP-5 was studied computationally by applying a QM/MM methodology at the ωB97X-D/6-31G(d,p):AMBER level. The results of our study clearly demonstrate the critical role of the catalytic triad Cys280-His166-Glu183 in S-depalmitoylation activity. This is evidenced through a comparison of the pathways catalyzed by the Cys280-His166-Glu183 triad and those with only Cys280 involved. The computed reaction barriers are in agreement with the catalytic efficiency. The calculated Gibb's free-energy profile suggests that S-depalmitoylation is a rate-limiting step compared to the preceding S-palmitoylation, with barriers of 26.1 and 25.3 kcal/mol, respectively. The energetics were complemented by monitoring the fluctuations in the electron density distribution through NBO charges and bond strength alterations via local mode stretching force constants during the catalytic pathways. This comprehensive protocol led to a more holistic picture of the reaction mechanism at the atomic level. It forms the foundation for future studies on the effects of gene mutations on both the S-palmitoylation and S-depalmitoylation steps, providing valuable data for the further development of enzyme replacement therapy, which is currently the only FDA-approved therapy for childhood neurodegenerative diseases, including Batten disease.
Collapse
Affiliation(s)
- Yuvraj Dangat
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Marek Freindorf
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
7
|
Gao S, Guan H, Bloomer H, Wich D, Song D, Khirallah J, Ye Z, Zhao Y, Chen M, Xu C, Liu L, Xu Q. Harnessing non-Watson-Crick's base pairing to enhance CRISPR effectors cleavage activities and enable gene editing in mammalian cells. Proc Natl Acad Sci U S A 2024; 121:e2308415120. [PMID: 38150477 PMCID: PMC10786293 DOI: 10.1073/pnas.2308415120] [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: 05/19/2023] [Accepted: 11/21/2023] [Indexed: 12/29/2023] Open
Abstract
Genomic DNA of the cyanophage S-2L virus is composed of 2-aminoadenine (Z), thymine (T), guanine (G), and cytosine (C), forming the genetic alphabet ZTGC, which violates Watson-Crick base pairing rules. The Z-base has an extra amino group on the two position that allows the formation of a third hydrogen bond with thymine in DNA strands. Here, we explored and expanded applications of this non-Watson-Crick base pairing in protein expression and gene editing. Both ZTGC-DNA (Z-DNA) and ZUGC-RNA (Z-RNA) produced in vitro show detectable compatibility and can be decoded in mammalian cells, including Homo sapiens cells. Z-crRNA can guide CRISPR-effectors SpCas9 and LbCas12a to cleave specific DNA through non-Watson-Crick base pairing and boost cleavage activities compared to A-crRNA. Z-crRNA can also allow for efficient gene and base editing in human cells. Together, our results help pave the way for potential strategies for optimizing DNA or RNA payloads for gene editing therapeutics and give insights to understanding the natural Z-DNA genome.
Collapse
Affiliation(s)
- Shuliang Gao
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Huiwen Guan
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Hanan Bloomer
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Douglas Wich
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Donghui Song
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Jennifer Khirallah
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Zhongfeng Ye
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Yu Zhao
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Mengting Chen
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Chutian Xu
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Lihan Liu
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| | - Qiaobing Xu
- Department of Biomedical Engineering, Tufts University, Medford, MA02155
| |
Collapse
|
8
|
Drobotenko MI, Svidlov AA, Dorohova АA, Baryshev MG, Dzhimak SS. Medium viscosity influence on the open states genesis in a DNA molecule. J Biomol Struct Dyn 2023:1-9. [PMID: 38102872 DOI: 10.1080/07391102.2023.2294178] [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: 08/10/2023] [Accepted: 12/01/2023] [Indexed: 12/17/2023]
Abstract
The studies were carried out by the mathematical modeling of DNA mechanical deformations. Numerical calculations done for the interferon alpha 17 gene, which consists of 980 base pairs. It has been established that the genesis and dynamics of open states in the DNA molecule depends on the magnitude of the external influence (torque) and on the viscosity of the environment. In addition, it is shown that the dynamics of open states zones can have a jump-like character with a small change in the magnitude of the torque. When torque is applied to all 980 base pairs of the gene, the following effect is observed: an increase in the viscosity of the medium leads to an increase in the value of the torque necessary for the occurrence of OS and DNA unwinding, i.e. viscosity plays an important stabilizing role in DNA dynamics. Under the influence of a localized torque on different (by the content of A-T and G-C pairs and location) regions of the interferon alpha 17 gene, it was found that the magnitude of the external torque necessary for the occurrence of open states at all calculated values of viscosity depends on the nucleotide composition. The dependence of the torque magnitude required for the open states occurrence on viscosity is observed when the torque is applied to areas close to the gene boundaries. At the same time, the significance of the end effect, which weakens DNA, decreased with increasing viscosity of the medium.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Mikhail I Drobotenko
- Department of Radiophysics and Nanothechnology, Physics Faculty, Kuban State University, Krasnodar, Russia
| | - Alexandr A Svidlov
- Department of Radiophysics and Nanothechnology, Physics Faculty, Kuban State University, Krasnodar, Russia
- Laboratory of Problems of Stable Isotope Spreading in Living Systems, Federal Research Center the Southern Scientific Center of the Russian Academy of Sciences, Rostov-on-Don, Russia
| | - Аnna A Dorohova
- Department of Radiophysics and Nanothechnology, Physics Faculty, Kuban State University, Krasnodar, Russia
- Laboratory of Problems of Stable Isotope Spreading in Living Systems, Federal Research Center the Southern Scientific Center of the Russian Academy of Sciences, Rostov-on-Don, Russia
| | - Mikhail G Baryshev
- Department of Safety and Productivity of agroecosystems, All-Russian Research Institute of Phytopathology, Russia
| | - Stepan S Dzhimak
- Department of Radiophysics and Nanothechnology, Physics Faculty, Kuban State University, Krasnodar, Russia
- Laboratory of Problems of Stable Isotope Spreading in Living Systems, Federal Research Center the Southern Scientific Center of the Russian Academy of Sciences, Rostov-on-Don, Russia
| |
Collapse
|
9
|
Derewenda ZS. C-H Groups as Donors in Hydrogen Bonds: A Historical Overview and Occurrence in Proteins and Nucleic Acids. Int J Mol Sci 2023; 24:13165. [PMID: 37685972 PMCID: PMC10488043 DOI: 10.3390/ijms241713165] [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: 08/08/2023] [Revised: 08/21/2023] [Accepted: 08/23/2023] [Indexed: 09/10/2023] Open
Abstract
Hydrogen bonds constitute a unique type of non-covalent interaction, with a critical role in biology. Until fairly recently, the canonical view held that these bonds occur between electronegative atoms, typically O and N, and that they are mostly electrostatic in nature. However, it is now understood that polarized C-H groups may also act as hydrogen bond donors in many systems, including biological macromolecules. First recognized from physical chemistry studies, C-H…X bonds were visualized with X-ray crystallography sixty years ago, although their true significance has only been recognized in the last few decades. This review traces the origins of the field and describes the occurrence and significance of the most important C-H…O bonds in proteins and nucleic acids.
Collapse
Affiliation(s)
- Zygmunt Stanislaw Derewenda
- Department of Molecular Physiology and Biological Physics, School of Medicine, University of Virginia, Charlottesville, VA 22903-2628, USA
| |
Collapse
|
10
|
Freindorf M, McCutcheon M, Beiranvand N, Kraka E. Dihydrogen Bonding-Seen through the Eyes of Vibrational Spectroscopy. MOLECULES (BASEL, SWITZERLAND) 2022; 28:molecules28010263. [PMID: 36615456 PMCID: PMC9822382 DOI: 10.3390/molecules28010263] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Revised: 12/20/2022] [Accepted: 12/20/2022] [Indexed: 12/30/2022]
Abstract
In this work, we analyzed five groups of different dihydrogen bonding interactions and hydrogen clusters with an H3+ kernel utilizing the local vibrational mode theory, developed by our group, complemented with the Quantum Theory of Atoms-in-Molecules analysis to assess the strength and nature of the dihydrogen bonds in these systems. We could show that the intrinsic strength of the dihydrogen bonds investigated is primarily related to the protonic bond as opposed to the hydridic bond; thus, this should be the region of focus when designing dihydrogen bonded complexes with a particular strength. We could also show that the popular discussion of the blue/red shifts of dihydrogen bonding based on the normal mode frequencies is hampered from mode-mode coupling and that a blue/red shift discussion based on local mode frequencies is more meaningful. Based on the bond analysis of the H3+(H2)n systems, we conclude that the bond strength in these crystal-like structures makes them interesting for potential hydrogen storage applications.
Collapse
|
11
|
Madushanka A, Verma N, Freindorf M, Kraka E. Papaya Leaf Extracts as Potential Dengue Treatment: An In-Silico Study. Int J Mol Sci 2022; 23:12310. [PMID: 36293162 PMCID: PMC9610845 DOI: 10.3390/ijms232012310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 10/04/2022] [Accepted: 10/06/2022] [Indexed: 11/05/2022] Open
Abstract
Dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) cause serious public health problems, with nearly 390 million people affected and 20,000 deaths per year in tropical and subtropical countries. Despite numerous attempts, no antiviral drug or vaccine is currently available to combat the manifestation. The challenge of discovering an efficient vaccine is enhanced by the surplus presence of efficient vectors and drug resistance from the virus. For centuries, papaya (Carica papaya) extracts have been traditionally used to treat DF, DHF, and DSS. In the present study, we systematically investigated seven compounds isolated from papaya leaf extract with regard to their potential as inhibitors for non-structural (NS) proteins, NS3 and NS5, which play a crucial role in viral RNA replication. The computational tools applied stretched across classical molecular docking, molecular dynamics (MD) simulations and SwissADME used to calculate binding affinities; binding free energies; Absorption, Distribution, Metabolism, and Excretion (ADME); and drug-likeness properties, thus, identifying Kaempferol, Chlorogenic acid, and Quercetin as potential candidates, with Kaempferol and Quercetin scoring best. Therefore, for the Kaempferol and Quercetin complexes, hybrid quantum mechanical/molecular mechanical (QM/MM) geometry and frequency calculations were performed, followed by the local mode analysis developed in our group to quantify Kaempferol-NS and Quercetin-NS hydrogen bonding. Given the non-toxic nature and the wide availability of the Kaempferol and Quercetin papaya extract in almost all of the susceptible regions, and our results showing high NS3 and NS5 binding affinities and energies, strong hydrogen bonding with both NS3 and NS5, and excellent ADME properties, we suggest Kaempferol and Quercetin as a strong NS3 and NS5 inhibitor to be further investigated in vitro.
Collapse
Affiliation(s)
| | | | | | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, P.O. Box 750314, Dallas, TX 75275, USA
| |
Collapse
|
12
|
Freindorf M, Delgado AAA, Kraka E. CO bonding in hexa‐ and pentacoordinate carboxy‐neuroglobin: A quantum mechanics/molecular mechanics and local vibrational mode study. J Comput Chem 2022. [DOI: 10.1002/jcc.26973] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Marek Freindorf
- Department of Chemistry Southern Methodist University Dallas Texas USA
| | | | - Elfi Kraka
- Department of Chemistry Southern Methodist University Dallas Texas USA
| |
Collapse
|
13
|
Mechanistic Details of the Sharpless Epoxidation of Allylic Alcohols—A Combined URVA and Local Mode Study. Catalysts 2022. [DOI: 10.3390/catal12070789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Abstract
In this work, we investigated the catalytic effects of a Sharpless dimeric titanium (IV)–tartrate–diester catalyst on the epoxidation of allylalcohol with methyl–hydroperoxide considering four different orientations of the reacting species coordinated at the titanium atom (reactions R1–R4) as well as a model for the non-catalyzed reaction (reaction R0). As major analysis tools, we applied the URVA (Unified Reaction Valley Approach) and LMA (Local Mode Analysis), both being based on vibrational spectroscopy and complemented by a QTAIM analysis of the electron density calculated at the DFT level of theory. The energetics of each reaction were recalculated at the DLPNO-CCSD(T) level of theory. The URVA curvature profiles identified the important chemical events of all five reactions as peroxide OO bond cleavage taking place before the TS (i.e., accounting for the energy barrier) and epoxide CO bond formation together with rehybridization of the carbon atoms of the targeted CC double bond after the TS. The energy decomposition into reaction phase contribution phases showed that the major effect of the catalyst is the weakening of the OO bond to be broken and replacement of OH bond breakage in the non-catalyzed reaction by an energetically more favorable TiO bond breakage. LMA performed at all stationary points rounded up the investigation (i) quantifying OO bond weakening of the oxidizing peroxide upon coordination at the metal atom, (ii) showing that a more synchronous formation of the new CO epoxide bonds correlates with smaller bond strength differences between these bonds, and (iii) elucidating the different roles of the three TiO bonds formed between catalyst and reactants and their interplay as orchestrated by the Sharpless catalyst. We hope that this article will inspire the computational community to use URVA complemented with LMA in the future as an efficient mechanistic tool for the optimization and fine-tuning of current Sharpless catalysts and for the design new of catalysts for epoxidation reactions.
Collapse
|
14
|
Margaret McCutcheon M, Freindorf M, Kraka E. Bonding in Nitrile Photo-dissociating Ruthenium Drug Candidates --A Local Vibrational Mode Study. J Chem Phys 2022; 157:014301. [DOI: 10.1063/5.0094567] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we investigated bonding features 15 ruthenium complexes of the type [Ru(tpy)(L)-(CH3CN)]n+, containing the tridentate tpy ligand (tpy = 2,2':6',2'--terpyridine) and various bidentate ancillary ligands, 12 compounds originally synthesized by Loftus et al. (J. Phys. Chem. C 123, 10291-10299 (2019)) complemented with three additional complexes. The main focus of our work was to relate these local features to the experimental data of Loftus et al. which assess the efficiency of nitrile release in an indirect way via observed quantum yields for ruthenium water association after nitrile release. As a tool to quantitatively assess Ru-NC and Ru-L bonding we utilized the local vibrational mode analysis complemented by the topological analysis of the electron density and the natural bond orbital analysis. Interestingly, the stronger Ru-NC bonds have the greater observed quantum yields, leading to the conclusion that the observed quantum yields are a result of a complex interplay of several processes excluding a direct relationship between QY and Ru-NC or Ru-L bond strengths. We identified the ST splitting as one of the key players and not the Ru-NC bond strength, as one may have thought. In summary, this work has presented a modern computational tool set for the investigation of bonding features applied to nitrile photo-dissociating ruthenium drug candidates forming a valuable basis for future design and fine tuning of nitrile releasing ruthenium compounds, as well as for the understanding of how local properties affect overall experimental outcomes.
Collapse
Affiliation(s)
| | | | - Elfi Kraka
- Chemistry, Southern Methodist University, United States of America
| |
Collapse
|
15
|
Costa Peluzo BMT, Kraka E. Uranium: The Nuclear Fuel Cycle and Beyond. Int J Mol Sci 2022; 23:ijms23094655. [PMID: 35563047 PMCID: PMC9101921 DOI: 10.3390/ijms23094655] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Revised: 04/13/2022] [Accepted: 04/18/2022] [Indexed: 01/01/2023] Open
Abstract
This review summarizes the recent developments regarding the use of uranium as nuclear fuel, including recycling and health aspects, elucidated from a chemical point of view, i.e., emphasizing the rich uranium coordination chemistry, which has also raised interest in using uranium compounds in synthesis and catalysis. A number of novel uranium coordination features are addressed, such the emerging number of U(II) complexes and uranium nitride complexes as a promising class of materials for more efficient and safer nuclear fuels. The current discussion about uranium triple bonds is addressed by quantum chemical investigations using local vibrational mode force constants as quantitative bond strength descriptors based on vibrational spectroscopy. The local mode analysis of selected uranium nitrides, N≡U≡N, U≡N, N≡U=NH and N≡U=O, could confirm and quantify, for the first time, that these molecules exhibit a UN triple bond as hypothesized in the literature. We hope that this review will inspire the community interested in uranium chemistry and will serve as an incubator for fruitful collaborations between theory and experimentation in exploring the wealth of uranium chemistry.
Collapse
|
16
|
BF3–Catalyzed Diels–Alder Reaction between Butadiene and Methyl Acrylate in Aqueous Solution—An URVA and Local Vibrational Mode Study. Catalysts 2022. [DOI: 10.3390/catal12040415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
In this study we investigate the Diels–Alder reaction between methyl acrylate and butadiene, which is catalyzed by BF3 Lewis acid in explicit water solution, using URVA and Local Mode Analysis as major tools complemented with NBO, electron density and ring puckering analyses. We considered four different starting orientations of methyl acrylate and butadiene, which led to 16 DA reactions in total. In order to isolate the catalytic effects of the BF3 catalyst and those of the water environment and exploring how these effects are synchronized, we systematically compared the non-catalyzed reaction in gas phase and aqueous solution with the catalyzed reaction in gas phase and aqueous solution. Gas phase studies were performed at the B3LYP/6-311+G(2d,p) level of theory and studies in aqueous solution were performed utilizing a QM/MM approach at the B3LYP/6-311+G(2d,p)/AMBER level of theory. The URVA results revealed reaction path curvature profiles with an overall similar pattern for all 16 reactions showing the same sequence of CC single bond formation for all of them. In contrast to the parent DA reaction with symmetric substrates causing a synchronous bond formation process, here, first the new CC single bond on the CH2 side of methyl acrylate is formed followed by the CC bond at the ester side. As for the parent DA reaction, both bond formation events occur after the TS, i.e., they do not contribute to the energy barrier. What determines the barrier is the preparation process for CC bond formation, including the approach diene and dienophile, CC bond length changes and, in particular, rehybridization of the carbon atoms involved in the formation of the cyclohexene ring. This process is modified by both the BF3 catalyst and the water environment, where both work in a hand-in-hand fashion leading to the lowest energy barrier of 9.06 kcal/mol found for the catalyzed reaction R1 in aqueous solution compared to the highest energy barrier of 20.68 kcal/mol found for the non-catalyzed reaction R1 in the gas phase. The major effect of the BF3 catalyst is the increased mutual polarization and the increased charge transfer between methyl acrylate and butadiene, facilitating the approach of diene and dienophile and the pyramidalization of the CC atoms involved in the ring formation, which leads to a lowering of the activation energy. The catalytic effect of water solution is threefold. The polar environment leads also to increased polarization and charge transfer between the reacting species, similar as in the case of the BF3 catalyst, although to a smaller extend. More important is the formation of hydrogen bonds with the reaction complex, which are stronger for the TS than for the reactant, thus stabilizing the TS which leads to a further reduction of the activation energy. As shown by the ring puckering analysis, the third effect of water is space confinement of the reacting partners, conserving the boat form of the six-member ring from the entrance to the exit reaction channel. In summary, URVA combined with LMA has led to a clearer picture on how both BF3 catalyst and aqueous environment in a synchronized effort lower the reaction barrier. These new insights will serve to further fine-tune the DA reaction of methyl acrylate and butadiene and DA reactions in general.
Collapse
|
17
|
Fu TH, Lin MY, Fu CB, Yu XF, Xiao B, Cheng JB, Li Q. The role of nitro group on the excited-state relaxation mechanism of P-Z base pair. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2022; 267:120549. [PMID: 34810098 DOI: 10.1016/j.saa.2021.120549] [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: 05/25/2021] [Revised: 09/09/2021] [Accepted: 10/24/2021] [Indexed: 06/13/2023]
Abstract
DNAs' photostability is significant to the normal function of organisms. P-Z is a hydrogen bonded artificial DNA base pair, where P and Z represent 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)one and 6-amino-5nitro-2(1H)-pyridone, respectively. The excited-state relaxation mechanism of P-Z pair is investigated using static TDDFT calculations combined with the non-adiabatic dynamic simulations at TDDFT level. The roles of nitro rotation, nitro out-of-plane deformation, and single proton transfer (SPT) along hydrogen bond are revealed. The results of potential energy profile calculations demonstrate that the SPT processes along the hydrogen bonds are unfavorable to occur statically, which is in great contrast to the natural base pair. The non-adiabatic dynamic simulations show that the excited-state nitro rotation and nitro out-of-plane deformation are the two important relaxation channels which lead to the fast internal conversion to S0 state. The SPT from Z to P is also observed, followed by distortion on P, inducing the fast internal conversion to S0 state. However, this channel (decay via SPT process) plays minor roles on the excited-state relaxation mechanism statistically. This work shows the great differences of the excited-state relaxation mechanism between the natural base pairs and artificial base pair, also sheds new light into the role of hydrogen bond and nitro group in P-Z base pair.
Collapse
Affiliation(s)
- Ting-He Fu
- The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai 264005, People's Republic of China
| | - Man-Yu Lin
- The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai 264005, People's Republic of China
| | - Cheng-Bin Fu
- The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai 264005, People's Republic of China
| | - Xue-Fang Yu
- The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai 264005, People's Republic of China.
| | - Bo Xiao
- The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai 264005, People's Republic of China
| | - Jian-Bo Cheng
- The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai 264005, People's Republic of China
| | - Qingzhong Li
- The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai 264005, People's Republic of China
| |
Collapse
|
18
|
Nanayakkara S, Tao Y, Kraka E. Capturing Individual Hydrogen Bond Strengths in Ices via Periodic Local Vibrational Mode Theory: Beyond the Lattice Energy Picture. J Chem Theory Comput 2021; 18:562-579. [PMID: 34928619 DOI: 10.1021/acs.jctc.1c00357] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Local stretching force constants derived from periodic local vibrational modes at the vdW-DF2 density functional level have been employed to quantify the intrinsic hydrogen bond strength of 16 ice polymorphs, ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI, XII, XIII, XIV, XV, XVII, and XIX, that are stable under ambient to elevated pressures. Based on this characterization on 1820 hydrogen bonds, relationships between local stretching force constants and structural parameters such as hydrogen bond length and angle were identified. Moreover, different bond strength distributions, from uniform to inhomogeneous, were observed for the 16 ices and could be explained in relation to different local structural elements within ices, that is, rings, that consist of different hydrogen bond types. In addition, criteria for the classification of hydrogen bonds as strong, intermediate, and weak were introduced. The latter was used to explore a different dimension of the water-ice phase diagram. These findings will provide important guidelines for assessing the credibility of new ice structures.
Collapse
Affiliation(s)
- Sadisha Nanayakkara
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Yunwen Tao
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
19
|
Krüger R, Larroza A, Fronza MG, Tisoco I, Savegnago L, Reis JS, Back DF, Iglesias BA, Alves D. Bis-triazolylchalcogenium-Functionalized Benzothiadiazole Derivatives as Light-up Sensors for DNA and BSA. J Org Chem 2021; 86:17866-17883. [PMID: 34843245 DOI: 10.1021/acs.joc.1c02153] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
A range of bis-triazolylchalcogenium-BTD 3 was synthesized by a copper-catalyzed azide-alkyne cycloaddition of azido arylchalcogenides 1 and 4,7-diethynylbenzo[c][1,2,5]thiadiazole 2. Eight new compounds were obtained in moderate to good yields using 1 mol % of copper(II) acetate monohydrate under mild reaction conditions. In addition, the synthesized bis-triazolylchalcogenium-BTD 3a-3h were investigated regarding their photophysical, electrochemical, and biomolecule binding properties in solution. In general, compounds presented strong absorption bands at the 250-450 nm region and cyan to green emission properties. The redox process attributed to the chalcogen atom was observed by electrochemical analysis (CV techniques). In addition, spectroscopic studies by UV-vis, steady-state emission fluorescence, and molecular docking calculations evidenced the ability of each derivative to establish interactions with calf-thymus DNA (CT-DNA) and bovine serum albumin (BSA). The behavior presented for this new class of compounds makes them a promising tool as optical sensors for biomolecules.
Collapse
Affiliation(s)
- Roberta Krüger
- Clean Organic Synthesis Laboratory - LASOL - CCQFA - Federal University of Pelotas - UFPel, CEP, 96010-900 Pelotas, RS, Brazil
| | - Allya Larroza
- Clean Organic Synthesis Laboratory - LASOL - CCQFA - Federal University of Pelotas - UFPel, CEP, 96010-900 Pelotas, RS, Brazil
| | - Mariana G Fronza
- Neurobiotechnology Research Group - GPN, CDTec, Federal University of Pelotas, UFPel, CEP, 96010-900 Pelotas, RS, Brazil
| | - Isadora Tisoco
- Chemistry Department, Bioinorganics and Porphyrinic Materials Laboratory, Federal University of Santa Maria, UFSM, CEP, 97105-900 Santa Maria, RS, Brazil
| | - Lucielli Savegnago
- Neurobiotechnology Research Group - GPN, CDTec, Federal University of Pelotas, UFPel, CEP, 96010-900 Pelotas, RS, Brazil
| | - Joel S Reis
- Institute of Chemistry, Fluminense Federal University, UFF, CEP, 24020-141 Niteroi, RJ, Brazil
| | - Davi F Back
- Chemistry Department, Inorganic Materials Laboratory, Federal University of Santa Maria, UFSM, CEP, 97105-900 Santa Maria, RS, Brazil
| | - Bernardo A Iglesias
- Chemistry Department, Bioinorganics and Porphyrinic Materials Laboratory, Federal University of Santa Maria, UFSM, CEP, 97105-900 Santa Maria, RS, Brazil
| | - Diego Alves
- Clean Organic Synthesis Laboratory - LASOL - CCQFA - Federal University of Pelotas - UFPel, CEP, 96010-900 Pelotas, RS, Brazil
| |
Collapse
|
20
|
Freindorf M, Beiranvand N, Delgado AAA, Tao Y, Kraka E. On the formation of CN bonds in Titan's atmosphere-a unified reaction valley approach study. J Mol Model 2021; 27:320. [PMID: 34633543 DOI: 10.1007/s00894-021-04917-8] [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: 07/01/2021] [Accepted: 09/15/2021] [Indexed: 08/30/2023]
Abstract
In this work, we investigated the formation of protonated hydrogen cyanide HCNH+ and methylene amine cation CH[Formula: see text] (both identified in Titan's upper atmosphere) from three different pathways which stem from the interaction between CH4 and N+(3P). As a mechanistic tool, we used the Unified Reaction Valley Approach (URVA) complemented with the Local Mode Analysis (LMA) assessing the strength of the CN bonds formed in these reactions. Our URVA studies could provide a comprehensive overview on bond formation/cleavage processes relevant to the specific mechanism of eight reactions R1- R8 that occur across the three pathways. In addition, we could explain the formation of CH[Formula: see text] and the appearance of HCNH+ and CHNH[Formula: see text] along these paths. Although only smaller molecules are involved in these reactions including isomerization, hydrogen atom abstraction, and hydrogen molecule capture, we found a number of interesting features, such as roaming in reaction R3 or the primary interaction of H2 with the carbon atom in HCNH+ in reaction R8 followed by migration of one of the H2 hydrogen atoms to the nitrogen which is more cost effective than breaking the HH bond first; a feature often found in catalysis. In all cases, charge transfer between carbon and nitrogen could be identified as a driving force for the CN bond formation. As revealed by LMA, the CN bonds formed in reactions R1-R8 cover a broad bond strength range from very weak to very strong, with the CN bond in protonated hydrogen cyanide HCNH+ identified as the strongest of all molecules investigated in this work. Our study demonstrates the large potential of both URVA and LMA to shed new light into these extraterrestrial reactions to help better understand prebiotic processes as well as develop guidelines for future investigations involving areas of complex interstellar chemistry. In particular, the formation of CN bonds as a precursor to the extraterrestrial formation of amino acids will be the focus of future investigations. Formation of CN bonds in Titan's atmosphere visualized via the reaction path curvature.
Collapse
Affiliation(s)
- Marek Freindorf
- Chemistry Department, SMU, Fondren Science Building, Dallas, 75275-0314, TX, USA
| | - Nassim Beiranvand
- Chemistry Department, SMU, Fondren Science Building, Dallas, 75275-0314, TX, USA
| | - Alexis A A Delgado
- Chemistry Department, SMU, Fondren Science Building, Dallas, 75275-0314, TX, USA
| | - Yunwen Tao
- Chemistry Department, SMU, Fondren Science Building, Dallas, 75275-0314, TX, USA
| | - Elfi Kraka
- Chemistry Department, SMU, Fondren Science Building, Dallas, 75275-0314, TX, USA.
| |
Collapse
|
21
|
Oliveira VP, Marcial BL, Machado FBC, Kraka E. Relating Bond Strength and Nature to the Thermodynamic Stability of Hypervalent Togni-Type Iodine Compounds. Chempluschem 2021; 86:1199-1210. [PMID: 34437775 DOI: 10.1002/cplu.202100285] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2021] [Revised: 08/09/2021] [Indexed: 11/08/2022]
Abstract
The bond strength and nature of a set of 32 Togni-like reagents have been investigated at the M062X/def2-TZVP(D) level of theory in acetonitrile described with the SMD continuum solvent model, to rationalize the main factors responsible for their thermodynamic stability in different conformations, and trifluoromethylation capabilities. For the assessment of bond strength, we utilized local stretching force constants and associated bond strength orders, complemented with local features of the electron density to access the nature of the bonds. Bond dissociation energies varied from 31.6 to 79.9 kcal/mol depending on the polarizing power of the ligand trans to CF3 . Based on the analysis of the Laplacian of the density, we propose that the charge-shift bond character plays an important role in the stability of the molecules studied, especially for those containing I-O bonds. New insights on the trans influence and on possible ways to fine-tune the stability of these reagents are provided.
Collapse
Affiliation(s)
- Vytor Pinheiro Oliveira
- Instituto Tecnológico de Aeronáutica (ITA), Departamento de Química, São José dos Campos, 12228-900, São Paulo, Brazil
| | - Bruna Luana Marcial
- Instituto Federal Goiano (IF Goiano), Núcleo de Química, Campus Morrinhos, Goiás, Brazil
| | - Francisco B C Machado
- Instituto Tecnológico de Aeronáutica (ITA), Departamento de Química, São José dos Campos, 12228-900, São Paulo, Brazil
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas, 75275-0314, USA
| |
Collapse
|