<i>Pcetk</i>: A pDynamo-based Toolkit for Protonation State Calculations in Proteins

<i>Pcetk</i> (a pDynamo-based continuum electrostatic toolkit) is an open-source, object-oriented toolkit for the calculation of proton binding energetics in proteins. The toolkit is a module of the pDynamo software library, combining the versatility of the Python scripting language and the efficiency of the compiled languages, C and Cython. In the toolkit, we have connected pDynamo to the external Poisson–Boltzmann solver, extended-MEAD. Our goal was to provide a modern and extensible environment for the calculation of protonation states, electrostatic energies, titration curves, and other electrostatic-dependent properties of proteins. <i>Pcetk</i> is freely available under the CeCILL license, which is compatible with the GNU General Public License. The toolkit can be found on the Web at the address http://github.com/mfx9/pcetk. The calculation of protonation states in proteins requires a knowledge of p<i>K</i>_a values of protonatable groups in aqueous solution. However, for some groups, such as protonatable ligands bound to protein, the p<i>K</i>_a^aq values are often difficult to obtain from experiment. As a complement to <i>Pcetk</i>, we revisit an earlier computational method for the estimation of p<i>K</i>_a^aq values that has an accuracy of ±0.5 p<i>K</i>_a-units or better. Finally, we verify the <i>Pcetk</i> module and the method for estimating p<i>K</i>_a^aq values with different model cases

Glycerol Dehydratation by the <i>B</i>₁₂-Independent Enzyme May Not Involve the Migration of a Hydroxyl Group: A Computational Study

A combination of continuum electrostatic and density functional calculations has been employed to study the mechanism of the <i>B</i>₁₂-independent glycerol dehydratase, a novel glycyl-radical enzyme involved in the microbial conversion of glycerol to 3-hydroxylpropionaldehyde. The calculations indicate that the dehydratation of glycerol by the <i>B</i>₁₂-independent enzyme does not need to involve a mechanistically complicated migration of the middle hydroxyl group to one of the two terminal positions of a molecule, as previously suggested. Instead, the reaction can proceed in three elementary steps. First, a radical transfer from the catalytically active Cys433 to the ligand generates a substrate-related intermediate. Second, a hydroxyl group splits off at the middle position of the ligand and is protonated by the neighboring His164 to form a water molecule. The other active site residue Glu435 accepts a proton from one of the terminal hydroxyl groups of the ligand and a CO double bond is created. Third, the reaction is completed by a radical back transfer from the product-related intermediate to Cys433. On the basis of our calculations, the catalytic functions of the active site residues have been suggested. Cys433 is a radical relay site; His164 and Glu435 make up a proton accepting/donating system; Asn156, His281, and Asp447 form a network of hydrogen bonds responsible for the electrostatic stabilization of the transition state. A synergistic participation of these residues in the reaction seems to be crucial for the catalysis

Catalytic mechanism of peptidoglycan deacetylase: A computational study

Bacterial peptidoglycan deacetylase enzymes are potentially important targets for the design of new drugs. In pathogenic bacteria, they modify cell-wall peptidoglycan by removing the acetyl group, which makes the bacteria more resistant to the host's immune response and other forms of attack, such as degradation by lysozyme. In this study, we have investigated the mechanism of reaction of acetyl removal from a model substrate, the N-acetylglucosamine/N-acetylmuramic acid dimer, by peptidogylcan deacetylase from Helicobacter pylori. For this, we employed a range of computational approaches, including molecular docking, Poisson-Boltzmann electrostatic pKacalculations, molecular dynamics simulations, and hybrid quantum chemical/molecular mechanical potential calculations, in conjunction with reaction-path-finding algorithms. The active site of this enzyme is in a region of highly negative electrostatic potential and contains a zinc dication with a bound water molecule. In the docked enzyme-substrate complex, our pKa calculations indicate that in the most stable protonation states of the active site the zinc-bound water molecule is in its hydroxide form and that the adjacent histidine residue, His247, is doubly protonated. In addition, there are one or two excess protons, with the neighboring aspartate residues, Asp12 and/or Asp199, being protonated. Overall, we find five classes of feasible reaction mechanisms, with the favored mechanism depending heavily on the protonation state of the active site. In the major one-excess-proton form, the mechanism with the lowest barrier (84 kJ mol-1) involves an initial protonation of the substrate nitrogen, followed by nucleophilic attack of the zinc-bound hydroxide and rupture of the substrate's carbon-nitrogen bond. However, in the minor two-excess-proton form, four mechanisms are almost equienergetic (83-86 kJ mol-1), comprising both those that start with nitrogen protonation and those in which nucleophilic attack by hydroxide occurs first

