How Can Hydrophobic Association Be Enthalpy Driven?

Capturing the Effects of Single Atom Substitutions on the Inhibition Efficiency of Glycogen Synthase Kinase-3β Inhibitors via Markov State Modeling and Experiments

Small modifications in the chemical structure of ligands are known to dramatically change their ability to inhibit the activity of a protein. Unraveling the mechanisms that govern these dramatic changes requires scrutinizing the dynamics of protein-ligand binding and unbinding at the atomic level. As an exemplary case, we have studied Glycogen Synthase Kinase-3β (GSK-3β), a multifunctional kinase that has been implicated in a host of pathological processes. As such, there is a keen interest in identifying ligands that inhibit GSK-3β activity. One family of compounds that are highly selective and potent inhibitors of GSK-3β is exemplified by a molecule termed COB-187. COB-187 consists of a five-member heterocyclic ring with a thione at C2, a pyridine substituted methyl at N3, and a hydroxyl and phenyl at C4. We have studied the inhibition of GSK-3β by COB-187-related ligands that differ in a single heavy atom from each other (either in the location of nitrogen in their pyridine ring, or with the pyridine ring replaced by a phenyl ring), or in the length of the alkyl group joining the pyridine and the N3. The inhibition experiments show a large range of half-maximal inhibitory concentration (IC50) values from 10 nM to 10 μM, implying that these ligands exhibit vastly different propensities to inhibit GSK-3β. To explain these differences, we perform Markov State Modeling (MSM) using fully atomistic simulations. Our MSM results are in excellent agreement with the experiments in that they accurately capture differences in the binding propensities of the ligands. The simulations show that the binding propensities are related to the ligands' ability to attain a compact conformation where their two aromatic rings are spatially close. We rationalize this result by sampling numerous binding and unbinding events via funnel metadynamics simulations, which show that indeed while approaching the bound state, the ligands prefer to be in their compact conformation. We find that the presence of nitrogen in the aromatic ring increases the probability of attaining the compact conformation. Protein-ligand binding is understood to be dictated by the energetics of interactions and entropic factors, like the release of bound water from the binding pockets. This work shows that changes in the conformational distribution of ligands due to atom-level modifications in the structure play an important role in protein-ligand binding.

Unveiling the dynamic and thermodynamic interactions of hydrocortisone with β-cyclodextrin and its methylated derivatives through insights from molecular dynamics simulations

Cyclodextrins (CDs) can enhance the stability and bioavailability of pharmaceutical compounds by encapsulating them within their cavities. This study utilized molecular dynamics simulations to investigate the interaction mechanisms between hydrocortisone (HC) and various methylated CD derivatives. The results reveal that the loading of HC into CD cavities follows different mechanisms depending on the degree and position of methylation. Loading into βCD and 6-MeβCD was more complete, with the hydroxyl groups of HC facing the primary hydroxyl rim (PHR) and the ketone side facing the secondary hydroxyl rim (SHR). In contrast, 2,3-D-MeβCD and 2,6-D-MeβCD showed a different loading mechanism, with the ketone side facing the PHR and the hydroxyl groups facing the SHR. The root mean square fluctuation (RMSF) analysis demonstrated that methylation increases the flexibility of CD heavy atoms, with 3-MeβCD and 2,3-D-MeβCD exhibiting the highest flexibility. However, upon inclusion of HC, 3-MeβCD, 2,3-D-MeβCD, 2-MeβCD, and 6-MeβCD showed a significant reduction in flexibility, suggesting a more rigid structure that effectively retains HC within their cavities. The radial distribution function revealed a significant reduction in the number of water molecules within the innermost layer of the methylated CD cavities, particularly in TMeβCD, indicating a decrease in polarity. The presence of HC led to the release of high-energy water molecules, creating more favorable conditions for HC loading. Conformational analysis showed that methylation caused a partial decrease in the area of the PHR, a significant decrease in the area of the middle rim, and a notable decrease in the area of the SHR. The loading of HC increased the area of the PHR in most derivatives, with the most pronounced increase observed in 2,6-D-MeβCD and 6-MeβCD. The analysis of interaction energies and binding free energies demonstrated that the binding of HC to methylated CD derivatives is thermodynamically more favorable than to βCD, with the strongest association observed for 6-MeβCD, 2-MeβCD, and 2,3-D-MeβCD.

Solvent effects in anion recognition

Anion recognition is pertinent to a range of environmental, medicinal and industrial applications. Recent progress in the field has relied on advances in synthetic host design to afford a broad range of potent recognition motifs and novel supramolecular structures capable of effective binding both in solution and at derived molecular films. However, performance in aqueous media remains a critical challenge. Understanding the effects of bulk and local solvent on anion recognition by host scaffolds is imperative if effective and selective detection in real-world media is to be viable. This Review seeks to provide a framework within which these effects can be considered both experimentally and theoretically. We highlight proposed models for solvation effects on anion binding and discuss approaches to retain strong anion binding in highly competitive (polar) solvents. The synthetic design principles for exploiting the aforementioned solvent effects are explored.

Protein structure-function continuum model: Emerging nexuses between specificity, evolution, and structure

The rationale for replacing the old binary of structure-function with the trinity of structure, disorder, and function has gained considerable ground in recent years. A continuum model based on the expanded form of the existing paradigm can now subsume importance of both conformational flexibility and intrinsic disorder in protein function. The disorder is actually critical for understanding the protein-protein interactions in many regulatory processes, formation of membrane-less organelles, and our revised notions of specificity as amply illustrated by moonlighting proteins. While its importance in formation of amyloids and function of prions is often discussed, the roles of intrinsic disorder in infectious diseases and protein function under extreme conditions are also becoming clear. This review is an attempt to discuss how our current understanding of protein function, specificity, and evolution fit better with the continuum model. This integration of structure and disorder under a single model may bring greater clarity in our continuing quest for understanding proteins and molecular mechanisms of their functionality.

A Prebiotic Precursor to Life's Phosphate Transfer System with an ATP Analog and Histidyl Peptide Organocatalysts

