Widom Line in Supercritical Water in Terms of Changes in Local Structure: Theoretical Perspective.

Performing molecular dynamics simulations with the TIP4P/2005 water model along 9 isobars (from 175 to 375 bar) in the temperature range between 300 and 1100 K, we have found that the loci of the extrema in the rate of change of specific structural properties can be used to define purely structure-based Widom lines. We have examined several parameters that describe the local structure of water, such as the tetrahedral arrangement, nearest neighbor distance, local density around the molecules, and the size of the largest dense domain. The last two parameters were determined using the Voronoi polyhedral and density-based spatial clustering of applications with noise methods, respectively. By analyzing the moments of the associated distributions, we show that along a given isobar, the temperature at which we observe a maximum in the fluctuation, the rate of change of the average values, or in the skewness values unambiguously determines the Widom line that is in agreement with the experimentally detected, thermodynamic response function-based ones.

Impact of the polymer dispersity on the properties of curcumin/polyvinylpyrrolidone amorphous solid dispersions.

Amorphous solid dispersions (ASD) are known to enhance the absorption of poorly water-soluble drugs. In this work we synthesise well-defined Polyvinylpyrrolidone (PVP) to establish the impact of dispersity and chain-end functionality on the physical properties of Curcumin (CUR)/PVP ASD. Thermodynamic characterisation of synthesised PVP emphasises a strong effect of the dispersity on the glass transition temperature (Tg), 50 °C higher for synthesised PVP than for commercial PVP K12 of same molar mass. This increase of Tg affects the thermodynamic properties of CUR/PVP ASD successfully formulated up to 70 wt% of CUR by milling or solvent evaporation. The evolution of both the Tg and CUR solubility values versus CUR content points out the development of fairly strong CUR-PVP interactions that strengthen the antiplasticising effect of PVP on the Tg of ASD. However, for ASD formulated with commercial PVP this effect is counterbalanced at low CUR content by a plasticising effect due to the shortest PVP chains. Moreover, the overlay of the phase and state diagrams highlights the strong impact of the polymer dispersity on the stability of CUR/PVP ASD. ASD formulated with low dispersity PVP are stable on larger temperature and concentration ranges than those formulated with PVP K12.

Using synchrotron high-resolution powder X-ray diffraction for the structure determination of a new cocrystal formed by two active principle ingredients.

The crystal structure of a new 1:1 cocrystal of carbamazepine and S-naproxen (C15H12N2O·C14H14O3) was solved from powder X-ray diffraction (PXRD). The PXRD pattern was measured at the high-resolution beamline CRISTAL at synchrotron SOLEIL (France). The structure was solved using Monte Carlo simulated annealing, then refined with Rietveld refinement. The positions of the H atoms were obtained from density functional theory (DFT) ground-state calculations. The symmetry is orthorhombic with the space group P212121 (No. 19) and the following lattice parameters: a = 33.5486 (9), b = 26.4223 (6), c = 5.3651 (10) Šand V = 4755.83 (19) Å3.

Evidence of Strong Guest-Host Interactions in Simvastatin Loaded in Mesoporous Silica MCM-41.

A rational design of drug delivery systems requires in-depth knowledge not only of the drug itself, in terms of physical state and molecular mobility, but also of how it is distributed among a carrier and its interactions with the host matrix. In this context, this work reports the behavior of simvastatin (SIM) loaded in mesoporous silica MCM-41 matrix (average pore diameter ~3.5 nm) accessed by a set of experimental techniques, evidencing that it exists in an amorphous state (X-ray diffraction, ssNMR, ATR-FTIR, and DSC). The most significant fraction of SIM molecules corresponds to a high thermal resistant population, as shown by thermogravimetry, and which interacts strongly with the MCM silanol groups, as revealed by ATR-FTIR analysis. These findings are supported by Molecular Dynamics (MD) simulations predicting that SIM molecules anchor to the inner pore wall through multiple hydrogen bonds. This anchored molecular fraction lacks a calorimetric and dielectric signature corresponding to a dynamically rigid population. Furthermore, differential scanning calorimetry showed a weak glass transition that is shifted to lower temperatures compared to bulk amorphous SIM. This accelerated molecular population is coherent with an in-pore fraction of molecules distinct from bulklike SIM, as highlighted by MD simulations. MCM-41 loading proved to be a suitable strategy for a long-term stabilization (at least three years) of simvastatin in the amorphous form, whose unanchored population releases at a much higher rate compared to the crystalline drug dissolution. Oppositely, the surface-attached molecules are kept entrapped inside pores even after long-term release assays.

