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
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):
- 1846408
- NSF-PAR ID:
- 10440375
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 6
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
null (Ed.)Although unsaturated organotrifluoroborates are common synthons in metal–organic chemistry, their transition metal complexes have received little attention. [CH 2 (3,5-(CH 3 ) 2 Pz) 2 ]Cu(CH 2 CHBF 3 ), (SIPr)Cu(MeCN)(CH 2 CHBF 3 ) and [CH 2 (3,5-(CH 3 ) 2 Pz) 2 ]Ag(CH 2 CHBF 3 ) represent rare, isolable molecules featuring a vinyltrifluoroborate ligand on coinage metals. The X-ray crystal structures show the presence of three-coordinate metal sites in these complexes. The vinyltrifluoroborate group binds asymmetrically to the metal site in [CH 2 (3,5-(CH 3 ) 2 Pz) 2 ]M(CH 2 CHBF 3 ) (M = Cu, Ag) with relatively closer M–C(H) 2 distances. The computed structures of [CH 2 (3,5-(CH 3 ) 2 Pz) 2 ]M(CH 2 CHBF 3 ) and M(CH 2 CHBF 3 ), however, have shorter M–C(H)BF 3 distances than M–C(H) 2 . These molecules feature various inter- or intra-molecular contacts involving fluorine of the BF 3 group, possibly affecting these M–C distances. The binding energies of [CH 2 CHBF 3 ] − to Cu + , Ag + and Au + have been calculated at the wB97XD/def2-TZVP level of theory, in the presence and absence of the supporting ligand CH 2 (3,5-(CH 3 ) 2 Pz) 2 . The calculation shows that Au + has the strongest binding to the [CH 2 CHBF 3 ] − ligand, followed by Cu + and Ag + , irrespective of the presence of the supporting ligand. However, in all three metals, the supporting ligand weakens the binding of olefin to the metal. The same trends were also found from the analysis of the σ-donation and π-backbonding interactions between the metal fragment and the π and π* orbitals of [CH 2 CHBF 3 ] − .more » « less
-
Cleavage of aromatic carbon–chlorine bonds is critical for the degradation of toxic industrial compounds. Here, we solved the X-ray crystal structure of chlorothalonil dehalogenase (Chd) from Pseudomonas sp. CTN-3, with 15 of its N-terminal residues truncated (Chd T ), using single-wavelength anomalous dispersion refined to 1.96 Å resolution. Chd has low sequence identity (<15%) compared with all other proteins whose structures are currently available, and to the best of our knowledge, we present the first structure of a Zn(II)-dependent aromatic dehalogenase that does not require a coenzyme. Chd T forms a “head-to-tail” homodimer, formed between two α-helices from each monomer, with three Zn(II)-binding sites, two of which occupy the active sites, whereas the third anchors a structural site at the homodimer interface. The catalytic Zn(II) ions are solvent-accessible via a large hydrophobic (8.5 × 17.8 Å) opening to bulk solvent and two hydrophilic branched channels. Each active-site Zn(II) ion resides in a distorted trigonal bipyramid geometry with His 117 , His 257 , Asp 116 , Asn 216 , and a water/hydroxide as ligands. A conserved His residue, His 114 , is hydrogen-bonded to the Zn(II)-bound water/hydroxide and likely functions as the general acid-base. We examined substrate binding by docking chlorothalonil (2,4,5,6-tetrachloroisophtalonitrile, TPN) into the hydrophobic channel and observed that the most energetically favorable pose includes a TPN orientation that coordinates to the active-site Zn(II) ions via a CN and that maximizes a π–π interaction with Trp 227 . On the basis of these results, along with previously reported kinetics data, we propose a refined catalytic mechanism for Chd-mediated TPN dehalogenation.more » « less
-
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