Structural Determinants of Improved Fluorescence in a Family of Bacteriophytochrome-Based Infrared Fluorescent Proteins: Insights from Continuum Electrostatic Calculations and Molecular Dynamics Simulations

International audienceUsing X-ray crystallography, continuum electrostatic calculations, and molecular dynamics simulations, we have studied the structure, protonation behavior, and dynamics of the biliverdin chromophore and its molecular environment in a series of genetically engineered infrared fluorescent proteins (IFPs) based on the chromophore-binding domain of the Deinococcus radiodurans bacteriophytochrome. Our study suggests that the experimentally observed enhancement of fluorescent properties results from the improved rigidity and planarity of the biliverdin chromophore, in particular of the first two pyrrole rings neighboring the covalent linkage to the protein. We propose that the increases in the levels of both motion and bending of the chromophore out of planarity favor the decrease in fluorescence. The chromophore-binding pocket in some of the studied proteins, in particular the weakly fluorescent parent protein, is shown to be readily accessible to water molecules from the solvent. These waters entering the chromophore region form hydrogen bond networks that affect the otherwise planar conformation of the first three rings of the chromophore. On the basis of our simulations, the enhancement of fluorescence in IFPs can be achieved either by reducing the mobility of water molecules in the vicinity of the chromophore or by limiting the interactions of the nearby protein residues with the chromophore. Finally, simulations performed at both low and neutral pH values highlight differences in the dynamics of the chromophore and shed light on the mechanism of fluorescence loss at low pH

Catalytic Mechanism of the Glycyl Radical Enzyme 4‑Hydroxyphenylacetate Decarboxylase from Continuum Electrostatic and QC/MM Calculations

Using continuum electrostatics and QC/MM calculations, we investigate the catalytic cycle of the glycyl radical enzyme 4-hydroxyphenylacetate decarboxylase, an enzyme involved in the fermentative production of <i>p</i>-cresol from tyrosine in clostridia. On the basis of our calculations, we propose a five-step mechanism for the reaction. In the first step, the substrate 4-hydroxyphenylacetate is activated by an unusual concerted abstraction of an electron and a proton. Namely, Cys503 radical abstracts an electron from the substrate and Glu637 abstracts a proton. Thus in total, a hydrogen atom is abstracted from the substrate. In the second step, the carboxylic group readily splits off from the phenoxy-acetate radical anion to give carbon dioxide. This decarboxylation step is coupled to a proton transfer from Glu637 back to the phenolic hydroxyl group which leads to a <i>p</i>-hydroxybenzyl radical. The remaining steps of the reaction involve a rotation of the Cys503 side chain followed by a proton transfer from Glu505 to Cys503 and a hydrogen atom transfer from Cys503 to the <i>p</i>-hydroxybenzyl radical to form <i>p</i>-cresol. The calculated mechanism agrees with experimental data suggesting that both Cys503 and Glu637 are essential for the catalytic function of 4-hydroxyphenylacetate decarboxylase and that the substrate requires a hydroxyl group in <i>para</i>-position to the acetate moiety

Catalytic Mechanism of the Glycyl Radical Enzyme 4‑Hydroxyphenylacetate Decarboxylase from Continuum Electrostatic and QC/MM Calculations

