How Can Hydrophobic Association Be Enthalpy Driven?

Hydrophobic residues of melittin mediate its binding to αA-crystallin

The molecular chaperone αA-crystallin, mainly localized in the human ocular lens, is believed to protect the lens from opacification and cataract, by suppressing the aggregation of the other lens proteins. The present study provides structural and thermodynamic insights into the ability of human αA-crystallin (HAA) to bind to its partially unfolded clients in the lens, using a small peptide, melittin from bee venom, as a model client. We characterized the thermodynamic parameters of the binding process between melittin and HAA through isothermal titration calorimetry (ITC), and found the binding to be endothermic and entropy-driven. We identified the amino acids in melittin important for binding to HAA by saturation-transfer difference (STD) nuclear magnetic resonance (NMR) experiments, and analysis of NMR line broadening upon titration of melittin with HAA. Our results suggest that hydrophobic residues Ile17 and Ile20 on the C-terminal region of melittin are in close contact with HAA in the melittin-HAA complex. Information obtained from NMR experiments was used to generate structural models of the melittin-HAA complex by molecular docking with high-ambiguity driven docking (HADDOCK). Structural models of the melittin-HAA complex reveal important principles underlying the interaction of HAA with its clients.

Identification of Potential Inhibitors of Cathepsin-B using Shape & Pharmacophore-based Virtual Screening, Molecular Docking and Explicit Water Thermodynamics

Lysosome has been long understood as a vital digestive organelle. Increasing reports indicate that the lysosome also plays a crucial role in the pathogenesis of a variety of neurodegenerative diseases, including Huntington's disease, Alzheimer's disease, Parkinson's disease and amyotrophic lateral sclerosis. Abnormal protein degradation and deposition stimulated by lysosomal dysfunction may cause age-related neurodegeneration. Enormous efforts have been devoted to the development of effective therapeutics against Alzheimer's disease, the most debilitating neurodegenerative disease. Endopeptidase activity of the Cathepsin-B is associated with the pathological processes. Work presented here focuses on identification of new inhibitors against Cathepsin-B protein using diverse computational approaches together. The inhibitors identified were further tested for in-vitro activity using enzyme based assay method. The identified inhibitors provided interesting understanding on how the water thermodynamic properties along with hydrophobic, steric, electronic, and structural requirements contribute to cathepsin-B inhibitory activity. These water thermodynamic studies, may further be used in computer aided drug discovery pipeline to design and predict more potent derivatives of various scaffolds as cathepsin-B inhibitors.

Role of Displacing Confined Solvent in the Conformational Equilibrium of β-Cyclodextrin

This study investigates the role of hydration and its relationship to the conformational equilibrium of the host molecule β-cyclodextrin. Molecular dynamics simulations indicate that the unbound β-cyclodextrin exhibits two state behavior in explicit solvent due to the opening and closing of its cavity. In implicit solvent, these transitions are not observed, and there is one dominant conformation of β-cyclodextrin with an open cavity. Based on these observations, we investigate the hypothesis that the expulsion of thermodynamically unfavorable water molecules into the bulk plays an important role in controlling the accessibility of the closed macrostate at room temperature. We compare the results of the molecular mechanics analytical generalized Born plus nonpolar solvation approach to those obtained through grid inhomogeneous solvation theory analysis with explicit solvation to elucidate the thermodynamic forces at play. The work illustrates the use of continuum solvent models to tease out solvation effects related to the inhomogeneity and the molecular nature of water and demonstrates the key role of the thermodynamics of enclosed hydration in driving the conformational equilibrium of molecules in solution.

Inclusion of enclosed hydration effects in the binding free energy estimation of dopamine D3 receptor complexes

Confined hydration and conformational flexibility are some of the challenges encountered for the rational design of selective antagonists of G-protein coupled receptors. We present a set of C3-substituted (-)-stepholidine derivatives as potent binders of the dopamine D3 receptor. The compounds are characterized biochemically, as well as by computer modeling using a novel molecular dynamics-based alchemical binding free energy approach which incorporates the effect of the displacement of enclosed water molecules from the binding site. The free energy of displacement of specific hydration sites is obtained using the Hydration Site Analysis method with explicit solvation. This work underscores the critical role of confined hydration and conformational reorganization in the molecular recognition mechanism of dopamine receptors and illustrates the potential of binding free energy models to represent these key phenomena.

Role of Water Hydrogen Bonding on Transport of Small Molecules inside Hydrophobic Channels