The role of the diffusion in the predictions of the classical nucleation theory for quasi-real systems differ in dipole moment value.

In this paper, we examine the crystallization tendency for two quasi-real systems, which differ exclusively in the dipole moment's value. The main advantage of the studied system is the fact that despite that their structures are entirely identical, they exhibit different physical properties. Hence, the results obtained for one of the proposed model systems cannot be scaled to reproduce the results for another corresponding system, as it can be done for simple model systems, where structural differences are modeled by the different parameters of the intermolecular interactions. Our results show that both examined systems exhibit similar stability behavior below the melting temperature. This finding is contrary to the predictions of the classical nucleation theory, which suggests a significantly higher crystallization tendency for a more polar system. Our studies indicate that the noted discrepancies are caused by the kinetic aspect of the classical nucleation theory, which overestimates the role of diffusion in the nucleation process.

Pressure Dependence of the Crystallization Rate for the S-Enantiomer and a Racemic Mixture of Ibuprofen.

This paper examines the pressure effect on the crystallization rate of the pharmaceutically active enantiomerically pure S-enantiomer and the racemic mixture of the well-known drug ibuprofen. Performed experimental studies revealed that at ambient pressure S-ibuprofen crystallizes faster than the racemic mixture. When the pressure increases, the crystallization rate slows down for both systems, but interestingly it is more apparent in the case of the S-enantiomer. It is found that this experimentally observed trend can be understood based on the predictions of the classical nucleation theory. We suggest that the solid-liquid interfacial free energy is the main reason for the observed variations in S- and RS-ibuprofen's stability behaviors. Employing a special method of computational studies, i.e., the capillary fluctuation method, we show that the increase in pressure affects the solid-liquid interfacial free energy for S- and RS-ibuprofen in an entirely different way. Importantly, the detected differences correspond to the experimentally observed variations in the overall crystallization rates.

Structure determination of riboflavin by synchrotron high-resolution powder X-ray diffraction.

The crystal structure of the stable form of vitamin B2 or riboflavin (C17H20N4O6) was solved using high-resolution powder X-ray diffraction (PXRD). The high-resolution PXRD pattern of riboflavin was recorded at room temperature at the European Synchrotron Radiation Facility (Grenoble, France). The starting structural model was generated using a Monte Carlo simulated annealing method. The final structure was obtained through Rietveld refinement. The positions of the H atoms belonging to hydroxy groups were estimated from computational energy minimizations. The symmetry is orthorhombic with the space group P212121 and the following lattice parameters: a = 20.01308, b = 15.07337 and c = 5.31565 Å.

Impact of chirality on the amorphous state of conglomerate forming systems: a case study of N-acetyl-α-methylbenzylamine.

The present work aims at addressing the issue of molecular handedness in glassy and liquid states and its impact on heterogeneous equilibrium. For this purpose, we evaluated the glass forming ability (GFA), crystallization propensity, molecular mobility and hydrogen bonding structure of a chiral conglomerate forming system, N-acetyl-α-methylbenzylamine (Nac-MBA), at various enantiomeric excesses (ees) using experimental and computational techniques. We revealed that the rich relaxational landscape (Debye (D), α, ßJG and ϒ) and the temperature dependence of the time scale of each process were insensitive to chirality. The most remarkable impact of chirality was expressed on the GFA and the recrystallization of heterochiral arrangements. In fact the GFA increases with decreasing ee, while the crystallization propensity increases with increasing ee. The counter enantiomer acted as a disruptor of crystallization and favored the glass formation upon cooling. The molecular dynamics simulation (MDS) results on the architecture of chiral sequences showed that homochiral sequences were more favorable when compared to heterochiral ones in the liquid state. However, this predisposition to form homochiral sequences in the liquid state was not the precursor of the future crystalline structure, since the liquid or the glassy system recrystallizes as heterochiral sequences. As per our understanding the crystallization was mostly controlled by the mean free migration path of an enantiomer to build homochiral or heterochiral sequences. In the present case, it seems that the mean free migration path achieved by an enantiomer for heterochiral sequences is shorter compared to homochiral arrangements in such a way that the crystallization of the metastable racemic compound is kinetically more favorable.

