Interaction of human 3-phosphoglycerate kinase with its two substrates: is substrate antagonism a kinetic advantage?

Substrate antagonism has been described for a variety of enzymes with more than one substrate and is characterized by a lowering of the affinity of one substrate in the presence of the other(s). 3-Phosphoglycerate kinase (PGK) catalyzes phosphotransfer from 1,3-bisphosphoglycerate (bPG) to ADP to give 3-phosphoglycerate (PG) and ATP, and is subject to substrate antagonism. Because of the instability of bPG, antagonism has only been described between PG and ATP or ADP. Here, we show that antagonism also occurs between bPG and ADP. Using the stopped-flow method, we show that the dissociation constant for one substrate increases in the presence of the other, and that this decrease in affinity is mainly due to an increase in the dissociation rate constant. As a consequence, there is an increase in the overall interaction kinetics. Interestingly, in the presence of the mirror image of natural d-ADP, l-ADP (a good substrate for PGK), antagonism is absent. Using rapid-quench-flow, we studied the kinetics of ATP formation. The time courses present the following: (1) a lag with l-ADP, but not with d-ADP, the kinetics of which were similar to the interaction kinetics measured by stopped-flow; (2) a burst that is directed by the phosphotransfer; and (3) a steady-state that is rate limited by the release of product kinetics. Structural explanations for these results are proposed by analyzing the crystallographic structure of the fully closed conformation of PGK in complex with l-ADP, PG, and the transition-state analogue AlF₄⁻ compared to previously determined structures.     

Differences in the transient kinetics of the binding of D-ADP and its mirror image L-ADP to human 3-phosphoglycerate kinase revealed by the presence of 3-phosphoglycerate

L-Nucleosides comprise a new class of antiviral and anticancer agents that are converted in vivo by a cascade of kinases to pharmacologically active nucleoside triphosphates. The last step of the cascade may be catalyzed by 3-phosphoglycerate kinase (PGK), an enzyme that has low specificity for nucleoside diphosphate (NDP): NDP + 1,3-bisphosphoglycerate <--> NTP + 3-phosphoglycerate. Here we compared the kinetics of the formation of the complexes of human PGK with d- and its mirror image l-ADP and the effect of 3-phosphoglycerate (PG) on these by exploiting the fluorescence signal of PGK that occurs upon its interaction with nucleotide substrate. Two types of experiment were carried out: equilibrium (estimation of dissociation constants) and stopped-flow (transient kinetics of the interactions). We show that under our experimental conditions (buffer containing 30% methanol, 4 degrees C) PGK binds d- and l-ADP with similar kinetics. However, whereas PG increased the dissociation rate constant for d-ADP by a factor of 8-which is a kinetic explanation for "substrate antagonism"-PG had little effect on this constant for l-ADP. We explain this difference by a molecular modeling study that showed that the beta-phosphates of d- and l-ADP have different orientations when bound to the active site of human PGK. The difference is unexpected because l-ADP is almost as catalytically competent as d-ADP [ Varga, A. et al. (2008) Biochem. Biophys. Res. Commun. 366, 994-1000].     

The mechanism of inhibition of the cyclin-dependent kinase-2 as revealed by the molecular dynamics study on the complex CDK2 with the peptide substrate HHASPRK

Molecular dynamics (MD) simulations were used to explain structural details of cyclin-dependent kinase-2 (CDK2) inhibition by phosphorylation at T14 and/or Y15 located in the glycine-rich loop (G-loop). Ten-nanosecond-long simulations of fully active CDK2 in a complex with a short peptide (HHASPRK) substrate and of CDK2 inhibited by phosphorylation of T14 and/or Y15 were produced. The inhibitory phosphorylations at T14 and/or Y15 show namely an ATP misalignment and a G-loop shift (~5 A) causing the opening of the substrate binding box. The biological functions of the G-loop and GxGxxG motif evolutionary conservation in protein kinases are discussed. The position of the ATP gamma-phosphate relative to the phosphorylation site (S/T) of the peptide substrate in the active CDK2 is described and compared with inhibited forms of CDK2. The MD results clearly provide an explanation previously not known as to why a basic residue (R/K) is preferred at the P(2) position in phosphorylated S/T peptide substrates.     

A comparative study of HIV-1 and HTLV-I protease structure and dynamics reveals a conserved residue interaction network

