EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization.

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.

Oxygen vacancies as active sites for water dissociation on rutile TiO2(110)

Through an interplay between scanning tunneling microscopy experiments and density functional theory calculations, we determine unambiguously the active surface site responsible for the dissociation of water molecules adsorbed on rutile TiO2(110). Oxygen vacancies in the surface layer are shown to dissociate H2O through the transfer of one proton to a nearby oxygen atom, forming two hydroxyl groups for every vacancy. The amount of water dissociation is limited by the density of oxygen vacancies present on the clean surface exclusively. The dissociation process sets in as soon as molecular water is able to diffuse to the active site.

Electrochemical Properties of Oxo-Manganese Complex Biomimicking Enzyme Active Sites and Its Electrocatalytic Application for Dopamine Determination

This work describes the characterization of the [Mn2 IV,IVO2(terpy)2(H2O)2]4+ complex in aqueous solution by UV-vis spectrophotometry, cyclic voltammetry, and linear sweep voltammetry with a rotating disk electrode. The pH effect, potential scan rate, effect of perfluorosulfonate polymer, and anion of supporting electrode on the electrochemical behavior of the modified electrode for better performance were investigated. The potential peak of the modified electrode was linearly dependent upon the ratio [ionic charge]/[ionic radius]. The modified electrode exerted an electrocatalytic effect on dopamine oxidation in aqueous solution with a decrease in the overpotential compared with the unmodified glassy carbon electrode. This way, the modified electrode showed an enzymatic biomimicking behavior. Tafel plot analyses were used to elucidate the kinetics and mechanism of dopamine oxidation. © 2013 Springer Science+Business Media New York.

Theoretical studies of mononuclear non-heme iron active sites

The quantum chemical investigations presented in this thesis use hybrid density functional theory to shed light on the catalytic mechanisms of mononuclear non-heme iron oxygenases, accommodating a ferrous ion in their active sites. More specifically, the dioxygen activation process and the subsequent oxidative reactions in the following enzymes were studied: tetrahydrobiopterin-dependent hydroxylases, naphthalene 1,2-dioxygenase and α-ketoglutarate-dependent enzymes. In light of many experimental efforts devoted to the functional mimics of non-heme iron oxygenases, the reactivity of functional analogues was also examined. The computed energetics and the available experimental data served to assess the feasibility of the reaction mechanisms investigated. Dioxygen activation in tetrahydrobiopterin- and α-ketoglutarate-dependent enzymes were found to involve a high-valent iron-oxo species, which was then capable of substrate hydroxylation. In the case of naphthalene 1,2-dioxygenase, the reactivity of an iron(III)-hydroxperoxo species toward the substrate was investigated and compared to the biomimetic counterpart.

Identification of active sites in amidase: Evolutionary relationship between amide bond- and peptide bond-cleaving enzymes

Mainly based on various inhibitor studies previously performed, amidases came to be regarded as sulfhydryl enzymes. Not completely satisfied with this generally accepted interpretation, we performed a series of site-directed mutagenesis studies on one particular amidase of Rhodococcus rhodochrous J1 that was involved in its nitrile metabolism. For these experiments, the recombinant amidase was produced as the inclusion body in Escherichia coli to greatly facilitate its recovery and subsequent purification. With regard to the presumptive active site residue Cys203, a Cys203 → Ala mutant enzyme still retained 11.5% of the original specific activity. In sharp contrast, substitutions in certain other positions in the neighborhood of Cys203 had a far more dramatic effect on the amidase. Glutamic acid substitution of Asp191 reduced the specific activity of the mutant enzyme to 1.33% of the wild-type activity. Furthermore, Asp191 → Asn substitution as well as Ser195 → Ala substitution completely abolished the specific activity. It would thus appear that, among various conserved residues residing within the so-called signature sequence common to all amidases, the real active site residues are Asp191 and Ser195 rather than Cys203. Inasmuch as an amide bond (CO-NH2) in the amide substrate is not too far structurally removed from a peptide bond (CO-NH-), the signature sequences of various amidases were compared with the active site sequences of various types of proteases. It was found that aspartic acid and serine residues corresponding to Asp191 and Ser195 of the Rhodococcus amidase are present within the active site sequences of aspartic proteinases, thus suggesting the evolutionary relationship between the two.

Phosphohistidyl active sites in polyphosphate kinase of Escherichia coli