Biochemistry is dependent upon enzyme catalysts accelerating key reactions. At the origin of life, prebiotic chemistry must have incorporated catalytic reactions. While this would have yielded much needed amplification of certain reaction products, it would come at the possible cost of rapidly depleting the high energy molecules that acted as chemical fuels. Biochemistry solves this problem by combining kinetically stable and thermodynamically activated molecules (e.g., ATP) with enzyme catalysts. Here, we demonstrate a prebiotic phosphate transfer system involving an ATP analog (imidazole phosphate) and histidyl peptides, which function as organocatalytic enzyme analogs. We demonstrate that histidyl peptides catalyze phosphorylations via a phosphorylated histidyl intermediate. We integrate these histidyl-catalyzed phosphorylations into a complete prebiotic scenario whereby inorganic phosphate is incorporated into organic compounds though physicochemical wet-dry cycles. Our work demonstrates a plausible system for the catalyzed production of phosphorylated compounds on the early Earth and how organocatalytic peptides, as enzyme precursors, could have played an important role in this.

Identifying Potential Ligand Binding Sites on Glycogen Synthase Kinase 3 Using Atomistic Cosolvent Simulations

Glycogen synthase kinase 3 β (GSK3β) is a serine/threonine kinase that phosphorylates several protein substrates in crucial cell signaling pathways. Owing to its therapeutic importance, there is a need to develop GSK3β inhibitors with high specificity and potency. One approach is to find small molecules that can allosterically bind to the GSK3β protein surface. We have employed fully atomistic mixed-solvent molecular dynamics (MixMD) simulations to identify three plausible allosteric sites on GSK3β that can facilitate the search for allosteric inhibitors. Our MixMD simulations narrow down the allosteric sites to precise regions on the GSK3β surface, thereby improving upon the previous predictions of the locations of these sites.

Thermodynamically Optimized Machine-Learned Reaction Coordinates for Hydrophobic Ligand Dissociation

Ligand unbinding is mediated by its free energy change, which has intertwined contributions from both energy and entropy. It is important, but not easy, to quantify their individual contributions to the free energy profile. We model hydrophobic ligand unbinding for two systems, a methane particle and a C60 fullerene, both unbinding from hydrophobic pockets in all-atom water. Using a modified deep learning framework, we learn a thermodynamically optimized reaction coordinate to describe the hydrophobic ligand dissociation for both systems. Interpretation of these reaction coordinates reveals the roles of entropic and enthalpic forces as the ligand and pocket sizes change. In both cases, we observe that the free-energy barrier to unbinding is dominated by entropy considerations. Furthermore, the process of methane unbinding is driven by methane solvation, while fullerene unbinding is driven first by pocket wetting and then fullerene wetting. For both solutes, the direct importance of the distance from the binding pocket to the learned reaction coordinate is present, but low. Our framework and subsequent feature important analysis thus give useful thermodynamic insight into hydrophobic ligand dissociation problems that are otherwise difficult to glean.

Exploring the role of cyclodextrins as a cholesterol scavenger: a molecular dynamics investigation of conformational changes and thermodynamics

This study presents a comprehensive analysis of the cholesterol binding mechanism and conformational changes in cyclodextrin (CD) carriers, namely βCD, 2HPβCD, and MβCD. The results revealed that the binding of cholesterol to CDs was spontaneous and thermodynamically favorable, with van der Waals interactions playing a dominant role, while Coulombic interactions have a negligible contribution. The solubility of cholesterol/βCD and cholesterol/MβCD complexes was lower compared to cholesterol/2HPβCD complex due to stronger vdW and Coulombic repulsion between water and CDs. Hydrogen bonding was found to have a minor role in the binding process. The investigation of mechanisms and kinetics of binding demonstrated that cholesterol permeates into the CD cavities completely. Replicas consideration indicated that while the binding to 2HPβCD occurred perpendicularly and solely through positioning cholesterol's oxygen toward the primary hydroxyl rim (PHR), the mechanism of cholesterol binding to βCD and MβCD could take place with the orientation of oxygen towards both rims. Functionalization resulted in decreased cavity polarity, increased constriction tendency, and altered solubility and configuration of the carrier. Upon cholesterol binding, the CDs expanded, increasing the cavity volume in cholesterol-containing systems. The effects of cholesterol on the relative shape anisotropy (κ2) and asphericity parameter (b) in cyclodextrins were investigated. βCD exhibited a spherical structure regardless of cholesterol presence, while 2HPβCD and MβCD displayed more pronounced non-sphericity in the absence of cholesterol. Loading cholesterol transformed 2HPβCD and MβCD into more spherical shapes, with increased probabilities of higher κ2. MβCD showed a higher maximum peak of κ2 compared to 2HPβCD after cholesterol loading, while 2HPβCD maintained a significant maximum peak at 0.2 for b.

Biomedical applications of supramolecular hydrogels with enhanced mechanical properties

Supramolecular hydrogels bound by hydrogen bonding, host-guest, hydrophobic, and other non-covalent interactions are among the most attractive biomaterials available. Supramolecular hydrogels have attracted extensive attention due to their inherent dynamic reversibility, self-healing, stimuli-response, excellent biocompatibility, and near-physiological environment. However, the inherent contradiction between non-covalent interactions and mechanical strength makes the practical application of supramolecular hydrogels a great challenge. This review describes the mechanical strength of hydrogels mediated by supramolecular interactions, and focuses on the potential strategies for enhancing the mechanical strength of supramolecular hydrogels and illustrates their applications in related fields, such as flexible electronic sensors, wound dressings, and three-dimensional (3D) scaffolds. Finally, the current problems and future research prospects of supramolecular hydrogels are discussed. This review is expected to provide insights that will motivate more advanced research on supramolecular hydrogels.

Structural and interaction analysis of human lipocalin-type prostaglandin D synthase with the poorly water-soluble drug NBQX

Lipocalin-type prostaglandin D synthase (L-PGDS) is a secretory lipid-transporter protein that was shown to bind a wide variety of hydrophobic ligands in vitro. Exploiting this function, we previously examined the feasibility of using L-PGDS as a novel delivery vehicle for poorly water-soluble drugs. However, the mechanism by which human L-PGDS binds to poorly water-soluble drugs is unclear. In this study, we determined the solution structure of human L-PGDS and investigated the mechanism of L-PGDS binding to 6-nitro-7-sulfamoyl-benzo[f]quinoxalin-2,3-dione (NBQX), an α-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid receptor antagonist. NMR experiments showed that human L-PGDS has an eight-stranded antiparallel β-barrel structure that forms a central cavity, a short 310 -helix and two α-helices. Titration with NBQX was monitored using 1 H-15 N HSQC spectroscopy. At higher NBQX concentrations, some cross-peaks of the protein exhibited fast-exchanging shifts with a curvature, indicating at least two binding sites. These residues were located in the upper portion of the cavity. Singular value decomposition analysis revealed that human L-PGDS has two NBQX binding sites. Large chemical shift changes were observed in the H2-helix and A-, B-, C-, D-, H- and I-strands and H2-helix upon NBQX binding. Calorimetric experiments revealed that human L-PGDS binds two NBQX molecules with dissociation constants of 46.7 μm for primary binding and 185.0 μm for secondary binding. Molecular docking simulations indicated that these NBQX binding sites are located within the β-barrel. These results provide new insights into the interaction between poorly water-soluble drugs and human L-PGDS as a drug carrier.