Impact of Low Concentration of Strongly Hydrogen-Bonded Water Molecules on the Dynamics of Amorphous Terfenadine: Insights from Molecular Dynamics Simulations and Dielectric Relaxation Spectroscopy.

The impact of low water concentration of strongly hydrogen-bonded water molecules on the dynamical properties of amorphous terfenadine (TFD) is investigated through complementary molecular dynamics (MD) simulations and dielectric relaxation spectroscopy (DRS) experiments. In this article, we especially highlight the important role played by some residual water molecules in the concentration of 1-2% (w/w) trapped in the TFD glassy matrix, which are particularly difficult to remove experimentally without a specific heating/drying process. From MD computations and analyses of the hydrogen bonding (HB) interactions, different categories of water molecules are revealed and particularly the presence of strongly HB water molecules. These latter localize themselves in small pockets in empty spaces existing in between the TFD molecules due to the poor packing of the glassy state and preferentially interact with the polar groups close to the flexible central part of the TFD molecules. We present a simple model which rationalizes at the molecular scale the effect of these strongly HB water molecules on dynamics and how they give rise to a supplementary relaxation process (namely process S) which is detected for the first time in the glassy state of TFD annealed at room temperature while this process is completely absent in a non-annealed glass. It also explains how this supplementary relaxation is coupled with the intramolecular motion (namely process γ) of the very flexible central part of the TFD molecule. The present findings help to understand more generally the microscopic origin of the secondary relaxations often detected by DRS in the glassy states of molecular compounds for which the exact nature is still debated.

Structure determination of a new cocrystal of carbamazepine and DL-tartaric acid by synchrotron powder X-ray diffraction.

The crystal structure of a new cocrystal of carbamazepine (systematic name: 5H-dibenzo[b,f]azepine-5-carboxamide, C15H12N2O) and DL-tartaric acid (C4H6O6), obtained by liquid-assisted grinding, was solved by powder X-ray diffraction (PXRD). The high-resolution PXRD pattern of this new phase was recorded at room temperature thanks to synchrotron experiments at the European Synchrotron Radiation Facility (Grenoble, France). The starting structural model was generated by a Monte-Carlo simulated annealing method. The final structure was obtained through Rietveld refinement and an energy minimization simulation was used to estimate the H-atom positions. The stability of the proposed structure as a function of temperature was also assessed from molecular dynamics simulations. The symmetry is monoclinic (space group P21/c) and contains eight molecules per unit cell, namely, four DL-tartaric acid and four carbamazepine molecules.

Insight From Molecular Dynamics Simulations on the Crystallization Tendency of Indomethacin Polymorphs in the Undercooled Liquid State.

The crystallization tendency of 2 crystalline polymorphs of indomethacin (Iα, Iγ) in the undercooled melt has been investigated using molecular dynamics simulations. The main thermodynamical and dynamical physical parameters involved in the nucleation and growth processes have been determined. A careful attention has been given to the crystal-liquid interfacial free energy which remains really challenging to determine from experiments. The present work particularly sheds the light on the importance of the interplay between the solid-liquid interfacial free energy and the driving force for crystallization. The nucleation and the growth rates have been also estimated in the framework of the classical nucleation theory and some growth modes (normal mode, two-dimensional nucleation, and screw dislocation).

Impact of Amorphization Methods on the Physicochemical Properties of Amorphous Lactulose.

The influence of the amorphization technique on the physicochemical properties of amorphous lactulose was investigated. Four different amorphization techniques were used: quenching of the melt, milling, spray-drying, and freeze-drying, and amorphous samples were analyzed by differential scanning calorimetry, NMR spectroscopy, and powder X-ray diffraction analysis. Special attention was paid to the tautomeric composition and to the glass transition of amorphized materials. It was found that the tautomeric composition of the starting physical state (crystal, liquid, or solution) is preserved during the amorphization process and has a strong repercussion on the glass transition of the material. The correlation between these two properties as well as the plasticizing effect of the different tautomers was clarified by molecular dynamics simulations.

Morphological and Structural Properties of Amorphous Lactulose Studied by Scanning Electron Microscopy, Polarized Neutron Scattering, and Molecular Dynamics Simulations.