We conducted a systematic analysis of water networking inside smooth hyperboloid hydrophobic structures (cylindrical, barrel, and hourglass shapes) to elucidate the role of water hydrogen bonding on the transport of small hydrophobic molecules (ligands). Through a series of molecular dynamics simulations, we established that a hydrogen-bonded network forming along the centerline resulted in a water exclusion zone adjacent to the walls. The size of the exclusion zone is a function of the geometry and the nonbonded interaction strength, defining the effective hydrophobicity of the structure. Exclusion of water molecules from this zone results in lower apparent viscosity, leading to acceleration of ligand transport up to 7 times faster than that measured in the bulk. Transport of ligands into and out of the hydrophobic structures was shown to be controlled by a single water molecule that capped the narrow regions in the structure. This mechanism provides physical insights into the behavior and role of water in the bottleneck regions of hydrophobic enzyme channels. These findings were then used in a sister publication [ Escalante , D. E. , Comput. Struct. Biotechnol. J. 2019 17 757 760 ] to develop a model that can accurately predict the transport of ligands along nanochannels of broad-substrate specificity enzymes.

Biomolecular Solvation Structure Revealed by Molecular Dynamics Simulations

To compare ordered water positions from experiment with those from molecular dynamics (MD) simulations, a number of MD models of water structure in crystalline endoglucanase were calculated. The starting MD model was derived from a joint X-ray and neutron diffraction crystal structure, enabling the use of experimentally assigned protonation states. Simulations were performed in the crystalline state, using a periodic 2 × 2 × 2 supercell with explicit solvent. Water X-ray and neutron scattering density maps were computed from MD trajectories using standard macromolecular crystallography methods. In one set of simulations, harmonic restraints were applied to bias the protein structure toward the crystal structure. For these simulations, the recall of crystallographic waters using strong peaks in the MD water electron density was very good, and there also was substantial visual agreement between the boomerang-like wings of the neutron scattering density and the crystalline water hydrogen positions. An unrestrained simulation also was performed. For this simulation, the recall of crystallographic waters was much lower. For both restrained and unrestrained simulations, the strongest water density peaks were associated with crystallographic waters. The results demonstrate that it is now possible to recover crystallographic water structure using restrained MD simulations but that it is not yet reasonable to expect unrestrained MD simulations to do the same. Further development and generalization of MD water models for force-field development, macromolecular crystallography, and medicinal chemistry applications is now warranted. In particular, the combination of room-temperature crystallography, neutron diffraction, and crystalline MD simulations promises to substantially advance modeling of biomolecular solvation.

Enthalpy-Entropy Interplay in π-Stacking Interaction of Benzene Dimer in Water

Aromatic groups can engage in an interesting class of noncovalent interactions termed π-π interactions, which play a pivotal role in stabilizing a variety of molecular architectures, including nucleic acids, proteins, and supramolecular assemblies. When the aromatic compounds interact with each other in an aqueous environment, their association is facilitated by the hydrophobic effect-the trend of nonpolar solutes to aggregate in a polar solution. To develop an in-depth understanding of hydrophobic association, we investigate in the present work π-π interactions in water, employing as a paradigm the benzene dimer. Using DFT-CES, a mean-field QM/MM method recently developed by our group, we describe the benzene solute at a quantum-mechanical level. Full consideration of detailed solute-electron density enables an optimal description of the solute-solvent interactions, leading to an accurate prediction of hydration free energies. In π-stacking of benzene, we find an entropic stabilization associated with the shrinkage of the solvent-excluded volume, which agrees with the theory of hydrophobic effect at subnanoscales. However, at the equilibrium binding distance of the benzene dimer, we find that the entropic stabilization nearly cancels out due to the enthalpic cost required for dewetting the internal space. Such an enthalpy-entropy compensation leads the association free energy to be predominantly dictated by the solute-solute interaction enthalpy. The present work offers new insight into the mechanistic role of water and the primary thermodynamic driving force of hydrophobic association.

Transfer Free Energies and Partitioning of Small Molecules in Collapsed PNIPAM Polymers

A central quantity in the design of functional hydrogels used as nanocarrier systems, for instance, for drug delivery or adaptive nanocatalysis, is the partition ratio, which quantifies the uptake of a molecular substance by the polymer matrix. By employing all-atom molecular dynamics simulations, we study the solvation and partitioning (with respect to bulk water) of small subnanometer-sized solutes in a dense matrix of collapsed poly( N-isopropylacrylamide) (PNIPAM) polymers above the lower critical solution temperature in aqueous solution. We examine the roles of the solute's polarity and its size on the solubility properties in the thermoresponsive polymer. We show that the transfer free energies of nonpolar solutes from bulk water into the polymer are favorable and scale in a good approximation with the solute's surface area. Even for small solute size variation, partitioning can vary over orders of magnitude. A polar nature of the solute, on the other hand, generally opposes the transfer, at least for alkyl solutes. Finally, we find a strong correlation between the transfer free energies in the gel and the adsorption free energies on a single extended polymer chain, which enables us to relate the partition ratios in the swollen and collapsed state of a PNIPAM gel.

Molecular Dynamics Simulations of Water-Mediated Cholesterol Capture within an Open-Ended Single-Walled Carbon Nanotube

