From structure and dynamics to biomolecular functions: The ubiquitous role of solvent in biology

Hydrogen-Bonded Network of Water in Phase-Separated Biomolecular Condensates

Biomolecular condensates formed via phase separation of intrinsically disordered proteins/regions (IDPs/IDRs) and nucleic acids are associated with cell physiology and disease. Water makes up for ∼60-70% of the condensate volume and is thought to influence the complex interplay of chain-chain and chain-solvent interactions, modulating the mesoscale properties of condensates. The behavior of water in condensates and the key roles of protein hydration water in driving the phase separation remain elusive. Here, we employ single-droplet vibrational Raman spectroscopy to illuminate the structural redistribution within protein hydration water during the phase separation of neuronal IDPs. Our Raman measurements reveal the changes in the water hydrogen bonding network during homotypic and heterotypic phase separation governed by various molecular drivers. Such single-droplet water Raman measurements offer a potent generic tool to unmask the intriguing interplay of sequence-encoded chain-chain and chain-solvent interactions governing macromolecular phase separation into membraneless organelles, synthetic condensates, and protocells.

Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

Entropy Tug-of-War Determines Solvent Effects in the Liquid-Liquid Phase Separation of a Globular Protein

Liquid-liquid phase separation (LLPS) plays a key role in the compartmentalization of cells via the formation of biomolecular condensates. Here, we combined atomistic molecular dynamics (MD) simulations and terahertz (THz) spectroscopy to determine the solvent entropy contribution to the formation of condensates of the human eye lens protein γD-Crystallin. The MD simulations reveal an entropy tug-of-war between water molecules that are released from the protein droplets and those that are retained within the condensates, two categories of water molecules that were also assigned spectroscopically. A recently developed THz-calorimetry method enables quantitative comparison of the experimental and computational entropy changes of the released water molecules. The strong correlation mutually validates the two approaches and opens the way to a detailed atomic-level understanding of the different driving forces underlying the LLPS.

Force-induced unzipping of DNA in the presence of solvent molecules

The melting of double-stranded DNA (dsDNA) in the presence of solvent molecules is a fundamental process with significant implications for understanding the thermal and mechanical behavior of DNA and its interactions with the surrounding environment. The solvents play an essential role in the structural transformation of DNA subjected to a pulling force. In this study, we simulate the thermal and force induced denaturation of dsDNA and elucidate the solvent dependent melting behavior, identifying key factors that influence the stability of DNA melting in presence of solvent molecules. Using a statistical model, we first find the melting profile of short heterogeneous DNA molecules in the presence of solvent molecules in Force ensemble. We also investigate the effect of solvent's strengths on the melting profile of DNA. In the force ensemble, we consider two homogeneous DNA chains and apply the force on different locations along the chain in the presence of solvent molecules. Different pathways manifest the melting of the molecule in both ensembles, and we found several interesting features of melting DNA in a constant force ensemble, such as lower critical force when the chain is pulled from the base pair close to a solvent molecule. The results provide new insights into the force-induced unzipping of DNA and could be used to develop new methods for controlling the unzipping process. By providing a better understanding of melting and unzipping of dsDNA in the presence of solvent molecules, this study provides valuable guidelines for predicting DNA thermodynamic quantities and for designing DNA nanostructures.

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

Gastric cancer specific drug delivery with hydrophilic peptide probe conjugation

Cancer-specific diagnosis is challenging. Phage display is an approach that could contribute to finding new specific biomarkers. In this study, we developed a new peptide probe specific for gastric cancer and validated it for gastric cancer-specific theranostics. We isolated linear peptides by screening a combinatorial phage library for a cancer stem cell marker, LGR5 protein. Among these, peptides with high selectivity against gastric cancer cells were selected and examined for therapeutic poteintial in vitro as well as in vivo. Through leucine-rich G protein-coupled receptor 5 (LGR5) protein-based phage display, we obtained a hydrophilic 7-mer peptide sequence (STCTRSR, named STC). Both the STC-peptide-conjugated fluorescent dye and chlorin e6 (Ce6) displayed a significantly higher intensity in gastric cancer cells compared to that in healthy cells. In mice with gastric cancer, the fluorescence in the tumors was 3.4× more intense when treated with the Ce6-STC conjugate compared to that with free Ce6 and conferred higher phototoxicity after single laser irradiation. Repeated photodynamic therapy could further reduce the tumor volume after treating these mice with the Ce6-STC conjugate. The treatment with the Ce6-STC conjugate exhibited a significantly lower fluorescence in the liver than that with free Ce6. In conclusion, we confirmed that the STC peptide is a gastric cancer-specific probe that could be useful in gastric cancer theranostics. In conclusion, considering its targeting ability and hydrophilicity, various hydrophobic chemotherapeutic agents could be revisited for gastric cancer treatment in combination with the probe described in this study.

