Water in cavity-ligand recognition.
Ontology highlight
ABSTRACT:
SUBMITTER: Baron R
PROVIDER: S-EPMC2933114 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A key factor influencing a drug's efficacy is its residence time in the binding pocket of the host protein. Using atomistic computer simulation to predict this residence time and the associated dissociation process is a desirable but extremely difficult task due to the long timescales involved. This gets further complicated by the presence of biophysical factors such as steric and solvation effects. In this work, we perform molecular dynamics (MD) simulations of the unbinding of a popular prototypical hydrophobic cavity-ligand system using a metadynamics-based approach that allows direct assessment of kinetic pathways and parameters. When constrained to move in an axial manner, the unbinding time is found to be on the order of 4,000 s. In accordance with previous studies, we find that the cavity must pass through a region of sharp wetting transition manifested by sudden and high fluctuations in solvent density. When we remove the steric constraints on ligand, the unbinding happens predominantly by an alternate pathway, where the unbinding becomes 20 times faster, and the sharp wetting transition instead becomes continuous. We validate the unbinding timescales from metadynamics through a Poisson analysis, and by comparison through detailed balance to binding timescale estimates from unbiased MD. This work demonstrates that enhanced sampling can be used to perform explicit solvent MD studies at timescales previously unattainable, to our knowledge, obtaining direct and reliable pictures of the underlying physiochemical factors including free energies and rate constants.
| S-EPMC4593114 | biostudies-literature
Project description:Simplified model binding sites allow one to isolate entangled terms in molecular energy functions. Here, we investigate the effects on ligand recognition of the introduction of a histidine into a hydrophobic cavity in lysozyme. We docked 656040 molecules and tested 26 highly and nine poorly ranked. Twenty-one highly ranked molecules bound and five were false positives, while three poorly ranked molecules were false negatives. In the 16 X-ray complexes now known, the docking predictions overlaid well with the crystallographic results. Although ligand enrichment was high, the false negatives, the false positives, and the inability to rank order illuminated weaknesses in our scoring, particularly overweighed apolar and underweighted polar terms. Adjusting these led to new problems, reflecting the entangled nature of docking scoring functions. Changes in ligand affinity relative to other lysozyme cavities speak to the subtleties of molecular recognition even in these simple sites and to their relevance for testing different models of recognition.
| S-EPMC3624796 | biostudies-literature
Project description:Targeting the minor groove of DNA through binding to a small molecule has long been considered an important molecular-recognition strategy in biology. A wide range of synthetic heterocyclic molecules bind noncovalently in the minor groove of the double helix and are also effective against a number of human and animal diseases. A classic structural concept, the isohelicity principle, has guided much of this work: such heterocyclic molecules require a shape that complements the convex surface of the minor groove. Researchers have used this principle to design molecules that can read DNA sequences. This principle also predicts that molecules that lack the complementary shape requirement would only bind weakly to DNA. Recently, however, researchers have unexpectedly found that some essentially linear compounds, which do not have this feature, can have high DNA affinity. In this Account, we discuss an alternative recognition concept based on these new findings. We demonstrate that highly structured water molecules can play a key role in mediating between the ligand and DNA minor groove without loss of binding affinity. Combined structural and thermodynamic approaches to understanding the behavior of these molecules have shown that there are different categories of bound water in their DNA complexes. For example, application of this water-bridging concept to the phenylamidine platform has resulted in the discovery of molecules with high levels of biological activity and low nonspecific toxicity. Some of these molecules are now in advanced clinical trials.
| S-EPMC2668807 | biostudies-literature
Project description:Water plays a crucial part in virtually all protein-ligand binding processes in and out of equilibrium. Here, we investigate the role of water in the binding kinetics of a ligand to a prototypical hydrophobic pocket by explicit-water molecular dynamics (MD) simulations and implicit diffusional approaches. The concave pocket in the unbound state exhibits wet/dry hydration oscillations whose magnitude and time scale are significantly amplified by the approaching ligand. In turn, the ligand's stochastic motion intimately couples to the slow hydration fluctuations, leading to a sixfold-enhanced friction in the vicinity of the pocket entrance. The increased friction considerably decelerates association in the otherwise barrierless system, indicating the importance of molecular-scale hydrodynamic effects in cavity-ligand binding arising due to capillary fluctuations. We derive and analyze the diffusivity profile and show that the mean first passage time distribution from the MD simulation can be accurately reproduced by a standard Brownian dynamics simulation if the appropriate position-dependent friction profile is included. However, long-time decays in the water-ligand (random) force autocorrelation demonstrate violation of the Markovian assumption, challenging standard diffusive approaches for rate prediction. Remarkably, the static friction profile derived from the force correlations strongly resembles the profile derived on the Markovian assumption apart from a simple shift in space, which can be rationalized by a time-space retardation in the ligand's downhill dynamics toward the pocket. The observed spatiotemporal hydrodynamic coupling may be of biological importance providing the time needed for conformational receptor-ligand adjustments, typical of the induced-fit paradigm.
| S-EPMC3557044 | biostudies-literature
Project description:Developing a ligand with high affinity for a specific protein target is essential for drug design, and water molecules are well known to play a key role in protein-drug recognition. However, predicting the role of particularly ordered water molecules in drug binding remains challenging. Furthermore, hydration free energy contributed from the water network, including the second shell of water molecules, is far from being well studied. In this research we focused on these aspects to accurately and efficiently evaluate water effects in protein-ligand binding affinity. We developed a new strategy using a free-energy calculation method, VM2. We successfully predicted the stable ordered water molecules in a number of protein systems: PDE 10a, HSP90, tryptophan synthase (TRPS), CDK2 and Factor Xa. In some of these, the second shell of water molecules appeared to be critical in protein-ligand binding. We also applied the strategy to largely improve binding free energy calculation using the MM/PBSA method. When applying MM/PBSA alone for two systems, CDK2 and Factor Xa, the computed binding free energy resulted in poor to moderate R2 values with experimental data. However, including water free energy correction greatly improved the free energy calculation. Furthermore, our work helped to explain how xk263 is a 1000 times faster binder to HIVp than ritonavir, a potentially useful tool for investigating binding kinetics. Our studies reveal the importance of fully considering water effects in therapeutic developments in pharmaceutical and biotechnology industries and for fundamental research in protein-ligand recognition.
| S-EPMC9970846 | biostudies-literature
Project description:Glutamate racemase (GR) is a cofactor independent amino acid racemase that has recently garnered increasing attention as an antimicrobial drug target. There are numerous high resolution crystal structures of GR, yet these are invariably bound to either D-glutamate or very weakly bound oxygen-based salts. Recent in silico screens have identified a number of new competitive inhibitor scaffolds, which are not based on D-Glu, but exploit many of the same hydrogen bond donor positions. In silico studies on 1-H-benzimidazole-2-sulfonic acid (BISA) show that the sulfonic acid points to the back of the GR active site, in the most buried region, analogous to the C2-carboxylate binding position in the GR-d-glutamate complex. Furthermore, BISA has been shown to be the strongest nonamino acid competitive inhibitor. Previously published computational studies have suggested that a portion of this binding strength is derived from complexation with a more closed active site, relative to weaker ligands, and in which the internal water network is more isolated from the bulk solvent. In order to validate key contacts between the buried sulfonate moiety of BISA and moieties in the back of the enzyme active site, as well as to probe the energetic importance of the potentially large number of interstitial waters contacted by the BISA scaffold, we have designed several mutants of Asn75. GR-N75A removes a key hydrogen bond donor to the sulfonate of BISA, but also serves to introduce an additional interstitial water, due to the newly created space of the mutation. GR- N75L should also show the loss of a hydrogen bond donor to the sulfonate of BISA, but does not (a priori) seem to permit an additional interstitial water contact. In order to investigate the dynamics, structure, and energies of this water-mediated complexation, we have employed the extended linear response (ELR) approach for the calculation of binding free energies to GR, using the YASARA2 knowledge based force field on a set of ten GR complexes, and yielding an R-squared value of 0.85 and a RMSE of 2.0 kJ/mol. Surprisingly, the inhibitor set produces a uniformly large interstitial water contribution to the electrostatic interaction energy (<V(el)>), ranging from 30 to >50%, except for the natural substrate (D-glutamate), which has only a 7% contribution of <V(el)> from water. The broader implications for predicting and exploiting significant interstitial water contacts in ligand-enzyme complexation are discussed.
| S-EPMC3782002 | biostudies-literature
Project description:Ligands usually bind to proteins by displacing water from the binding site. The affinity and kinetics of binding therefore depend on the hydration characteristics of the site. Here, we show that the extreme case of a completely dehydrated free binding site is realized for the large nonpolar binding cavity in bovine beta-lactoglobulin. Because spatially delocalized water molecules may escape detection by x-ray diffraction, we use water (17)O and (2)H magnetic relaxation dispersion (MRD), (13)C NMR spectroscopy, molecular dynamics simulations, and free energy calculations to establish the absence of water from the binding cavity. Whereas carbon nanotubes of the same diameter are filled by a hydrogen-bonded water chain, the MRD data show that the binding pore in the apo protein is either empty or contains water molecules with subnanosecond residence times. However, the latter possibility is ruled out by the computed hydration free energies, so we conclude that the 315 A(3) binding pore is completely empty. The apo protein is thus poised for efficient binding of fatty acids and other nonpolar ligands. The qualitatively different hydration of the beta-lactoglobulin pore and carbon nanotubes is caused by subtle differences in water-wall interactions and water entropy.
| S-EPMC2359778 | biostudies-literature
Project description:Large-scale and high-efficient water collection of microfibers with long-term durability still remains challenging. Here we present well-controlled, bioinspired spindle-knot microfibers with cavity knots (named cavity-microfiber), precisely fabricated via a simple gas-in-water microfluidic method, to address this challenge. The cavity-microfiber is endowed with unique surface roughness, mechanical strength, and long-term durability due to the design of cavity as well as polymer composition, thus enabling an outstanding performance of water collection. The maximum water volume collected on a single knot is almost 495 times than that of the knot on the cavity-microfiber. Moreover, the spider-web-like networks assembled controllably by cavity-microfibers demonstrate excellent large-scale and high-efficient water collection. To maximize the water-collecting capacity, nodes/intersections should be designed on the topology of the network as many as possible. Our light-weighted yet tough, low-cost microfibers with high efficiency in directional water transportation offers promising opportunities for large-scale water collection in water-deficient areas.
| S-EPMC5714965 | biostudies-literature
Project description:A prominent aqueous cavity is formed by the junction of three identical subunits in the excitatory amino acid transporter (EAAT) family. To investigate the effect of this structure on the interaction of ligands with the transporter, we recorded currents in voltage-clamped Xenopus oocytes expressing EAATs and used concentration jumps to measure binding and unbinding rates of a high-affinity aspartate analog that competitively blocks transport (?-2-fluorenyl-aspartylamide; 2-FAA). The binding rates of the blocker were approximately one order of magnitude slower than l-Glu and were not significantly different for EAAT1, EAAT2, or EAAT3, but 2-FAA exhibited higher affinity for the neuronal transporter EAAT3 as a result of a slower dissociation rate. Unexpectedly, the rate of recovery from block was increased by l-Glu in a saturable and concentration-dependent manner, ruling out a first-order mechanism and suggesting that following unbinding, there is a significant probability of ligand rebinding to the same or neighboring subunits within a trimer. Consistent with such a mechanism, coexpression of wild-type subunits with mutant (R447C) subunits that do not bind glutamate or 2-FAA also increased the unblocking rate. The data suggest that electrostatic and steric factors result in an effective dissociation rate that is approximately sevenfold slower than the microscopic subunit unbinding rate. The quaternary structure, which has been conserved through evolution, is expected to increase the transporters' capture efficiency by increasing the probability that following unbinding, a ligand will rebind as opposed to being lost to diffusion.
| S-EPMC3169134 | biostudies-literature
Project description:Internal water molecules are important to protein structure and function, but positional disorder and low occupancies can obscure their detection by X-ray crystallography. Here, we show that water can be detected within the distal cavities of myoglobin mutants by subtle changes in the absorbance spectrum of pentacoordinate heme, even when the presence of solvent is not readily observed in the corresponding crystal structures. A well-defined, noncoordinated water molecule hydrogen bonded to the distal histidine (His64) is seen within the distal heme pocket in the crystal structure of wild type (wt) deoxymyoglobin. Displacement of this water decreases the rate of ligand entry into wt Mb, and we have shown previously that the entry of this water is readily detected optically after laser photolysis of MbCO complexes. However, for L29F and V68L Mb no discrete positions for solvent molecules are seen in the electron density maps of the crystal structures even though His64 is still present and slow rates of ligand binding indicative of internal water are observed. In contrast, time-resolved perturbations of the visible absorption bands of L29F and V68L deoxyMb generated after laser photolysis detect the entry and significant occupancy of water within the distal pockets of these variants. Thus, the spectral perturbation of pentacoordinate heme offers a potentially robust system for measuring nonspecific hydration of the active sites of heme proteins.
| S-EPMC2774273 | biostudies-literature