A Study Showcasing Crystal Medley of Metformin Free Base: Virtual and Experimental Insights

Objective: The present work reports the existence of polymorphism in metformin free base applying virtual and experimental approach. Methods: The experimental screening of various polymorphs of metformin free base has been coupled with virtual approach resulting in isolation of three distinct crystal phases from different solvents/ or solvent mixtures. The virtual screening of polymorphs was performed using Polymorph Predictor module of BIOVIA Material Studio (MS) software. A state-of-the-art approach was undertaken utilizing integrated molecular modeling and Monte Carlo simulation technology. The sophisticated tools were used to characterize the experimentally isolated forms such as Differential Scanning Calorimetry (DSC), Fourier Transform Infra-Red Spectroscopy (FTIR), Powder X Ray Diffraction (PXRD) and optical microscopy. Lattice energy landscape was examined to look for these observed forms. Crystal morphology study was conducted to analyse Morphologically Important (MI) facets. Results: A list of potential polymorphs of the drug molecule was generated with varied lattice constants using MS. Crystal energy landscape obtained from the study provided systematic energy ranking to the polymorphs. DSC of these forms exhibited endothermic peaks at 134.01, 111.24 and 79.84°C respectively indicating them to be different forms of metformin base. Structure determination from powder patterns showed that although all the forms, form I, II and III crystallize out in same triclinic space group, P-1 but they possessed different crystallographic parameters and hydrogen bonding patterns. The predicted results were quite complacent with experimental observations with respect to lattice parameters. The three observed forms have been found as local minima in the lattice energy landscape of potential polymorphs. Conclusion: The study explained the practical relevance of Crystal Structure Prediction (CSP) in studying the potential polymorphs of a drug molecule by successfully predicting various polymorphs of metformin free base. Out of all predicted polymorphs, three polymorphs were actually isolated experimentally. The information generated about the existence of multiple forms of metformin free base provides an opportunity to select and scale-up the desired crystal form well suited on bioavailability and stability grounds. The process can be escalated further to preclinical study, to emerge as an invaluable technology in future in increasing the efficiency of drug development process.


Introduction
Multifariousness in crystal forms of a drug molecule has become interesting topic of discussion in the field of pharmaceuticals. This inherent tendency of the compounds to exist in various crystal forms is called polymorphism. In general polymorphism is classified into two major types, viz., monotropism and enantiotropism. This general behavior is based on their temperature linked stability. In the monotropic system, only one polymorph is stable below the melting point of the drug substance and in enantiotropic system, the substance exists in different crystal forms which are interconvertible based on their transition temperature (the temperature at which two polymorphic forms exhibit same free energy). Besides this polymorphism can be categorized into two types based on their arrangements in solid state, i.e., packing polymorphism and conformational polymorphism. In the former the differences are due to molecular packing in the unit cell and in the later case, the substances differ due to the presence of different conformers in their solid-state arrangements. Moreover, literature is flooded with reports and reviews on this interesting phenomenon of polymorphism [1][2][3][4][5][6][7][8][9][10][11][12]. From regulatory perspectives, a new polymorphic form is designated as a "pharmaceutical alternative" to the original form. Exploring various polymorphic forms of a drug molecule has become a critical part of research and development as phase transformations have a major impact on product quality and stability [1]. It is also quite useful from intellectual proprietary benefits for the pharmaceutical industry [3,4]. It is very well stated that different solid forms of a drug substance exhibit different physicochemical and biopharmaceutical properties such as solubility, dissolution, stability, manufacturability including flowability, tablet compaction, hygroscopicity, etc. This will ultimately affect the bioavailability and bioequivalence of that substance. Furthermore, different solid forms of the active pharmaceutical ingredients can possess variable organoleptic properties such as taste, color etc., [5][6][7][8][9][10][11]. Thus, in-depth knowledge of all the crystal forms should be acquired in advance in order to ensure the efficient and successful development of a robust product [12]. The present study focuses on the screening of various crystal forms of metformin free base. Metformin hydrochloride is a first line treatment of Type 2 diabetes [13]. Additionally, modern metformin has received a new attention due to its anti-ageing properties by dramatically increasing the health span of individuals making people healthy, active and alert years longer. This miracle pill favorably modulates the major biological outcomes associated with longevity such as inflammation, stress defence, cellular survival, autophagy and protein synthesis [14][15][16]. It has countless of other health benefits such as improvement of cognitive function, cardioprotection, anticancer properties, etc., [17]. The present study thus selected this important molecule to play with its solid-state properties. The study focuses on the virtual and experimental screening of various crystal forms of metformin free base. Metformin is a strong base and belongs to BCS class III drug category [17] and is commercialized in its hydrochloride salt form. Unfortunately, its unpleasant taste and high dosing requirements could lead to compliance problems. Moreover, hydrochloride salts could present problems in handling especially in high humidity conditions and has high compaction susceptibility. This salt during manufacturing also generates corrosive action on tablet tooling [18][19][20][21][22][23][24][25]. Furthermore, controlled release formulation of this salt form is difficult to prepare due to its very high aqueous solubility parameter [26]. Thus, in order to deal with the disadvantages of most frequently used hydrochloride salt form of drug, the present work was initiated to explore the various polymorphs of metformin free base with the aim of selection of suitable crystal form for future development. The literature survey has revealed that not much work has been done on polymorphic concepts of metformin free base. The present study is a well-designed screening strategy which is focused on development timelines for the product in pharmaceutical sector. Biovia Material Studio in silico strategy involving Crystal Structure Prediction (CSP) is combined with experimental screening of polymorphic forms. All the observed forms were characterized using sophisticated analytical techniques. The approach followed in the present work is quite efficient in detecting and interpreting various crystal forms of metformin base. The Crystal Structure Determination (CSD) using powder X-ray diffraction scans of various observed forms is then undertaken using Biovia Material Studio. The final outcome of CSP, crystal energy landscape is looked upon for the experimentally observed crystal forms. This gives the idea about the lattice energy associated with these forms which in turn expresses the stability of these solid forms. Furthermore, morphology simulation of observed forms of metformin free base was performed. The final outcome, crystal habit generated from morphology study has assisted CSP in controlling the existence of various polymorphs of drug substance.

