Simulation of pH-Dependent Properties of Proteins Using Mesoscopic Models. In: Lipkowitz KB, Boyd DB, editors. Reviews in Computational Chemistry. Hoboken: John Wiley & Sons, Inc.; 1999. pp. 249-311. [DOI: 10.1002/9780470125908.ch5]

Searching for Hydrodynamic Orienting Effects in the Association of Tri-N-acetylglucosamine with Hen Egg-White Lysozyme

Using stopped-flow fluorometry, we determined rate constants for the formation of diffusional encounter complexes of tri-N-acetylglucosamine (NAG₃) with hen egg-white lysozyme (k_a^WT) and its double mutant Asp48Asn/Lys116Gln (k_a^MT). We defined binding anisotropy, κ ≡ (k_a^WT – k_a^MT)/(k_a^WT + k_a^MT), and determined its ionic strength dependence. Our goal was to check if this ionic strength dependence provides information about the orienting hydrodynamic effects in the ligand-binding process. We also computed ionic strength dependence of the binding anisotropy from Brownian dynamics simulations using simple models of the lysozyme–NAG₃ system. The results of our experiments indicate that in the case of lysozyme and NAG₃ such hydrodynamic orienting effects are rather negligible. On the other hand, the results of our Brownian dynamics simulations prove that there exist molecular systems for which such orienting effects are substantial. However, the ionic strength dependence of the rate constants for the wild-type and modified systems do not exhibit any qualitative features that would allow us to conclude the presence of hydrodynamic orienting effects from stopped-flow experiments alone. Nevertheless, the results of our simulations suggest the presence of hydrodynamic orienting effects in the receptor–ligand association when the anisotropy of binding depends on the solvent viscosity.

Poisson-Boltzmann model analysis of binding mRNA cap analogues to the translation initiation factor eIF4E

The electrostatic free energy of binding of two analogues of the 5'-mRNA cap, differing in size and electric charge, to the wild type and mutated eukaryotic initiation factor eIF4E was computed using the finite difference solutions to the Poisson-Boltzmann equation. Two definitions of the solute-solvent dielectric boundary were used: van der Waals model, solvent exclusion (SE) model. The computed electrostatic energies were supplemented by estimations of the non polar and entropic contributions. A comparison with experimental data for the investigated systems was done. It appears that the SE model with additional contribution fits experimental findings better than the van der Waals model does.

Experimental and calculated shift in pK(a) upon binding of phosphotyrosine peptide to the SH2 domain of p56(lck)

The pH dependence of the affinity of a 11-mer phosphotyrosine (pY) peptide (EPQpYEEIPIYL-NH2) for the SH2 domain of the tyrosine kinase p56(lck) was investigated with surface plasmon resonance (SPR). From SPR competition experiments the affinity in solution was obtained. The pH dependence of the affinity in solution can be well described by a proton linkage model with a single pK(a) shift upon binding, from 6.1 to 4.7. This shift is ascribed to the transition from the -2 to the -1 ionisation state of the tyrosine phosphate group. Based on the X-ray structure for the complex with Lck SH2, a pK(a) value of 5.3 for the bound pY peptide was computed, modelling the solvated protein as a system of point charges in a continuum. With the phosphate in the -2 state the binding energy is 1.8 kcal/mol more favourable than for the -1 state, corresponding to a 20-fold higher affinity. A proper charge is relevant in the design of potential therapeutic Lck SH2 ligands with mimics for the metabolically unstable tyrosine phosphate group.

Revisiting the structural flexibility of the complex p21(ras)-GTP: the catalytic conformation of the molecular switch II

The hydrolysis of GTP in p21(ras) triggers conformational changes that regulate the ras/ERK signaling pathway. An important active site residue is Gln61, which has been found to be mutated in 30% of human tumors. The dynamics of the active site conformation is studied by using molecular dynamics simulation of two independent structures of the GTP-bound uncomplexed enzyme. Two distinct conformations of the enzyme are observed, in which the side-chain residue Gln61 is in different orientations. Essential dynamics analysis is used to describe the essential motions in the transition between the two conformations. Results are compared with earlier simulations of p21(ras) and its complex with GTPase activating protein p21-GAP.

Ionization state and molecular docking studies for the macrophage migration inhibitory factor: the role of lysine 32 in the catalytic mechanism

The macrophage migration inhibitory factor (MIF) is a cytokine that is structurally similar to certain isomerases and for which multiple immune and catalytic roles have been proposed. Different catalytic activities have been reported for MIF, yet the exact mechanism by which MIF acts is not completely known. As a tautomerase, the enzyme uses a general acid-base mechanism of proton transfer in which the amino-terminal proline has been shown to function as the catalytic base. We report the results of molecular docking simulations of macrophage migration inhibitory factor with three substrates, D-dopachrome, L-dopachrome methyl ester and p-(hydroxyphenyl)pyruvate. Electrostatic pK(a) predictions were also performed for the free and complexed forms of the enzyme. The predicted binding mode of p-(hydroxyphenyl)pyruvate is in agreement with the recently published X-ray structure. A model for the binding mode of D-dopachrome and L-dopachrome methyl ester to MIF is proposed which offers insights into the catalytic mechanism of D-dopachrome tautomerase activity of MIF. The proposed catalytic mechanism is further supported by the pK(a) predictions, which suggest that residue Lys32 acts as the general acid for the enzymatic catalysis of D-dopachrome.