Energetics of nucleophile activation in a protein tyrosine phosphatase

The exposure of OPFRs in fish from aquaculture area: Backward tracing of the ecological risk regulation

In this study, we backward traced and controlled the pollution of organophosphorus flame retardants (OPFRs) in aquaculture areas from the standpoints of terminal treatment, migration and transformation resistance, and source molecular substitution technology. A regulatory plan to considerably reduce the combined biotoxicity of fish exposed to OPFRs in aquaculture areas and significantly improves the biodegradation of sewage treatment and the efficiency of soil plant-microorganism combined remediation was formulated. Environmentally friendly alternatives of OPFRs were designed. The supplementation scheme of aquatic feed significantly alleviates the toxicity risk of fish exposure to OPFRs in aquafarm (reduced by 121.02%). The regulatory scheme of external stimulus to enhance the biodegradation of OPFRs in wastewater treatment process included an H₂O₂ concentration of 400 mg/L, voltage gradient of 1.5 V/m, and pH of 6.5 can improve the degradation capacity of OPFRs molecules by 88.86%. The degradation of OPFRs can be enhanced by plant-microorganism combined remediation (up to 98.64%) by growing plants whose primary function is phytoextraction in soils dominated by Sphingopyxis sp. and Rhodococcus sp. A 3D-QSAR pharmacophore model based on apoptosis toxicity, mitochondrial dysfunction, oxidative stress response, reproductive, neurotoxicity, gill-inhalation combined toxicity of fish exposed to OPFRs in aquafarm was fabricated. The recommended aquatic feed scheme and the control scheme of enhanced degradation of OPFRs by sewage treatment and soil environment had better applicability for the new-designed OPFRs substitution molecules (the maximum combined toxicity/degradation is reduced/increased by 75.46% and 63.24%, respectively). In this paper, a technical scheme of OPFRs terminal treatment, process regulation, and source control was applied as a cradle-to-grave approach to reduce the ecological toxicity risk of fish exposed to OPFRs in aquaculture areas providing theoretical support for the realization of OPFRs environmental pollution control.

Fine-Tuning of Sequence Specificity by Near Attack Conformations in Enzyme-Catalyzed Peptide Hydrolysis

The catalytic role of near attack conformations (NACs), molecular states that lie on the pathway between the ground state (GS) and transition state (TS) of a chemical reaction, is not understood completely. Using a computational approach that combines Bürgi–Dunitz theory with all-atom molecular dynamics simulations, the role of NACs in catalyzing the first stages of HIV-1 protease peptide hydrolysis was previously investigated using a substrate that represents the recognized SP1-NC cleavage site of the HIV-1 Gag polyprotein. NACs were found to confer no catalytic effect over the uncatalyzed reaction there ( Δ Δ G N ‡ ∼ 0 kcal/mol). Here, using the same approach, the role of NACs across multiple substrates that each represent a further recognized cleavage site is investigated. Overall rate enhancement varies by | Δ Δ G ‡ | ∼ 12–15 kcal/mol across this set, and although NACs contribute a small and approximately constant barrier to the uncatalyzed reaction (< Δ G N ‡ u > = 4.3 ± 0.3 kcal/mol), they are found to contribute little significant catalytic effect ( | Δ Δ G N ‡ | ∼ 0–2 kcal/mol). Furthermore, no correlation is exhibited between NAC contributions and the overall energy barrier ( R 2 = 0.01). However, these small differences in catalyzed NAC contributions enable rates to match those required for the kinetic order of processing. Therefore, NACs may offer an alternative and subtle mode compared to non-NAC contributions for fine-tuning reaction rates during complex evolutionary sequence selection processes—in this case across cleavable polyproteins whose constituents exhibit multiple functions during the virus life-cycle.

AMBER and CHARMM Force Fields Inconsistently Portray the Microscopic Details of Phosphorylation

Phosphorylation of serine, threonine, and tyrosine is one of the most frequently occurring and crucial post-translational modifications of proteins often associated with important structural and functional changes. We investigated the direct effect of phosphorylation on the intrinsic conformational preferences of amino acids as a potential trigger of larger structural events. We conducted a comparative study of force fields on terminally capped amino acids (dipeptides) as the simplest model for phosphorylation. Our bias-exchange metadynamics simulations revealed that all model dipeptides sampled a great heterogeneity of ensembles affected by introduction of mono- and dianionic phosphate groups. However, the detected changes in populations of backbone conformers and side-chain rotamers did not reveal a strong discriminatory shift in preferences, as could be anticipated for the bulky, charged phosphate group. Furthermore, the AMBER and CHARMM force fields provided inconsistent populations of individual conformers as well as net structural trends upon phosphorylation. Detailed analysis of ensembles revealed competition between hydration and formation of internal hydrogen bonds involving amide hydrogens and the phosphate group. The observed difference in hydration free energy and potential for hydrogen bonding in individual force fields could be attributed to the different partial atomic charges used in each force field and, hence, the different parametrization strategies. Nevertheless, conformational propensities and net structural changes upon phosphorylation are difficult to extract from experimental measurements, and existing experimental data provide limited guidance for force field assessment and further development.

Vibrio cholerae LMWPTP-2 display unique surface charge and grooves around the active site: Indicative of distinctive substrate specificity and scope to design specific inhibitor