Virtual screening of polymorphs
The Crystal Structure Prediction (CSP) was implemented using Polymorph module of Biovia Material Studio (MS) software suite, version 7.0. The initial step of energy minimization involves geometry optimization using Dmol3 program and COMPASS forcefield technology of MS based on quantum mechanical methods, starting from the input molecular structure of metformin. It creates an optimized input structure for polymorph prediction. The hunt for the polymorphs was executed within most commonly occurring space groups i.e., P-1, Pbca, P21/c, P212121, P2/c, C2/c and P21 [27] and was repeated thrice till all predicted crystal structures converges to equal energy minima. The Monte Carlo method of MS involves a systematic and stepwise approach towards polymorph prediction involving various steps such as packing, pre-clustering, optimization and clustering. The eventual outcome of polymorph prediction is crystal lattice energy landscape (energy versus density plot). It ranks the crystal structures based upon their lattice energy and density parameters and generates a list of lattice parameters of predicted structures within selected space groups [28][29][30]. This is further looked up for local minima (metastable crystal structures), global minima and thermodynamically most stable form.

Experimental screening of polymorphs
Metformin free base was isolated from commercial sample of Metformin HCl (procured from Thermo Fischer Scientific) as stated in United States patent US 6,031,004 and US 8,076,377 [25,26]. This free base form of metformin was further explored for its polymorphic tendency and success was achieved with slow evaporation method. Different solid forms were successfully isolated by slow evaporation technique using various solvents/ or solvent mixture. The saturated solutions of metformin base in various solvents and solvent mixtures were heated below the boiling point of these solvent/ or solvent mixtures and then kept for evaporation at room temperature. Form I crystallized out from benzaldehyde in 9 days. Form II from 2-propanol appeared in 2 days and form III crystallized out from a solvent mixture of cyclohexane and toluene in 21 days. Crystals observed were filtered and stored in desiccator under controlled humid conditions till further characterization.

Characterization of polymorphs
Thermal analysis: DSC thermograms of isolated forms were interpreted using Q20 DSC instrument by placing the sample (3-5 mg) into sealed aluminium pans and scanning (50-300°; 10°C/min.; nitrogen purge of 50 mL/min.) Data was interpreted using Universal analysis 2000 software (TA Q series advantage software).

Fourier Transform Infra-Red Spectroscopy (FTIR):
The spectral data of the observed forms of metformin free base pressed into pellet (containing 2-4 mg sample dispersed in approx. 20 mg KBr) was collected over the range of 450-4000cm -1 using Perkin-Elmer RX-1 FTIR spectrophotometer (UK). For metformin, the free base spectra exhibited a notable band at 1670 cm -1 (C=N imine stretch).

Powder X-Ray Diffraction (PXRD):
The samples of all the collected forms (about 200 mg each) were loaded on X′PERT-PRO diffractometer using Cu-Kα radiation (1.54°C) and scanned at 2θ values ranging from 3.5° to 49° with a step size of 0.0170°, scan step time set to 25.1996 s and slit size adjusted to 0.48°C. The diffractograms were recorded and analyzed using XPERT high score software.
Optical Microscopy: The overall macroscopic crystal habit of the different polymorphs was studied by observing optical micrographs of all the samples taken using a Leica DM3000 upright optical microscope with polarizer at an image resolution of 25X.

