On the calculation of solubilities via direct coexistence simulations: Investigation of NaCl aqueous solutions and Lennard-Jones binary mixtures.

Direct coexistence molecular dynamics simulations of NaCl solutions and Lennard-Jones binary mixtures were performed to explore the origin of reported discrepancies between solubilities obtained by direct interfacial simulations and values obtained from the chemical potentials of the crystal and solution phases. We find that the key cause of these discrepancies is the use of crystal slabs of insufficient width to eliminate finite-size effects. We observe that for NaCl crystal slabs thicker than 4 nm (in the direction perpendicular to the interface), the same solubility values are obtained from the direct coexistence and chemical potential routes, namely, 3.7 ± 0.2 molal at T = 298.15 K and p = 1 bar for the JC-SPC/E model. Such finite-size effects are absent in the Lennard-Jones system and are likely caused by surface dipoles present in the salt crystals. We confirmed that µs-long molecular dynamics runs are required to obtain reliable solubility values from direct coexistence calculations, provided that the initial solution conditions are near the equilibrium solubility values; even longer runs are needed for equilibration of significantly different concentrations. We do not observe any effects of the exposed crystal face on the solubility values or equilibration times. For both the NaCl and Lennard-Jones systems, the use of a spherical crystallite embedded in the solution leads to significantly higher apparent solubility values relative to the flat-interface direct coexistence calculations and the chemical potential values. Our results have broad implications for the determination of solubilities of molecular models of ionic systems.

Computational studies of pressure, temperature, and surface effects on the structure and thermodynamics of confined water.

The behavior of water confined on nanometer length scales is important in a diverse set of technical and scientific contexts, ranging from the performance of fuel cells and biological molecular machines to the design of self-assembling nanoscale materials. Here, we review recent insights into the structure and thermodynamics of confined water that have been elucidated primarily by computer simulation studies. We emphasize investigations in which interfacial chemistry and molecular topography are varied systematically and in which a wide range of thermodynamic conditions of temperature and pressure are explored. We consider homogeneous interfaces ranging from the simplest hard wall to chemically realistic, but structurally ideal, hydrophobic and hydrophilic surfaces, and the continuous scale of surface polarity is investigated. Features associated with interface heterogeneities arising from chemical patterning or from the natural characteristics of protein surfaces are discussed. Finally, we provide our thoughts on important directions for further studies.

Hydrogen bond strength and network structure effects on hydration of non-polar molecules.

We measure the solvation free energy, Δµ*, for hard spheres and Lennard-Jones particles in a number of artificial liquids made from modified water models. These liquids have reduced hydrogen bond strengths or altered bond angles. By measuring Δµ* for a number of state points at P = 1 bar and different temperatures, we obtain solvation entropies and enthalpies, which are related to the temperature dependence of the solubilities. By resolving the solvation entropy into the sum of the direct solute-solvent interaction and a term depending on the solvent reorganisation enthalpy we show that, although the hydrophobic effect in water at 300 K arises mainly from the small molecular size, its temperature dependence is anomalously low because the reorganisation enthalpy of liquid water is unusually small. We attribute this to the strong tetrahedral network which results from both the molecular geometry and the hydrogen bond strength.

Test of nonequilibrium thermodynamics in glassy systems: the soft-sphere case.

The scaling properties of the soft-sphere potential allow the derivation of an exact expression for the pressure of a frozen liquid, i.e., the pressure corresponding to configurations which are local minima in its multidimensional potential energy landscape. The existence of such a relation offers the unique possibility for testing the recently proposed extension of the liquid free energy to glassy out-of-equilibrium conditions and the associated expression for the temperature of the configurational degrees of freedom. We demonstrate that the nonequilibrium free energy provides an exact description of the soft-sphere pressure in glass states.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling.

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

Supercooled liquids and the glass transition.

Glasses are disordered materials that lack the periodicity of crystals but behave mechanically like solids. The most common way of making a glass is by cooling a viscous liquid fast enough to avoid crystallization. Although this route to the vitreous state-supercooling-has been known for millennia, the molecular processes by which liquids acquire amorphous rigidity upon cooling are not fully understood. Here we discuss current theoretical knowledge of the manner in which intermolecular forces give rise to complex behaviour in supercooled liquids and glasses. An intriguing aspect of this behaviour is the apparent connection between dynamics and thermodynamics. The multidimensional potential energy surface as a function of particle coordinates (the energy landscape) offers a convenient viewpoint for the analysis and interpretation of supercooling and glass-formation phenomena. That much of this analysis is at present largely qualitative reflects the fact that precise computations of how viscous liquids sample their landscape have become possible only recently.

