Drugs that are poorly soluble in water can be solubilized by the addition of hydrotropes. Albeit known for almost a century, how they work at a molecular basis is still controversial due to the lack of a rigorous theoretical basis. To clear up this situation, a combination of experimental data and Fluctuation Theory of Solutions (FTS) has been employed; information on the interactions between all the molecular species present in the solution has been evaluated directly. FTS has identified two major factors of hydrotrope-induced solubilization: preferential hydrotrope-solute interaction and water activity depression. The former is dominated by hydrotrope-solute association, and the latter is enhanced by ionic dissociation and hindered by the self-aggregation of the hydrotropes. Moreover, in stark contrast to previous hypotheses, neither the change of solute hydration nor the water structure accounts for hydrotropy. Indeed, the rigorous FTS poses serious doubts over the other common hypothesis: self-aggregation of the hydrotrope hinders, rather than promotes, solubilization.

Hydrophobic drugs can often be solubilized by the addition of hydrotropes. We have previously shown that preferential drug-hydrotrope association is one of the major factors of increased solubility (but not "hydrotrope clustering" or changes in "water structure"). How, then, can we understand this drug-hydrotrope interaction at a molecular level? Thermodynamic models based upon stoichiometric solute-water and solute-hydrotrope binding have long been used to understand solubilization microscopically. Such binding models have shown that the solvation numbers or coordination numbers of the water and hydrotrope molecules around the drug solute is the key quantity for solute-water and solute-hydrotrope interaction. However, we show that a rigorous statistical thermodynamic theory (the fluctuation solution theory originated by Kirkwood and Buff) requires the total reconsideration of such a paradigm. Here we show that (i) the excess solvation number (the net increase or decrease, relative to the bulk, of the solvent molecules around the solute), not the coordination number, is the key quantity for describing the solute-hydrotrope interaction; (ii) solute-hydrotrope binding is beyond the reach of the stoichiometric models because long-range solvation structure plays an important role.

Nicotinamide is an effective non-micellar hydrotrope (solubilizer) for drugs with low aqueous solubility. To clarify the molecular basis of nicotinamide's hydrotropic effectiveness, we present here a rigorous statistical thermodynamic theory, based on the Kirkwood-Buff theory of solutions, and our recent application of it to hydrotropy. We have shown that (i) nicotinamide self-association reduces solubilization efficiency, contrary to the previous hypothesis which claimed that self-association drives solubilization and (ii) the minimum hydrotrope concentration (MHC), namely, the threshold concentration above which solubility suddenly increases, is caused not by the bulk-phase self-association of nicotinamides as has been postulated previously, but by the enhancement of nicotinamide-nicotinamide interaction around the drug molecules. We have thus established a new view of hydrotropy - it is nicotinamide's non-stoichiometric accumulation around the drug that is the basis of solubility increase above MHC.

We all know that to enhance solubility using greener chemistry we should harness sound principles of molecular-based thermodynamics. The problem is that even for simple systems it can be hard to know how to use fundamental tools for formulation benefit, and for the more complex systems that we must often use, calculations required for molecular thermodynamics can often be quite involved. In this paper we show that a fundamental, assumption-free statistical thermodynamics approach, the Kirkwood-Buff theory, can be used in practical, complex aqueous systems to provide the insights we need to optimise formulations. The theory itself is not that difficult, but its implementation, which requires many steps of thermodynamic calculations, has up to now not been straightforward. Taking full advantage of an interactive approach, here we review what the Kirkwood-Buff theory can provide for formulators; we use the power of modern web browsers to provide open-source, user-friendly, responsive-design apps to do the hard work of data analysis, leaving formulators to focus on the interpretation of the results for their specific optimisation task. Indeed the apps are intended to be used by researchers and formulators for specific systems of interest to them.

In this paper, we briefly review the boxed molecular dynamics (BXD) method which allows analysis of thermodynamics and kinetics in complicated molecular systems. BXD is a multiscale technique, in which thermodynamics and long-time dynamics are recovered from a set of short-time simulations. In this paper, we review previous applications of BXD to peptide cyclization, solution phase organic reaction dynamics and desorption of ions from self-assembled monolayers (SAMs). We also report preliminary results of simulations of diamond etching mechanisms and protein unfolding in atomic force microscopy experiments. The latter demonstrate a correlation between the protein's structural motifs and its potential of mean force. Simulations of these processes by standard molecular dynamics (MD) is typically not possible, because the experimental time scales are very long. However, BXD yields well-converged and physically meaningful results. Compared with other methods of accelerated MD, our BXD approach is very simple; it is easy to implement, and it provides an integrated approach for simultaneously obtaining both thermodynamics and kinetics. It also provides a strategy for obtaining statistically meaningful dynamical results in regions of configuration space that standard MD approaches would visit only very rarely.