The two retroviruses human T-lymphotropic virus type I (HTLV-I) and human immunodeficiency virus type 1 (HIV-1) are the causative agents of severe and fatal diseases including adult T-cell leukemia and the acquired immune deficiency syndrome (AIDS). Both viruses code for a protease that is essential for replication and therefore represents a key target for drugs interfering with viral infection. The retroviral proteases from HIV-1 and HTLV-I share 31% sequence identity and high structural similarities. Yet, their substrate specificities and inhibition profiles differ substantially. In this study, we performed all-atom molecular dynamics (MD) simulations for both enzymes in their ligand-free states and in complex with model substrates in order to compare their dynamic behaviors and enhance our understanding of the correlation between sequence, structure, and dynamics in this protein family. We found extensive similarities in both local and overall protein dynamics, as well as in the energetics of their interactions with model substrates. Interestingly, those residues that are important for strong ligand binding are frequently not conserved in sequence, thereby offering an explanation for the differences in binding specificity. Moreover, we identified an interaction network of contacts between conserved residues that interconnects secondary structure elements and serves as a scaffold for the protein fold. This interaction network is conformationally stable over time and may provide an explanation for the highly similar dynamic behavior of the two retroviral proteases, even in the light of their rather low overall sequence identity.     

Molecular mechanism of substrate specificity in the bacterial neutral amino acid transporter LeuT

The recently published X-ray structure of LeuT, a Na(+)/Cl(-)-dependent neurotransmitter transporter, has provided fresh impetus to efforts directed at understanding the molecular principles governing specific neurotransmitter transport. The combination of the LeuT crystal structure with the results of molecular simulations enables the functional data on specific binding and transport to be related to molecular structure. All-atom FEP and molecular dynamics (MD) simulations of LeuT embedded in an explicit membrane were performed alongside a decomposition analysis to dissect the molecular determinants of the substrate specificity of LeuT. It was found that the ligand must be in a zwitterionic (ZW) form to bind tightly to the transporter. The theoretical results on the absolute binding-free energies for leucine, alanine, and glycine show that alanine can be a potent substrate for LeuT, although leucine is preferred, which is consistent with the recent experimental data (Singh et al., Nature 2007;448:952-956). Furthermore, LeuT displays robust specificity for leucine over glycine. Interestingly, the ability of LeuT to discriminate between substrates relies on the dynamics of residues that form its binding pocket (e.g., F253 and Q250) and the charged side chains (R30-D404) from a second coordination shell. The water-mediated R30-D404 salt bridge is thought to be part of the extracellular (EC) gate of LeuT. The introduction of a polar ligand such as glycine to the water-depleted binding pocket of LeuT gives rise to structural rearrangements of the R30-D404-Q250 hydrogen-bonding network and leads to increased hydration of the binding pocket. Conformational changes associated with the broken hydrogen bond between Q250 and R30 are shown to be important for tight and selective ligand binding to LeuT.     

Coupling of Calcium and Substrate Binding through Loop Alignment in the Outer-Membrane Transporter BtuB

In Gram-negative bacteria, TonB-dependent outer-membrane transporters bind large, scarce organometallic substrates with high affinity preceding active transport. The cobalamin transporter BtuB requires the additional binding of two Ca²⁺ ions before substrate binding can occur, but the underlying molecular mechanism is unknown. Using the crystallographic structures available for different bound states of BtuB, we have carried out extended molecular dynamics simulations of multiple functional states of BtuB to address the role of Ca²⁺ in substrate recruitment. We find that Ca²⁺ binding both stabilizes and repositions key extracellular loops of BtuB, optimizing interactions with the substrate. Interestingly, replacement by Mg²⁺ abolishes this effect, in accordance with experiments. Using a set of new force-field parameters developed for cyanocobalamin, we also simulated the substrate-bound form of BtuB, where we observed interactions not seen in the crystal structure between the substrate and loops previously found to be important for binding and transport. Based on our results, we suggest that the large size of cobalamin compared to other TonB-dependent transporter substrates explains the requirement of Ca²⁺ binding for high-affinity substrate recruitment in BtuB.     

An Allosteric Signaling Pathway of Human 3-Phosphoglycerate Kinase from Force Distribution Analysis

