Computational active site analysis of molecular pathways to improve functional classification of enzymes

This study describes a method to computationally assess the function of homologous enzymes through small molecule binding interaction energy. Three experimentally determined X-ray structures and four enzyme models from ornithine cyclo-deaminase, alanine dehydrogenase, and mu-crystallin were used in combination with nine small molecules to derive a function score (FS) for each enzyme-model combination. While energy values varied for a single molecule-enzyme combination due to differences in the active sites, we observe that the binding energies for the entire pathway were proportional for each set of small molecules investigated. This proportionality of energies for a reaction pathway appears to be dependent on the amino acids in the active site and their direct interactions with the small molecules, which allows a function score (FS) to be calculated to assess the specificity of each enzyme. Potential of mean force (PMF) calculations were used to obtain the energies, and the resulting FS values demonstrate that a measurement of function may be obtained using differences between these PMF values. Additionally, limitations of this method are discussed based on: (a) larger substrates with significant conformational flexibility; (b) low homology enzymes; and (c) open active sites. This method should be useful in accurately predicting specificity for single enzymes that have multiple steps in their reactions and in high throughput computational methods to accurately annotate uncharacterized proteins based on active site interaction analysis.     

Using a neural network and spatial clustering to predict the location of active sites in enzymes

Structural genomics projects aim to provide a sharp increase in the number of structures of functionally unannotated, and largely unstudied, proteins. Algorithms and tools capable of deriving information about the nature, and location, of functional sites within a structure are increasingly useful therefore. Here, a neural network is trained to identify the catalytic residues found in enzymes, based on an analysis of the structure and sequence. The neural network output, and spatial clustering of the highly scoring residues are then used to predict the location of the active site.A comparison of the performance of differently trained neural networks is presented that shows how information from sequence and structure come together to improve the prediction accuracy of the network. Spatial clustering of the network results provides a reliable way of finding likely active sites. In over 69% of the test cases the active site is correctly predicted, and a further 25% are partially correctly predicted. The failures are generally due to the poor quality of the automatically generated sequence alignments.We also present predictions identifying the active site, and potential functional residues in five recently solved enzyme structures, not used in developing the method. The method correctly identifies the putative active site in each case. In most cases the likely functional residues are identified correctly, as well as some potentially novel functional groups.     

Combining structural genomics and enzymology: completing the picture in metabolic pathways and enzyme active sites

An important goal of structural genomics is to complete the structural analysis of all the enzymes in metabolic pathways and to understand the structural similarities and differences. A preliminary glimpse of this type of analysis was achieved before structural genomics efforts with the glycolytic pathway and efforts are underway for many other pathways, including that of catecholamine metabolism. Structural enzymology necessitates a complete structural characterization, even for highly homologous proteins (greater than 80% sequence homology), as every active site has distinct structural features and it is these active site differences that distinguish one enzyme from another. Short cuts with homology modeling cannot be taken with our current knowledge base. Each enzyme structure in a pathway needs to be determined, including structures containing bound substrates, cofactors, products and transition state analogs, in order to obtain a complete structural and functional understanding of pathway-related enzymes.     

Structural and mechanistic analysis of engineered trichodiene synthase enzymes from Trichoderma harzianum: towards higher catalytic activities empowering sustainable agriculture

Trichoderma spp. are well-known bioagents for the plant growth promotion and pathogen suppression. The beneficial activities of the fungus Trichoderma spp. are attributed to their ability to produce and secrete certain secondary metabolites such as trichodermin that belongs to trichothecene family of molecules. The initial steps of trichodermin biosynthetic pathway in Trichoderma are similar to the trichothecenes from Fusarium sporotrichioides. Trichodiene synthase (TS) encoded by tri5 gene in Trichoderma catalyses the conversion of farnesyl pyrophosphate to trichodiene as reported earlier. In this study, we have carried out a comprehensive comparative sequence and structural analysis of the TS, which revealed the conserved residues involved in catalytic activity of the protein. In silico, modelled tertiary structure of TS protein showed stable structural behaviour during simulations. Two single-substitution mutants, i.e. D109E, D248Y and one double-substitution mutant (D109E and D248Y) of TS with potentially higher activities are screened out. The mutant proteins showed more stability than the wild type, an increased number of electrostatic interactions and better binding energies with the ligand, which further elucidates the amino acid residues involved in the reaction mechanism. These results will lead to devise strategies for higher TS activity to ultimately enhance the trichodermin production by Trichoderma spp. for its better exploitation in the sustainable agricultural practices.     