In the synthesis of inorganic polyphosphate (polyP) from ATP by polyphosphate kinase (PPK; EC 2.7.4.1) of Escherichia coli, an N—P-linked phosphoenzyme was previously identified as the intermediate. The phosphate is presumed to be linked to N3 of the histidine residue because of its chemical stabilities and its resemblance to other enzymes known to contain N3-phosphohistidine. Tryptic digests of [32P]PPK contain a predominant 32P-labeled peptide that includes His-441. Of the 16 histidine residues in PPK of E. coli, 4 are conserved among several bacterial species. Mutagenesis of these 4 histidines shows that two (His-430 and His-598) are unaffected in function when mutated to glutamine, whereas two others (His-441 and His-460) mutated to glutamine or alanine fail to be phosphorylated, show no enzymatic activities, and fail to support polyP accumulation in cells bearing these mutant enzymes.

In situ characterization of porous VPO catalysts with fibrous structure: Identifying the redox behavior and the stability of active sites

Two VPO materials with fibrillar morphology have been prepared by the aid of electrospinning technique. One is a VPO carbon-supported material (VCF200) with fibrous morphology and very high surface area that is stable under oxidizing conditions up to 350C. The other material is a bulk mixed VPO oxide (VPO500) with fibrous structure obtained after optimizing the calcination of the carbon support in VCF200. Despite it is a bulk oxide material, this material exhibits a high surface area (> 60 m2/g). The redox behavior of both samples was monitored by in situ Raman spectroscopy under oxidation/reduction cycles. For the dehydrated supported sample (VCF200), the pyrophosphate phase (VO)2P2O7 (Raman ~930 cm-1) is detected, which has been described as the active phase (see Figure (a) below). This phase is quite stable since it does not disappear upon subsequent oxidation/reduction cycles. Under reduction conditions at 125C, in consecutive cycles, additional Raman bands appear at ~1090 cm-1 that are characteristic of the αII-VOPO4 phase. On the other hand, the bulk phases show a reversible behavior under redox cycles (Figure (b)). Under reducing conditions, a Raman band appears at ~980 cm-1 (β-VPO phase), whereas under oxidation conditions some segregation to VOx oxides occurs. Nevertheless, this segregation is reversible and the β-VPO phase forms again under reducing conditions. Thus, these results demonstrate that the active VPO phases of these fibrous catalysts are quite stable, and that their structure is reversible under several redox cycles, which make them suitable as oxidation catalysts.

Active site-directed protein regulation

Regulation of protein function is vital for the control of cellular processes. Proteins are often regulated by allosteric mechanisms, in which effecters bind to regulatory sites distinct from the active sites and alter protein function. Intrasteric regulation, directed at the active site and thus the counterpart of allosteric control, is now emerging as an important regulatory mechanism.

Functional analysis, using in vitro mutagenesis, of amino acids located in the phenylalanine hydroxylase active site

The 3-dimensionaI structure determination of rat phenylalanine hydroxylase (PAH) has identified potentially important amino acids lining the active site cleft with the majority of these having hydrophobic side-chains including several with aromatic side chains. Here we have analyzed the effect on rat PAH enzyme kinetics of in vitro mutagenesis of a number of these amino acids lining the PAH active site. Mutation of F299, Y324, F331, and Y343 caused a significant decrease in enzyme activity but no change in the K-m for substrate or cofactor. me conclude that these aromatic residues are essential for activity but are not significantly involved in binding of the substrate or cofactor. in contrast the PAH mutant, S349T, showed an 18-fold increase in K-m for phenylalanine, showing the first functional evidence that this residue was binding at or near the phenylalanine binding site. This confirms the recently published model for the binding of phenylalanine to the PAH active site that postulated S349 interacts with the amino group on the main chain of the phenylalanine molecule. This result differs with that found for the equivalent mutation (S395T), in the closely related tyrosine hydroxylase, which had no effect on substrate K-m, showing that while the architecture of the two active sites are very similar the amino acids that bind to the respective substrates are different. (C) 2000 Academic Press.

Identification of the yeast 20S proteasome catalytic centers and subunit interactions required for active-site formation

The proteasome is responsible for degradation of substrates of the ubiquitin pathway. 20S proteasomes are cylindrical particles with subunits arranged in a stack of four heptameric rings. The outer rings are composed of α subunits, and the inner rings are composed of β subunits. A well-characterized archaeal proteasome has a single type of each subunit, and the N-terminal threonine of the β subunit is the active-site nucleophile. Yeast proteasomes have seven different β subunits and exhibit several distinct peptidase activities, which were proposed to derive from disparate active sites. We show that mutating the N-terminal threonine in the yeast Pup1 β subunit eliminates cleavage after basic residues in peptide substrates, and mutating the corresponding threonine of Pre3 prevents cleavage after acidic residues. Surprisingly, neither mutation has a strong effect on cell growth, and they have at most minor effects on ubiquitin-dependent proteolysis. We show that Pup1 interacts with Pup3 in each β subunit ring. Our data reveal that different proteasome active sites contribute very differently to protein breakdown in vivo, that contacts between particular subunits in each β subunit ring are critical for active-site formation, and that active sites in archaea and different eukaryotes are highly similar.

