Homology modeling of histidine-containing phosphocarrier protein and eosinophil-derived neurotoxin: construction of models and comparison with experiment.

Homology modeling methods have been used to construct models of two proteins--the histidine-containing phosphocarrier protein (HPr) from Mycoplasma capricolum and human eosinophil-derived neurotoxin (EDN). Comparison of the models with the subsequently determined X-ray crystal structures indicates that the core regions of both proteins are reasonably well reproduced, although the template structures are closer to the X-ray structures in these regions--possible enhancements are discussed. The conformations of most of the side chains in the core of HPr are well reproduced in the modeled structure. As expected, the conformations of surface side chains in this protein differ significantly from the X-ray structure. The loop regions of EDN were incorrectly modeled--reasons for this and possible enhancements are discussed.

On achieving better than 1-A accuracy in a simulation of a large protein: Streptomyces griseus protease A.

Computational methods are frequently used to simulate the properties of proteins. In these studies accuracy is clearly important, and the improvement of accuracy of protein simulation methodology is one of the major challenges in the application of theoretical methods, such as molecular dynamics, to structural studies of biological molecules. Much effort is being devoted to such improvements. Here, we present an analysis of a 187-ps molecular dynamics simulation of the serine protease Streptomyces griseus protease A in its crystal environment. The reproduction of the experimental structure is considerably better than has been achieved in earlier simulations--the root mean square deviation of the simulated structure from the x-ray structure being less than 1 A, a significant step toward the goal of simulating proteins to within experimental error. The use of a longer cutoff with truncation rather than a switching function, inclusion of all crystalline water and the counterions in the crystallization medium, and use of the consistent valence force field characterize the differences in this calculation.

Theoretical studies on the dihydrofolate reductase mechanism: electronic polarization of bound substrates.

We have applied local density functional theory, an ab initio quantum mechanical method, to study the shift in the spatial electron density of the substrate dihydrofolate that accompanies binding to the enzyme dihydrofolate reductase. The results shed light on fundamental electronic effects due to the enzyme that may contribute to catalysis. In particular, the enzyme induces a long-range polarization of the substrate that perturbs its electron density distribution in a specific and selective way in the vicinity of the bond that is reduced by the enzyme. Examination of the electron density changes that occur in folate reveals that a similar effect is seen but this time specifically at the bond that is reduced in this substrate. This suggests that the polarization effect may be implicated in the reaction mechanism and may play a role in determining the sequence whereby the 7,8-bond in folate is reduced first, followed by reduction of the 5,6-bond in the resulting dihydro compound.

The electrostatic potential of Escherichia coli dihydrofolate reductase.

Escherichia coli dihydrofolate reductase (DHFR) carries a net charge of -10 electrons yet it binds ligands with net charges of -4 (NADPH) and -2 (folate or dihydrofolate). Evaluation and analysis of the electrostatic potential of the enzyme give insight as to how this is accomplished. The results show that the enzyme is covered by an overall negative potential (as expected) except for the ligand binding sites, which are located inside "pockets" of positive potential that enable the enzyme to bind the negatively charged ligands. The electrostatic potential can be related to the asymmetric distribution of charged residues in the enzyme. The asymmetric charge distribution, along with the dielectric boundary that occurs at the solvent-protein interface, is analogous to the situation occurring in superoxide dismutase. Thus DHFR is another case where the shape of the active site focuses electric fields out into solution. The positive electrostatic potential at the entrance of the ligand binding site in E. coli DHFR is shown to be a direct consequence of the presence of three positively charged residues at positions 32, 52, and 57--residues which have also been shown recently to contribute significantly to electronic polarization of the ligand folate. The latter has been postulated to be involved in the catalytic process. A similar structural motif of three positively charged amino acids that gives rise to a positive potential at the entrance to the active site is also found in DHFR from chicken liver, and is suggested to be a common feature in DHFRs from many species. It is noted that, although the net charges of DHFRs from different species vary from +3 to -10, the enzymes are able to bind the same negatively charged ligands, and perform the same catalytic function.

Predicted three-dimensional structure of the protease inhibitor domain of the Alzheimer's disease beta-amyloid precursor.

Alzheimer's disease is characterized by the deposition of amyloid beta-protein as plaques and tangles in the brains of its victims. The amyloid precursor can be expressed with or without the inclusion of a protease inhibitor domain, the potential role of which in amyloidogenesis has prompted the generation of a model of its three-dimensional structure based on the known structure of a related inhibitor. The model structure predicts that the mutated residues are almost entirely on the surface of the inhibitor domain, while conserved residues constitute the hydrophobic core. In addition, several pairs of structurally complementary, or concerted, mutations are seen. These structural features provide strong evidence for the validity of the modeled structure, and it is suggested that the presence of complementary mutations may be used as a criterion for evaluating protein structures built by homology, in addition to the (spatial) location of the mutations. The terminal residues delimiting the domain are among those furthest from the protease binding site and are in close proximity to one another, thus suggesting the ability of the domain to function as a structural cassette within the context of a larger protein. The electrostatic potentials of the inhibitor and of the related bovine pancreatic trypsin inhibitor reveal how two inhibitors with very different net charges can bind with approximately the same binding constant to trypsin and suggest a mutation of trypsin that might selectively enhance the binding of the amyloid inhibitor domain. The model provides a structural basis for understanding the functional roles of residues in the domain and for designing simpler molecules to test as pharmacologic agents for intervention in Alzheimer's disease.

