Characterisation and processing of amorphous binary mixtures with low glass transition temperature

The number of drug substances in formulation development in the pharmaceutical industry is increasing. Some of these are amorphous drugs and have glass transition below ambient temperature, and thus they are usually difficult to formulate and handle. One reason for this is the reduced viscosity, related to the stickiness of the drug, that makes them complicated to handle in unit operations. Thus, the aim in this thesis was to develop a new processing method for a sticky amorphous model material. Furthermore, model materials were characterised before and after formulation, using several characterisation methods, to understand more precisely the prerequisites for physical stability of amorphous state against crystallisation. The model materials used were monoclinic paracetamol and citric acid anhydrate. Amorphous materials were prepared by melt quenching or by ethanol evaporation methods. The melt blends were found to have slightly higher viscosity than the ethanol evaporated materials. However, melt produced materials crystallised more easily upon consecutive shearing than ethanol evaporated materials. The only material that did not crystallise during shearing was a 50/50 (w/w, %) blend regardless of the preparation method and it was physically stable at least two years in dry conditions. Shearing at varying temperatures was established to measure the physical stability of amorphous materials in processing and storage conditions. The actual physical stability of the blends was better than the pure amorphous materials at ambient temperature. Molecular mobility was not related to the physical stability of the amorphous blends, observed as crystallisation. Molecular mobility of the 50/50 blend derived from a spectral linewidth as a function of temperature using solid state NMR correlated better with the molecular mobility derived from a rheometer than that of differential scanning calorimetry data. Based on the results obtained, the effect of molecular interactions, thermodynamic driving force and miscibility of the blends are discussed as the key factors to stabilise the blends. The stickiness was found to be affected glass transition and viscosity. Ultrasound extrusion and cutting were successfully tested to increase the processability of sticky material. Furthermore, it was found to be possible to process the physically stable 50/50 blend in a supercooled liquid state instead of a glassy state. The method was not found to accelerate the crystallisation. This may open up new possibilities to process amorphous materials that are otherwise impossible to manufacture into solid dosage forms.

Stabilization of Phospholipid Coatings in Capillary Electrophoresis

The development of a simple method of coating a semi-permanent phospholipid layer onto a capillary for electrochromatography use was the focus of this study. The work involved finding good coating conditions, stabilizing the phospholipid coating, and examining the effect of adding divalent cations, cetyltrimethylammonium bromide, and polyethylene glycol (PEG)-lipids on the stability of the coating. Since a further purpose was to move toward more biological membrane coatings, the capillaries were also coated with cholesterol-containing liposomes and liposomes of red blood cell ghost lipids. Liposomes were prepared by extrusion, and large unilamellar vesicles with a diameter of about 100 nm were obtained. Zwitterionic phosphatidylcholine (PC) was used as a basic component, mainly 1-palmitoyl-2-oleyl-sn-glycero-3-phosphocholine (POPC) but also eggPC and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). Different amounts of sphingomyelin, bovine brain phosphatidylserine, and cholesterol were added to the PC. The stability of the coating in 40 mM N-(2-hydroxyethyl)piperazine-N’-(2-ethanesulfonic acid) (HEPES) solution at pH 7.4 was studied by measuring the electroosmotic flow and by separating neutral steroids, basic proteins, and low-molar-mass drugs. The presence of PC in the coating solution was found to be essential to achieving a coating. The stability of the coating was improved by the addition of negative phosphatidylserine, cholesterol, divalent cations, or PEGylated lipids, and by working in the gel-state region of the phospholipid. Study of the effect on the PC coating of divalent metal ions calcium, magnesium, and zinc showed a molar ratio of 1:3 PC/Ca2+ or PC/Mg2+ to give increased rigidity to the membrane and the best coating stability. The PEGylated lipids used in the study were sterically stabilized commercial lipids with covalently attached PEG chains. The vesicle size generally decreased when PEGylated lipids of higher molar mass were present in the vesicle. The predominance of discoidal micelles over liposomes increased PEG chain length and the average size of the vesicles thus decreased. In the capillary electrophoresis (CE) measurements a highly stable electroosmotic flow was achieved with 20% PEGylated lipid in the POPC coating dispersion, the best results being obtained for disteroyl PEG (3000) conjugates. The results suggest that smaller particles (discoidal micelles) result in tighter packing and better shielding of silanol groups on the silica wall. The effect of temperature on the coating stability was investigated by using DPPC liposomes at temperatures above (45 C) and below (25 C) the main phase transition temperature. Better results were obtained with DPPC in the more rigid gel state than in the fluid state: the electroosmotic flow was heavily suppressed and the PC coating was stabilized. Also dispersions of DPPC with 0−30 mol% of cholesterol and sphingomyelin in different ratios, which more closely resemble natural membranes, resulted in stable coatings. Finally, the CE measurements revealed that a stable coating is formed when capillaries are coated with liposomes of red blood cell ghost lipids.

