Molecular modeling and its experimental verification for the catalytic mechanism of Candida antarctica lipase B

Quantum mechanical and molecular dynamics simulation analysis has been performed on the model system for CALB (Candida antarctica lipase B) with esters to study the reaction mechanism and conformational preference of catalytic hydrolysis and the esterification reaction. Using quantum mechanical analysis, the ping-pong bi-bi mechanism was applied and energies and 3-dimensional binding configurations of the whole reaction pathways were calculated. Further molecular dynamics simulation analysis was performed on the basis of the transition state obtained from quantum mechanical study to observe the effect of structures of,the substrates. Calculation results using substrates of different chain length and chiral configurations were compared for conformational preference. The calculated results showed very small influence on chain length, whereas chiral conformation showed big differences. Calculated results from molecular modeling studies have been compared qualitatively with the experimental data using racemic mixtures of (+/-)-cis-4-acetamido-cyclopent-2-ene-1-ethyl acetate as substrates.     

Refinement of the NMR structures for acyl carrier protein with scalar coupling data

Structure determination of small proteins using NMR data is most commonly pursued by combining NOE derived distance constraints with inherent constraints based on chemical bonding. Ideally, one would make use of a variety of experimental observations, not just distance constraints. Here, coupling constant constraints have been added to molecular mechanics and molecular dynamics protocols for structure determination in the form of a psuedoenergy function that is minimized in a search for an optimum molecular conformation. Application is made to refinement of a structure for a 77 amino acid protein involved in fatty acid synthesis, Escherichia coli acyl carrier protein (ACP). 54 3JHN alpha coupling constants, 12 coupling constants for stereospecifically assigned side chain protons, and 450 NOE distance constraints were used to calculate the 3-D structure of ACP. A three-step protocol for a molecular dynamics calculation is described, in analogy to the protocol previously used in molecular mechanics calculations. The structures calculated with the molecular mechanics approach and the molecular dynamics approach using a rigid model for the protein show similar molecular energies and similar agreement with experimental distance and coupling constant constraints. The molecular dynamics approach shows some advantage in overcoming local minimum problems, but only when a two-state averaging model for the protein was used, did molecular energies drop significantly.     

Modeling of substrate and inhibitor binding to phospholipase A2.

Molecular graphics and molecular mechanics techniques have been used to study the mode of ligand binding and mechanism of action of the enzyme phospholipase A2. A substrate-enzyme complex was constructed based on the crystal structure of the apoenzyme. The complex was minimized to relieve initial strain, and the structural and energetic features of the resultant complex analyzed in detail, at the molecular and residue level. The minimized complex was then used as a basis for examining the action of the enzyme on modified substrates, binding of inhibitors to the enzyme, and possible reaction intermediate complexes. The model is compatible with the suggested mechanism of hydrolysis and with experimental data about stereoselectivity, efficiency of hydrolysis of modified substrates, and inhibitor potency. In conclusion, the model can be used as a tool in evaluating new ligands as possible substrates and in the rational design of inhibitors, for the therapeutic treatment of diseases such as rheumatoid arthritis, atherosclerosis, and asthma.     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

The action of alpha chymotrypsin upon mixtures of esters and proteins at enzyme-saturating concentrations

The behavior of α-chymotrypsin has been studied in the simultaneous presence of two different substrates, each present in the reaction mixture at its saturation level. Mixtures of two esters were hydrolyzed at rates intermediate between the rates of hydrolysis of each ester when present alone, suggesting, in this case, competitive hydrolysis. In contrast, the rates of hydrolysis in mixtures of casein with gelatin or of either protein with an ester were equal to the sum of the rates of hydrolysis of the separate substrates, indicating in these cases independent hydrolysis. The activity of the α-chymotrypsin preparation used could not be attributed to contamination with other enzymes. Studies of the effect of soy bean inhibitor on chymotrypsin indicate that the mechanism of inhibition with protein substrates differs from that when esters are used, providing further evidence that α-chymotrypsin reacts differently with esters and proteins. These results indicate that if chymotrypsin forms specific complexes with its substrates, it must possess at least three distinct active sites. However there is independent chemical evidence that the proteolytic and esterolytic activities of this enzyme reside in the same active center. If this is true, the experimental observations reported here cannot be explained unless it is supposed that this enzyme does not form specific Michaelis complexes with its substrates.     