The Role of Carbohydrate Intake on the Gut Microbiome: A Weight of Evidence Systematic Review

(1) Background: Carbohydrates are the most important source of nutritional energy for the human body. Carbohydrate digestion, metabolism, and their role in the gut microbiota modulation are the focus of multiple studies. The objective of this weight of evidence systematic review is to investigate the potential relationship between ingested carbohydrates and the gut microbiota composition at different taxonomic levels. (2) Methods: Weight of evidence and information value techniques were used to evaluate the relationship between dietary carbohydrates and the relative abundance of different bacterial taxa in the gut microbiota. (3) Results: The obtained results show that the types of carbohydrates that have a high information value are: soluble fiber with Bacteroides increase, insoluble fiber with Bacteroides and Actinobacteria increase, and Firmicutes decrease. Oligosaccharides with Lactobacillus increase and Enterococcus decrease. Gelatinized starches with Prevotella increase. Starches and resistant starches with Blautia decrease and Firmicutes increase. (4) Conclusions: This work provides, for the first time, an integrative review of the subject by using statistical techniques that have not been previously employed in microbiota reviews.

Insights from molecular dynamics and DFT calculations into the interaction of 1,4-benzodiazepines with 2-hydroxypropyl-βCD in a theoretical study

This study delves into the interaction between benzodiazepine (BZD) drugs and 2-hydroxypropyl-β-cyclodextrin (2HPβCD), a cyclodextrin (CD) known to improve drug delivery and enhance therapeutic outcomes. We find that the 2HPβCD's atoms become more rigid in the presence of chlordiazepoxide (CDP), clonazepam (CLZ), and diazepam (DZM), whereas they become more flexible in the presence of nordazepam (NDM) and nitrazepam (NZP). We also investigated the structure of 2HPβCD and found that loading these drugs increases both the area and volume of the 2HPβCD cavity, making it more suitable for drug delivery. Moreover, this research found that all drugs exhibited negative values for the binding free energy, indicating thermodynamic favorability and improved solubility. The binding free energy order of the BZDs was consistent in both molecular dynamics and Monte Carlo methods, with CDP and DZM having the highest affinity for binding. We also analyzed the contribution of different interaction energies in binding between the carrier and the drugs and found that Van der Waals energy is the primary component. Our results indicate that the number of hydrogen bonds between 2HPβCD/water slightly decreases in the presence of BZDs, but the hydrogen bond's quality remains constant.

Elucidation of a Thermodynamical Feature Attributed to Substrate Binding to the Prokaryotic H+/Oligopeptide Cotransporter YdgR with Calorimetric Analysis: The Substrate Binding Driven by the Change in Entropy Implies the Release of Bound Water Molecules from the Binding Pocket

Here, we have elucidated the substrate recognition mechanism by a prokaryotic H⁺/oligopeptide cotransporter, YdgR, using isothermal titration calorimetry. Under acidic conditions (pH 6.0), the binding of a dipeptide, Val-Ala, to YdgR elicited endothermic enthalpy, which compensated for the increase in entropy due to dipeptide binding. A series of dipeptides were used in the binding titration. The dipeptides represent Val-X and X-Val, where X is Ala, Ser, Val, Tyr, or Phe. Most dipeptides revealed endothermic enthalpy, which was completely compensated by the increase in entropy due to dipeptide binding. The change in enthalpy due to binding correlated well with the change in entropy, whereas the Gibbs free energy involved in the binding of the dipeptide to YdgR remained unchanged irrespective of dipeptide sequences, implying that the binding reaction was driven by entropy, that is, the release of bound water molecules in the binding pocket. It is also important to clarify that, based on the prediction of water molecules in the ligand-binding pocket of YdgR, the release of three bound water molecules in the putative substrate binding pocket occurred through binding to YdgR. In the comparison of Val-X and X-Val dipeptides, the N-terminal region of the binding pocket might contain more bound water molecules than the C-terminal region. In light of these findings, we suggest that bound water molecules might play an important role in substrate recognition and binding by YdgR.

Determination of partition coefficients of phthalic acid esters between polydimethylsiloxane and water and its field application to surface waters

Phthalic acid esters (PAEs) or phthalates are endocrine-disrupting chemicals and among the most frequently detected hydrophobic organic pollutants, which can be gradually released from consumer products into the environment (e.g., water). This study measured the equilibrium partition coefficients for 10 selected PAEs, with a wide range of logarithms of the octanol-water partition coefficient (log K_ow) from 1.60 to 9.37, between poly(dimethylsiloxane) (PDMS) and water (K_PDMSw) using the kinetic permeation method. The desorption rate constant (k_d) and K_PDMSw for each PAEs were calculated from kinetic data. The experimental log K_PDMSw for the PAEs ranges from 0.8 to 5.9, which is linearly correlated with log K_ow values up to 8 from the literature (R² > 0.94); however, it slightly deviated for the PAEs with log K_ow values greater than 8. In addition, K_PDMSw decreased with the temperature and enthalpy for PAEs partitioning in PDMS-water in an exothermic manner. Furthermore, the effects of dissolved organic matter and ionic strength on the partitioning of PAEs in PDMS were investigated. PDMS was used as a passive sampler to determine the aqueous concentration of plasticizers in river surface water. The results of this study can be used to evaluate the bioavailability and risk of phthalates in real environmental samples.

Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization

Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying "magic bullets" amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand-target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural-energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Physical insight into the entropy-driven ion association

