Calculation of the relative change in binding free energy of a protein-inhibitor complex.

By means of a thermodynamic perturbation method implemented with molecular dynamics, the relative free energy of binding was calculated for the enzyme thermolysin complexed with a pair of phosphonamidate and phosphonate ester inhibitors. The calculated difference in free energy of binding was 4.21 +/- 0.54 kilocalories per mole. This compares well with the experimental value of 4.1 kilocalories per mole. The method is general and can be used to determine a change or "mutation" in any system that can be suitably represented. It is likely to prove useful for protein and drug design.

Free energy calculations by computer simulation.

A fundamental problem in chemistry and biochemistry is understanding the role of solvation in determining molecular properties. Recent advances in statistical mechanical theory and molecular dynamics methodology can be used to solve this problem with the aid of supercomputers. By using these advances the free energies of solvation of all the chemical classes of amino acid side chains, four nucleic acid bases and other organic molecules can be calculated. The effect of a site-specific mutation on the stability of trypsin is predicted. The results are in good agreement with available experiments.

Van der waals surfaces in molecular modeling: implementation with real-time computer graphics.

A method is described for generating van der Waals molecular surfaces with a real-time interactive calligraphic color display system. These surfaces maintain their proper representation during bond rotation and global transformations, and an interior atom removal method yields a comprehensible picture of the molecular surface for large molecules. Both algorithms are faster than previous methods. This combination provides a powerful tool for real-time interactive molecular modeling.

A phylogenetic approach to target selection for structural genomics: solution structure of YciH.

Structural genomics presents an enormous challenge with up to 100 000 protein targets in the human genome alone. At current rates of structure deter-mination, judicious selection of targets is necessary. Here, a phylogenetic approach to target selection is described which makes use of the National Center for Biotechnology Information database of Clusters of Orthologous Groups (COGS). The strategy is designed so that each new protein structure is likely to provide novel sequence-fold information. To demonstrate this approach, the NMR solution structure of YciH (COG0023), a putative translation initiation factor from Escherichia coli, has been determined and its fold classified. YciH is an ortholog of eIF-1/SUI1, an integral component of the translation initiation complex in eukaryotes. The structure consists of two antiparallel alpha-helices packed against the same side of a five-stranded beta-sheet. The first 31 residues of the 11.5 kDa protein are unstructured in solution. Comparative analysis indicates that the folded portion of YciH resembles a number of structures with the alpha-beta plait topology, though its sequence is not homologous to any of them. Thus, the phylogenetic approach to target selection described here was used successfully to identify a new homologous superfamily within this topology.

Computational enzymology.

Numerical simulations of enzyme reaction mechanisms are beginning to provide quantitative as well as qualitative insights. Methods based on hybrid quantum mechanical/molecular mechanical technique permit the natural inclusion of protein solvation effects. Coupled with modern experimental techniques, the numerical simulations are providing details at the atomic level about how enzyme structure influences its function.

Combining Laue diffraction and molecular dynamics to study enzyme intermediates.

Two separate techniques, Laue diffraction and computational molecular dynamics (MD) simulations, have been independently developed to allow the visualization and assessment of transient structural states. Recent studies on isocitrate dehydrogenase show that computational MD simulations of an enzymatic Michaelis complex are consistent with difference Fourier electron density maps of the same structure from a Laue experiment. The use of independent MD studies during crystallographic refinement has allowed us to assign with confidence a number of additional contacts and features important for hydride transfer. We find that unrestrained independent MD simulations provides a very useful method of cross-validation for highly mobile atoms in regions of experimental density that are poorly defined. Likewise, information from Laue difference maps provides information about substrate conformation and interactions that greatly facilitate MD simulations.

Distribution of protein folds in the three superkingdoms of life.

A sensitive protein-fold recognition procedure was developed on the basis of iterative database search using the PSI-BLAST program. A collection of 1193 position-dependent weight matrices that can be used as fold identifiers was produced. In the completely sequenced genomes, folds could be automatically identified for 20%-30% of the proteins, with 3%-6% more detectable by additional analysis of conserved motifs. The distribution of the most common folds is very similar in bacteria and archaea but distinct in eukaryotes. Within the bacteria, this distribution differs between parasitic and free-living species. In all analyzed genomes, the P-loop NTPases are the most abundant fold. In bacteria and archaea, the next most common folds are ferredoxin-like domains, TIM-barrels, and methyltransferases, whereas in eukaryotes, the second to fourth places belong to protein kinases, beta-propellers and TIM-barrels. The observed diversity of protein folds in different proteomes is approximately twice as high as it would be expected from a simple stochastic model describing a proteome as a finite sample from an infinite pool of proteins with an exponential distribution of the fold fractions. Distribution of the number of domains with different folds in one protein fits the geometric model, which is compatible with the evolution of multidomain proteins by random combination of domains. [Fold predictions for proteins from 14 proteomes are available on the World Wide Web at. The FIDs are available by anonymous ftp at the same location.]

Simulation analysis of triose phosphate isomerase: conformational transition and catalysis.