Crystal structure determination
This was done by collecting the PXRD patterns of observed forms and subjecting them to Reflex Plus module of Material Studio. It was performed following four basic steps i.e., indexing, Pawley fitting/ refinement, structure solution and Rietveld refinement [31]. X-cell (Newmann, 2003) indexing program was used for peak indexing to obtain various solutions of crystal lattices. The unit cell possessing highest Figure of Merit (FOM) was then subjected to refinement for attaining accurate lattice constants and lattice parameters. The optimization of unit cell was done from 10 cycles followed with a search for space groups. The optimized molecular structure was imported into the refined empty unit cell. The structure solution was finally carried out by exposing this whole assembly to Powder Solve program of Reflex Plus module of MS [31,32]. Rietveld refinement of the solved structure was executed to generate R wp (weighted Rietveld parameter) values which further confirms the similarity between a simulated and an experimental diffraction pattern.
The final optimization of the structure was achieved using Forcite module of MS.

Crystal morphology simulation
The morphology prediction of the observed polymorphic modifications of was carried out using the Morphology module of MS [33,34]. The maximum h,k,l values (Miller indices) for the faces were set to 5 each with an upper limit of 100 growing faces. The minimum inter-planar spacing (d hkl ) for face list generation was adjusted to 1.0 A°. Growth morphology algorithm working on the hypothesis that crystal face growth rate is proportional to its attachment energy [35][36][37] was then applied to observed crystal structures. The outcome generated is a list of growth faces along with their attachment energy (E att ). The habit properties including aspect ratio was calculated.

Solubility study
The equilibrium solubility study of observed forms of metformin free base was carried out by shaking an excess of drug (approx. 20 mg) in 10 ml of distilled water using MSW-275 (Macro scientific works, New Delhi) water bath shaker preset at 37°C and 200 rpm for 24 hrs. The aliquots were passed through 0.45µm membrane filter and analyzed spectrophotometrically at λ max = at 232 nm, E 1% 1 cm = 16,398. Concentration (C) of various forms was determined using the formula:

Virtual analysis from Crystal Structure Prediction (CSP)
The virtual analysis of polymorph screening using Polymorph Predictor module of Biovia Material Studio was carried out starting from molecular structure of metformin ( Figure 1). The simulation result was retrieved in form of a table that contains data on space groups, lattice energy and lattice constants corresponding to various predicted polymorphs. Out of the various generated structures, 60 unique sets were extracted corresponding based on low energy and high-density criteria (Table 1).

Characterization
The practically isolated three forms of metformin free base were analysed using sophisticated tools such as DSC, FTIR spectroscopy, PXRD and optical microscopy.

FTIR Spectroscopy
Variations in the characteristic spectral frequencies of observed forms further demonstrates the existence of polymorphism in metformin free base. The information obtained from FTIR spectra ( Table 3, Figure 3(a-c)) of form I, II and III indicated the differences in the intermolecular interactions within these modifications.

PXRD Study
The various peaks observed in PXRD scans ( Figure 4) of three forms of metformin are given in table 4. Sharp differences in their 2θ values clearly characterize them as different polymorphic forms.   Expt.

Morphology study
Crystal morphology and habit of the developed forms of metformin was studied using growth morphology algorithm of morphology module of MS starting from their respective crystal structures. It is based upon the hypothesis that the growth rate of a crystal face is proportional to its attachment energy (E att ). E att is a habit controlling factor which is defined as energy per mol of the molecule liberated when a new layer having thickness d hkl (hkl indicates the Miller indices assigning a set of faces that are equivalent by the symmetry of the crystal) is joined to the surface of the crystal). Thus, the crystal faces with lowest attachment energies (in terms of magnitude) are the slowest growing and, considered as morphologically most important facets [36,37]. Now, if the attachment energy of the slowest developing or morphologically important face is taken

Solubility study
The solubility study performed in distilled water at 37°C (Table 7) indicated that the form with highest lattice energy or stability (form III) was least soluble (271 ± 2.4). The solubility profile was thus found in accord with their lattice energy data.

Form Solubility in distilled water at 37°C (mg/ ml) ± SD
Form I 329 ± 1.2 Form II 307 ± 2.5 Form III 201 ± 3.2 Table 7: Solubility data of observed forms of metformin.