On the functional ordering of biomembranes : Implications for drug binding

Cells are packed with membrane structures, defining the inside and outside, and the different subcellular compartments. These membranes consisting mainly of phospholipids have a variety of functions in addition to providing a permeability barrier for various compounds. These functions involve cellular signaling, where lipids can act as second messengers, or direct regulation of membrane associating proteins. The first part of this study focuses on relating some of the physicochemical properties of membrane lipids to the association of drug compounds to membranes. A fluorescence based method is described allowing for determination of the membrane association of drugs. This method was subsequently applied to a novel drug, siramesine, previously shown to have anti-cancer activity. Siramesine was found to associate with anionic lipids. Especially interesting is its strong affinity for a second messenger lipid phosphatidic acid. This is the first example of a small molecule drug compound specifically interacting with a cellular lipid. Phosphatidic acid in cells is required for the activation of many signaling pathways mediating growth and proliferation. This provides an intriguing possibility for a simple molecular mechanism of the observed anti-cancer activity of siramesine. In the second part the thermal behavior and self assembly of charged and uncharged membrane assemblies was studied. Strong inter-lamellar co-operativity was observed for multilamellar DPPC vesicles using fluorescence techniques together with calorimetry. The commonly used membrane models, large unilamellar vesicles (LUV) and multilamellar vesicles (MLV) were found to possess different biophysical properties as interlamellar interactions of MLVs drive segregation of a pyrene labeled lipid analogue into clusters. The effect of a counter-ion lattice on the self assembly of a cationic gemini surfactant was studied. The presence of NaCl strongly influenced the thermal phase behavior of M-1 vesicles, causing formation of giant vesicles upon exceeding a phase transition temperature, followed by a subsequent transition into a more homogenous dispersion. Understanding the underlying biophysical aspects of cellular membranes is of fundamental importance as the complex picture of the structure and function of cells is evolving. Many of the cellular reactions take place on membranes and membranes are known to regulate the activity of many peripheral and intergral membrane associating proteins. From the point of view of drug design and gene technology, membranes can provide an interesting target for future development of drugs, but also a vehicle sensitive for environmental changes allowing for encapsulating drugs and targeting them to the desired site of action.

Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Dimensional Reduction Near the Deconfinement Transition

