Engineering Cocrystals of Poorly Water-Soluble Drugs to Enhance Dissolution in Aqueous Medium

Biopharmaceutics Classification System (BCS) Class II and IV drugs suffer from poor aqueous solubility and hence low bioavailability. Most of these drugs are hydrophobic and cannot be developed into a pharmaceutical formulation due to their poor aqueous solubility. One of the ways to enhance the aqueous solubility of poorlywater-soluble drugs is to use the principles of crystal engineering to formulate cocrystals of these molecules with water-soluble molecules (which are generally called coformers). Many researchers have shown that the cocrystals significantly enhance the aqueous solubility of poorly water-soluble drugs. In this review, we present a consolidated account of reports available in the literature related to the cocrystallization of poorly water-soluble drugs. The current practice to formulate new drug cocrystals with enhanced solubility involves a lot of empiricism. Therefore, in this work, attempts have been made to understand a general framework involved in successful (and unsuccessful) cocrystallization events which can yield different solid forms such as cocrystals, cocrystal polymorphs, cocrystal hydrates/solvates, salts, coamorphous solids, eutectics and solid solutions. The rationale behind screening suitable coformers for cocrystallization has been explained based on the rules of five i.e., hydrogen bonding, halogen bonding (and in general non-covalent bonding), length of carbon chain, molecular recognition points and coformer aqueous solubility. Different techniques to screen coformers for effective cocrystallization and methods to synthesize cocrystals have been discussed. Recent advances in technologies for continuous and solvent-free production of cocrystals have also been discussed. Furthermore, mechanisms involved in solubilization of these solid forms and the parameters influencing dissolution and stability of specific solid forms have been discussed. Overall, this review provides a consolidated account of the rationale for design of cocrystals, past efforts, recent developments and future perspectives for cocrystallization research which will be extremely useful for researchers working in pharmaceutical formulation development.


Introduction
Biopharmaceutics Classification System (BCS) of drugs classifies drugs into four major categories ( Figure 1) based on their solubility and permeability behavior [1]. BCS Class II and Class IV drugs suffer from poor aqueous solubility. Poor aqueous solubility of hydrophobic drugs can result in poor absorption, low bioavailability and poses challenges for drug development process [2]. Enhancing bioavailability of poorly water-soluble BCS class II and BCS class IV drugs therefore becomes necessary to improve drug's efficacy. In addition to BCS of drugs, Developability Classification System (DCS) of drugs also plays a significant role in determining the development of pharmaceutical formulations, especially the oral formulations based on its solubility in biorelevant media such as FaSSIF (Fast State Simulated Intestinal Fluid) and FeSSIF (Fed State Simulated Intestinal Fluid) rather than its solubility in buffers [3][4][5]. Very few reports are available in the literature where researchers have determined the dissolution rate of cocrystals of poorly water-soluble drugs in biorelevant media [6,7]. The studies illustrate that cocrystals exhibited enhanced dissolution rate in biorelevant media and buffer as well indicating that the DCS serves as a highly relevant tool in determining developability of cocrystals of poorly water-soluble APIs.
Enhancing aqueous solubility of poorly water-soluble drugs without compromising on stability is one of the major challenges faced by the pharmaceutical industries during drug discovery and development processes [8][9][10][11][12]. Crystal Engineering is a tool which can be used to tailor the physicochemical properties of Active Pharmaceutical Ingredients (APIs) [13] such as melting point, dissolution rate, aqueous solubility, refractive index, surface activity, habit, density, electrostatic, mechanical and optical properties [14]. Cocrystallization is one of the crystal engineering approaches adopted to prepare multicomponent pharmaceutical crystals to enhance the dissolution rates of poorly water-soluble APIs without affecting their intrinsic properties [15].
Tremendous research is now being conducted in cocrystallization. Discovery of new cocrystals is now considered as one of the patentable invention [16]. Cocrystal technology is being considered as an advanced technology used for improving the drug product by modifying the molecular conformations and intermolecular interactions [17]. Combination drug therapy (combining two drugs in a solid form and administering in a single dose) is considered to be another important advantage of the cocrystal technology [18,19]. However, the number of marketed cocrystal products available till date is very low [20,21]. Examples of such cocrystals available in market and approved by the Food and Drug Administration (FDA) are Entresto, Lexapro [20] and Depakote [21].
In this review, we provide a consolidated account of information available in the literature on cocrystal engineering. We have discussed different aspects of cocrystallization technology such as the synthesis techniques, type of intermolecular interactions involved in cocrystallization, different types of solids obtained during cocrystallization, techniques to characterize cocrystals, different types of cocrystals and dissolution of cocrystals in aqueous medium. In addition to BCS of drugs, Developability Classification System (DCS) of drugs also plays a significant role in determining the development of pharmaceutical formulations, especially the oral formulations based on its solubility in biorelevant media such as FaSSIF (Fast State Simulated Intestinal Fluid) and FeSSIF (Fed State Simulated Intestinal Fluid) rather than its solubility in buffers [3][4][5]. Very few reports are available in the literature where researchers have determined the dissolution rate of cocrystals of poorly water-soluble drugs in biorelevant media [6,7]. The studies illustrate that cocrystals exhibited enhanced dissolution rate in biorelevant media and buffer as well indicating that the DCS serves as a highly relevant tool in determining developability of cocrystals of poorly water-soluble APIs.
Enhancing aqueous solubility of poorly water-soluble drugs without compromising on stability is one of the major challenges faced by the pharmaceutical industries during drug discovery and development processes [8][9][10][11][12]. Crystal Engineering is a tool which can be used to tailor the physicochemical properties of Active Pharmaceutical Ingredients (APIs) [13] such as melting point, dissolution rate, aqueous solubility, refractive index, surface activity, habit, density, electrostatic, mechanical and optical properties [14]. Cocrystallization is one of the crystal engineering approaches adopted to prepare multicomponent pharmaceutical crystals to enhance the dissolution rates of poorly water-soluble APIs without affecting their intrinsic properties [15].
Tremendous research is now being conducted in cocrystallization. Discovery of new cocrystals is now considered as one of the patentable invention [16]. Cocrystal technology is being considered as an advanced technology used for improving the drug product by modifying the molecular conformations and intermolecular interactions [17]. Combination drug therapy (combining two drugs in a solid form and administering in a single dose) is considered to be another important advantage of the cocrystal technology [18,19]. However, the number of marketed cocrystal products available till date is very low [20,21]. Examples of such cocrystals available in market and approved by the Food and Drug Administration (FDA) are Entresto, Lexapro [20] and Depakote [21].
In this review, we provide a consolidated account of information available in the literature on cocrystal engineering. We have discussed different aspects of cocrystallization technology such as the synthesis techniques, type of intermolecular interactions involved in cocrystallization, different types of solids obtained during cocrystallization, techniques to characterize cocrystals, different types of cocrystals and dissolution of cocrystals in aqueous medium.

Cocrystals
Cocrystals are defined as crystalline materials with two or more different molecules (i.e., drug and coformers) in the same crystal lattice [47]. In pharmaceutical industries, cocrystals have gained a As new forms of pharmaceutical solids are being produced with different type of guest molecules (solvents/water/solids), assigning an exact nomenclature for each of the solid was a hotly debated topic for many years [44][45][46]. Inspired by the work of Aitipamula et al. [45], Grothe et al. [46] developed a straightforward system for classification of multicomponent solid forms. The classification system proposed by Grothe et al. [46] is shown in Figure 3.

