Photoswitching mechanism of a fluorescent protein revealed by time-resolved crystallography and transient absorption spectroscopy.

Reversibly switchable fluorescent proteins (RSFPs) serve as markers in advanced fluorescence imaging. Photoswitching from a non-fluorescent off-state to a fluorescent on-state involves trans-to-cis chromophore isomerization and proton transfer. Whereas excited-state events on the ps timescale have been structurally characterized, conformational changes on slower timescales remain elusive. Here we describe the off-to-on photoswitching mechanism in the RSFP rsEGFP2 by using a combination of time-resolved serial crystallography at an X-ray free-electron laser and ns-resolved pump-probe UV-visible spectroscopy. Ten ns after photoexcitation, the crystal structure features a chromophore that isomerized from trans to cis but the surrounding pocket features conformational differences compared to the final on-state. Spectroscopy identifies the chromophore in this ground-state photo-intermediate as being protonated. Deprotonation then occurs on the µs timescale and correlates with a conformational change of the conserved neighbouring histidine. Together with a previous excited-state study, our data allow establishing a detailed mechanism of off-to-on photoswitching in rsEGFP2.

On the control of the proton current in the voltage-gated proton channel Hv1.

The nature of the action of voltage-activated proton transport proteins is a conundrum of great current interest. Here we approach this issue by exploring the action of Hv1, a voltage-gated proton channel found in different cells in humans and other organisms. Our study focuses on evaluating the free energy of transporting a proton through the channel, as well as the effect of the proton transfer through D112, in both the closed and open channel conformations. It is found that D112 allows a transported proton to bypass the electrostatic barrier of the open channel, while not being able to help in passing the barrier in the closed form. This reflects the change in position of the gating arginine residues relative to D112, upon voltage activation. Significantly, the effect of D112 accounts for the observed trend in selectivity by overcoming the electrostatic barrier at its highest point. Thus, the calculations provide a structure/function correlation for the Hv1 system. The present work also clarifies that the action of Hv1 is not controlled by a Grotthuss mechanism but, as is always the case, by the protein electrostatic potential at the rate-limiting barriers.

Chromophore twisting in the excited state of a photoswitchable fluorescent protein captured by time-resolved serial femtosecond crystallography.

Chromophores absorb light in photosensitive proteins and thereby initiate fundamental biological processes such as photosynthesis, vision and biofluorescence. An important goal in their understanding is the provision of detailed structural descriptions of the ultrafast photochemical events that they undergo, in particular of the excited states that connect chemistry to biological function. Here we report on the structures of two excited states in the reversibly photoswitchable fluorescent protein rsEGFP2. We populated the states through femtosecond illumination of rsEGFP2 in its non-fluorescent off state and observed their build-up (within less than one picosecond) and decay (on the several picosecond timescale). Using an X-ray free-electron laser, we performed picosecond time-resolved crystallography and show that the hydroxybenzylidene imidazolinone chromophore in one of the excited states assumes a near-canonical twisted configuration halfway between the trans and cis isomers. This is in line with excited-state quantum mechanics/molecular mechanics and classical molecular dynamics simulations. Our new understanding of the structure around the twisted chromophore enabled the design of a mutant that displays a twofold increase in its off-to-on photoswitching quantum yield.

Predicting substituent effects on activation energy changes by static catalytic fields.

Catalytic fields illustrate topology of the optimal charge distribution of a molecular environment reducing the activation energy for any process involving barrier crossing, like chemical reaction, bond rotation etc. Until now, this technique has been successfully applied to predict catalytic effects resulting from intermolecular interactions with individual water molecules constituting the first hydration shell, aminoacid mutations in enzymes or Si→Al substitutions in zeolites. In this contribution, hydrogen to fluorine (H→F) substitution effects for two model reactions have been examined indicating qualitative applicability of the catalytic field concept in the case of systems involving intramolecular interactions. Graphical abstract Hydrogen to fluorine (H→F) substitution effects on activation energy in [kcal/mol].

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 3'-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-1445. © 2017 Wiley Periodicals, Inc.

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 pKa calculations, 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.

Using 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.

Pcetk: A pDynamo-based Toolkit for Protonation State Calculations in Proteins.

