The nature of enzyme catalysis in trypsin

Optically computed phase microscopy for quantitative dynamic imaging of label-free cells and nanoparticles

Numerous drug delivery systems based on nanoparticles have been developed, such as those used in BioNTech/Pfizer's and Moderna's Covid vaccines. Knowledge on mechanical interactions between cells and nanoparticles is critical to advance the efficiency and safety of these drug delivery systems. To quantitatively track the motion of cell (transparent) and nanoparticles (nontransparent) with nanometer displacement sensitivity, we investigate a novel imaging technology, optically computed phase microscopy (OCPM) that processes 3D spatial-spectral data through optical computation. We demonstrate that OCPM has the capability to image the motion of cells and magnetic nanoparticles that are mechanically excited by an external magnetic field, quantitatively and in the en face plane.

Bioactive Peptides from Liquid Milk Protein Concentrate by Sequential Tryptic and Microbial Hydrolysis

Recently, bioactive peptides as a health-promoting agent have come to the forefront of health research; however, industrial production is limited, possibly due to the lack of the required technological knowledge. The objective of the investigation was to prepare bioactive peptides with hypoallergenic properties from liquid milk protein concentrate (LMPC), through sequential enzymatic and microbial hydrolysis. LMPC was produced from ultra-heat-treated (UHT) skimmed cow’s milk using a nanofiltration membrane. The effect of the concentration of trypsin (0.008–0.032 g·L−1) on the hydrolysis of LMPC was studied. Subsequently, the hydrolysis of tryptic-hydrolyzed LMPC (LMPC-T) with lactic acid bacteria was performed, and the effect of glucose in microbial hydrolysis was studied. Aquaphotomic analysis of the hydrolysis of LMPC was performed using the spectral range of 1300–1600 nm (near-infrared spectra). Changes in antioxidant capacity, anti-angiotensin-converting enzyme activity, and antibacterial activity against Bacillus cereus, Staphylococcus aureus and Listeria monocytogenes were noted after the sequential tryptic and microbial hydrolysis of LMPC. Allergenicity in LMPC was reduced, due to sequential hydrolysis with 0.016 g·L−1 of trypsin and lacteal acid bacteria. According to the aquaphotomic analysis result, there was a dissociation of hydrogen bonds in compounds during the initial period of fermentation and, subsequently, the formation of compounds with hydrogen bonds. The formation of compounds with a hydrogen bond was more noticeable when microbial hydrolysis was performed with glucose. This may support the belief that the results of the present investigation will be useful to scale up the process in the food and biopharmaceutical industries.

Optimization and Pretreatment for Hot Water Extraction of Korean Deer (Cervus canadensis Erxleben) Velvet Antlers

Velvet antler (VA) is a historically traditional medicinal supplement and is well known in Asian countries for its pharmaceutical and health benefits. The objectives for this study were to optimize the hot water extraction (HWE) of VA for the Korean VA industry, and to determine the most effective pretreatment method among microwave (MW), ultrasonication (US), and enzymatic (EZ) techniques. Using response surface methodology, optimum extraction temperatures and times were determined by central composite design configuration based on extraction yield and sialic acid content. Various quality parameters of VA extract including yield, soluble solid, protein, and sialic acid contents were also compared with the conjunction of HWE and pretreatment. The yield and sialic acid content of VA extract were determined to be 40% and 0.73 mg/g, respectively, under an optimum temperature of 100°C at 24 h of extraction time. The yields from VA extracts pretreated with MW, US, and EZ were 17.42%, 19.73%, and 29.15%, respectively. Among the tested commercial enzymes, pepsin was the most effective proteolytic enzyme and led to the highest yield (47.65%), soluble solids (4.03 °brix), protein (1.12 mg/ml), and sialic acid (3.04 mg/ml) contents from VA extract.

Catalysis by solvation rather than the desolvation effect: exploring the catalytic efficiency of SAM-dependent chlorinase

Chlorinase SalL active sites provide electrostatic stabilization of the transition state which is the origin of its catalytic effect.

Improved Production of Active Streptomyces griseus Trypsin with a Novel Auto-Catalyzed Strategy

