Water-mediated Interactions Research Articles

Interactions that influence the binding of synthetic heparan sulfate based disaccharides to fibroblast growth factor-2.

Heparan sulfate (HS) is a linear sulfated polysaccharide that mediates protein activities at the cell-extracellular interface. Its interactions with proteins depend on the complex patterns of sulfonations and sugar residues. Previously, we synthesized all 48 potential disaccharides found in HS and used them for affinity screening and X-ray structural analysis with fibroblast growth factor-1 (FGF1). Herein, we evaluated the affinities of the same sugars against FGF2 and determined the crystal structures of FGF2 in complex with three disaccharides carrying N-sulfonated glucosamine and 2-O-sulfonated iduronic acid as basic backbones. The crystal structures show that water molecules mediate different interactions between the 3-O-sulfonate group and Lys125. Moreover, the 6-O-sulfonate group forms intermolecular interactions with another FGF2 unit apart from the main binding site. These findings suggest that the water-mediated interactions and the intermolecular interactions influence the binding affinity of different disaccharides with FGF2, correlating with their respective dissociation constants in solution.

Glutamine synthetase mutations that cause glutamine deficiency: mechanistic insights from molecular dynamics simulations

Glutamine synthetase mutations that cause glutamine deficiency: mechanistic insights from molecular dynamics simulations

Polypharmacology rescored: Protein–ligand interaction profiles for remote binding site similarity assessment

Detection of remote binding site similarity in proteins plays an important role for drug repositioning and off-target effect prediction. Various non-covalent interactions such as hydrogen bonds and van-der-Waals forces drive ligands' molecular recognition by binding sites in proteins. The increasing amount of available structures of protein-small molecule complexes enabled the development of comparative approaches. Several methods have been developed to characterize and compare protein–ligand interaction patterns. Usually implemented as fingerprints, these are mainly used for post processing docking scores and (off-)target prediction. In the latter application, interaction profiles detect similarities in the bound interactions of different ligands and thus identify essential interactions between a protein and its small molecule ligands. Interaction pattern similarity correlates with binding site similarity and is thus contributing to a higher precision in binding site similarity assessment of proteins with distinct global structure. This renders it valuable for existing drug repositioning approaches in structural bioinformatics.Current methods to characterize and compare structure-based interaction patterns – both for protein-small-molecule and protein–protein interactions – as well as their potential in target prediction will be reviewed in this article. The question of how the set of interaction types, flexibility or water-mediated interactions, influence the comparison of interaction patterns will be discussed. Due to the wealth of protein–ligand structures available today, predicted targets can be ranked by comparing their ligand interaction pattern to patterns of the known target. Such knowledge-based methods offer high precision in comparison to methods comparing whole binding sites based on shape and amino acid physicochemical similarity.

Water-mediated ion-ion interactions are enhanced at the water vapor-liquid interface.

There is overwhelming evidence that ions are present near the vapor-liquid interface of aqueous salt solutions. Charged groups can also be driven to interfaces by attaching them to hydrophobic moieties. Despite their importance in many self-assembly phenomena, how ion-ion interactions are affected by interfaces is not understood. We use molecular simulations to show that the effective forces between small ions change character dramatically near the water vapor-liquid interface. Specifically, the water-mediated attraction between oppositely charged ions is enhanced relative to that in bulk water. Further, the repulsion between like-charged ions is weaker than that expected from a continuum dielectric description and can even become attractive as the ions are drawn to the vapor side. We show that thermodynamics of ion association are governed by a delicate balance of ion hydration, interfacial tension, and restriction of capillary fluctuations at the interface, leading to nonintuitive phenomena, such as water-mediated like charge attraction. "Sticky" electrostatic interactions may have important consequences on biomolecular structure, assembly, and aggregation at soft liquid interfaces. We demonstrate this by studying an interfacially active model peptide that changes its structure from α-helical to a hairpin-turn-like one in response to charging of its ends.

Structures of darunavir-resistant HIV-1 protease mutant reveal atypical binding of darunavir to wide open flaps.