Pcetk (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. Pcetk 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 pKa values of protonatable groups in aqueous solution. However, for some groups, such as protonatable ligands bound to protein, the pKa aq values are often difficult to obtain from experiment. As a complement to Pcetk, we revisit an earlier computational method for the estimation of pKa aq values that has an accuracy of ± 0.5 pKa-units or better. Finally, we verify the Pcetk module and the method for estimating pKa aq values with different model cases.

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 p-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 p-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 p-hydroxybenzyl radical to form p-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 para-position to the acetate moiety.

Low cost prediction of relative stabilities of hydrogen bonded complexes from atomic multipole moments for overly short intermolecular distances.

The relative stability of biologically relevant, hydrogen bonded complexes with shortened distances can be assessed at low cost by the electrostatic multipole term alone more successfully than by ab initio methods. These results imply that atomic multipole moments may help improve ligand-receptor ranking predictions, particularly in cases where accurate structural data are not available.

Glycerol dehydratation by the B12-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 B(12)-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 B(12)-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.

Structural basis for a Kolbe-type decarboxylation catalyzed by a glycyl radical enzyme.

4-Hydroxyphenylacetate decarboxylase is a [4Fe-4S] cluster containing glycyl radical enzyme proposed to use a glycyl/thiyl radical dyad to catalyze the last step of tyrosine fermentation in clostridia. The decarboxylation product p-cresol (4-methylphenol) is a virulence factor of the human pathogen Clostridium difficile . Here we describe the crystal structures at 1.75 and 1.81 Å resolution of substrate-free and substrate-bound 4-hydroxyphenylacetate decarboxylase from the related Clostridium scatologenes . The structures show a (ßγ)(4) tetramer of heterodimers composed of a catalytic ß-subunit harboring the putative glycyl/thiyl dyad and a distinct small γ-subunit with two [4Fe-4S] clusters at 40 Å distance from the active site. The γ-subunit comprises two domains displaying pseudo-2-fold symmetry that are structurally related to the [4Fe-4S] cluster-binding scaffold of high-potential iron-sulfur proteins. The N-terminal domain coordinates one cluster with one histidine and three cysteines, and the C-terminal domain coordinates the second cluster with four cysteines. Whereas the C-terminal cluster is buried in the ßγ heterodimer interface, the N-terminal cluster is not part of the interface. The previously postulated decarboxylation mechanism required the substrate's hydroxyl group in the vicinity of the active cysteine residue. In contrast to expectation, the substrate-bound state shows a direct interaction between the substrate's carboxyl group and the active site Cys503, while His536 and Glu637 at the opposite side of the active site pocket anchor the hydroxyl group. This state captures a possible catalytically competent complex and suggests a Kolbe-type decarboxylation for p-cresol formation.

The molecular basis of urokinase inhibition: from the nonempirical analysis of intermolecular interactions to the prediction of binding affinity.

Urokinase-type plasminogen activator (uPA) is a trypsin-like serine protease that plays a crucial role in angiogenesis process. In addition to its physiological role in healthy organisms, angiogenesis is extremely important in cancer growth and metastasis, resulting in numerous attempts to understand its control and to develop new approaches to anticancer therapy. The alpha-aminoalkylphosphonate diphenyl esters are well known as highly efficient serine protease inhibitors. However, their mode of binding has not been verified experimentally in details. For a group of average and potent phosphonic inhibitors of urokinase, flexible docking calculations were performed to gain an insight into the active site interactions responsible for observed enzyme inhibition. The docking results are consistent with the previously suggested mode of inhibitors binding. Subsequently, rigorous ab initio study of binding energy was carried out, followed by its decomposition according to the variation-perturbation procedure to reveal stabilization energy constituents with clear physical meaning. Availability of the experimental inhibitory activities and comparison with theoretical binding energy allows for the validation of theoretical models of inhibition, as well as estimation of the possible potential for binding affinity prediction. Since the docking results accompanied by molecular mechanics optimization suggested that several crucial active site contacts were too short, the optimal distances corresponding to the minimum ab initio interaction energy were also evaluated. Despite the deficiencies of force field-optimized enzyme-inhibitor structures, satisfactory agreement with experimental inhibitory activity was obtained for the electrostatic interaction energy, suggesting its possible application in the binding affinity prediction.

Oxime-induced reactivation of sarin-inhibited AChE: a theoretical mechanisms study.

Oximes (especially oximate anions) are used as potential reactivators of OP-inhibited AChE due to their unique alpha-effect nucleophilic reactivity. In the present study, by applying the DFT approach at the B3LYP/6-311G(d,p) level and the Møller-Plesset perturbation theory at the MP2/6-311G(d,p) level, the formoximate-induced reactivation patterns of the sarin-AChE adduct and the corresponding reaction mechanism have been investigated. The potential energy surface along the pathway of the reactivation reaction of sarin-inhibited AChE by oxime reveals that the reaction can occur quickly due to the relatively low energy barriers. A two-step process is a major pathway proposed for the studied reactivation reaction. Through the nucleophilic attack, the oximate first binds to the sarin-AChE adduct to form a relatively stable phosphorus complex. The regeneration of the serine takes place subsequently through an elimination step, which is expected to be competitive with the nucleophilic attacking process. The polarizable continuum model (PCM) has been applied to evaluate the solvate effects on the pathway. It is concluded that the reaction energy barriers are also low enough for the reaction to easily occur in solvent. The results derived from both the gas-phase model and the aqueous solvation model suggest that the studied oximate anion is an efficient antidote reagent for sarin-inhibited AChE.