Electron redistribution on binding of a substrate to an enzyme: folate and dihydrofolate reductase.

The migration of electron density of a substrate (folate) on binding to an enzyme (dihydrofolate reductase) is studied by a quantum-mechanical method originally developed in solid state physics. A significant polarization of the substrate is induced by the enzyme, toward the transition state of the enzymatic reaction, at the same time giving rise to "electronic strain energy" in the substrate and enhanced protein-ligand interactions. The spatial arrangement of protein charges that induces the polarization is identified and found to be structurally conserved for bacterial and vertebrate dihydrofolate reductases.

Changes in the electron density of the cofactor NADPH on binding to E. coli dihydrofolate reductase.

Quantum-mechanical electron density calculations reveal that a significant polarization is induced in the cofactor NADPH (reduced nicotinamide adenine dinucleotide phosphate) on binding to the enzyme dihydrofolate reductase. The calculations indicate that electron density corresponding to approximately 0.7 electron charges is shifted within the molecule, extending over more than 20 A. Further calculations on proposed enzyme mutants show that the polarization of NADPH on binding to DHFR is, in large part, induced by a motif of three positively charged residues. This motif was also identified to be directly responsible for the positive electrostatic potential surrounding the cofactor binding site in the enzyme. The possibility of this long-range polarization of NADPH was originally proposed based on a previous study of ligand binding to DHFR where a conserved structural motif of three positively charged residues was found to play a major role in polarizing the substrate folate over its entire length of 18 A.

Molecular dynamics study of the structure and dynamics of a protein molecule in a crystalline ionic environment, Streptomyces griseus protease A.

A large-scale molecular dynamics simulation of the behavior of a serine protease (Streptomyces griseus protease A) in a crystalline environment has been performed. All atoms (including hydrogens) of two protein molecules and the surrounding solvent of crystallization, consisting of both water and salt ions, were explicitly represented, and a relatively long range of interactions (up to 15 A) were included. The simulation is the longest so far reported for a protein in such an environment (60 ps). The use of the full crystalline environment allows a direct comparison of the structure and dynamic properties of the protein and surrounding solvent to be made with the experimental X-ray structure. Here we report the comparison of the protein structures and analyze the energetics of the system, including interaction with the aqueous environment. Subsequent papers will deal with other aspects of the simulation. The overall root mean square differences between the time-averaged molecular dynamics structure and that from crystallography, for all well-ordered, non-hydrogen atoms, are 1.67 and 1.25 A for the two molecules taken as the asymmetric unit. An extensive analysis of the conformation of substructural elements and individual residues and their deviation from experiment has revealed a strong influence of the ionic medium on their behavior. Implications of the results for free energy calculations and for future directions are also discussed.

Catalysis of a rotational transition in a peptide by crystal forces.

Detailed examination of the dynamics trajectories of the isolated cyclic peptide cyclo-(Ala-Pro-D-Phe)2 and of the molecule in its crystalline environment led to the unexpected observation that the methyl groups of the alanine residues rotated more frequently during a simulation in the crystal environment than in a simulation of the isolated peptide. In effect, the crystal environment is "catalyzing" the rotational isomerization of the methyl groups. In order to understand how the crystal forces increase the rate of this rotation, and to explore any possible analogy to the inducing of strained conformations of ligands by enzymes, the barriers to rotation in the two environments were studied by using the torsion angle forcing method. The crystal forces induce a different, higher energy, conformation of the peptide than is found for the isolated molecule, and the different rates of rotation have been explained in terms of the resulting specific intramolecular interactions that, it turns out, give rise to the lower rotational barrier. Molecular dynamics simulations of the peptide were also run at higher temperatures in order to calculate the barriers to rotation through the use of Arrhenius plots. The barriers obtained in this way agree well with the barriers obtained through an adiabatic reaction path derived by rotating the methyl through the barrier while minimizing all remaining degrees of freedom. The rates of rotation calculated from these adiabatic barriers also agree well with the rates observed during the 300 K simulations.

Theoretical studies of the structure and molecular dynamics of a peptide crystal.

Energy minimizations and molecular dynamics simulations have been performed on the cyclic peptide cyclo-(Ala-Pro-D-Phe)2 in both the isolated and crystal states. The results of these calculations have been analyzed, both to investigate our ability to reproduce experimental data (structure and vibrational and NMR spectra) and to investigate the effects of environment on the energy, structure, and dynamics of peptides. Comparison of the minimized and time-averaged crystal systems with the experimental peptide structure shows that the calculations have closely reproduced the experimental structure. Molecular dynamics of the isolated molecule has led to a new conformation, which is approximately equal to 8.5 kcal/mol more stable than the conformation that exists in the crystal, the latter conformation being stabilized by intermolecular (packing) forces. This illustrates the considerable effect that environment can have on the conformation of peptides. The crystal environment has also been shown to significantly reduce the dynamic conformational fluctuations seen for the isolated molecule. The behavior of the peptide during the isolated simulation also supports the experimental NMR observation of a symmetric structure that differs from the asymmetric, instantaneous structures which characterize the molecule during the dynamics. Calculations of vibrational frequencies of the peptide in the crystal and isolated states show the expected shifts in bond-stretching frequencies due to intermolecular interactions. Finally, we have calculated NMR coupling constants from the dynamics simulation of the isolated peptide and have compared these with the experimental values. This has led to a possible reinterpretation of the experimental data.