Investigation of molecular and mesoscale clusters in undersaturated glycine aqueous solutions

The influence of solute concentration and filtration on mesoscale clusters of flufenamic acid, a typical pharmaceutical compound, in ethanol

HYPOTHESIS

It is hypothesised in this work that mesoscale clusters will be present in both undersaturated and supersaturated solutions of organic pharmaceutical molecules. These clusters, being loose aggregates, could be sensitive to shear forces experienced during filtration. Thus, comparing the behaviour of these clusters alongside nanoparticles during filtration-an important sample treatment parameter during crystallization-will elucidate qualitative differences from solid, crystalline nanoparticles of similar size.

EXPERIMENTS

The impact of filtration with different pore sizes and membranes on (i) mesoscale clusters of flufenamic acid (FFA) ethanol solutions and (ii) aqueous FFA nanosuspensions was studied with dynamic light scattering and nanoparticle tracking analysis.

FINDINGS

FFA solutions, ranging from undersaturated to supersaturated, were found to form mesoscale clusters, where the cluster size and number concentration were independent of solute concentration. Under filtration stress, irrespective of pore size and membrane used, the mesoscale cluster peak disappeared from the size distribution with no detectable change in concentration. In contrast, similarly sized FFA nanoparticles were removed by filtration, causing a significant change in solute concentration and size distribution. Mesoscale clusters of FFA in ethanol constitute only a tiny fraction of the total solute concentration and possess poor light scattering properties, lower mass density than solid particles of similar size, and no clear phase boundary.

Amorphous aggregates with a very wide size distribution play a central role in crystal nucleation

There is mounting evidence that crystal nucleation from supersaturated solution involves the formation and reorganization of prenucleation clusters, contradicting classical nucleation theory. One of the key unresolved issues pertains to the origin, composition, and structure of these clusters. Here, a range of amino acids and peptides is investigated using light scattering, mass spectrometry, and in situ terahertz Raman spectroscopy, showing that the presence of amorphous aggregates is a general phenomenon in supersaturated solutions. Significantly, these aggregates are found on a vast range of length scales from dimers to 30-mers to the nanometre and even micrometre scale, implying a continuous distribution throughout this range. Larger amorphous aggregates are sites of spontaneous crystal nucleation and act as intermediates for laser-induced crystal nucleation. These results are shown to be consistent with a nonclassical nucleation model in which barrierless (homogeneous) nucleation of amorphous aggregates is followed by the nucleation of crystals from solute-enriched aggregates. This provides a novel perspective on crystal nucleation and the role of nonclassical pathways.

Mesoscale Clusters in the Crystallisation of Organic Molecules

The early stages of the molecular self-assembly pathway leading to crystal nucleation have a significant influence on the properties and purity of organic materials. This mini review collates the work on organic mesoscale clusters and discusses their importance in nucleation processes, with a particular focus on their critical properties and susceptibility to sample treatment parameters. This is accomplished by a review of detection methods, including dynamic light scattering, nanoparticle tracking analysis, small angle X-ray scattering, and transmission electron microscopy. Considering the challenges associated with crystallisation of flexible and large-molecule active pharmaceutical ingredients, the dynamic nature of mesoscale clusters has the potential to expand the discovery of novel crystal forms. By collating literature on mesoscale clusters for organic molecules, a more comprehensive understanding of their role in nucleation will evolve and can guide further research efforts.

Amino acid chiral amplification using Monte Carlo dynamic

This study investigates the stability of chiral-molecule solution phases, with a specific focus on amino acids. The model framework is based on a two-dimensional square lattice model, where individual sites may be occupied by oriented chiral molecules or structureless solvent particles. Utilizing the Glauber dynamics and statistical mechanical formalism, as previously introduced and examined by Lombardo et al., we explore the influence of temperature, amino acid concentration, enantiomeric excess, and homochiral interaction strength on nucleation mechanisms, equilibrium phase behavior, and crystal composition. Our findings offer thermodynamic insights into the chiral amplification process of amino acids, contributing to a deeper understanding of the underlying processes.

Molecular Dynamics Investigation of Giant Clustering in Small-Molecule Solutions: The Case of Aqueous PEHA