3-Phosphogycerate kinase (PGK) is a two domain enzyme, which transfers a phosphate group between its two substrates, 1,3-bisphosphoglycerate bound to the N-domain and ADP bound to the C-domain. Indispensable for the phosphoryl transfer reaction is a large conformational change from an inactive open to an active closed conformation via a hinge motion that should bring substrates into close proximity. The allosteric pathway resulting in the active closed conformation has only been partially uncovered. Using Molecular Dynamics simulations combined with Force Distribution Analysis (FDA), we describe an allosteric pathway, which connects the substrate binding sites to the interdomain hinge region. Glu192 of alpha-helix 7 and Gly394 of loop L14 act as hinge points, at which these two secondary structure elements straighten, thereby moving the substrate-binding domains towards each other. The long-range allosteric pathway regulating hPGK catalytic activity, which is partially validated and can be further tested by mutagenesis, highlights the virtue of monitoring internal forces to reveal signal propagation, even if only minor conformational distortions, such as helix bending, initiate the large functional rearrangement of the macromolecule.     

Molecular dynamics simulation of the substitution of serine for the conserved glycine in the G-loop in the cdc28-srm yeast mutant using the crystal lattice of human CDK2 kinase

Two-nanosecond molecular dynamics simulations of the crystal lattice of an active complex of pT160-CDK2 kinase/cyclin A/ATP-Mg²⁺/substrate were performed. The simulations showed that the structures of the wild-type CDK2 complex and the mutant CDK2 complex involving the substitution G16S-CDK2 corresponding to the yeast substitution G20S-CDC28 differ noticeably and the differences between the structural conformations are most pronounced in the regions that play a key role in the kinase functioning. The results of the computer calculations were used to consider the structural elements that may affect the kinase activity, the regulatory phosphorylation, and the binding of protein kinase with cyclins and substrates.     

Long‐range molecular dynamics show that inactive forms of Protein Kinase A are more dynamic than active forms

Many protein kinases are characterized by at least two structural forms corresponding to the highest level of activity (active) and low or no activity, (inactive). Further, protein dynamics is an important consideration in understanding the molecular and mechanistic basis of enzyme function. In this work, we use protein kinase A (PKA) as the model system and perform microsecond range molecular dynamics (MD) simulations on six variants which differ from one another in terms of active and inactive form, with or without bound ligands, C‐terminal tail and phosphorylation at the activation loop. We find that the root mean square fluctuations in the MD simulations are generally higher for the inactive forms than the active forms. This difference is statistically significant. The higher dynamics of inactive states has significant contributions from ATP binding loop, catalytic loop, and αG helix. Simulations with and without C‐terminal tail show this differential dynamics as well, with lower dynamics both in the active and inactive forms if C‐terminal tail is present. Similarly, the dynamics associated with the inactive form is higher irrespective of the phosphorylation status of Thr 197. A relatively stable stature of active kinases may be better suited for binding of substrates and detachment of the product. Also, phosphoryl group transfer from ATP to the phosphosite on the substrate requires precise transient coordination of chemical entities from three different molecules, which may be facilitated by the higher stability of the active state.     

基于弹性网络模型的蛋白质变构路径与关键残基识别研究

目的 变构效应在蛋白质生物学功能执行过程中发挥着重要的调控作用，如何基于蛋白质空间结构，有效识别变构信号的传播路径和关键的残基位点是蛋白质结构-功能关系研究领域的热点科学问题。方法 本研究利用基于弹性网络模型（elastic network model，ENM）的力分布计算方法，通过分析蛋白质对外力的响应过程，来识别体系的变构路径以及变构过程中的关键残基。在该方法中，对蛋白质的关键变构位点施加外力，通过对体系形变以及内力分布情况的分析，有效识别与外力承载区域形变相耦合的关键残基，从而得到力信号在蛋白质结构内的传播路径。结果 利用该方法研究了人类磷酸甘油酸激酶（human phosphoglycerate kinase，hPGK）和蛋白质酪氨酸磷酸酶（protein tyrosine phosphatase，PTP）PDZ2结构域的变构调控路径和关键残基。对于hPGK，识别出从底物结合位点到铰链区的两条变构信号传导路径。对于PTP PDZ2，也成功识别出从配体结合位点传递到蛋白质远端的两条长程变构调控路径。计算结果与实验和分子动力学（molecular dynamics，MD）模拟得到的结果一致。结论 本研究为蛋白质体系关键残基识别及变构路径研究提供了有效的分析方法。     

Substrate Binding Specifically Modulates Domain Arrangements in Adenylate Kinase