Low molecular weight protein tyrosine phosphatases (LMWPTPs) are ubiquitously found as small cytoplasmic enzymes which act on phospho-tyrosine containing proteins that are engaged in various cellular functions. Vibrio cholerae O395 contains two LMWPTPs having widely different sequence. Phylogenetic analysis based on a non redundant set of 124 LMWPTP sequences, designate that LMWPTP-2 from Vibrio choleraeO395 (VcLMWPTP-2) is a single taxon. We have determined the crystal structure of VcLMWPTP-2 at 2.6 Å with MOPS bound in the active site. Tertiary structure analysis indicates that VcLMWPTP-2 forms dimer. Studies in solution state also confirm exclusive presence of a dimeric form. Kinetic studies demonstrate that VcLMWPTP-2 dimer is catalytically active while inactivation through oligomerisation was reported as one of the regulatory mechanism in case of mammalian LMWPTP viz., Bos taurus LMWPTP, BPTP. Kinetic studies using p-nitrophenyl phosphate (p-NPP) as a substrate demonstrate active participation of both the P-loop cysteine in catalysis. Vicinal Cys17, in addition plays a role of protecting the catalytic Cys12 under oxidative stress. Structural analysis and MD simulations allowed us to propose the role of several conserved residues around the active site. Distribution of surface charges and grooves around the active site delineates unique features of VcLMWPTP-2 which could be utilized to design specific inhibitor.

Crystal structure of undecaprenyl-pyrophosphate phosphatase and its role in peptidoglycan biosynthesis

As a protective envelope surrounding the bacterial cell, the peptidoglycan sacculus is a site of vulnerability and an antibiotic target. Peptidoglycan components, assembled in the cytoplasm, are shuttled across the membrane in a cycle that uses undecaprenyl-phosphate. A product of peptidoglycan synthesis, undecaprenyl-pyrophosphate, is converted to undecaprenyl-phosphate for reuse in the cycle by the membrane integral pyrophosphatase, BacA. To understand how BacA functions, we determine its crystal structure at 2.6 Å resolution. The enzyme is open to the periplasm and to the periplasmic leaflet via a pocket that extends into the membrane. Conserved residues map to the pocket where pyrophosphorolysis occurs. BacA incorporates an interdigitated inverted topology repeat, a topology type thus far only reported in transporters and channels. This unique topology raises issues regarding the ancestry of BacA, the possibility that BacA has alternate active sites on either side of the membrane and its possible function as a flippase.

Bacterial cell wall components are assembled in a transmembrane cycle that involves the membrane integral pyrophosphorylase, BacA. Here the authors solve the crystal structure of BacA which shows an interdigitated inverted topology repeat that hints towards a flippase function for BacA.

Fortilin binds IRE1α and prevents ER stress from signaling apoptotic cell death

The endoplasmic reticulum, the cytoplasmic organelle that matures a massive amount of nascent secretory polypeptides, is particularly sensitive to stress. Endoplasmic reticulum stress causes unfolded proteins to populate the organelle, eliciting the unfolded protein response. During the unfolded protein response, GRP78—an endoplasmic reticulum master stress regulator—detaches from three endoplasmic reticulum stress sensors (IRE1α, PERK, and ATF6) and allows them to activate the apoptotic signaling pathway. Fortilin, a pro-survival molecule, is known to inhibit apoptosis by binding and inhibiting p53, but its role in endoplasmic reticulum stress-induced apoptosis remains unknown. Here, we report that fortilin directly interacts with the cytoplasmic domain of IRE1α, inhibits both kinase and endoribonuclease (RNase) activities of the stress sensor, and protects cells against apoptotic cell death at both cellular and whole animal levels. Our data support a role of fortilin in the unfolded protein response and its potential participation in human diseases caused by unfolded protein response.

IRE1α is an ER stress sensor, whose activity induces apoptosis. Here, the authors report that fortilin, a pro-survival factor, with yet unknown roles in ER stress, interacts with active IRE1α, inhibits both its kinase end RNase activities, and protects cells from apoptosis both in vitro and in vivo.

Low molecular weight protein tyrosine phosphatase: Multifaceted functions of an evolutionarily conserved enzyme

Originally identified as a low molecular weight acid phosphatase, LMW-PTP is actually a protein tyrosine phosphatase that acts on many phosphotyrosine-containing cellular proteins that are primarily involved in signal transduction. Differences in sequence, structure, and substrate recognition as well as in subcellular localization in different organisms enable LMW-PTP to exert many different functions. In fact, during evolution, the LMW-PTP structure adapted to perform different catalytic actions depending on the organism type. In bacteria, this enzyme is involved in the biosynthesis of group 1 and 4 capsules, but it is also a virulence factor in pathogenic strains. In yeast, LMW-PTPs dephosphorylate immunophilin Fpr3, a peptidyl-prolyl-cis-trans isomerase member of the protein chaperone family. In humans, LMW-PTP is encoded by the ACP1 gene, which is composed of three different alleles, each encoding two active enzymes produced by alternative RNA splicing. In animals, LMW-PTP dephosphorylates a number of growth factor receptors and modulates their signalling processes. The involvement of LMW-PTP in cancer progression and in insulin receptor regulation as well as its actions as a virulence factor in a number of pathogenic bacterial strains may promote the search for potent, selective and bioavailable LMW-PTP inhibitors.

Hydroxyapatite Growth Inhibition Effect of Pellicle Statherin Peptides