Discussion
The results of the prediction study compiled in a tabular form revealed the significance of CSP in acquiring the data on possible polymorphs of a drug molecule. The knowledge generated from lattice energy landscape would give an idea of prevalence of multiple solid forms of a drug molecule based on crystal packing space featuring low energy and high-density crystal structures. Polymorph Predictor module and userfriendly algorithms of Biovia MS has shaped CSP as a powerful tool to predict the crystal structures of a drug substance.
The sharp endothermic peaks of observed forms at different positions in DSC curves clearly gives a clue about the existence of different solid form of drug.
The shift in FTIR bands depicting changes in vibrational frequencies support the existence of different solid phases due to variable hydrogen bonding patterns. This technique helps in identifying functional groups involved in hydrogen bonding.
The unique differences in PXRD patterns of three forms has confirmed the prevalence of polymorphism in metformin free base. The marked appearance of new peaks and disappearance of certain prominent peaks supports the evolution of new solid phases.   The simulated and experimental PXRD scans complement each other as interpreted from R wp value which further proclaims the neatness of crystal structure determination using powder technology (Figure 12(a-c)). R wp is a basic measure of similarity between a simulated and an experimental diffraction pattern of each crystal structure and is found out to be below 15%.   . 13 .
Structures corresponding to predicted and observed polymorphs of metformin free base found in crystal energy landscape are featured in table 5. The predicted and experimental results are found to be quite complacent with each other. It is clearly seen that developed crystal structures exists in the list of simulated results within an energy difference of 7.5 kcal/mol from the global minima. The simulated crystal structure nos. 6, 11 and 20 complement well with experimentally observed forms. The lattice energy diagram was generated from the predicted structures ( Figure 13) depicting colorfully highlighted experimentally matched forms (mentioned under legend). The interpretation of crystal energy landscape has very well positioned the experimentally observed solid forms of metformin as local minima within an energy difference of approximately 7.0 kcal/mol from the global minimum. This further indicates the existence of multiple metastable forms of the drug.
On interpretation of lattice energy landscape, form III occupies sixth position in the list of predicted structures generated form CSP runs and is nearest to the global minimum. It is high density form having approximately 2.5 kcal/mol higher energy than that of global minimum which further confirms it to be the most stable crystalline form. Form II exists as eleventh local minima in the lattice energy landscape. The relatively low lattice energy and highest density value provided it good stability as compared to other possible forms. Form I is found as the twentieth local minimum predicted in the study.
The morphology prediction has well demonstrated that lattice energy model has ranked the observed crystal forms in correct manner with respect to their stability and observing tendency.
The morphology prediction of various forms of metformin and their ranking in the lattice energy landscape are in good agreement with each other as far as the stability is concerned (Table 1). It was concluded after careful examination that form I is most stable having highest E att of its M.I. facet (100), followed by form II (facet (010)), and form III (facet (1-10)). This also expresses the growing speed of crystal and their observing tendency during experimentation.
It is likely to understand that many low energies predicted structures have not assimilated as polymorphs as they may not be favored kinetically and thus may not be appear experimental crystallization temperatures. The higher aspect ratios (>7) can also be the other parameter for not observing these hypothetical low energy forms including the global minimum. Thus, morphology prediction further assists the CSP by filtering out some of the hypothetical structures based on the consideration of growth rates, using the attachmentenergy model. The knowledge of aspect ratios indicates that crystals growing as thin plates so slowly that they are unlikely to be found in competition with equally thermodynamically stable structures.
Visual examination by optical microscopy, generates information on the crystal topography and habit suggesting the prevalence of different crystal phases of metformin free base.
Solubility study also demonstrates the existence of multiple forms of the drug.

Conclusion
The findings illustrate the implications of CSP in understanding the concepts of polymorphism in metformin free base. Out of the three isolated forms, form I of metformin free base has been found to be the most stable which can be further considered for development. Thus, the varied options of polymorphic forms with desired characteristics has been provided by lattice energy landscape of drug molecule. The efficiency of experimental work can be facilitated by utilizing a CSP which can reduce the amount of time and costs of the experimental part of research and development. Metformin free base has a strong tendency to exist in variable crystal forms as shown by CSP and illustrated in the crystal energy landscape. Onset of multiple phases of drug molecule is depicted from thermograms and PXRD scans which is also backed by spectral analysis using FTIR. Crystal structure determination from PXRD patterns confirmed the development of three forms of metformin. The variation in solubility values further helps in selection of suitable form for isolation. The study insists in scrutinizing other parameters such as solvent effects on crystallization for more detailed vision on the polymorphic aspects of a drug molecule. The study plays a constructive role in accumulating the structurally diversified data of a drug molecule which further permits to select the desired crystal form for future development.