Determinants of reactivity and selectivity in soluble epoxide hydrolase from quantum mechanics/molecular mechanics modeling

Soluble epoxide hydrolase (sEH) is an enzyme involved in drug metabolism that catalyzes the hydrolysis of epoxides to form their corresponding diols. sEH has a broad substrate range and shows high regio- and enantioselectivity for nucleophilic ring opening by Asp333. Epoxide hydrolases therefore have potential synthetic applications. We have used combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations (at the AM1/CHARMM22 level) and high-level ab initio (SCS-MP2) QM/MM calculations to analyze the reactions, and determinants of selectivity, for two substrates: trans-stilbene oxide (t-SO) and trans-diphenylpropene oxide (t-DPPO). The calculated free energy barriers from the QM/MM (AM1/CHARMM22) umbrella sampling MD simulations show a lower barrier for phenyl attack in t-DPPO, compared with that for benzylic attack, in agreement with experiment. Activation barriers in agreement with experimental rate constants are obtained only with the highest level of QM theory (SCS-MP2) used. Our results show that the selectivity of the ring-opening reaction is influenced by several factors, including proximity to the nucleophile, electronic stabilization of the transition state, and hydrogen bonding to two active site tyrosine residues. The protonation state of His523 during nucleophilic attack has also been investigated, and our results show that the protonated form is most consistent with experimental findings. The work presented here illustrates how determinants of selectivity can be identified from QM/MM simulations. These insights may also provide useful information for the design of novel catalysts for use in the synthesis of enantiopure compounds.     

Comparative studies of the specificities of α-chymotrypsin and subtilisin BPN′. Studies with flexible and `locked'' substrates

Subtilisin BPN' hydrolysed N-acetyl-l-3-(2-naphthyl)-alanine methyl ester, N-acetyl-l-leucine methyl ester and N-acetyl-l-valine methyl ester, faster than alpha-chymotrypsin. Of eight ;locked' substrates tested, only methyl 5,6-benzindan-2-carboxylate was hydrolysed faster by subtilisin, whereas the other esters were better substrates for chymotrypsin. Compared with the values for chymotrypsin, the stereospecific ratios during the hydrolysis of the optically active locked substrates by subtilisin were decreased by one and two orders of magnitude for bi- and tri-cyclic substrates respectively. The polar groups adjacent to the alpha-carbon atom of locked substrates did not contribute significantly to the reactivity of the more active optical isomers, but had a detrimental effect on the less active antipodes during hydrolysis by both the enzymes. These studies show that the binding site of subtilisin BPN' is longer and broader than that of alpha-chymotrypsin.     

A steady-state kinetic model of butyrylcholinesterase from horse plasma

The steady-state kinetics of the butyrylcholinesterase-catalysed hydrolysis of butyrylthiocholine and thiophenyl acetate were shown to deviate from Michaelis–Menten kinetics. The `best' empirical rate law was selected by fitting different rate equations to the experimental data by non-linear regression methods. The results were analysed in view of two alternative interpretations: (1) the reaction is catalysed by a mixture of enzymes, or (2) the activity is due to a single enzyme displaying deviations from Michaelis–Menten kinetics. It was concluded that the second alternative applies, and this conclusion was further supported by experiments involving simultaneous hydrolysis of alternative thiol ester substrates (butyrylthiocholine/thiophenyl acetate) as well as alternative thiol ester and oxygen ester substrates (butyrylthiocholine/phenyl acetate; thiophenyl acetate/butyrylcholine; acetylthiocholine/phenyl acetate). On the basis of the conclusion that a single enzyme is responsible for the activity, a molecular model is proposed. This model involves an acylated enzyme, and implies binding to the enzyme of one acyl group and one ester molecule, but not two ester molecules at the same time. Thus butyrylcholinesterase, which is structurally a tetramer, behaves functionally as a co-operative dimer, an interpretation in accordance with available data from active-site titrations.     