The excess concentration of cholesterol in the bloodstream can be brought down to a safer level by utilizing a potential cholesterol-binding agent such as a carbon nanotube (CNT). Here, we have probed solvent-mediated interactions between cholesterol and CNT by performing molecular dynamics simulations and potential-of-mean force (PMF) calculations. Simulations predict favorable interactions between water-mediated cholesterol and CNT owing to strong mutual interactions between them, whereas water plays an opposing role in the association. The breakdown of PMF into its enthalpic and entropic contributions indicates that contrary to traditional entropy-driven hydrophobic association, the cholesterol encapsulation within a CNT is primarily driven by enthalpy.

Water access and ligand dissociation at the binding site of proteins

Although water is undoubtedly an essential mediator of protein-ligand interactions, whether or not such water molecules are critical for the progress of ligand dissociation remains unclear. To gain a more complete understanding, molecular dynamics simulations are performed with two molecular systems, rigid model binding sites and trypsin-benzamidine. Free-energy landscapes are calculated with a suitably chosen solvent coordinate, which well describes water access to the ligand binding site. The results of free energy provided clear description of water-ligand exchange process, where two different mechanisms appear depending on whether the binding site is buried or not. As the site is more buried, water access is more difficult. When water does not access the site, ligand dissociation produces a large energy barrier, i.e., slow dissociation kinetics. This indicates that control of ligand dissociation kinetics becomes possible with burying the binding site. However, the results also showed that appropriate burying is important because burying reduces not only water access but also ligand binding. The role of the protein structural change is also discussed; it likely plays a similar role to water access because during ligand dissociation, it can make new coordination with the ligand binding site like water. These results contribute to the future pharmaceutical drug design and will be useful for fundamental exploration of various molecular events.

Heterogeneous Solvation in Distinctive Protein-Protein Interfaces Revealed by Molecular Dynamics Simulations

Water, despite being a driving force in biochemical processes, has an elusively complex microscopic behavior. While water can increase its local density near amphiphilic protein surfaces, water is also thought to evaporate from hydrophobic surfaces and cavities, an effect known as "dewetting". The existence and extent of dewetting effects remains elusive due to the difficulty in observing clear "drying" transitions in experiments or simulations. Here, we use explicit solvent molecular dynamics (MD) simulations to study the molecular solvation at the binding interfaces of two distinctive molecular complexes: the highly hydrophilic barnase-barstar and the highly hydrophobic MDM2-p53. Our simulations, in conjunction with simple volumetric analyses, reveal a strikingly different water behavior at the binding interfaces of these two molecular complexes. In both complexes, we observe significant changes in the water local density as the two proteins approach, supporting the existence of a clear dewetting transition in the case of MDM2-p53, with an onset distance of 5.6-7.6 Å. Furthermore, the solvation analysis reported herein is a valuable tool to capture and quantify persistent or transient dewetting events in future explicit solvent MD simulations.

How Ligand Protonation State Controls Water in Protein-Ligand Binding

The role of water in protein-ligand binding has been an intensely studied topic in recent years; however, how ligand protonation state change perturbs water has not been considered. Here we show that water dynamics and interactions can be controlled by the protonation state of ligand using continuous constant pH molecular dynamics simulations of two closely related model systems, β-secretase 1 and 2 (BACE1 and BACE2), in complex with a small-molecule inhibitor. Simulations revealed that, upon binding, the inhibitor pyrimidine ring remains deprotonated in BACE1 but becomes protonated in BACE2. Pyrimidine protonation results in water displacement, rigidification of the binding pocket, and shift in the ligand binding mode from water-mediated to direct hydrogen bonding. These findings not only support but also rationalize the most recent structure-selectivity data in BACE1 drug design. Binding-induced protonation state changes are likely common; our work offers a glimpse at how modeling protein-ligand binding while allowing ligand titration can further advance the understanding of water and structure-based drug design.

Identification of a Water-Coordinating HER2 Inhibitor by Virtual Screening Using Similarity-Based Scoring

Human epidermal growth factor receptor 2 (HER2) is a validated breast cancer drug target for small molecule inhibitors that target the ATP-binding pocket of the kinase domain. In this work, a large-scale virtual screen was performed to a novel homology model of HER2, in a hypothesized “fully active” state, that considered water-mediated interactions during the prioritization of compounds for experimental testing. This screen led to the identification of a new inhibitor with micro molar affinity and potency (K_d = 7.0 μM, IC₅₀ = 4.6 μM). Accompanying molecular dynamics simulations showed that inhibitor binding likely involves water coordination through an important water-mediated network previously identified in our laboratory. The predicted binding geometry also showed a remarkable overlap with the crystallographic poses for two previously reported inhibitors of the related Chk1 kinase. Concurrent with the HER2 studies, we developed formalized computational protocols that leverage solvated footprints (per-residue interaction maps that include bridging waters) to identify ligands that can “coordinate” or “displace” key binding site waters. Proof-of-concept screens targeting HIVPR and PARP1 demonstrate that molecules with high footprint overlap can be effectively identified in terms of their coordination or displacement patterns relative to a known reference. Overall, the procedures developed as a result of this study should be useful for researchers targeting HER2 and, more generally, for any protein in which the identification of compounds that exploit binding site waters is desirable.