Relationship between structural order and the anomalies of liquid water.

In contrast to crystalline solids--for which a precise framework exists for describing structure--quantifying structural order in liquids and glasses has proved more difficult because even though such systems possess short-range order, they lack long-range crystalline order. Some progress has been made using model systems of hard spheres, but it remains difficult to describe accurately liquids such as water, where directional attractions (hydrogen bonds) combine with short-range repulsions to determine the relative orientation of neighbouring molecules as well as their instantaneous separation. This difficulty is particularly relevant when discussing the anomalous kinetic and thermodynamic properties of water, which have long been interpreted qualitatively in terms of underlying structural causes. Here we attempt to gain a quantitative understanding of these structure-property relationships through the study of translational and orientational order in a models of water. Using molecular dynamics simulations, we identify a structurally anomalous region--bounded by loci of maximum orientational order (at low densities) and minimum translational order (at high densities)--in which order decreases on compression, and where orientational and translational order are strongly coupled. This region encloses the entire range of temperatures and densities for which the anomalous diffusivity and thermal expansion coefficient of water are observed, and enables us to quantify the degree of structural order needed for these anomalies to occur. We also find that these structural, kinetic and thermodynamic anomalies constitute a cascade: they occur consecutively as the degree of order is increased.

Protein purification with vapor-phase carbon dioxide.

Gaseous CO2 was used as an antisolvent to induce the fractional precipitation of alkaline phosphatase, insulin, lysozyme, ribonuclease, trypsin, and their mixtures from dimethylsulfoxide (DMSO). Compressed CO2 was added continuously and isothermally to stationary DMSO solutions (gaseous antisolvent, GAS). Dissolution of CO2 was accompanied by a pronounced, pressure-dependent volumetric expansion of DMSO and a consequent reduction in solvent strength of DMSO towards dissolved proteins. View cell experiments were conducted to determine the pressures at which various proteins precipitate from DMSO. The solubility of each protein in CO2-expanded DMSO was different, illustrating the potential to separate and purify proteins using gaseous antisolvents. Polyacrylamide gel electrophoresis in sodium dodecyl sulfate (SDS-PAGE) was used to quantify the separation of lysozyme from ribonuclease, alkaline phosphatase from insulin, and trypsin from catalase. Lysozyme biological activity assays were also performed to determine the composition of precipitates from DMSO initially containing lysozyme and ribonuclease. SDS-PAGE characterizations suggest that the composition and purity of solid-phase precipitated from a solution containing multiple proteins may be accurately controlled through the antisolvent's pressure. Insulin, lysozyme, ribonuclease, and trypsin precipitates recovered substantial amounts of biological activity upon redissolution in aqueous media. Alkaline phosphatase, however, was irreversibly denaturated. Vapor-phase antisolvents, which are easily separated and recovered from proteins and liquid solvents upon depressurization, appear to be a reliable and effective means of selectively precipitating proteins.

Long-term and high-temperature storage of supercritically-processed microparticulate protein powders.

PURPOSE: The long-term and high-temperature storage of dry, micron-sized particles of lysozyme, trypsin, and insulin was investigated. Subsequent to using supercritical carbon dioxide as an antisolvent to induce their precipitation from a dimethylsulfoxide solution, protein microparticles were stored in sealed containers at -25, -15, 0, 3, 20, 22, and 60 degrees C. The purpose of this study was to investigate the suitability of supercritical antisolvent precipitation as a finishing step in protein processing. METHODS: Karl Fisher titrations were used to determine the residual moisture content of commercial and supercritically-processed protein powders. The secondary structure of the dry protein particles was determined periodically during storage using Raman spectroscopy. The proteins were also redissolved periodically in aqueous buffers and assayed spectrophotometrically for biological activity and by circular dichroism for structural conformation in solution. RESULTS: Amide I band Raman spectra indicate that the secondary structure of the protein particles, while perturbed from that of the solution state, remained constant in time, regardless of the storage temperature. The recoverable biological activity upon reconstitution for the supercritically-processed lysozyme and trypsin microparticles was also preserved and found to be independent of storage temperature. Far UV circular dichroism spectra support the bioactivity assays and further suggest that adverse structural changes, with potential to hinder renaturation upon redissolution, do not take place during storage. CONCLUSIONS: The present study suggests that protein precipitation using supercritical fluids may yield particles suitable for long-term storage at ambient conditions.

