Computational studies on the catalytic potential of the double active site for enzyme engineering
Ontology highlight
ABSTRACT:
SUBMITTER: Krishna N
PROVIDER: S-EPMC11297320 | biostudies-literature | 2024 Jan
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A loop closure-based sequential algorithm, PRODA_MATCH, was developed to match catalytic residues onto a scaffold for enzyme design in silico. The computational complexity of this algorithm is polynomial with respect to the number of active sites, the number of catalytic residues, and the maximal iteration number of cyclic coordinate descent steps. This matching algorithm is independent of a rotamer library that enables the catalytic residue to take any required conformation during the reaction coordinate. The catalytic geometric parameters defined between functional groups of transition state (TS) and the catalytic residues are continuously optimized to identify the accurate position of the TS. Pseudo-spheres are introduced for surrounding residues, which make the algorithm take binding into account as early as during the matching process. Recapitulation of native catalytic residue sites was used as a benchmark to evaluate the novel algorithm. The calculation results for the test set show that the native catalytic residue sites were successfully identified and ranked within the top 10 designs for 7 of the 10 chemical reactions. This indicates that the matching algorithm has the potential to be used for designing industrial enzymes for desired reactions.
| S-EPMC3190151 | biostudies-literature
Project description:Recent studies reveal that the core sequences of many proteins were nearly optimized for stability by natural evolution. Surface residues, by contrast, are not so optimized, presumably because protein function is mediated through surface interactions with other molecules. Here, we sought to determine the extent to which the sequences of protein ligand-binding and enzyme active sites could be predicted by optimization of scoring functions based on protein ligand-binding affinity rather than structural stability. Optimization of binding affinity under constraints on the folding free energy correctly predicted 83% of amino acid residues (94% similar) in the binding sites of two model receptor-ligand complexes, streptavidin-biotin and glucose-binding protein. To explore the applicability of this methodology to enzymes, we applied an identical algorithm to the active sites of diverse enzymes from the peptidase, beta-gal, and nucleotide synthase families. Although simple optimization of binding affinity reproduced the sequences of some enzyme active sites with high precision, imposition of additional, geometric constraints on side-chain conformations based on the catalytic mechanism was required in other cases. With these modifications, our sequence optimization algorithm correctly predicted 78% of residues from all of the enzymes, with 83% similar to native (90% correct, with 95% similar, excluding residues with high variability in multiple sequence alignments). Furthermore, the conformations of the selected side chains were often correctly predicted within crystallographic error. These findings suggest that simple selection pressures may have played a predominant role in determining the sequences of ligand-binding and active sites in proteins.
| S-EPMC1177389 | biostudies-literature
Project description:The ubiquitin proteasome system (UPS) signals for degradation of proteins through attachment of K48-linked polyubiquitin chains, or alterations in protein-protein recognition through attachment of K63-linked chains. Target proteins are ubiquitinated in three sequential chemical steps by a three-component enzyme system. Ubiquitination, or E2 enzymes, catalyze the central step by facilitating reaction of a target protein lysine with the C-terminus of Ub that is attached to the active site cysteine of the E2 through a thioester bond. E2 reactivity is modulated by dynamics of an active site gate, whose central residue packs against the active site cysteine in a closed conformation. Interestingly, for the E2 Ubc13, which specifically catalyzes K63-linked ubiquitination, the central gate residue adopts an open conformation. We set out to determine if active site gate dynamics play a role in catalysis for E2-25K, which adopts the canonical, closed gate conformation, and which selectively synthesizes K48-linked ubiquitin chains. Gate dynamics were characterized using mutagenesis of key residues, combined with enzyme kinetics measurements, and main chain NMR relaxation. The experimental data were interpreted with all atom MD simulations. The data indicate that active site gate opening and closing rates for E2-25K are precisely balanced.
| S-EPMC5934386 | biostudies-literature
Project description:Berberine bridge enzyme (BBE) is a paradigm for the class of bicovalently flavinylated oxidases, which catalyzes the oxidative cyclization of (S)-reticuline to (S)-scoulerine. His174 was identified as an important active site residue because of its role in the stabilization of the reduced state of the flavin cofactor. It is also strictly conserved in the family of BBE-like oxidases. Here, we present a detailed biochemical and structural characterization of a His174Ala variant supporting its importance during catalysis and for the structural organization of the active site. Substantial changes in all kinetic parameters and a decrease in midpoint potential were observed for the BBE His174Ala variant protein. Moreover, the crystal structure of the BBE His174Ala variant showed significant structural rearrangements compared to wild-type enzyme. On the basis of our findings, we propose that His174 is part of a hydrogen bonding network that stabilizes the negative charge at the N1-C2=O locus via interaction with the hydroxyl group at C2' of the ribityl side chain of the flavin cofactor. Hence, replacement of this residue with alanine reduces the stabilizing effect for the transiently formed negative charge and results in drastically decreased kinetic parameters as well as a lower midpoint redox potential.
| S-EPMC3413249 | biostudies-literature
Project description:Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.
| S-EPMC10095239 | biostudies-literature
Project description:Rational enzyme design presents a major challenge that has not been overcome by computational approaches. One of the key challenges is the difficulty in assessing the magnitude of the maximum possible catalytic activity. In an attempt to overcome this challenge, we introduce a strategy that takes an active enzyme (assuming that its activity is close to the maximum possible activity), design mutations that reduce the catalytic activity, and then try to restore that catalysis by mutating other residues. Here we take as a test case the enzyme haloalkane dehalogenase (DhlA), with a 1,2-dichloroethane substrate. We start by demonstrating our ability to reproduce the results of single mutations. Next, we design mutations that reduce the enzyme activity and finally design double mutations that are aimed at restoring the activity. Using the computational predictions as a guide, we conduct an experimental study that confirms our prediction in one case and leads to inconclusive results in another case with 1,2-dichloroethane as substrate. Interestingly, one of our predicted double mutants catalyzes dehalogenation of 1,2-dibromoethane more efficiently than the wild-type enzyme.
| S-EPMC6329970 | biostudies-literature
Project description:M-CSA (Mechanism and Catalytic Site Atlas) is a database of enzyme active sites and reaction mechanisms that can be accessed at www.ebi.ac.uk/thornton-srv/m-csa. Our objectives with M-CSA are to provide an open data resource for the community to browse known enzyme reaction mechanisms and catalytic sites, and to use the dataset to understand enzyme function and evolution. M-CSA results from the merging of two existing databases, MACiE (Mechanism, Annotation and Classification in Enzymes), a database of enzyme mechanisms, and CSA (Catalytic Site Atlas), a database of catalytic sites of enzymes. We are releasing M-CSA as a new website and underlying database architecture. At the moment, M-CSA contains 961 entries, 423 of these with detailed mechanism information, and 538 with information on the catalytic site residues only. In total, these cover 81% (195/241) of third level EC numbers with a PDB structure, and 30% (840/2793) of fourth level EC numbers with a PDB structure, out of 6028 in total. By searching for close homologues, we are able to extend M-CSA coverage of PDB and UniProtKB to 51 993 structures and to over five million sequences, respectively, of which about 40% and 30% have a conserved active site.
| S-EPMC5753290 | biostudies-literature
Project description:How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme's reaction cycle. To examine positioning and motions and test catalytic proposals, we collected "room temperature" X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Ensemble analyses indicated limited change through KSI's reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI's general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure-function relationships, providing the basis for a new era of "ensemble-function" interrogation of enzymes.
| S-EPMC7776779 | biostudies-literature
Project description:Enzymological paradigms have shifted recently to acknowledge the biological importance of catalytic promiscuity. However, catalytic promiscuity is a poorly understood property, and no thermodynamic treatment has described the conformational landscape of promiscuous versus substrate-specific enzymes. Here, two structurally similar glutathione transferase (GST, glutathione S-transferase) isoforms with high specificity or high promiscuity are compared. Differential scanning calorimetry (DSC) indicates a reversible low temperature transition for the promiscuous GSTA1-1 that is not observed with substrate-specific GSTA4-4. This transition is assigned to rearrangement of the C terminus at the active site of GSTA1-1 based on the effects of ligands and mutations. Near-UV and far-UV circular dichroism indicate that this transition is due to repacking of tertiary contacts with the remainder of the subunit, rather than "unfolding" of the C terminus per se. Analysis of the DSC data using a modified Landau theory indicates that the local conformational landscape of the active site of GSTA1-1 is smooth, with barrierless transitions between states. The partition function of the C-terminal states is a broad unimodal distribution at all temperatures within this DSC transition. In contrast, the remainder of the GSTA1-1 subunit and the GSTA4-4 protein exhibit folded and unfolded macrostates with a significant energy barrier separating them. Their partition function includes a sharp unimodal distribution of states only at temperatures that yield either folded or unfolded macrostates. At intermediate temperatures the partition function includes a bimodal distribution. The barrierless rearrangement of the GSTA1-1 active site within a local smooth energy landscape suggests a thermodynamic basis for catalytic promiscuity.
| S-EPMC3234982 | biostudies-literature
Project description:The positioning of catalytic groups within proteins plays an important role in enzyme catalysis, and here we investigate the positioning of the general base in the enzyme ketosteroid isomerase (KSI). The oxygen atoms of Asp38, the general base in KSI, were previously shown to be involved in anion-aromatic interactions with two neighboring Phe residues. Here we ask whether those interactions are sufficient, within the overall protein architecture, to position Asp38 for catalysis or whether the side chains that pack against Asp38 and/or the residues of the structured loop that is capped by Asp38 are necessary to achieve optimal positioning for catalysis. To test positioning, we mutated each of the aforementioned residues, alone and in combinations, in a background with the native Asp general base and in a D38E mutant background, as Glu at position 38 was previously shown to be mispositioned for general base catalysis. These double-mutant cycles reveal positioning effects as large as 10(3)-fold, indicating that structural features in addition to the overall protein architecture and the Phe residues neighboring the carboxylate oxygen atoms play roles in positioning. X-ray crystallography and molecular dynamics simulations suggest that the functional effects arise from both restricting dynamic fluctuations and disfavoring potential mispositioned states. Whereas it may have been anticipated that multiple interactions would be necessary for optimal general base positioning, the energetic contributions from positioning and the nonadditive nature of these interactions are not revealed by structural inspection and require functional dissection. Recognizing the extent, type, and energetic interconnectivity of interactions that contribute to positioning catalytic groups has implications for enzyme evolution and may help reveal the nature and extent of interactions required to design enzymes that rival those found in biology.
| S-EPMC4004248 | biostudies-literature