In our recent studies, we have shown that in vivo–acquired enamel pellicle is a sophisticated biological structure containing a significant portion of naturally occurring salivary peptides. From a functional aspect, the identification of peptides in the acquired enamel pellicle is of interest because many salivary proteins exhibit functional domains that maintain the activities of the native protein. Among the in vivo–acquired enamel pellicle peptides that have been newly identified, 5 peptides are derived from statherin. Here, we assessed the ability of these statherin pellicle peptides to inhibit hydroxyapatite crystal growth. In addition, atomistic molecular dynamics (MD) simulations were performed to better understand the underlying physical mechanisms of hydroxyapatite growth inhibition. A microplate colorimetric assay was used to quantify hydroxyapatite growth. Statherin protein, 5 statherin-derived peptides, and a peptide lacking phosphate at residues 2 and 3 were analyzed. Statherin peptide phosphorylated on residues 2 and 3 indicated a significant inhibitory effect when compared with the 5 other peptides ( P < 0.05). MD simulations showed a strong affinity and fast adsorption to hydroxyapatite for phosphopeptides, whereas unphosphorylated peptides interacted weakly with the hydroxyapatite. Our data suggest that the presence of a covalently linked phosphate group (at residues 2 and 3) in statherin peptides modulates the effect of hydroxyapatite growth inhibition. This study provides a mechanism to account for the composition and function of acquired enamel pellicle statherin peptides that will contribute as a base for the development of biologically stable and functional synthetic peptides for therapeutic use against dental caries and/or periodontal disease.

Mutations Impairing GSK3-Mediated MAF Phosphorylation Cause Cataract, Deafness, Intellectual Disability, Seizures, and a Down Syndrome-like Facies

Transcription factors operate in developmental processes to mediate inductive events and cell competence, and perturbation of their function or regulation can dramatically affect morphogenesis, organogenesis, and growth. We report that a narrow spectrum of amino-acid substitutions within the transactivation domain of the v-maf avian musculoaponeurotic fibrosarcoma oncogene homolog (MAF), a leucine zipper-containing transcription factor of the AP1 superfamily, profoundly affect development. Seven different de novo missense mutations involving conserved residues of the four GSK3 phosphorylation motifs were identified in eight unrelated individuals. The distinctive clinical phenotype, for which we propose the eponym Aymé-Gripp syndrome, is not limited to lens and eye defects as previously reported for MAF/Maf loss of function but includes sensorineural deafness, intellectual disability, seizures, brachycephaly, distinctive flat facial appearance, skeletal anomalies, mammary gland hypoplasia, and reduced growth. Disease-causing mutations were demonstrated to impair proper MAF phosphorylation, ubiquitination and proteasomal degradation, perturbed gene expression in primary skin fibroblasts, and induced neurodevelopmental defects in an in vivo model. Our findings nosologically and clinically delineate a previously poorly understood recognizable multisystem disorder, provide evidence for MAF governing a wider range of developmental programs than previously appreciated, and describe a novel instance of protein dosage effect severely perturbing development.

Structures, dynamics, and water permeation free energy across bilayers of Lipid A and its analog studied with molecular dynamics simulation

Fundamental studies of the supramolecular layer structures, dynamics and water permeation free energy of hexa-acyl-chain Lipid A and its analogue of tetra-acyl chains would be useful for polymer membranes design for endotoxin removal in water treatment, drug delivery and other biotechnologies. In this work, we studied their supramolecular bilayer by using molecular dynamics simulations and efficient free energy computations. Our simulation accuracy was verified by the agreement between the bilayer structural properties (structure factor, bilayer thickness, and the area per lipid) and lateral diffusion coefficient in our simulation and experimental measurements. More importantly, our simulation for the first time illustrated hexagonal compact packing of the hydrocarbon acyl chains within a leaflet of Lipid A membrane (at 298 K and water content of 40 wt %), which is consistent with experiments. In contrast, Lipid A analogue is found with less ordered ripple structures at the same condition. Our study also demonstrated slower dynamics and larger and broader free energy barrier (∼23 kJ/mol) for water permeation for Lipid A, compared with that of Lipid A analogue. Moreover, the analysis of dynamics showed that highly hydrated hydrophilic diglucosamine backbone is structurally stable, whereas the interdigitated hydrophobic acyl chain tails inside the membrane with faster dynamics screen the aqueous environment from the lipid interior and also reinforce the membrane's structural stability.

Binding of disordered proteins to a protein hub

A small number of proteins, called hubs, have high connectivity and are essential for interactome functionality and integrity. Keap1 is a crucial hub in the oxidative stress response and apoptosis. The Kelch domain of Keap1 preferentially binds to disordered regions of its partners, which share similar binding motifs, but have a wide range of binding affinities. Isothermal titration calorimetry (ITC) and multi-microsecond molecular dynamics (MD) simulations were used to determine the factors that govern the affinity of all currently known disordered binding partners to Kelch. Our results show that the affinities to this hub are largely determined by the extent of preformed bound state-like conformation in the free state structures of these disordered targets. Based on our findings, we have designed a high-affinity peptide that can specifically disrupt the Keap1-NRF2 interaction and has the potential for therapeutic applications.

Identification of multiple phosphorylation sites on maize endosperm starch branching enzyme IIb, a key enzyme in amylopectin biosynthesis

Starch branching enzyme IIb (SBEIIb) plays a crucial role in amylopectin biosynthesis in maize endosperm by defining the structural and functional properties of storage starch and is regulated by protein phosphorylation. Native and recombinant maize SBEIIb were used as substrates for amyloplast protein kinases to identify phosphorylation sites on the protein. A multidisciplinary approach involving bioinformatics, site-directed mutagenesis, and mass spectrometry identified three phosphorylation sites at Ser residues: Ser(649), Ser(286), and Ser(297). Two Ca(2+)-dependent protein kinase activities were partially purified from amyloplasts, termed K1, responsible for Ser(649) and Ser(286) phosphorylation, and K2, responsible for Ser(649) and Ser(297) phosphorylation. The Ser(286) and Ser(297) phosphorylation sites are conserved in all plant branching enzymes and are located at opposite openings of the 8-stranded parallel β-barrel of the active site, which is involved with substrate binding and catalysis. Molecular dynamics simulation analysis indicates that phospho-Ser(297) forms a stable salt bridge with Arg(665), part of a conserved Cys-containing domain in plant branching enzymes. Ser(649) conservation appears confined to the enzyme in cereals and is not universal, and is presumably associated with functions specific to seed storage. The implications of SBEIIb phosphorylation are considered in terms of the role of the enzyme and the importance of starch biosynthesis for yield and biotechnological application.

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPDclosed ) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPDopen ) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPDclosed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPDopen conformation. When the active site water molecule network that is a part of the free WPDclosed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