The enzyme adenylate kinase (ADK) features two substrate binding domains that undergo large-scale motions during catalysis. In the apo state, the enzyme preferentially adopts a globally open state with accessible binding sites. Binding of two substrate molecules (AMP + ATP or ADP + ADP) results in a closed domain conformation, allowing efficient phosphoryl-transfer catalysis. We employed molecular dynamics simulations to systematically investigate how the individual domain motions are modulated by the binding of substrates. Two-dimensional free-energy landscapes were calculated along the opening of the two flexible lid domains for apo and holo ADK as well as for all single natural substrates bound to one of the two binding sites of ADK. The simulations reveal a strong dependence of the conformational ensembles on type and binding position of the bound substrates and a nonsymmetric behavior of the lid domains. Altogether, the ensembles suggest that, upon initial substrate binding to the corresponding lid site, the opposing lid is maintained open and accessible for subsequent substrate binding. In contrast, ATP binding to the AMP-lid induces global domain closing, preventing further substrate binding to the ATP-lid site. This might constitute a mechanism by which the enzyme avoids the formation of a stable but enzymatically unproductive state.     

Perturbation of yeast 3-phosphoglycerate kinase reaction mixtures with ADP: transient kinetics of formation of ATP from bound 1,3-bisphosphoglycerate

3-Phosphoglycerate kinase (PGK) is the first ATP-producing enzyme in glycolysis: ADP + 1,3-bisphosphoglycerate (bPG) <--> ATP + 3-phosphoglycerate (PG). Whereas extensive studies have been carried out on its structure, there is less information about its reaction pathway, which is usually studied in the reverse direction because of the instability of bPG. We studied the transients of the PGK reaction by chemical sampling in a rapid quench flow apparatus, using [gamma-(32)P]ATP, in 30% methanol at 4 degrees C to decrease k(cat). There were two types of experiment, both at low PG concentrations to prevent bPG release. In the first, reaction mixtures were quenched in acid at different times (from 4 ms) and the bPG concentrations were determined. This type gave information about the ATP binding and phospho-transfer steps. In the second, PGK reaction mixtures at equilibrium were perturbed by the injection of ADP, the new mixtures aged for different times and quenched in acid, and the bPG concentrations were determined. This gave information about the kinetics of the binding of ADP to a PGK intermediate. The data from the two types of experiments were fitted to simple schemes and then treated together by a global fitting procedure using a five-step pathway, deduced from previous structural studies. Under our conditions, it appears that (1) a binary PGK.bPG complex is an important intermediate on the reaction pathway, i.e., that ADP is released before bPG, (2) ADP binds to a "closed" conformation in the PGK.bPG complex, and (3) the PGK reaction can be studied in the physiologically important direction without having to handle bPG.     

The Molecular Mechanism of Ion-Dependent Gating in Secondary Transporters

LeuT-like fold Na-dependent secondary active transporters form a large family of integral membrane proteins that transport various substrates against their concentration gradient across lipid membranes, using the free energy stored in the downhill concentration gradient of sodium ions. These transporters play an active role in synaptic transmission, the delivery of key nutrients, and the maintenance of osmotic pressure inside the cell. It is generally believed that binding of an ion and/or a substrate drives the conformational dynamics of the transporter. However, the exact mechanism for converting ion binding into useful work has yet to be established. Using a multi-dimensional path sampling (string-method) followed by all-atom free energy simulations, we established the principal thermodynamic and kinetic components governing the ion-dependent conformational dynamics of a LeuT-like fold transporter, the sodium/benzyl-hydantoin symporter Mhp1, for an entire conformational cycle. We found that inward-facing and outward-facing states of Mhp1 display nearly the same free energies with an ion absent from the Na2 site conserved across the LeuT-like fold transporters. The barrier separating an apo-state from inward-facing or outward-facing states of the transporter is very low, suggesting stochastic gating in the absence of ion/substrate bound. In contrast, the binding of a Na2 ion shifts the free energy stabilizing the outward-facing state and promoting substrate binding. Our results indicate that ion binding to the Na2 site may also play a key role in the intracellular thin gate dynamics modulation by altering its interactions with the transmembrane helix 5 (TM5). The Potential of Mean Force (PMF) computations for a substrate entrance displays two energy minima that correspond to the locations of the main binding site S1 and proposed allosteric S2 binding site. However, it was found that substrate''s binds to the site S1 ∼5 kcal/mol more favorable than that to the site S2 for all studied bound combinations of ions and a substrate.     

Small molecules containing hetero-bicyclic ring systems compete with UDP-Glc for binding to WaaG glycosyltransferase

