1
|
Synthesis, cytotoxic activity and quantum chemical calculations of new 7-thioxopyrazolo[1,5-f]pyrimidin-2-one derivatives. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2019.127261] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
2
|
Abstract
To accurately predict the free energy barrier for urea elimination in aqueous solution, we examined the reaction coordinates for the direct and water-assisted elimination pathways, and evaluated the corresponding free energy barriers by using the surface and volume polarization for electrostatics (SVPE) model-based first-principles electronic-structure calculations. Based on the computational results, the water-assisted elimination pathway is dominant for urea elimination in aqueous solution, and the corresponding free energy barrier is 25.3 kcal/mol. The free energy barrier of 25.3 kcal/mol predicted for the dominant reaction pathway of urea elimination in aqueous solution is in good agreement with available experimental kinetic data.
Collapse
|
3
|
Yao M, Tu W, Chen X, Zhan CG. Reaction pathways and free energy profiles for spontaneous hydrolysis of urea and tetramethylurea: unexpected substituent effects. Org Biomol Chem 2014; 11:7595-605. [PMID: 24097048 DOI: 10.1039/c3ob41055b] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
It has been difficult to directly measure the spontaneous hydrolysis rate of urea and, thus, 1,1,3,3-tetramethylurea (Me4U) was used as a model to determine the "experimental" rate constant for urea hydrolysis. The use of Me4U was based on an assumption that the rate of urea hydrolysis should be 2.8 times that of Me4U hydrolysis because the rate of acetamide hydrolysis is 2.8 times that of N,N-dimethyl-acetamide hydrolysis. The present first-principles electronic-structure calculations on the competing non-enzymatic hydrolysis pathways have demonstrated that the dominant pathway is the neutral hydrolysis via the CN addition for both urea (when pH < ~11.6) and Me4U (regardless of pH), unlike the non-enzymatic hydrolysis of amides where alkaline hydrolysis is dominant. Based on the computational data, the substituent shift of the free energy barrier calculated for the neutral hydrolysis is remarkably different from that for the alkaline hydrolysis, and the rate constant for the urea hydrolysis should be ~1.3 × 10(9)-fold lower than that (4.2 × 10(-12) s(-1)) measured for the Me4U hydrolysis. As a result, the rate enhancement and catalytic proficiency of urease should be 1.2 × 10(25) and 3 × 10(27) M(-1), respectively, suggesting that urease surpasses proteases and all other enzymes in its power to enhance the rate of reaction. All of the computational results are consistent with available experimental data for Me4U, suggesting that the computational prediction for urea is reliable.
Collapse
Affiliation(s)
- Min Yao
- Key Laboratory of Pesticide and Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P. R. China
| | | | | | | |
Collapse
|
4
|
Yang J, Wang G, Wang F, Gong X. An improved approach for predicting the density of azido compounds. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.824570] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
5
|
DFT study of the interaction between the conjugated fluorescein and dabcyl system, using fluorescene quenching method. J Mol Model 2012; 18:4113-20. [DOI: 10.1007/s00894-012-1413-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2012] [Accepted: 03/22/2012] [Indexed: 10/28/2022]
|
6
|
|
7
|
Li Q, Wei Y, Guo H, Zhan CG. Syntheses, structural characterizations and electronic absorption spectra simulation of three phenylimido substituted hexamolybdates incorporating a remote chloro group. Inorganica Chim Acta 2008. [DOI: 10.1016/j.ica.2007.11.031] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
8
|
Lu H, Chen X, Zhan CG. First-principles calculation of pKa for cocaine, nicotine, neurotransmitters, and anilines in aqueous solution. J Phys Chem B 2007; 111:10599-605. [PMID: 17691837 PMCID: PMC2882246 DOI: 10.1021/jp072917r] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The absolute pKa values of 24 representative amine compounds, including cocaine, nicotine, 10 neurotransmitters, and 12 anilines, in aqueous solution were calculated by performing first-principles electronic structure calculations that account for the solvent effects using four different solvation models, i.e., the surface and volume polarization for electrostatic interaction (SVPE) model, the standard polarizable continuum model (PCM), the integral equation formalism for the polarizable continuum model (IEFPCM), and the conductor-like screening solvation model (COSMO). Within the examined computational methods, the calculations using the SVPE model lead to the absolute pKa values with the smallest root-mean-square-deviation (rmsd) value (1.18). When the SVPE model was replaced by the PCM, IEFPCM, and COSMO, the rmsd value of the calculated absolute pKa values became 3.21, 2.72, and 3.08, respectively. All types of calculated pKa values linearly correlate with the experimental pKa values very well. With the empirical corrections using the linear correlation relationships, the theoretical pKa values are much closer to the corresponding experimental data and the rmsd values become 0.51-0.83. The smallest rmsd value (0.51) is also associated with the SVPE model. All of the results suggest that the first-principles electronic structure calculations using the SVPE model are a reliable approach to the pKa prediction for the amine compounds.
Collapse
Affiliation(s)
- Haiting Lu
- Key Laboratory of Pesticide & Chemical Biology of the Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, 430079, P. R. China
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536
| | - Xi Chen
- Key Laboratory of Pesticide & Chemical Biology of the Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, 430079, P. R. China
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536
| |
Collapse
|
9
|
Wu W, Li W, Wang L, Zhang P, Zhang J. Density functional theory study on lactides: Geometries, IR, NMR and electronic spectra. ACTA ACUST UNITED AC 2007. [DOI: 10.1016/j.theochem.2007.03.033] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
10
|
Xiong Y, Zhan CG. Theoretical studies of the transition-state structures and free energy barriers for base-catalyzed hydrolysis of amides. J Phys Chem A 2007; 110:12644-52. [PMID: 17107116 PMCID: PMC2892839 DOI: 10.1021/jp063140p] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The transition-state structures and free energy barriers for the rate-determining step (i.e. the formation of a tetrahedral intermediate) of base-catalyzed hydrolysis of a series of amides in aqueous solution have been studied by performing first-principle electronic structure calculations using a hybrid supermolecule-polarizable continuum approach. The calculated results and a revisit of recently reported experimental proton inventory data reveal that the favorable transition-state structure optimized for the tetrahedral intermediate formation of hydroxide ion-catalyzed hydrolysis of formamide may have three solvating water molecules remaining on the attacking hydroxide oxygen and two additional water molecules attached to the carbonyl oxygen of formamide. The calculated results have also demonstrated interesting substituent effects on the optimized transition-state geometries, on the transition-state stabilization, and on the calculated free energy barriers for the base-catalyzed hydrolysis of amides. When some or all of the hydrogen atoms of formamide are replaced by methyl groups, the total number of water molecules hydrogen-bonding with the attacking hydroxide in the transition state decreases from three for formamide to two for N-methylacetamide, N,N-dimethylformamide (DMF), and N,N-dimethylacetamide (DMA). The larger substituents of the amide hinder the solvent water molecules approaching the attacking hydroxide oxygen in the transition state and, therefore, destabilize the transition-state structure and increase the free energy barrier. By using the optimized most favorable transition-state structures, the calculated free energy barriers, i.e., 21.6 (or 21.7), 22.7, 23.1, and 26.0 kcal/mol for formamide, N-methylacetamide, DMF, and DMA, respectively, are in good agreement with the available experimental free energy barriers, i.e., 21.2, 21.5, 22.6, and 24.1 kcal/mol for formamide, N-methylacetamide, DMF, and DMA, respectively.
Collapse
Affiliation(s)
| | - Chang-Guo Zhan
- Corresponding author. Correspondence: Chang-Guo Zhan, Ph.D., Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536, TEL: 859-323-3943, FAX: 859-323-3575,
| |
Collapse
|
11
|
Huang X, Zheng F, Crooks PA, Dwoskin L, Zhan CG. Modeling multiple species of nicotine and deschloroepibatidine interacting with alpha4beta2 nicotinic acetylcholine receptor: from microscopic binding to phenomenological binding affinity. J Am Chem Soc 2005; 127:14401-14. [PMID: 16218635 PMCID: PMC3182463 DOI: 10.1021/ja052681+] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A variety of molecular modeling, molecular docking, and first-principles electronic structure calculations were performed to study how the alpha4beta2 nicotinic acetylcholine receptor (nAChR) binds with different species of two typical agonists, (S)-(-)-nicotine and (R)-(-)-deschloroepibatidine, each of which is distinguished by different free bases and protonation states. On the basis of these results, predictions were made regarding the corresponding microscopic binding free energies. Hydrogen-bonding and cation-pi interactions between the receptor and the respective ligands were found to be the dominant factors differentiating the binding strengths of different microscopic binding species. The calculated results and analyses demonstrate that, for each agonist, all the species are interchangeable and can quickly achieve a thermodynamic equilibrium in solution and at the nAChR binding site. This allows quantitation of the equilibrium concentration distributions of the free ligand species and the corresponding microscopic ligand-receptor binding species, their pH dependence, and their contributions to the phenomenological binding affinity. The predicted equilibrium concentration distributions, pK(a) values, absolute phenomenological binding affinities, and their pH dependence are all in good agreement with available experimental data, suggesting that the computational strategy from the microscopic binding species and affinities to the phenomenological binding affinity is reliable for studying alpha4beta2 nAChR-ligand binding. This should provide valuable information for future rational design of drugs targeting nAChRs. The general strategy of the "from-microscopic-to-phenomenological" approach for studying interactions of alpha4beta2 nAChRs with (S)-(-)-nicotine and (R)-(-)-deschloroepibatidine may also be useful in studying other types of ligand-protein interactions involving multiple molecular species of a ligand and in associated rational drug design.
Collapse
Affiliation(s)
- Xiaoqin Huang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, Kentucky 40536
| | - Fang Zheng
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, Kentucky 40536
| | - Peter A. Crooks
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, Kentucky 40536
| | - Linda Dwoskin
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, Kentucky 40536
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, Kentucky 40536
| |
Collapse
|
12
|
Zhan CG, Deng SX, Skiba JG, Hayes BA, Tschampel SM, Shields GC, Landry DW. First-principle studies of intermolecular and intramolecular catalysis of protonated cocaine. J Comput Chem 2005; 26:980-6. [PMID: 15880781 PMCID: PMC2875688 DOI: 10.1002/jcc.20241] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
We have performed a series of first-principles electronic structure calculations to examine the reaction pathways and the corresponding free energy barriers for the ester hydrolysis of protonated cocaine in its chair and boat conformations. The calculated free energy barriers for the benzoyl ester hydrolysis of protonated chair cocaine are close to the corresponding barriers calculated for the benzoyl ester hydrolysis of neutral cocaine. However, the free energy barrier calculated for the methyl ester hydrolysis of protonated cocaine in its chair conformation is significantly lower than for the methyl ester hydrolysis of neutral cocaine and for the dominant pathway of the benzoyl ester hydrolysis of protonated cocaine. The significant decrease of the free energy barrier, approximately 4 kcal/mol, is attributed to the intramolecular acid catalysis of the methyl ester hydrolysis of protonated cocaine, because the transition state structure is stabilized by the strong hydrogen bond between the carbonyl oxygen of the methyl ester moiety and the protonated tropane N. The relative magnitudes of the free energy barriers calculated for different pathways of the ester hydrolysis of protonated chair cocaine are consistent with the experimental kinetic data for cocaine hydrolysis under physiologic conditions. Similar intramolecular acid catalysis also occurs for the benzoyl ester hydrolysis of (protonated) boat cocaine in the physiologic condition, although the contribution of the intramolecular hydrogen bonding to transition state stabilization is negligible. Nonetheless, the predictability of the intramolecular hydrogen bonding could be useful in generating antibody-based catalysts that recruit cocaine to the boat conformation and an analog that elicited antibodies to approximate the protonated tropane N and the benzoyl O more closely than the natural boat conformer might increase the contribution from hydrogen bonding. Such a stable analog of the transition state for intramolecular catalysis of cocaine benzoyl-ester hydrolysis was synthesized and used to successfully elicit a number of anticocaine catalytic antibodies.
Collapse
Affiliation(s)
- Chang-Guo Zhan
- Division of Clinical Pharmacology and Experimental Therapeutics, College of Physicians & Surgeons, Columbia University, 630 West 168th Street, New York, New York 10032, USA.
| | | | | | | | | | | | | |
Collapse
|
13
|
Xiong Y, Zhan CG. Reaction Pathways and Free Energy Barriers for Alkaline Hydrolysis of Insecticide 2-Trimethylammonioethyl Methylphosphonofluoridate and Related Organophosphorus Compounds: Electrostatic and Steric Effects. J Org Chem 2004; 69:8451-8. [PMID: 15549820 DOI: 10.1021/jo0487597] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Reaction pathways and free energy barriers for alkaline hydrolysis of the highly neurotoxic insecticide 2-trimethylammonioethyl methylphosphonofluoridate and related organophosphorus compounds were studied by performing first-principles electronic structure calculations on representative methylphosphonofluoridates, (RO)CH3P(O)F, in which R = CH2CH2N+(CH3)3, CH3, CH2CH2C(CH3)3, CH2CH2CH(CH3)2, CH(CH3)CH2N+(CH3)3, and CH(CH3)CH2N(CH3)2. The dominant reaction pathway was found to be associated with a transition state in which the attacking nucleophile OH- and the leaving group F- are positioned on opposite sides of the plane formed by the three remaining atoms attached to the phosphorus in order to minimize the electrostatic repulsion between these two groups. The free energy barriers calculated for the rate-determining step of the dominant pathway are 12.5 kcal/mol when R = CH2CH2N+(CH3)3, 15.5 kcal/mol when R = CH3, 17.9 kcal/mol when R = CH2CH2C(CH3)3, 16.5 kcal/mol when R = CH2CH2CH(CH3)2, 13.4 kcal/mol when R = CH(CH3)CH2N+(CH3)3, and 18.7 kcal/mol when R = CH(CH(3))CH(2)N(CH(3))(2). The calculated free energy barriers are in good agreement with available experimentally derived activation free energies, i.e. 14.7 kcal/mol when R = CH(3), 13.4 kcal/mol when R = CH2CH2N+(CH3)3, and 13.9 kcal/mol when R = CH(CH3)CH2N+(CH3)3. A detailed analysis of the calculated energetic results and available experimental data suggests that the net charge of the molecule (M) being hydrolyzed is a prominent factor affecting the free energy barrier (DeltaG) for the alkaline hydrolysis of phosphodiesters, phosphonofluoridates, and related organophosphorus compounds. The electrostatic interactions between the attacking nucleophile OH- and the molecule M being hydrolyzed favor such an order of the free energy barrier: DeltaG(M(+)+OH-) < DeltaG(M0+OH-) < DeltaG(M(-)+OH-), where M+, M0, and M- represent the cationic, neutral, and anionic molecules, respectively. The change of the substituent R in (RO)CH(3)P(O)F from CH3 to CH2CH2N+(CH3)3 is associated with both the electrostatic and steric effects on the free energy barrier, but the electrostatic effect dominates the substituent shift of the free energy barrier. This helps to better understand why the alkaline hydrolysis of (RO)CH3P(O)F with R = CH2CH2N+(CH3)3 and CH(CH3)CH2N+(CH3)3 is significantly faster than that with R = CH3. The effect of electrostatic interaction also helps to understand why the rate constants for the alkaline hydrolysis of phosphodiesters, such as intramolecular second messenger adenosine 3',5'-phosphate (cAMP), are generally smaller than those for the alkaline hydrolysis of the phosphonofluoridates and related phosphotriesters.
Collapse
Affiliation(s)
- Ying Xiong
- Key Laboratory of Pesticide & Chemical Biology of the Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P. R. China
| | | |
Collapse
|