When ordinary nuclear matter is heated to a high temperature of ~ 10^12 K, it undergoes a deconfinement transition to a new phase, strongly interacting quark-gluon plasma. While the color charged fundamental constituents of the nuclei, the quarks and gluons, are at low temperatures permanently confined inside color neutral hadrons, in the plasma the color degrees of freedom become dominant over nuclear, rather than merely nucleonic, volumes. Quantum Chromodynamics (QCD) is the accepted theory of the strong interactions, and confines quarks and gluons inside hadrons. The theory was formulated in early seventies, but deriving first principles predictions from it still remains a challenge, and novel methods of studying it are needed. One such method is dimensional reduction, in which the high temperature dynamics of static observables of the full four-dimensional theory are described using a simpler three-dimensional effective theory, having only the static modes of the various fields as its degrees of freedom. A perturbatively constructed effective theory is known to provide a good description of the plasma at high temperatures, where asymptotic freedom makes the gauge coupling small. In addition to this, numerical lattice simulations have, however, shown that the perturbatively constructed theory gives a surprisingly good description of the plasma all the way down to temperatures a few times the transition temperature. Near the critical temperature, the effective theory, however, ceases to give a valid description of the physics, since it fails to respect the approximate center symmetry of the full theory. The symmetry plays a key role in the dynamics near the phase transition, and thus one expects that the regime of validity of the dimensionally reduced theories can be significantly extended towards the deconfinement transition by incorporating the center symmetry in them. In the introductory part of the thesis, the status of dimensionally reduced effective theories of high temperature QCD is reviewed, placing emphasis on the phase structure of the theories. In the first research paper included in the thesis, the non-perturbative input required in computing the g^6 term in the weak coupling expansion of the pressure of QCD is computed in the effective theory framework at an arbitrary number of colors. The two last papers on the other hand focus on the construction of the center-symmetric effective theories, and subsequently the first non-perturbative studies of these theories are presented. Non-perturbative lattice simulations of a center-symmetric effective theory for SU(2) Yang-Mills theory show --- in sharp contrast to the perturbative setup --- that the effective theory accommodates a phase transition in the correct universality class of the full theory. This transition is seen to take place at a value of the effective theory coupling constant that is consistent with the full theory coupling at the critical temperature.

Monte Carlo simulations of molecular clusters in nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

Polymeereillä stabiloitu amorfinen formulaatio

Nowadays growing number of new active pharmaceutical ingredients (API) have large molecular weight and are hydrophobic. The energy of their crystal lattice is bigger and polarity has decreased. This leads to weakened solubility and dissolution rate of the drug. These properties can be enhanced for example by amorphization. Amorphous form has the best dissolution rate in the solid state. In the amorphous form drug molecules are randomly arranged, so the energy required to dissolve molecules is lower compared to the crystalline counterpart. The disadvantage of amorphous form is that it is unstable. Amorphous form tends to crystallize. Stability of amorphous form can be enhanced by adding an adjuvant to drug product. Adjuvant is usually a polymer. Polymers prevent crystallization both by forming bonds with API molecules and by steric hindrance. The key thing in stabilizing amorphous form is good miscibility between API and polymer. They have to be mixed in a molecular level so that the polymer is able to prevent crystallization. The aim of this work was to study miscibility of drug and polymer and stability of their dispersion with different analytical methods. Amorphous dispersions were made by rotary evaporator and freeze dryer. Amorphicity was confirmed with X-ray powder diffraction (XRPD) right after preparation. Itraconazole and theophylline were the chosen molecules to be stabilized. Itraconazole was expected to be easier and theophylline more difficult to stabilize. Itraconazole was stabilized with HPMC and theophylline was stabilized with PVP. Miscibility was studied with XRPD and differential scanning calorimetry (DSC). In addition it was studied with polarized light microscope if miscibility was possible to see visually. Dispersions were kept in stressed conditions and the crystallization was analyzed with XRPD. Stability was also examined with isothermal microcalorimetry (IMC). The dispersion of itraconazole and theophylline 40/60 (w/w) was completely miscible. It was proved by linear combination of XRPD results and single glass transition temperature in DSC. Homogenic well mixed film was observed with light microscope. Phase separation was observed with other compositions. Dispersions of theophylline and PVP mixed only partly. Stability of itraconazole dispersions were better than theophylline dispersions which were mixed poorer. So miscibility was important thing considering stability. The results from isothermal microcalorimetry were similar to results from conventional stability studies. Complementary analytical methods should be used when studying miscibility so that the results are more reliable. Light microscope is one method in addition to mostly used XRPD and DSC. Analyzing light microscope photos is quite subjective but it gives an idea of miscibility. Isothermal microcalorimetry can be one option for conventional stability studies. If right conditions can be made where the crystallization is not too fast, it may be possible to predict stability with isothermal microcalorimetry.