The ion association is widely believed to be dominated by the favorable entropy change arising from the release of water molecules from ion hydration shells. However, no direct thermodynamic evidence exists to validate the reliability and suitability of this view. Herein, we employ complicated free energy calculations to rigorously split the free energy including its entropic and enthalpic components into the water-induced contributions and ion-ion interaction terms for several ion pairs from monatomic to polyatomic ions, spanning the size range from small kosmotropes to large chaotropes (Na⁺ , Cs⁺ , Ca²⁺ , F^- , I^- , CO₃ ^2- , and HPO₄ ^2- ). Our results successfully reveal that though ion associations are indeed determined by a delicate balance between the favorable entropy variation and the repulsive enthalpy change, the entropy gain dominated by the solvent occurs only for the monatomic ion pairing. The water-induced entropic contribution significantly goes against the ion pairing between polyatomic anion and cation, which is, alternatively, dominated by the favorable entropy from the ion-ion interaction term, due to the configurational arrangement of polyatomic anions involved in ion association. The structural and dynamic analysis demonstrates that the entropy penalty from the water phase is primarily ascribed to the enhanced stability of water molecules around the cation imposed by the incoming anion. Our study successfully provides a fundamental understanding of water-mediated ion associations and highlights disparate lengthscale dependencies of the dehydration thermodynamics on the specific types of ions.

Solubility of Hydrophobes into Macrocyclic Hosts

The intrinsic nature of macrocyclic molecules to preferentially absorb a specific solute has been opening up supramolecular chemistry. Nevertheless, the determinant factor with molecular perspectives in promoting host-guest complexations remains inconclusive, due to the lack of rigorous thermodynamic examination on the guest solubility inside the host. Here, we quantify the solute-solvent energetic and entropic contributions between the end states and on the docking route during inclusion of noble gases in cucurbit[5]uril, cucurbit[6]uril, and α-cyclodextrin, using molecular dynamics simulations in combination with the potential distribution theorem. Results show that in all of the pairs examined both the solute-solvent energy and entropy favor the inclusion, while the former is rather dominant. The frequency of interior drying, which pertains to the entropic contribution, differs between the hosts and is controlled by the existence of lid water at portal and the flexibility of host framework. Moreover, the hosts exhibit various types of absorption manners, involving non-, single-, and double-free-energy barriers.

Enhancing Sampling of Water Rehydration on Ligand Binding: A Comparison of Techniques

Water often plays a key role in protein structure, molecular recognition, and mediating protein-ligand interactions. Thus, free energy calculations must adequately sample water motions, which often proves challenging in typical MD simulation time scales. Thus, the accuracy of methods relying on MD simulations ends up limited by slow water sampling. Particularly, as a ligand is removed or modified, bulk water may not have time to fill or rearrange in the binding site. In this work, we focus on several molecular dynamics (MD) simulation-based methods attempting to help rehydrate buried water sites: BLUES, using nonequilibrium candidate Monte Carlo (NCMC); grand, using grand canonical Monte Carlo (GCMC); and normal MD. We assess the accuracy and efficiency of these methods in rehydrating target water sites. We selected a range of systems with varying numbers of waters in the binding site, as well as those where water occupancy is coupled to the identity or binding mode of the ligand. We analyzed the rehydration of buried water sites in binding pockets using both clustering of trajectories and direct analysis of electron density maps. Our results suggest both BLUES and grand enhance water sampling relative to normal MD and grand is more robust than BLUES, but also that water sampling remains a major challenge for all of the methods tested. The lessons we learned for these methods and systems are discussed.

Cholesterol's Condensing Effect: Unpacking a Century-Old Mystery

The ability of cholesterol to uncoil (i.e., condense) the acyl chains of phospholipids has been known for a century. Despite extensive studies of the interactions between cholesterol and phospholipids, a molecular-level understanding of this uncoiling phenomenon has remained elusive. Equally unclear has been whether cholesterol's two different faces (i.e., its relatively smooth α face and its relatively rough β face) contribute to its condensing power. Because cholesterol's condensing effect is believed to play a major role in controlling the fluidity, structure, and functioning of all animal cell membranes, its biological importance cannot be overstated. This Perspective focuses on experimental evidence that addresses (i) the credibility of a popular "umbrella" mechanism that has been used to account for cholesterol's condensing effect, (ii) the credibility of an alternate "template" mechanism, (iii) the importance of cholesterol two-faced character with respect to its condensing power, and (iv) the viability of a surface occupancy model.

Fruit Bromelain-Derived Peptide Potentially Restrains the Attachment of SARS-CoV-2 Variants to hACE2: A Pharmacoinformatics Approach

Before entering the cell, the SARS-CoV-2 spike glycoprotein receptor-binding domain (RBD) binds to the human angiotensin-converting enzyme 2 (hACE2) receptor. Hence, this RBD is a critical target for the development of antiviral agents. Recent studies have discovered that SARS-CoV-2 variants with mutations in the RBD have spread globally. The purpose of this in silico study was to determine the potential of a fruit bromelain-derived peptide. DYGAVNEVK. to inhibit the entry of various SARS-CoV-2 variants into human cells by targeting the hACE binding site within the RBD. Molecular docking analysis revealed that DYGAVNEVK interacts with several critical RBD binding residues responsible for the adhesion of the RBD to hACE2. Moreover, 100 ns MD simulations revealed stable interactions between DYGAVNEVK and RBD variants derived from the trajectory of root-mean-square deviation (RMSD), radius of gyration (Rg), and root-mean-square fluctuation (RMSF) analysis, as well as free binding energy calculations. Overall, our computational results indicate that DYGAVNEVK warrants further investigation as a candidate for preventing SARS-CoV-2 due to its interaction with the RBD of SARS-CoV-2 variants.

Systematic Evaluation of Fluorination as Modification for Peptide-Based Fusion Inhibitors against HIV-1 Infection

With the emergence of novel viruses, the development of new antivirals is more urgent than ever. A key step in human immunodeficiency virus type 1 (HIV-1) infection is six-helix bundle formation within the envelope protein subunit gp41. Selective disruption of bundle formation by peptides has been shown to be effective; however, these drugs, exemplified by T20, are prone to rapid clearance from the patient. The incorporation of non-natural amino acids is known to improve these pharmacokinetic properties. Here, we evaluate a peptide inhibitor in which a critical Ile residue is replaced by fluorinated analogues. We characterized the influence of the fluorinated analogues on the biophysical properties of the peptide. Furthermore, we show that the fluorinated peptides can block HIV-1 infection of target cells at nanomolar levels. These findings demonstrate that fluorinated amino acids are appropriate tools for the development of novel peptide therapeutics.