The importance of the formation of giant clusters in solution, in nature and industry, is increasingly recognized. However, relatively little attention has been paid to the formation of giant clusters in solutions of small, relatively soluble but nonamphiphilic molecules. In this work, we present a general methodology based on molecular dynamics that can be used to investigate such systems. As a case study, we focus on the formation of apparently stable clusters of pentaethylenehexamine (PEHA) in water. These clusters have been used as templates for the construction of bioinspired silica nanoparticles. To better understand clustering in this system, we study the effect of PEHA protonation state (neutral, +1, and +2) and counterion type (chloride or acetate) on PEHA clustering in dilute aqueous solutions (200 and 400 mM) using large-scale classical molecular dynamics. We find that large stable clusters are formed by singly charged PEHA with chloride or acetate as the counterion, although it is not clear for the case with acetate whether bulk phase separation, that might lead to precipitation, would eventually occur. Large clusters also appear to be stable for doubly charged PEHA with acetate, the less soluble counterion. We attribute this behavior to a form of complex coacervation, observed here for relatively small and highly soluble molecules (PEHA + counterion) rather than the large polyions usually found to form such coacervates. We discuss whether this behavior might also be described by an effective SALR (short-range attraction, long-range repulsion) interaction. This work might help future studies of additives for the design of novel bioinspired templated nanomaterials and of giant clustering in small-molecule solutions more generally.

Tuning Interfacial Concentration Enhancement through Dispersion Interactions to Facilitate Heterogeneous Nucleation

Classical molecular dynamics simulations were used to investigate how dispersion (van der Waals) interactions between non-polar, hydrophobic surfaces and aqueous glycine solutions affect the solution composition, molecular orientation, and dynamics at the interface. Simulations revealed that dispersion interactions lead to a major increase in the concentration of glycine at the interface in comparison with the bulk solution, resulting from a competition between solute and solvent molecules to be or not to be near the interface. This can then lead to kinetic and/or structural effects facilitating heterogeneous nucleation of glycine at non-polar surfaces, in agreement with recent observations for tridecane, graphene, and polytetrafluoroethylene. A novel parameterization process was developed to map a model surface with tunable dispersion interactions to heptane, tridecane, and graphite materials. The model surface was capable of reproducing the solution structure observed in fully atomistic simulations with excellent agreement and also provided good agreement for dynamic properties, at a significantly reduced computational cost. This approach can be used as an effective tool for screening materials for heterogeneous nucleation enhancement or suppression, based on non-specific dispersion interactions based on bulk material molecular properties, rather than interfacial functional groups, templating or confinement effects.

Molecular Dynamics Investigation of Clustering in Aqueous Glycine Solutions

Recent experiments with undersaturated aqueous glycine solutions have repeatedly exhibited the presence of giant liquid-like clusters or nanodroplets around 100 nm in diameter. These nanodroplets re-appear even after careful efforts for their removal and purification of the glycine solution. The composition of these clusters is not clear, although it has been suggested that they are mainly composed of glycine, a small and very soluble amino acid. To gain insights into this phenomenon, we study the aggregation of glycine in aqueous solutions at concentrations below the experimental solubility limit using large-scale molecular dynamics simulations under ambient conditions. Three protonation states of glycine (zwitterion = GLZ, anion = GLA, and cation = GLC) are simulated using molecular force fields based on the 1.14*CM1A partial charge scheme, which incorporates the OPLS all-atom force field and TIP3P water. When initiated from dispersed states, we find that giant clusters do not form in our simulations unless salt impurities are present. Moreover, if simulations are initiated from giant cluster states, we find that they tend to dissolve in the absence of salt impurities. Therefore, the simulation results provide little support for the possibility that the giant clusters seen in experiments are composed purely of glycine (and water). Considering that strenuous efforts are made in experiments to remove impurities such as salt, we propose that the giant clusters observed might instead result from the aggregation of reaction products of aqueous glycine, such as diketopiperazine or other oligoglycines which may be difficult to separate from glycine using conventional methods, or their co-aggregation with glycine.

Fluorescence of KCl Aqueous Solution: A Possible Spectroscopic Signature of Nucleation

Ion pairing in water solutions alters both the water hydrogen-bond network and ion solvation, modifying the dynamics and properties of electrolyte water solutions. Here, we report an anomalous intrinsic fluorescence of KCl aqueous solution at room temperature and show that its intensity increases with the salt concentration. From the ab initio density functional theory (DFT) and time-dependent DFT modeling, we propose that the fluorescence emission could originate from the stiffening of the hydrogen bond network in the hydration shell of solvated ion-pairs that suppresses the fast nonradiative decay and allows the slower radiative channel to become a possible decay pathway. Because computations suggest that the fluorophores are the local ion-water structures present in the prenucleation phase, this band could be the signature of the incoming salt precipitation.

Mesoscale clusters of organic solutes in solution and their role in crystal nucleation

It is becoming evident that primary nucleation of crystals of organic molecules from solution is often anything but ‘classical’ in its complexity. It is also becoming increasingly clear that mesoscopic...