N-terminal sequences play crucial roles in regulating expression, translation, activation and enzymatic properties of proteins. To reduce cell toxicity of intracellular trypsin and increase secretory expression, we developed a novel auto-catalyzed strategy to produce recombinant trypsin by engineering the N-terminus of mature Streptomyces griseus trypsin (SGT). The engineered N-terminal peptide of SGT was composed of the thioredoxin, glycine-serine linker, His₆-tag and the partial bovine trypsinogen pro-peptide (DDDDK). Furthermore, we constructed a variant TLEI with insertion of the artificial peptide at N-terminus and site-directed mutagenesis of the autolysis residue R145. In fed-batch fermentation, the production of extracellular trypsin activity was significantly improved to 47.4 ± 1.2 U·ml⁻¹ (amidase activity, 8532 ± 142.2 U·ml⁻¹ BAEE activity) with a productivity of 0.49 U·ml⁻¹·h⁻¹, which was 329% greater than that of parent strain Pichia pastoris GS115-SGT. This work has significant potential to be scaled-up for microbial production of SGT. In addition, the N-terminal peptide engineering strategy can be extended to improve heterologous expression of other toxic enzymes.

Evolution of catalytic antibody repertoire in autoimmune mice

We have attempted to efficiently obtain catalytic antibodies (catAbs) with amidase/esterase activity in the expanded sequence space of the antibody repertoire. In doing so, we used an autoimmune mouse strain, MRL/lpr, that is known to produce enhanced levels of autoantibodies. We applied different types of haptens, such as, and, that are supposed to mimic the transition state of the substrate in the ester/amide hydrolysis. Among them, hapten (2) could not be used, as it was readily broken down after synthesis. Upon immunization with hapten (1), catAbs preferentially evolved in MRL/lpr mice, but this did not happen upon immunization with haptens (3) and (4). Independently, immunization to MRL/lpr mice with successfully elicited the catAbs with the ability to activate vitamin B(6) prodrugs. The common observation seen in these two cases is that most of the catAbs derived from MRL/lpr mice by hapten (1) and half of them by hapten (5) had a Lys at H95, which is at the junctional N region between the V(H) and J(H) gene segments. Despite the conservation of Lys (H95), analyses of the N-region and utilization of the D gene segment in the heavy chain gene showed that these catAbs were from several independent clones of the same family. Studies of site-directed mutagenesis suggest that, in the catAbs elicited from hapten (1), a Lys (H95) and a His (L91) are involved in the catalytic function. Both residues are known to interact with the phosphonate moiety of hapten (1). Such studies also suggest that, in the catAbs elicited from hapten (5), a Lys (H95) and a His (H35) are involved in the catalytic function. These basic amino acids seem to be important for binding to the phosphonate hapten, as they were not changed even after extensive evolution following multiple mutations. By contrast, in normal BALB/c mice, immunization of hapten (1) resulted in eliciting catAbs in lower yield and the majority were the non-catAbs, whose sequences were quite different from those of the catAbs from MRL/lpr mice. They were clonally related to one another and most of them originated from a single clone. The positions of the interacting key residues in the CDRs that interact with the phosphorus moiety strongly differ between our catAbs and other reported catAbs with esterase/amidase activity, which were elicited by the phosphonate/phosphonamidate haptens from normal mice. Further comparison of antibodies elicited by the phosphorus haptens, such as DNA, RNA, phosphocholine, and phosphotyrosine, indicated that none of them had sequence similarity in the basic amino acids and their positions in the CDRs, except for one example, which is anti-DNA antibody elicited from C3H-lpr mice. Analysis based on the classification of canonical structures of the antibodies again suggested that our catAbs derived from MRL/lpr mice belong to an unusual class that is not listed in the literature. Taken together, the above evidence suggests that the unique catalytic subsets that existed in the initial repertoire in the MRL/lpr mice could effectively be captured by the phosphonate haptens through the interaction with the Lys at H95. In the BALB/c mice, however, another noncatalytic subset with an ability to bind only to a moiety other than the phosphonate moiety alternatively evolved, because of the lowest abundance or elimination of the catalytic subsets.

Roles of electrostatic interaction in proteins

