Dielectric properties of proteins from simulation: the effects of solvent, ligands, pH, and temperature.

We have used a standard Fröhlich-Kirkwood dipole moment fluctuation model to calculate the static dielectric permittivity, epsilon(0), for four different proteins, each of which was simulated under at least two different conditions of pH, temperature, solvation, or ligand binding. For the range of proteins and conditions studied, we calculate values for epsilon(0) between 15 and 40. Our results show, in agreement with prior work, that the behavior of charged residues is the primary determinant of the effective permittivity. Furthermore, only environmental changes that alter the properties of charged residues exert a significant effect on epsilon. In contrast, buried water molecules or ligands have little or no effect on protein dielectric properties.

Simulations of the estrogen receptor ligand-binding domain: affinity of natural ligands and xenoestrogens.

We have carried out molecular dynamics (MD) simulations and free energy calculations on the alpha-subtype of the human estrogen receptor ligand-binding domain (ERalpha LBD) complexed with a number of known agonists and putative xenoestrogens. Our dynamical simulations of ligand-receptor complexes underscore the highly structured nature of the complex and offer some interesting insights into the structure-activity relationship (SAR) for these ligands. With traditional thermodynamic integration (TI) calculations, we calculate relative binding free energies for three known agonists, in good agreement with experimental values. The sheer number of possible xenoestrogenic compounds makes an approach using traditional free energy calculations unfeasible. Instead, we have made use of a single-step perturbation methodology that allows the calculation of relative free energies for a large number of related polyaromatic hydrocarbons (PAHs) from a single simulation. Our results show good (maximum deviation 3.3 kJ mol(-1)) agreement with experimental data, suggesting the possibility of large-scale xenoestrogen screening in silico to obtain strongly estrogenic compounds for subsequent experimental testing.

Exhaustive mutagenesis in silico: multicoordinate free energy calculations on proteins and peptides.

We have extended and applied a multicoordinate free energy method, chemical Monte Carlo/Molecular Dynamics (CMC/MD), to calculate the relative free energies of different amino acid side-chains. CMC/MD allows the calculation of the relative free energies for many chemical species from a single free energy calculation. We have previously shown its utility in host:guest chemistry (Pitera and Kollman, J Am Chem Soc 1998;120:7557-7567)1 and ligand design (Eriksson et al., J Med Chem 1999;42:868-881)2, and here demonstrate its utility in calculations of amino acid properties and protein stability. We first study the relative solvation free energies of N-methylated and acetylated alanine, valine, and serine amino acids. With careful inclusion of rotameric states, internal energies, and both the solution and vacuum states of the calculation, we calculate relative solvation free energies in good agreement with thermodynamic integration (TI) calculations. Interestingly, we find that a significant amount of the unfavorable solvation of valine seen in prior work (Sun et al., J Am Chem Soc 1992;114:6798-6801)3 is caused by restraining the backbone in an extended conformation. In contrast, the solvation free energy of serine is calculated to be less favorable than expected from experiment, due to the formation of a favorable intramolecular hydrogen bond in the vacuum state. These monomer calculations emphasize the need to accurately consider all significant conformations of flexible molecules in free energy calculations. This development of the CMC/MD method paves the way for computations of protein stability analogous to the biochemical technique of "exhaustive mutagenesis." We have carried out just such a calculation at position 133 of T4 lysozyme, where we use CMC/MD to calculate the relative stability of eight different side-chain mutants in a single free energy calculation. Our T4 calculations show good agreement with the prior free energy calculations of Veenstra et al. (Prot Eng 1997;10:789-807)4 and excellent agreement with the experiments of Mendel et al. (Science 1992;256:1798-1802).

Understanding substrate specificity in human and parasite phosphoribosyltransferases through calculation and experiment.

We present molecular dynamics (MD) simulations on two enzymes: a human hypoxanthine-guanine-phosphoribosyltransferase (HGPRTase) and its analogue in the protozoan parasite Tritrichomonas foetus. The parasite enzyme has an additional ability to process xanthine as a substrate, making it a hypoxanthine-guanine-xanthine phosphoribosyltransferase (HGXPRTase) [Chin, M. S., and Wang, C. C. (1994) Mol. Biochem. Parasitol. 63 (2), 221-229 (1)]. X-ray crystal structures of both enzymes complexed to guanine monoribosyl phosphate (GMP) have been solved, and show only subtle differences in the two active sites [Eads et al. (1994) Cell 78 (2), 325-334 (2); Somoza et al. (1996) Biochemistry 35 (22), 7032-7040 (3)]. Most of the direct contacts with the base region of the substrate are made by the protein backbone, complicating the identification of residues significantly associated with xanthine recognition. Our calculations suggest that the broader specificity of the parasite enzyme is due to a significantly more flexible base-binding region, and rationalize the effect of two mutations, R155E and D163N, that alter substrate specificity [Munagala, N. R., and Wang, C. C. (1998) Biochemistry 37 (47), 16612-16619 (4)]. In addition, our simulations suggested a double mutant (D106E/D163N) that might rescue the D163N mutant. This double mutant was expressed and assayed, and its catalytic activity was confirmed. Our molecular dynamics trajectories were also used with a structure-based design program, Pictorial Representation Of Free Energy Changes (PROFEC), to suggest parasite-selective derivatives of GMP. Our calculations here successfully rationalize the parasite-selectivity of two novel inhibitors derived from the computer-aided design of Somoza et al. (5) and demonstrate the utility of PROFEC in the design of species-selective inhibitors.

Biochemical characterization and crystallographic structure of an Escherichia coli protein from the phosphotriesterase gene family.

Phosphotriesterase homology protein (PHP) is a member of a recently discovered family of proteins related to phosphotriesterase, a hydrolytic, bacterial enzyme with an unusual substrate specificity for synthetic organophosphate triesters and phosphorofluoridates, which are common constituents of chemical warfare agents and agricultural pesticides. No natural substrate has been identified for phosphotriesterase, and it has been suggested that the enzyme may have evolved the ability to hydrolyze synthetic compounds in bacteria under selective pressure to meet nutritional needs. PHP, which has 28% sequence identity with phosphotriesterase, may belong to the family of proteins from which phosphotriesterase evolved. Here we report the cloning, expression, initial characterization, and high-resolution X-ray crystallographic structure of PHP. Biochemical analysis shows that PHP is monomeric and binds two zinc ions per monomer. Unlike phosphotriesterase, PHP does not catalyze the hydrolysis of nonspecific phosphotriesters. The structure, similar to that of phosphotriesterase, consists of a long, elliptical alpha/beta barrel and has a binuclear zinc center in a cleft at the carboxy end of the barrel at the location of the presumptive active site.