Comparison of Secondary Structure Formation Using 10 Different Force Fields in Microsecond Molecular Dynamics Simulations

We have compared molecular dynamics (MD) simulations of a β-hairpin forming peptide derived from the protein Nrf2 with 10 biomolecular force fields using trajectories of at least 1 μs. The total simulation time was 37.2 μs. Previous studies have shown that different force fields, water models, simulation methods, and parameters can affect simulation outcomes. The MD simulations were done in explicit solvent with a 16-mer Nrf2 β-hairpin forming peptide using Amber ff99SB-ILDN, Amber ff99SB*-ILDN, Amber ff99SB, Amber ff99SB*, Amber ff03, Amber ff03*, GROMOS96 43a1p, GROMOS96 53a6, CHARMM27, and OPLS-AA/L force fields. The effects of charge-groups, terminal capping, and phosphorylation on the peptide folding were also examined. Despite using identical starting structures and simulation parameters, we observed clear differences among the various force fields and even between replicates using the same force field. Our simulations show that the uncapped peptide folds into a native-like β-hairpin structure at 310 K when Amber ff99SB-ILDN, Amber ff99SB*-ILDN, Amber ff99SB, Amber ff99SB*, Amber ff03, Amber ff03*, GROMOS96 43a1p, or GROMOS96 53a6 were used. The CHARMM27 simulations were able to form native hairpins in some of the elevated temperature simulations, while the OPLS-AA/L simulations did not yield native hairpin structures at any temperatures tested. Simulations that used charge-groups or peptide capping groups were not largely different from their uncapped counterparts with single atom charge-groups. On the other hand, phosphorylation of the threonine residue located at the β-turn significantly affected the hairpin formation. To our knowledge, this is the first study comparing such a large set of force fields with respect to β-hairpin folding. Such a comprehensive comparison will offer useful guidance to others conducting similar types of simulations.

Detailed regulatory mechanism of the interaction between ZO-1 PDZ2 and connexin43 revealed by MD simulations

The gap junction protein connexin43 (Cx43) binds to the second PDZ domain of Zonula occludens-1 (ZO-1) through its C-terminal tail, mediating the regulation of gap junction plaque size and dynamics. Biochemical study demonstrated that the very C-terminal 12 residues of Cx43 are necessary and sufficient for ZO-1 PDZ2 binding and phosphorylation at residues Ser (-9) and Ser (-10) of the peptide can disrupt the association. However, only a crystal structure of ZO-1 PDZ2 in complex with a shorter 9 aa peptide of connexin43 was solved experimentally. Here, the interactions between ZO-1 PDZ2 and the short, long and phosphorylated Cx43 peptides were studied using molecular dynamics (MD) simulations and free energy calculation. The short peptide bound to PDZ2 exhibits large structural variations, while the extension of three upstream residues stabilizes the peptide conformation and enhanced the interaction. Phosphorylation at Ser(-9) significantly weakens the binding and results in conformational flexibility of the peptide. Glu210 of ZO-1 PDZ2 was found to be a key regulatory point in Cx43 binding and phosphorylation induced dissociation.

Insights into the reaction of protein-tyrosine phosphatase 1B: crystal structures for transition state analogs of both catalytic steps

Catalysis by protein-tyrosine phosphatase 1B (PTP1B) occurs through a two-step mechanism involving a phosphocysteine intermediate. We have solved crystal structures for the transition state analogs for both steps. Together with previously reported crystal structures of apo-PTP1B, the Michaelis complex of an inactive mutant, the phosphoenzyme intermediate, and the product complex, a full picture of all catalytic steps can now be depicted. The transition state analog for the first catalytic step comprises a ternary complex between the catalytic cysteine of PTP1B, vanadate, and the peptide DADEYL, a fragment of a physiological substrate. The equatorial vanadate oxygen atoms bind to the P-loop, and the apical positions are occupied by the peptide tyrosine oxygen and by the PTP1B cysteine sulfur atom. The vanadate assumes a trigonal bipyramidal geometry in both transition state analog structures, with very similar apical O-O distances, denoting similar transition states for both phosphoryl transfer steps. Detailed interactions between the flanking peptide and the enzyme are discussed.

Enzyme millisecond conformational dynamics do not catalyze the chemical step

The idea that enzymes catalyze reactions by dynamical coupling between the conformational motions and the chemical coordinates has recently attracted major experimental and theoretical interest. However, experimental studies have not directly established that the conformational motions transfer energy to the chemical coordinate, and simulating enzyme catalysis on the relevant timescales has been impractical. Here, we introduce a renormalization approach that transforms the energetics and dynamics of the enzyme to an equivalent low-dimensional system, and allows us to simulate the dynamical coupling on a ms timescale. The simulations establish, by means of several independent approaches, that the conformational dynamics is not remembered during the chemical step and does not contribute significantly to catalysis. Nevertheless, the precise nature of this coupling is a question of great importance.

New insights into the activation of Escherichia coli tyrosine kinase revealed by molecular dynamics simulation and biochemical analysis

