The effect of cosolvents on biomolecular equilibria has traditionally been rationalized using simple binding models. More recently, a renewed interest in the use of Kirkwood-Buff (KB) theory to analyze solution mixtures has provided new information on the effects of osmolytes and denaturants and their interactions with biomolecules. Here we review the status of KB theory as applied to biological systems. In particular, the existing models of denaturation are analyzed in terms of KB theory, and the use of KB theory to interpret computer simulation data for these systems is discussed.

We introduce and validate a new precision oncology framework for the systematic prioritization of drugs targeting mechanistic tumor dependencies in individual patients. Compounds are prioritized on the basis of their ability to invert the concerted activity of master regulator proteins that mechanistically regulate tumor cell state, as assessed from systematic drug perturbation assays. We validated the approach on a cohort of 212 gastroenteropancreatic neuroendocrine tumors (GEP-NETs), a rare malignancy originating in the pancreas and gastrointestinal tract. The analysis identified several master regulator proteins, including key regulators of neuroendocrine lineage progenitor state and immunoevasion, whose role as critical tumor dependencies was experimentally confirmed. Transcriptome analysis of GEP-NET-derived cells, perturbed with a library of 107 compounds, identified the HDAC class I inhibitor entinostat as a potent inhibitor of master regulator activity for 42% of metastatic GEP-NET patients, abrogating tumor growth in vivo. This approach may thus complement current efforts in precision oncology.

A molecular dynamics study of the effects of salt on the conformational preferences of leucine enkephalin (LE) in water has been performed to outline a general approach for quantifying the effects of a cosolvent on a biomolecular equilibrium in solution using molecular simulation. The simulations of LE suggest that the peptide exists as a mixture of folded and unfolded conformations in pure water, and a potential of mean force calculation indicates the population of the folded form to be 53±3%. The addition of 2 M salt reduces the population of the folded form to 12±5%. A combination of the simulation data and Kirkwood−Buff theory is used to determine the preferential interactions of the cosolvent with the peptide and to reproduce the population changes indicated by the potential of mean force calculations. The results demonstrate the potential of a combined simulation and Kirkwood−Buff approach for quantifying simulation data in an effort to understand the effects of cosolvents on the thermodynamics of biomolecules in solution, especially for systems where potential of mean force calculations are unfeasible.

Opioid receptors, like other members of the G protein-coupled receptor (GPCR) family, have been shown to associate to form dimers and/or oligomers at the plasma membrane. Whether this association is stable or transient is not known. Recent compelling evidence suggests that at least some GPCRs rapidly associate and dissociate. We have recently calculated binding affinities from free energy estimates to predict transient association between mouse delta opioid receptor (DOR) protomers at a symmetric interface involving the fourth transmembrane (TM4) helix (herein termed “4” dimer). Here we present disulfide cross-linking experiments with DOR constructs with cysteines substituted at the extracellular ends of TM4 or TM5 that confirm the formation of DOR complexes involving these helices. Our results are consistent with the involvement of TM4 and/or TM5 at the DOR homodimer interface, but possibly with differing association propensities. Coarse-grained (CG) well-tempered metadynamics simulations of two different dimeric arrangements of DOR involving TM4 alone or with TM5 (herein termed “4/5” dimer) in an explicit lipid−water environment confirmed the presence of two structurally and energetically similar configurations of the 4 dimer, as previously assessed by umbrella sampling calculations, and revealed a single energetic minimum of the 4/5 dimer. Additional CG umbrella sampling simulations of the 4/5 dimer indicated that the strength of association between DOR protomers varies depending on the protein region at the interface, with the 4 dimer being more stable than the 4/5 dimer.