Calculation of protein-ligand binding free energy by using a polarizable potential.

The binding of charged ligands benzamidine and diazamidine to trypsin was investigated by using a polarizable potential energy function and explicit-water molecular dynamics simulations. The binding free energies were computed from the difference between the free energies of decoupling the ligand from water and protein environments. Both the absolute and the relative free energies from the perturbation simulations agree with experimental measurements to within 0.5 kcal.mol(-1). Comparison of free-energy components sampled from different thermodynamic paths indicates that electrostatics is the main driving force behind benzamidine recognition of trypsin. The contribution of electronic polarization to binding appears to be crucial. By computing the free-energy contribution caused by the polarization between the ligand and its surroundings, we found that polarization has the opposite effect in dissimilar environments. Although polarization favors ligand solvation in water, it weakens the protein-ligand attraction by screening the electrostatic interaction between trypsin and benzamidine. We also examined the relative binding free energies of a benzamidine analog diazamidine to trypsin. The changes in free energy on benzamidine-diazamidine substitution were tens of kilocalories in both water and trypsin environments; however, the change in the total binding free energy is <2 kcal.mol(-1) because of cancellation, consistent with the experimental results. Overall, our results suggest that the use of a polarizable force field, given adequate sampling, is capable of achieving chemical accuracy in molecular simulations of protein-ligand recognition.

A transferable coarse-grained model for hydrogen-bonding liquids.

We present here a recent development of a generalized coarse-grained model for use in molecular simulations. In this model, interactions between coarse-grained particles consist of both van der Waals and explicit electrostatic components. As a result, the coarse-grained model offers the transferability that is lacked by most current effective-potential based approaches. The previous center-of-mass framework (P. A. Golubkov and P. Ren, J. Chem. Phys., 2006, 125, 64103) is generalized here to include arbitrary off-center interaction sites for both Gay-Berne and multipoles. The new model has been applied to molecular dynamic simulations of neat methanol liquid. By placing a single point multipole at the oxygen atom rather than at the center of mass of methanol, there is a significant improvement in the ability to capture hydrogen-bonding. The critical issue of transferability of the coarse-grained model is verified on methanol-water mixtures, using parameters derived from neat liquids without any modification. The mixture density and internal energy from coarse-grained molecular dynamics simulations show good agreement with experimental measurements, on a par with what has been obtained from more detailed atomic models. By mapping the dynamics trajectory from the coarse-grained simulation into the all-atom counterpart, we are able to investigate atomic-level structure and interaction. Atomic radial distribution functions of neat methanol, neat water and mixtures compare favorably to experimental measurements. Furthermore, hydrogen-bonded 6- and 7-molecule chains of water and methanol observed in the mixture are in agreement with previous atomic simulations.

Substrate specificity of human kallikreins 1 and 6 determined by phage display.

The human tissue kallikrein (KLK) family contains 15 secreted serine proteases that are expressed in a wide range of tissues and have been implicated in different physiological functions and disease states. Of these, KLK1 has been shown to be involved in the regulation of multiple physiological processes such as blood pressure, smooth muscle contraction, and vascular cell growth. KLK6 is overexpressed in breast and ovarian cancer tissues and has been shown to cleave peptide derived from human myelin protein and Abeta amyloid peptide in vitro. Here we analyzed the substrate specificity of KLK1 and KLK6, by substrate phage display using a random octapeptide library. Consistent with earlier biochemical data, KLK1 was shown to exhibit both trypsin- and chymotrypsin-like selectivities with Tyr/Arg preferred at site P1, Ser/Arg strongly preferred at P1', and Phe/Leu at P2. KLK6 displayed trypsin-like activity, with the P1 position occupied only by Arg and a strong preference for Ser in P1'. Docking simulations of consensus peptide provide information on the identity of the enzyme residues that are responsible for substrate binding. Bioinformatic analysis suggested several putative KLK6 protein substrates, such as ionotropic glutamate receptor (GluR) and synphilin.

Generalized coarse-grained model based on point multipole and Gay-Berne potentials.

This paper presents a general coarse-grained molecular mechanics model based on electric point multipole expansion and Gay-Berne [J. Chem. Phys. 74, 3316 (1981)] potential. Coarse graining of van der Waals potential is achieved by treating molecules as soft uniaxial ellipsoids interacting via a generalized anisotropic Gay-Berne function. The charge distribution is represented by point multipole expansion, including point charge, dipole, and quadrupole moments placed at the center of mass. The Gay-Berne and point multipole potentials are combined in the local reference frame defined by the inertial frame of the all-atom counterpart. The coarse-grained model has been applied to rigid-body molecular dynamics simulations of molecular liquids including benzene and methanol. The computational efficiency is improved by several orders of magnitude, while the results are in reasonable agreement with all-atom models and experimental data. We also discuss the implications of using point multipole for polar molecules capable of hydrogen bonding and the applicability of this model to a broad range of molecular systems including highly charged biopolymers.

Inactivation of the phenylpyruvate tautomerase activity of macrophage migration inhibitory factor by 2-oxo-4-phenyl-3-butynoate.

Macrophage migration inhibitory factor (MIF) is an important immunoregulatory protein that has been implicated in several inflammatory diseases. MIF also has a phenylpyruvate tautomerase (PPT) activity, the role of which remains elusive in these biological activities. The acetylene compound, 2-oxo-4-phenyl-3-butynoate (2-OPB), has been synthesized and tested as a potential irreversible inhibitor of its enzymatic activity. Incubation of the compound with MIF results in the rapid and irreversible loss of the PPT activity. Mass spectral analysis established that the amino-terminal proline, previously implicated as a catalytic base in the PPT-catalyzed reaction, is the site of covalent modification. Inactivation of the PPT activity likely occurs by a Michael addition of Pro-1 to C-4 of the inhibitor. Attempts to crystallize the inactivated complex to confirm the structure of the adduct on the covalently modified Pro-1 by X-ray crystallography were not successful. Nor was it possible to unambiguously interpret electron density observed in the active sites of the native crystals soaked with the inhibitor. This may be due to crystal packing in that the side chain of Glu-16 from an adjacent trimer occupies one active site. However, this crystal contact may be partially responsible for the high-resolution quality of these MIF crystals. Nonetheless, because MIF is a member of the tautomerase superfamily, a group of structurally homologous proteins that share a beta-alpha-beta structural motif and a catalytic Pro-1, 2-OPB may find general use as a probe of tautomerase superfamily members that function as PPTs.