A computational method for the analysis and prediction of protein:phosphopeptide-binding sites

Phosphopeptide-binding domains, including the FHA, SH2, WW, WD40, MH2, and Polo-box domains, as well as the 14-3-3 proteins, exert control functions in important processes such as cell growth, division, differentiation, and apoptosis. Structures and mechanisms of phosphopeptide binding are generally diverse, revealing few general principles. A computational method for analysis of phosphopeptide-binding domains was therefore developed to elucidate the physical and chemical nature of phosphopeptide binding, given this lack of structural similarity. The surfaces of nine phosphopeptide-binding proteins, representing seven distinct classes of phosphopeptide-binding modules, were discretized, and encoded with information about amino acid identity, surface curvature, and electrostatic potential at every point on the surface in order to identify local surface properties enriched in phosphoresidue contact sites. Cross-validation indicated that propensities corresponding to this enrichment calculated from a subset of the training data could be used to predict the phosphoresidue contact site on proteins not used in training with no false negative results, and with few unconfirmed positive predictions. The locations of phosphoresidue contact sites were then predicted on the surfaces of the checkpoint kinase Chk1 and the BRCA1 BRCT repeat domain, and these predictions are consistent with recent experimental evidence.     

Analysis of crystal structures of aspartic proteinases: On the role of amino acid residues adjacent to the catalytic site of pepsin-like enzymes

To elucidate the role of amino acid residues adjacent to the catalytic site of pepsin-like enzymes, we analyzed and compared the crystal structures of these enzymes, their complexes with inhibitors, and zymogens in the active site area (a total of 82 structures). In addition to the water molecule (W1) located between the active carboxyls and playing a role of the nucleophile during catalytic reaction, another water molecule (W2) at the vicinity of the active groups was found to be completely conserved. This water molecule plays an essential role in formation of a chain of hydrogen-bonded residues between the active site flap and the active carboxyls on ligand binding. These data suggest a new approach to understanding the role of residues around the catalytic site, which can assist the development of the catalytic reaction. The influence of groups adjacent to the active carboxyls is manifested by pepsin activity at pH 1.0. Some features of pepsin-like enzymes and their mutants are discussed in the framework of the approach.     

Systematic development of computational models for the catalytic site in galactose oxidase: impact of outer-sphere residues on the geometric and electronic structures

A systematic in silico approach has been employed to generate sound, experimentally validated active-site models for galactose oxidase (GO) using a hybrid density functional, B(38HF)P86. GO displays three distinct oxidation states: oxidized [Cu(II)-Y*]; semireduced [Cu(II)-Y]; and reduced [Cu(I)-Y]. Only the [Cu(II)-Y*] and the [Cu(I)-Y] states are assumed to be involved in the catalytic cycle, but their structures have not yet been determined. We have developed several models (1-7) for the [Cu(II)-Y*] state that were evaluated by comparison of our computational results with experimental data. An extended model system (6) that includes solvent molecules and second coordination sphere residues (R330, Y405, and W290) is essential to obtain an experimentally correct electronic structure of the active site. The optimized structure of 6 resulted in a five-coordinate Cu site with a protein radical centered on the Tyr-Cys cofactor. We further validated our converged model with the largest model (7) that included additional outer-sphere residues (Q406, H334, Y329, G513, and T580) and water molecules. Adding these residues did not affect significantly the active site's electronic and geometric structures. Using both 6 and 7, we explored the redox dependence of the active-site structure. We obtained four- and three-coordinate Cu sites for [Cu(II)-Y] and [Cu(I)-Y] states, respectively, that corroborate well with the experimental data. The relative energies of these states were validated by a comparison with experimental redox potentials. Collectively, our computational GO models well reproduce the physicochemical characteristics of the individual states, including their redox behaviors.     