Functions of FKBP12 and Mitochondrial Cyclophilin Active Site Residues In Vitro and In Vivo in Saccharomyces cerevisiae

Cyclophilin and FK506 binding protein (FKBP) accelerate cis–trans peptidyl-prolyl isomerization and bind to and mediate the effects of the immunosuppressants cyclosporin A and FK506. The normal cellular functions of these proteins, however, are unknown. We altered the active sites of FKBP12 and mitochondrial cyclophilin from the yeast Saccharomyces cerevisiae by introducing mutations previously reported to inactivate these enzymes. Surprisingly, most of these mutant enzymes were biologically active in vivo. In accord with previous reports, all of the mutant enzymes had little or no detectable prolyl isomerase activity in the standard peptide substrate-chymotrypsin coupled in vitro assay. However, in a variation of this assay in which the protease is omitted, the mutant enzymes exhibited substantial levels of prolyl isomerase activity (5–20% of wild-type), revealing that these mutations confer sensitivity to protease digestion and that the classic in vitro assay for prolyl isomerase activity may be misleading. In addition, the mutant enzymes exhibited near wild-type activity with two protein substrates, dihydrofolate reductase and ribonuclease T1, whose folding is accelerated by prolyl isomerases. Thus, a number of cyclophilin and FKBP12 “active-site” mutants previously identified are largely active but protease sensitive, in accord with our findings that these mutants display wild-type functions in vivo. One mitochondrial cyclophilin mutant (R73A), and also the wild-type human FKBP12 enzyme, catalyze protein folding in vitro but lack biological activity in vivo in yeast. Our findings provide evidence that both prolyl isomerase activity and other structural features are linked to FKBP and cyclophilin in vivo functions and suggest caution in the use of these active-site mutations to study FKBP and cyclophilin functions.

DNA polymerase active site is highly mutable: Evolutionary consequences

DNA polymerases contain active sites that are structurally superimposable and highly conserved in sequence. To assess the significance of this preservation and to determine the mutational burden that active sites can tolerate, we randomly mutated a stretch of 13 amino acids within the polymerase catalytic site (motif A) of Thermus aquaticus DNA polymerase I. After selection, by using genetic complementation, we obtained a library of approximately 8,000 active mutant DNA polymerases, of which 350 were sequenced and analyzed. This is the largest collection of physiologically active polymerase mutants. We find that all residues of motif A, except one (Asp-610), are mutable while preserving wild-type activity. A wide variety of amino acid substitutions were obtained at sites that are evolutionarily maintained, and conservative substitutions predominate at regions that stabilize tertiary structures. Several mutants exhibit unique properties, including DNA polymerase activity higher than the wild-type enzyme or the ability to incorporate ribonucleotide analogs. Bacteria dependent on these mutated polymerases for survival are fit to replicate repetitively. The high mutability of the polymerase active site in vivo and the ability to evolve altered enzymes may be required for survival in environments that demand increased mutagenesis. The inherent substitutability of the polymerase active site must be addressed relative to the constancy of nucleotide sequence found in nature.

Swiveling-domain mechanism for enzymatic phosphotransfer between remote reaction sites.

The crystal structure of pyruvate phosphate dikinase, a histidyl multiphosphotransfer enzyme that synthesizes adenosine triphosphate, reveals a three-domain molecule in which the phosphohistidine domain is flanked by the nucleotide and the phosphoenolpyruvate/pyruvate domains, with the two substrate binding sites approximately 45 angstroms apart. The modes of substrate binding have been deduced by analogy to D-Ala-D-Ala ligase and to pyruvate kinase. Coupling between the two remote active sites is facilitated by two conformational states of the phosphohistidine domain. While the crystal structure represents the state of interaction with the nucleotide, the second state is achieved by swiveling around two flexible peptide linkers. This dramatic conformational transition brings the phosphocarrier residue in close proximity to phosphoenolpyruvate/pyruvate. The swiveling-domain paradigm provides an effective mechanism for communication in complex multidomain/multiactive site proteins.