In Vitro and In Vivo Trypanocidal Efficacy of Nitrofuryl- and Nitrothienylazines

African trypanosomiasis is a vector-borne disease of animals and humans in the tsetse fly belt of Africa. Trypanosoma congolense ("nagana") is the most pathogenic trypanosome in livestock and causes high morbidity and mortality rates among cattle. In the absence of effective preventative vaccines, the management of trypanosomiasis relies on chemoprophylaxis and/or -therapy. However, the trypanocides in clinical use exhibit poor oral bioavailability and toxicity, and therapeutic failures occur because of resistant strains. Because nitrofurantoin displayed, in addition to its clinical use, promising antiparasitic activity, the current study was conducted to evaluate the in vitro trypanocidal activity and preliminary in vivo treatment efficacy of previously synthesized nitrofuranylazines. The trypanocidal activity of these nitrofuran derivatives varied among the evaluated trypanosome species; however, T. congolense strain IL3000 was more susceptible than other animal and human trypanosomes. The nitrofurylazines 4a (IC50 0.04 μM; SI > 7761) and 7a (IC50 0.03 μM; SI > 9542) as well as the nitrothienylazine 8b (IC50 0.04 μM; SI 232), with nanomolar IC50 values, were revealed as early antitrypanosomal leads. Although these derivatives showed strong trypanocidal activity in vitro, no in vivo treatment efficacy was observed in T. congolense IL3000 infected mice after both oral and intraperitoneal administration in a preliminary study. This was attributed to the poor solubility of the test compounds in the in vivo testing media. Indeed, a challenge in drug discovery is finding a balance between the physicochemical properties of a drug candidate, particularly lipophilicity and water solubility, and maintaining adequate potency to provide an effective dose. Hence, future chemical modifications may be required to generate lead-like to lead-like nitrofuranylazines that possess optimal physicochemical and pharmacokinetic properties while retaining in vitro and, ultimately, in vivo trypanocidal efficacy.

Thermodynamic forces from protein and water govern condensate formation of an intrinsically disordered protein domain

Liquid-liquid phase separation (LLPS) can drive a multitude of cellular processes by compartmentalizing biological cells via the formation of dense liquid biomolecular condensates, which can function as membraneless organelles. Despite its importance, the molecular-level understanding of the underlying thermodynamics of this process remains incomplete. In this study, we use atomistic molecular dynamics simulations of the low complexity domain (LCD) of human fused in sarcoma (FUS) protein to investigate the contributions of water and protein molecules to the free energy changes that govern LLPS. Both protein and water components are found to have comparably sizeable thermodynamic contributions to the formation of FUS condensates. Moreover, we quantify the counteracting effects of water molecules that are released into the bulk upon condensate formation and the waters retained within the protein droplets. Among the various factors considered, solvation entropy and protein interaction enthalpy are identified as the most important contributions, while solvation enthalpy and protein entropy changes are smaller. These results provide detailed molecular insights on the intricate thermodynamic interplay between protein- and solvation-related forces underlying the formation of biomolecular condensates.

Dynamics of Hydrogen Bonds between Water and Intrinsically Disordered and Structured Regions of Proteins

Recent studies indicate more restricted dynamics of water around intrinsically disordered proteins (IDPs) than structured proteins. We examine here the dynamics of hydrogen bonds between water molecules and two proteins, small ubiquitin-related modifier-1 (SUMO-1) and ubiquitin-conjugating enzyme E2I (UBC9), which we compare around intrinsically disordered regions (IDRs) and structured regions of these proteins. It has been recognized since some time that excluded volume effects, which influence access of water molecules to hydrogen-bonding sites, and the strength of hydrogen bonds between water and protein affect hydrogen bond lifetimes. While we find those two properties to mediate lifetimes of hydrogen bonds between water and protein residues in this study, we also find that the lifetimes are affected by the concentration of charged groups on other nearby residues. These factors are more important in determining the hydrogen bond lifetimes than whether a residue hydrogen bonding with water belongs to an IDR or to a structured region.

Curcuminoids as Anticancer Drugs: Pleiotropic Effects, Potential for Metabolic Reprogramming and Prospects for the Future