Escherichia coli tyrosine kinase (Etk) regulates the export of pathogenic capsular polysaccharide (CPS) by intermolecularly autophosphorylating its C-terminal tyrosine cluster. The kinase Etk, however, needs to be first activated by the intramolecular phosphorylation of a tyrosine residue, Y574, next to the active site. The recently determined structure of Etk shows that dephosphorylated Y574 blocks the active site and prevents substrate access. After phosphorylation, the negatively charged P-Y574 side chain was previously postulated to flip out to associate with a positively charged R614, unblocking the active site. This proposed activation is unique among protein kinases; however, there is no direct structural evidence in support of this hypothesis. In this paper, we carried out molecular dynamics simulation, mutagenesis, and biochemical analysis to study the activation mechanism of Etk. Our simulation results are in excellent agreement with the proposed molecular switch involving P-Y574 and R614 in the activation of Etk. Further, we show that a previously unidentified residue, R572, modulates the rotation of the P-Y574 side chain through electrostatic interaction, slowing down the opening of the active site. Our enzymatic assays demonstrate that the R572A mutant of Etk possesses significantly increased kinase activity, providing direct experimental support for the unique activation mechanism of Etk. In addition, the simulation of the Etk Y574F mutant predicted short periods of unblocked active site by Y574F, in good agreement with the low kinase activity of this mutant. The C-terminal substrate peptide and the nucleotide cofactor were also docked into the active site, and their implications are discussed.

Inclusion of ionization states of ligands in affinity calculations

When estimating binding affinities of a ligand, which can exists in multiple forms, for a target molecule, one must consider all possible competing equilibria. Here, a method is presented that estimates the contribution of the protonation equilibria of a ligand in solution to the measured or calculated binding affinity. The method yields a correction to binding constants that are based on the total concentration of inhibitor (the sum of all ionized forms of the inhibitor in solution) to account for the complexed form of the inhibitor only. The method is applied to the calculation of the difference in binding affinity of two inhibitors, 2-phosphoglycolate (PGA) and its phoshonate analog 3-phosphonopropionate (3PP), for the glycolytic enzyme triosephosphate isomerase. Both inhibitors have three titrating sites and exist in solution as a mixture of different forms. In this case the form that actually binds to the enzyme is present at relative low concentrations. The contributions of the alternative forms to the difference in binding energies is estimated by means of molecular dynamics simulations and corrections. The inhibitors undergo a pK(a) shift upon binding that is estimated by ab initio calculations. An interesting finding is that the affinity difference of the two inhibitors is not due to different interactions in the active site of the enzyme, but rather due to the difference in the solvation properties of the inhibitors.

The empirical valence bond as an effective strategy for computer-aided enzyme design

The ability of the empirical valence bond (EVB) to be used in screening active site residues in enzyme design is explored in a preliminary way. This validation is done by comparing the ability of this approach to evaluate the catalytic contributions of various residues in chorismate mutase. It is demonstrated that the EVB model can serve as an accurate tool in the final stages of computer-aided enzyme design (CAED). The ability of the model to predict quantitatively the catalytic power of enzymes should augment the capacity of current approaches for enzyme design.

Toward accurate screening in computer-aided enzyme design

The ability to design effective enzymes is one of the most fundamental challenges in biotechnology and in some respects in biochemistry. In fact, such ability would be one of the most convincing manifestations of a full understanding of the origin of enzyme catalysis. In this work, we explore the reliability of different simulation approaches, in terms of their ability to rank different possible active site constructs. This validation is done by comparing the ability of different approaches to evaluate the catalytic contributions of various residues in chorismate mutase. It is demonstrated that the empirical valence bond (EVB) model can serve as a practical yet accurate tool in the final stages of computer-aided enzyme design (CAED). Other approaches for fast screening are also examined and found to be less accurate and mainly useful for qualitative screening of ionized residues. It is pointed out that accurate ranking of different options for enzyme design cannot be accomplished by approaches that cannot capture the electrostatic preorganization effect. This is in particular true with regard to current design approaches that use gas phase or small cluster calculations and then estimate the interaction between the enzyme and the transition state (TS) model rather than the TS binding free energy or the relevant activation free energy. The ability of the EVB model to provide a tool for quantitative ranking in the final stage of CAED may help in progressing toward the design of enzymes whose catalytic power is closer to that of native enzymes than to that of the current generation of designer enzymes.

Theoretical study of specific hydrogen-bonding effects on the bridging P-OR bond strength of phosphate monoester dianions

It has been proposed that the driving force for the initial phosphoryl transfer step of protein tyrosine phosphatases (PTPases) could be activation of the substrate ROPO32- by means of an enforced hydrogen-bonding interaction between an aspartic general acid and the bridging oxygen atom O (Zhang et al. Biochemistry 1995, 34, 16088-16096). The potential catalytic effect of this type of interaction, with regard to P-OR bond cleavage, was investigated computationally through simple model systems in which an efficient intramolecular hydrogen bond can take place between a H-bond donor group and the bridging oxygen atom of the dianionic phosphate. The dielectric effect of the environment (epsilon = 1, 4, and 78) was also explored. The results indicate that this interaction causes significant lengthenings of the scissile P-OR bond in all media but with more extreme effects observed in the low dielectric fields epsilon = 1 and epsilon = 4. It is interesting that, in all cases examined, this interaction actually contributes to stabilize the reactant state while causing its P-OR bond to lengthen. Overall, our results support the idea that this specific hydrogen-bonding situation might well be used by PTPases as an important driving force for promoting phosphoryl transfer reactions through highly dissociative transition states.

The planar conformation of a strained proline ring: a QM/MM study