Using continuum electrostatics and QC/MM calculations, we investigate the catalytic cycle of the glycyl radical enzyme 4-hydroxyphenylacetate decarboxylase, an enzyme involved in the fermentative production of <i>p</i>-cresol from tyrosine in clostridia. On the basis of our calculations, we propose a five-step mechanism for the reaction. In the first step, the substrate 4-hydroxyphenylacetate is activated by an unusual concerted abstraction of an electron and a proton. Namely, Cys503 radical abstracts an electron from the substrate and Glu637 abstracts a proton. Thus in total, a hydrogen atom is abstracted from the substrate. In the second step, the carboxylic group readily splits off from the phenoxy-acetate radical anion to give carbon dioxide. This decarboxylation step is coupled to a proton transfer from Glu637 back to the phenolic hydroxyl group which leads to a <i>p</i>-hydroxybenzyl radical. The remaining steps of the reaction involve a rotation of the Cys503 side chain followed by a proton transfer from Glu505 to Cys503 and a hydrogen atom transfer from Cys503 to the <i>p</i>-hydroxybenzyl radical to form <i>p</i>-cresol. The calculated mechanism agrees with experimental data suggesting that both Cys503 and Glu637 are essential for the catalytic function of 4-hydroxyphenylacetate decarboxylase and that the substrate requires a hydroxyl group in <i>para</i>-position to the acetate moiety

Catalytic Mechanism of the Glycyl Radical Enzyme 4‑Hydroxyphenylacetate Decarboxylase from Continuum Electrostatic and QC/MM Calculations

Using continuum electrostatics and QC/MM calculations, we investigate the catalytic cycle of the glycyl radical enzyme 4-hydroxyphenylacetate decarboxylase, an enzyme involved in the fermentative production of <i>p</i>-cresol from tyrosine in clostridia. On the basis of our calculations, we propose a five-step mechanism for the reaction. In the first step, the substrate 4-hydroxyphenylacetate is activated by an unusual concerted abstraction of an electron and a proton. Namely, Cys503 radical abstracts an electron from the substrate and Glu637 abstracts a proton. Thus in total, a hydrogen atom is abstracted from the substrate. In the second step, the carboxylic group readily splits off from the phenoxy-acetate radical anion to give carbon dioxide. This decarboxylation step is coupled to a proton transfer from Glu637 back to the phenolic hydroxyl group which leads to a <i>p</i>-hydroxybenzyl radical. The remaining steps of the reaction involve a rotation of the Cys503 side chain followed by a proton transfer from Glu505 to Cys503 and a hydrogen atom transfer from Cys503 to the <i>p</i>-hydroxybenzyl radical to form <i>p</i>-cresol. The calculated mechanism agrees with experimental data suggesting that both Cys503 and Glu637 are essential for the catalytic function of 4-hydroxyphenylacetate decarboxylase and that the substrate requires a hydroxyl group in <i>para</i>-position to the acetate moiety

Norovirus RNA-dependent RNA polymerase: A computational study of metal-binding preferences

Norovirus (NV) RNA-dependent RNA polymerase (RdRP) is essential for replicating the genome of the virus, which makes this enzyme a key target for the development of antiviral agents against NV gastroenteritis. In this work, a complex of NV RdRP bound to manganese ions and an RNA primer-template duplex was investigated using X-ray crystallography and hybrid quantum chemical/molecular mechanical simulations. Experimentally, the complex crystallized in a tetragonal crystal form. The nature of the primer/template duplex binding in the resulting structure indicates that the complex is a closed back-tracked state of the enzyme, in which the (Formula presented.) -end of the primer occupies the position expected for the post-incorporated nucleotide before translocation. Computationally, it is found that the complex can accept a range of divalent metal cations without marked distortions in the active site structure. The highest binding energy is for copper, followed closely by manganese and iron, and then by zinc, nickel, and cobalt. Proteins 2017; 85:1435\ue2\u80\u931445. \uc2\ua9 2017 Wiley Periodicals, Inc