Morphological and structural properties of amorphous disaccharide lactulose (C12H22O11), obtained by four different amorphization methods (milling, quenching of the melt form, spray-drying, and freeze-drying), are investigated by scanning electron microscopy, polarized neutron scattering, and molecular dynamics simulations. While major differences on the morphology of the different amorphous samples are revealed by scanning electron microscopy images, only subtle structural differences have been found by polarized neutron scattering. Microstructure of the milled sample appears slightly different from the other amorphized materials with the presence of remaining crystalline germs which are not detected by X-ray diffraction. Quantitative phase analysis shows that these remaining crystallites are present in a ratio between 1 and 4%, and their size remains between 20 and 30 nm despite a long milling time of about 8 h. The impact of the change in tautomeric concentrations on the physical properties of lactulose in the amorphous state has been investigated from molecular dynamics simulations. It is suggested that chemical differences between lactulose tautomers could be at the origin of small structural differences detected by polarized neutron scattering.

Molecular Mobility of Terfenadine: Investigation by Dielectric Relaxation Spectroscopy and Molecular Dynamics Simulation.

The molecular mobility of an amorphous active pharmaceutical ingredient, terfenadine, was carefully investigated by dielectric relaxation spectroscopy and molecular dynamics simulation for the first time. Comprehensive characterization on a wide frequency (10-2 to 109 Hz) and temperature (300 K) range highlights the fragile nature of this good glass-former (m = 112) and the relatively large nonexponentiality of the main relaxation (ßKWW = 0.53 ± 0.01). In the glassy state, a particularly broad secondary relaxation of intramolecular origin is evidenced. Terfenadine is a flexible molecule, and from molecular dynamics simulation, a clear link is established between the flexibility of the central part of the molecule (carrying, on the one side, the nitrogen group, and on the other side, the OH group) and the distribution of dipole moments, which explains that broadness. Terfenadine is one of the very few cases for which the molecular mobility of the glass obtained by the quench of the melt or by milling can be compared. From the present study, no major difference in terms of molecular mobility is found between these two glasses. However, terfenadine amorphized by milling (for 1-20 h) clearly shows a lower stability than the quenched liquid as we observed its recrystallization upon heating. Interestingly, it is shown that this recrystallization upon heating is not complete and that the 1-2% of the remaining amorphous phase has an original behavior. Indeed, it exhibits an enhanced main mobility induced by an autoconfinement effect created by the surrounding crystalline phase.

Lactulose: A Model System to Investigate Solid State Amorphization Induced by Milling.

In this article, we show that crystalline lactulose can be amorphized directly in the solid state by mechanical milling. Moreover, compared to similar materials, the amorphization kinetics of lactulose is found to be very rapid and the amorphous state thus obtained appears to be very stable against recrystallization on heating. These features make lactulose a model compound for this type of solid state transformation. The ease of crystalline lactulose to be amorphized on milling is explained by comparing elastic constants of lactulose with those of several other disaccharides. These constants have been determined by molecular dynamics simulations. The article also shows how isothermal dissolution calorimetry can be used effectively for the determination of amorphization kinetics during grinding when the usual characterization techniques (differential scanning calorimetry and powder X-ray diffraction) fail.

Molecular mobility of amorphous N-acetyl-α-methylbenzylamine and Debye relaxation evidenced by dielectric relaxation spectroscopy and molecular dynamics simulations.

The present work focusses on the molecular mobility characterization of amorphous N-acetyl-α-methylbenzylamine (Nac-MBA) by Broadband Dielectric Relaxation Spectroscopy (DRS) coupled with Fast Scanning Calorimetry (FSC) and Molecular Dynamics (MD) simulations covering over 12 decades in the frequency range. This study reveals another example of a secondary amide that shows a very intense Debye-like contribution (almost 90% of the global dielectric intensity) in addition to the structural α-relaxation and secondary Johari-Goldstein ß-relaxation. The D- and α-relaxations are separated by about one decade (in frequency) and their relaxation times follow a near parallel temperature evolution (Vogel-Fulcher-Tammann-Hesse). The micro-structure of Nac-MBA has been investigated from MD simulations. It is shown that the intense Debye-like process emanates from the formation of linear intermolecular H-bonding aggregates (precursors of the crystalline structure) generating super-dipole moments.

Impact of chirality on peculiar ibuprofen molecular dynamics: hydrogen bonding organization and syn vs. anti carboxylic group conformations.

