Prediction of Carbamazepine-Succinic Acid Co-Crystal Dissolution in Ethanolic Solution using a Computational Molecular Dynamic Simulation Technique

  • Authors

    • Nornizar Anuar
    • Noor Sazmira Pauzi
    • Siti Nurul ‘Ain Yusop
    • Nik Salwani Md Azmi
    • Syarifah Abd Rahim
    • Muhammad Fitri Othman
    2018-11-27
    https://doi.org/10.14419/ijet.v7i4.18.21884
  • Carbamazepine, co-crystallization, dissolution, molecular dynamic simulation, morphology, solubility, succinic acid, surface chemistry.

  • Abstract

    Carbamazepine (CBZ) is an anti-epileptic Class II drug according to the biopharmaceutical classification system and it forms a co-crystal with succinic acid (SA). The physicochemical properties of co-crystal are of interest since they control how the co-crystal behaves in different environment. In this study, the morphology of carbamazepine-succinic acid (CBZ-SA) co-crystal was predicted by applying the atomic charge calculated from MOPAC and DMol3 while the dissolution behavior of CBZ-SA in ethanol  was investigated using dynamic simulation by considering the transport properties of both co-crystal and solvent. The predicted CBZ-SA morphology shows a plate-like shape with the main crystal facet of (1 0 -1), (1 0 1), (1 1 0), (0 1 1), (0 1 0), (2 0 0) and (0 0 2). Attachment energy calculation shows that facet (1 0 -1) is the slowest growing facet having the lowest attachment energy of -27.91 kcal/mol, while facet (0 1 1) is the fastest growing facet with the highest value of attachment energy of -96.74 kcal/mol. The dissolution behavior of CBZ-SA in ethanol was also assessed through the mean square displacement (MSD) and diffusion coefficient calculation, and the result shows that the dissolution of co-crystal first occurs at facet (1 0 1) and last to dissolve at facet (1 1 0). Because the molecular configuration of the molecules at facet (1 0 1) leaves a large gap between molecules that makes it easier for the solvent molecules to penetrate the surface layer and the hydrophilic part of CBZ and succinic acid molecules is exposed to the surface forming hydrogen bond. Calculation of binding energy shows that the interaction of ethanol was more favorable at facet (1 1 0) with the value of -745.631 kcal/mol. The analysis carried out in this work showed that there is a promising relationship between RDF and diffusion coefficient in predicting the diffused molecules from the surface of the facets, as both methods predicted the facets (0 0 2), (2 0 0) and (1 0 1) as among the top three facets dissolved into the ethanol  environment. Both MSD and diffusion coefficient also predict that facets (0 1 0) and (1 1 0) are among the last facets dissolved. However, the binding energy calculation does not show any apparent relationship to neither RDF nor diffusion coefficient.

     

     

  • References

    1. [1] Ali S, Velaga S, & Jouyban A (2014), Solubility of carbamazepine, nicotinamide and carbamazepine–nicotinamide cocrystal in ethanol–water mixtures. Fluid Phase Equilibria 363, 97–105.

      [2] Tian F, Sandler N, Aaltonen J, Lang C, Saville DJ, Gordon KC, Strachan CJ, Rantanen J & Rades T (2017), Influence of polymorphic form, morphology, and excipient interactions on the dissolution of carbamazepine compacts. Journal of Pharmaceutical Sciences 96, 584 – 594.

      [3] Vadlamudhi MK & Dhanaraj S (2017), Significance of excipients to enhance the bioavailability of poorly water-soluble drugs in oral solid dosage forms: a review. IOP Conference Series: Materials Scince and Engineering, 263, 022023

      [4] Cascone, S (2017), Modeling and comparison of release profiles: effect of the dissolution method. European Journal of Pharmaceutical Sciences 106, 352–361

      [5] Ab Rahman F, Abd Rahim S, Tan CC, Low SH & Ramle NA (2017), Carbamazepine-fumaric acid and carbamazepine-succinic acid co-crystal screening using solution based method. International Journal of Chemical Engineering and Applications 8, 136 – 140

      [6] Liu M, Hong C, Yao Y, Shen H, Ji G, Li G, & Xie Y (2016), Development of a pharmaceutical cocrystal with solution crystallization technology: Preparation, characterization, and evaluation of myricetin-proline cocrystals. European Journal of Pharmaceutics and Biopharmaceutics 107, 151–159

      [7] Rahim SA, Rahman FA, Nasir ENEM & Ramle NA (2015), Carbamazepine co-crystal screening with dicarboxylic acids co-crystal formers. International Journal of Chemical and Molecular Engineering 9(5), 442-445.

      [8] Fuliaş A, Vlase G, Vlase T, Şuta L-M, Şoica C & Ledeţi I (2015), Screening and characterization of cocrystal formation between carbamazepine and succinic acid. Journal of Thermal Analysis and Calorimetry 121, 1081 – 1086

      [9] Ullah M, Hussain I & Sun CC (2016), The development of carbamazepine-succinic acid cocrystal tablet formulations with improved in vitro and in vivo performance. Drug Development and industrial Pharmacy 42, 969-976

      [10] Gao Y & Olsen KW (2013), Molecular dynamics of drug crystal dissolution: simulation of acetaminophen form I in water. Molecular Pharmaceutics 10, 905 – 917

      [11] Chan T & Ouyang D (2018), Investigating the molecular dissolution process of binary solid dispersions by molecular dynamics simulations. Asian Journal of Pharmaceutical Sciences 13, 248 – 254

      [12] Song L, Chen L, Wang J, Chen F & Lan G (2017), Prediction of crystal morphology of 3, 4-Dinitro-1H-pyrazole (DNP) in different solvents. Journal of Molecular Graphics and Modelling 75, 62 – 70

      [13] Shi W, Chu Y, Xia M, Lei W & Wang F (2016), Crystal morphology prediction of 1,3,3-trinitroazetidine in ethanol solvent by molecular dynamics simulation. Journal of Molecular Graphics and Modelling 64, 94 – 100

      [14] Rohl AL (2003), Computer prediction of crystal morphology. Current Opinion in Solid State and Materials Science 7, 21 – 26

      [15] Shim HM, Kim JK, Kim HS & Koo KK (2014), Molecular dynamics simulation on nucleation of ammonium perchlorate from an aqueous solution. Crystal Growth & Design 14, 5897-5903

      [16] Zeng, JP., Bai Y-S., Chen S & Ma C-A (2013), Molecular dynamics simulation of diffusion of nitrobenzene in 3-methylimidazolium hexa fluoro phosphate ionic liquids. Journal of Molecular Liquids 183, 1-7

      [17] Sun Q, Liu Q, Xiao J (2014), Molecular dynamics simulation of interface interactions and mechanical properties of CL-20/HMX cocrystal and its based PBXs. Acta Chim. Sinica 72, 1036 –1042.

      [18] Clydesdale G, Roberts KJ, Telfer GB & Grant DJW (1997), Modelling the crystal morphology of α-lactose monohydrate. Journal of Pharmaceutical Sciences 86, 135 – 141

      [19] Schmidt C, Yϋrϋdϋ C, Wachsmuth A & Ulrich J (2011), modeling the morphology of benzoic acid crystals grown from aqueous solution. Crystal Engineering Communication 13, 1159 – 1169

      [20] Coombes DS, Richard C, Catlow A, Gale JD, Rohl AL & Price SL (2005), Calculation of attachment energies and relative volume growth rates as an aid to polymorph prediction. Crystal Growth & Design 5, 879 – 885

      [21] Berry DJ & Steed JW (2017), Pharmaceutical cocrystals, salts and multicomponent systems; intermolecular interactions and property-based design. Advanced Drug Delivery Reviews 117, 3 – 24.

      [22] Zeng J, Aimin W, Xuedong G, Jingwen C, Song C & Feng X (2012), Molecular dynamics simulation of diffusion of vitamin c in water solution. Chinese Journal of Chemistry 30, 115 – 120

      [23] Gholizadeh R., Wang Y & Yu Y-X (2017), Molecular dynamics simulations of stability at the early stages of silica materials preparation. Journal of Molecular Structure 1138, 198 – 207

      [24] Shi W, Xia M, Lei W & Wang F (2014), Solvent effect on the crystal morphology of 2,6- diamino-3,5-dinitropyridine-1-oxide: a molecular dynamics simulation study. Journal of Molecular Graphics and Modelling 50, 71–77

  • Downloads

  • How to Cite

    Anuar, N., Sazmira Pauzi, N., Nurul ‘Ain Yusop, S., Salwani Md Azmi, N., Abd Rahim, S., & Fitri Othman, M. (2018). Prediction of Carbamazepine-Succinic Acid Co-Crystal Dissolution in Ethanolic Solution using a Computational Molecular Dynamic Simulation Technique. International Journal of Engineering & Technology, 7(4.18), 122-128. https://doi.org/10.14419/ijet.v7i4.18.21884

    Received date: 2018-11-27

    Accepted date: 2018-11-27

    Published date: 2018-11-27