1
|
Shu Q, Huang P, Dong Z, Wang W. Molecular dynamics investigation on synthesis of a pH- and temperature-sensitive carbon nanotube loaded with doxorubicin. iScience 2024; 27:108812. [PMID: 38303688 PMCID: PMC10831279 DOI: 10.1016/j.isci.2024.108812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Revised: 12/21/2023] [Accepted: 01/02/2024] [Indexed: 02/03/2024] Open
Abstract
The many exotic properties of carbon nanotubes (CNTs) make them a powerful attraction in the field of drug delivery systems (DDS). In this work, based on quantum chemical calculation and molecular simulation techniques, polyacrylic acid (PAA) and N-isopropyl acrylamide (NIP) are selected and acted simultaneously on the CNT to form a stable system (FCNT). As a potential DDS, FCNT captures the dispersed doxorubicin (DOX) molecules around it and maintains a stable configuration. In these processes, electrostatic and van der Waals forces act synergistically, with van der Waals forces dominating. Compared to NIP, PAA molecules exhibit stronger adhesion and encapsulation efficiency to CNT and stronger adsorption capacity to DOX. This study reveals the mechanism of action among PAA, NIP, CNT, and DOX, providing feasibility verification and prospective guidance for the experimental synthesis of PAA-NIP-CNT-type multifunctional DDS, and also broadening the idea for exploring more efficient DDS suitable for DOX.
Collapse
Affiliation(s)
- Qijiang Shu
- Institute of Information, Yunnan University of Chinese Medicine, Kunming, Yunnan 650500, China
- Yunnan Key Laboratory of Southern Medicinal Utilization, Yunnan University of Chinese Medicine, Kunming, Yunnan 650500, China
- Yunnan Traditional Chinese Medicine Prevention and Treatment Engineering Research Center, Yunnan University of Chinese Medicine, Kunming, Yunnan 650500, China
| | - Pengru Huang
- Guangxi Key Laboratory of Information Materials and Guangxi Collaborative Innovation Center of Structure and Property for New Energy and Materials, School of Material Science & Engineering, Guilin University of Electronic Technology, Guilin, Guangxi 541004, China
| | - Zhi Dong
- College of Chinese Materia Medica, Yunnan University of Chinese Medicine, Kunming, Yunnan 650500, China
| | - Wenping Wang
- College of Chinese Materia Medica, Yunnan University of Chinese Medicine, Kunming, Yunnan 650500, China
| |
Collapse
|
2
|
Karthikeyan M, Dhinesh Kumar M, Kaniraja G, Karunakaran C. Theoretical investigations of free energy of binding and chiral recognition studies of (R)- and (S)-Noradrenaline towards β-cyclodextrin. J Mol Graph Model 2023; 124:108552. [PMID: 37379759 DOI: 10.1016/j.jmgm.2023.108552] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Revised: 06/08/2023] [Accepted: 06/13/2023] [Indexed: 06/30/2023]
Abstract
Noradrenaline (NA), one of the important excitatory catecholamine neurotransmitters, is used as a medication for Parkinson's Disease (PD). The β-cyclodextrin (β-CD) is one of the most effective drug carrier & also used in chiral separation. So, in this theoretical investigation, the R/S-Noradrenaline (R/S-NA) forms binding & chiral recognition mechanisms and energies with β-CD were explored. Using the AutoDock, R/S forms were first docked into the cavity of β-CD giving host-guest complexes with the free energy of binding for S-NA (-4.81 kcal/mol) larger than R-NA (-4.53 kcal/mol). The host-guest inclusion 1:1 complexes between R/S-NA and β-CD have been also modeled and optimized with ONIOM2 (B3LYP/6-31g++DP: PM6) method by using the Gaussian software. Further, frequency calculations were carried out to obtain the free energies. In comparison to the R-NA (-54.59 kcal/mol), it was observed that the S-NA (-56.48 kcal/mol) with β-CD is more stable. Furthermore, the H-bond results from molecular dynamics simulation revealed that S-NA/β-CD was more stable than R-NA/β-CD. In addition, the thermodynamic properties, vibrational analysis (IR), HOMO-LUMO band gap energy, inter molecular hydrogen bond interactions, and conformational analysis were investigated for both the R/S forms to support & compare the stability of the inclusion complex. These inclusion & high stability of S-NA/β-CD and in turn its theoretical chiral recognition behavior observed agreeing well with the reported NMR experimental data have implications in drug delivery and chiral separation research.
Collapse
Affiliation(s)
- Murugesan Karthikeyan
- Biomedical Research Laboratory, Department of Chemistry, Virudhunagar Hindu Nadars' Senthikumara Nadar College (Autonomous & Affiliated to Madurai Kamaraj University), Virudhunagar- 626 001, Tamil Nadu, India
| | - Marimuthu Dhinesh Kumar
- Biomedical Research Laboratory, Department of Chemistry, Virudhunagar Hindu Nadars' Senthikumara Nadar College (Autonomous & Affiliated to Madurai Kamaraj University), Virudhunagar- 626 001, Tamil Nadu, India
| | - Ganesan Kaniraja
- Biomedical Research Laboratory, Department of Chemistry, Virudhunagar Hindu Nadars' Senthikumara Nadar College (Autonomous & Affiliated to Madurai Kamaraj University), Virudhunagar- 626 001, Tamil Nadu, India
| | - Chandran Karunakaran
- Biomedical Research Laboratory, Department of Chemistry, Virudhunagar Hindu Nadars' Senthikumara Nadar College (Autonomous & Affiliated to Madurai Kamaraj University), Virudhunagar- 626 001, Tamil Nadu, India.
| |
Collapse
|
3
|
Qiu S, Zhao X, Li Y, Ding W, Huang J. Simulation of Internal Defects in TKX-50 Crystals. MATERIALS (BASEL, SWITZERLAND) 2023; 16:ma16114063. [PMID: 37297197 DOI: 10.3390/ma16114063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Revised: 05/28/2023] [Accepted: 05/29/2023] [Indexed: 06/12/2023]
Abstract
1,1'-Dihydroxy-5,5'-bi-tetrazolium dihydroxylamine salt (TKX-50) is a new type of high-energy low-sense explosive with great application value, but TKX-50 made directly from the reaction has problems such as irregular crystal morphology and relatively large length-diameter, and these factors seriously affect the sensitivity of TKX-50 and limit its large-scale application. The internal defects of TKX-50 crystals have a great influence on their weakness, and studying its related properties is of great theoretical significance and application value. To further investigate the microscopic properties of TKX-50 crystals and to explore the connection between microscopic parameters and macroscopic susceptibility, this paper reports the use of molecular dynamics simulations to construct TKX-50 crystal scaling models with three types of defects-vacancy, dislocation and doping-and conducts molecular dynamics simulations. The influence of TKX-50 crystal defects on the initiation bond length, density, bonding diatomic interaction energy, and cohesive energy density of the crystal was obtained. The simulation results show that the models with a higher bond length of the initiator bond and higher percentage activated the initiator's N-N bond and lowered the bond-linked diatomic energy, cohesive energy density, and density corresponding to higher crystal sensitivities. This led to a preliminary connection between TKX-50 microscopic model parameters and macroscopic susceptibility. The results of the study can provide a reference for the design of subsequent experiments, and the research method can be extended to the research work on other energy-containing materials.
Collapse
Affiliation(s)
- Siqi Qiu
- State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
| | - Xue Zhao
- State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
| | - Yuanyuan Li
- State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
| | - Wenyuan Ding
- State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
| | - Junrui Huang
- State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
| |
Collapse
|
4
|
Saleh OA, Badawey AM, Aboul-Enein HY, Fouad MA. Enantioseparation, quantification, molecular docking and molecular dynamics study of five β-adrenergic blockers on Lux-Cellulose-2 column. BMC Chem 2023; 17:22. [PMID: 36927568 PMCID: PMC10018884 DOI: 10.1186/s13065-023-00925-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2022] [Accepted: 02/25/2023] [Indexed: 03/18/2023] Open
Abstract
Enantioseparation of five β-adrenergic blockers was studied using two mobile phases on a cellulose tris(3-chloro-4-methylphenylcarbamate) (Lux-Cellulose-2) chiral column in normal phase mode. The first mobile phase composed of n-hexane: ethanol: diethylamine 60: 40: 0.1 by volume has successfully resolved the chromatographic peaks of three pairs of β-adrenergic blockers namely, bisoprolol, carvedilol and atenolol. A mixture of n-hexane: ethanol: diethyl amine 75: 25: 0.1 by volume was used as the second mobile phase to separate the four pairs of enantiomers, metoprolol, carvedilol, nebivolol and atenolol with high resolution values. The mobile phases were pumped at a flow rate 1 mL/min with column temperature 25 °C using a UV detector at 230 nm. Molecular docking simulations of the five pairs of enantiomers was carried out in the cavities of the chiral stationary phase to gain a better understanding of the interaction between analyte enantiomers and chiral stationary phase and to better understand the mechanism of chiral recognition. According to the results, hydrogen bond interactions and π-π- interactions were the main types of interaction involved in the chiral recognition. Molecular dynamics simulation was performed to investigate the solvent effect on the interaction of the five pair of enantiomers in the chiral stationary phase cavity under dynamic conditions.
Collapse
Affiliation(s)
- Ola Ahmed Saleh
- Medicinal and Pharmaceutical Chemistry Department, Pharmaceutical and Drug Industries Research Division, National Research Centre (ID: 60014618), P.O. Box 12622, Giza, Egypt.
| | - Amr Mohamed Badawey
- Analytical Chemistry Department, Faculty of Pharmacy, Cairo University, Kasr El-Aini St., Cairo, 11562, Egypt
| | - Hassan Y Aboul-Enein
- Medicinal and Pharmaceutical Chemistry Department, Pharmaceutical and Drug Industries Research Division, National Research Centre (ID: 60014618), P.O. Box 12622, Giza, Egypt.
| | - Marwa Ahmed Fouad
- Pharmaceutical Chemistry Department, Faculty of Pharmacy, Cairo University, Kasr El-Aini St., Cairo, 11562, Egypt.,Pharmaceutical Chemistry Department, School of Pharmacy, NewGiza University, NewGiza, Km 22 Cairo-Alexandria Desert Road, Cairo, Egypt
| |
Collapse
|
5
|
Li J, Gao H, Ye Z, Deng J, Ouyang D. In silico formulation prediction of drug/cyclodextrin/polymer ternary complexes by machine learning and molecular modeling techniques. Carbohydr Polym 2022; 275:118712. [PMID: 34742437 DOI: 10.1016/j.carbpol.2021.118712] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2021] [Revised: 09/13/2021] [Accepted: 09/20/2021] [Indexed: 01/17/2023]
Abstract
Ternary cyclodextrin (CD) complexes (drug/CD/polymer) can effectively improve the solubility of water-insoluble drugs with large size than binary CD formulations. However, ternary formulations are screened by a trial-and-error approach, which is laborious and material-wasting. Current research aims to develop a prediction model for ternary CD formulations by combined machine learning and molecular modeling. 596 ternary formulations data were collected to build a prediction model by machine learning. The random forest model achieved good performance with R2 = 0.887 in ST prediction and R2 = 0.815 in ST/SB prediction. Two ternary formulations (Hydrocortisone/β-CD/HPMC and dovitinib/γ-CD/CMC) were used to validate the prediction model. Molecular modeling results showed that HPMC not only warped around hydrocortisone but also prevented CD molecules from self-aggregation to increase solubility. In conclusion, a prediction model for the ternary CD formulations was successfully developed, which will significantly accelerate the formulation screening process to benefit the formulation development of water-insoluble drugs.
Collapse
Affiliation(s)
- Junjun Li
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China.
| | - Hanlu Gao
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China.
| | - Zhuyifan Ye
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China.
| | - Jiayin Deng
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China.
| | - Defang Ouyang
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China.
| |
Collapse
|
6
|
Why 2,6-di-methyl-β-cyclodextrin can encapsulate OH-substituted naphthalenes better than β-cyclodextrin: Binding pose, non-covalent interaction and solvent effect. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113496] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
7
|
Dai X, Chen L, Liao Y, Sheng M, Qu Q, Shi Y, Shi X. Formulation design and mechanism study of hydrogel based on computational pharmaceutics theories. J Mol Graph Model 2021; 110:108051. [PMID: 34715467 DOI: 10.1016/j.jmgm.2021.108051] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2021] [Revised: 10/10/2021] [Accepted: 10/11/2021] [Indexed: 12/01/2022]
Abstract
Formulation design and mechanism study of the drug delivery system (DDS) is an important but difficult subject in pharmaceutical research. The study of formulation factors is the most time- and labor-consuming work of formulation design. In this paper, a multiscale computational pharmaceutics strategy was developed to guide the systematic study of formulation factors of a typical polymer-based DDS, hydrogel, and further to guide the formulation design. According to the strategy, the combination of solubility parameter (δ) and diffusion coefficient (D) calculated by the AA-MD simulation was suggested as the general evaluation method for the matrix screening of the hydrogels at the pre-formulation stage. At the formulation design stage, the CG-MD simulation method was suggested to predict the morphology and drug-releasing behavior of the hydrogels under different formulation factors. The influence mechanism can be explained by the combination of multiple parameters, such as the microstructure diagram, the radius of gyration (Rg), the radial distribution function (RDF), and the free diffusion volume (Vdiffusion). The simulation results are in good agreement with the in vitro release experiment, indicating that the strategy has good applicability.
Collapse
Affiliation(s)
- Xingxing Dai
- Beijing University of Chinese Medicine, No. 11 of North 3rd Ring East Road, Chaoyang District, Beijing, 100029, China; Key Laboratory for Production Process Control and Quality Evaluation of Traditional Chinese Medicine, Beijing Municipal Science & Technology Commission, Beijing, 100029, China
| | - Liping Chen
- Beijing University of Chinese Medicine, No. 11 of North 3rd Ring East Road, Chaoyang District, Beijing, 100029, China; Key Laboratory for Production Process Control and Quality Evaluation of Traditional Chinese Medicine, Beijing Municipal Science & Technology Commission, Beijing, 100029, China
| | - Yuyao Liao
- Beijing University of Chinese Medicine, No. 11 of North 3rd Ring East Road, Chaoyang District, Beijing, 100029, China; Key Laboratory for Production Process Control and Quality Evaluation of Traditional Chinese Medicine, Beijing Municipal Science & Technology Commission, Beijing, 100029, China
| | - Mengke Sheng
- Beijing University of Chinese Medicine, No. 11 of North 3rd Ring East Road, Chaoyang District, Beijing, 100029, China; Key Laboratory for Production Process Control and Quality Evaluation of Traditional Chinese Medicine, Beijing Municipal Science & Technology Commission, Beijing, 100029, China
| | - Qingsong Qu
- Beijing University of Chinese Medicine, No. 11 of North 3rd Ring East Road, Chaoyang District, Beijing, 100029, China; Key Laboratory for Production Process Control and Quality Evaluation of Traditional Chinese Medicine, Beijing Municipal Science & Technology Commission, Beijing, 100029, China
| | - Yanshuang Shi
- Beijing University of Chinese Medicine, No. 11 of North 3rd Ring East Road, Chaoyang District, Beijing, 100029, China; Key Laboratory for Production Process Control and Quality Evaluation of Traditional Chinese Medicine, Beijing Municipal Science & Technology Commission, Beijing, 100029, China
| | - Xinyuan Shi
- Beijing University of Chinese Medicine, No. 11 of North 3rd Ring East Road, Chaoyang District, Beijing, 100029, China; Key Laboratory for Production Process Control and Quality Evaluation of Traditional Chinese Medicine, Beijing Municipal Science & Technology Commission, Beijing, 100029, China.
| |
Collapse
|
8
|
Mazurek AH, Szeleszczuk Ł, Gubica T. Application of Molecular Dynamics Simulations in the Analysis of Cyclodextrin Complexes. Int J Mol Sci 2021; 22:9422. [PMID: 34502331 PMCID: PMC8431145 DOI: 10.3390/ijms22179422] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2021] [Revised: 08/27/2021] [Accepted: 08/27/2021] [Indexed: 11/17/2022] Open
Abstract
Cyclodextrins (CDs) are highly respected for their ability to form inclusion complexes via host-guest noncovalent interactions and, thus, ensofance other molecular properties. Various molecular modeling methods have found their applications in the analysis of those complexes. However, as showed in this review, molecular dynamics (MD) simulations could provide the information unobtainable by any other means. It is therefore not surprising that published works on MD simulations used in this field have rapidly increased since the early 2010s. This review provides an overview of the successful applications of MD simulations in the studies on CD complexes. Information that is crucial for MD simulations, such as application of force fields, the length of the simulation, or solvent treatment method, are thoroughly discussed. Therefore, this work can serve as a guide to properly set up such calculations and analyze their results.
Collapse
Affiliation(s)
- Anna Helena Mazurek
- Department of Physical Chemistry, Chair of Physical Pharmacy and Bioanalysis, Faculty of Pharmacy, Doctoral School, Medical University of Warsaw, Banacha 1 Street, 02-093 Warsaw, Poland;
| | - Łukasz Szeleszczuk
- Department of Physical Chemistry, Chair of Physical Pharmacy and Bioanalysis, Faculty of Pharmacy, Medical University of Warsaw, Banacha 1 Street, 02-093 Warsaw, Poland;
| | - Tomasz Gubica
- Department of Physical Chemistry, Chair of Physical Pharmacy and Bioanalysis, Faculty of Pharmacy, Medical University of Warsaw, Banacha 1 Street, 02-093 Warsaw, Poland;
| |
Collapse
|
9
|
Shi JH, Lin ZY, Kou SB, Wang BL, Jiang SL. Enantioseparation of mandelic acid and substituted derivatives by high-performance liquid chromatography with hydroxypropyl-β-cyclodextrin as chiral mobile additive and evaluation of inclusion complexes by molecular dynamics. Chirality 2021; 33:675-684. [PMID: 34390038 DOI: 10.1002/chir.23348] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2021] [Revised: 07/15/2021] [Accepted: 07/21/2021] [Indexed: 12/21/2022]
Abstract
The enantioseparation and resolution mechanism of mandelic acid (MA), 4-methoxymandelic acid (MMA), and 4-propoxymandelic acid (PMA) were investigated by reversed-phase high-performance liquid chromatography (HPLC) with 2-hydroxypropyl-β-cyclodextrin (HP-β-CD) as a chiral mobile-phase additive and molecular dynamics simulation. The suitable chromatographic conditions for the enantioseparation of MA, MMA, and PMA were obtained. Under the selected chromatographic conditions, these enantiomers could achieve baseline separation. The results of thermodynamic parameter analysis revealed that the main driven forces for the enantioseparation of MA, MMA, and PMA could be van der Waals forces and hydrogen-bonding interactions and the chromatographic retention of these chiral compounds was an enthalpy-driven process. The results of the molecular simulation revealed that their chiral resolution mechanism on HP-β-CD was responsible for the formation of inclusion complexes of enantiomers with HP-β-CD with different conformations and binding energies. And the binding energy of HP-β-CD with (S)-isomer was larger than that with (R)-isomer, which is consistent with the experimental results of the first elution of (S)-isomer. Additionally, it is also confirmed that the interaction energies included the van der Waals energy (∆Evdw ), electrostatic energy (∆Eelec ), polar solvation energy, and SASA energy (∆Esasa ), and the separation factor (α) was closely connected with the disparity in the binding energies of optical isomers and HP-β-CD complexes. Meanwhile, from molecular dynamics simulation, it can be found that the ∆(∆Ebinding ), (∆(∆Ebinding ) = ∆Ebinding,R - ∆Ebinding,S ) value was in order of MA-HP-β-CD complex > MMA-HP-β-CD complex > PMA-HP-β-CD complex, which was consistent with the order of Δ(ΔG) values obtained from van't Hoff plot. This indicated that the molecular dynamics simulation has predictive function for chiral resolution.
Collapse
Affiliation(s)
- Jie-Hua Shi
- College of Pharmaceutical Science, Zhejiang University of Technology, Hangzhou, P.R. China
| | - Zhen-Yi Lin
- College of Pharmaceutical Science, Zhejiang University of Technology, Hangzhou, P.R. China
| | - Song-Bo Kou
- College of Pharmaceutical Science, Zhejiang University of Technology, Hangzhou, P.R. China
| | - Bao-Li Wang
- College of Pharmaceutical Science, Zhejiang University of Technology, Hangzhou, P.R. China
| | - Shao-Liang Jiang
- College of Pharmaceutical Science, Zhejiang University of Technology, Hangzhou, P.R. China
| |
Collapse
|
10
|
Zhao L, Sun Q, Pu H, Tang P, Liu Y, Li M, Ren X, Li H. Experimental and computer simulation investigations of ethyl red with modified β-cyclodextrins: Inclusion mechanism and structure characterization. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137725] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
11
|
|