Volumetric measurements in binary solvents: theory to experiment

Volumetric Properties of Four-Stranded DNA Structures

Four-stranded non-canonical DNA structures including G-quadruplexes and i-motifs have been found in the genome and are thought to be involved in regulation of biological function. These structures have been implicated in telomere biology, genomic instability, and regulation of transcription and translation events. To gain an understanding of the molecular determinants underlying the biological role of four-stranded DNA structures, their biophysical properties have been extensively studied. The limited libraries on volume, expansibility, and compressibility accumulated to date have begun to provide insights into the molecular origins of helix-to-coil and helix-to-helix conformational transitions involving four-stranded DNA structures. In this article, we review the recent progress in volumetric investigations of G-quadruplexes and i-motifs, emphasizing how such data can be used to characterize intra-and intermolecular interactions, including solvation. We describe how volumetric data can be interpreted at the molecular level to yield a better understanding of the role that solute-solvent interactions play in modulating the stability and recognition events of nucleic acids. Taken together, volumetric studies facilitate unveiling the molecular determinants of biological events involving biopolymers, including G-quadruplexes and i-motifs, by providing one more piece to the thermodynamic puzzle describing the energetics of cellular processes in vitro and, by extension, in vivo.

Toward Biotherapeutics Formulation Composition Engineering using Site-Identification by Ligand Competitive Saturation (SILCS)

Formulation of protein-based therapeutics employ advanced formulation and analytical technologies for screening various parameters such as buffer, pH, and excipients. At a molecular level, physico-chemical properties of a protein formulation depend on self-interaction between protein molecules, protein-solvent and protein-excipient interactions. This work describes a novel in silico approach, SILCS-Biologics, for structure-based modeling of protein formulations. SILCS Biologics is based on the Site-Identification by Ligand Competitive Saturation (SILCS) technology and enables modeling of interactions among different components of a formulation at an atomistic level while accounting for protein flexibility. It predicts potential hotspot regions on the protein surface for protein-protein and protein-excipient interactions. Here we apply SILCS-Biologics on a Fab domain of a monoclonal antibody (mAbN) to model Fab-Fab interactions and interactions with three amino acid excipients, namely, arginine HCl, proline and lysine HCl. Experiments on 100 mg/ml formulations of mAbN showed that arginine increased, lysine reduced, and proline did not impact viscosity. We use SILCS-Biologics modeling to explore a structure-based hypothesis for the viscosity modulating effect of these excipients. Current efforts are aimed at further validation of this novel computational framework and expanding the scope to model full mAb and other protein therapeutics.

Duplex-tetraplex equilibria in guanine- and cytosine-rich DNA

Noncanonical four-stranded DNA structures, including G-quadruplexes and i-motifs, have been discovered in the cell and are implicated in a variety of genomic regulatory functions. The tendency of a specific guanine- and cytosine-rich region of genomic DNA to adopt a four-stranded conformation depends on its ability to overcome the constraints of duplex base-pairing by undergoing consecutive duplex-to-coil and coil-to-tetraplex transitions. The latter ability is determined by the balance between the free energies of participating ordered and disordered structures. In this review, we present an overview of the literature on the stability of G-quadruplex and i-motif structures and discuss the extent of duplex-tetraplex competition as a function of the sequence context of the DNA and environmental conditions including temperature, pH, salt, molecular crowding, and the presence of G-quadruplex-binding ligands. We outline how the results of in vitro studies can be expanded to understanding duplex-tetraplex equilibria in vivo.

Effect of cosolvent on protein stability: a theoretical investigation

We developed a statistical thermodynamic algorithm for analyzing solvent-induced folding/unfolding transitions of proteins. The energetics of protein transitions is governed by the interplay between the cavity formation contribution and the term reflecting direct solute-cosolvent interactions. The latter is viewed as an exchange reaction in which the binding of a cosolvent to a solute is accompanied by release of waters of hydration to the bulk. Our model clearly differentiates between the stoichiometric and non-stoichiometric interactions of solvent or co-solvent molecules with a solute. We analyzed the urea- and glycine betaine (GB)-induced conformational transitions of model proteins of varying size which are geometrically approximated by a sphere in their native state and a spherocylinder in their unfolded state. The free energy of cavity formation and its changes accompanying protein transitions were computed based on the concepts of scaled particle theory. The free energy of direct solute-cosolvent interactions were analyzed using empirical parameters previously determined for urea and GB interactions with low molecular weight model compounds. Our computations correctly capture the mode of action of urea and GB and yield realistic numbers for (∂ΔG°/∂a3)T,P which are related to the m-values of protein denaturation. Urea is characterized by negative values of (∂ΔG°/∂a3)T,P within the entire range of urea concentrations analyzed. At concentrations below ∼1 M, GB exhibits positive values of (∂ΔG°/∂a3)T,P which turn positive at higher GB concentrations. The balance between the thermodynamic contributions of cavity formation and direct solute-cosolvent interactions that, ultimately, defines the mode of cosolvent action is extremely subtle. A 20% increase or decrease in the equilibrium constant for solute-cosolvent binding may change the sign of (∂ΔG°/∂a3)T,P thereby altering the mode of cosolvent action (stabilizing to destabilizing or vice versa).

Caffeine dimerization: effects of sugar, salts, and water structure

Sugars and salts strongly affect the dimerization of caffeine in water. Such a change of dimerization, considered to be crucial for bitter taste suppression, has long been rationalized by the change of "water structure" induced by the additives; "kosmotropic" (water structure enhancing) salts and sugars promote dimerization, whereas "chaotropic" (water structure breaking) salts suppress dimerization. Based on statistical thermodynamics, here we challenge this consensus; we combine the rigorous Kirkwood-Buff theory of solution with the classical isodesmic model of caffeine association. Instead of the change of water structure, we show that the enhancement of caffeine dimerization is due to the exclusion of additives from caffeine, and that the weakening of dimerization is due to the binding of additives on caffeine.

Hydrotrope accumulation around the drug: the driving force for solubilization and minimum hydrotrope concentration for nicotinamide and urea

A rigorous statistical thermodynamic theory explains how urea and nicotinamide can solubilize hydrophobic drugs in water.

Infinitely dilute partial molar properties of proteins from computer simulation

A detailed understanding of temperature and pressure effects on an infinitely dilute protein's conformational equilibrium requires knowledge of the corresponding infinitely dilute partial molar properties. Established molecular dynamics methodologies generally have not provided a way to calculate these properties without either a loss of thermodynamic rigor, the introduction of nonunique parameters, or a loss of information about which solute conformations specifically contributed to the output values. Here we implement a simple method that is thermodynamically rigorous and possesses none of the above disadvantages, and we report on the method's feasibility and computational demands. We calculate infinitely dilute partial molar properties for two proteins and attempt to distinguish the thermodynamic differences between a native and a denatured conformation of a designed miniprotein. We conclude that simple ensemble average properties can be calculated with very reasonable amounts of computational power. In contrast, properties corresponding to fluctuating quantities are computationally demanding to calculate precisely, although they can be obtained more easily by following the temperature and/or pressure dependence of the corresponding ensemble averages.

Mechanism of hydrophobic drug solubilization by small molecule hydrotropes

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.