More than 60 years after the analyses by Linderstrom-Lang and Kirkwood of their hypothetical 'protein' structures, we have now a plethora of experimental evidence and computational estimates of the electrostatic forces in proteins, with very many protein 3D structures at atomic resolution. In the mean time, there were in the beginning, many arguments and suggestions about the roles of electrostatics, mainly from empirical findings and tendencies. A few experimental results indicated that the electrostatic contribution is of the order of several kcal/mol, which was theoretically difficult to reproduce correctly, because a large opposing reaction field should be subtracted from a large, direct Coulombic field. Although the importance of the reaction field was recognized even 70 years ago, appropriate applications to protein molecules were made only in this decade, with the development of numerical computation. Now, an electrostatic molecular surface is one of the most popular pictures in journals of structural biology, indicating that the electrostatic force is one of the important components contributing to molecular recognition, which is a major focus of current biology and biochemistry. The development of NMR techniques has made it possible to observe the individual ionizations of ionizable groups in a protein, in addition to the determination of the 3D structure. Since it does not require any additional probe, each charge state can report the very local and heterogeneous electrostatic potentials working in the protein, without disturbing the original field. From the pKa values, the contributions of electrostatic interactions, ion pairs, charge-dipole interactions, and hydrogen bonds to protein stability have been correctly evaluated. Protein engineering also provides much more information than that obtainable from the native proteins, as the residues concerned can now be easily substituted with other amino acid residues having electrostatically different characteristics. Those experimental results have revealed smaller contributions than previously expected, probably because we underestimated the reaction field effects. Especially, a single ion pair stabilizes a protein only slightly, although a cooperative salt-bridge network can contribute significantly to protein stability. Marginal stabilities of proteins arise from small difference between many factors with driving and opposing forces. In spite of the small contribution of each single electrostatic interaction to the protein stability, the sum of their actions works to maintain the specific 3D structure of the protein. The 'negative' roles of electrostatics, which might destabilize protein conformation, should be pointed out. Unpaired buried charges are energetically too expensive to exit in the hydrophobic core. Isolated hydrogen bond donors and acceptors also exert negative effects, but they are not as expensive as the unpaired buried charges, with costs of a few kcal/mol. Therefore, statistical analyses of protein 3D structures reveal only rare instances of isolated hydrogen bond donors and acceptors. This must be the main reason why alpha-helices and beta-sheets are only observed in protein cores as the backbone structures. Such secondary structures do not leave any backbone hydrogen donors or acceptors unpaired, because of their intrinsically regular packing. Otherwise, it might be very difficult to construct a backbone structure, in which all the backbone amide and carbonyl groups had their own hydrogen bond partners in the protein core. There are two theoretical approaches to protein electrostatics, the macroscopic or continuum model, and the microscopic or molecular model. As described in this article, the macroscopic model has inherent problems because the protein-solvent system is very hetergeneous from the physical point of view...

Model assembly study of the ligand binding by p-hydroxybenzoate hydroxylase: correlation between the calculated binding energies and the experimental dissociation constants

The energies of binding of seven ligands by p-hydroxybenzoate hydroxylase (PHBH) were calculated theoretically. Direct enzyme-ligand interaction energies were calculated using the ab initio quantum mechanical model assembly of the active site at the 3-21G level. Solvation energies of the ligands needed in the evaluation of the binding energies were calculated with the semiempirical AM1-SM2 method and the long-range electrostatic interaction energies between the ligands and the protein matrix classically using the static charge distributions of the ligands and the protein. Energies for proton-transfer between the ligands' OH or SH substituent at position 4 and the active-site tyrosine within the ab initio model assemblies were calculated and compared to the corresponding pKas in aqueous solution. Excluding 3,4-dihydroxybenzoate, the natural product of PHBH, a linear relationship between the calculated binding energies and the experimental binding free energies was found with a correlation coefficient of 0.90. Contributions of the direct enzyme-ligand interaction energies, solvation energies and the long-range electrostatic interaction energies to the calculated binding energies were analyzed. The proton-transfer energies of the ligands with substituents ortho to the ionized OH were found to be perturbed less in the model calculations than the energies of their meta isomers as deduced from the corresponding pKas.

GREEN: a program package for docking studies in rational drug design

