skip to main content

Title: Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
The Journal of Chemical Physics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Glycoside hydrolase enzymes are important for hydrolyzing the β-1,4 glycosidic bond in polysaccharides for deconstruction of carbohydrates. The two-step retaining reaction mechanism of Glycoside Hydrolase Family 7 (GH7) was explored with different sized QM-cluster models built by the Residue Interaction Network ResidUe Selector (RINRUS) software using both the wild-type protein and its E217Q mutant. The first step is the glycosylation, in which the acidic residue 217 donates a proton to the glycosidic oxygen leading to bond cleavage. In the subsequent deglycosylation step, one water molecule migrates into the active site and attacks the anomeric carbon. Residue interaction-based QM-cluster models lead to reliable structural and energetic results for proposed glycoside hydrolase mechanisms. The free energies of activation for glycosylation in the largest QM-cluster models were predicted to be 19.5 and 31.4 kcal mol −1 for the wild-type protein and its E217Q mutant, which agree with experimental trends that mutation of the acidic residue Glu217 to Gln will slow down the reaction; and are higher in free energy than the deglycosylation transition states (13.8 and 25.5 kcal mol −1 for the wild-type protein and its mutant, respectively). For the mutated protein, glycosylation led to a low-energy product. This thermodynamic sink may correspond to the intermediate state which was isolated in the X-ray crystal structure. Hence, the glycosylation is validated to be the rate-limiting step in both the wild-type and mutated enzyme. 
    more » « less
  2. Saccharomyces cerevisiae OYE 3 shares 80% sequence identity with the well-studied Saccharomyces pastorianus OYE 1; however, wild-type OYE 3 shows different stereoselectivities toward some alkene substrates. Site-saturation mutagenesis of Trp 116 in OYE 3 followed by substrate profiling showed that the mutations had relatively little effect, opposite to that observed previously for OYE 1. The X-ray crystal structures of unliganded and phenol-bound OYE 3 were solved to 1.8 and 1.9 Å resolution, respectively. Both structures were nearly identical to that of OYE 1, with only a single amino acid difference in the active site region (Ser 296 versus Phe 296, part of loop 6). Despite their essentially identical static X-ray structures, molecular dynamics (MD) simulations revealed that loop 6 conformations differed significantly in solution between OYE 3 and OYE 1. In OYE 3, loop 6 remained nearly as open as observed in the crystal structure; by contrast, loop 6 closed over the active site of OYE 1 by ca. 4 Å. Loop closure likely generates a greater number of active site protein contacts for substrate bound to OYE 1 as compared to OYE 3. These differences provide an explanation for the differing stereoselectivities of OYE 3 and OYE 1, despite their nearly identical X-ray crystal structures. 
    more » « less
  3. Although the 3D structures of active and inactive cannabinoid receptors type 2 (CB2) are available, neither the X-ray crystal nor the cryo-EM structure of CB2-orthosteric ligand-modulator has been resolved, prohibiting the drug discovery and development of CB2 allosteric modulators (AMs). In the present work, we mainly focused on investigating the potential allosteric binding site(s) of CB2. We applied different algorithms or tools to predict the potential allosteric binding sites of CB2 with the existing agonists. Seven potential allosteric sites can be observed for either CB2-CP55940 or CB2-WIN 55,212-2 complex, among which sites B, C, G and K are supported by the reported 3D structures of Class A GPCRs coupled with AMs. Applying our novel algorithm toolset-MCCS, we docked three known AMs of CB2 including Ec2la (C-2), trans-β-caryophyllene (TBC) and cannabidiol (CBD) to each site for further comparisons and quantified the potential binding residues in each allosteric binding site. Sequentially, we selected the most promising binding pose of C-2 in five allosteric sites to conduct the molecular dynamics (MD) simulations. Based on the results of docking studies and MD simulations, we suggest that site H is the most promising allosteric binding site. We plan to conduct bio-assay validations in the future. 
    more » « less
  4. Ethylenediaminetetraacetic acid (EDTA), which has two amine and four carboxylate protonation sites, forms stable complexes with lanthanide ions. This work analyzes the coordination structure, in atomic resolution, of the Eu3+ ion complexed with EDTA in all its protonation states in aqueous solution. Eu-EDTA complexes were modeled using classical molecular dynamics (MD) simulations using force field parameters optimized with ab initio molecular dynamics (AIMD) simulations. Structures from the MD simulations were used to predict extended X-ray absorption fine structure (EXAFS) spectra and compared with EXAFS measurements of the Eu3+ aqua ion and Eu-EDTA complexes at pH 3 and 11. This work details how Eu-EDTA complex coordination structures change with increasing protonation of the EDTA ligand in the complex, from the tightly bound unprotonated complex to the unbinding of the fully protonated EDTA ligand from the Eu3+ ion as both become solvated by water. Agreement between predicted and measured EXAFS spectra supports the findings from simulation. 
    more » « less
  5. Room-temperature macromolecular crystallography allows protein structures to be determined under close-to-physiological conditions, permits dynamic freedom in protein motions and enables time-resolved studies. In the case of metalloenzymes that are highly sensitive to radiation damage, such room-temperature experiments can present challenges, including increased rates of X-ray reduction of metal centres and site-specific radiation-damage artefacts, as well as in devising appropriate sample-delivery and data-collection methods. It can also be problematic to compare structures measured using different crystal sizes and light sources. In this study, structures of a multifunctional globin, dehaloperoxidase B (DHP-B), obtained using several methods of room-temperature crystallographic structure determination are described and compared. Here, data were measured from large single crystals and multiple microcrystals using neutrons, X-ray free-electron laser pulses, monochromatic synchrotron radiation and polychromatic (Laue) radiation light sources. These approaches span a range of 18 orders of magnitude in measurement time per diffraction pattern and four orders of magnitude in crystal volume. The first room-temperature neutron structures of DHP-B are also presented, allowing the explicit identification of the hydrogen positions. The neutron data proved to be complementary to the serial femtosecond crystallography data, with both methods providing structures free of the effects of X-ray radiation damage when compared with standard cryo-crystallography. Comparison of these room-temperature methods demonstrated the large differences in sample requirements, data-collection time and the potential for radiation damage between them. With regard to the structure and function of DHP-B, despite the results being partly limited by differences in the underlying structures, new information was gained on the protonation states of active-site residues which may guide future studies of DHP-B. 
    more » « less