Steady-state and pre-steady-state kinetic evaluation of severe acute respiratory syndrome coronavirus (SARS-CoV) 3CLpro cysteine protease: development of an ion-pair model for catalysis

Severe acute respiratory syndrome (SARS) was a worldwide epidemic caused by a coronavirus that has a cysteine protease (3CLpro) essential to its life cycle. Steady-state and pre-steady-state kinetic methods were used with highly active 3CLpro to characterize the reaction mechanism. We show that 3CLpro has mechanistic features common and disparate to the archetypical proteases papain and chymotrypsin. The kinetic mechanism for 3CLpro-mediated ester hydrolysis, including the individual rate constants, is consistent with a simple double displacement mechanism. The pre-steady-state burst rate was independent of ester substrate concentration indicating a high commitment to catalysis. When homologous peptidic amide and ester substrates were compared, a series of interesting observations emerged. Despite a 2000-fold difference in nonenzymatic reactivity, highly related amide and ester substrates were found to have similar kinetic parameters in both the steady-state and pre-steady-state. Steady-state solvent isotope effect (SIE) studies showed an inverse SIE for the amide but not ester substrates. Evaluation of the SIE in the pre-steady-state revealed normal SIEs for both amide and ester burst rates. Proton inventory (PI) studies on amide peptide hydrolysis were consistent with two proton-transfer reactions in the transition state while the ester data was consistent with a single proton-transfer reaction. Finally, the pH-inactivation profile of 3CLpro with iodoacetamide is indicative of an ion-pair mechanism. Taken together, the data are consistent with a 3CLpro mechanism that utilizes an "electrostatic" trigger to initiate the acylation reaction, a cysteine-histidine catalytic dyad ion pair, an enzyme-facilitated release of P1, and a general base-catalyzed deacylation reaction.     

The retinal conformation and its environment in rhodopsin in light of a new 2.2 A crystal structure

A new high-resolution structure is reported for bovine rhodopsin, the visual pigment in rod photoreceptor cells. Substantial improvement of the resolution limit to 2.2 A has been achieved by new crystallization conditions, which also reduce significantly the probability of merohedral twinning in the crystals. The new structure completely resolves the polypeptide chain and provides further details of the chromophore binding site including the configuration about the C6-C7 single bond of the 11-cis-retinal Schiff base. Based on both an earlier structure and the new improved model of the protein, a theoretical study of the chromophore geometry has been carried out using combined quantum mechanics/force field molecular dynamics. The consistency between the experimental and calculated chromophore structures is found to be significantly improved for the 2.2 A model, including the angle of the negatively twisted 6-s-cis-bond. Importantly, the new crystal structure refinement reveals significant negative pre-twist of the C11-C12 double bond and this is also supported by the theoretical calculation although the latter converges to a smaller value. Bond alternation along the unsaturated chain is significant, but weaker in the calculated structure than the one obtained from the X-ray data. Other differences between the experimental and theoretical structures in the chromophore binding site are discussed with respect to the unique spectral properties and excited state reactivity of the chromophore.     

Sensitivity of the folding/unfolding transition state ensemble of chymotrypsin inhibitor 2 to changes in temperature and solvent