A program package, GREEN, has been developed that enables docking studies between ligand molecules and a protein molecule. Based on the structure of the protein molecule, the physical and chemical environment of the ligand-binding site is expressed as three-dimensional grid-point data. The grid-point data are used for the real-time evaluation of the protein-ligand interaction energy, as well as for the graphical representation of the binding-site environment. The interactive docking operation is facilitated by various built-in functions, such as energy minimization, energy contribution analysis and logging of the manipulation trajectory. Interactive modeling functions are incorporated for designing new ligand molecules while considering the binding-site environment and the protein-ligand interaction. As an example of the application of GREEN, a docking study is presented on the complex between trypsin and a synthetic trypsin inhibitor. The program package will be useful for rational drug design, based on the 3D structure of the target protein.

Effect of specificity on ligand conformation in acyl-chymotrypsins

I.r. difference spectroscopy combined with 13C and 18O double-isotope substitution was used to examine the ester acyl carbonyl stretching vibration of hydrocinnamoyl-chymotrypsin. A single acyl carbonyl stretching band was observed at 1731 cm-1. This contrasts with previous i.r. and resonance Raman spectroscopic studies of a number of trans-3-arylacryloyl-chymotrypsins which showed two acyl carbonyl stretching bands in the region of 1700 cm-1, which were proposed to represent productive and non-productive conformations of the acyl-enzyme. The single acyl carbonyl band for hydrocinnamoyl-chymotrypsin suggests only a single conformation, and the comparatively high frequency of this band implies little or no hydrogen-bonding to this carbonyl group. Enzymic hydrogen-bonding to the acyl carbonyl is believed to give bond polarization and thereby catalytic-rate acceleration. Thus, in view of the apparent lack of such hydrogen-bonding in hydrocinnamoyl-chymotrypsin, it should be the case that this acyl-chymotrypsin is less specific than trans-3-arylacryloyl-chymotrypsins, whereas the opposite is true. It is therefore proposed that there may be a productive acyl carbonyl population of lower stretching frequency for hydrocinnamoyl-chymotrypsin, but that this is too small to be discerned because of either a relatively high deacylation rate or an unfavourable conformational equilibrium. The single acyl carbonyl band for hydrocinnamoyl-chymotrypsin is significantly broader than those for trans-3-arylacryloyl-chymotrypsins, indicating that this group is more conformationally mobile and dispersed in the former. This can be correlated with the absence of acyl carbonyl hydrogen-bonding in hydrocinnamoyl-chymotrypsin, and with the much greater flexibility of the saturated hydrocinnamoyl group than unsaturated trans-3-arylacryloyl. This flexibility is presumably the reason why hydrocinnamoyl-chymotrypsin is more specific than trans-3-arylacryloyl-chymotrypsins. Resonance Raman spectroscopy is limited to the non-specific trans-3-arylacryloyl-chymotrypsins because of its chromophoric requirement, whereas i.r. may be used to examine non-chromophoric more specific acyl-enzymes such as hydrocinnamoyl-chymotrypsin. The results presented in this paper suggest that trans-3-arylacryloyl-chymotrypsins are atypical.

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.

Active site dynamics of acyl-chymotrypsin

The motions of water molecules, the acyl moiety, the catalytic triad, and the oxyanion binding site of acyl-chymotrypsin were studied by means of a stochastic boundary molecular dynamics simulation. A water molecule that could provide the nucleophilic OH- for the deacylation stage of the catalysis was found to be trapped between the imidazole ring of His-57 and the carbonyl carbon of the acyl group. It makes a hydrogen bond with the N epsilon 2 of His-57 and is held in place through a network of hydrogen-bonded water molecules in the active site. The water molecule was found as close as 2.8 A to the carbonyl carbon. This appears to be due to the constraints imposed by nonbonded interaction in the active site. Configurations were found in which one hydrogen of the trapped water shared a bifurcated hydrogen bond with His-57-N epsilon 2 and Ser-195-O gamma, with the water oxygen very close to the carbonyl carbon. The existence of such a water molecule suggests that large movement of the His-57 imidazole ring between positions suitable for providing general-base catalyzed assistance and for providing general-acid catalyzed assistance may not be required during the reaction. The simulation indicates that the side chains of residues involved in catalysis (i.e., His-57, Ser-195, and Asp-102) are significantly less flexible than other side chains in the protein. The 40% reduction in rms fluctuations is consistent with a comparable reduction calculated from the temperature factors obtained in the X-ray crystallographic data of gamma-chymotrypsin. The greater rigidity of active site residues seems to result from interconnected hydrogen bonding networks among the residues and between the residues and the solvent water in the active site.