The α-1,3-glucosyltransferase WaaG is involved in the synthesis of the core region of lipopolysaccharides in E. coli. A fragment-based screening for inhibitors of the WaaG glycosyltrasferase donor site has been performed using NMR spectroscopy. Docking simulations were performed for three of the compounds of the fragment library that had shown binding activity towards WaaG and yielded 3D models for the respective complexes. The three ligands share a hetero-bicyclic ring system as a common structural motif and they compete with UDP-Glc for binding. Interestingly, one of the compounds promoted binding of uridine to WaaG, as seen from STD NMR titrations, suggesting a different binding mode for this ligand. We propose these compounds as scaffolds for the design of selective high-affinity inhibitors of WaaG. Binding of natural substrates, enzymatic activity and donor substrate selectivity were also investigated by NMR spectroscopy. Molecular dynamics simulations of WaaG were carried out with and without bound UDP and revealed structural changes compared to the crystal structure and also variations in flexibility for some amino acid residues between the two WaaG systems studied.     

Phosphorylation-induced conformational changes in a mitogen-activated protein kinase substrate. Implications for tyrosine hydroxylase activation

Mitogen-activated protein (MAP) kinase-mediated phosphorylation of specific residues in tyrosine hydroxylase leads to an increase in enzyme activity. However, the mechanism whereby phosphorylation affects enzyme turnover is not well understood. We used a combination of fluorescence resonance energy transfer (FRET) measurements and molecular dynamics simulations to explore the conformational free energy landscape of a 10-residue MAP kinase substrate found near the N terminus of the enzyme. This region is believed to be part of an autoregulatory sequence that overlies the active site of the enzyme. FRET was used to measure the effect of phosphorylation on the ensemble of peptide conformations, and molecular dynamics simulations generated free energy profiles for both the unphosphorylated and phosphorylated peptides. We demonstrate how FRET transfer efficiencies can be calculated from molecular dynamics simulations. For both the unphosphorylated and phosphorylated peptides, the calculated FRET efficiencies are in excellent agreement with the experimentally determined values. Moreover, the FRET measurements and molecular simulations suggest that phosphorylation causes the peptide backbone to change direction and fold into a compact structure relative to the unphosphorylated state. These results are consistent with a model of enzyme activation where phosphorylation of the MAP kinase substrate causes the N-terminal region to adopt a compact structure away from the active site. The methods we employ provide a general framework for analyzing the accessible conformational states of peptides and small molecules. Therefore, they are expected to be applicable to a variety of different systems.     

Enantioselectivity of human AMP, dTMP and UMP-CMP kinases

L-nucleoside analogues such as lamivudine are active for treating viral infections. Like D-nucleosides, the biological activity of the L-enantiomers requires their stepwise phosphorylation by cellular or viral kinases to give the triphosphate. The enantioselectivity of NMP kinases has not been thoroughly studied, unlike that of deoxyribonucleoside kinases. We have therefore investigated the capacity of L-enantiomers of some natural (d)NMP to act as substrates for the recombinant forms of human uridylate-cytidylate kinase, thymidylate kinase and adenylate kinases 1 and 2. Both cytosolic and mitochondrial adenylate kinases were strictly enantioselective, as they phosphorylated only D-(d)AMP. L-dTMP was a substrate for thymidylate kinase, but with an efficiency 150-fold less than D-dTMP. Both L-dUMP and L-(d)CMP were phosphorylated by UMP-CMP kinase although much less efficiently than their natural counterparts. The stereopreference was conserved with the 2'-azido derivatives of dUMP and dUMP while, unexpectedly, the 2'-azido-D-dCMP was a 4-fold better substrate for UMP-CMP kinase than was CMP. Docking simulations showed that the small differences in the binding of D-(d)NMP to their respective kinases could account for the differences in interactions of the L-isomers with the enzymes. This in vitro information was then used to develop the in vivo activation pathway for L-dT.     

Molecular basis for the lack of enantioselectivity of human 3-phosphoglycerate kinase

Non-natural L-nucleoside analogues are increasingly used as therapeutic agents to treat cancer and viral infections. To be active, L-nucleosides need to be phosphorylated to their respective triphosphate metabolites. This stepwise phosphorylation relies on human enzymes capable of processing L-nucleoside enantiomers. We used crystallographic analysis to reveal the molecular basis for the low enantioselectivity and the broad specificity of human 3-phosphoglycerate kinase (hPGK), an enzyme responsible for the last step of phosphorylation of many nucleotide derivatives. Based on structures of hPGK in the absence of nucleotides, and bound to L and d forms of MgADP and MgCDP, we show that a non-specific hydrophobic clamp to the nucleotide base, as well as a water-filled cavity behind it, allows high flexibility in the interaction between PGK and the bases. This, combined with the dispensability of hydrogen bonds to the sugar moiety, and ionic interactions with the phosphate groups, results in the positioning of different nucleotides so to expose their diphosphate group in a position competent for catalysis. Since the third phosphorylation step is often rate limiting, our results are expected to alleviate in silico tailoring of L-type prodrugs to assure their efficient metabolic processing.     