The molecular basis for high resistance to clinical inhibitors of HIV-1 protease (PR) was examined for the variant designated PRP51 that was selected for resistance to darunavir (DRV). High resolution crystal structures of PRP51 with the active site D25N mutation revealed a ligand-free form and an inhibitor-bound form showing a unique binding site and orientation for DRV. This inactivating mutation is known to increase the dimer dissociation constant and decrease DRV affinity of PR. The PRP51-D25N dimers were in the open conformation with widely separated flaps, as reported for other highly resistant variants. PRP51-D25N dimer bound two DRV molecules and showed larger separation of 8.7 Å between the closest atoms of the two flaps compared with 4.4 Å for the ligand-free structure of this mutant. The ligand-free structure, however, lacked van der Waals contacts between Ile50 and Pro81′ from the other subunit in the dimer, unlike the majority of PR structures. DRV is bound inside the active site cavity; however, the inhibitor is oriented almost perpendicular to its typical position and exhibits only 2 direct hydrogen bond and two water-mediated interactions with atoms of PRP51-D25N compared with 11 hydrogen bond interactions seen for DRV bound in the typical position in wild-type enzyme. The atypical location of DRV may provide opportunities for design of novel inhibitors targeting the open conformation of PR drug-resistant mutants.

Structural Basis for PI(4)P-Specific Membrane Recruitment of the Legionella pneumophila Effector DrrA/SidM

Recruitment of the Legionella pneumophila effector DrrA to the Legionella-containing vacuole, where it activates and AMPylates Rab1, is mediated by a P4M domain that binds phosphatidylinositol 4-phosphate [PI(4)P] with high affinity and specificity. Despite the importance of PI(4)P in Golgi trafficking and its manipulation by pathogens, the structural bases for PI(4)P-dependent membrane recruitment remain poorly defined. Here, we determined the crystal structure of a DrrA fragment including the P4M domain in complex with dibutyl PI(4)P and investigated the determinants of phosphoinositide recognition and membrane targeting. Headgroup recognition involves an elaborate network of direct and water-mediated interactions with basic and polar residues in the context of a deep, constrictive binding pocket. An adjacent hydrophobic helical element packs against the acyl chains and inserts robustly into PI(4)P-containing monolayers. The structural, biochemical, and biophysical data reported here support a detailed structural mechanism for PI(4)P-dependent membrane targeting by DrrA.

Heterogeneous Hydration of p53/MDM2 Complex.

Water-mediated interactions play critical roles in biomolecular recognition processes. Explicit solvent molecular dynamics (MD) simulations and the variational implicit-solvent model (VISM) are used to study those hydration properties during binding for the biologically important p53/MDM2 complex. Unlike simple model solutes, in such a realistic and heterogeneous solute–solvent system with both geometrical and chemical complexity, the local water distribution sensitively depends on nearby amino acid properties and the geometric shape of the protein. We show that the VISM can accurately describe the locations of high and low density solvation shells identified by the MD simulations and can explain them by a local coupling balance of solvent–solute interaction potentials and curvature. In particular, capillary transitions between local dry and wet hydration states in the binding pocket are captured for interdomain distance between 4 to 6 Å, right at the onset of binding. The underlying physical connection between geometry and polarity is illustrated and quantified. Our study offers a microscopic and physical insight into the heterogeneous hydration behavior of the biologically highly relevant p53/MDM2 system and demonstrates the fundamental importance of hydrophobic effects for biological binding processes. We hope our study can help to establish new design rules for drugs and medical substances.

Direct and indirect interactions in the recognition between a cross‐neutralizing antibody and the four serotypes of dengue virus

Dengue fever is the most important vector-borne viral disease. Four serotypes of dengue virus, DENV1 to DENV4, coexist. Secondary infection by a different serotype is a risk factor for severe dengue. Monoclonal antibody mAb4E11 neutralizes the four serotypes of DENV with varying efficacies by recognizing an epitope located within domain-III (ED3) of the viral envelope (E) protein. To better understand the cross-reactivities between mAb4E11 and the four serotypes of DENV, we constructed mutations in both Fab4E11 fragment and ED3, and we searched for indirect interactions in the crystal structures of the four complexes. According to the serotype, 7 to 12 interactions are mediated by one water molecule, 1 to 10 by two water molecules, and several of these interactions are conserved between serotypes. Most interfacial water molecules make hydrogen bonds with both antibody and antigen. Some residues or atomic groups are engaged in both direct and water-mediated interactions. The doubly-indirect interactions are more numerous in the complex of lowest affinity. The third complementarity determining region of the light chain (L-CDR3) of mAb4E11 does not contact ED3. The structures and double-mutant thermodynamic cycles showed that the effects of (hyper)-mutations in L-CDR3 on affinity were caused by conformational changes and indirect interactions with ED3 through other CDRs. Exchanges of residues between ED3 serotypes showed that their effects on affinity were context dependent. Thus, conformational changes, structural context, and indirect interactions should be included when studying cross-reactivity between antibodies and different serotypes of viral antigens for a better design of diagnostics, vaccine, and therapeutic tools against DENV and other Flaviviruses.

Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations

Water is the natural medium of molecules in the cell and plays an important role in protein structure, function and interaction with small molecule ligands. However, the widely used molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method for binding energy calculation does not explicitly take account of water molecules that mediate key protein-ligand interactions. We have developed a protocol to include water molecules that mediate ligand-protein interactions as part of the protein structure in calculation of MM/PBSA binding energies (a method we refer to as water-MM/PBSA) for a series of JNK3 kinase inhibitors. Improved correlation between water-MM/PBSA binding energies and experimental IC50 values was obtained compared to that obtained from classical MM/PBSA binding energy. This improved correlation was further validated using sets of neuraminidase and avidin inhibitors. The observed improvement, however, appears to be limited to systems in which there are water-mediated ligand-protein hydrogen bond interactions. We conclude that the water-MM/PBSA method performs better than classical MM/PBSA in predicting binding affinities when water molecules play a direct role in mediating ligand-protein hydrogen bond interactions.

Bridging a Gap Between Cytochrome Bc1 Complex Structure and Function

Cytochrome (cyt) bc1 complex has two active sites where energy-sustaining electron/proton transfer reactions occur. To determine the exact transfer steps involving the substrate quinone, the cyt bc1 complex of Rhodobacter capsulatus was studied using atomistic molecular dynamics simulations. The computations give clear-cut ideas how the proton transfer might take place at both of the active sites. The simulations strongly suggest that water-mediated interactions are needed to achieve quinol oxidation and to prevent short-circuit suppression at the Qo-site. In the case of the Qi-site, the simulations emphasize the role of negatively charged cardiolipin in quinone reduction. Together with site-directed mutagenesis this in silico approach provides long-awaited information regarding the structure-function relationship of the cyt bc1 complex.

Antigen-Antibody Interaction Database (AgAbDb): a compendium of antigen-antibody interactions.

Antigen-Antibody Interaction Database (AgAbDb) is an immunoinformatics resource developed at the Bioinformatics Centre, University of Pune, and is available online at http://bioinfo.net.in/AgAbDb.htm. Antigen-antibody interactions are a special class of protein-protein interactions that are characterized by high affinity and strict specificity of antibodies towards their antigens. Several co-crystal structures of antigen-antibody complexes have been solved and are available in the Protein Data Bank (PDB). AgAbDb is a derived knowledgebase developed with an objective to compile, curate, and analyze determinants of interactions between the respective antigen-antibody molecules. AgAbDb lists not only the residues of binding sites of antigens and antibodies, but also interacting residue pairs. It also helps in the identification of interacting residues and buried residues that constitute antibody-binding sites of protein and peptide antigens. The Antigen-Antibody Interaction Finder (AAIF), a program developed in-house, is used to compile the molecular interactions, viz. van der Waals interactions, salt bridges, and hydrogen bonds. A module for curating water-mediated interactions has also been developed. In addition, various residue-level features, viz. accessible surface area, data on epitope segment, and secondary structural state of binding site residues, are also compiled. Apart from the PDB numbering, Wu-Kabat numbering and explicit definitions of complementarity-determining regions are provided for residues of antibodies. The molecular interactions can be visualized using the program Jmol. AgAbDb can be used as a benchmark dataset to validate algorithms for prediction of B-cell epitopes. It can as well be used to improve accuracy of existing algorithms and to design new algorithms. AgAbDb can also be used to design mimotopes representing antigens as well as aid in designing processes leading to humanization of antibodies.

Simulation optimization of spherical non-polar guest recognition by deep-cavity cavitands

Biomimetic deep-cavity cavitand hosts possess unique recognition and encapsulation properties that make them capable of selectively binding a range of non-polar guests within their hydrophobic pocket. Adamantane based derivatives which snuggly fit within the pocket of octa-acid deep cavity cavitands exhibit some of the strongest host binding. Here we explore the roles of guest size and attractiveness on optimizing guest binding to form 1:1 complexes with octa-acid cavitands in water. Specifically we simulate the water-mediated interactions of the cavitand with adamantane and a range of simple Lennard-Jones guests of varying diameter and attractive well-depth. Initial simulations performed with methane indicate hydrated methanes preferentially reside within the host pocket, although these guests frequently trade places with water and other methanes in bulk solution. The interaction strength of hydrophobic guests increases with increasing size from sizes slightly smaller than methane to Lennard-Jones guests comparable in size to adamantane. Over this guest size range the preferential guest binding location migrates from the bottom of the host pocket upwards. For guests larger than adamantane, however, binding becomes less favorable as the minimum in the potential-of-mean force shifts to the cavitand face around the portal. For a fixed guest diameter, the Lennard-Jones well-depth is found to systematically shift the guest-host potential-of-mean force to lower free energies, however, the optimal guest size is found to be insensitive to increasing well-depth. Ultimately our simulations show that adamantane lies within the optimal range of guest sizes with significant attractive interactions to match the most tightly bound Lennard-Jones guests studied.

