Calculating pKa values in the cAMP-dependent protein kinase: the effect of conformational change and ligand binding

Protonation and pK changes in protein-ligand binding

Formation of protein-ligand complexes causes various changes in both the receptor and the ligand. This review focuses on changes in pK and protonation states of ionizable groups that accompany protein-ligand binding. Physical origins of these effects are outlined, followed by a brief overview of the computational methods to predict them and the associated corrections to receptor-ligand binding affinities. Statistical prevalence, magnitude and spatial distribution of the pK and protonation state changes in protein-ligand binding are discussed in detail, based on both experimental and theoretical studies. While there is no doubt that these changes occur, they do not occur all the time; the estimated prevalence varies, both between individual complexes and by method. The changes occur not only in the immediate vicinity of the interface but also sometimes far away. When receptor-ligand binding is associated with protonation state change at particular pH, the binding becomes pH dependent: we review the interplay between sub-cellular characteristic pH and optimum pH of receptor-ligand binding. It is pointed out that there is a tendency for protonation state changes upon binding to be minimal at physiologically relevant pH for each complex (no net proton uptake/release), suggesting that native receptor-ligand interactions have evolved to reduce the energy cost associated with ionization changes. As a result, previously reported statistical prevalence of these changes - typically computed at the same pH for all complexes - may be higher than what may be expected at optimum pH specific to each complex. We also discuss whether proper account of protonation state changes appears to improve practical docking and scoring outcomes relevant to structure-based drug design. An overview of some of the existing challenges in the field is provided in conclusion.

Small-molecule screening identifies modulators of aquaporin-2 trafficking

In the principal cells of the renal collecting duct, arginine vasopressin (AVP) stimulates the synthesis of cAMP, leading to signaling events that culminate in the phosphorylation of aquaporin-2 water channels and their redistribution from intracellular domains to the plasma membrane via vesicular trafficking. The molecular mechanisms that control aquaporin-2 trafficking and the consequent water reabsorption, however, are not completely understood. Here, we used a cell-based assay and automated immunofluorescence microscopy to screen 17,700 small molecules for inhibitors of the cAMP-dependent redistribution of aquaporin-2. This approach identified 17 inhibitors, including 4-acetyldiphyllin, a selective blocker of vacuolar H(+)-ATPase that increases the pH of intracellular vesicles and causes accumulation of aquaporin-2 in the Golgi compartment. Although 4-acetyldiphyllin did not inhibit forskolin-induced increases in cAMP formation and downstream activation of protein kinase A (PKA), it did prevent cAMP/PKA-dependent phosphorylation at serine 256 of aquaporin-2, which triggers the redistribution to the plasma membrane. It did not, however, prevent cAMP-induced changes to the phosphorylation status at serines 261 or 269. Last, we identified the fungicide fluconazole as an inhibitor of cAMP-mediated redistribution of aquaporin-2, but its target in this pathway remains unknown. In conclusion, our screening approach provides a method to begin dissecting molecular mechanisms underlying AVP-mediated water reabsorption, evidenced by our identification of 4-acetyldiphyllin as a modulator of aquaporin-2 trafficking.

pKa modulation of the acid/base catalyst within GH32 and GH68: a role in substrate/inhibitor specificity?

Glycoside hydrolases of families 32 (GH32) and 68 (GH68) belong to clan GH-J, containing hydrolytic enzymes (sucrose/fructans as donor substrates) and fructosyltransferases (sucrose/fructans as donor and acceptor substrates). In GH32 members, some of the sugar substrates can also function as inhibitors, this regulatory aspect further adding to the complexity in enzyme functionalities within this family. Although 3D structural information becomes increasingly available within this clan and huge progress has been made on structure-function relationships, it is not clear why some sugars bind as inhibitors without being catalyzed. Conserved aspartate and glutamate residues are well known to act as nucleophile and acid/bases within this clan. Based on the available 3D structures of enzymes and enzyme-ligand complexes as well as docking simulations, we calculated the pKa of the acid-base before and after substrate binding. The obtained results strongly suggest that most GH-J members show an acid-base catalyst that is not sufficiently protonated before ligand entrance, while the acid-base can be fully protonated when a substrate, but not an inhibitor, enters the catalytic pocket. This provides a new mechanistic insight aiming at understanding the complex substrate and inhibitor specificities observed within the GH-J clan. Moreover, besides the effect of substrate entrance on its own, we strongly suggest that a highly conserved arginine residue (in the RDP motif) rather than the previously proposed Tyr motif (not conserved) provides the proton to increase the pKa of the acid-base catalyst.

Predicting the open conformations of protein kinases using molecular dynamics simulations

Protein kinases (PK) control phosphorylation in eukaryotic cells, and thereby regulate metabolic pathways, cell cycle progression, apoptosis, and transcription. Consequently, there is significant interest in manipulating PK activity and treat diseases by using small-molecule drugs. All PK catalytic domains undergo large conformational changes as a result of substrate binding and phosphorylation. The "closed" state of a PK catalytic domain is the only state able to phosphorylate the target substrate, which makes the two other observed states (the "open" and the "intermediate" states) interesting drug targets. We investigate whether molecular dynamics (MD) simulations starting from the closed state of the catalytic domain of protein kinase A (C-PKA) can be used to produce realistic structures representing the intermediate and/or open conformation of C-PKA, because this would allow for drug docking calculations and drug design using MD snapshots. We perform 36 ten-nanosecond MD simulations starting from the closed conformation [PDB ID: ATP] of C-PKA in various liganded and phosphorylated states. The results show that MD simulations are capable of reproducing the open conformation of C-PKA with good accuracy within 1 ns of simulation as measured by Cα root mean square deviations (RMSDs) and RMSDs of atoms defining the ATP-binding pocket. Importantly, we are able to show that even without knowledge of the structure of the open form of C-PKA, we can identify the MD snapshots resembling the open conformation most using the open structure of a different PK displaying only 23% sequence identity to C-PKA.