Chitosan-Based Biocompatible Copolymers for Thermoresponsive Drug Delivery Systems: On the Development of a Standardization System

Chitosan is a natural polysaccharide that is considered to be biocompatible, biodegradable and non-toxic. The polymer has been used in drug delivery applications for its positive charge, which allows for adhesion with and recognition of biological tissues via non-covalent interactions. In recent times, chitosan has been used for the preparation of graft copolymers with thermoresponsive polymers such as poly-N-vinylcaprolactam (PNVCL) and poly-N-isopropylamide (PNIPAM), allowing the combination of the biodegradability of the natural polymer with the ability to respond to changes in temperature. Due to the growing interest in the utilization of thermoresponsive polymers in the biological context, it is necessary to increase the knowledge of the key principles of thermoresponsivity in order to obtain comparable results between different studies or applications. In the present review, we provide an overview of the basic principles of thermoresponsivity, as well as a description of the main polysaccharides and thermoresponsive materials, with a special focus on chitosan and poly-N-Vinyl caprolactam (PNVCL) and their biomedical applications.

Near-Instantaneously Self-Healing Coating toward Stable and Durable Electromagnetic Interference Shielding

Durable electromagnetic interference (EMI) shielding is highly desired, as electromagnetic pollution is a great concern for electronics' stable performance and human health. Although a superhydrophobic surface can extend the service lifespan of EMI shielding materials, degradation of its protection capability and insufficient self-healing are troublesome issues due to unavoidable physical/chemical damages under long-term application conditions. Here, we report, for the first time, an instantaneously self-healing approach via microwave heating to achieve durable shielding performance. First, a hydrophobic 1H,1H,2H,2H-perfluorooctyltriethoxysilane (POTS) layer was coated on a polypyrrole (PPy)-modified fabric (PPy@POTS), enabling protection against the invasion of water, salt solution, and corrosive acidic and basic solutions. Moreover, after being damaged, the POTS layer can, for the first time, be instantaneously self-healed via microwave heating for a very short time, i.e., 4 s, benefiting from the intense thermal energy generated by PPy under electromagnetic wave radiation. This self-healing ability is also repeatable even after intentionally severe plasma etching, which highlights the great potential to achieve robust and durable EMI shielding applications. Significantly, this approach can be extended to other EMI shielding materials where heat is a triggering stimulus for healing thin protection layers. We envision that this work could provide insights into fabricating EMI shielding materials with durable performance for portable and wearable devices, as well as for human health care.

Native and substituted cyclodextrins as chiral selectors for capillary electrophoresis enantioseparations: Structures, features, application, and molecular modeling

CDs are cyclic oligosaccharides consisting of α-d-glucopyranosyl units linked through 1,4-linkages, which are obtained from enzymatic degradation of starch. The coexistence of hydrophilic and hydrophobic regions in the same structure makes these macrocycles extremely versatile as complexing host with application in food, cosmetics, environmental, agriculture, textile, pharmaceutical, and chemical industries. Due to their inherent chirality, CDs have been also successfully used as chiral selectors in enantioseparation science, in particular, for CE enantioseparations. In the last decades, multidisciplinary approaches based on CE, NMR spectroscopy, X-ray crystallography, microcalorimetry, and molecular modeling have shed light on some aspects of recognition mechanisms underlying enantiodiscrimination. With the ever growing improvement of computer facilities, hardware and software, computational techniques have become a useful tool to model at molecular level the dynamics of diastereomeric associate formation to sample low-energy conformations, the binding energies between the enantiomer and the CD, and to profile noncovalent interactions contributing to the stability of CD/enantiomer association. On this basis, the aim of this review is to provide the reader with a critical overview on the applications of CDs in CE. In particular, the contemporary theory of the electrophoretic technique and the main structural features of CDs are described, with a specific focus on techniques, methods, and approaches to model CE enantioseparations promoted by native and substituted CDs. A systematic compilation of all published literature has not been attempted.

Cavitand Complexes in Aqueous Solution: Collaborative Experimental and Computational Studies of the Wetting, Assembly, and Function of Nanoscopic Bowls in Water

Water is the dominant liquid on Earth. Despite this, the main focus of supramolecular chemistry research has been on binding and assembly events in organic solvents. This arose because it is more straightforward to synthesize organic-media-soluble hosts and because of the relative simplicity of organic solvents compared to water. Nature, however, relies on water as a solvent, and spurred by this fact, supramolecular chemists have recently been making forays into the aqueous domain to understand water-mediated non-covalent interactions. These studies can benefit from the substantial understanding of the hydrophobic effect and electrostatic interactions developed by physical chemists. Nearly 20 years ago, the Gibb group first synthesized a class of water-soluble host molecules, the deep-cavity cavitands, that possess non-polar pockets that readily bind non-polar moieties in aqueous solution and are capable of assembling into a wide range of complexes with distinct stoichiometries. As such, these amphipathic host species are ideal platforms for studying the role of negatively curved features on guest complexation and the structural requirements for guided assembly processes driven by the hydrophobic effect. Here we review the collaborative experimental and computational investigations between Gibb and Ashbaugh over the past 10 years exploring questions including the following: How does water wet/solvate the non-polar surfaces of non-polar pockets? How does this wetting control the binding of non-polar guests? How does wetting affect the binding of anionic species? How does the nature and size of a guest size impact the assembly of cavitand hosts into multimeric capsular complexes? What are the conformational motifs of guests packed within the confines of capsular complexes? How might the electrostatic environment engendered by hosts impact the properties and reactivity of internalized guests?

Protein flexibility reduces solvent-mediated friction barriers of ligand binding to a hydrophobic surface patch

Solvent fluctuations have been explored in detail for idealized and rigid hydrophobic model systems, but so far it has remained unclear how internal protein motions and their coupling to the surrounding solvent affect the dynamics of ligand binding to biomolecular surfaces. Here, molecular dynamics simulations were used to elucidate the solvent-mediated binding of a model ligand to the hydrophobic surface patch of ubiquitin. The ligand's friction profiles reveal pronounced long-time correlations and enhanced friction in the vicinity of the protein, similar to idealized hydrophobic surfaces. Interestingly, these effects are shaped by internal protein motions. Protein flexibility modulates water density fluctuations near the hydrophobic surface patch and smooths out the friction profile of ligand binding.