QM and QM/MM energy calculations have been carried out on an atomic resolution structure of liganded triosephosphate isomerase (TIM) that has an active site proline (Pro168) in a planar conformation. The origin of the planarity of this proline has been identified. Steric interactions between the atoms of the proline ring and a tyrosine ring (Tyr166) on one side of the proline prevent the ring from adopting the up pucker (chi1 is approximately -30 degrees), while the side chain of a nearby alanine (Ala171) forbids the down pucker (chi1 is approximately +30 degrees). To obtain a proline conformation that is in agreement with the experimentally observed planar state, a quantum system of sufficient size is required and should at least include the nearby side chains of Tyr166, Ala171, and Glu129 to provide enough stabilization. It is argued that the current force fields for structure optimization do not describe strained protein fragments correctly. The proline is part of a catalytic loop that closes upon ligand binding. Comparison of the proline conformation in different TIM X-ray structures, indicates that in the closed conformation of TIM the proline is planar or nearly planar, while in the open conformation it is down puckered. This suggests that the planarity possibly plays a role in the overall catalytic cycle of TIM, presumable acting as a reservoir of energy that becomes available upon loop opening.

Conformational changes in protein loops and helices induced by post-translational phosphorylation

Post-translational phosphorylation is a ubiquitous mechanism for modulating protein activity and protein-protein interactions. In this work, we examine how phosphorylation can modulate the conformation of a protein by changing the energy landscape. We present a molecular mechanics method in which we phosphorylate proteins in silico and then predict how the conformation of the protein will change in response to phosphorylation. We apply this method to a test set comprised of proteins with both phosphorylated and non-phosphorylated crystal structures, and demonstrate that it is possible to predict localized phosphorylation-induced conformational changes, or the absence of conformational changes, with near-atomic accuracy in most cases. Examples of proteins used for testing our methods include kinases and prokaryotic response regulators. Through a detailed case study of cyclin-dependent kinase 2, we also illustrate how the computational methods can be used to provide new understanding of how phosphorylation drives conformational change, why substituting Glu or Asp for a phosphorylated amino acid does not always mimic the effects of phosphorylation, and how a phosphatase can “capture” a phosphorylated amino acid. This work illustrates how computational methods can be used to elucidate principles and mechanisms of post-translational phosphorylation, which can ultimately help to bridge the gap between the number of known sites of phosphorylation and the number of structures of phosphorylated proteins.

Many proteins are chemically modified after they are synthesized in the cell. These post-translational modifications can modulate the ability of a protein to perform chemical reactions and to interact with other proteins. At the cellular level, for example, these chemical modifications are critical for allowing the cell to respond to its environment and control its division. One of the most common mechanisms by which proteins can be modified is by phosphorylation—the addition of a phosphate group to an amino acid side chain of the protein. Thousands of proteins are known to be modified by phosphorylation, but only for a small minority of these do we have any detailed understanding of how the chemical modification regulates the function of the protein. The authors describe a computational method that can make testable predictions about the structural changes that occur in a protein induced by post-translational phosphorylation. Their results show that the method can produce structural models of the phosphorylated proteins with near-atomic accuracy, and provide insight into the energetics of conformational switches driven by phosphorylation. As such, the computational method complements experiments aimed at understanding the mechanisms of protein regulation by phosphorylation.

Incorporating the effect of ionic strength in free energy calculations using explicit ions

The incorporation of explicit ions to mimic the effect of ionic strength or to neutralize the overall charge on a system in free energy calculations using molecular dynamics simulations is investigated. The difference in the free energy of hydration between two triosephosphate isomerase inhibitors calculated at five different ion concentrations is used as an example. We show that the free energy difference can be highly sensitive to the presence of explicit ions even in cases where the mutation itself does not involve a change in the overall charge. The effect is most significant if the molecule carries a net charge close to the site mutated. Furthermore, it is shown that the introduction of a small number of ions can lead to very severe sampling problems suggesting that in practical calculations convergence can best be achieved by incorporating either no counterions or by simulating at high ionic strength to ensure sufficient sampling of the ion distribution.

14-3-3 Protein C-terminal Stretch Occupies Ligand Binding Groove and Is Displaced by Phosphopeptide Binding

14-3-3 proteins are important regulators of numerous cellular signaling circuits. They bind to phosphorylated protein ligands and regulate their functions by a number of different mechanisms. The C-terminal part of the 14-3-3 protein is known to be involved in the regulation of 14-3-3 binding properties. The structure of this region is unknown; however, a possible location of the C-terminal stretch within the ligand binding groove of the 14-3-3 protein has been suggested. To fully understand the role of the C-terminal stretch in the regulation of the 14-3-3 protein binding properties, we investigated the physical location of the C-terminal stretch and its changes upon the ligand binding. For this purpose, we have used Forster resonance energy transfer (FRET) measurements and molecular dynamics simulation. FRET measurements between Trp242 located at the end of the C-terminal stretch and a dansyl group attached at two different cysteine residues (Cys25 or Cys189) indicated that in the absence of the ligand, the C-terminal stretch occupies the ligand binding groove of 14-3-3 protein. Our data also showed that phosphopeptide binding displaces the C-terminal stretch from the ligand binding groove. Intramolecular distances calculated from FRET measurements fit well with distances obtained from molecular dynamics simulation of full-length 14-3-3zeta protein.

How enzymes work: analysis by modern rate theory and computer simulations

Advances in transition state theory and computer simulations are providing new insights into the sources of enzyme catalysis. Both lowering of the activation free energy and changes in the generalized transmission coefficient (recrossing of the transition state, tunneling, and nonequilibrium contributions) can play a role. A framework for understanding these effects is presented, and the contributions of the different factors, as illustrated by specific enzymes, are identified and quantified by computer simulations. The resulting understanding of enzyme catalysis is used to comment on alternative proposals of how enzymes work.