Reproduction of correct electrostatic field by charges and dipoles on a closed surface

A general algorithm based on the Green function theorem has been developed to correctly reproduce electrostatic fields inside a closed space by point charges and point dipoles on the surface surrounding the space. For actual computations, limited numbers of point charges, including charge pairs replacing point dipoles, are enough to approximate the inner fields. As examples, reaction fields were reproduced by the current surface charges and dipoles for the dielectric models, where a monopole, dipole, or quadrupole was individually set at the center in a vacuum sphere surrounded by high dielectric continuum. The potentials due to those reaction fields agree well with the analytical ones. As an application of this method to the analysis of the electronic structure of the active site of a protein, a combination of the continuum dielectric model and ab initio molecular orbital calculation was carried out. Other applications to molecular dynamics and quantum mechanical calculations are also discussed.

A priori evaluation of aqueous polarization effects through Monte Carlo QM-MM simulations

A Monte Carlo quantum mechanical-molecular mechanical (QM-MM) simulation method was used to determine the contributions of the solvent polarization effect to the total interaction energies between solute and solvent for amino acid side chains and nucleotide bases in aqueous solution. In the present AM1-TIP3P approach, the solute molecule is characterized by valence electrons and nucleus cores with Hartree-Fock theory incorporating explicit solvent effects into the total Hamiltonian, while the solvent is approximated by the three-point charge TIP3P model. The polarization energy contributes 10 to 20 percent of the total electrostatic energy in these systems. The performance of the hybrid AM1-TIP3P model was further validated by consideration of bimolecular complexes with water and by computation of the free energies of solvation of organic molecules using statistical perturbation theory. Excellent agreement with ab initio 6-31G(d) results and experimental solvation free energies was obtained.

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.

Possible coupling of chemical to structural dynamics in subtilisin BPN' catalyzed hydrolysis

The viscosity dependence of enzymatic catalysis was examined in subtilisin BPN' catalyzed hydrolysis of N-succinyl-Ala-Ala-Pro-Phe-p-nitroanilide and N-succinyl-Ala-Ala-Pro-Phe-thiobenzyl ester. The viscosity of the reaction medium was varied by added glycerol, ethylene glycol, sucrose, glucose, fructose, poly(ethylene glycol) and Ficoll-400. Responses of the Michaelis-Menten parameters associated with hydrolysis were calculated from data obtained by spectrophotometric techniques. The reactions with these two substrates have catalytic rates well below the diffusion-controlled limit and thus enable us to study the viscosity effects on catalytic steps of non-transport nature. It was found that the Km values for both amide and ester reactions remained relatively independent of cosolvents. On the other hand, while the kcat values for amide were insensitive to cosolvents, those for ester were substantially attenuated except in the case of poly(ethylene glycol). The observed rate attenuations cannot be explained by changes in proton activity, water activity, dielectric constant of the reaction medium or shifts of any kinetically important pKa. Instead, the results can be adequately described by microviscosity effects on the unimolecular deacylation step with a coupling constant of 0.65 +/- 0.11. In addition, the different viscosity dependence in the acylation vs deacylation step can be rationalized in terms of fluctuation-dependent chemical dynamics of proton transfers in the context of the Bogris-Hynes model.

Hydrogen-bonding in enzyme catalysis. Fourier-transform infrared detection of ground-state electronic strain in acyl-chymotrypsins and analysis of the kinetic consequences