Time-dependent conformational analysis of ALK5-lumican complex in presence of graphene and graphene oxide employing molecular dynamics and MMPBSA calculation

Lumican, an extracellular matrix protein avails wound healing by binding to ALK5 membrane receptor (TGF-beta receptor I). Their interaction enables epithelialization and substantiates rejuvenation of injured tissue. To enrich permanence of ALK5-lumican interaction, we employed graphene and graphene oxide co-factors. Herein, this study explicates concomitancy of graphene and graphene oxide with ALK5-lumican. We performed an in silico approach involving molecular modelling, molecular docking, molecular dynamics for 200 ns, DSSP analysis and MMPBSA calculations. Results of molecular dynamics indicate cofactors influential in altering bioactive site of lumican than ALK5. Similarly, MMPBSA calculations unveiled binding energy of apoenzyme as -108.09 kcal/mol, holoenzyme (G) as -79.20 kcal/mol and holoenzyme (GO) as -114.33 kcal/mol. This concludes graphene oxide lucrative in enhancing binding energy of ALK5-lumican in holoenzyme (GO) via coil formation of Lum C₁₃ domain. In contrast, graphene reduced binding energy of ALK5-lumican in holoenzyme (G) modifying Lum C₁₃ into beta sheets. MMPBSA residual contribution analysis of Lum C₁₃ residues revealed binding energy of -13.9 kcal/mol for apoenzyme, -6.8 kcal/mol for holoenzyme (G) and -19.5 kcal/mol for holoenzyme (GO). This supports coil formation propitious for better ALK5-Lum interaction. Highest SASA energy of -21.05 kcal/mol of holoenzyme (G) assures graphene reasonable for improved ALK5-lumican hydrophobicity. As per the motive of the study, graphene oxide enriches permanence of ALK5-lumican. This provides counsel for plausible exploitation of lumican and graphene oxide as targeted/nano drug delivery system to reinstate acute wounds, chronic wounds, corneal wounds, hypertrophic scars and keloids in near future. Communicated by Ramaswamy H. Sarma.

Water-Mediated Interactions Determine Helix Formation of Peptides in Open Nanotubes

Water-mediated interactions (WMIs) play diverse roles in molecular biology. They are particularly relevant in geometrically confined spaces such as the interior of the chaperonin, at the interface between ligands and their binding partners, and in the ribosome tunnel. Inspired in part by the geometry of the ribosome tunnel, we consider confinement effects on the stability of peptides. We describe results from replica exchange molecular dynamics simulations of a system containing a 23-alanine or 23-serine polypeptide confined to nonpolar and polar nanotubes in the gas phase and when open to a water reservoir. We quantify the effect of water in determining the preferred conformational states of these polypeptides by calculating the difference in the solvation free energy for the helix and coil states in the open nanotube in the two phases. Our simulations reveal several possibilities. We find that nanoscopic confinement preferentially stabilizes the helical state of polypeptides with hydrophobic side chains, which is explained by the entropic stabilization mechanism proposed on the basis of polymer physics. Polypeptide chains with hydrophilic side chains can adopt helical structures within nanotubes, but helix formation is sensitive to the nature of the nanotube due to WMIs. We elaborate on the potential implications of our findings to the stability of peptides in the ribosome tunnel.

Tuning Local Hydration Enables a Deeper Understanding of Protein-Ligand Binding: The PP1-Src Kinase Case

Water plays a key role in biomolecular recognition and binding. Despite the development of several computational and experimental approaches, it is still challenging to comprehensively characterize water-mediated effects on the binding process. Here, we investigate how water affects the binding of Src kinase to one of its inhibitors, PP1. Src kinase is a target for treating several diseases, including cancer. We use biased molecular dynamics simulations, where the hydration of predetermined regions is tuned at will. This computational technique efficiently accelerates the SRC-PP1 binding simulation and allows us to identify several key and yet unexplored aspects of the solvent's role. This study provides a further perspective on the binding phenomenon, which may advance the current drug design approaches for the development of new kinase inhibitors.

Selective Affimers Recognise the BCL-2 Family Proteins BCL-xL and MCL-1 through Noncanonical Structural Motifs*

The BCL-2 family is a challenging group of proteins to target selectively due to sequence and structural homologies across the family. Selective ligands for the BCL-2 family regulators of apoptosis are useful as probes to understand cell biology and apoptotic signalling pathways, and as starting points for inhibitor design. We have used phage display to isolate Affimer reagents (non-antibody-binding proteins based on a conserved scaffold) to identify ligands for MCL-1, BCL-xL , BCL-2, BAK and BAX, then used multiple biophysical characterisation methods to probe the interactions. We established that purified Affimers elicit selective recognition of their target BCL-2 protein. For anti-apoptotic targets BCL-xL and MCL-1, competitive inhibition of their canonical protein-protein interactions is demonstrated. Co-crystal structures reveal an unprecedented mode of molecular recognition; where a BH3 helix is normally bound, flexible loops from the Affimer dock into the BH3 binding cleft. Moreover, the Affimers induce a change in the target proteins towards a desirable drug-bound-like conformation. These proof-of-concept studies indicate that Affimers could be used as alternative templates to inspire the design of selective BCL-2 family modulators and more generally other protein-protein interaction inhibitors.

Molecular View into the Cyclodextrin Cavity: Structure and Hydration

We find, through atomistic molecular dynamics simulation of native cyclodextrins (CDs) in water, that although the outer surface of a CD appears like a truncated cone, the inner cavity resembles a conical hourglass because of the inward protrusion of the glycosidic oxygens. Furthermore, the conformations of the constituent α-glucose molecules are found to differ significantly from a free monomeric α-glucose molecule. This is the first computational study that maps the conformational change to the preferential hydrogen bond donating capacity of one of the secondary hydroxyl groups of CD, in consensus with an NMR experiment. We have developed a simple and novel geometry-based technique to identify water molecules occupying the nonspherical CD cavity, and the computed water occupancies are in close agreement with the experimental and density functional theory studies. Our analysis reveals that a water molecule in CD cavity loses out about two hydrogen bonds and remains energetically frustrated but possesses higher orientational degree of freedom compared to bulk water. In the context of CD-drug complexation, these imply a nonclassical, that is, enthalpically driven hydrophobic association of a drug in CD cavity.

Entropy-Enthalpy Compensation in Peptide Adsorption on Solid Surfaces: Dependence on Surface Hydration

