1
|
Li J, Liang J, Wang Z, Ptaszek AL, Liu X, Ganoe B, Head-Gordon M, Head-Gordon T. Highly Accurate Prediction of NMR Chemical Shifts from Low-Level Quantum Mechanics Calculations Using Machine Learning. J Chem Theory Comput 2024; 20:2152-2166. [PMID: 38331423 DOI: 10.1021/acs.jctc.3c01256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Theoretical predictions of NMR chemical shifts from first-principles can greatly facilitate experimental interpretation and structure identification of molecules in gas, solution, and solid-state phases. However, accurate prediction of chemical shifts using the gold-standard coupled cluster with singles, doubles, and perturbative triple excitations [CCSD(T)] method with a complete basis set (CBS) can be prohibitively expensive. By contrast, machine learning (ML) methods offer inexpensive alternatives for chemical shift predictions but are hampered by generalization to molecules outside the original training set. Here, we propose several new ideas in ML of the chemical shift prediction for H, C, N, and O that first introduce a novel feature representation, based on the atomic chemical shielding tensors within a molecular environment using an inexpensive quantum mechanics (QM) method, and train it to predict NMR chemical shieldings of a high-level composite theory that approaches the accuracy of CCSD(T)/CBS. In addition, we train the ML model through a new progressive active learning workflow that reduces the total number of expensive high-level composite calculations required while allowing the model to continuously improve on unseen data. Furthermore, the algorithm provides an error estimation, signaling potential unreliability in predictions if the error is large. Finally, we introduce a novel approach to keep the rotational invariance of the features using tensor environment vectors (TEVs) that yields a ML model with the highest accuracy compared to a similar model using data augmentation. We illustrate the predictive capacity of the resulting inexpensive shift machine learning (iShiftML) models across several benchmarks, including unseen molecules in the NS372 data set, gas-phase experimental chemical shifts for small organic molecules, and much larger and more complex natural products in which we can accurately differentiate between subtle diastereomers based on chemical shift assignments.
Collapse
Affiliation(s)
- Jie Li
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Jiashu Liang
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Zhe Wang
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Aleksandra L Ptaszek
- Christian Doppler Laboratory for High-Content Structural Biology and Biotechnology, Department of Structural and Computational Biology, Max Perutz Laboratories, University of Vienna, Campus Vienna Biocenter 5, Vienna 1030, Austria
- Laboratory for Computer-Aided Molecular Design, Division of Medicinal Chemistry, Otto Loewi Research Center, Medical University Graz, Neue Stiftingtalstrasse 6/III, Graz 8010, Austria
| | - Xiao Liu
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Brad Ganoe
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Ma L, Geng Y, Zhang G, Hu Z, James TD, Wang X, Wang Z. Near-Infrared Bodipy-Based Molecular Rotors for β-Amyloid Imaging In Vivo. Adv Healthc Mater 2023; 12:e2300733. [PMID: 37523149 DOI: 10.1002/adhm.202300733] [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/09/2023] [Revised: 07/12/2023] [Indexed: 08/01/2023]
Abstract
β-amyloid (Aβ) is one of the important biomarkers for diagnosing Alzheimer's disease (AD). Many near-infrared probes based on the donor-π-acceptor structure have been developed to detect Aβ. Most reported Aβ probes are based on the N,N-dimethylamino group as the ideal donor, which is a widely accepted binding unit. As such, the development of fluorescent probes with improved binding units to detect Aβ is urgently required. Therefore, with this research three anchoring molecular rotor electron donors consisting of cyclic amines of different ring sizes are developed, namely five-membered ring (TPyr), six-membered ring (TPip), and seven-membered ring (THAI). These new anchored molecular rotors are connected to a 4,4-difluoro-4-bora-3a,4a-diaza-s-indacene (BODIPY) and named TPyrBDP, TPipBDP, and THAIBDP. These probes exhibit high affinities (from 28 to 54 nm) for Aβ1-42 aggregates. The six-membered ring dye TPipBDP exhibits the highest signal-to-noise (75.5-fold) and higher affinity (28.30 ± 5.94 nm). TPipBDP can cross the blood-brain barrier and exhibits higher fluorescence enhancement with APP/PS1 (AD) double transgenic (Tg) mice than with wild-type (WT) mice.
Collapse
Affiliation(s)
- Lijun Ma
- State Key Laboratory of Chemical Resource Engineering, College of Chemistry, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Yujie Geng
- State Key Laboratory of Chemical Resource Engineering, College of Chemistry, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Guoyang Zhang
- State Key Laboratory of Chemical Resource Engineering, College of Chemistry, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Ziwei Hu
- State Key Laboratory of Chemical Resource Engineering, College of Chemistry, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
| | - Tony D James
- Department of Chemistry, University of Bath, Bath, BA2 7AY, UK
- School of Chemistry and Chemical Engineering, Henan Normal University, Xinxiang, 453007, China
| | - Xuefei Wang
- School of Chemistry and Chemical Engineering, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Zhuo Wang
- State Key Laboratory of Chemical Resource Engineering, College of Chemistry, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, 100029, China
- Department of Chemistry, University of Bath, Bath, BA2 7AY, UK
| |
Collapse
|
3
|
Cohen RD, Wood JS, Lam YH, Buevich AV, Sherer EC, Reibarkh M, Williamson RT, Martin GE. DELTA50: A Highly Accurate Database of Experimental 1H and 13C NMR Chemical Shifts Applied to DFT Benchmarking. Molecules 2023; 28:molecules28062449. [PMID: 36985422 PMCID: PMC10051451 DOI: 10.3390/molecules28062449] [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/25/2023] [Revised: 02/23/2023] [Accepted: 02/28/2023] [Indexed: 03/30/2023] Open
Abstract
Density functional theory (DFT) benchmark studies of 1H and 13C NMR chemical shifts often yield differing conclusions, likely due to non-optimal test molecules and non-standardized data acquisition. To address this issue, we carefully selected and measured 1H and 13C NMR chemical shifts for 50 structurally diverse small organic molecules containing atoms from only the first two rows of the periodic table. Our NMR dataset, DELTA50, was used to calculate linear scaling factors and to evaluate the accuracy of 73 density functionals, 40 basis sets, 3 solvent models, and 3 gauge-referencing schemes. The best performing DFT methodologies for 1H and 13C NMR chemical shift predictions were WP04/6-311++G(2d,p) and ωB97X-D/def2-SVP, respectively, when combined with the polarizable continuum solvent model (PCM) and gauge-independent atomic orbital (GIAO) method. Geometries should be optimized at the B3LYP-D3/6-311G(d,p) level including the PCM solvent model for the best accuracy. Predictions of 20 organic compounds and natural products from a separate probe set had root-mean-square deviations (RMSD) of 0.07 to 0.19 for 1H and 0.5 to 2.9 for 13C. Maximum deviations were less than 0.5 and 6.5 ppm for 1H and 13C, respectively.
Collapse
Affiliation(s)
- Ryan D Cohen
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
- Department of Chemistry and Biochemistry, Seton Hall University, South Orange, NJ 07079, USA
| | - Jared S Wood
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
- Department of Chemistry and Biochemistry, University of North Carolina Wilmington, Wilmington, NC 28409, USA
| | - Yu-Hong Lam
- Department of Computational and Structural Chemistry, Merck & Co., Inc., Rahway, NJ 07065, USA
| | - Alexei V Buevich
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
| | - Edward C Sherer
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
| | - Mikhail Reibarkh
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
| | - R Thomas Williamson
- Department of Chemistry and Biochemistry, University of North Carolina Wilmington, Wilmington, NC 28409, USA
| | - Gary E Martin
- Department of Chemistry and Biochemistry, Seton Hall University, South Orange, NJ 07079, USA
| |
Collapse
|
4
|
Pérez-Conesa S, Keeler EG, Zhang D, Delemotte L, McDermott AE. Informing NMR experiments with molecular dynamics simulations to characterize the dominant activated state of the KcsA ion channel. J Chem Phys 2021; 154:165102. [PMID: 33940802 PMCID: PMC9250420 DOI: 10.1063/5.0040649] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Accepted: 02/26/2021] [Indexed: 11/14/2022] Open
Abstract
As the first potassium channel with an x-ray structure determined, and given its homology to eukaryotic channels, the pH-gated prokaryotic channel KcsA has been extensively studied. Nevertheless, questions related, in particular, to the allosteric coupling between its gates remain open. The many currently available x-ray crystallography structures appear to correspond to various stages of activation and inactivation, offering insights into the molecular basis of these mechanisms. Since these studies have required mutations, complexation with antibodies, and substitution of detergents in place of lipids, examining the channel under more native conditions is desirable. Solid-state nuclear magnetic resonance (SSNMR) can be used to study the wild-type protein under activating conditions (low pH), at room temperature, and in bacteriomimetic liposomes. In this work, we sought to structurally assign the activated state present in SSNMR experiments. We used a combination of molecular dynamics (MD) simulations, chemical shift prediction algorithms, and Bayesian inference techniques to determine which of the most plausible x-ray structures resolved to date best represents the activated state captured in SSNMR. We first identified specific nuclei with simulated NMR chemical shifts that differed significantly when comparing partially open vs fully open ensembles from MD simulations. The simulated NMR chemical shifts for those specific nuclei were then compared to experimental ones, revealing that the simulation of the partially open state was in good agreement with the SSNMR data. Nuclei that discriminate effectively between partially and fully open states belong to residues spread over the sequence and provide a molecular level description of the conformational change.
Collapse
Affiliation(s)
- Sergio Pérez-Conesa
- KTH Royal Institute of Technology, Science for Life Laboratory, Stockholm, Sweden
| | - Eric G. Keeler
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Dongyu Zhang
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Lucie Delemotte
- KTH Royal Institute of Technology, Science for Life Laboratory, Stockholm, Sweden
| | - Ann E. McDermott
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
5
|
Cárdenas R, Martínez-Seoane J, Amero C. Combining Experimental Data and Computational Methods for the Non-Computer Specialist. Molecules 2020; 25:E4783. [PMID: 33081072 PMCID: PMC7594097 DOI: 10.3390/molecules25204783] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Revised: 08/25/2020] [Accepted: 08/28/2020] [Indexed: 01/01/2023] Open
Abstract
Experimental methods are indispensable for the study of the function of biological macromolecules, not just as static structures, but as dynamic systems that change conformation, bind partners, perform reactions, and respond to different stimulus. However, providing a detailed structural interpretation of the results is often a very challenging task. While experimental and computational methods are often considered as two different and separate approaches, the power and utility of combining both is undeniable. The integration of the experimental data with computational techniques can assist and enrich the interpretation, providing new detailed molecular understanding of the systems. Here, we briefly describe the basic principles of how experimental data can be combined with computational methods to obtain insights into the molecular mechanism and expand the interpretation through the generation of detailed models.
Collapse
Affiliation(s)
| | | | - Carlos Amero
- Laboratorio de Bioquímica y Resonancia Magnética Nuclear, Centro de Investigaciones Químicas, Instituto de Investigación en Ciencias Básicas y Aplicadas, Universidad Autónoma del Estado de Morelos, Cuernavaca, Morelos 62209, Mexico; (R.C.); (J.M.-S.)
| |
Collapse
|
6
|
Wang D, Huang M, Gao D, Chen K, Xinxie, Xu W, Li X. Screening anti-TMV agents targeting tobacco mosaic virus helicase protein. PESTICIDE BIOCHEMISTRY AND PHYSIOLOGY 2020; 166:104449. [PMID: 32448412 DOI: 10.1016/j.pestbp.2019.07.017] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/03/2019] [Revised: 06/29/2019] [Accepted: 07/13/2019] [Indexed: 06/11/2023]
Abstract
Tobacco mosaic virus helicase (TMV-Hel) plays important roles in viral multiplication. TMV-Hel is a potential target of anti-TMV agents. Our previous studies expressed and purified TMV-Hel as target protein for cytosinpeptidemycin. In this study, we preform molecular docking to study the binding sites of commercial antiviral agents with TMV-Hel. Then we verify the interactions between the potential anti-TMV agents and TMV-Hel in vitro using Microscale Thermophoresis experiment and study the inhibiting expression of TMV-Hel with the potential anti-TMV agents in vivo using Western-blot (WB) method. The results showed that ribavirin bound to TMV-Hel with a dissociation constant of 1.55 μM by direct interaction with eight binding sites, which was consistent with the docking studies. Ribavirin inhibited the expression of TMV-Hel in Nicotiana benthamiana. Docking studies combined Microscale Thermophoresis and WB experiment provided a new method to screen anti-TMV agents targeting TMV-Hel.
Collapse
Affiliation(s)
- Dongmei Wang
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang 550025, PR China
| | - Maoxi Huang
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang 550025, PR China
| | - Di Gao
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang 550025, PR China
| | - Kai Chen
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang 550025, PR China
| | - Xinxie
- Department of Plant Pathology, College of Agriculture, Guizhou University, Guiyang 550025, PR China
| | - Weiming Xu
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang 550025, PR China.
| | - Xiangyang Li
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang 550025, PR China.
| |
Collapse
|
7
|
Jin X, Zhu T, Zhang JZH, He X. Automated Fragmentation QM/MM Calculation of NMR Chemical Shifts for Protein-Ligand Complexes. Front Chem 2018; 6:150. [PMID: 29868556 PMCID: PMC5952040 DOI: 10.3389/fchem.2018.00150] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Accepted: 04/16/2018] [Indexed: 01/13/2023] Open
Abstract
In this study, the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) method was applied for NMR chemical shift calculations of protein-ligand complexes. In the AF-QM/MM approach, the protein binding pocket is automatically divided into capped fragments (within ~200 atoms) for density functional theory (DFT) calculations of NMR chemical shifts. Meanwhile, the solvent effect was also included using the Poission-Boltzmann (PB) model, which properly accounts for the electrostatic polarization effect from the solvent for protein-ligand complexes. The NMR chemical shifts of neocarzinostatin (NCS)-chromophore binding complex calculated by AF-QM/MM accurately reproduce the large-sized system results. The 1H chemical shift perturbations (CSP) between apo-NCS and holo-NCS predicted by AF-QM/MM are also in excellent agreement with experimental results. Furthermore, the DFT calculated chemical shifts of the chromophore and residues in the NCS binding pocket can be utilized as molecular probes to identify the correct ligand binding conformation. By combining the CSP of the atoms in the binding pocket with the Glide scoring function, the new scoring function can accurately distinguish the native ligand pose from decoy structures. Therefore, the AF-QM/MM approach provides an accurate and efficient platform for protein-ligand binding structure prediction based on NMR derived information.
Collapse
Affiliation(s)
- Xinsheng Jin
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, East China Normal University, Shanghai, China
| | - Tong Zhu
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, East China Normal University, Shanghai, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China
| | - John Z. H. Zhang
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, East China Normal University, Shanghai, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China
- Department of Chemistry, New York University, New York, NY, United States
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy, School of Chemistry and Molecular Engineering, Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, East China Normal University, Shanghai, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China
- National Engineering Research Centre for Nanotechnology, Shanghai, China
| |
Collapse
|
8
|
Li X, Chen K, Gao D, Wang D, Huang M, Zhu H, Kang J. Binding studies between cytosinpeptidemycin and the superfamily 1 helicase protein of tobacco mosaic virus. RSC Adv 2018; 8:18952-18958. [PMID: 35539684 PMCID: PMC9080650 DOI: 10.1039/c8ra01466c] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2018] [Accepted: 05/14/2018] [Indexed: 11/21/2022] Open
Abstract
Tobacco mosaic virus (TMV) helicases play important roles in viral multiplication and interactions with host organisms. They can also be targeted by antiviral agents. Cytosinpeptidemycin has a good control effect against TMV. However, the mechanism of action is unclear. In this study, we expressed and purified TMV superfamily 1 helicase (TMV-Hel) and analyzed its three-dimensional structure. Furthermore, the binding interactions of TMV-Hel and cytosinpeptidemycin were studied. Microscale thermophoresis and isothermal titration calorimetry experiments showed that cytosinpeptidemycin bound to TMV-Hel with a dissociation constant of 0.24–0.44 μM. Docking studies provided further insights into the interaction of cytosinpeptidemycin with the His375 of TMV-Hel. Mutational and Microscale thermophoresis analyses showed that cytosinpeptidemycin bound to a TMV-Hel mutant (H375A) with a dissociation constant of 14.5 μM. Thus, His375 may be the important binding site for cytosinpeptidemycin. The data are important for designing and synthesizing new effective antiphytoviral agents. Tobacco mosaic virus (TMV) helicases play important roles in viral multiplication and interactions with host organisms.![]()
Collapse
Affiliation(s)
- Xiangyang Li
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering
- Key Laboratory of Green Pesticide and Agricultural Bioengineering
- Ministry of Education
- Guizhou University
- Guiyang 550025
| | - Kai Chen
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering
- Key Laboratory of Green Pesticide and Agricultural Bioengineering
- Ministry of Education
- Guizhou University
- Guiyang 550025
| | - Di Gao
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering
- Key Laboratory of Green Pesticide and Agricultural Bioengineering
- Ministry of Education
- Guizhou University
- Guiyang 550025
| | - Dongmei Wang
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering
- Key Laboratory of Green Pesticide and Agricultural Bioengineering
- Ministry of Education
- Guizhou University
- Guiyang 550025
| | - Maoxi Huang
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering
- Key Laboratory of Green Pesticide and Agricultural Bioengineering
- Ministry of Education
- Guizhou University
- Guiyang 550025
| | - Hengmin Zhu
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering
- Key Laboratory of Green Pesticide and Agricultural Bioengineering
- Ministry of Education
- Guizhou University
- Guiyang 550025
| | - Jinxin Kang
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering
- Key Laboratory of Green Pesticide and Agricultural Bioengineering
- Ministry of Education
- Guizhou University
- Guiyang 550025
| |
Collapse
|