I.r. difference spectra are presented for 3-(indol-3-yl)acryloyl-, cinnamoyl-, 3-(5-methylthien-2-yl)acryloyl-, dehydrocinnamoyl- and dihydrocinnamoyl-chymotrypsins at low pH, where the acyl-enzymes are catalytically inactive. At least two absorption bands are seen in each case in the ester carbonyl stretching region of the spectrum. Cinnamoyl-chymotrypsin substituted at the carbonyl carbon atom with 13C was prepared. A difference spectrum in which 13C-substituted acyl-enzyme was subtracted from [12C]acyl-enzyme shows two bands in the ester carbonyl region and thus confirms the assignment of the features to the single ester carbonyl group. The frequencies of the ester carbonyl bands are interpreted in terms of differential hydrogen-bonding. In each case a lower-frequency relatively narrow band is assigned to a productive potentially reactive binding mode in which the carbonyl oxygen atom is inserted in the oxyanion hole of the enzyme active centre. The higher-frequency band, which is broader, is assigned to a non-productive binding mode in each case, where a water molecule bridges from the carbonyl oxygen atom to His-57; this mode is equivalent to the crystallographically determined structure of 3-(indol-3-yl)acryloyl-chymotrypsin, i.e. the Henderson structure. A difference spectrum of dihydrocinnamoyl-chymotrypsin taken at higher pH shows resolution of a feature centred upon 1731 cm-1, which is assigned to a non-bonded conformer in which the carbonyl oxygen atom is not hydrogen-bonded. Perturbation of the protein spectrum in the presence of acyl groups is interpreted in terms of enhanced structural rigidity. It is reported that the ester carbonyl region of the difference spectrum of cinnamoyl-subtilisin is complicated by overlap of features that arise from protein perturbation. Measurements of carbonyl absorption frequencies in a number of solvents of the methyl esters of the acyl groups used to make acyl-enzymes have permitted determination of the apparent dielectric constants experienced by carbonyl groups in the enzyme active centre as well as a discussion of the effects of polarity. The ester carbonyl bond strengths of the various conformations were estimated by using simple harmonic oscillator theory and an empirical relation between the force constants and bond strengths. The fractional bond breaking induced by hydrogen-bonding was used to calculate rate enhancement factors by using absolute reaction rate theory.(ABSTRACT TRUNCATED AT 400 WORDS)

Simulations of dynamical properties of a Michaelis complex: porcine pancreatic elastase and the hexapeptide, Thr-Pro-n Val-Leu-Tyr-Thr

Molecular dynamic simulations (30ps) of the Michaelis complex of hexapeptide (Thr-Pro-nVal-Leu-Tyr-Thr) bound to porcine pancreatic elastase (PPE) hydrated by about 2000 water molecules have been performed using the AMBER 3.0 program package. Dynamical properties of the conformation of the active site have been examined. A comparison with previously reported simulations of native PPE shows that after the substrate is bound, the catalytically crucial H-bond between O gamma-H group of (Ser 195) and nitrogen N epsilon (His 57) is more readily formed. These results show, however, that the H-bond does not adopt the most favorable conformation. The O gamma-H group of Ser 195 has a statistical preference for an attractive interaction with the O = C carbonyl (Ser 214) rather than the nitrogen N epsilon (His 57).

Dynamic properties of the first enzymatic reaction steps of porcine pancreatic elastase. How rigid is the active site of the native enzyme? Molecular dynamics simulation

Two molecular dynamics simulations (100 and 50 ps) of native porcine pancreatic elastase i.e., without bound substrate and with the active site hydrated by a dome of water (630 molecules) have been performed. Dynamical properties of the catalytic tetrad have been examined. While relative conformations of the Asp 102, His 57, and Ser 214 are rather stable in time, the side chain of Ser 195 undergoes several conformational changes. No preferences are observed for the formation of a hydrogen bond between the O gamma-H group (Ser 195) and nitrogen N, (His 57). A cluster of ordered water molecules effectively competes with the H-O gamma group (Ser 195) and thereby prevents the formation of this H bond, which is generally agreed to be crucial for catalysis.

Enzymes work by solvation substitution rather than by desolvation