To better characterize the transition state for folding/unfolding and its sensitivity to environmental changes, we have run multiple molecular dynamics simulations of chymotrypsin inhibitor 2 (CI2) under varying solvent conditions and temperature. The transition state structures agree well with experiment, and are similar under all of the conditions investigated here. Increasing the temperature leads to some movement in the position of the transition state along several reaction coordinates, as measured by changes in properties of the transition state structures. These structural changes are in the direction of a more native-like transition state as denaturation conditions become more severe, as expected for a Hammond effect. These structural changes are not, however, reflected in the global structure as measured by the total number of contacts or the average S-values. These results suggest that the small changes in average Phi-values with temperature seen by experiment may be due to an increase in the sensitivity of the transition state to mutation rather than a change in the average structure of the transition state. A simple analysis of the rates of unfolding indicates that the free energy barrier to unfolding decreases with increasing temperature, but even in our very high temperature simulations there is a small free energy barrier.     

Role of the ribose residue of substrates in reactions catalyzed by ribonuclease

The rate constants and Km for the hydrolysis of the optically active nonglycosidic analogues of the CpA and C greater than p catalysed by RNase A and RNase BS-I were measured. The rate of hydrolysis of the model substrates in 10(5) and 10(3) slower that for the appropriate dinucleoside phosphate and nucleoside cyclophosphate. However, substitution of the relatively rigid ribofuranose ring with flexible alifatic chains is accompanied by little variation in binding constants. The analyses based on the single substrate system indicate that the observed difference in rate constants must be accounted for by a difference between the binding of the substrates in the transition state to the RNase active site. Consequently, the "rigidity" of the ribose rings in RNA leads to large decreases in the free energy of activation for the reactions catalysed by RNases.     

Free energy calculations on binding and catalysis by alpha-lytic protease: the role of substrate size in the P1 pocket

We present free energy calculations using molecular dynamics on different substrates of alpha-lytic protease in the gas phase, in solution, while forming a noncovalent Michaelis complex with the enzyme, and in a tetrahedral structure representing a transition state/intermediate for acylation by the enzyme. Various P1 substrates were studied, with P1 = Gly, Ala, Val, and Leu. In qualitative agreement with experiment, the enzyme was calculated to bind and catalyze most effectively substrates with P1 = Ala over those with P1 = Gly, Val or Leu. Also, the calculated relative solvation free energies of Gly----Ala and Ala----Val were in qualitative agreement with experimental values in corresponding model systems. However, the level of quantitative agreement with experiment achieved in our earlier study of relative binding and catalysis of native subtilisin and an Asn-155----Ala mutant was not achieved. We surmise that this is due to the greater difficulty in quantitatively simulating effects that are predominantly van der Waals and hydrophobic compared to those that are hydrogen bonding/electrostatic.     

Mapping the interactions present in the transition state for unfolding/folding of FKBP12.

The structure of the transition state for folding/unfolding of the immunophilin FKBP12 has been characterised using a combination of protein engineering techniques, unfolding kinetics, and molecular dynamics simulations. A total of 34 mutations were made at sites throughout the protein to probe the extent of secondary and tertiary structure in the transition state. The transition state for folding is compact compared with the unfolded state, with an approximately 30 % increase in the native solvent-accessible surface area. All of the interactions are substantially weaker in the transition state, as probed by both experiment and molecular dynamics simulations. In contrast to some other proteins of this size, no element of structure is fully formed in the transition state; instead, the transition state is similar to that found for smaller, single-domain proteins, such as chymotrypsin inhibitor 2 and the SH3 domain from alpha-spectrin. For FKBP12, the central three strands of the beta-sheet, beta-strand 2, beta-strand 4 and beta-strand 5, comprise the most structured region of the transition state. In particular Val101, which is one of the most highly buried residues and located in the middle of the central beta-strand, makes approximately 60 % of its native interactions. The outer beta-strands and the ends of the central beta-strands are formed to a lesser degree. The short alpha-helix is largely unstructured in the transition state, as are the loops. The data are consistent with a nucleation-condensation model of folding, the nucleus of which is formed by side-chains within beta-strands 2, 4 and 5, and the C terminus of the alpha-helix. The precise residues involved in the nucleus differ in the two simulated transition state ensembles, but the interacting regions of the protein are conserved. These residues are distant in the primary sequence, demonstrating the importance of tertiary interactions in the transition state. The two independently derived transition state ensembles are structurally similar, which is consistent with a Bronsted analysis confirming that the transition state is an ensemble of states close in structure.     