14-3-3ζ C-terminal Stretch Changes Its Conformation upon Ligand Binding and Phosphorylation at Thr232

14-3-3 proteins are abundant binding proteins involved in many biologically important processes. 14-3-3 proteins bind to other proteins in a phosphorylation-dependent manner and function as scaffold molecules modulating the activity of their binding partners. In this work, we studied the conformational changes of 14-3-3 C-terminal stretch, a region implicated in playing a role in the regulation of 14-3-3. Time-resolved fluorescence and molecular dynamics were used to investigate structural changes of the C-terminal stretch induced by phosphopeptide binding and phosphorylation at Thr232, a casein kinase I phosphorylation site located within this region. A tryptophan residue placed at position 242 was exploited as an intrinsic fluorescence probe of the C-terminal stretch dynamics. Other tryptophan residues were mutated to phenylalanine. Time-resolved fluorescence measurements revealed that phosphopeptide binding changes the conformation and increases the flexibility of 14-3-3zeta C-terminal stretch, demonstrating that this region is directly involved in ligand binding. Phosphorylation of 14-3-3zeta at Thr232 resulted in inhibition of phosphopeptide binding and suppression of 14-3-3-mediated enhancement of serotonin N-acetyltransferase activity. Time-resolved fluorescence of Trp242 also revealed that phosphorylation at Thr232 induces significant changes of the C-terminal stretch conformation. In addition, molecular dynamic simulations suggest that phosphorylation at Thr232 induces a more extended conformation of 14-3-3zeta C-terminal stretch and changes its interaction with the rest of the 14-3-3 molecule. These results indicate that the conformation of the C-terminal stretch plays an important role in the regulation of 14-3-3 binding properties.

Hydrolysis of cyclic phosphates by ribonuclease A: a computational study using a simplified ab initio quantum model

The second step in the enzyme-catalyzed hydrolysis of phosphate esters by ribonuclease A (RNase A) was studied using an ab initio quantum-based model of the active site including constrained parts of three critical residues, His-12, His-119, and Lys-41, and a small substrate. The competition between release of the cyclic phosphate intermediate and subsequent hydrolysis following transphosphorylation was explored to determine the electronic factors that contribute to preferential intermediate product release observed experimentally. The structural and energetic results obtained at both the RHF and MP2 levels reveal several contributing factors consistent with experimental observation. Although the intrinsic electronic effects tend to favor hydrolysis slightly with an overall activation free energy of approximately 70 kJ mol(-1), entropic and environmental effects favor release of the cyclic phosphate intermediate over hydrolysis. Exploration of the second, hydrolysis step also revealed interesting similarity with the transphosphorylation step, including the observation of autocatalysis by the substrate. Moreover, both steps of the overall RNase A reaction reveal multiple pathways involving proton transfers to sites of similar proton affinities. The anionic phosphate in both steps can act as a stable proton binding site as protons are moved around the active site throughout the progress of the reaction. These results suggest autocatalysis may be representative of more general behavior in enzymes containing highly charged substrates, especially phosphates.

Computer simulations of enzyme catalysis: methods, progress, and insights

Understanding the action of enzymes on an atomistic level is one of the important aims of modern biophysics. This review describes the state of the art in addressing this challenge by simulating enzymatic reactions. It considers different modeling methods including the empirical valence bond (EVB) and more standard molecular orbital quantum mechanics/molecular mechanics (QM/MM) methods. The importance of proper configurational averaging of QM/MM energies is emphasized, pointing out that at present such averages are performed most effectively by the EVB method. It is clarified that all properly conducted simulation studies have identified electrostatic preorganization effects as the source of enzyme catalysis. It is argued that the ability to simulate enzymatic reactions also provides the chance to examine the importance of nonelectrostatic contributions and the validity of the corresponding proposals. In fact, simulation studies have indicated that prominent proposals such as desolvation, steric strain, near attack conformation, entropy traps, and coherent dynamics do not account for a major part of the catalytic power of enzymes. Finally, it is pointed out that although some of the issues are likely to remain controversial for some time, computer modeling approaches can provide a powerful tool for understanding enzyme catalysis.

Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase

Understanding the chemical step in the catalytic reaction of DNA polymerases is essential for elucidating the molecular basis of the fidelity of DNA replication. The present work evaluates the free energy surface for the nucleotide transfer reaction of T7 polymerase by free energy perturbation/empirical valence bond (FEP/EVB) calculations. A key aspect of the enzyme simulation is a comparison of enzymatic free energy profiles with the corresponding reference reactions in water using the same computational methodology, thereby enabling a quantitative estimate for the free energy of the nucleotide insertion reaction. The reaction is driven by the FEP/EVB methodology between valence bond structures representing the reactant, pentacovalent intermediate, and the product states. This pathway corresponds to three microscopic chemical steps, deprotonation of the attacking group, a nucleophilic attack on the P(alpha) atom of the dNTP substrate, and departure of the leaving group. Three different mechanisms for the first microscopic step, the generation of the RO(-) nucleophile from the 3'-OH hydroxyl of the primer, are examined: (i) proton transfer to the bulk solvent, (ii) proton transfer to one of the ionic oxygens of the P(alpha) phosphate group, and (iii) proton transfer to the ionized Asp654 residue. The most favorable reaction mechanism in T7 pol is predicted to involve the proton transfer to Asp654. This finding sheds light on the long standing issue of the actual role of conserved aspartates. The structural preorganization that helps to catalyze the reaction is also considered and analyzed. The overall calculated mechanism consists of three subsequent steps with a similar activation free energy of about 12 kcal/mol. The similarity of the activation barriers of the three microscopic chemical steps indicates that the T7 polymerase may select against the incorrect dNTP substrate by raising any of these barriers. The relative height of these barriers comparing right and wrong dNTP substrates should therefore be a primary focus of future computational studies of the fidelity of DNA polymerases.