Considerable attention has recently been drawn to the hypothesis that enzymes catalyze their reactions by displacing solvent and creating an environment similar to the gas phase for the reacting substrates. This "desolvation hypothesis" is reexamined in this paper by defining a common reference energy for reactions in various environments. It is argued that consistent attempts to describe the actual energetics of enzymatic reactions, taking either gas phase or solution as a reference, would contradict the above hypothesis. That is, the enzyme does remove water molecules from its substrate, but substitutes these molecules for another polar environment (namely, its active site). By taking amide hydrolysis as an example, we use experimentally estimated solvation energies and analyze the reaction profile in the gas phase, in solution, and in enzyme active sites. We show that the gas-phase reaction is characterized by an enormous activation barrier (associated with forming the charged nucleophile from neutral fragments), although the nucleophilic attack is essentially barrierless. On the other hand, the enzyme and solution reactions are found to have similar reaction profiles, with a lower activation barrier for the enzymatic reaction. Presumably, the fact that previous analyses of this problem did not involve the construction of the relevant thermodynamic cycles (and quantitative estimates of the corresponding solvation energies) led to the desolvation hypothesis. Our conclusion is that enzyme active sites provide specific polar environments that do not resemble the gas phase but that are designed for electrostatic stabilization of ionic transition states and that "solvate" these states more than water does.

Calculations of free energy profiles for the staphylococcal nuclease catalyzed reaction

Calculations of the free energy profile for the first two (rate-limiting) steps of the staphylococcal nuclease catalyzed reaction are reported. The calculations are based on the empirical valence bond method in combination with free energy perturbation molecular dynamics simulations. The calculated activation free energy is in good agreement with experimental kinetic data, and the catalytic effect of the enzyme is reproduced without any arbitrary adjustment of parameters. The enormous reduction of the activation barrier (relative to the reference reaction in water) appears to be largely associated with the strong electrostatic effect of the Ca2+ ion and the two arginine residues in the active site. This favorable electrostatic environment reduces the cost of the general-base catalysis step by almost 15 kcal/mol (by stabilizing the OH- nucleophile) and then stabilizes the developing negative charge on the 5'-phosphate group in the second step of the reaction by about 19 kcal/mol. The basic features of the originally postulated enzyme mechanism (Cotton et al., 1979) are found to be compatible with the observed activation free energy. However, the proposed modification of the mechanism (Sepersu et al., 1987), in which Arg 87 interacts only with the pentacoordinated transition state, is supported by the simulations. Further calculations on the D21E mutant also give results in good agreement with kinetic data.

Mechanism of the chain extension step in the biosynthesis of fatty acids

The chain extension step in the enzymatic synthesis of fatty acids by fatty acid synthase, involving a formal Claisen condensation of thio esters, has been clarified by theoretical calculations for model systems, using the modified neglect of diatomic overlap and Austin Model 1 parametric self-consistent field molecular orbital procedures. The reaction involves a free carbanion, formed by decarboxylation of a malonate ion. Formation of the carbanion and condensation with the fatty acid thio ester are not concerted. The decarboxylation is strongly endothermic. It is brought about by electrostatic interaction (field effect) with an ammonium ion derived from an adjacent lysine residue, the ions being far enough apart to inhibit proton transfer between them. Proton transfer would lead to an enol that is predicted not to be able to undergo the Claisen condensation. The formation of the ammonium ion is considered in terms of the pKa of the relevant groups. The bearing of this work on a recent interpretation of the activity and selectivity of enzyme reactions is discussed, and some misunderstandings concerning this interpretation are clarified.

Modeling complex molecular interactions involving proteins and DNA

We have presented a perspective of progress in three areas of simulations of complex molecules: the development of force fields for molecular simulation; the application of computer graphics, molecular mechanics and molecular dynamics in simulations of DNA and DNA-drug complexes and the application of computer graphics, molecular mechanics and quantum mechanics in studies of enzyme substrate interactions. It is our perspective that improvements are being made in force fields, and these will allow a more accurate simulation of structures and energies of complex molecules. In the area of DNA molecular mechanics and dynamics, it is clear that the use of computer graphics model building combined with NMR NOE data is a potentially very powerful tool in accurately determining structures of drug-DNA complexes using molecular mechanics and dynamics. Finally, we are in a position to reasonably simulate structures and (qualitatively) energies for complete reaction pathways of enzymes using a combination of computer graphics, molecular mechanics and quantum mechanics. More accurate energies and pathways are sure to follow, using the combined molecular mechanics/quantum mechanics optimization developed by Singh and the free energy perturbation methods pioneered in Groningen and Houston.