Charge and hydration structure of dendritic polyelectrolytes: molecular simulations of polyglycerol sulphate

Macromolecules based on dendritic or hyperbranched polyelectrolytes have been emerging as high potential candidates for biomedical applications. Here we study the charge and solvation structure of dendritic polyglycerol sulphate (dPGS) of generations 0 to 3 in aqueous sodium chloride solution by explicit-solvent molecular dynamics computer simulations. We characterize dPGS by calculating several important properties such as relevant dPGS radii, molecular distributions, the solvent accessible surface area, and the partial molecular volume. In particular, as the dPGS exhibits high charge renormalization effects, we address the challenges of how to obtain a well-defined effective charge and surface potential of dPGS for practical applications. We compare implicit- and explicit-solvent approaches in our all-atom simulations with the coarse-grained simulations from our previous work. We find consistent values for the effective electrostatic size (i.e., the location of the effective charge of a Debye-Hückel sphere) within all the approaches, deviating at most by the size of a water molecule. Finally, the excess chemical potential of water insertion into dPGS and its thermodynamic signature are presented and rationalized.

Dimer-monomer equilibrium of human HSP27 is influenced by the in-cell macromolecular crowding environment and is controlled by fatty acids and heat

Small heat shock protein 27 (HSP27) is an essential element of the proteostasis network in human cells. The HSP27 monomer coexists with the dimer, which can bind unfolded client proteins. Here, we evaluated the in-cell dimer-monomer equilibrium and its relevance to the binding of client proteins in a normal human vascular endothelial cell line. When cells were treated with a membrane-permeable crosslinker, the protein existed primarily as a free monomer (27 kDa) with a markedly smaller percentage of dimer (54 kDa), hetero-conjugates, and minor smear-like bands. When the protein was crosslinked in a cell-free lysate, two of the hetero-conjugates that were crosslinked in live cells were also detected, but the dimer and other complexes were absent. However, when cells were pretreated with fatty acid (FA) and/or heat (42.5 °C), dissociation of the dimer was selectively prevented and two types of covalently linked dimers were increased. These changes occurred most prominently in cells treated with docosahexaenoic acid (DHA) and heat, which appeared to intensify the heat resistance of the cell. Both the formation of covalently linked dimers and heat resistance were prevented by N-acetylcysteine. By contrast, nearly all of the free monomers in the lysate converted to disulfide bond-linked dimers by a simple, long incubation at 4 °C. These results strongly suggest that the monomer-dimer equilibrium of HSP27 was inversed between the in-cell and cell-free systems. Temperature- and amphiphile-regulated dimerization was restricted probably due to the low hydration of the in-cell crowding environment.

Predicting Displaceable Water Sites Using Mixed-Solvent Molecular Dynamics

Water molecules are an important factor in protein-ligand binding. Upon binding of a ligand with a protein's surface, waters can either be displaced by the ligand or may be conserved and possibly bridge interactions between the protein and ligand. Depending on the specific interactions made by the ligand, displacing waters can yield a gain in binding affinity. The extent to which binding affinity may increase is difficult to predict, as the favorable displacement of a water molecule is dependent on the site-specific interactions made by the water and the potential ligand. Several methods have been developed to predict the location of water sites on a protein's surface, but the majority of methods are not able to take into account both protein dynamics and the interactions made by specific functional groups. Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique that explicitly accounts for the interaction of both water and small molecule probes with a protein's surface, allowing for their direct competition. This method has previously been shown to identify both active and allosteric sites on a protein's surface. Using a test set of eight systems, we have developed a method using MixMD to identify conserved and displaceable water sites. Conserved sites can be determined by an occupancy-based metric to identify sites which are consistently occupied by water even in the presence of probe molecules. Conversely, displaceable water sites can be found by considering the sites which preferentially bind probe molecules. Furthermore, the inclusion of six probe types allows the MixMD method to predict which functional groups are capable of displacing which water sites. The MixMD method consistently identifies sites which are likely to be nondisplaceable and predicts the favorable displacement of water sites that are known to be displaced upon ligand binding.

Tailoring the Variational Implicit Solvent Method for New Challenges: Biomolecular Recognition and Assembly

Predicting solvation free energies and describing the complex water behavior that plays an important role in essentially all biological processes is a major challenge from the computational standpoint. While an atomistic, explicit description of the solvent can turn out to be too expensive in large biomolecular systems, most implicit solvent methods fail to capture “dewetting” effects and heterogeneous hydration by relying on a pre-established (i.e., guessed) solvation interface. Here we focus on the Variational Implicit Solvent Method, an implicit solvent method that adds water “plasticity” back to the picture by formulating the solvation free energy as a functional of all possible solvation interfaces. We survey VISM's applications to the problem of molecular recognition and report some of the most recent efforts to tailor VISM for more challenging scenarios, with the ultimate goal of including thermal fluctuations into the framework. The advances reported herein pave the way to make VISM a uniquely successful approach to characterize complex solvation properties in the recognition and binding of large-scale biomolecular complexes.