The number of published studies on curcuminoids in cancer research, including its lead molecule curcumin and synthetic analogs, has been increasing substantially during the past two decades. Insights on the diversity of inhibitory effects they have produced on a multitude of pathways involved in carcinogenesis and tumor progression have been provided. As this wealth of data was obtained in settings of various experimental and clinical data, this review first aimed at presenting a chronology of discoveries and an update on their complex in vivo effects. Secondly, there are many interesting questions linked to their pleiotropic effects. One of them, a growing research topic, relates to their ability to modulate metabolic reprogramming. This review will also cover the use of curcuminoids as chemosensitizing molecules that can be combined with several anticancer drugs to reverse the phenomenon of multidrug resistance. Finally, current investigations in these three complementary research fields raise several important questions that will be put among the prospects for the future research related to the importance of these molecules in cancer research.

Bimodal 1/f Noise and Anticorrelation between DNA-Water and DNA-Ion Energy Fluctuations

The coupling between the conformational fluctuations of DNA and its surrounding environment, consisting of water and ions in solution, remains poorly understood and relatively less investigated as compared to proteins. Here, with the help of molecular dynamics simulations and statistical mechanical analyses, we explore the dynamical coupling among DNA, water, and counterions through correlations among respective energy fluctuations in both double- (ds-) and single-stranded (ss-) DNA solutions. Fluctuations in the collective DNA-water and DNA-ion interaction energies are found to be strongly anticorrelated across all the systems. The fluctuations of DNA self-energy, however, are weakly coupled to DNA-water and DNA-ion interactions in ds-DNA. An enhancement of the DNA-water coupling is observed in ss-DNA, where the system is less rigid. All the interaction energies exhibit 1/f noise in their energy power spectra with surprisingly prominent bimodality in the DNA-water and DNA-ion fluctuations. The nature of the energy spectra appears to be indifferent to the relative rigidity of the DNA. We discuss the role of the observed correlations in ion-water motions on a DNA duplex in the experimentally observed anomalous slow dielectric relaxation and solvation dynamics and in furthering our understanding of the DNA energy landscape.

Locating dynamic contributions to allostery via determining rates of vibrational energy transfer

Determining rates of energy transfer across non-covalent contacts for different states of a protein can provide information about dynamic and associated entropy changes during transitions between states. We investigate the relationship between rates of energy transfer across polar and nonpolar contacts and contact dynamics for the β₂-adrenergic receptor, a rhodopsin-like G-protein coupled receptor, in an antagonist-bound inactive state and agonist-bound active state. From structures sampled during molecular dynamics (MD) simulations, we find the active state to have, on average, a lower packing density, corresponding to generally more flexibility and greater entropy than the inactive state. Energy exchange networks (EENs) are computed for the inactive and active states from the results of the MD simulations. From the EENs, changes in the rates of energy transfer across polar and nonpolar contacts are found for contacts that remain largely intact during activation. Change in dynamics of the contact, and entropy associated with the dynamics, can be estimated from the change in rates of energy transfer across the contacts. Measurement of change in the rates of energy transfer before and after the transition between states thereby provides information about dynamic contributions to activation and allostery.

Semiclassical Theory of Multistage Nonequilibrium Electron Transfer in Macromolecular Compounds in Polar Media with Several Relaxation Timescales

Many specific features of ultrafast electron transfer (ET) reactions in macromolecular compounds can be attributed to nonequilibrium configurations of intramolecular vibrational degrees of freedom and the environment. In photoinduced ET, nonequilibrium nuclear configurations are often produced at the stage of optical excitation, but they can also be the result of electron tunneling itself, i.e., fast redistribution of charges within the macromolecule. A consistent theoretical description of ultrafast ET requires an explicit consideration of the nuclear subsystem, including its evolution between electron jumps. In this paper, the effect of the multi-timescale nuclear reorganization on ET transitions in macromolecular compounds is studied, and a general theory of ultrafast ET in non-Debye polar environments with a multi-component relaxation function is developed. Particular attention is paid to designing the multidimensional space of nonequilibrium nuclear configurations, as well as constructing the diabatic free energy surfaces for the ET states. The reorganization energies of individual ET transitions, the equilibrium energies of ET states, and the relaxation properties of the environment are used as input data for the theory. The effect of the system-environment interaction on the ET kinetics is discussed, and mechanisms for enhancing the efficiency of charge separation in macromolecular compounds are analyzed.