Density functional study of the mechanism of a tyrosine phosphatase: I. Intermediate formation

The first step in the catalytic mechanism of a protein tyrosine phosphatase, the transfer of a phosphate group from the phosphotyrosine substrate to a cysteine side chain of the protein to form a phosphoenzyme intermediate, has been studied by combining density functional calculations of an active-site cluster with continuum electrostatic descriptions of the solvent and the remainder of the protein. This approach provides the high level of quantum chemical methodology needed to adequately model phosphotransfer reactions with a reasonable description of the environment around the active site. Energy barriers and geometries along a reaction pathway are calculated. In the literature, mechanisms assuming both a monoanionic and a dianionic substrate have been proposed; this disagreement is addressed by performing calculations for both possibilities. For the dianionic substrate, a dissociative reaction pathway with early proton transfer to the leaving group and a 9 kcal/mol energy barrier is predicted (the experimental estimate is ca. 14 kcal/mol), while for the monoanionic substrate, an associative pathway with late proton transfer and a 22 kcal/mol energy barrier is predicted. These results, together with a review of experimental evidence, support the dianionic-substrate/dissociative-pathway alternative. The relationship between a dianionic or monoanionic substrate and a dissociative or associative pathway, respectively, can be understood in terms of classical organic chemical reaction pathways.

Quantum mechanical methods for enzyme kinetics

This review discusses methods for the incorporation of quantum mechanical effects into enzyme kinetics simulations in which the enzyme is an explicit part of the model. We emphasize three aspects: (a) use of quantum mechanical electronic structure methods such as molecular orbital theory and density functional theory, usually in conjunction with molecular mechanics; (b) treating vibrational motions quantum mechanically, either in an instantaneous harmonic approximation, or by path integrals, or by a three-dimensional wave function coupled to classical nuclear motion; (c) incorporation of multidimensional tunneling approximations into reaction rate calculations.

A novel dual specificity phosphatase SKRP1 interacts with the MAPK kinase MKK7 and inactivates the JNK MAPK pathway. Implication for the precise regulation of the particular MAPK pathway

Mitogen-activated protein kinases (MAPKs) are activated in response to various extracellular stimuli, and their activities are regulated by upstream activating kinases and protein phosphatases such as MAPK phosphatases (MKPs). We report the identification and characterization of a novel MKP termed SKRP1 (SAPK pathway-regulating phosphatase 1). It contains an extended active site sequence motif conserved in all MKPs but lacks a Cdc25 homology domain. Immunoblotting analysis revealed that SKRP1 is constitutively expressed, and its transcripts of 4.0 and 1.0 kb were detected in almost tissues examined. SKRP1 was highly specific for c-Jun N-terminal kinase (JNK) in vitro and effectively suppressed the JNK activation in response to tumor necrosis factor alpha or thapsigargin. Endogenous SKRP1 was present predominantly in the cytoplasm and co-localized with JNK. However, SKRP1 does not bind directly to its target JNK, but co-precipitation of SKRP1 with the MAPK kinase MKK7, a JNK activator, was found in vitro and in vivo. Furthermore, we found that SKRP1 did not interfere with the co-precipitation of MKK7 with JNK. Together, our findings indicate that SKRP1 interacts with its physiological substrate JNK through MKK7, thereby leading to the precise regulation of JNK activity in vivo.

The catalytic mechanism of protein tyrosine phosphatases revisited

Experimental and theoretical studies of the catalytic mechanism in protein tyrosine phosphatases and dual specific phosphatases are reviewed. The structural properties of these enzymes contributing to the efficient rate enhancement of phosphate monoester hydrolysis have been established during the last decade. There are, however, uncertainties in the interpretation of available experimental data that make the commonly assumed reaction mechanism somewhat doubtful. Theoretical calculations as well as analysis of crystal structures point towards an alternative interpretation of the ionisation state in the reactive complex.

*NO, RSNO, ONOO-, NO+, *NOO, NOx--dynamic regulation of oxidant scavenging, nitric oxide stores, and cyclic GMP-independent cell signaling

Following its release from nitric oxide synthase, nitric oxide seldom perfuses the cytosol; rather this reactive mediator quickly interacts with available target molecules proximate to its site of release. Within the cell, virtually every component, low-molecular-weight oxidants and reductants, proteins, lipids, sugars, and nucleic acids can be modified by nitrogen oxides thus acting as potential targets for reactive nitrogen oxides. Adducts formed by nitrogen oxides often modulate the cellular activities of the target molecules, and these modified molecules may be differentially metabolized or localized. The formation of nitrogen oxide adducts can be a reversible process, and the reactive nitrogen species released may be specifically oxidized or reduced during the process. Recently, numerous studies have demonstrated that reversible nitration of cellular proteins acts to transduce molecular signals regulating such diverse processes as muscle contraction, neurotransmission, protein metabolism, and apoptosis. The vast numbers of molecules that undergo biologically relevant interactions with nitrogen oxides imply that the cellular concentration of nitrosated and nitrated species may effectively comprise a reserve or cellular store. Potentially, these nitroso reserves function as critical components of the overall redox status of the intracellular environs. Understanding the dynamic regulation of nitric oxide/nitrogen oxides release from these stores is likely to provide clues important in resolving the complex pathophysiology of poorly understood multifactorial disorders, including neurodegeneration, multiorgan failure, cardiomyopathy, and septic shock.