Binding Thermodynamics and Kinetics Calculations Using Chemical Host and Guest: A Comprehensive Picture of Molecular Recognition

Understanding the fine balance between changes of entropy and enthalpy and the competition between a guest and water molecules in molecular binding is crucial in fundamental studies and practical applications. Experiments provide measurements. However, illustrating the binding/unbinding processes gives a complete picture of molecular recognition not directly available from experiments, and computational methods bridge the gaps. Here, we investigated guest association/dissociation with β-cyclodextrin (β-CD) by using microsecond-time-scale molecular dynamics (MD) simulations, postanalysis and numerical calculations. We computed association and dissociation rate constants, enthalpy, and solvent and solute entropy of binding. All the computed values of kon, koff, ΔH, ΔS, and ΔG using GAFF-CD and q4MD-CD force fields for β-CD could be compared with experimental data directly and agreed reasonably with experiment findings. In addition, our study further interprets experiments. Both force fields resulted in similar computed ΔG from independently computed kinetics rates, ΔG = -RT ln(kon·C0/koff), and thermodynamics properties, ΔG = ΔH - TΔS. The water entropy calculations show that the entropy gain of desolvating water molecules are a major driving force, and both force fields have the same strength of nonpolar attractions between solutes and β-CD as well. Water molecules play a crucial role in guest binding to β-CD. However, collective water/β-CD motions could contribute to different computed kon and ΔH values by different force fields, mainly because the parameters of β-CD provide different motions of β-CD, hydrogen-bond networks of water molecules in the cavity of free β-CD, and strength of desolvation penalty. As a result, q4MD-CD suggests that guest binding is mostly driven by enthalpy, while GAFF-CD shows that gaining entropy is the major driving force of binding. The study deepens our understanding of ligand-receptor recognition and suggests strategies for force field parametrization for accurately modeling molecular systems.

Enthalpy-driven micellization of oligocarbonate-fluorene end-functionalized Poly(ethylene glycol)☆

A fluorescent pyrene probe method was applied to measure the critical micelle concentration (CMC) of oligocarbonate-fluorene end-functionalized poly(ethylene glycol) (FmE445Fm) triblock copolymers in water. The CMC decreases with lower temperature and higher values of the hydrophobic block length, m. When analyzed by a closed-assembly micelle model, the estimated energetic parameters find a negative ΔH°mic and small positive ΔS°mic suggestive of enthalpy-driven micellization, which differs from entropy-driven oxyethylene/oxybutylene triblock copolymers and octaethylene glycol-n-alkyl ethers. The enthalpy-driven micellization of FmE445Fm may result from the limited hydration of individual hydrophobic F blocks that leads to few hydrogen-bonded waters released during F block association. The π-π stacking oligocarbonate-fluorene system also observed enthalpy-entropy compensation when compared to a series of published data on diblock and triblock copolymer systems. An anomalously low partition equilibrium constant for m = 15.3 implies a tightly-packed core that excludes pyrene intercalation into the fluorene core. This is discussed along with the possible limited applicability to estimate the CMC and potential model drug molecule insertions into the intercalated micelle core.

Molecular Shape and the Hydrophobic Effect

This review focuses on papers published since 2000 on the topic of the properties of solutes in water. More specifically, it evaluates the state of the art of our understanding of the complex relationship between the shape of a hydrophobe and the hydrophobic effect. To highlight this, we present a selection of references covering both empirical and molecular dynamics studies of small (molecular-scale) solutes. These include empirical studies of small molecules, synthetic hosts, crystalline monolayers, and proteins, as well as in silico investigations of entities such as idealized hard and soft spheres, small solutes, hydrophobic plates, artificial concavity, molecular hosts, carbon nanotubes and spheres, and proteins.

Functional details of the Mycobacterium tuberculosis VapBC26 toxin-antitoxin system based on a structural study: insights into unique binding and antibiotic peptides

Toxin-antitoxin (TA) systems are essential for bacterial persistence under stressful conditions. In particular, Mycobacterium tuberculosis express VapBC TA genes that encode the stable VapC toxin and the labile VapB antitoxin. Under normal conditions, these proteins interact to form a non-toxic TA complex, but the toxin is activated by release from the antitoxin in response to unfavorable conditions. Here, we present the crystal structure of the M. tuberculosis VapBC26 complex and show that the VapC26 toxin contains a pilus retraction protein (PilT) N-terminal (PIN) domain that is essential for ribonuclease activity and that, the VapB26 antitoxin folds into a ribbon-helix-helix DNA-binding motif at the N-terminus. The active site of VapC26 is sterically blocked by the flexible C-terminal region of VapB26. The C-terminal region of free VapB26 adopts an unfolded conformation but forms a helix upon binding to VapC26. The results of RNase activity assays show that Mg²⁺ and Mn²⁺ are essential for the ribonuclease activity of VapC26. As shown in the nuclear magnetic resonance spectra, several residues of VapB26 participate in the specific binding to the promoter region of the VapBC26 operon. In addition, toxin-mimicking peptides were designed that inhibit TA complex formation and thereby increase toxin activity, providing a novel approach to the development of new antibiotics.