Heterogeneous behavior of metalloproteins toward metal ion binding and selectivity: insights from molecular dynamics studies

About one-third of the existing proteins require metal ions as cofactors for their catalytic activities and structural complexities. While many of them bind only to a specific metal, others bind to multiple (different) metal ions. However, the exact mechanism of their metal preference has not been deduced to clarity. In this study, we used molecular dynamics (MD) simulations to investigate whether a cognate metal (bound to the structure) can be replaced with other similar metal ions. We have chosen seven different proteins (phospholipase A₂, sucrose phosphatase, pyrazinamidase, cysteine dioxygenase (CDO), plastocyanin, monoclonal anti-CD4 antibody Q425, and synaptotagmin 1 C₂B domain) bound to seven different divalent metal ions (Ca²⁺, Mg²⁺, Zn²⁺, Fe²⁺, Cu²⁺, Ba²⁺, and Sr²⁺, respectively). In total, 49 MD simulations each of 50 ns were performed and each trajectory was analyzed independently. Results demonstrate that in some cases, cognate metal ions can be exchanged with similar metal ions. On the contrary, some proteins show binding affinity specifically to their cognate metal ions. Surprisingly, two proteins CDO and plastocyanin which are known to bind Fe²⁺ and Cu²⁺, respectively, do not exhibit binding affinity to any metal ion. Furthermore, the study reveals that in some cases, the active site topology remains rigid even without cognate metals, whereas, some require them for their active site stability. Thus, it will be interesting to experimentally verify the accuracy of these observations obtained computationally. Moreover, the study can help in designing novel active sites for proteins to sequester metal ions particularly of toxic nature.     

A critical role for the histidine residues in the catalytic function of acyl-CoA:cholesterol acyltransferase catalysis: evidence for catalytic difference between ACAT1 and ACAT2

To investigate a role for histidine residues in the expression of normal acyl-CoA:cholesterol acyltransferase (ACAT) activity, the histidine residues located at five different positions in two isoenzymes were substituted by alanine, based on the sequence homology between ACAT1 and ACAT2. Among the 10 mutants generated by baculovirus expression technology, H386A-ACAT1, H460A-ACAT1, H360A-ACAT2, and H399A-ACAT2 lost their enzymatic activity completely. A reduction in catalytic activity is unlikely to result from structural changes in the substrate-binding pocket, because their substrate-binding affinities were normal. However, the enzymatic activity of H386A-ACAT1 was restored to <37% of the level of the wild-type activity when cholesterol was replaced by 25-hydroxycholesterol as substrate. H527A-ACAT1 and H501A-ACAT2, termed carboxyl end mutants, exhibit activities of ∼96% and ∼75% of that of the wild-type. Interestingly, H425A-ACAT1 showed 59% of the wild-type activity, in contrast to its equivalent mutant, H399A-ACAT2. These results demonstrate that the histidine residues located at the active site are very crucial both for the catalytic activity of the enzyme and for distinguishing ACAT1 from ACAT2 with respect to enzyme catalysis and substrate specificity.     

Interconverting the catalytic activities of (betaalpha)(8)-barrel enzymes from different metabolic pathways: sequence requirements and molecular analysis