Types of Solid Forms Obtained from Cocrystallization
While cocrystal formation is the main desired outcome, an unsuccessful cocrystallization event might result in formation of different solid forms such as eutectics [22][23][24][25][26], a drug polymorph [27,28], coamorphous solids [25,[28][29][30][31][32][33], physical mixtures [22,25,34], salts [35], solvates [35,36], hydrates [35][36][37] and solid solutions [38][39][40]. These solid forms can consist of a single component or multiple components. Single component solids include amorphous forms and polymorphs. Multicomponent solids can be cocrystals [41,42], salts [41], coamorphous solids [28], polymorphs of cocrystals [43], solvates/hydrates [36,37], and continuous solid solutions or discontinuous solid solutions/eutectics [34]. Figure 2 and Table 1 explains the differences between different solid forms of an API. As new forms of pharmaceutical solids are being produced with different type of guest molecules (solvents/water/solids), assigning an exact nomenclature for each of the solid was a hotly debated topic for many years [44][45][46]. Inspired by the work of Aitipamula et al. [45], Grothe et al. [46] developed a straightforward system for classification of multicomponent solid forms. The classification system proposed by Grothe et al. [46] is shown in Figure 3.  Grothe et al. [46] [Reprinted from [46] with permission. Copyright 2016 American Chemical Society].  [34,55] Melting point of the solid formed Lower melting than the parent components Mostly in between the melting points of parent molecules but may also be higher or lower than the parent molecules Exhibits a solidus-liquidus melting behavior [34,55] Binary phase diagram 'V'-shaped curve 'W'-shaped curve Unary phase diagram [34,55] Intermolecular interactions Short-range and weaker non-covalent adhesive interactions Stronger and non-covalent adhesive interactions (hydrogen bonding, halogen bonding, Π-Π interactions, etc.) Stronger and non-covalent cohesive interactions [34,55] Arrangement of molecules in crystal lattice Molecules are randomly arranged Molecules are well organized and well-packed Molecules are well organized [34,55] Predominant thermodynamic force of the system Entropy Enthalpy Enthalpy [34,55] Crystal structure of the solid formed No significant change from the parent components Characterized by a new crystal phase formation and hence possess a new crystal structure Characterized by a new crystal phase formation and therefore possess a new crystal structure [34,55] Thermodynamic stability of the solid formed Less stable More stable than eutectics More stable than eutectics [34]  The cocrystal phase showed 18 times increased intrinsic dissolution rate than the commercial API [58] AMG 517 A transient receptor potential vanilloid 1 antagonist (TRPV1)

Eutectics
Eutectics are another important class of multicomponent solids that has been gaining tremendous attention in pharmaceutical research in the recent years. When attempts to obtain a cocrystal fail, one may end up with a eutectic. Eutectics are multicomponent crystalline materials which exist in the form of discontinuous solid solutions (as shown in Figure 2). Cherukuvada and Nangia [34] defined eutectics as a conglomerate of solid solutions [34]. Unlike cocrystals, eutectic phases are stabilized by weaker adhesive interactions between the unlike molecules or stronger cohesive interactions between the molecules having similar structure [34,55]. Eutectic phases do not have a distinct or a specific crystal structure. Instead their crystalline nature resembles the combination of crystalline nature of the parent components. Eutectics possess a melting point less than the melting point of the pure components. Table 1 presents characteristic features of eutectics vis-à-vis' other solid forms.
Eutectics can be as good as cocrystals in fine-tuning the physicochemical properties of an API [55]. Duarte et al. [76] reported a eutectic dispersion of fenofibrate with a low molecular weight polymer, polyethylene glycol. This eutectic dispersion is commercially being available with the trade name 'Fenoglide' [76,77]. These eutectic dispersions possessed lower melting point and crystalline nature similar to the parent components [76]. Furthermore, Faeges [78] formulated liquefied form of a eutectic of aspirin with 2-3 parts of glycerin or propylene glycol (w/v) [78] as an ointment for topical applications. These eutectic formulations were also reported to enhance the shelf-life of aspirin by preventing its hydrolysis [78].
Several reports are available in the literature where eutectics were used for increasing the dissolution of poorly water-soluble drugs. The summary of reports available in the literature on various drug eutectics and their dissolution behavior are presented in Table 3. From Table 3, it is evident that the eutectic mixtures show significantly enhanced dissolution as compared to the raw drug. This ability of eutectic mixtures to exhibit enhanced dissolution rates than the raw drug is attributed to the randomized lattice arrangement which exists in the eutectic phase [24].
Cherukuvada and Row [55] employed hydrogen-bonding principles to understand and design generalized rules involved in cocrystal/eutectic/solid solution formation by investigating cocrystallization of 4,4 -bipyridine, isonicotinamide, isoniazid Fluoxetine hydrochloride drug systems as a model for their study [55]. Furthermore, Cherukuvada and Nangia [34] formulated ground rules based on the crystal engineering approach to explain the circumstances favoring cocrystal/eutectic/solid solution formation with respect to the structural similarity and intermolecular interactions existing in the binary system. While eutectic formation occurs with a pair of molecules having either similar or dissimilar structures [23,34] and sustained by weaker adhesive interactions, cocrystals are stabilized by strong adhesive interactions [34]. Table 1 presents the characteristic features of cocrystals vis-à-vis eutectics and solid solutions in comparison with cocrystals.

Solid Solutions
A solid solution is a homogeneous phase formed out of a solid-state reaction between the drug and coformer which are miscible into each other. When two molecules possess structural similarity (isomorphous and isostructural), then the resultant solid form is generally a solid solution [38][39][40]. Solid solutions are stabilized by strong cohesive interactions. Table 1 presents salient features of solid solutions and explains the difference in how different a solid solution is vis-à-vis other solid forms. Table 4 presents the summary of a few reports available in the literature on drug solid solutions.

Coamorphous Solids
Coamorphous solids are another interesting group of pharmaceutical solids. Coamorphous solids are solid forms in which the amorphous state is stabilized by weak intermolecular interactions [85,86] between the drug and coformer molecules. Coamorphous solids are used to enhance the bioavailability of hydrophobic drugs as their amorphous nature can improve aqueous solubility of an API. The term 'coamorphous' was first coined by Chieng et al. [85]. A drug and the polymer combination may result into 'Amorphous Solid Dispersions (ASDs)' whereas a mixture of amorphous drug and low molecular weight coformers can form a 'coamorphous solid' [87]. Weak intermolecular interactions (such as hydrogen bonding and/or Π-Π interactions [85,86]) dominate between drug and coformer in coamorphous solids. The coamorphous solids possess glass transition temperature (T g ) in between the glass transition temperatures of individual components [86]. There are also reports where coamorphous solids with no well-defined stoichiometry have been reported [30]. Table 5 presents a summary of few literature reports where coamorphous solids of APIs have been reported.
From the literature reports (as shown in Table 5), it is evident that the coamorphous solids showed enhanced dissolution rates than the raw drug molecules [25,[28][29][30][31]33,88]. This enhancement in dissolution rate mainly depends on two important factors namely: (i) capability of the drug to maintain higher supersaturation level and (ii) the strength of the API-coformer interaction to slow down the nucleation. In a study on Curcumin-artemisin coamorphous solid, Nangia and coworkers [28] have stated that coamorphization induces micronization of particles and in a way also leads to enhanced dissolution rates [28]. From the literature information available on coamorphous solids, one can conclude that coamorphous solids are one of the byproducts of a cocrystallization process which can enhance dissolution of hydrophobic drugs in aqueous medium [89], and in turn their bioavailability. Also, it was observed that the higher level of apparent solubility attained with the coamorphous solids during in vitro dissolution studies also resulted in enhanced in vivo bioavailability [89,90]. However, there is a need to investigate the long-term dissolution rates of pharmaceutical coamorphous solids. Curcumin-piperazine coamorphous phase showed lower dissolution than raw curcumin at temperature above T g of the coamorphous solid and exhibited higher dissolution than raw curcumin at a temperature below T g [29] Folic Acid The nanoconfined coamorphous solid showed enhanced dissolution than raw ibuprofen [33] Pharmaceutics 2018, 10, 108 13 of 74

Salts
Salt formation is one of the traditional methods employed to enhance the aqueous solubility of poorly water-soluble drugs. More than half of the medicines in market exist in the form of salt [90]. Salts are formed because of intermolecular hydrogen bonding due to proton transfer between the molecules with ionizable functional groups. Figure 4 presents the schematic representation of the difference between a salt and cocrystal [91]. Hence, salt formation is favored only when the API contains an ionizable site in it [15,92]. Formation of an API salt with a coformer occurs when there is a proton transfer from an acid to a base in the ionic state [91,93]. According to Sarma et al. [44], API molecules with nitrogenous functionality and -COOH functional groups in their chemical structure are more susceptible to salt formation [44] since these molecules can favor proton transfer to the coformer molecules. It has been also reported that presence of adequate number of counterions in an API and coformer molecule facilitates salt formation [94]. The presence of charge-assisted hydrogen bonding in salts enables coformer molecules to dissociate easily from the salt complexes resulting in peak solubility of drug in a dissolution medium within a few hours. Moreover, from thermodynamic point of view, salts possess higher enthalpy of hydration which also facilitates attainment of higher dissolution rate [95].

Salts
Salt formation is one of the traditional methods employed to enhance the aqueous solubility of poorly water-soluble drugs. More than half of the medicines in market exist in the form of salt [90]. Salts are formed because of intermolecular hydrogen bonding due to proton transfer between the molecules with ionizable functional groups. Figure 4 presents the schematic representation of the difference between a salt and cocrystal [91]. Hence, salt formation is favored only when the API contains an ionizable site in it [15,92]. Formation of an API salt with a coformer occurs when there is a proton transfer from an acid to a base in the ionic state [91,93]. According to Sarma et al. [44], API molecules with nitrogenous functionality and -COOH functional groups in their chemical structure are more susceptible to salt formation [44] since these molecules can favor proton transfer to the coformer molecules. It has been also reported that presence of adequate number of counterions in an API and coformer molecule facilitates salt formation [94]. The presence of charge-assisted hydrogen bonding in salts enables coformer molecules to dissociate easily from the salt complexes resulting in peak solubility of drug in a dissolution medium within a few hours. Moreover, from thermodynamic point of view, salts possess higher enthalpy of hydration which also facilitates attainment of higher dissolution rate [95].

Salt-Cocrystal Continuum
Salt-cocrystal continuum is other interesting subset of multicomponent pharmaceutical solids which falls under cocrystal/salt category. When a multicomponent solid form contains mixed ionization states (the extent of proton transfer from one molecule to the other is not predictable), it is difficult to understand whether the resultant solid form is a salt or cocrystal. This mainly occurs when the difference in pKa value of the drug and coformer lies between 0 and 3 [93] (Explained later in Section 3.1). Childs et al. [93] investigated the influence of crystal structure on the ionization states of the salt-cocrystal continuum. Jacobs and Noa [96] reported a hybrid salt-cocrystal methanol water solvate of p-Coumaric acid-Quinine with unexpected stoichiometry prepared by slow evaporation [96].

∆pKa Rule
∆pKa value has been used to assess cocrystal formation ability of a coformer with a given API [93,97]. pKa value (negative logarithm of dissociation constant) indicates the ability of an acid molecule to give up a proton [93]. When the difference between pKa value of API and coformer (∆pKa) ranges in negative values, there will be no proton transfer [93,98]. Therefore, one can possibly expect cocrystal formation in such a case [93,[97][98][99]. On the other hand, salt formation is observed when ∆pKa value is greater than 3 due to completion of proton transfer [93,98,99]. According to Berry and Steed [98], when ∆pKa value remains close to that of a base, then the system forms a salt and when it exists close to the acid, then the system forms a cocrystal [98]. da Silva et al. [71] designed and developed five 5-Fluorocytosine cocrystals with adipic, succinic, terephthalic,

Salt-Cocrystal Continuum
Salt-cocrystal continuum is other interesting subset of multicomponent pharmaceutical solids which falls under cocrystal/salt category. When a multicomponent solid form contains mixed ionization states (the extent of proton transfer from one molecule to the other is not predictable), it is difficult to understand whether the resultant solid form is a salt or cocrystal. This mainly occurs when the difference in pKa value of the drug and coformer lies between 0 and 3 [93] (Explained later in Section 3.1). Childs et al. [93] investigated the influence of crystal structure on the ionization states of the salt-cocrystal continuum. Jacobs and Noa [96] reported a hybrid salt-cocrystal methanol water solvate of p-Coumaric acid-Quinine with unexpected stoichiometry prepared by slow evaporation [96].

∆pKa Rule
∆pKa value has been used to assess cocrystal formation ability of a coformer with a given API [93,97]. pKa value (negative logarithm of dissociation constant) indicates the ability of an acid molecule to give up a proton [93]. When the difference between pKa value of API and coformer (∆pKa) ranges in negative values, there will be no proton transfer [93,98]. Therefore, one can possibly expect cocrystal formation in such a case [93,[97][98][99]. On the other hand, salt formation is observed when ∆pKa value is greater than 3 due to completion of proton transfer [93,98,99]. According to Berry and Steed [98], when ∆pKa value remains close to that of a base, then the system forms a salt and when it exists close to the acid, then the system forms a cocrystal [98]. da Silva et al. [71] designed and developed five 5-Fluorocytosine cocrystals with adipic, succinic, terephthalic, benzoic, and malic acid based on the ∆pKa values for API-coformer [71]. When the ∆pKa value ranges between 0 and 3 (partially ionized states), in such cases, these solid forms were referred to as Salt-cocrystal continuum [93,99,100]. Interestingly, Nangia and coworkers [100], while attempting to develop Clotrimazole (CLT) cocrystals with some carboxylic acid coformers, identified salt formation with maleic acid (MA) at a stoichiometric ratio of 1:0.5 (CLT:MA) whereas the calculated ∆pKa value for the system was 0.93 [101]. Thus, an ∆pKa value cannot always be used to predict/confirm the nature of a solid phase in all the cases and an experimental analysis is required for accuracy.

Hydrogen Bond Donors and Acceptors
The number of hydrogen bond donors and acceptors in a coformer and drug molecules also determines the extent of success in a cocrystallization event. Molecules that can form multiple hydrogen bonds are likely to form cocrystals with the coformer molecules [102]. Etter [103] and Donohue [104] framed HydrogenBond Rules to predict the circumstances under which hydrogen bond interactions that result into cocrystals [54,103,104]. These rules are as given below: The best proton donors and acceptors available after intramolecular hydrogenbond formation then participate in intermolecular hydrogen bonds d.
All acidic hydrogen atoms are included in hydrogen bonding in the crystal structure

Molecular Recognition Points
Almarsson and Zaworotko [105] pointed out that the API molecules contain certain functional group (or molecular recognition point) in their structure which interacts with the coformer and thereby create a supramolecular unit (or molecular recognition point) called supramolecular synthons [105]. The term 'Synthon' was first introduced by Corey in 1967 who defined Synthons as "Structural units within supermolecules which can be formed and/or assembled by known or conceivable synthetic operations involving intermolecular interactions" [106]. Desiraju [107] defined supramolecular synthons as spatial arrangement of intermolecular interactions which serves as a base for any supramolecular synthesis [107]. Thus, synthons are design elements in crystal engineering which are different from the term 'intermolecular interactions' [108]. Sometimes, a synthon can also be represented as a single interaction [107] (such as Cl. . . Cl and N. . . Br interactions) in a few supermolecular structures.

Flexibility of Synthon-Forming Functional Groups
In addition to molecular recognition points, the position of functional groups and the conformational flexibility of participating molecules play a significant role in determining success rate of a cocrystallization. For instance, Nangia and coworkers [22,72] identified that resorcinol can cocrystallize with curcumin whereas hydroquinone and catechol could not cocrystallize with curcumin though all the three molecules possessed same functional groups [22,72]. Such observations suggest that understanding the rationale behind formation of supramolecular synthons in the crystal lattice of the cocrystals is highly necessary for cocrystal design and development.
Aakeroy et al. [108] carried out an extensive study to understand how polymorphic compounds serve as good cocrystallizing agents/coformers and emphasized the significance of flexibility of synthon-forming functional groups of coformers during any cocrystallization [108]. They experimentally studied the cocrystal forming ability of three polymorphic compounds, isonicotinamide, 2-amino 3-nitropyridine, 4-chlorobenzamide and maleic hydrazide [108]. It was observed in their study that isonicotinamide, 2-amino 3-nitropyridine and 4-chlorobenzamide participated actively in intermolecular hydrogen bonding with a variety of aliphatic and aromatic carboxylic acids and thereby favored the formation of binary/ternary cocrystals whereas maleic hydrazide was not found suitable candidate for cocrystallization with aromatic or aliphatic compounds having acid, amide and oxime functional groups. This was attributed to the flexibility of functional groups in coformers in addition to their polymorphic nature. Isonicotinamide, 2-amino 3-nitropyridine and 4-chlorobenzamide exhibit hydrogen bonding between different functional groups in each of their polymorphs whereas all the three polymorphs of maleic hydrazide always exhibited the primary hydrogen-bonding interactions between same functional groups which decreases the possibility of formation of new hydrogen bond synthons with the other molecules.

Carbon Chain Length of Dicarboxylic Acid Coformers
Carboxylic acids are some of the most commonly used coformers for cocrystallization of many small molecules since they can form heterosynthons with molecules containing amide and pyridine functional groups and homosynthons with API molecules containing acid functional group. However, the cocrystal forming tendency of carboxylic acids also depends on the length of carbon chain in it. Shevchenko et al. [109] while investigating cocrystallization of itraconazole with different aliphatic dicarboxylic acids containing carbon chains of varying length observed that as the length of carbon chain in the coformer molecules increases, the packing of these molecules within the crystal lattice of drug molecules becomes increasingly difficult due to steric hindrance or incompatibility of large carbon chain to exactly fit into the crystal lattice of drug [109]. During their study, the researchers identified that itraconazole formed cocrystals with oxalic acid (C2), adipic acid (C6), malonic acid (C3), glutaric acid (C5) and pimelic acid (C7) and not with suberic acid (C8), azelaic acid (C9) and sebacic acid (C10). Based on this observation, it can be safely concluded that coformers containing longer carbon chains possibly are not suitable candidates for cocrystallization with API molecules where the probability of geometrical fitness of coformer into the API lattice is low [109].

Effect of Solvents
Solubility of the API and coformer in a solvent used for cocrystallization plays a significant role in determining the success of cocrystallization experiment. The solubility of the individual components must be determined a priori to cocrystallization experiments [110]. The Phase Solubility Diagram (PSD), also called as Ternary (API-Coformer-Solvent) Phase Diagram can be constructed using this data which then serves as a fundamental tool to identify region of cocrystal formation, understand the solution chemistry and solubility behavior of cocrystals [110]. The polarity of the solvent system, solubility of API and coformer, temperature and pH are the important parameters which determine the cocrystal forming zone in a ternary system. Robertson and his coworkers [111] observed that the polarity of the solvent determined the type of non-covalent interactions (hydrogen-bond or halogen-bond), and thereby controlling intermolecular interactions in the cocrystal phases [111]. It was summarized that the hydrogen-bonded cocrystals formation was favored by less polar solvents (such as toluene) whereas the more polar solvents (such as chloroform, dichloromethane, acetone, acetonitrile, nitromethane and 1-propanol) favored the formation of halogen-bonded cocrystals and in some cases, mixed halogen and hydrogen-bonded cocrystals [111]. This is mainly attributed to the influence of polarity of the different solvents on the strength of intermolecular interactions.

Screening Methods for Cocrystals
As cocrystallization is influenced widely by several important parameters, selecting a suitable coformer for cocrystallization requires an effective screening process. Screening of a suitable coformer can be carried out experimentally or computationally. Experimental methods are exhaustive and time-consuming. On the other hand, computational methods can serve as a rapidscreening tool for initial assessment of coformers that are suitable for cocrystallization process. Sections 4.1 and 4.2 presents the various computational and experimental methods that can be used for coformer screening.

Computational Methods
Most of the computational methods reported for coformer screening in literature till date are mainly thermodynamics-based methods. Issa et al. [112] used lattice energy calculations as an effective approach to screen coformers to form thermodynamically stable cocrystals [112]. If the lattice energy of the cocrystal is lower than the sum of lattice energies of individual components, then the cocrystal phase is said to be a thermodynamically stable phase [112]. Apart from calculation of lattice energies [112][113][114][115], several other parameters have also been employed for computational prediction of successful cocrystal formation or coformer screening. These include calculation of interaction energies [116,117], electrostatic potentials [118], molecular complementarity between API and coformers [119,120], solubility behavior [121][122][123], crystal energy landscapes of API-coformer pairs [124] and hydrogen bond propensities [125]. Table 6 presents the summary of various computational coformer screening methods  reported so far in the literature. Despite being less labor intensive, computational screening methods suffer from the requirement of large simulation times to perform molecular dynamic simulations.

Experimental Methods
Newman [133] has summarized various methods reported in the literature to obtain various solid forms such as polymorphs, salts, cocrystals and amorphous solid dispersions [133]. Though hydrogen-bonding motif structures play a significant role in stabilizing the crystal structures of cocrystals [134][135][136][137][138][139], in some cases, evaluation of crystal structures of some cocrystals synthesized using techniques such as Solution crystallization, mechanocrystallization and hot-stage microscopy showed that these cocrystals do not possess hydrogen-bonding motifs as given by Crystal Engineering principles [124]. Therefore, a suitable coformer screening technique is highly essential for a cocrystallization process. Given below is a brief account of different experimental techniques which could be used for screening coformers for cocrystallization.  Calculation of energy difference N,N -Diphenylureas and Triphenylphosphine Oxide Relative ∆∆E int between heterodimers and homodimers served as a good predictor of cocrystal formation in the investigated system [129] Free energy calculation Pentoxifylline with coformers, aspirin, salicylic acid, and benzoic acid FLEXCRYSTprogram was used to calculatethe free energy of experimental and hypothetical crystal structures and then correlated with each other [130] In silico screening Phenylpiperazine derivatives anddicarboxylic acids Computed values for the mixing enthalpies and solubility advantage factors determine the cocrystallization propensities [131] Ab initio screening Nicotinamide, isonicotinamide, picolinamide and two paracetamol cocrystals were screened for their stability Lattice energy calculations [114] Heat of formation distribution of components 492 pairs Calculation of Heat of formation (Hf) values [132] Fabian approach 974 cocrystal structures were investigated by Fabian approach Calculation ofQuantitative Structure-Activity Relationship (QSAR) molecular descriptors of cocrystallizing components [125] Utilization ofcomputed crystal energy landscapes Carbamazepine (CBZ)-Nicotinamide (NCT), Computed crystal energy landscapes of the drug-coformer pairs were compared and analyzed with the binary and ternary phase diagrams Molecular complementarity, hydrogen bond motifs and multicomponent hydrogenbond propensities were compared and analyzed [119] Knowledge-based hydrogen bond prediction Pyrimethamine drug with carbamazepine, theophylline, aspirin, α-ketoglutaric acid, saccharin, p-coumaric acid, succinimide and L-isoleucine as coformers Hydrogen bond propensity calculations were carried out to detect the formation of new molecular adducts [120] Multistage lattice energy minimization methodology

Differential Scanning Calorimetric (DSC) Analysis
Thermal analysis can be used for selection of a suitable coformer to form a cocrystal with a desired drug molecule. While one can detect the polymorphic transformation exhibited by the drug/coformer, a formation of a new phase (such as cocrystal), with the multicomponent reactant molecules (API and coformer) can also be easily detected. DSC has been reported to be a rapid thermal screening tool by Lu et al. [62]. Similarly, Yamashita et al. [140,141] conducted an extensive study using DSC as a coformer screening tool to identify formation of many pharmaceutical cocrystals [140,141]. Recently, Saganowska and Wesolowski [142] used DSC as a rapid screening tool and identified 15 different benzodiazepine cocrystals (at stoichiometric ratio of 1:1) [142].
DSC can be used as a tool to identify the coformers capable of forming a cocrystal with an API molecule using either of the following two approaches:

(a) Nature of DSC Thermograms
The nature of a DSC thermogram obtained for a binary mixture of API and coformer at a specific stoichiometric ratio can help in identifying the nature of the solid form as explained below: Physical mixture-If heating of a binary mixture at a specific heating rate yields two endotherms, each corresponding to the melting point of the individual components, in a DSC thermogram, then it suggests that the binary mixture remains as a physical mixture with no intermolecular interactions occurring between these molecules [25,34].
Eutectic-A eutectic is said to have formed when the DSC thermogram shows a single endotherm with melting temperature less than the melting points of either of the parent molecules (the drug and the coformer) [34,55,62,140,141].
Cocrystal-A binary system can be said to be capable of forming a cocrystal when the DSC thermogram exhibits a single endotherm with melting points lower or in between or greater than the melting points of the individual components [15,34,140,141]. DSC analysis along with single crystal XRD analysis can confirm the formation of a new cocrystal phase.
Furthermore, if a DSC thermogram shows consecutive multiple peaks, it possibly indicates the cocrystal formation. In case of presence of two consecutive peaks, the first melting peak corresponds to the eutectic melting, followed by the cocrystal melting (represented by the second endotherm). It is also possible that two or multiple peaks are an indication of polymorphic transformation of cocrystals [62,140].
In case of presence of three consecutive endotherms in a DSC thermogram, the first represents the eutectic melt, the second represents the melting of excess component (drug or coformer), followed by the final melting of cocrystal in the third endotherm [24,62,140].
Coamorphous solid-The DSC thermogram of coamorphous solids show a glasstransition phase change implying the formation of a coamorphous solid phase [25,28,143].
Solid solution-If a DSC thermogram exhibits a simple solidus-liquidus behavior, it indicates the formation of a solid solution [34]. When two components in a binary system form a complete solution in both the solid and liquid phases, then this behavior is termed a simple solidus-liquidus behavior. If the two components in the binary system possess similar crystal structure, then the intermolecular interactions in these compounds remains similar and the difference in the size of atoms remains small. Such isomorphous compounds can result into a solid solution.

(b) Binary Phase Diagrams
DSC thermograms obtained for different mixtures containing drug and coformer molecules in varying stoichiometric ratios can be used to construct a binary phase diagram for any drug-coformer pair. These binary phase diagrams can be used to detect the cocrystal formation zone for a given system. Section 4.2.2 explains this in detail.

Phase Diagrams
Phase diagrams are utilized as means to identifying different solid phases that can be formed between any drug-coformer pair. These phase diagrams can be generated either for two components (API-coformer) or for three components (API-coformer-solvent) as well.

(a) Binary Phase Diagrams
Binary phase diagrams are generally constructed with the data points obtained from thermal analysis methods such as DSC analysis [140]. The onset temperature of the first endotherm in a DSC thermogram is generally chosen as solidus point. The peak temperature of the second endotherm is chosen as liquidus point for constructing the phase diagram. The melting behavior of an API and coformer (congruent/incongruent melting) determines the solid solution/eutectic and cocrystal forming property for the investigated system. In general, the eutectic forming binary system adopts 'V'-shaped curve whereas the cocrystal forming system adopts 'W'-shaped curve indicating cocrystal formation between two eutectics. Reports are available in the literature where researchers have used binary phase diagram for determining the cocrystal formation zone for the investigated drug-coformer combination [24,25,34,50,73,140,141,144,145]. Sangster [80] has reported a set of binary phase diagrams for 60 different drug systems constructed with the help of computer-coupled phase diagram/thermodynamic analysis [80].
The commonly observed binary phase diagram for most of the drug-coformer pairs reported so far in the literature are shown in Figure 8.

(b) Ternary Phase Diagrams
Solution crystallization experiments can result in cocrystals. However, sometimes multicomponent single crystals are not formed during solution crystallization due to the incongruent solubility behavior of individual components in the solvent system. Therefore, the thermodynamic behavior of a system involving ternary components needs to be ascertained before attempting solution crystallization [146]. Ternary Phase Diagrams (Solute-solute-solvent phase diagram) help in determining the cocrystal formation region for a given system. Many researchers have employed ternary phase diagram as an optimization tool to determine suitable stoichiometric ratio of API and coformer for cocrystal formation [146][147][148][149][150]. When the two components manifest a similar/congruent solubility, the phase diagram exhibits a more symmetrical trend [148]. On the other hand, when the two components show dissimilar/incongruent solubility, then the phase diagram will be less symmetrical behavior [148]. Figure 9A,B presents schematics of a symmetric and an asymmetric ternary phase diagram.

Synthesis of Cocrystals
Cocrystals can be synthesized using several methods which could either be the batch processes or continuous processes. Given below is a brief account of different techniques which are being used to synthesize cocrystals.

Synthesis of Cocrystals
Cocrystals can be synthesized using several methods which could either be the batch processes or continuous processes. Given below is a brief account of different techniques which are being used to synthesize cocrystals.

Batch Processes
The batch mode cocrystallization techniques include mechanochemical methods such as Solid-State grinding [151,152], Liquid-Assisted grinding [151], Ion and Liquid-Assisted Grinding [153] and solvent-based techniques such as slow evaporation [154], ultrasound-assisted cocrystallization [155,156], slurry conversion [140,157] or solvent-mediated phase transformation [158], generation of cocrystals from moisture [159,160] and anti-solvent precipitation (Liquid Anti-Solvent precipitation) [161] and Gas Anti-Solvent precipitation [162]. Table 7 presents the principle involved behind these batch cocrystallization techniques [151,[153][154][155]159,[161][162][163]. Mechanochemical methods, ultrasound-assisted cocrystallization, cocrystal generation using moisture are discussed in the Section 5.1 in detail. Liquid-Assisted Grinding (LAG) Addition of solvent molecules during grinding of API-coformer mixture facilitates diffusion of low molecular weight coformer molecules into the crystal lattice of an API [151] Ion and Liquid-Assisted Grinding (ILAG) Grinding of multiple components along with liquid phase facilitates diffusion of liquid phase into the solid-state thereby facilitating mobility of molecules and exposing the hidden molecules to the surface [153] Slow evaporation Preparing a supersaturated solution and allowing solvent to evaporate slowly at ambient conditions induces primary nucleation, thereby leading to slow growth of crystals [154] Ultrasound-assisted cocrystallization The cavitation energy of ultrasonic waves induces primary nucleation of particles and leads to attainment of supersaturation for cocrystal growth [155] Solvent-Mediated Phase Transformation (SMPT) Keeping the binary mixture of API and coformer into a solvent/solvent mixture for a long-time enables phase transformation to a new phase (cocrystal) or conversion from one polymorphic form to another due to the activity of solvent molecules and their nature of interaction with solute molecules over a period [163] Cocrystal formation from moisture Uptake of moisture by a solid mixture, dissolution of reactant molecules in the mixture, attaining supersaturation leading to cocrystal nucleation and growth are the primary steps involved in cocrystal generation from moisture [159] Liquid Anti-Solvent precipitation (LAS) Supersaturation of molecules (drug + coformer mixture) is attained by mixing of solvent and anti-solvent. This results in increase in supersaturation and nucleation rate and thereby facilitates production of sub-micron cocrystal particles [161] Gas Anti-Solvent precipitation (GAS) In GAS process, gaseous phase is used as anti-solvent instead of a liquid to obtain a cocrystal from a solution of API and coformer from a liquid solution [162] (

a) Mechanochemical Methods (Grinding)
The API and coformer molecules can be cocrystallized by grinding together without solvent, which is termed as Solid-State Grinding (SSG) or Neat grinding or with a few drops of solvent [termed as Solvent-Drop grinding or Liquid-Assisted Grinding (LAG)]. During grinding, the low molecular weight components (solid/liquid/gas phase) diffuse easily into the API crystal lattice [152] which can result in a formation of intermediate phases such as eutectic or amorphous phase which can further result in formation of a cocrystal. According to Etter [103], conversion of individual parent molecules into hydrogen-bonded cocrystals basically occurs based on the rate and force of grinding, particle size of the components and vapor pressure of the phases [103].
Chadwick et al. [152] explained different mechanisms of cocrystal formation via grinding and the various thermodynamic parameters that influence cocrystal formation during grinding. Friscic [153] explored ion and liquid-assisted grinding (ILAG) in designing multicomponent cocrystals. It was proposed that grinding of multiple components along with liquid facilitates diffusion of liquid phase into solid-state thereby enabling mobility of molecules and exposing the hidden molecules to the surface.
Numerous reports are available in the literature where grinding (solid-state or liquid-assisted) methods were used extensively to produce cocrystals [164][165][166][167][168][169][170][171][172]. Table 8 presents a summary of a few reports available in the literature where the researchers have proposed scientific mechanisms behind solid-state or liquid-assisted grinding process [151,152,167,168,171,173,174].

(b) Ultrasound-assisted Cocrystallization
Ultrasound has been widely used for inducing nucleation in solution and cocrystallizing small molecules. Ultrasound-assisted cocrystallization has several advantages over conventional solution crystallization. The mechanical energy released during passage of ultrasonic waves induces primary nucleation at lower supersaturation levels, thereby reducing the induction time and metastable zone width [175,176]. Ultrasound can therefore induce crystallization easily from solution which otherwise is difficult to attain by conventional solution crystallization experiments [155]. Reports are available in the literature where ultrasound-assisted solution cocrystallization (USSC) [60,155,156,177] has been reported as a method to crystallize caffeine-maleic acid (2:1) [155] and caffeine-oxalic acid (2:1) [155,156].
Morrison et al. [178] proposed a bench-top inexpensive method called 'Solvent-Drop Floating Foam Rack/Sonic Bath Method' for synthesizing cocrystals and tested the feasibility of this method for synthesis of carbamazepine-saccharin (1:1) form II, carbamazepine-nicotinamide (1:1), AMG 517-benzoic acid (1:1) and AMG 517-malic acid (1:1) cocrystals. In this technique, the moist powders were subjected to sonication in a sonic bath for an optimized period (90 min was used for carbamazepine and AMG 517 drug systems since sonication for less than 5 min yielded only partial conversion of cocrystal for many systems) using a floating foam rack for inducing cocrystallization [178]. The time duration provided for sonication and a suitable solvent selection were found to be the two critical parameters which influence cocrystal formation by this technique [178].

(c) Cocrystal Generation Using Moisture
Rodriguez-Hornedo and coworkers [159] conducted a study with API molecules namely carbamazepine, theophylline and caffeine to understand the mechanism behind cocrystal formation in presence of deliquescent additive such as sucrose or fructose [159]. Figure 10 presents a schematic of cocrystal formation from moisture in presence of deliquescent additive. The study concluded that uptake of moisture by the solid reactants facilitated nucleation and growth of carbamazepine-nicotinamide, carbamazepine-saccharin, caffeine-maleic acid, caffeine-oxalic acid, caffeine-malonic acid, caffeine-glutaric acid, theophylline-maleic acid, theophylline-oxalic acid, theophylline-malonic acid and theophylline-glutaric acid cocrystals.
Furthermore, Good et al. [160] attempted to understand the mechanism behind carbamazepine (CBZ)-nicotinamide (NCT) cocrystal formation facilitated by moisture sorption by Poly Vinyl Pyrollidone K12 (PVP K12). As PVP K12 decreased the ratio of cocrystal to drug solubility, increase in PVP K12 concentration increased the stability of carbamazepine (CBZ)-nicotinamide (NCT) cocrystal [160]. Figure 11 represents the optical microscopic images showing carbamazepine-nicotinamide cocrystal formation by moisture sorption, dissolution, cocrystal nucleation and growth in presence of PVP K12 as a deliquescent polymer at 75% RH and 25 • C.

(c) Cocrystal Generation Using Moisture
Rodriguez-Hornedo and coworkers [159] conducted a study with API molecules namely carbamazepine, theophylline and caffeine to understand the mechanism behind cocrystal formation in presence of deliquescent additive such as sucrose or fructose [159]. Figure 10 presents a schematic of cocrystal formation from moisture in presence of deliquescent additive. The study concluded that uptake of moisture by the solid reactants facilitated nucleation and growth of carbamazepine-nicotinamide, carbamazepine-saccharin, caffeine-maleic acid, caffeine-oxalic acid, caffeine-malonic acid, caffeine-glutaric acid, theophylline-maleic acid, theophylline-oxalic acid, theophylline-malonic acid and theophylline-glutaric acid cocrystals.
Furthermore, Good et al. [160] attempted to understand the mechanism behind carbamazepine (CBZ)-nicotinamide (NCT) cocrystal formation facilitated by moisture sorption by Poly Vinyl Pyrollidone K12 (PVP K12). As PVP K12 decreased the ratio of cocrystal to drug solubility, increase in PVP K12 concentration increased the stability of carbamazepine (CBZ)-nicotinamide (NCT) cocrystal [160]. Figure 11 represents the optical microscopic images showing carbamazepine-nicotinamide cocrystal formation by moisture sorption, dissolution, cocrystal nucleation and growth in presence of PVP K12 as a deliquescent polymer at 75% RH and 25 °C.

(c) Cocrystal Generation Using Moisture
Rodriguez-Hornedo and coworkers [159] conducted a study with API molecules namely carbamazepine, theophylline and caffeine to understand the mechanism behind cocrystal formation in presence of deliquescent additive such as sucrose or fructose [159]. Figure 10 presents a schematic of cocrystal formation from moisture in presence of deliquescent additive. The study concluded that uptake of moisture by the solid reactants facilitated nucleation and growth of carbamazepine-nicotinamide, carbamazepine-saccharin, caffeine-maleic acid, caffeine-oxalic acid, caffeine-malonic acid, caffeine-glutaric acid, theophylline-maleic acid, theophylline-oxalic acid, theophylline-malonic acid and theophylline-glutaric acid cocrystals.
Furthermore, Good et al. [160] attempted to understand the mechanism behind carbamazepine (CBZ)-nicotinamide (NCT) cocrystal formation facilitated by moisture sorption by Poly Vinyl Pyrollidone K12 (PVP K12). As PVP K12 decreased the ratio of cocrystal to drug solubility, increase in PVP K12 concentration increased the stability of carbamazepine (CBZ)-nicotinamide (NCT) cocrystal [160]. Figure 11 represents the optical microscopic images showing carbamazepine-nicotinamide cocrystal formation by moisture sorption, dissolution, cocrystal nucleation and growth in presence of PVP K12 as a deliquescent polymer at 75% RH and 25 °C.

Continuous Production of Cocrystals
Continuous production of cocrystals has gained much attention from pharmaceutical sectors [179][180][181][182]. In addition to enabling large-scale production of cocrystals, continuous production methods also offer significant advantages over batch processes [182] as mentioned below: a.
Batchwise variation in cocrystals quality can be rectified by means of continuous production by adopting a Quality-based Design (QbD) approach b.
Continuous production of cocrystals is less labor intensive when compared to batch mode c.
As compared to batch process, the maintenance of uniform particle distribution throughout the entire production, monitoring of quality of cocrystals by Process Analytical Technology (PAT) and maintaining consistency in the quality of products obtained is easier for continuous production Examples of different continuous cocrystal production methods are Twin Screw Extrusion [183], Hot-Melt Extrusion [65,184], Co-Extrusion, Solvent-Free Continuous Cocrystallization [185], Spray Drying [186], Spray Flash evaporation process [187], Solid-State Shear Milling Technology [188], melt crystallization [48] and melt-assisted grinding [179]. Table 9 presents a summary of a few literature reports available on continuous cocrystallization production techniques. Solid-State Shear Milling Technology [188], melt crystallization [48] and supercritical fluid-based methods [189,190] are discussed in detail in the following Section 5.2 (a-c).
Though continuous cocrystallization production methods can be employed to give high yield of cocrystals, these methods have certain limitations in terms of reproducible cycling efficiency and accurate process control. Maintaining stability and quality of products when transferred from one unit operation to another (Quality assurance) to make these products suitable for further use can be a challenge.

(a) Solid-State Shear Milling (S3M) Technology
Recently, Korde et al. [188] employed 'Solid-State Shear Milling' technology to produce carbamazepine-salicylic acid cocrystals with 5-25 wt. % of poly(ethylene oxide) (PEO). This 'Solid-State Shear Milling' technology was reported to be a scalable and polymer-assisted technology for continuous production of cocrystals [188]. Figure 12 presents the schematic of 'Solid-State Shear Milling (S3M)' technology employed by Korde et al. [188]. In this technique, well-mixed API and coformer are blended along with a polymer and then subjected to fine grinding by S3M. The S3M consists of two inlaid pans, a rotor and a stator which are made up of wear-resistant materials. The applied shear force for milling along with the polymeric aid facilitates efficient milling due to generation of high stress fields for grinding [188].

Continuous Production of Cocrystals
Continuous production of cocrystals has gained much attention from pharmaceutical sectors [179][180][181][182]. In addition to enabling large-scale production of cocrystals, continuous production methods also offer significant advantages over batch processes [182] as mentioned below: a. Batchwise variation in cocrystals quality can be rectified by means of continuous production by adopting a Quality-based Design (QbD) approach b. Continuous production of cocrystals is less labor intensive when compared to batch mode c. As compared to batch process, the maintenance of uniform particle distribution throughout the entire production, monitoring of quality of cocrystals by Process Analytical Technology (PAT) and maintaining consistency in the quality of products obtained is easier for continuous production Examples of different continuous cocrystal production methods are Twin Screw Extrusion [183], Hot-Melt Extrusion [65,184], Co-Extrusion, Solvent-Free Continuous Cocrystallization [185], Spray Drying [186], Spray Flash evaporation process [187], Solid-State Shear Milling Technology [188], melt crystallization [48] and melt-assisted grinding [179]. Table 9 presents a summary of a few literature reports available on continuous cocrystallization production techniques. Solid-State Shear Milling Technology [188], melt crystallization [48] and supercritical fluid-based methods [189,190] are discussed in detail in the following Section 5.2 (a-c).
Though continuous cocrystallization production methods can be employed to give high yield of cocrystals, these methods have certain limitations in terms of reproducible cycling efficiency and accurate process control. Maintaining stability and quality of products when transferred from one unit operation to another (Quality assurance) to make these products suitable for further use can be a challenge.

(a) Solid-State Shear Milling (S3M) Technology
Recently, Korde et al. [188] employed 'Solid-State Shear Milling' technology to produce carbamazepine-salicylic acid cocrystals with 5-25 wt. % of poly(ethylene oxide) (PEO). This 'Solid-State Shear Milling' technology was reported to be a scalable and polymer-assisted technology for continuous production of cocrystals [188]. Figure 12 presents the schematic of 'Solid-State Shear Milling (S3M)' technology employed by Korde et al. [188]. In this technique, well-mixed API and coformer are blended along with a polymer and then subjected to fine grinding by S3M. The S3M consists of two inlaid pans, a rotor and a stator which are made up of wear-resistant materials. The applied shear force for milling along with the polymeric aid facilitates efficient milling due to generation of high stress fields for grinding [188].

(b) Melt Crystallization
Melt crystallization has been adopted to produce cocrystals by several researchers [144,191,192]. This is an indirect way in which the cocrystal phase is formed from the melting of eutectic phase. Melting of API and coformer beyond the eutectic temperature creates a eutectic melt and the cocrystal phase growth appears from nucleation in the eutectic melt. For the drug-coformer pairs which form cocrystals from the eutectic melt, melt crystallization can yield cocrystal without any phase impurity [191,192]. Melt crystallization is a solvent-free, scalable and continuous method for cocrystal production [48].
Douroumis and coworkers [179] have contributed an extensive review on various solvent-free continuous cocrystal manufacturing methods [179]. Recently, Crawford [193] has contributed a review on the applications of Twin Screw Extrusion in continuous production of organic compounds [193]. Also, Shastri and coworkers [181] have recently provided a detailed review on the various challenges associated with continuous cocrystal production methods and the prospects in this technology [181].

(c) Supercritical CO 2 -Based Methods
Utilization of supercritical CO 2 as solvent or an anti-solvent instead of using liquid solvents serves as an excellent mean for large-scale production of cocrystals [194]. Pando et al. [162] had presented an extensive review on preparation of pharmaceutical cocrystals using supercritical CO 2 [162]. Efforts have been made by many researchers to prepare cocrystals by Gas Anti-solvent Crystallization (GAS) [189,190,[194][195][196][197][198] and to compare its morphology, crystalline nature and dissolution rates with cocrystals prepared via Liquid Anti-Solvent (LAS) precipitation or traditional solution-based methods [194][195][196]. In addition to GAS, Supercritical fluid Enhanced Atomization (SEA) [190], Atomization and Anti-solvent Crystallization (AAS) [191] and Supercritical Anti-solvent crystallization (SAS) [189] are a few other techniques where supercritical CO 2 has been utilized as an anti-solvent for cocrystals production. Rapid Expansion of Supercritical Solutions (RESS) is an interesting technique where both micronization of particles and cocrystallization can be achieved simultaneously using supercritical CO 2 as a solvent [199]. RESS involves no toxic solvents and therefore serves as an ecofriendly technique for a large-scale production of cocrystals. While supercritical CO 2 methods can be used for continuous synthesis of cocrystals, requirement of high pressure and especially designed nozzles for atomization can limit the usefulness of such techniques for large-scale production.

Structural Analysis
The crystalline nature of the cocrystals can be characterized by Powder X-Ray Diffraction (PXRD) and single crystal X-Ray Diffraction (SCXRD). While SCXRD provides detailed structural information including the lattice parameters, space group, miller indices, crystal system, unit cell volume, inter and intramolecular interactions, the PXRD provides information about the crystallinity of the solid phase. However, when the crystals produced by solution crystallization are not of good quality to conduct SCXRD analysis, one can possibly extract the structural data from PXRD data using indexing programs such as TREOR90 [200], ITO [201] and AUTOX [202], DASH [203], Rex.Cell [204] and Rietveld refinement programs such as TOPAS [205] and EXPO [206]. International Union of Crystallography (IUCr) serves as a valuable resource, containing information about different crystallographic software used for determining single crystal data of cocrystals from their corresponding PXRD data [207]. Several reports are available in the literature where researchers have solved the crystal structure of cocrystals from Powder X-Ray Diffraction data [208,209].

Thermal Analysis
Hot-Stage Microscopy (HSM or Kofler method): A thermal microscopic method provides a visual understanding of phase transitions which occur while heating an API-coformer mixture to a certain temperature [210][211][212][213][214][215]. It is also called as Mixed-fusion method in which one of the components (first component) melts initially. It is then followed by the solubilization of a second component into the molten component. This complete recrystallization results in the formation of 'Zone of mixing' which represents cocrystal formation from eutectic melt or polymorphic transformation of cocrystals (as shown in Figure 13A). HSM was reported to be used as a cocrystal/cocrystal polymorphs screening technique in several reports [216,217].

Thermal Analysis
Hot-Stage Microscopy (HSM or Kofler method): A thermal microscopic method provides a visual understanding of phase transitions which occur while heating an API-coformer mixture to a certain temperature [210][211][212][213][214][215]. It is also called as Mixed-fusion method in which one of the components (first component) melts initially. It is then followed by the solubilization of a second component into the molten component. This complete recrystallization results in the formation of 'Zone of mixing' which represents cocrystal formation from eutectic melt or polymorphic transformation of cocrystals (as shown in Figure 13A). HSM was reported to be used as a cocrystal/cocrystal polymorphs screening technique in several reports [216,217].  Figure 13B and C show an example of a cocrystal formation between lamotrigine and caffeine by using HSM [211]. Figure 13B indicates formation of cocrystal phase at 130 °C [202] followed by melting of eutectic phase at 190 °C whereas the cocrystal phase melt is observed when the temperature is increased to 200 °C (as shown in Figure 13C) [211]. The results obtained from HSM were also compared with the DSC thermograms obtained for lamotrigine-caffeine (2:1) pair. The results obtained in HSM and DSC analysis showed a good correlation between these two techniques. Thus, HSM/DSC screening method was suggested to be an efficient method for cocrystal screening when the quantity of sample available for cocrystal screening is less [211].
In addition to Hot-Stage Microscopy, Differential Scanning Calorimetry also serves as an effective thermal analysis technique for characterizing the cocrystal phases (as explained in Section 4.2).

Spectroscopic Analysis
There are several spectroscopic techniques which are commonly used to characterize the cocrystal phases. Fourier Transform (FT-IR) Infrared Spectroscopy (a vibrational spectroscopy) is a non-destructive analysis technique used for identifying hydrogen bonds and molecular conformations in a cocrystal phase with respect to its pure components based on the differences in the vibrational modes exhibited by the samples due to absorption of light.
Fourier Transform Near-Infra Red (FTNIR) Spectroscopy is a non-destructive vibrational spectroscopic technique which records the vibrational changes in the near-infra red region due to absorption of light by vibrating molecules. It is used to detect hydrogen bonds in a cocrystal phase based on the spectral changes seen in a spectrum of a cocrystal phase in comparison to near-infrared spectrum obtained for the parent molecules [218].
Raman spectroscopy is also a non-destructive spectroscopic technique used to detect the vibrational changes in the spectrum of a cocrystal phase which occurs due to scattering of incident  Figure 13B,C shows an example of a cocrystal formation between lamotrigine and caffeine by using HSM [211]. Figure 13B indicates formation of cocrystal phase at 130 • C [202] followed by melting of eutectic phase at 190 • C whereas the cocrystal phase melt is observed when the temperature is increased to 200 • C (as shown in Figure 13C) [211]. The results obtained from HSM were also compared with the DSC thermograms obtained for lamotrigine-caffeine (2:1) pair. The results obtained in HSM and DSC analysis showed a good correlation between these two techniques. Thus, HSM/DSC screening method was suggested to be an efficient method for cocrystal screening when the quantity of sample available for cocrystal screening is less [211].
In addition to Hot-Stage Microscopy, Differential Scanning Calorimetry also serves as an effective thermal analysis technique for characterizing the cocrystal phases (as explained in Section 4.2).

Spectroscopic Analysis
There are several spectroscopic techniques which are commonly used to characterize the cocrystal phases. Fourier Transform (FT-IR) Infrared Spectroscopy (a vibrational spectroscopy) is a non-destructive analysis technique used for identifying hydrogen bonds and molecular conformations in a cocrystal phase with respect to its pure components based on the differences in the vibrational modes exhibited by the samples due to absorption of light.
Fourier Transform Near-Infra Red (FTNIR) Spectroscopy is a non-destructive vibrational spectroscopic technique which records the vibrational changes in the near-infra red region due to absorption of light by vibrating molecules. It is used to detect hydrogen bonds in a cocrystal phase based on the spectral changes seen in a spectrum of a cocrystal phase in comparison to near-infrared spectrum obtained for the parent molecules [218].
Raman spectroscopy is also a non-destructive spectroscopic technique used to detect the vibrational changes in the spectrum of a cocrystal phase which occurs due to scattering of incident light. Furthermore, through solid-state Nuclear Magnetic Resonance (ssNMR) spectroscopy, one can obtain detailed structural information of pharmaceutical solid phases such as cocrystals, eutectics and amorphous solids [219]. ssNMR spectroscopy provides information regarding molecular mobility, differences in hydrogen bonding and crystallinity of solids. The cocrystal powders which are difficult to be characterized by single crystal X-Ray diffraction analysis (due to poor crystal quality) to obtain a crystal structure can be characterized by solid-state Nuclear Magnetic Resonance (ssNMR) spectroscopy [219,220].
Proton NMR spectroscopy is a nuclear magnetic resonance spectroscopy which can be used for determining stoichiometric ratio of cocrystal phases [221].
Furthermore, Variable Temperature X-ray Diffraction (VTXRD) analysis is used to determine phase transformation events which occur at a specific temperature during heating of a cocrystal phase. Through VTXRD, in addition to identification of polymorphic phase transformations, one can also determine the crystalline nature of the solid phase at different temperatures. Reports are available in the literature were VTXRD has been used to determine phase transformation events observed for different [222,223]. Vangala et al. [222] reported that dimorphs of caffeine-glutaric acid (1:1) cocrystal are enantiotropically related with the transition temperature of 79 • C [222] using VTXRD studies.
In addition to the characterization techniques discussed above, some of the well-advanced techniques such as pair-distribution function PXRD [224], Terahertz (THz) spectroscopic imaging [225], Probe-Type Low-Frequency Raman Spectroscopy [226], combined near IR and Raman spectroscopies [227], in situ Raman spectroscopy, in situthermomicroscopy [192], calorimetric analysis [192] and PAT [75] are being used widely by many researchers nowadays in order to understand the phase transition events that possibly occur during a cocrystallization process.

Morphological Characterization
Polarized optical microscopy [228], Fluorescence microscopy [229] and Scanning Electron Microscopy (SEM) [24,147,217] are the various techniques which can be used to characterize the morphology of cocrystals.

Drug Release Testing
According to United States Food and Drug Administration (US FDA), to ensure that pure drug is released from the cocrystal in vivo, evaluation of the performance of cocrystals by in vitro dissolution study or solubility testing is highly mandatory [230,231]. Determination of cocrystal solubility, determination of powder and intrinsic dissolution rates using United States Pharmacopeia apparatus [61,65,68,72,101] and evaluating cocrystal diffusivity by dialysis and Franz-diffusion cell method [232] are some of the in vitro drug release testing methods used for cocrystal formulations. Achieving in vitro-in vivo correlation (IVIVC) of a pure drug from its cocrystal form is an important criterion that must be satisfied by cocrystals for their use in pharmaceutical applications.

Influence of Intermolecular Interactions on Cocrystallization
The type and strength of intermolecular interactions that occur between the drug and coformer molecules determine the nature and the stability of the multicomponent solids (cocrystals/salts/coamorphous solids/eutectics/solid solutions/physical mixture) formed during cocrystallization. There are several types of intermolecular interactions (hydrogen bonding, secondary bonding, pi-pi interactions, ionic interactions, halogen bonding, vanderwaal's forces and dipolar interactions) which can possibly exist in multicomponent systems. However, most of the definitions for cocrystals are primarily based on considering hydrogen bonding as a predominant intermolecular interaction stabilizing the cocrystal lattice [233][234][235]. In 1940, Glasstone called cocrystals as 'lattice compounds' in his book, Textbook of Physical Chemistry and defined them as 'substances formed between stoichiometric amounts of "two" molecular species, which owe their stability to packing in the crystal lattice, and not to ordinary valency forces' [235]. Later, Dunitz [236] proposed that the word 'two' in Glasstone's definition should be changed to 'two or more' [236]. However, some of the definitions given for cocrystals by many researchers were solely based on the nature of components (gas phase/ionic molecules/solid components) which constitute the crystal lattice. Later, an unambiguous definition which covers all the characteristics of cocrystal had become necessary. Attempts were made by Lara-Ochoa and Espinosa-Perez [237] to propose a universal definition for cocrystals by considering three non-ordinary valence forces (π-π interactions and vanderwaal's forces, π-π stacking interactions and halogen bonding) stabilizing the cocrystal structure [237].
The following are the different types of intermolecular interactions reported to exist in the crystal structure of different cocrystals:

Hydrogen Bonding
Hydrogen bonding is a non-covalent interaction that occurs between a hydrogen atom attached to any electronegative atom such as oxygen or nitrogen and another electronegative atom. It is commonly represented as X-H-Y where X and Y are the electronegative atoms. Hydrogen-bonding interaction can be an intermolecular or intramolecular interaction. Hydrogen bonding is a specific type of dipolar interaction that serves as a key parameter in crystal engineering approach. Most of the cocrystal structures reported in Cambridge Structural Database (CSD) are stabilized by hydrogen-bonding interactions. Hydrogen bonds are individually weak bonds but collectively they are stronger enough to form supramolecular synthons in multicomponent systems and stabilize their structures [99,238,239]. Figure 14 presents the different types of hydrogen bond synthons which were reported to be observed in theophylline-cinnamic acid (1:1) cocrystals [240]. The crystal structure obtained for theophylline-cinnamic acid (1:1) cocrystal ( Figure 14) showed that cinnamic acid forms a carboxylic acid two-point synthon with carbonyl O of the pyrimidine ring and with the N-H hydrogen of the imidazole ring of theophylline to stabilize the crystal structure. This kind of stronger interaction prevents the formation of cocrystal hydrate [240].
The following are the different types of intermolecular interactions reported to exist in t ystal structure of different cocrystals:

Hydrogen Bonding
Hydrogen bonding is a non-covalent interaction that occurs between a hydrogen atom attache any electronegative atom such as oxygen or nitrogen and another electronegative atom. It mmonly represented as X-H---Y where X and Y are the electronegative atoms. Hydrogen-bondin teraction can be an intermolecular or intramolecular interaction. Hydrogen bonding is a speci pe of dipolar interaction that serves as a key parameter in crystal engineering approach. Most e cocrystal structures reported in Cambridge Structural Database (CSD) are stabilized b drogen-bonding interactions. Hydrogen bonds are individually weak bonds but collectively th e stronger enough to form supramolecular synthons in multicomponent systems and stabili eir structures [99,238,239]. Figure 14 presents the different types of hydrogen bond synthons which were reported to served in theophylline-cinnamic acid (1:1) cocrystals [240]. The crystal structure obtained f eophylline-cinnamic acid (1:1) cocrystal ( Figure 14) showed that cinnamic acid forms a carboxy id two-point synthon with carbonyl O of the pyrimidine ring and with the N-H hydrogen of th idazole ring of theophylline to stabilize the crystal structure. This kind of stronger interactio events the formation of cocrystal hydrate [240].

Halogen Bonding
In recent years, halogen-bonded cocrystals have gained an extensive attention. Figure   Figure 14.

Halogen Bonding
In recent years, halogen-bonded cocrystals have gained an extensive attention. Figure 15 presents the schematic of halogen bonding interaction between two molecules. Halogen bonds are more directional than hydrogen bonds [242] but are analogous to hydrogen bonds [243]. Novick et al. [244] stated that halogen bonds are 'anti-hydrogen bonds' [245] while Alkorta et al. [245], Desiraju and Steiner [246] called halogen bonds as 'inverse hydrogen bonds' [245,246] to emphasize the difference between the hydrogen and halogen bonding. Since the halogen atom attached to an electronegative atom (a molecule with 'R' group) contains both nucleophilic and electrophilic end, it can interact with neighboring molecule to form non-covalent bonding, in turn contributing to the building of a new supramolecular structure (as shown in Figure 15).

Ionic Interactions
Ionic cocrystals are another interesting class of cocrystals. These cocrystals are formed out of organic molecules and inorganic salts [248,249]. Braga and his coworkers [250] first coined the term 'ionic cocrystals'. Ionic interactions along with non-covalent interactions such as hydrogen bonding and dipole-bonding interactions stabilize their crystal structure [248]. However, United States Food and Drug Administration (US FDA) and EMA has declared that pharmaceutical cocrystals are crystalline materials composed of two or more components interconnected by non-ionic interactions in a crystal lattice [48][49][50]. Hence, it appears that ionic cocrystals are exception to the definition given by US FDA and EMA on pharmaceutical cocrystals. Ionic cocrystals are also referred as 'Salt cocrystals' when the organic molecules cocrystallize with organic salts [249]. The salt entity in ionic cocrystals can be either organic [57] or inorganic [248,250,251]. The schematic shown in Figure 17 explains a difference between an organic cocrystal and ionic cocrystal [252].

Ionic Interactions
Ionic cocrystals are another interesting class of cocrystals. These cocrystals are formed out of organic molecules and inorganic salts [248,249]. Braga and his coworkers [250] first coined the term 'ionic cocrystals'. Ionic interactions along with non-covalent interactions such as hydrogen bonding and dipole-bonding interactions stabilize their crystal structure [248]. However, United States Food and Drug Administration (US FDA) and EMA has declared that pharmaceutical cocrystals are crystalline materials composed of two or more components interconnected by non-ionic interactions in a crystal lattice [48][49][50]. Hence, it appears that ionic cocrystals are exception to the definition given by US FDA and EMA on pharmaceutical cocrystals. Ionic cocrystals are also referred as 'Salt cocrystals' when the organic molecules cocrystallize with organic salts [249]. The salt entity in ionic cocrystals can be either organic [57] or inorganic [248,250,251]. The schematic shown in Figure 17 explains a difference between an organic cocrystal and ionic cocrystal [252]. crystalline materials composed of two or more components interconnected by non-ionic interactions in a crystal lattice [48][49][50]. Hence, it appears that ionic cocrystals are exception to the definition given by US FDA and EMA on pharmaceutical cocrystals. Ionic cocrystals are also referred as 'Salt cocrystals' when the organic molecules cocrystallize with organic salts [249]. The salt entity in ionic cocrystals can be either organic [57] or inorganic [248,250,251]. The schematic shown in Figure 17 explains a difference between an organic cocrystal and ionic cocrystal [252].  Recently, Song et al. [223] synthesized ionic cocrystal of piracetam with calcium chloride by optimizing solvent parameters for solution crystallization through construction of Ternary Phase Diagram. The resultant ionic cocrystal had molecular formula of piracetam 2 ·CaCl 2 ·2H 2 O which is the piracetam dihydrate ionic cocrystal [223]. The coordination around Ca 2+ cation in piracetam CaCl 2 2H 2 O ionic cocrystal and the 2D layer of Piracetam CaCl 2 2H 2 O ionic cocrystal have been shown in Figure 18 [223]. The piracetam CaCl 2 .2H 2 O crystallized in in a monoclinic space group, P2 1 /n. A unit cell of the ionic cocrystal consisted of two molecules of piracetam, one CaCl 2 molecule and two H 2 O molecules [223]. Table 10 presents the summary of ionic cocrystals reported till date in the literature. Recently, Song et al. [223] synthesized ionic cocrystal of piracetam with calcium chloride by optimizing solvent parameters for solution crystallization through construction of Ternary Phase Diagram. The resultant ionic cocrystal had molecular formula of piracetam2·CaCl2·2H2O which is the piracetam dihydrate ionic cocrystal [223]. The coordination around Ca 2+ cation in piracetam CaCl2 2H2O ionic cocrystal and the 2D layer of Piracetam CaCl2 2H2O ionic cocrystal have been shown in Figure 18 [223]. The piracetam CaCl2.2H2O crystallized in in a monoclinic space group, P21/n. A unit cell of the ionic cocrystal consisted of two molecules of piracetam, one CaCl2 molecule and two H2O molecules [223]. Table 10 presents the summary of ionic cocrystals reported till date in the literature.   There was no significant difference in the IDR of raw PIR and the ionic cocrystal [251] Trimesic acid (H 3 TMA) 2,6-bis(4-pyridylmethylene)cyclohexanone Nil Crystallization Not reported [254] Carbamazepine (CBZ) Sodium iodide (NaI) Acetyl chloride Crystallization in methanol Not reported [255] Sodium iodide (NaI) and Hydrobromic acid (HBr) Acridinium I 2 X species Crystallization in methanol Not reported [255] Benzoic acid (BA) Phenoxy acetic acid Nil Slow evaporation technique Not reported [249] 7.4. Π. . . Π Stacking Interactions Π. . . Π stacking interactions are a type of attractive non-covalent interactions which occur between aromatic rings (as they contain Π bonds). The probability of occurrence of Π. . . Π stacking interactions increases when the number of hydrogen poor aromatic residues is more in a molecule [256]. Several reports are available in the literature where the role of hydrogen bonding in altering physiochemical properties of an API has been reported [15,95]. However, very limited reports are available in the literature where Π. . . Π stacking interactions were utilized for modulating the physiochemical properties of pharmaceutical cocrystals [17,257,258].
Recently, Bora et al. [17] reported that Π. . . Π interactions (with T-shaped motif structures) and C−H· · · π interactions can modulate the physiochemical properties of cocrystals such as pH-dependent solubility, crystal packing and permeability while exploring cocrystallization of acridine with isomeric hydroxybenzoic acids. Since acridine actively participates in intermolecular hydrogen bonding and Π-stacking interactions (namely face-to-face π-stacking interactions and edge to face T-shaped C-H. . . Π interactions), it was chosen for their study [17]. It was concluded that packing of the cocrystal lattice stabilized by π-stacking interactions and C-H. . . Π interactions affects the solubility and membrane permeability of cocrystals. Figure 19 provides a pictorial representation of π-stacking and C-H. . . Π interactions observed in these cocrystal structures.
Sangtani et al. [257] observed concomitant color polymorphism in Furosemide-4,4 -bipyridine (2:1) cocrystals (Forms I and II). Form I cocrystal appeared as pale yellow needles whereas form II cocrystal appeared as orange blocks [257]. Form I cocrystal was observed to crystallize at faster rate and with higher yield than the form II cocrystal. The authors proposed that the color polymorphism in furosemide-4,4 -bipyridine (2:1) cocrystals can be attributed to the difference in Π. . . Π * separation (as shown in Figure 20) between the benzene ring of furosemide and pyridine ring of 4,4 -bipyridine in the sandwich motifs of form I and form II cocrystals [257]. This was further supported by Dynamic Functional Theory (DFT) calculations. From DFT calculations, it was evident that the Highest Occupied Molecular Orbital (HOMO)-Lowest Unoccupied Molecular Orbital (LUMO) gap for form I pale yellow needles was more when compared with form II orange blocks [257].

Π…Π Stacking Interactions
Π…Π stacking interactions are a type of attractive non-covalent interactions which occur between aromatic rings (as they contain Π bonds). The probability of occurrence of Π…Π stacking interactions increases when the number of hydrogen poor aromatic residues is more in a molecule [256]. Several reports are available in the literature where the role of hydrogen bonding in altering physiochemical properties of an API has been reported [15,95]. However, very limited reports are available in the literature where Π…Π stacking interactions were utilized for modulating the physiochemical properties of pharmaceutical cocrystals [17,257,258].
Recently, Bora et al. [17] reported that Π…Π interactions (with T-shaped motif structures) and C−H⋯π interactions can modulate the physiochemical properties of cocrystals such as pH-dependent solubility, crystal packing and permeability while exploring cocrystallization of acridine with isomeric hydroxybenzoic acids. Since acridine actively participates in intermolecular hydrogen bonding and Π-stacking interactions (namely face-to-face π-stacking interactions and edge to face T-shaped C-H…Π interactions), it was chosen for their study [17]. It was concluded that packing of the cocrystal lattice stabilized by π-stacking interactions and C-H…Π interactions affects the solubility and membrane permeability of cocrystals. Figure 19 provides a pictorial representation of π-stacking and C-H…Π interactions observed in these cocrystal structures.
Sangtani et al. [257] observed concomitant color polymorphism in Furosemide-4,4′-bipyridine (2:1) cocrystals (Forms I and II). Form I cocrystal appeared as pale yellow needles whereas form II cocrystal appeared as orange blocks [257]. Form I cocrystal was observed to crystallize at faster rate and with higher yield than the form II cocrystal. The authors proposed that the color polymorphism in furosemide-4,4′-bipyridine (2:1) cocrystals can be attributed to the difference in Π…Π * separation (as shown in Figure 20) between the benzene ring of furosemide and pyridine ring of 4,4′-bipyridine in the sandwich motifs of form I and form II cocrystals [257]. This was further supported by Dynamic Functional Theory (DFT) calculations. From DFT calculations, it was evident that the Highest Occupied Molecular Orbital (HOMO)-Lowest Unoccupied Molecular Orbital (LUMO) gap for form I pale yellow needles was more when compared with form II orange blocks [257].  Thus, Π…Π interactions can significantly contribute to modification of physiochemical properties of an API molecule. Further research is needed to explore the role of Π…Π interactions in fine-tuning the solid-state properties (especially the dissolution rate, solubility, physical and chemical stability) of poorly water-soluble drugs.

Vanderwaal's Interactions
Vanderwaal's interactions are yet other type of intermolecular interactions which contribute to the packing of the cocrystal lattice [259]. It is the weakest of all intermolecular interactions. However, collectively these interactions are strong. These interactions were also observed to influence the formation of coamorphous solid dispersions to a certain extent [260].

Ternary and Quaternary Cocrystals
Ternary cocrystals and quaternary cocrystals have recently gained attention from several researchers mainly for poorly water-soluble drugs to enhance dissolution rates and aqueous solubility [261,262]. Ternary (three-component) and quaternary (four-component) cocrystals are supramolecular structures stabilized by robust synthons which build the entire crystal structure. Synthesizing ternary cocrystals is not easy when compared with the synthesis of binary cocrystals as maintenance of congruent solubility of drug and the two coformers is difficult and requires a lot of screening experiments. There is a need to understand the basic principles involved in formation of ternary cocrystals and intermolecular interactions involved in stabilizing their crystal structures [263][264][265]. Several approaches have been used by researchers to design ternary cocrystals based on crystal engineering and synthon engineering principles [265][266][267][268]. Some of the approaches have been explained below: Thus, Π. . . Π interactions can significantly contribute to modification of physiochemical properties of an API molecule. Further research is needed to explore the role of Π. . . Π interactions in fine-tuning the solid-state properties (especially the dissolution rate, solubility, physical and chemical stability) of poorly water-soluble drugs.

Vanderwaal's Interactions
Vanderwaal's interactions are yet other type of intermolecular interactions which contribute to the packing of the cocrystal lattice [259]. It is the weakest of all intermolecular interactions. However, collectively these interactions are strong. These interactions were also observed to influence the formation of coamorphous solid dispersions to a certain extent [260].

Ternary and Quaternary Cocrystals
Ternary cocrystals and quaternary cocrystals have recently gained attention from several researchers mainly for poorly water-soluble drugs to enhance dissolution rates and aqueous solubility [261,262]. Ternary (three-component) and quaternary (four-component) cocrystals are supramolecular structures stabilized by robust synthons which build the entire crystal structure. Synthesizing ternary cocrystals is not easy when compared with the synthesis of binary cocrystals as maintenance of congruent solubility of drug and the two coformers is difficult and requires a lot of screening experiments. There is a need to understand the basic principles involved in formation of ternary cocrystals and intermolecular interactions involved in stabilizing their crystal structures [263][264][265]. Several approaches have been used by researchers to design ternary cocrystals based on crystal engineering and synthon engineering principles [265][266][267][268]. Some of the approaches have been explained below:

A. Hierarchical Fashion of Hydrogen Bond Formation
Aakeroy et al. [266] developed an approach to design ternary cocrystals based on the two principle rules: a.
Hydrogen bond formation take place in a hierarchical manner (best donor forms hydrogen with the best acceptor, the second-best donor with the second-best acceptor) [103] b.

B. Long-range Synthon Aufbau Module (LSAM)
LSAM has its origin from the 'Aufbau principle' proposed by Kitaigorodski in the year 1961 [270]. According to LSAM, the molecules which are closely-packed in one-dimensional chains, arrange themselves to form two-dimensional sheets. These two-dimensional sheets assemble and pack closely to form three-dimensional supramolecular structural crystal units. Reports are available in the literature where researchers have reported formation of binary (two-component) cocrystals [265], ternary (three-component) cocrystals [265,267,271,272] and quaternary (four-component) cocrystals [267] based on LSAM.

B. Long-range Synthon Aufbau Module (LSAM)
LSAM has its origin from the 'Aufbau principle' proposed by Kitaigorodski in the year 1961 [270]. According to LSAM, the molecules which are closely-packed in one-dimensional chains, arrange themselves to form two-dimensional sheets. These two-dimensional sheets assemble and pack closely to form three-dimensional supramolecular structural crystal units. Reports are available in the literature where researchers have reported formation of binary (two-component) cocrystals [265], ternary (three-component) cocrystals [265,267,271,272] and quaternary (four-component) cocrystals [267] based on LSAM.

C. Drug-Bridge-Drug Ternary Cocrystallization Strategy
Recently, Liu et al. [261] developed a 'Drug-Bridge-Drug Ternary Cocrystallization strategy' for cocrystallizing two anti-tuberculosis drugs, isoniazid and pyrazinamide using trans-fumaric acid as a bridge. Pyrazinamide has aqueous solubility of 15,000 mg/L [273] whereas isoniazid has aqueous solubility of 140,000 mg/L at 25 °C [274]. Figure 23 presents (a) the chemical structures of the API molecules, isoniazid and pyrazinamide, trans and cis fumaric acid and (b) the schematic how fumaric acid act as a bridge in interconnecting isoniazid and pyrazinamide molecules.
Fumaric acid was chosen as a bridge to link isoniazid and pyrazinamide as the carboxylic acid functional groups in fumaric acid can interact effectively with N-heterocycles of isoniazid and pyrazinamide to form stronger heterosynthon [261]. Moreover, the amide functional group of pyrazinamide and hydrazide functional group in isoniazid has been reported to form hydrogen-bonding motifs with dicarboxylic acids and tricarboxylic acids in several studies [275,276]. The selection of suitable dicarboxylic acid plays an important role in determining the success of this 'drug-bridge-drug ternary cocrystallization strategy'.  Long-range Synthon Assembly Module (LSAM) observed in acetazolomidenicotinamide-hydroxypyridine (1:1:1) ternary cocrystal [Reprinted from [267] with permission of International Union of Crystallography].

C. Drug-Bridge-Drug Ternary Cocrystallization Strategy
Recently, Liu et al. [261] developed a 'Drug-Bridge-Drug Ternary Cocrystallization strategy' for cocrystallizing two anti-tuberculosis drugs, isoniazid and pyrazinamide using trans-fumaric acid as a bridge. Pyrazinamide has aqueous solubility of 15,000 mg/L [273] whereas isoniazid has aqueous solubility of 140,000 mg/L at 25 • C [274]. Figure 23 presents (a) the chemical structures of the API molecules, isoniazid and pyrazinamide, trans and cis fumaric acid and (b) the schematic how fumaric acid act as a bridge in interconnecting isoniazid and pyrazinamide molecules.
Fumaric acid was chosen as a bridge to link isoniazid and pyrazinamide as the carboxylic acid functional groups in fumaric acid can interact effectively with N-heterocycles of isoniazid and pyrazinamide to form stronger heterosynthon [261]. Moreover, the amide functional group of pyrazinamide and hydrazide functional group in isoniazid has been reported to form hydrogen-bonding motifs with dicarboxylic acids and tricarboxylic acids in several studies [275,276]. The selection of suitable dicarboxylic acid plays an important role in determining the success of this 'drug-bridge-drug ternary cocrystallization strategy'.

C. Drug-Bridge-Drug Ternary Cocrystallization Strategy
Recently, Liu et al. [261] developed a 'Drug-Bridge-Drug Ternary Cocrystallization strategy' for cocrystallizing two anti-tuberculosis drugs, isoniazid and pyrazinamide using trans-fumaric acid as a bridge. Pyrazinamide has aqueous solubility of 15,000 mg/L [273] whereas isoniazid has aqueous solubility of 140,000 mg/L at 25 °C [274]. Figure 23 presents (a) the chemical structures of the API molecules, isoniazid and pyrazinamide, trans and cis fumaric acid and (b) the schematic how fumaric acid act as a bridge in interconnecting isoniazid and pyrazinamide molecules.
Fumaric acid was chosen as a bridge to link isoniazid and pyrazinamide as the carboxylic acid functional groups in fumaric acid can interact effectively with N-heterocycles of isoniazid and pyrazinamide to form stronger heterosynthon [261]. Moreover, the amide functional group of pyrazinamide and hydrazide functional group in isoniazid has been reported to form hydrogen-bonding motifs with dicarboxylic acids and tricarboxylic acids in several studies [275,276]. The selection of suitable dicarboxylic acid plays an important role in determining the success of this 'drug-bridge-drug ternary cocrystallization strategy'.

D. Ditopic Hydrogen Bond Donors and Acceptors Combination
Etter's rule of hydrogen bonding [103] suggests that design of a cocrystal mainly depends on hydrogen bond acceptors and hydrogen bond donors which participate in cocrystallization (discussed in Section 3.2). In general, ternary cocrystals are stabilized by strong intermolecular interactions that exist between three components. Later, Aakeroy et al. [266] proposed that these intermolecular hydrogen-bonding interactions in ternary cocrystals can occur as per any one of the patterns illustrated in Figure 24.

D. Ditopic Hydrogen Bond Donors and Acceptors Combination
Etter's rule of hydrogen bonding [103] suggests that design of a cocrystal mainly depends on hydrogen bond acceptors and hydrogen bond donors which participate in cocrystallization (discussed in Section 3.2). In general, ternary cocrystals are stabilized by strong intermolecular interactions that exist between three components. Later, Aakeroy et al. [266] proposed that these intermolecular hydrogen-bonding interactions in ternary cocrystals can occur as per any one of the patterns illustrated in Figure 24. Adsmond et al. [264] adopted strategy 2 (shown in Figure 25) to design ternary cocrystals [264] of acridine, 3-hydroxybenzoic acid and 2-amino-4,6-dimethylpyridine using carboxyphenol (3-hydroxybenzoic acid) as a ditopic hydrogen donors molecule. According to the proposed strategy, a molecule with ditopic hydrogenbond donors (such as carboxy phenols) can be combined with two different molecules containing hydrogenbond acceptors in them to generate a new three-component crystal [264]. The schematic representation of the strategy used by Adsmond et al. [264] in formation of acridine-3-hydroxybenzoic acid-2-amino-4,6-dimethylpyridine (1:1:1) cocrystal is shown in Figure 25. Also, Bhogala and Nangia [277] synthesized ternary and quaternary cocrystals of 1, 3cis,5cis-cyclohexanetricarboxylic acids with 4,4′-bipyridine bases connected by methylene and alkene chains [277]. It was concluded from the study that utilization of triacid and differentiated bipyridine bases is responsible for the self-assembly of the different components through O-H…N hydrogen bonds in the ternary and quaternary systems. Tilborg et al. [278] reported a three-component cocrystal of 1:1:1 L-proline-D-proline-fumaric acid with a racemic compound, DL-proline and fumaric acid [278]. Interestingly, a 1:1 eutectic mixture of Pyrazinamide-Isoniazid Adsmond et al. [264] adopted strategy 2 (shown in Figure 25) to design ternary cocrystals [264] of acridine, 3-hydroxybenzoic acid and 2-amino-4,6-dimethylpyridine using carboxyphenol (3-hydroxybenzoic acid) as a ditopic hydrogen donors molecule. According to the proposed strategy, a molecule with ditopic hydrogenbond donors (such as carboxy phenols) can be combined with two different molecules containing hydrogenbond acceptors in them to generate a new three-component crystal [264]. The schematic representation of the strategy used by Adsmond et al. [264] in formation of acridine-3-hydroxybenzoic acid-2-amino-4,6-dimethylpyridine (1:1:1) cocrystal is shown in Figure 25.

D. Ditopic Hydrogen Bond Donors and Acceptors Combination
Etter's rule of hydrogen bonding [103] suggests that design of a cocrystal mainly depends on hydrogen bond acceptors and hydrogen bond donors which participate in cocrystallization (discussed in Section 3.2). In general, ternary cocrystals are stabilized by strong intermolecular interactions that exist between three components. Later, Aakeroy et al. [266] proposed that these intermolecular hydrogen-bonding interactions in ternary cocrystals can occur as per any one of the patterns illustrated in Figure 24. Adsmond et al. [264] adopted strategy 2 (shown in Figure 25) to design ternary cocrystals [264] of acridine, 3-hydroxybenzoic acid and 2-amino-4,6-dimethylpyridine using carboxyphenol (3-hydroxybenzoic acid) as a ditopic hydrogen donors molecule. According to the proposed strategy, a molecule with ditopic hydrogenbond donors (such as carboxy phenols) can be combined with two different molecules containing hydrogenbond acceptors in them to generate a new three-component crystal [264]. The schematic representation of the strategy used by Adsmond et al. [264] in formation of acridine-3-hydroxybenzoic acid-2-amino-4,6-dimethylpyridine (1:1:1) cocrystal is shown in Figure 25. Also, Bhogala and Nangia [277] synthesized ternary and quaternary cocrystals of 1, 3cis,5cis-cyclohexanetricarboxylic acids with 4,4′-bipyridine bases connected by methylene and alkene chains [277]. It was concluded from the study that utilization of triacid and differentiated bipyridine bases is responsible for the self-assembly of the different components through O-H…N hydrogen bonds in the ternary and quaternary systems. Tilborg et al. [278] reported a three-component cocrystal of 1:1:1 L-proline-D-proline-fumaric acid with a racemic compound, DL-proline and fumaric acid [278]. Interestingly, a 1:1 eutectic mixture of Pyrazinamide-Isoniazid Also, Bhogala and Nangia [277] synthesized ternary and quaternary cocrystals of 1, 3cis,5cis-cyclohexanetricarboxylic acids with 4,4 -bipyridine bases connected by methylene and alkene chains [277]. It was concluded from the study that utilization of triacid and differentiated bipyridine bases is responsible for the self-assembly of the different components through O-H. . . N hydrogen bonds in the ternary and quaternary systems. Tilborg et al. [278] reported a three-component cocrystal of 1:1:1 L-proline-D-proline-fumaric acid with a racemic compound, DL-proline and fumaric acid [278]. Interestingly, a 1:1 eutectic mixture of Pyrazinamide-Isoniazid [18] yielded two different ternary cocrystals with succinic acid and fumaric acid as coformers (each with a stoichiometric ratio of 1:1:1) [23]. Also, Cheung et al. [279] reported that solid-state 13 C NMR can serve as a direct method to determine the structure of ternary cocrystal of benzoquinone (BQ), racemic bis-β-Naphthol (BN) and anthracene (AN) in a stoichiometric ratio of 1:1:0.5, prepared by solid-state grinding [279]. Recently, ternary cocrystals of Bumetanide with isonicotinamide-2-picolinic acid, isonicotinamide-vanillic acid, isonicotinamide-para aminosalicylic acid and 2-hydroxypyridone-2-picolinic acid pairs were reported by Allu et al. [280]. Aitipamula et al. [262] prepared ternary cocrystals of isoniazid with enhanced dissolution rates as compared to the commercially available isoniazid [262].

Advantages of Ternary/Quaternary Cocrystals
Ternary/quaternary cocrystals are highly useful in enhancing the efficacy of drugs when it is not achieved by formulating binary cocrystals. Furthermore, when it is not possible to cocrystallize two different APIs (to enhance their efficacy) due to their shape incompatibility; another molecule(s) can be introduced between the two drugs which can serve as a linker, propagate the growth unit and enable cocrystallization, and form a supramolecular structure [261]. Like binary cocrystals, three-or four-component cocrystals enhance dissolution rate [261,262], bioavailability and physical stability of poorly water-soluble drugs

Disadvantages of Ternary/Quaternary Cocrystals
Though ternary/quaternary cocrystals possess advantages such as enhanced bioavailability, solubility, dissolution rates and physical stability, the synthesis of ternary cocrystals is difficult as the intermolecular interactions in the ternary cocrystals should be very specific and balanced. Therefore, the synthesis of ternary cocrystals involves a lot of efforts in screening and synthesis and poses difficulty in terms of understanding of the chemistry.

Polymorphism in Cocrystals
Polymorphism can be defined as the ability of a drug to exists in more than one crystalline phases with variation in arrangements or variation in the conformation of drug molecules in a crystal lattice [102,281,282]. Drugs existing in more than one crystal form can have different physicochemical properties such as dissolution rate, solubility, morphology, mechanical properties and physicochemical stability [283,284]. Polymorphism in cocrystals can influence aqueous solubility, dissolution property and bioavailability of a drug [35].
Cocrystals exhibit different types of polymorphic behaviors such as synthon polymorphism [285,286], concomitant polymorphism [257,287,288], conformational polymorphism, packing polymorphism and pseudopolymorphism [285]. Each of these types has been discussed in brief in sections below:

Concomitant Polymorphism
Concomitant polymorphs are polymorphs which crystallize simultaneously in identical crystallization conditions in the same batch. Reports are available in the literature where researchers have reported formation of concomitant cocrystal polymorphs [285,287,288,[300][301][302]. Color polymorphism is another interesting variation of concomitant polymorphism in which cocrystal forms of more than one color appear during crystallization. A few findings (See Section 7.4) on colored polymorphism of cocrystal have been reported by Gonnade and coworkers [257,258].

Concomitant Polymorphism
Concomitant polymorphs are polymorphs which crystallize simultaneously in identical crystallization conditions in the same batch. Reports are available in the literature where researchers have reported formation of concomitant cocrystal polymorphs [285,287,288,[300][301][302]. Color polymorphism is another interesting variation of concomitant polymorphism in which cocrystal forms of more than one color appear during crystallization. A few findings (See Section 7.4) on colored polymorphism of cocrystal have been reported by Gonnade and coworkers [257,258].
Tothadi [306] identified an interesting packing polymorphism in urea-4,4 -bipyridine (1:1) cocrystals [306]. Figure 29 presents the packing polymorphs (Form I and II) of urea-4,4 -bipyridine (1:1) cocrystals. Occurrence of urea tape in the crystal structure of urea-containing cocrystals is one of the rare observations in crystallography. Till date, out of 194 urea-containing cocrystal structures only 5 cocrystal structures were identified to possess urea tape structures [308]. Interestingly, both the polymorphic forms contained urea tapes in their crystal structures [306]. Though the urea tape arrangements are common in both the polymorphs, the 3D packing significantly differed from each other. In the Form 1 cocrystal, the consecutive planes are arranged in an ABA'B ABA B fashion whereas in the Form 2 cocrystal, the consecutive planes are arranged in an ABAB fashion.
of the rare observations in crystallography. Till date, out of 194 urea-containing cocrystal structures only 5 cocrystal structures were identified to possess urea tape structures [308]. Interestingly, both the polymorphic forms contained urea tapes in their crystal structures [306]. Though the urea tape arrangements are common in both the polymorphs, the 3D packing significantly differed from each other. In the Form 1 cocrystal, the consecutive planes are arranged in an ABA'B′ABA′B′ fashion whereas in the Form 2 cocrystal, the consecutive planes are arranged in an ABAB fashion. Recently, Surov et al. [307] addressed that weak interactions can cause packing polymorphism in pharmaceutical cocrystals while investigating packing polymorphism (form I and II) in salicylamide-oxalic acid (2:1) cocrystal [307]. The crystal structures of form I and form II cocrystal polymorphs revealed that both the polymorphs consisted of conformationally identical molecules. However, packing of the two polymorphs differed in terms of arrangement of the neighboring ribbons [307]. In form I cocrystal, the adjusted ribbons were packed in zig-zag fashion at an angle of ~76.2° while it is ~48.6° in case of form II cocrystal. Crystal structure analysis and Hirshfeld surface analysis revealed that the distribution of weak intermolecular interactions in polymorphic forms influenced the packing of the cocrystal lattice to a larger extent [307]. Figure 30 represents the molecular packing projections of salicylamide-oxalic acid (2:1) form I and form II cocrystal along an axis.  Recently, Surov et al. [307] addressed that weak interactions can cause packing polymorphism in pharmaceutical cocrystals while investigating packing polymorphism (form I and II) in salicylamide-oxalic acid (2:1) cocrystal [307]. The crystal structures of form I and form II cocrystal polymorphs revealed that both the polymorphs consisted of conformationally identical molecules. However, packing of the two polymorphs differed in terms of arrangement of the neighboring ribbons [307]. In form I cocrystal, the adjusted ribbons were packed in zig-zag fashion at an angle of~76.2 • while it is~48.6 • in case of form II cocrystal. Crystal structure analysis and Hirshfeld surface analysis revealed that the distribution of weak intermolecular interactions in polymorphic forms influenced the packing of the cocrystal lattice to a larger extent [307]. Figure 30 represents the molecular packing projections of salicylamide-oxalic acid (2:1) form I and form II cocrystal along an axis.
other. In the Form 1 cocrystal, the consecutive planes are arranged in an ABA'B′ABA′B′ fashion whereas in the Form 2 cocrystal, the consecutive planes are arranged in an ABAB fashion. Recently, Surov et al. [307] addressed that weak interactions can cause packing polymorphism in pharmaceutical cocrystals while investigating packing polymorphism (form I and II) in salicylamide-oxalic acid (2:1) cocrystal [307]. The crystal structures of form I and form II cocrystal polymorphs revealed that both the polymorphs consisted of conformationally identical molecules. However, packing of the two polymorphs differed in terms of arrangement of the neighboring ribbons [307]. In form I cocrystal, the adjusted ribbons were packed in zig-zag fashion at an angle of ~76.2° while it is ~48.6° in case of form II cocrystal. Crystal structure analysis and Hirshfeld surface analysis revealed that the distribution of weak intermolecular interactions in polymorphic forms influenced the packing of the cocrystal lattice to a larger extent [307]. Figure 30 represents the molecular packing projections of salicylamide-oxalic acid (2:1) form I and form II cocrystal along an axis.
Row and coworkers [288] reported gallic acid-succinimide cocrystal solvates/hydrates of various stoichiometric ratios formed with different types of solvents (such as 1,4-dioxane, water, tetrahydrofuran, ethyl acetate, acetone) and water inclusions as pseudopolymorphs in their study [288]. Figures 31 and 32 present the staircase network of gallic acid-succinimide-tetrahydrofuran (2:2:1) cocrystal and gallic acid-succinimide-water (1:1:1) cocrystal made up of carboxylic acid and carboxamide homodimers with tetrahydrofuran and water molecules as inclusions in it. The solvent molecules occupied the voids in cocrystal lattice through C-H…O interactions (in case of tetrahydrofuran) and through O-H…O interactions (in case of water) [288].

Impact of Cocrystal Polymorphism on Solid-State Properties of API
Similar to the single component solids, polymorphism in multicomponent solids also plays a significant role in determining quality, efficacy and safety of an API molecule. Polymorphism has been known to influence aqueous solubility [310] and physical stability [56,230,262,310,311] of cocrystals. Paradkar and coworkers [310] synthesized dimorphic forms of carbamazepine-saccharin (1:1) cocrystal (Form I and Form II) by slow evaporative solution crystallization method and characterized their thermodynamic interrelationship [310]. During the study, it was observed from the van 't Hoff plot and DSC thermograms that Form I carbamazepine-saccharin (1:1) remains as a stable form whereas the Form II carbamazepine-saccharin (1:1) exists as a metastable form.
Row and coworkers [288] reported gallic acid-succinimide cocrystal solvates/hydrates of various stoichiometric ratios formed with different types of solvents (such as 1,4-dioxane, water, tetrahydrofuran, ethyl acetate, acetone) and water inclusions as pseudopolymorphs in their study [288]. Figures 31 and 32 present the staircase network of gallic acid-succinimide-tetrahydrofuran (2:2:1) cocrystal and gallic acid-succinimide-water (1:1:1) cocrystal made up of carboxylic acid and carboxamide homodimers with tetrahydrofuran and water molecules as inclusions in it. The solvent molecules occupied the voids in cocrystal lattice through C-H…O interactions (in case of tetrahydrofuran) and through O-H…O interactions (in case of water) [288].

Impact of Cocrystal Polymorphism on Solid-State Properties of API
Similar to the single component solids, polymorphism in multicomponent solids also plays a significant role in determining quality, efficacy and safety of an API molecule. Polymorphism has been known to influence aqueous solubility [310] and physical stability [56,230,262,310,311] of cocrystals. Paradkar and coworkers [310] synthesized dimorphic forms of carbamazepine-saccharin (1:1) cocrystal (Form I and Form II) by slow evaporative solution crystallization method and characterized their thermodynamic interrelationship [310]. During the study, it was observed from the van 't Hoff plot and DSC thermograms that Form I carbamazepine-saccharin (1:1) remains as a stable form whereas the Form II carbamazepine-saccharin (1:1) exists as a metastable form.

Impact of Cocrystal Polymorphism on Solid-State Properties of API
Similar to the single component solids, polymorphism in multicomponent solids also plays a significant role in determining quality, efficacy and safety of an API molecule. Polymorphism has been known to influence aqueous solubility [310] and physical stability [56,230,262,310,311] of cocrystals. Paradkar and coworkers [310] synthesized dimorphic forms of carbamazepine-saccharin (1:1) cocrystal (Form I and Form II) by slow evaporative solution crystallization method and characterized their thermodynamic interrelationship [310]. During the study, it was observed from the van 't Hoff plot and DSC thermograms that Form I carbamazepine-saccharin (1:1) remains as a stable form whereas the Form II carbamazepine-saccharin (1:1) exists as a metastable form. Moreover, Form II polymorph showed enhanced solubility than form I polymorph with respect to different temperatures in deionized water owing to its less stability than form I [310].

Solubility and Dissolution Enhancement by Cocrystals
Bioavailability of a cocrystal is determined by its dissolution rate, aqueous solubility and permeability. Solubility is defined as a maximum amount of drug that can be dissolved in a solution and is nothing but the thermodynamic equilibrium attained by a solute between a solid and liquid phase [312]. Cocrystal dissolution refers to the amount of solute which dissolves in an aqueous medium at a specific time interval. The rate at which the solute dissolves in the aqueous medium to reach its equilibrium state is called as the dissolution rate [313]. Hence, solubility is a thermodynamic process whereas dissolution is a kinetic process. Given below is a brief account of information available on solubility and dissolution of cocrystals.

Cocrystal Solubilization
Many researchers have attempted to understand the mechanism behind the solubility of cocrystals [9,95,[314][315][316]. The insight obtained from the review of such literature reports indicates that the solubility of cocrystals depends on two important parameters namely the strength of intermolecular interactions in the crystal lattice and solvation of cocrystal components [9,10,95,[314][315][316]. Maheshwari et al. [316] reported that solubilization of a cocrystal in dissolution medium involves two main steps: (1) Release of the solute molecules from the crystal lattice of the cocrystal and (2) the solvation of the released molecules [316] (as shown in Figure 33). Therefore, free energy of cocrystal solubilization depends on the free energy associated with the release of solute molecules from the cocrystal lattice and the free energy associated with the solvation barrier of the cocrystal as given in Equation (1) [316], where ∆G solution is the Gibb's free energy associated with the solubilization process, ∆G lattice is the Gibb's free energy associated with the cocrystal lattice and ∆G solvation is the Gibb's free energy associated with the solvation barrier. When the free energy associated with the lattice interactions and the free energy associated with the solvation barrier becomes negligible, enhancement in cocrystal dissolution is achieved due to decrease in free energy change for solubilization [316]. The cocrystal formed with a highly water-soluble coformer, was found to possess high solubility [9,64]. Therefore, Maheshwari et al. [316] suggested that solvation is the most important barrier for the solubilization of cocrystals for hydrophobic drugs [316].

Phase Solubility Diagram (PSD)
PSDs are used to identify the regions for drug and conformer stability in terms of drug and conformer concentration [9]. For a cocrystal [A α B β ], drug (A) concentration can be expressed as a function of conformer (B) concentration using following equation, where K sp is the solubility product for dissolution of a cocrystal. Figure 34 is a typical PSD for two different cocrystals namely stable and metastable cocrystal [9]. In this figure, the dotted line represents stoichiometric concentrations of the drug and coformer and its intersection with the cocrystal solubility curves (denoted by filled circles) gives the maximum drug concentration corresponding to the cocrystal solubility. It can be observed from Figure 34 that the stable cocrystal exhibits lower solubility than the drug solubility in a given solvent/medium. On the other hand, a metastable cocrystal exhibits solubility higher than the drug solubility [9]. The invariant points (characterized by zero degrees of freedom and represented by x marks) shown in this figure indicate the transition points or eutectic points where solid drug, cocrystal and solution containing drug and coformer are in equilibrium [9]. Furthermore, it can be noted from Figure 34 that the stoichiometric composition of the cocrystal affects the cocrystal solubility.

Factors Influencing Solubility of Cocrystals
During cocrystal dissolution in the aqueous medium, factors such as ionization of parent components [317], pH of the aqueous medium [317], drug-solubilizing agents [317] and coformer concentration [317] influence the solubility of cocrystals.

(a) Ionization of Parent Components and pH of Aqueous Medium
Cocrystal can be synthesized by cocrystallizing a neutral drug molecule with acidic or amphoteric coformer molecules or a zwitterionic drug with acidic coformer or a basic drug with

Factors Influencing Solubility of Cocrystals
During cocrystal dissolution in the aqueous medium, factors such as ionization of parent components [317], pH of the aqueous medium [317], drug-solubilizing agents [317] and coformer concentration [317] influence the solubility of cocrystals.

(a) Ionization of Parent Components and pH of Aqueous Medium
Cocrystal can be synthesized by cocrystallizing a neutral drug molecule with acidic or amphoteric coformer molecules or a zwitterionic drug with acidic coformer or a basic drug with

Factors Influencing Solubility of Cocrystals
During cocrystal dissolution in the aqueous medium, factors such as ionization of parent components [317], pH of the aqueous medium [317], drug-solubilizing agents [317] and coformer concentration [317] influence the solubility of cocrystals.

(a) Ionization of Parent Components and pH of Aqueous Medium
Cocrystal can be synthesized by cocrystallizing a neutral drug molecule with acidic or amphoteric coformer molecules or a zwitterionic drug with acidic coformer or a basic drug with acidic coformer molecule. As a result, the properties of a cocrystal designed for the same API but with different coformers also exhibit different solubility behavior mainly due to variation in ionization behavior of coformer molecules. Since ionization is highly dependent on pH, pH is another parameter which significantly affects the cocrystal solubility.
Maheshwari et al. [316] used Gabapentin-lactam (GBPL) as a model drug and synthesized its cocrystals with carboxylic acids, fumaric acid (FA), 4-hydrobenzoic acid (4HBA), genitisic acid (GA), 4-aminobenzoic acid (4ABA) and benzoic acid (BA) as coformers by reaction crystallization. The aim was to understand the influence of coformers with low aqueous solubility on the solubility of cocrystals [303]. From their study, it was identified that the so formed GBPL 2 -FA (2:1), GBPL-4HBA (1:1), GBPL-GA (1:1), GBPL-4ABA (1:1) and GBPL-BA (1:1) cocrystals showed lower aqueous solubility than raw GBPL [316] and exhibited pH-dependent solubilities. Figure 35 presents variation in solubility of gabapentin-lactam cocrystals as a function of pH. The pH value at which the solubility curve of the cocrystal and drug intersects is called as pH max . Figure 35 shows the pH max values exhibited by five different gabapentin-lactam cocrystals [316]. Below pH max value, the five GBPL cocrystals showed lower solubilitythan the drug. Therefore, this point is called as 'the Eutectic point' or 'transition point' [9,318]. The pH max is also called as Gibb's pH as it defines the thermodynamic stability of the cocrystal and above this pH the cocrystal is unstable [316]. ionization behavior of coformer molecules. Since ionization is highly dependent on pH, pH is another parameter which significantly affects the cocrystal solubility. Maheshwari et al. [316] used Gabapentin-lactam (GBPL) as a model drug and synthesized its cocrystals with carboxylic acids, fumaric acid (FA), 4-hydrobenzoic acid (4HBA), genitisic acid (GA), 4-aminobenzoic acid (4ABA) and benzoic acid (BA) as coformers by reaction crystallization. The aim was to understand the influence of coformers with low aqueous solubility on the solubility of cocrystals [303]. From their study, it was identified that the so formed GBPL2-FA (2:1), GBPL-4HBA (1:1), GBPL-GA (1:1), GBPL-4ABA (1:1) and GBPL-BA (1:1) cocrystals showed lower aqueous solubility than raw GBPL [316] and exhibited pH-dependent solubilities. Figure 35 presents variation in solubility of gabapentin-lactam cocrystals as a function of pH. The pH value at which the solubility curve of the cocrystal and drug intersects is called as pHmax. Figure 35 shows the pHmax values exhibited by five different gabapentin-lactam cocrystals [316]. Below pHmax value, the five GBPL cocrystals showed lower solubilitythan the drug. Therefore, this point is called as 'the Eutectic point' or 'transition point' [9,318]. The pHmax is also called as Gibb's pH as it defines the thermodynamic stability of the cocrystal and above this pH the cocrystal is unstable [316]. Thakuria et al. [314] has provided a detailed review on how the ionization of components affects the pH solubility profile of different drug pairs reported in the literature. Table 11 presents generalized observations reported in the literature about the solubility behavior of cocrystals. Thakuria et al. [314] has provided a detailed review on how the ionization of components affects the pH solubility profile of different drug pairs reported in the literature. Table 11 presents generalized observations reported in the literature about the solubility behavior of cocrystals. Gabapentin-lactam Non-ionizable Gentisic acid Acidic Increase in solubility [316] Gabapentin-lactam Non-ionizable Benzoic acid Acidic Increase in solubility [316] Gabapentin-lactam Non-ionizable 4-aminobenzoic acid Monoprotic acid U-shaped trend in which solubility reached minimum with an increase in pH [316] Gabapentin-lactam Non-ionizable 4-hydroxybenzoic acid Acidic Increase in solubility [316] Gabapentin-lactam Non-ionizable Fumaric acid Diprotic acid Increase in solubility [316] Gabapentin Zwitterionic 3-hydroxybenzoic acid Acidic U-shaped trend leading to increase in solubility [321] (

b) Effect of Solubilizing Agent
Though drug-solubilizing agents such as surfactants can increase the solubility of a cocrystal in an aqueous medium [322], at a certain concentration level it also reduces the solubility of a cocrystal. The concentration of the surfactant [M] therefore determines the solubility and stability of the cocrystals in anaqueous medium. The cocrystal solubility is directly proportional to the square root of [M] whereas the drug solubility is directly proportional to [M] [314]. The concentration of surfactant at which solubility of a cocrystal and drug becomes identical is called as Critical Stabilization Concentration (CSC). The cocrystal becomes thermodynamically unstable below the CSC in the dissolution medium (See Figure 36). On the other hand, it remains in equilibrium with pure drug at CSC (this point is known as 'the eutectic point") and remains thermodynamically stable above CSC. Therefore, CSC serves as an important parameter in influencing the solubility of the cocrystals. Figure 37 shows typical dissolution profiles of insoluble drugs [85] which can be characterized by "spring and parachute model" [323]. Babu and Nangia [85] stated that dissolution of cocrystals can also be explained using 'Spring and parachute model'.

Dissolution of Cocrystals
During dissolution, cocrystals exhibit maximum peak value in the drug concentration which is attained in a shorter time (such as less than 30 min) (spring effect). This high concentration is maintained for a longer period before ultimately decreasing to the equilibrium solubility level (parachute effect) [95]. The spring and parachute effect can be explained using a mechanism proposed by Babu and Nangia [95]. As per the proposed mechanism, three main steps are involved during dissolution of cocrystals (as shown in Figure 38): (i) Dissociation of the cocrystal into amorphous or nanocrystalline drug clusters, which is represented as 'spring' effect ( Figures 37 and 38) (ii) Transformation of the amorphous or nanocrystalline drug clusters into a stable form through formation of a metastable phase by adopting Ostwald's Law of Stages (iii) Attainment of higher apparent solubility and maintenance of optimal drug concentration in the aqueous medium, which is represented as 'parachute' effect ( Figures 37 and 38).
Pharmaceutics 2018, 10, x FOR PEER REVIEW 53 of 74 (iii) Attainment of higher apparent solubility and maintenance of optimal drug concentration in the aqueous medium, which is represented as 'parachute' effect ( Figures 37 and 38).   (iii) Attainment of higher apparent solubility and maintenance of optimal drug concentration in the aqueous medium, which is represented as 'parachute' effect ( Figures 37 and 38).

Factors Influencing Dissolution of Cocrystals
Dissolution of cocrystals can be influenced by several factors such as aqueous solubility of coformers (as discussed in Section 10.2.2), intermolecular interactions in cocrystals (as discussed in Section 10.2.2), crystal habit of cocrystals (as discussed below) and pH of the dissolution medium (as discussed below).

(a) Crystal Habit
Crystal habit of cocrystals can also influence its dissolution in dissolution medium. A drug can cocrystallize with a coformer molecule in various sizes and shapes depending on different crystallization conditions. A crystallization event can change crystal properties such as habit, polymorphism and size [324]. The term 'crystal habit' is used to describe the general shape of a crystal. Modification of a drug crystal's habit during crystallization can alter its dissolution behavior due to a change in the nature of crystal faces exposed to the dissolution medium. However, studies related to the modification of cocrystal habits and understanding the effect of cocrystal habits on its dissolution properties are limited [325]. Sulfadimidine (SDM)-4-aminosalicylic acid (4-ASA) cocrystal (in 1:1 molar ratio) [polymorph I and polymorph II] were cocrystallized in four different crystal habits. These habits were obtained by solvent evaporation using ethanol (habit I) and acetone (habit II), solvent evaporation followed by grinding (habit III) and spray drying (habit IV) [See Table  12] [325]. It was observed that cocrystals prepared by milling showed highest dissolution than the cocrystal powders prepared via solvent evaporation and spray drying (as shown in Table 12). Milling produced cocrystals of very fine (smaller) particle size thereby increasing its surface area, increasing its flowability and enhancing its dissolution. On the other hand, powders prepared by solvent evaporation exhibited comparatively lesser dissolution than the former. Despite smaller particle size and higher surface area, cocrystals prepared by spray drying process showed lower dissolution which was attributed to the agglomeration of particles [325]. Therefore, Serrano and coworkers [325] suggested that crystal habits of cocrystals with poor pharmaceutical characteristics can be engineered to alter dissolution, flowability and compaction behavior [325].

(b) pH
Cao et al. [326] showed that the dissolution behavior of cocrystals with ionizable components is highly dependent on the interfacial pH at the dissolving solid-liquid interface [326]. When a cocrystal containing non-ionizable drug and ionizable coformer is dissolved in a dissolution medium, the dissociation reaction can modify pH at the solid-liquid interface. Authors studied the dissolution behavior of carbamazepine-saccharin and carbamazepine-salicylic acid cocrystals [326] in a dissolution medium (dissolution medium was prepared by dissolving Sodium Lauryl Sulphate (SLS) in water) with varying pH such as 1.27, 2.16, 3.02, 4.03, 5.97 and 7.66. The dissolution rate of the cocrystals was found to increase with increase in pH [326]. The mechanism of cocrystal dissolution was explained by the process of interfacial mass transport. Interfacial equilibrium model and surface saturation model were used to analyze the mass transport process. The dissolution rates of these cocrystals were found to be greatly influenced by interfacial pH and the bulk pH was not found to adequately explain the dissolution behavior [326].

Cocrystals with Low Dissolution Rates
While cocrystallization can be used to enhance solubility of drugs, it can also be used to reduce the aqueous solubility of drugs [327]. Though enhanced aqueous solubility of poorly water-soluble drugs by formulating as cocrystals has been well-documented in the literature, the number of reports on cocrystals with lower dissolution rates than the raw drug is limited. The following are the case studies of such pharmaceutical cocrystals with lower dissolution rates reported in the literature:

(a) Sulfacetamide Cocrystals
Sulfacetamide is a topical antibiotic used to treat conjunctivitis. Despite its therapeutic efficacy, its pharmaceutical use is limited because of physiological constraints such as tear flow, reflex blinking and drug loss. This drawback can possibly be eradicated by frequent dosing. However, frequent dosing in turn can lead to excess drug loading in patients. Therefore, reducing the aqueous solubility of sulfacetamide can help in getting rid of excess drug loading in patients without reducing the frequency of drug administration. Therefore, Nangia and coworkers [327] adopted cocrystallization approach to bring down the aqueous solubility of sulfacetamide [327]. Sulfacetamide-caffeine (1:1) and sulfacetamide-isonicotinamide (1:1) cocrystals exhibited 0.68 and 0.64 times lower dissolution rates than raw sulfacetamide whereas sulfacetamide-theophylline (1:1) cocrystal exhibited dissolution rate equivalent to that of raw sulfacetamide [327]. It was proposed that modification of intermolecular interactions (hydrogen bonding) existing in the parent drug molecule by means of cocrystallization resulted in better crystal packing, stable crystal lattice and increase in density of the cocrystals. These modifications in the crystal lattice were responsible for lower dissolution of the synthesized sulfacetamide cocrystals, which in turn can increase the residence time of sulfacetamide drug at the site of action [327].

(b) Fluoxetine HCl-Benzoic Acid (1:1) Cocrystal
Fluoxetine HCl is an antidepressant drug available in the market with the trade name, Prozac. During a study carried out by Childs et al. [57] to prepare cocrystals of amine hydrochlorides with organic acids by means of crystal engineering approach, it was found that Fluoxetine HCl-Benzoic acid (1:1) cocrystal exhibited lower powder and intrinsic dissolution rates than that of commercial Fluoxetine HCl [57]. However, fluoxetine HCl cocrystals prepared with other acids such as fumaric acid and succinic acid showed enhanced powder dissolution rates. Thus, it is evident that one can increase or decrease the dissolution rates of drugs by appropriate selection of coformers for cocrystallization. The reason behind the influence of the structural and thermodynamic parameters on the dissolution rates of these cocrystals remain unclear [57]. However, it is to be noticed that the dissolution rates obtained for the cocrystals were observed to be in good correlation with the aqueous solubility of coformers [57].

(c) Curcumin-phloroglucinol (1:1) Cocrystal
Curcumin is a natural phenolic ingredient with anticancer activity. Its therapeutic use is limited due to its poor aqueous solubility, and hence its poor bioavailability. Efforts have been made by several researchers to enhance its aqueous solubility by forming cocrystals [24,72,328]. However, cocrystallization with phloroglucinol formed curcumin-phloroglucinol (1:1) cocrystal with lower dissolution rates [73] though phloroglucinol has higher aqueous solubility.

(d) Lamotrigine-phenobarbital (1:1) Cocrystal
Recently, Kaur et al. [329] reported that lamotrigine-phenobarbital cocrystal (of 1:1 stoichiometric ratio) exhibited lower dissolution rate and poor aqueous solubility than raw phenobarbital and lamotrigine. Pure lamotrigine crystal is stabilized by N-H. . . N amide/pyridine homodimer intermolecular interactions whereas the crystal structure of pure phenobarbital is stabilized by N-H. . . O amine/carbonyl homodimer intermolecular interactions (as shown in Figure 39). In case of lamotrigine-phenobarbital (1:1) cocrystal, these lamotrigine and phenobarbital homodimers are held together by stronger intermolecular interactions (via formation of one N-H. . . N bond and two N-H. . . O bonds between the molecules of lamotrigine and phenobarbital) forming heterodimers (as shown in Figure 39). It was proposed by the authors that the stronger heterodimer interactions in cocrystal led to poor dissolution (as shown in Figure 40) [329]. Fluoxetine HCl is an antidepressant drug available in the market with the trade name, Prozac. During a study carried out by Childs et al. [57] to prepare cocrystals of amine hydrochlorides with organic acids by means of crystal engineering approach, it was found that Fluoxetine HCl-Benzoic acid (1:1) cocrystal exhibited lower powder and intrinsic dissolution rates than that of commercial Fluoxetine HCl [57]. However, fluoxetine HCl cocrystals prepared with other acids such as fumaric acid and succinic acid showed enhanced powder dissolution rates. Thus, it is evident that one can increase or decrease the dissolution rates of drugs by appropriate selection of coformers for cocrystallization. The reason behind the influence of the structural and thermodynamic parameters on the dissolution rates of these cocrystals remain unclear [57]. However, it is to be noticed that the dissolution rates obtained for the cocrystals were observed to be in good correlation with the aqueous solubility of coformers [57].

(c) Curcumin-phloroglucinol (1:1) Cocrystal
Curcumin is a natural phenolic ingredient with anticancer activity. Its therapeutic use is limited due to its poor aqueous solubility, and hence its poor bioavailability. Efforts have been made by several researchers to enhance its aqueous solubility by forming cocrystals [24,72,328]. However, cocrystallization with phloroglucinol formed curcumin-phloroglucinol (1:1) cocrystal with lower dissolution rates [73] though phloroglucinol has higher aqueous solubility.

(d) Lamotrigine-phenobarbital (1:1) Cocrystal
Recently, Kaur et al. [329] reported that lamotrigine-phenobarbital cocrystal (of 1:1 stoichiometric ratio) exhibited lower dissolution rate and poor aqueous solubility than raw phenobarbital and lamotrigine. Pure lamotrigine crystal is stabilized by N-H…N amide/pyridine homodimer intermolecular interactions whereas the crystal structure of pure phenobarbital is stabilized by N-H…O amine/carbonyl homodimer intermolecular interactions (as shown in Figure  39). In case of lamotrigine-phenobarbital (1:1) cocrystal, these lamotrigine and phenobarbital homodimers are held together by stronger intermolecular interactions (via formation of one N-H…N bond and two N-H…O bonds between the molecules of lamotrigine and phenobarbital) forming heterodimers (as shown in Figure 39). It was proposed by the authors that the stronger heterodimer interactions in cocrystal led to poor dissolution (as shown in Figure 40) [329].

In Vitro and In Vivo Studies
Manypharmaceutical cocrystals have been synthesized by several researchers to enhance the aqueous solubility of poorly water-soluble drugs (especially the BCS Class II and class IV drugs). However, the number of cocrystals being used as a regular drug with the approval of FDA is very low. Also, the bioavailability studies for the synthesized cocrystals have been rare. At present, Entresto (used in the treatment of chronic heart failure), Lexapro (used as an antidepressant) and Depakote (used in the treatment of seizure disorders and manic depression) are the three pharmaceutical cocrystals being approved by FDA for clinical use. Table 13 presents a summary of reports available in literature on the bioavailability assessment of pharmaceutical cocrystals. In these reports (listed in Table 13), the pharmacokinetic behavior and therapeutic effect of different solid forms have been studied. It is evident from this table that the pharmaceutical cocrystals help in fine-tuning the dissolution behavior and aqueous solubility of API molecules and in turn enhances their bioavailability. Therefore, there is a need to conduct bioavailability studies for the new solid phases (cocrystals/eutectics/coamorphous solids) prepared by cocrystallization and evaluate their efficacy in terms of their aqueous solubility under biological pH conditions, in vivo permeability and bioavailability. This would pave the way for these cocrystals to enter the next step of clinical trials for necessary approval by FDA.

In Vitro and In Vivo Studies
Manypharmaceutical cocrystals have been synthesized by several researchers to enhance the aqueous solubility of poorly water-soluble drugs (especially the BCS Class II and class IV drugs). However, the number of cocrystals being used as a regular drug with the approval of FDA is very low. Also, the bioavailability studies for the synthesized cocrystals have been rare. At present, Entresto (used in the treatment of chronic heart failure), Lexapro (used as an antidepressant) and Depakote (used in the treatment of seizure disorders and manic depression) are the three pharmaceutical cocrystals being approved by FDA for clinical use. Table 13 presents a summary of reports available in literature on the bioavailability assessment of pharmaceutical cocrystals. In these reports (listed in Table 13), the pharmacokinetic behavior and therapeutic effect of different solid forms have been studied. It is evident from this table that the pharmaceutical cocrystals help in fine-tuning the dissolution behavior and aqueous solubility of API molecules and in turn enhances their bioavailability. Therefore, there is a need to conduct bioavailability studies for the new solid phases (cocrystals/eutectics/coamorphous solids) prepared by cocrystallization and evaluate their efficacy in terms of their aqueous solubility under biological pH conditions, in vivo permeability and bioavailability. This would pave the way for these cocrystals to enter the next step of clinical trials for necessary approval by FDA. The coamorphous phase exhibited greater dissolution and pharmacokinetic profile than raw curcumin. It also showed higher bioavailability and therapeutic effect than raw curcumin. The coamorphous solid was non-toxic even at a dose of 10 times higher dose at 2000 mg/kg in xenograft model [89] Hesperetin

Challenges and Future Perspectives
Cocrystallization of poorly water-soluble drugs is one of the novel ways to improve their aqueous solubility. A lot of research efforts are now focused on synthesizing cocrystals of poorly water-soluble drugs with appropriate coformers. The physicochemical properties of cocrystals such as melting point, aqueous solubility and hence their bioavailability depends upon the type of coformer used. Also, as mentioned earlier, the cocrystallization attempts do not always lead to a successful formation of cocrystals. At times, a eutectic or even a coamorphous solid is obtained. Therefore, choosing a correct coformer is of utmost importance. However, at present, conformers are either chosen based on empirical understanding or based on cumbersome methods requiring detailed analysis and calculations. Therefore, development of a new and fast coformer screening tool is necessary to screen coformers suitable for cocrystallization. Furthermore, efforts are also needed to develop a generalized understanding of intermolecular interactions that influence the cocrystallization outcome by employing supramolecular chemistry and crystal engineering principles. While a rationale design of a cocrystal can lead to a successful outcome at the end of cocrystallization, it is equally important to develop solvent-free cocrystal production methods such as melt crystallization. At present researchers are trying to develop continuous techniques for cocrystal synthesis. Further efforts in this direction will enable a large-scale production of cocrystals. Additionally, future research also needs to focus on stability of cocrystals. At present, very little information is available on the aspects related to the cocrystal stability such as conversion of cocrystals to drug polymorphs or degradation of cocrystals upon exposure to aqueous environment, pH or temperature, etc. It is also evident from the literature reports on cocrystal/eutectic/coamorphous/solid solution formulations reviewed so far that most of the solids prepared by cocrystallization have not been studied for their bioavailability, pharmacokinetics and pharmacodynamics. Mere production of cocrystals without investigating its efficacy, pharmacokinetic or pharmacodynamic behavior in biological system would end up in vain. Moreover, preclinical trials followed by clinical trials have needed to develop cocrystals into marketed product.

Conflicts of Interest:
The authors declare no conflict of interest.