Rational design of drug-like compounds targeting Mycobacterium marinum MelF protein

The mycobacterial mel2 locus (mycobacterial enhanced infection locus, Rv1936-1941) is Mycobacterium marinum and M. tuberculosis specific, which can withstand reactive oxygen species (ROS) and reactive nitrogen species (RNS) induced stress. A library of over a million compounds was screened using in silico virtual ligand screening (VLS) to identify inhibitors against the modeled structure of MelF protein expressed by melF of mel2 locus so that M. marinum’s ability to withstand ROS/RNS stress could be reduced. The top ranked 1000 compounds were further screened to identify 178 compounds to maximize the scaffold diversity by manually evaluating the interaction of each compound with the target site. M. marinum melF was cloned, expressed and purified as maltose binding protein (MBP)-tagged recombinant protein in Escherichia coli. After establishing the flavin dependent oxidoreductase activity of MelF (~ 84 kDa), the inhibitors were screened for the inhibition of enzyme activity of whole cell lysate (WCL) and the purified MelF. Amongst these, 16 compounds could significantly inhibit the enzyme activity of purified MelF. For the six best inhibitory compounds, the minimal inhibitory concentration (MIC) was determined to be 3.4–19.4 μM and 13.5–38.8 μM for M. marinum and M. tuberculosis, respectively. Similarly, the minimal bactericidal concentration (MBC) was determined to be 6.8–38.8 μM and 27–38.8 μM against M. marinum and M. tuberculosis, respectively. One compound each in combination with isoniazid (INH) also showed synergistic inhibitory effect against M. marinum and M. tuberculosis with no cytotoxicity in HeLa cells. Interestingly, these inhibitors did not display any non-specific protein-structure destabilizing effect. Such inhibitors targeting the anti-ROS/RNS machinery may facilitate the efficient killing of replicating and nonreplicating mycobacteria inside the host cells.

Molecular Dynamics Studies of Overbased Detergents on a Water Surface

Molecular dynamics (MD) simulations are reported of model overbased detergent nanoparticles on a model water surface which mimic their behavior on a Langmuir trough or large water droplet in engine oil. The simulations predict that the structure of the nanoparticle on a water surface is different to when it is immersed in a bulk hydrophobic solvent. The surfactant tails are partly directed out of the water, while the carbonate core maximizes its extent of contact with the water. Umbrella sampling calculations of the potential of mean force between two particles showed that they are associated with varying degrees with a maximum binding free energy of ca. 10 k_BT for the salicylate stabilized particle, ca. 8 k_BT for a sulfurized alkyl phenate stabilized particle, and ca. 5 k_BT for a sulfonate stabilized particle. The differences in the strength of attraction depend on the proximity of nearest approach and the energy penalty associated with the disruption of the hydration shell of water molecules around the calcium carbonate core when the two particles approach. This is greatest for the sulfonate particle, which partially loses the surfactant ions to the solution, and least for the salicylate, which forms the weakest water "cage". The particles are separated by a water hydration layer, even at the point of closest approach.

Water-Restructuring Mutations Can Reverse the Thermodynamic Signature of Ligand Binding to Human Carbonic Anhydrase

This study uses mutants of human carbonic anhydrase (HCAII) to examine how changes in the organization of water within a binding pocket can alter the thermodynamics of protein-ligand association. Results from calorimetric, crystallographic, and theoretical analyses suggest that most mutations strengthen networks of water-mediated hydrogen bonds and reduce binding affinity by increasing the enthalpic cost and, to a lesser extent, the entropic benefit of rearranging those networks during binding. The organization of water within a binding pocket can thus determine whether the hydrophobic interactions in which it engages are enthalpy-driven or entropy-driven. Our findings highlight a possible asymmetry in protein-ligand association by suggesting that, within the confines of the binding pocket of HCAII, binding events associated with enthalpically favorable rearrangements of water are stronger than those associated with entropically favorable ones.

Interaction of human serum albumin with uremic toxins: a thermodynamic study

Interaction of uremic toxins with HSA is studied by ITC and understood in terms of thermodynamic driving forces.

Solvent effects on ligand binding to a serine protease

ITC experiments and MD simulations reveal the mechanism behind enthalpy/entropy compensation upon trypsin-benzamidine binding at different solvation conditions.