Although protein adsorption at the solid-water interface is of immense importance, understanding the crucial role of the water phase in mediating protein-surface interactions is lacking, particularly due to the lack of fundamental thermodynamic data. Herein, we have performed complicated free energy calculations and successfully extracted the entropy and enthalpy changes of molecular adsorption on solids. Using the gold and graphene as the surface models with distinct affinities to the water phase, we successfully unravel the sharply opposite manners of entropy-enthalpy compensation in driving water and tripeptide adsorptions on two surfaces. Though the thermodynamic features of water adsorption on surface are enthalpically dominated based on the positions of free energy barriers and minima, the favorable entropy term significantly decreases the free energy barrier and further stabilizes the adsorbate at the adsorption site on the graphene surface. For the peptide, the shape of the adsorption free energy profile is jointly determined by the enthalpy and entropy changes, which, however, alternatively act the driving force to promote the peptide adsorption on the Au surface and graphene surface. The distinct structural and dynamic properties of solid-liquid interfaces account for the special role of the interfacial water phase in regulating the competitive relationship between the entropy and enthalpy variations.

Supertough Lignin Hydrogels with Multienergy Dissipative Structures and Ultrahigh Antioxidative Activities

Hydrogels derived from lignin are typically weak and contain only a small amount of lignin, which limits their broad application prospects. In the present work, a novel lignin/poly(N,N-dimethylacrylamide) (PDMA) hydrogel with a high lignin content, superb toughness, and ultrahigh antioxidative performance is constructed by employing a facile dissolve-dry-swell solvent exchange method. Through this process, lignin and PDMA are self-assembled into a multienergy dissipative structure containing rigid lignin-rich domains. Precisely, the PDMA chains both interpenetrated inside and adhered on the surface of these domains through hydrophobic associations. This structure enables the lignin hydrogels to dissipate energy efficiently during the fracture process. At an optimized ultrahigh lignin content of 58% (dry weight basis), the prepared lignin hydrogel exhibited remarkable mechanical properties, such as a high elastic modulus (2.5 MPa), tensile strength (2.5 MPa), and super tensile strain (11.3), and an extremely high fracture energy above 16 000 J m^-2. In addition, the tough lignin hydrogel exhibited a commendable antioxidant property and nontoxicity. All these advantageous properties provide the lignin/PDMA hydrogels with the potential for use in biomedical materials applications.

Distinct binding of cetirizine enantiomers to human serum albumin and the human histamine receptor H1

Cetirizine, a major metabolite of hydroxyzine, became a marketed second-generation H₁ antihistamine that is orally active and has a rapid onset of action, long duration of effects and a very good safety record at recommended doses. The approved drug is a racemic mixture of (S)-cetirizine and (R)-cetirizine, the latter being the levorotary enantiomer that also exists in the market as a third-generation, non-sedating and highly selective antihistamine. Both enantiomers bind tightly to the human histamine H₁ receptor (hH₁R) and behave as inverse agonists but the affinity and residence time of (R)-cetirizine are greater than those of (S)-cetirizine. In blood plasma, cetirizine exists in the zwitterionic form and more than 90% of the circulating drug is bound to human serum albumin (HSA), which acts as an inactive reservoir. Independent X-ray crystallographic work has solved the structure of the hH₁R:doxepin complex and has identified two drug-binding sites for cetirizine on equine serum albumin (ESA). Given this background, we decided to model a membrane-embedded hH₁R in complex with either (R)- or (S)-cetirizine and also the complexes of both ESA and HSA with these two enantiomeric drugs to analyze possible differences in binding modes between enantiomers and also among targets. The ensuing molecular dynamics simulations in explicit solvent and additional computational chemistry calculations provided structural and energetic information about all of these complexes that is normally beyond current experimental possibilities. Overall, we found very good agreement between our binding energy estimates and extant biochemical and pharmacological evidence. A much higher degree of solvent exposure in the cetirizine-binding site(s) of HSA and ESA relative to the more occluded orthosteric binding site in hH₁R is translated into larger positional fluctuations and considerably lower affinities for these two nonspecific targets. Whereas it is demonstrated that the two known pockets in ESA provide enough stability for cetirizine binding, only one such site does so in HSA due to a number of amino acid replacements. At the histamine-binding site in hH₁R, the distinct interactions established between the phenyl and chlorophenyl moieties of the two enantiomers with the amino acids lining up the pocket and between their free carboxylates and Lys179 in the second extracellular loop account for the improved pharmacological profile of (R)-cetirizine.

From Specific γ-CD/[Nb6 Cl12 (H2 O)6 ]2+ Recognition to Biological Activity Tuning

Specific molecular recognition of γ-cyclodextrin (γ-CD) by the cationic hexanuclear niobium [Nb₆ Cl₁₂ (H₂ O)₆ ]²⁺ cluster complex in aqueous solutions results in a 1:1 supramolecular assembly {[Nb₆ Cl₁₂ (H₂ O)₆ ]@γ-CD}²⁺ . NMR spectroscopy, isothermal titration calorimetry (ITC), and ESI-MS were used to study the interaction between the inorganic cluster and the organic macrocycle. Such molecular association affects the biological activity of [Nb₆ Cl₁₂ (H₂ O)₆ ]²⁺ , decreasing its cytotoxicity despite enhanced cellular uptake. The 1:1 stoichiometry is maintained in solution over a large window of the reagents' ratio, but crystallization by slow evaporation produces a 1:2 host-guest complex [Nb₆ Cl₁₂ (H₂ O)₆ @(γ-CD)₂ ]Cl₂ ⋅20 H₂ O featuring the cluster encapsulated between two molecules of γ-CD. The 1:2 complex was characterized by XRD, elemental analysis, IR spectroscopy, and thermogravimetric analysis (TGA). Quantum chemical calculations were performed to model host-guest interaction.

Pressure Induced Wetting and Dewetting of the Nonpolar Pocket of Deep-Cavity Cavitands in Water