Studies of the impact of chirality on amorphous states are scarce. Here, we present combined dielectric relaxation spectroscopy (DRS) experiments and molecular dynamics (MD) simulation investigations of homochiral and racemic ibuprofen in the liquid, undercooled liquid and glassy states. The influence of chirality is particularly investigated on the syn and anti conformations of the -COOH moiety of the ibuprofen molecule and its link to the peculiar Debye-like dynamical process detected in this compound. Most of the studied properties are found to be nearly identical in the homochiral and racemic systems. But the polarity and intensity of the Debye-like process are clearly found to be more intense in the racemic mixture than in the enantiomerically pure ibuprofen. The difference is explained by the higher population of the anti conformation (with the higher dipole moment) and the lower population of hydrogen bonded cyclic dimers that can be transiently formed in the racemic mixture.

Molecular Packing, Hydrogen Bonding, and Fast Dynamics in Lysozyme/Trehalose/Glycerol and Trehalose/Glycerol Glasses at Low Hydration.

Water and glycerol are well-known to facilitate the structural relaxation of amorphous protein matrices. However, several studies evidenced that they may also limit fast (∼picosecond-nanosecond, ps-ns) and small-amplitude (∼Å) motions of proteins, which govern their stability in freeze-dried sugar mixtures. To determine how they interact with proteins and sugars in glassy matrices and, thereby, modulate their fast dynamics, we performed molecular dynamics (MD) simulations of lysozyme/trehalose/glycerol (LTG) and trehalose/glycerol (TG) mixtures at low glycerol and water concentrations. Upon addition of glycerol and/or water, the glass transition temperature, Tg, of LTG and TG mixtures decreases, the molecular packing of glasses is improved, and the mean-square displacements (MSDs) of lysozyme and trehalose either decrease or increase, depending on the time scale and on the temperature considered. A detailed analysis of the hydrogen bonds (HBs) formed between species reveals that water and glycerol may antiplasticize the fast dynamics of lysozyme and trehalose by increasing the total number and/or the strength of the HBs they form in glassy matrices.

How strongly does trehalose interact with lysozyme in the solid state? Insights from molecular dynamics simulation and inelastic neutron scattering.

Therapeutic proteins are usually conserved in glassy matrixes composed of stabilizing excipients and a small amount of water, which both control their long-term stability, and thus their potential use in medical treatments. To shed some light on the protein-matrix interactions in such systems, we performed molecular dynamics (MD) simulations on matrixes of (i) the model globular protein lysozyme (L), (ii) the well-known bioprotectant trehalose (T), and (iii) the 1:1 (in weight) lysozyme/trehalose mixture (LT), at hydration levels h of 0.0, 0.075, and 0.15 (in g of water/g of protein or sugar). We also supplemented these simulations with complementary inelastic neutron scattering (INS) experiments on the L, T, and LT lyophilized (freeze-dried) samples. The densities and free volume distributions indicate that trehalose improves the molecular packing of the LT glass with respect to the L one. Accordingly, the low-frequency vibrational densities of states (VDOS) and the mean square displacements (MSDs) of lysozyme reveal that it is less flexible-and thus less likely to unfold-in the presence of trehalose. Furthermore, at low contents (h = 0.075), water systematically stiffens the vibrational motions of lysozyme and trehalose, whereas it increases their MSDs on the nanosecond (ns) time scale. This stems from the hydrogen bonds (HBs) that lysozyme and trehalose form with water, which, interestingly, are stronger than the ones they form with each other but which, nonetheless, relax faster on the ns time scale, given the larger mobility of water. Moreover, lysozyme interacts preferentially with water in the hydrated LT mixtures, and trehalose appears to slow down significantly the relaxation of lysozyme-water HBs. Overall, our results suggest that the stabilizing efficiency of trehalose arises from its ability to (i) increase the number of HBs formed by proteins in the dry state and (ii) make the HBs formed by water with proteins stable on long (>ns) time scales.

Structure determination of L-arabinitol by powder X-ray diffraction.

Powder X-ray diffraction patterns of the commercial phase of L-arabinitol were recorded with a laboratory diffractometer. The starting structural model was found by a Monte-Carlo simulated annealing method. The final structure was obtained through Rietveld refinements with soft restraints on the interatomic bond lengths and bond angles. H atoms of hydroxyl groups were localized by minimization of the crystalline energy. The cell is triclinic with the space group P1 and contains two molecules. The crystalline cohesion is achieved by an important network of O-H···O hydrogen bonds.