Unified elucidation of the entropy-driven and -opposed hydrophobic effects

The association of nonpolar solutes is generally believed to be entropy driven, which has been shown to be true for the contact of small molecules, ellipsoids, and plates.

A combined treatment of hydration and dynamical effects for the modeling of host-guest binding thermodynamics: the SAMPL5 blinded challenge

As part of the SAMPL5 blinded experiment, we computed the absolute binding free energies of 22 host-guest complexes employing a novel approach based on the BEDAM single-decoupling alchemical free energy protocol with parallel replica exchange conformational sampling and the AGBNP2 implicit solvation model specifically customized to treat the effect of water displacement as modeled by the Hydration Site Analysis method with explicit solvation. Initial predictions were affected by the lack of treatment of ionic charge screening, which is very significant for these highly charged hosts, and resulted in poor relative ranking of negatively versus positively charged guests. Binding free energies obtained with Debye-Hückel treatment of salt effects were in good agreement with experimental measurements. Water displacement effects contributed favorably and very significantly to the observed binding affinities; without it, the modeling predictions would have grossly underestimated binding. The work validates the implicit/explicit solvation approach employed here and it shows that comprehensive physical models can be effective at predicting binding affinities of molecular complexes requiring accurate treatment of conformational dynamics and hydration.

Kinetics and equilibrium partitioning of dissolved BTEX in PDMS and POM sheets

Passive sampling of volatile organic chemicals from soil and groundwater is primarily important in assessing the status of environmental contamination. A group of low molecular weight pollutants usually found in petroleum fuels, benzene, toluene, ethylbenzene, and xylenes (BTEX) was studied for its kinetics and equilibrium partitioning with single-phase passive samplers using polydimethylsiloxane (PDMS) and polyoxymethylene (POM) as sorbing phase. PDMS (1 mm) and POM (0.076 mm) sheets were used for sorption of BTEX and concentrations were analyzed using GC-FID. The equilibrium absorption and desorption of PDMS in water was achieved after 120 min while POM sheets absorbed up to 35 days and desorbed in 7 days. The kinetic rate constants in PDMS is higher than in POM up to 3 orders of magnitude. Logarithms of partition coefficient were determined to be in the range of 1.6-2.8 for PDMS and 2.1-3.1 for POM. The results indicate that POM is a stronger sorbent for BTEX and has slower equilibration time than PDMS. The partitioning process for both polymers was found to be enthalpy-driven by measurement of K d values at varying temperatures. K d values increase at low temperature and high ionic strength conditions. Presence of other gasoline components, as well as dissolved organic matter, did not significantly affect equilibrium partitioning. A good 1:1 correlation between the measured and the predicted concentrations was established on testing the potential application of the constructed PDMS sampler on natural soils and artificial soils spiked with gasoline-contaminated water.

Recent developments in the theoretical, simulational, and experimental studies of the role of water hydrogen bonding in hydrophobic phenomena

Hydrophobic effects (hydrophobic hydration and hydrophobic interaction) constitute an important element of a wide variety of phenomena relevant to biological, physical, chemical, environmental, engineering, and pharmaceutical sciences, such as the immiscibility of oil and water, self-assembly of amphiphiles leading to micelle and membrane formation, folding and stability and unfolding of the native structure of a biologically active protein, gating of ion channels, wetting, froth floatation, and adhesion. On the other hand, the hydrogen bonding ability of water plays a major (if not crucial) role in hydrophobic phenomena. We present a review of most important and relatively recent experimental, simulational, and theoretical research on hydrophobic phenomena in various systems. With a particular interest we survey investigations clarifying the role of water hydrogen bonding therein, because it has been the main object of our own recent research. We have developed a probabilistic hydrogen bond (PHB) model that allows one to obtain an analytic expression for the number of bonds per water molecule as a function of its distance to a hydrophobe, hydrophobe radius, and temperature. Knowing that function, one can explicitly identify a water hydrogen bond contribution to the external potential whereto a water molecule is subjected near a hydrophobe. Combining the PHB model with the classical density functional theory (DFT), one can examine the contribution of water hydrogen bonding to the temperature and lengthscale effects on the hydration of particles and on their solvent-mediated interactions over the entire low-to-high temperature and small-to-large lengthscale ranges. We applied the combined DFT/PHB model to study a variety of hydrophobic phenomena such as (liquid) water in contact with a hydrophobic plate, solvation of spherical solutes of various radii in associated and non-associated liquids at various temperatures, the solvent-mediated interaction of spherical solutes and its temperature dependence, interaction of C60 fullerenes in water, temperature effect on the evaporation lengthscale of water confined between two hydrophobes, temperature dependence of the effective width of the solute-solvent transition layer and average density therein. These applications demonstrated that the DFT/PHB model can serve as a valuable tool in studying hydrophobic phenomena because it constitutes a balanced combination of simplicity, accuracy, and detail. The predictions of the combined DFT/PHB approach for the solvent density profiles and thermodynamic aspects of hydrophobic phenomena are generally in good agreement with experiments and simulations. For example, it predicts the small-to-large crossover lengthscale of its mechanism to be approximately in the range from 1nm to 4nm, and decreasing with increasing temperature. It also suggests that, in terms of the average fluid density in the solute-solvent transition layer, the transition layer for small hydrophobes (of radii ≲2 nm) becomes enriched with rather than depleted of fluid when both the solvent-solute affinity and hb-energy alteration ratio become large enough. The boundary values of these parameters, needed for the depletion-to-enrichment crossover, are predicted to decrease with increasing temperature.