Interaction of human 3-phosphoglycerate kinase with L-ADP, the mirror image of D-ADP

l-Nucleoside-analogues, mirror images of the natural d-nucleosides, are a new class of antiviral and anticancer agents. In the cell they have to be phosphorylated to pharmacologically active triphosphate forms, the last step seems to involve human 3-phosphoglycerate kinase (hPGK). Here we present a steady state kinetic and biophysical study of the interaction of the model compound l-MgADP with hPGK. l-MgADP is a good substrate with k_cat and K_m values of 685 s⁻¹ and 0.27 mM, respectively. Double inhibition studies suggest that l-MgADP binds to the specific adenosine-binding site and protects the conformation of hPGK molecule against heat denaturation, as detected by microcalorimetry. Structural details of the interaction in the enzyme active site are different for the d- and l-enantiomers (e.g. the effect of Mg²⁺), but these differences do not prevent the occurrence of the catalytic cycle, which is accompanied by the hinge-bending domain closure, as indicated by SAXS measurements.     

Dynamics of Preferential Substrate Recognition in HIV-1 Protease: Redefining the Substrate Envelope

Human immunodeficiency virus type 1 (HIV-1) protease (PR) permits viral maturation by processing the gag and gag-pro-pol polyproteins. HIV-1 PR inhibitors (PIs) are used in combination antiviral therapy but the emergence of drug resistance has limited their efficacy. The rapid evolution of HIV-1 necessitates consideration of drug resistance in novel drug design. Drug-resistant HIV-1 PR variants no longer inhibited efficiently, continue to hydrolyze the natural viral substrates. Though highly diverse in sequence, the HIV-1 PR substrates bind in a conserved three-dimensional shape we termed the substrate envelope. Earlier, we showed that resistance mutations arise where PIs protrude beyond the substrate envelope, because these regions are crucial for drug binding but not for substrate recognition. We extend this model by considering the role of protein dynamics in the interaction of HIV-1 PR with its substrates. We simulated the molecular dynamics of seven PR-substrate complexes to estimate the conformational flexibility of the bound substrates. Interdependence of substrate-protease interactions might compensate for variations in cleavage-site sequences and explain how a diverse set of sequences are recognized as substrates by the same enzyme. This diversity might be essential for regulating sequential processing of substrates. We define a dynamic substrate envelope as a more accurate representation of PR-substrate interactions. This dynamic substrate envelope, described by a probability distribution function, is a powerful tool for drug design efforts targeting ensembles of resistant HIV-1 PR variants with the aim of developing drugs that are less susceptible to resistance.     

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. II. substrate-enzyme interactions and dynamics

Molecular dynamics simulations of protein tyrosine phosphatase 1B (PTP1B) complexed with the phosphorylated peptide substrate DADEpYL and the free substrate have been conducted to investigate 1) the physical forces involved in substrate-protein interactions, 2) the importance of enzyme and substrate flexibility for binding, 3) the electrostatic properties of the enzyme, and 4) the contribution from solvation. The simulations were performed for 1 ns, using explicit water molecules. The last 700 ps of the trajectories was used for analysis determining enthalpic and entropic contributions to substrate binding. Based on essential dynamics analysis of the PTP1B/DADEpYL trajectory, it is shown that internal motions in the binding pocket occur in a subspace of only a few degrees of freedom. In particular, relatively large flexibilities are observed along several eigenvectors in the segments: Arg(24)-Ser(28), Pro(38)-Arg(47), and Glu(115)-Gly(117). These motions are correlated to the C- and N-terminal motions of the substrate. Relatively small fluctuations are observed in the region of the consensus active site motif (H/V)CX(5)R(S/T) and in the region of the WPD loop, which contains the general acid for catalysis. Analysis of the individual enzyme-substrate interaction energies revealed that mainly electrostatic forces contribute to binding. Indeed, calculation of the electrostatic field of the enzyme reveals that only the field surrounding the binding pocket is positive, while the remaining protein surface is characterized by a predominantly negative electrostatic field. This positive electrostatic field attracts negatively charged substrates and could explain the experimentally observed preference of PTP1B for negatively charged substrates like the DADEpYL peptide.