Precipitation of proteins in supercritical carbon dioxide.

Supercritical CO2 was used as an antisolvent to form protein particles that exhibited minimal loss of activity upon reconstitution. Organic protein solutions were sprayed under a variety of operating conditions into the supercritical fluid, causing precipitation of dry, microparticulate (1-5 microns) protein powders. Three proteins were studied: trypsin, lysozyme, and insulin. Amide I band Raman spectra were used to estimate the alpha-helix and beta-sheet structural contents of native and precipitate powders of each protein. Analysis of the Raman spectral revealed minimal (lysozyme), intermediate (trypsin), and appreciable (insulin) changes in secondary structure with respect to the commercial starting materials. The perturbations in secondary structure suggest that the most significant event during supercritical fluid-induced precipitation involved the formation of beta-sheet structures with concomitant decreases of alpha-helix. Amide I band Raman and Fourier-transform infrared (FTIR) spectra indicate that higher operating temperatures and pressures lead to more extensive beta-sheet-mediated intermolecular interactions in the precipitates. Raman and FTIR spectra of redissolved precipitates are similar to those of aqueous commercial proteins, indicating that conformational changes were reversible upon reconstitution. These results suggest that protein precipitation in supercritical fluids can be used to form particles suitable for controlled release, direct aerosol delivery to the lungs, and long-term storage at ambient conditions.

Freeze crystallization of imipenem.

This work is a study of the freeze crystallization process developed for the production of a rapidly soluble and stable crystalline form of imipenem, an antibiotic. The objective is to understand the relationship between process conditions and the product crystallinity during the freeze crystallization. Solutions of imipenem and sodium bicarbonate, a drug stabilizer, were crystallized in acetone-water solvent systems at various ratios and temperatures. The degree of crystallinity of the resulting products was measured using an X-ray powder diffraction analysis technique. The acetone-water S-L phase diagram was used to correlate the resulting degrees of crystallinity of imipenem products and equilibrium properties of the system. It was discovered that the conversion of amorphous imipenem to crystalline products was directly related to the percentage of equilibrium liquid during the freeze crystallization. Solubility data indicated that imipenem is virtually insoluble in these equilibrium liquid phases. This suggests that the phase transition from the amorphous state to the crystalline state is mediated by the presence of the liquid phase.

Formation of bioerodible polymeric microspheres and microparticles by rapid expansion of supercritical solutions.

Polyhydroxy acids [poly(L-lactic acid) (L-PLA), poly(D,L-lactic acid) (DL-PLA), and poly-(glycolic acid) (PGA)], biocompatible and bioerodible polymers that are being investigated for controlled delivery of pharmaceuticals and are approved by the Food and Drug Administration for in vivo sutures and bone repair implants, have been dissolved in supercritical CO2 and precipitated by rapid expansion of the resulting supercritical solutions (RESS). The formation of these microparticles and microspheres is a first step toward the goal of producing, in a single processing step, drug-loaded polymeric microspheres for use in controlled release applications. Nucleation of poly(L-lactic acid) from CO2 and CO2-acetone mixtures produced microparticles and microspheres ranging from 4 to 25 microns. Microspheres (2-20 microns) were also obtained with chlorotrifluoromethane as solvent. Commercial L-PLA precipitated after extraction of low molecular weight oligomers showed degradation kinetics similar to that of the starting material. The precipitation of DL-PLA from CO2 produced irregular-sized particles (10-20 microns). PGA, a polymer insoluble in most organic solvents, was found to be soluble in supercritical CO2. Nucleation of PGA from CO2 produced both regular-sized particles and needles of 10-40-microns length. The total solubility of commercial L-PLA in supercritical CO2 at 250 bar and 55 degrees C decreased from 0.14 wt % to less than 0.05 wt % and then leveled off as the cumulative flow of CO2 per unit mass of L-PLA loaded in the extractor increased beyond 20 standard L of CO2/g of L-PLA. Use of acetone (1 wt %) as a cosolvent increased L-PLA solubility by approximately 500%.