The (betaalpha)(8)-barrel enzymes N'-[(5'-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide isomerase (tHisA) and imidazole glycerol phosphate synthase (tHisF) from Thermotoga maritima catalyze two successive reactions in the biosynthesis of histidine. In both enzymes, aspartate residues at the C-terminal end of beta-strand 1 (Asp8 in tHisA and Asp11 in tHisF) and beta-strand 5 (Asp127 in tHisA and Asp130 in tHisF) are essential for catalytic activity. It was demonstrated earlier that in tHisA the substitution of Asp127 by valine (tHisA-D127V) generates phosphoribosylanthranilate isomerase (TrpF) activity, a related (betaalpha)(8)-barrel enzyme participating in tryptophan biosynthesis. It is shown here that in tHisF the corresponding substitution of Asp130 by valine (tHisF-D130V) also generates TrpF activity. To determine the effectiveness of individual amino acid exchanges in these conversions, each of the 20 standard amino acid residues was introduced at position 127 of tHisA and 130 of tHisF by saturation random mutagenesis. The tHisA-D127X and tHisF-D130X variants with TrpF activity were identified by selection in vivo, and the proteins purified and characterized. The results obtained show that removal of the negatively charged carboxylate side-chain at the C-terminal end of beta-strand 5 is sufficient to establish TrpF activity in tHisA and tHisF, presumably because it allows the binding of the negatively charged TrpF substrate, phosphoribosylanthranilate. In contrast, the double mutants tHisA-D8N+D127V and tHisF-D11N+D130V did not show detectable activity, demonstrating that the aspartate residues at the C-terminal end of beta-strand 1 are essential for catalysis of the TrpF reaction. The ease with which TrpF activity can be established on both the tHisA and tHisF scaffolds supports the evolutionary relationship of these three enzymes and highlights the functional plasticity of the (betaalpha)(8)-barrel enzyme fold.     

Crystal structure of Stefin A in complex with cathepsin H: N-terminal residues of inhibitors can adapt to the active sites of endo- and exopeptidases

Binding of cystatin-type inhibitors to papain-like exopeptidases cannot be explained by the stefin B-papain complex. The crystal structure of human stefin A bound to an aminopeptidase, porcine cathepsin H, has been determined in monoclinic and orthorhombic crystal forms at 2.8A and 2.4A resolutions, respectively. The asymmetric unit of each form contains four complexes. The structures are similar to the stefin B-papain complex, but with a few distinct differences. On binding, the N-terminal residues of stefin A adopt the form of a hook, which pushes away cathepsin H mini-chain residues and distorts the structure of the short four residue insertion (Lys155A-Asp155D) unique to cathepsin H. Comparison with the structure of isolated cathepsin H shows that the rims of the cathepsin H structure are slightly displaced (up to 1A) from their position in the free enzyme. Furthermore, comparison with the stefin B-papain complex showed that molecules of stefin A bind about 0.8A deeper into the active site cleft of cathepsin H than stefin B into papain. The approach of stefin A to cathepsin H induces structural changes along the interaction surface of both molecules, whereas no such changes were observed in the stefin B-papain complex. Carboxymethylation of papain seems to have prevented the formation of the genuine binding geometry between a papain-like enzyme and a cystatin-type inhibitor as we observe it in the structure presented here.     

Structural and dynamical properties of manganese catalase and the synthetic protein DF1 and their implication for reactivity from classical molecular dynamics calculations

There is a pressing need for accurate force fields to assist in metalloprotein analysis and design. Recent work on the design of mimics of dimetal proteins highlights the requirements for activity. DF1 is a de novo designed protein, which mimics the overall fold and active site geometry of a series of diiron and dimanganese proteins. Specifically, the dimanganese form of DF1 is a mimic of the natural enzyme manganese catalase, which catalyzes the dismutation reaction of hydrogen peroxide into water and oxygen. During catalytic turnover, the active site has to accommodate both the reduced and the oxidized state of the dimanganese core. The biomimetic protein DF1 is only stable in the reduced form and thus not active. Furthermore, the synthetic protein features an additional bridging glutamate sidechain, which occupies the substrate binding site. The goal of this study is to develop classical force fields appropriate for design of such important dimanganese proteins. To this aim, we use a nonbonded model to represent the metal-ligand interactions, which implicitly takes into account charge transfer and local polarization effects between the metal and its ligands. To calibrate this approach, we compare and contrast geometric and dynamical properties of manganese catalase and DF1. Having demonstrated a good correspondence with experimental structural data, we examine the effect of mutating the bridging glutamate to aspartate (M1) and serine (M2). Classical MD based on the refined forcefield shows that these point mutations affect not only the immediate coordination sphere of the manganese ions, but also the relative position of the helices, improving the similarity to Mn-catalase, especially in case of M2. On the basis of these findings, classical molecular dynamics calculations with the active site parameterization scheme introduced herein seem to be a promising addition to the protein design toolbox.     

