Search

Water molecules decrease the potential of mean force of a hydrogen bond (H-bond), as well as modulate (de)solvation forces, but exactly how much has not been easy to determine. Crystallographic water molecules provide snapshots of optimal solutions for the role of solvent in protein interactions, information that is often ignored by implicit solvent models. Motivated by high-resolution crystal structures, we describe a simple quantitative approach to explicitly incorporate the role of molecular water in protein interactions. Applications to protein-DNA interactions show that the accuracy of binding free-energy estimates improves significantly if a distinction is made between H-bonds that are desolvated (or only contact crystal waters), solvated by mobile waters trapped at the binding interface, or partially solvated through connections to bulk water. These different environments are modeled by a unique "water" scaling factor that decreases or increases the strength of hydrogen bonds depending on whether water contacts the acceptor or donor atoms or the bond is fully desolvated, respectively. Our empirical energies are fully consistent with mobile water molecules having a strong polarization effect in direct intermolecular interactions., (© 2010 Wiley-Liss, Inc.) more...

The dynamics of adenylate kinase of Escherichia coli (AKeco) and its complex with the inhibitor AP(5)A, are characterized by correlating the theoretical results obtained with the Gaussian Network Model (GNM) and the anisotropic network model (ANM) with the order parameters and correlation times obtained with Slowly Relaxing Local Structure (SRLS) analysis of (15)N-NMR relaxation data. The AMPbd and LID domains of AKeco execute in solution large amplitude motions associated with the catalytic reaction Mg(+2)*ATP + AMP --> Mg(+2)*ADP + ADP. Two sets of correlation times and order parameters were determined by NMR/SRLS for AKeco, attributed to slow (nanoseconds) motions with correlation time tau( perpendicular) and low order parameters, and fast (picoseconds) motions with correlation time tau( parallel) and high order parameters. The structural connotation of these patterns is examined herein by subjecting AKeco and AKeco*AP(5)A to GNM analysis, which yields the dynamic spectrum in terms of slow and fast modes. The low/high NMR order parameters correlate with the slow/fast modes of the backbone elucidated with GNM. Likewise, tau( parallel) and tau( perpendicular) are associated with fast and slow GNM modes, respectively. Catalysis-related domain motion of AMPbd and LID in AKeco, occurring per NMR with correlation time tau( perpendicular), is associated with the first and second collective slow (global) GNM modes. The ANM-predicted deformations of the unliganded enzyme conform to the functional reconfiguration induced by ligand-binding, indicating the structural disposition (or potential) of the enzyme to bind its substrates. It is shown that NMR/SRLS and GNM/ANM analyses can be advantageously synthesized to provide insights into the molecular mechanisms that control biological function., ((c) 2004 Wiley-Liss, Inc.) more...

Understanding the molecular mechanisms of HIV-1 reverse transcriptase (RT) action and drug inhibition is essential for designing effective antiretroviral therapies. Although comparisons of the different crystal forms of RT give insights into the flexibility of different domains, a direct computational assessment of the effect of inhibitor binding on the collective dynamics of RT is lacking. A structure-based approach is used here for exploring the dynamics of RT in unliganded and inhibitor-bound forms. Non-nucleoside RT inhibitors (NNRTI) are shown to interfere directly with the global hinge-bending mechanism that controls the cooperative motions of the p66 fingers and thumb subdomains. The net effect of nevirapine binding is to change the direction of domain movements rather than suppress their mobilities. The second generation NNRTI, efavirenz, on the other hand, shows the stronger effect of simultaneously reorienting domain motions and obstructing the p66 thumb fluctuations. A second hinge site controlling the global rotational reorientations of the RNase H domain is identified, which could serve as a target for potential inhibitors of RNase H activity., (Copyright 2002 Wiley-Liss, Inc.) more...