Docking of ATP to Ca-ATPase: considering protein domain motions.

Most standard molecular docking algorithms take into account only ligand flexibility, while numerous studies demonstrate that receptor flexibility may be also important. While some efficient methods have been proposed to take into account local flexibility of protein side chains, the influence of large-scale domain motions on the docking results still represents a challenge for computational methods. In this work we compared the results of ATP docking to different models of Ca-ATPase: crystallographic apo- and holo-forms of the enzyme as well as "flexible" target models generated via molecular dynamics (MD) simulations in water. MD simulations were performed for two different apo-forms and one holo-form of Ca2+-ATPase and reveal large-scale domain motions of type "closure", which is consistent with experimental structures. Docking to a set of MD-conformers yielded correct solutions with ATP bound in both domains regardless of the starting Ca2+-ATPase structure. Also, special attention was paid to proper ranking of docking solutions and some particular features of different scoring functions and their applicability for the model of "flexible" receptor. Particularly, the results of docking ATP were ranked by a scoring criterion specially designed to estimate ATP-protein interactions. This criterion includes stacking and hydrophobic interactions characteristic of ATP-protein complexes. The performance of this ligand-specific scoring function was considerably better than that of a standard scoring function used in the docking algorithm.

Complementarity of hydrophobic properties in ATP-protein binding: a new criterion to rank docking solutions.

ATP is an important substrate of numerous biochemical reactions in living cells. Molecular recognition of this ligand by proteins is very important for understanding enzymatic mechanisms. Considerable insight into the problem may be gained via molecular docking simulations. At the same time, standard docking protocols are often insufficient to predict correct conformations for protein-ATP complexes. Thus, in most cases the native-like solutions can be found among the docking poses, but current scoring functions have only limited ability to discriminate them from false positives. To improve the selection of correct docking solutions obtained with the GOLD software, we developed a new ranking criterion specific for ATP-protein binding. The method is based on detailed analysis of the intermolecular interactions in 40 high-resolution 3D structures of ATP-protein complexes (the training set). We found that the most important factors governing this recognition are hydrogen-bonding, stacking between adenine and aromatic protein residues, and hydrophobic contacts between adenine and protein residues. To address the latter, we applied the formalism of 3D molecular hydrophobicity potential. The results obtained were used to construct an ATP-oriented scoring criterion as a linear combination of the terms describing these intermolecular interactions. The criterion was then validated using the test set of 10 additional ATP-protein complexes. As compared with the standard scoring functions, the new ranking criterion significantly improved the selection of correct docking solutions in both sets and allowed considerable enrichment at the top of the list containing docking poses with correct solutions.

Development of the force field parameters for phosphoimidazole and phosphohistidine.

Phosphorylation of histidine-containing proteins is a key step in the mechanism of many phosphate transfer enzymes (kinases, phosphatases) and is the first stage in a wide variety of signal transduction cascades in bacteria, yeast, higher plants, and mammals. Studies of structural and dynamical aspects of such enzymes in the phosphorylated intermediate states are important for understanding the intimate molecular mechanisms of their functioning. Such information may be obtained via molecular dynamics and/or docking simulations, but in this case appropriate force field parameters for phosphohistidine should be explicitly defined. In the present article we describe development of the GROMOS96 force field parameters for phosphoimidazole molecule--a realistic model of the phosphohistidine side chain. The parameterization is based on the results of ab initio quantum chemical calculations with subsequent refinement and testing using molecular mechanics and molecular dynamics simulations. The set of force constants and equilibrium geometry is employed to derive force field for the phosphohistidine moiety. Resulting parameters and topology are incorporated into the molecular modeling package GROMACS and used in molecular dynamics simulations of a phosphohistidine-containing protein in explicit solvent.

The distinct functional properties of the nucleotide-binding domain of ATP7B, the human copper-transporting ATPase: analysis of the Wilson disease mutations E1064A, H1069Q, R1151H, and C1104F.

Copper transport by the P(1)-ATPase ATP7B, or Wilson disease protein (WNDP),1 is essential for human metabolism. Perturbation of WNDP function causes intracellular copper accumulation and severe pathology, known as Wilson disease (WD). Several WD mutations are clustered within the WNDP nucleotide-binding domain (N-domain), where they are predicted to disrupt ATP binding. The mechanism by which the N-domain coordinates ATP is presently unknown, because residues important for nucleotide binding in the better characterized P(2)-ATPases are not conserved within the P(1)-ATPase subfamily. To gain insight into nucleotide binding under normal and disease conditions, we generated the recombinant WNDP N-domain and several WD mutants. Using isothermal titration calorimetry, we demonstrate that the N-domain binds ATP in a Mg(2+)-independent manner with a relatively high affinity of 75 microm, compared with millimolar affinities observed for the P(2)-ATPase N-domains. The WNDP N-domain shows minimal discrimination between ATP, ADP, and AMP, yet discriminates well between ATP and GTP. Similar results were obtained for the N-domain of ATP7A, another P(1)-ATPase. Mutations of the invariant WNDP residues E1064A and H1069Q drastically reduce nucleotide affinities, pointing to the likely role of these residues in nucleotide coordination. In contrast, the R1151H mutant exhibits only a 1.3-fold reduction in affinity for ATP. The C1104F mutation significantly alters protein folding, whereas C1104A does not affect the structure or function of the N-domain. Together, the results directly demonstrate the phenotypic diversity of WD mutations within the N-domain and indicate that the nucleotide-binding properties of the P(1)-ATPases are distinct from those of the P(2)-ATPases.

Molecular modelling of the nucleotide-binding domain of Wilson's disease protein: location of the ATP-binding site, domain dynamics and potential effects of the major disease mutations.

WNDP (Wilson's disease protein) is a copper-transporting ATPase that plays an essential role in human physiology. Mutations in WNDP result in copper accumulation in tissues and cause a severe hepato-neurological disorder known as Wilson's disease. Several mutations were surmised to affect the nucleotide binding and hydrolysis by WNDP; however, how the nucleotides bind to normal and mutated WNDP remains unknown. To aid such studies, we performed the molecular modelling of the spatial structure and dynamics of the ATP-binding domain of WNDP and its interactions with ATP. The three-dimensional models of this domain in two conformations were built using the X-ray structures of the Ca2+-ATPase in the E1 and E2 states. To study the functional aspects of the models, they were subjected to long-term molecular dynamics simulations in an explicit solvent; similar calculations were performed for the ATP-binding domain of Ca2+-ATPase. In both cases, we found large-scale motions that lead to significant changes of distances between several functionally important residues. The ATP docking revealed two possible modes of ATP binding: via adenosine buried in the cleft near residues H1069, R1151 and D1164, and via phosphate moiety 'anchored' by H-bonds with residues in the vicinity of catalytic D1027. Furthermore, interaction of ATP with both sites occurs if they are spatially close to each other. This may be achieved after relative domain motions of the 'closure' type observed in molecular dynamics simulations. The results provide a framework for analysis of disease mutations and for future mutagenesis studies.