A theoretical approach is employed to study the catalysis of the dihydroxyacetone phosphate (DHAP) to D-glyceraldehyde 3-phosphate (GAP) reaction by the enzyme triose phosphate isomerase (TIM). The conformational change in a loop involved in protecting the active site from solvent is examined by use of X-ray data and molecular dynamics simulations. A mixed quantum-mechanics and molecular mechanics potential is used to determine the energy surface along the reaction path. The calculations address the role of the enzyme in lowering the barrier to reaction and provide a decomposition into specific residue contributions. To obtain a clearer understanding of the electronic effects, the polarization of the substrate carbonyl group by the active site residues is examined and compared with FTIR measurements on the wild-type and mutant forms of the enzyme.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Free energy perturbation calculations on binding and catalysis after mutating Asn 155 in subtilisin.

Site-directed mutagenesis is a very powerful approach to altering the biological functions of proteins, the structural stability of proteins and the interactions of proteins with other molecules. Several experimental studies in recent years have been directed at estimating the changes in catalytic properties, (rates of binding and catalysis) in site-directed mutants of enzymes compared to the native enzymes. Simulation approaches to the study of complex molecules have also become more powerful, in no small measure owing to the increase in computer power. These simulations have often allowed results of experiments to be rationalized and understood mechanistically. A new approach called the free-energy pertubation method, which uses statistical mechanics and molecular dynamics can often be used for quantitative calculation of free energy differences. We have applied such a technique to calculate the differential free energy of binding and free energy of activation for catalysis of a tripeptide substrate by native subtilisin and a subtilisin mutant (Asn 155----Ala 155). Our studies lead to a calculated difference in free energy of binding which is relatively small, but a calculated change in free energy of catalysis which is substantial. These energies are very close to those determined experimentally (J. A. Wells and D. A. Estell, personal communication), which were not known to us until the simulations were completed. This demonstrates the predictive power and utility of theoretical simulation methods in studies of the effects of site-specific mutagenesis on both enzyme binding and catalysis.

The importance of reactant positioning in enzyme catalysis: a hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase.

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the S(N)2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2-4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. S(N)2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

Structure of the triosephosphate isomerase-phosphoglycolohydroxamate complex: an analogue of the intermediate on the reaction pathway.

The glycolytic enzyme triosephosphate isomerase (TIM) catalyzes the interconversion of the three-carbon sugars dihydroxyacetone phosphate (DHAP) and D-glyceraldehyde 3-phosphate (GAP) at a rate limited by the diffusion of substrate to the enzyme. We have solved the three-dimensional structure of TIM complexed with a reactive intermediate analogue, phosphoglycolohydroxamate (PGH), at 1.9-A resolution and have refined the structure to an R-factor of 18%. Analysis of the refined structure reveals the geometry of the active-site residues and the interactions they make with the inhibitor and, by analogy, the substrates. The structure is consistent with an acid-base mechanism in which the carboxylate of Glu-165 abstracts a proton from carbon while His-95 donates a proton to oxygen to form an enediol (or enediolate) intermediate. The conformation of the bound substrate stereoelectronically favors proton transfer from substrate carbon to the syn orbital of Glu-165. The crystal structure suggests that His-95 is neutral rather than cationic in the ground state and therefore would have to function as an imidazole acid instead of the usual imidazolium. Lys-12 is oriented so as to polarize the substrate oxygens by hydrogen bonding and/or electrostatic interaction, providing stabilization for the charged transition state. Asn-10 may play a similar role.

Computer simulation and analysis of the reaction pathway of triosephosphate isomerase.

A theoretical approach designed for chemical reactions in the condensed phase is used to determine the energy along the reaction path of the enzyme triosephosphate isomerase. The calculations address the role of the enzyme in lowering the barrier to reaction and provide a decomposition into specific residue contributions. The results suggest that, although Lys-12 is most important, many other residues within 16 A of the substrate contribute and that histidine-95 as the imidazole/imidazolate pair could act as an acid/base catalyst.

Simulation of the enzyme reaction mechanism of malate dehydrogenase.

A hybrid numerical method, which employs molecular mechanics to describe the bulk of the solvent-protein matrix and a semiempirical quantum-mechanical treatment for atoms near the reactive site, was utilized to simulate the minimum energy surface and reaction pathway for the interconversion of malate and oxaloacetate catalyzed by the enzyme malate dehydrogenase (MDH). A reaction mechanism for proton and hydride transfers associated with MDH and cofactor nicotinamide adenine dinucleotide (NAD) is deduced from the topology of the calculated energy surface. The proposed mechanism consists of (1) a sequential reaction with proton transfer preceding hydride transfer (malate to oxaloacetate direction), (2) the existence of two transition states with energy barriers of approximately 7 and 15 kcal/mol for the proton and hydride transfers, respectively, and (3) reactant (malate) and product (oxaloacetate) states that are nearly isoenergetic. Simulation analysis of the calculated energy profile shows that solvent effects due to the protein matrix dramatically alter the intrinsic reactivity of the functional groups involved in the MDH reaction, resulting in energetics similar to that found in aqueous solution. An energy decomposition analysis indicates that specific MDH residues (Arg-81, Arg-87, Asn-119, Asp-150, and Arg-153) in the vicinity of the substrate make significant energetic contributions to the stabilization of proton transfer and destabilization of hydride transfer. This suggests that these amino acids play an important role in the catalytic properties of MDH.