Hydrophobic interactions drive the binding of nonpolar ligands to the oily pockets of proteins and supramolecular species in aqueous solution. As such, the wetting of host pockets is expected to play a critical role in determining the thermodynamics of guest binding. Here we use molecular simulations to examine the impact of pressure on the wetting and dewetting of the nonpolar pockets of a series of deep-cavity cavitands in water. The portals to the cavitand pockets are functionalized with both nonpolar (methyl) and polar (hydroxyl) groups oriented pointing either upward or inward toward the pocket. We find wetting of the pocket is favored by the hydroxyl groups and dewetting is favored by the methyl groups. The distribution of waters in the pocket is found to exhibit a two-state-like equilibrium between wet and dry states with a free energy barrier between the two states. Moreover, we demonstrate that the pocket hydration of the cavitands can be collapsed onto a unified adsorption isotherm by assuming the effective pressures within each cavitand pocket differ by a shift pressure that depends on the chemical identity and number of functional groups placed about the portal. These observations support the development of a two-state capillary evaporation model that accurately describes the equilibrium between states and naturally gives rise to the effective shift pressures observed from simulation. This work demonstrates that the hydration of host pockets can be tuned following simple design rules that in turn are expected to impact the thermodynamics of guest complexation.

Spontaneous drying of non-polar deep-cavity cavitand pockets in aqueous solution

There are many open questions regarding the hydration of solvent-exposed non-polar tracts and pockets in proteins. Although water is predicted to de-wet purely repulsive surfaces and evacuate crevices, the extent of de-wetting is unclear when ubiquitous van der Waals interactions are in play. The structural simplicity of synthetic supramolecular hosts imbues them with considerable potential to address this issue. To this end, here we detail a combination of densimetry and molecular dynamics simulations of three cavitands, coupled with calorimetric studies of their complexes with short-chain carboxylates. Our results reveal the range of wettability possible within the ostensibly identical cavitand pockets-which differ only in the presence and/or position of the methyl groups that encircle the portal to their non-polar pockets. The results demonstrate the ability of macrocycles to template water cavitation within their binding sites and show how the orientation of methyl groups can trigger the drying of non-polar pockets in liquid water, which suggests new avenues to control guest complexation.

Structure-thermodynamics-relationships of hepatitis C viral NS3 protease inhibitors

Thermodynamic parameters were determined for structurally-related inhibitors of HCV NS3 protease to assess how binding entropies and enthalpies vary with incremental changes at the P2 and P3 inhibitor subsites. Changing the heterocyclic substituent at P2 from a pyridyl to a 7-methoxy-2-phenyl-4-quinolyl group leads to a 710-fold increase in affinity. Annelating a benzene ring onto a pyridine ring leads to quinoline-derived inhibitors having higher affinities, but the individual enthalpy and entropy contributions are markedly different for each ligand pair. Introducing a phenyl group at C2 of the heterocyclic ring at P2 uniformly leads to higher affinity analogs with more favorable binding entropies, while adding a methoxy group at C7 of the quinoline ring at P2 provides derivatives with more favorable binding enthalpies. Significant enthalpy/entropy compensation is observed for structural changes made to inhibitors lacking a 2-phenyl substituent, whereas favorable changes in both binding enthalpies and entropies accompany structural modifications when a 2-phenyl group is present. Overall, binding energetics of inhibitors having a 2-phenyl-4-quinolyl group at P2 are dominated by entropic effects, whereas binding of the corresponding norphenyl analogs are primarily enthalpy driven. Notably, the reversal from an entropy driven association to an enthalpy driven one for this set of inhibitors also correlates with alternate binding modes. When the steric bulk of the side chain at P3 is increased from a hydrogen atom to a tert-butyl group, there is a 770-fold improvement in affinity. The 30-fold increase resulting from the first methyl group is solely the consequence of a more favorable change in entropy, whereas subsequent additions of methyl groups leads to modest increases in affinity that arise primarily from incremental improvements in binding enthalpies accompanied with smaller favorable entropic contributions.

Proximal charge effects on guest binding to a non-polar pocket

Science still does not have the ability to accurately predict the affinity that ligands have for proteins. In an attempt to address this, the Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) series of blind predictive challenges is a community-wide exercise aimed at advancing computational techniques as standard predictive tools in rational drug design. In each cycle, a range of biologically relevant systems of different levels of complexity are selected to test the latest modeling methods. As part of this on-going exercise, and as a step towards understanding the important factors in context dependent guest binding, we challenged the computational community to determine the affinity of a series of negatively and positively charged guests to two constitutionally isomeric cavitand hosts: octa-acid 1, and exo-octa acid 2. Our affinity determinations, combined with molecular dynamics simulations, reveal asymmetries in affinities between host-guest pairs that cannot alone be explained by simple coulombic interactions, but also point to the importance of host-water interactions. Our work reveals the key facets of molecular recognition in water, emphasizes where improvements need to be made in modelling, and shed light on the complex problem of ligand-protein binding in the aqueous realm.

Water structure in solution and crystal molecular dynamics simulations compared to protein crystal structures

The function of proteins is influenced not only by the atomic structure but also by the detailed structure of the solvent surrounding it. Computational studies of protein structure also critically depend on the water structure around the protein. Herein we compare the water structure obtained from molecular dynamics (MD) simulations of galectin-3 in complex with two ligands to crystallographic water molecules observed in the corresponding crystal structures. We computed MD trajectories both in a water box, which mimics a protein in solution, and in a crystallographic unit cell, which mimics a protein in a crystal. The calculations were compared to crystal structures obtained at both cryogenic and room temperature. Two types of analyses of the MD simulations were performed. First, the positions of the crystallographic water molecules were compared to peaks in the MD density after alignment of the protein in each snapshot. The results of this analysis indicate that all simulations reproduce the crystallographic water structure rather poorly. However, if we define the crystallographic water sites based on their distances to nearby protein atoms and follow these sites throughout the simulations, the MD simulations reproduce the crystallographic water sites much better. This shows that the failure of MD simulations to reproduce the water structure around proteins in crystal structures observed both in this and previous studies is caused by the problem of identifying water sites for a flexible and dynamic protein (traditionally done by overlaying the structures). Our local clustering approach solves the problem and shows that the MD simulations reasonably reproduce the water structure observed in crystals. Furthermore, analysis of the crystal MD simulations indicates a few water molecules that are close to unmodeled electron density peaks in the crystal structures, suggesting that crystal MD could be used as a complementary tool for identifying and modelling water in protein crystallography.

Highly stretchable, self-healing, and 3D printing prefabricatable hydrophobic association hydrogels with the assistance of electrostatic interaction

Self-healing and 3D printing prefabricatable physically crosslinked hydrogels were prepared by copolymerization of butyl acrylate, 2-(dimethylamino)ethyl methacrylate, and methacrylic acid, followed by soaking in water.