Enthalpic Breakdown of Water Structure on Protein Active-Site Surfaces

The principles underlying water reorganization around simple nonpolar solutes are well understood and provide the framework for the classical hydrophobic effect, whereby water molecules structure themselves around solutes so that they maintain favorable energetic contacts with both the solute and the other water molecules. However, for certain solute surface topographies, water molecules, due to their geometry and size, are unable to simultaneously maintain favorable energetic contacts with both the surface and neighboring water molecules. In this study, we analyze the solvation of ligand-binding sites for six structurally diverse proteins using hydration site analysis and measures of local water structure, in order to identify surfaces at which water molecules are unable to structure themselves in a way that maintains favorable enthalpy relative to bulk water. These surfaces are characterized by a high degree of enclosure, weak solute-water interactions, and surface constraints that induce unfavorable pair interactions between neighboring water molecules. Additionally, we find that the solvation of charged side chains in an active site generally results in favorable enthalpy but can also lead to pair interactions between neighboring water molecules that are significantly unfavorable relative to bulk water. We find that frustrated local structure can occur not only in apolar and weakly polar pockets, where overall enthalpy tends to be unfavorable, but also in charged pockets, where overall water enthalpy tends to be favorable. The characterization of local water structure in these terms may prove useful for evaluating the displacement of water from diverse protein active-site environments.

The effect of desolvation on the binding of inhibitors to HIV-1 protease and cyclin-dependent kinases: Causes of resistance

Studies of the cyclin-dependent kinase inhibitors and HIV-1 protease inhibitors have confirmed that ligand-protein binding is dependent on desolvation effects. It has been found that a four parameter linear model incorporating desolvation energy, lipophilicity, dipole moment and molecular volume of the ligands is a good model to describe the binding between ligands and kinases or proteases. The resistance shown by MDR proteases to the anti-viral drugs is multi-faceted involving varying changes in desolvation, lipophilicity and dipole moment interaction compared to the non-resistant protease. Desolvation has been shown to be the dominant factor influencing the effect of inhibitors against the cyclin-dependent kinases, but lipophilicity and dipole moment are also significant factors. The model can differentiate between the inhibitory activity of CDK2/cycE, CDK1/cycB and CDK4/cycD enzymes.

Force fields for simulating the interaction of surfaces with biological molecules

The interaction of biomolecules with solid interfaces is of fundamental importance to several emerging biotechnologies such as medical implants, anti-fouling coatings and novel diagnostic devices. Many of these technologies rely on the binding of peptides to a solid surface, but a full understanding of the mechanism of binding, as well as the effect on the conformation of adsorbed peptides, is beyond the resolution of current experimental techniques. Nanoscale simulations using molecular mechanics offer potential insights into these processes. However, most models at this scale have been developed for aqueous peptide and protein simulation, and there are no proven models for describing biointerfaces. In this review, we detail the current research towards developing a non-polarizable molecular model for peptide-surface interactions, with a particular focus on fitting the model parameters as well as validation by choice of appropriate experimental data.

Spatial Decomposition of Translational Water-Water Correlation Entropy in Binding Pockets

A number of computational tools available today compute the thermodynamic properties of water at surfaces and in binding pockets by using inhomogeneous solvation theory (IST) to analyze explicit-solvent simulations. Such methods enable qualitative spatial mappings of both energy and entropy around a solute of interest and can also be applied quantitatively. However, the entropy estimates of existing methods have, to date, been almost entirely limited to the first-order terms in the IST’s entropy expansion. These first-order terms account for localization and orientation of water molecules in the field of the solute but not for the modification of water–water correlations by the solute. Here, we present an extension of the Grid Inhomogeneous Solvation Theory (GIST) approach which accounts for water–water translational correlations. The method involves rewriting the two-point density of water in terms of a conditional density and utilizes the efficient nearest-neighbor entropy estimation approach. Spatial maps of this second order term, for water in and around the synthetic host cucurbit[7]uril and in the binding pocket of the enzyme Factor Xa, reveal mainly negative contributions, indicating solute-induced water–water correlations relative to bulk water; particularly strong signals are obtained for sites at the entrances of cavities or pockets. This second-order term thus enters with the same, negative, sign as the first order translational and orientational terms. Numerical and convergence properties of the methodology are examined.