Insight of Endo-1,4-Xylanase II from Trichoderma reesei: Conserved Water-Mediated H-Bond and Ion Pairs Interactions

Endo-1,4-Xylanase II is an enzyme which degrades the linear polysaccharide beta-1,4-xylan into xylose. This enzyme shows highest enzyme activity around 55°C, even without being stabilized by the disulphide bridges. A set of nine high resolution crystal structures of Xylanase II (1.11-1.80Å) from Trichoderma reesei were selected and analyzed in order to identify the invariant water molecules, ion pairs and water-mediated ionic interactions. The crystal structure (PDB-id: 2DFB) solved at highest resolution (1.11Å) was chosen as the reference and the remaining structures were treated as mobile molecules. These structures were then superimposed with the reference molecule to observe the invariant water molecules using 3-dimensional structural superposition server. A total of 37 water molecules were identified to be invariant molecules in all the crystal structures, of which 26 invariant molecules have hydrogen bond interactions with the back bone of residues and 21 invariant water molecules have interactions with side chain residues. The structural and functional roles of these water molecules and ion pairs have been discussed. The results show that the invariant water molecules and ion pairs may be involved in maintaining the structural architecture, dynamics and function of the Endo-1,4-Xylanase II.

Structure-guided design of potent diazobenzene inhibitors for the BET bromodomains.

BRD4, characterized by two acetyl-lysine binding bromodomains and an extra-terminal (ET) domain, is a key chromatin organizer that directs gene activation in chromatin through transcription factor recruitment, enhancer assembly, and pause release of the RNA polymerase II complex for transcription elongation. BRD4 has been recently validated as a new epigenetic drug target for cancer and inflammation. Our current knowledge of the functional differences of the two bromodomains of BRD4, however, is limited and is hindered by the lack of selective inhibitors. Here, we report our structure-guided development of diazobenzene-based small-molecule inhibitors for the BRD4 bromodomains that have over 90% sequence identity at the acetyl-lysine binding site. Our lead compound, MS436, through a set of water-mediated interactions, exhibits low nanomolar affinity (estimated Ki of 30-50 nM), with preference for the first bromodomain over the second. We demonstrated that MS436 effectively inhibits BRD4 activity in NF-κB-directed production of nitric oxide and proinflammatory cytokine interleukin-6 in murine macrophages. MS436 represents a new class of bromodomain inhibitors and will facilitate further investigation of the biological functions of the two bromodomains of BRD4 in gene expression.

Molecular docking of thiamine reveals similarity in binding properties between the prion protein and other thiamine-binding proteins

Prion-induced diseases are a global health concern. The lack of effective therapy and 100% mortality rates for such diseases have made the prion protein an important target for drug discovery. Previous NMR experimental work revealed that thiamine and its derivatives bind the prion protein in a pocket near the N-terminal loop of helix 1, and conserved intermolecular interactions were noted between thiamine and other thiamine-binding proteins. Furthermore, water-mediated interactions were observed in all of the X-ray crystallographic structures of thiamine-binding proteins, but were not observed in the thiamine-prion NMR study. To better understand the potential role of water in thiamine-prion binding, a docking study was employed using structural X-ray solvent. Before energy minimization, docked thiamine assumed a "V" shape similar to some of the known thiamine-dependent proteins. Following minimization with NMR-derived restraints, the "F" conformation was observed. Our findings confirmed that water is involved in ligand stabilization and phosphate group interaction. The resulting refined structure of thiamine bound to the prion protein allowed the 4-aminopyrimidine ring of thiamine to π-stack with Tyr150, and facilitated hydrogen bonding between Asp147 and the amino group of 4-aminopyrimidine. Investigation of the π-stacking interaction through mutation of the tyrosine residue further revealed its importance in ligand placement. The resulting refined structure is in good agreement with previous experimental restraints, and is consistent with the pharmacophore model of thiamine-binding proteins.

Probing the Binding of Syzygium-Derived α-Glucosidase Inhibitors with N- and C-Terminal Human Maltase Glucoamylase by Docking and Molecular Dynamics Simulation