Molecular modeling,dynamics studies and density functional theory approaches to identify potential inhibitors of SIRT4 protein from Homo sapiens : a novel target for the treatment of type 2 diabetes

Type 2 diabetes is one of the biggest health challenges in the world and WHO projects it to be the 7th leading cause of death in 2030. It is a chronic condition affecting the way our body metabolizes sugar. Insulin resistance is high risk factor marked by expression of Lipoprotein Lipases and Peroxisome Proliferator-Activated Receptor that predisposes to type 2 diabetes. AMP-dependent protein kinase in AMPK signaling pathway is a central sensor of energy status. Deregulation of AMPK signaling leads to inflammation, oxidative stress, and deactivation of autophagy which are implicated in pathogenesis of insulin resistance. SIRT4 protein deactivates AMPK as well as directly inhibits insulin secretion. SIRT4 overexpression leads to dyslipidimeia, decreased fatty acid oxidation, and lipogenesis which are the characteristic features of insulin resistance promoting type 2 diabetes. This makes SIRT4 a novel therapeutic target to control type 2 diabetes. Virtual screening and molecular docking studies were performed to obtain potential ligands. To further optimize the geometry of protein–ligand complexes Quantum Polarized Ligand Docking was performed. Binding Free Energy was calculated for the top three ligand molecules. In view of exploring the stereoelectronic features of the ligand, density functional theory approach was implemented at B3LYP/6-31G* level. 30 ns MD simulation studies of the protein–ligand complexes were done. The present research work proposes ZINC12421989 as potential inhibitor of SIRT4 with docking score (?7.54 kcal/mol), docking energy (?51.34 kcal/mol), binding free energy (?70.21 kcal/mol), and comparatively low energy gap (?0.1786 eV) for HOMO and LUMO indicating reactivity of the lead molecule.     

Conformation and microenvironment of the active site of a low molecular weight 1,4-beta-D-glucan glucanohydrolase from an alkalothermophilic Thermomonospora sp.: involvement of lysine and cysteine residues

Conformation and microenvironment at the active site of 1,4-beta-D-glucan glucanohydrolase was probed with fluorescent chemo-affinity labeling using o-phthalaldehyde. OPTA has been known to form a fluorescent isoindole derivative by cross-linking the proximal thiol and amino groups of cysteine and lysine. Modification of lysine of the enzyme by TNBS and of cysteine residue by PHMB abolished the ability of the enzyme to form an isoindole derivative with OPTA. Kinetic analysis of the TNBS and PHMB-modified enzyme suggested the presence of essential lysine and cysteine residues, respectively, at the active site of the enzyme. The substrate protection of the enzyme with carboxymethylcellulose (CMC) confirmed the involvement of lysine and cysteine residues in the active site of the enzyme. Multiple sequence alignment of peptides obtained by tryptic digestion of the enzyme showed cysteine is one of the conserved amino acids corroborating the chemical modification studies.     

A method for activity staining after native polyacrylamide gel electrophoresis using a coupled enzyme assay and fluorescence detection: application to the analysis of several glycolytic enzymes.

We describe a method for the detection of isoforms of several glycolytic enzymes by activity staining after native PAGE. The staining is based on coupled enzyme assays carried out on the gel after electrophoresis and is linked to the disappearance of NADH, which is visualized by fluorescence. This method offers reliable and sensitive detection for phosphoenolpyruvate carboxylase, PPi-dependent phosphofructokinase, and pyruvate kinase from plant tissues. It can be applied to the detection of all enzymes which are normally detected spectrophotometrically using coupled enzyme assays consuming NAD(P)H.