Segmented Transition Pathway of the Signaling Protein Nitrogen Regulatory Protein C

Recent advances in experimental methods provide increasing evidence that proteins sample the conformational substates that are important for function in the absence of their ligands. An example is the receiver domain of nitrogen regulatory protein C, a member of the phosphorylation-mediated signaling family of “two-component systems.” The receiver domain of nitrogen regulatory protein C samples both inactive conformation and the active conformation before phosphorylation. Here we determine a possible pathway of interconversion between the active state and the inactive state by targeted molecular dynamics simulations and quasi-harmonic analysis; these methods are used because the experimental conversion rate is in the high microsecond range, longer than those that are easily accessible to atomistic molecular dynamics simulations. The calculated pathway is found to be composed of four consecutive stages described by different progress variables. The lowest quasi-harmonic principal components from unbiased molecular dynamics simulations on the active state correspond to the first stage, but not to the subsequent stages of the transition. The targeted molecular dynamics pathway suggests that several transient nonnative hydrogen bonds may facilitate the transition.     

Molecular dynamics simulations from putative transition states of alpha-spectrin SH3 domain

A series of molecular dynamics simulations in explicit solvent were started from nine structural models of the transition state of the SH3 domain of alpha-spectrin, which were generated by Lindorff-Larsen et al. (Nat Struct Mol Biol 2004;11:443-449) using molecular dynamics simulations in which experimental Phi - values were incorporated as restraints. Two of the nine models were simulated 10 times for 200 ns and the remaining models simulated two times for 200 ns. Complete folding was observed in one case, while in the other simulations partial folding or unfolding events were observed, which were characterized by a regularization of elements of secondary structure. These results are consistent with recent experimental evidence that the folding of SH3 domains involves low populated intermediate states.     

Transition state affinity jump chromatography. A double selection method for isolating catalytically active enzymes and other molecules

A double selection method for isolating active enzyme molecules, using substrate analog affinity chromatography and elution with transition state analogs, is described. To demonstrate the principle, a mixture containing native chymotrypsin and [3H]deoxychymotrypsin, in which the active site serine had been converted to [3H]alanine, was applied to a column containing immobilized D-tryptophan methyl ester. Both forms of chymotrypsin were retained. Catalytically active enzyme was selectively desorbed with the peptide aldehyde chymostatin, leaving catalytically inactive deoxychymotrypsin bound to the substrate analog affinity column. This affinity technique may afford a simple and general method for separating enzymes and other catalysts according to their molecular turnover numbers.     

Simulations of the p97 complex suggest novel conformational states of hydrolysis intermediates

The vitally important AAA (ATPases associated with various cellular activities) protein p97 is involved in cellular functions ranging from replication to degradation of misfolded proteins and has recently been proposed as a novel chemotherapeutic target. p97 is a large molecular machine that has been shown to hexamerize in vitro, with each monomer consisting of an N domain responsible for binding to effector proteins and two AAA repeats (D1 and D2). However, structural studies are inconclusive or in disagreement with one another on several important features such as the locations of the N domains, the relative orientations of the D1 and D2 rings, and the dimensions of the central pore. Here, we present atomic-scale simulations of the p97 hexamer in the prehydrolysis, transition, and post-hydrolysis states. To improve the agreement between low- and high-resolution experimental studies, we first use a biased simulation technique, molecular dynamics flexible fitting (MDFF), to improve the correlation between the structures described in these experiments. We follow this with extended, classical molecular dynamics simulations, which not only show that structures generated in the MDFF phase are stable, but reveal insights into the dynamics important to each state. Simulation results suggest a hybrid model for hydrolysis, in which the N and D2 domains are dynamic while the D1 domains are relatively static, salt bridges stabilize the position of the N domains in the pre-hydrolysis state, and the rings formed by D1 and D2 rotate relative to one another.