Human maltase glucoamylase (MGAM) is a potent molecular target for controlling post prandial glucose surplus in type 2 diabetes. Binding of small molecules from Syzygium sp. with α-glucosidase inhibitory potential in MGAM has been investigated in silico. Our results suggest that myricetin was the most potent inhibitor with high binding affinity for both N- and C-terminals of MGAM. Molecular dynamics revealed that myricetin interacts in its stretched conformation through water-mediated interactions with C-terminal of MGAM and by normal hydrogen bonding with the N-terminal. W1369 of the extended 21 amino acid residue helical loop of C-terminal plays a major role in myricetin binding. Owing to its additional sugar sites, overall binding of small molecules favours C-terminal MGAM.

Unintended consequences? Water molecules at biological and crystallographic protein–protein interfaces

The importance of protein–protein interactions (PPIs) is becoming increasingly appreciated, as these interactions lie at the core of virtually every biological process. Small molecule modulators that target PPIs are under exploration as new therapies. One of the greatest obstacles faced in crystallographically determining the 3D structures of proteins is coaxing the proteins to form “artificial” PPIs that lead to uniform crystals suitable for X-ray diffraction. This work compares interactions formed naturally, i.e., “biological”, with those artificially formed under crystallization conditions or “non-biological”. In particular, a detailed analysis of water molecules at the interfaces of high-resolution (≤2.30Å) X-ray crystal structures of protein–protein complexes, where 140 are biological protein–protein complex structures and 112 include non-biological protein–protein interfaces, was carried out using modeling tools based on the HINT forcefield. Surprisingly few and relatively subtle differences were observed between the two types of interfaces: (i) non-biological interfaces are more polar than biological interfaces, yet there is better organized hydrogen bonding at the latter; (ii) biological associations rely more on water-mediated interactions with backbone atoms compared to non-biological associations; (iii) aromatic/planar residues play a larger role in biological associations with respect to water, and (iv) Lys has a particularly large role at non-biological interfaces. A support vector machines (SVMs) classifier using descriptors from this study was devised that was able to correctly classify 84% of the two interface types.

Variational Implicit-Solvent Modeling of Host–Guest Binding: A Case Study on Cucurbit[7]uril|

The synthetic host cucurbit[7]uril (CB[7]) binds aromatic guests or metal complexes with ultrahigh affinity compared with that typically displayed in protein–ligand binding. Due to its small size, CB[7] serves as an ideal receptor–ligand system for developing computational methods for molecular recognition. Here, we apply the recently developed variational implicit-solvent model (VISM), numerically evaluated by the level-set method, to study hydration effects in the high-affinity binding of the B2 bicyclo[2.2.2]octane derivative to CB[7]. For the unbound host, we find that the host cavity favors the hydrated state over the dry state due to electrostatic effects. For the guest binding, we find reasonable agreement to experimental binding affinities. Dissection of the individual VISM free-energy contributions shows that the major driving forces are water-mediated hydrophobic interactions and the intrinsic (vacuum) host–guest van der Waals interactions. These findings are in line with recent experiments and molecular dynamics simulations with explicit solvent. It is expected that the level-set VISM, with further refinement on the electrostatic descriptions, can efficiently predict molecular binding and recognition in a wide range of future applications.

Altered Backbone and Side-Chain Interactions Result in Route Heterogeneity during the Folding of Interleukin-1β (IL-1β)

Deletion of the β-bulge trigger-loop results in both a switch in the preferred folding route, from the functional loop packing folding route to barrel closure, as well as conversion of the agonist activity of IL-1β into antagonist activity. Conversely, circular permutations of IL-1β conserve the functional folding route as well as the agonist activity. These two extremes in the folding-functional interplay beg the question of whether mutations in IL-1β would result in changes in the populations of heterogeneous folding routes and the signaling activity. A series of topologically equivalent water-mediated β-strand bridging interactions within the pseudosymmetric β-trefoil fold of IL-1β highlight the backbone water interactions that stabilize the secondary and tertiary structure of the protein. Additionally, conserved aromatic residues lining the central cavity appear to be essential for both stability and folding. Here, we probe these protein backbone-water molecule and side chain-side chain interactions and the role they play in the folding mechanism of this geometrically stressed molecule. We used folding simulations with structure-based models, as well as a series of folding kinetic experiments to examine the effects of the F42W core mutation on the folding landscape of IL-1β. This mutation alters water-mediated backbone interactions essential for maintaining the trefoil fold. Our results clearly indicate that this perturbation in the primary structure alters a structural water interaction and consequently modulates the population of folding routes accessed during folding and signaling activity.

Strategic targeting of multiple water-mediated interactions: a concise and rational structure-based design approach to potent and selective MMP-13 inhibitors.

As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.