Alpha7 helix plays an important role in the conformational stability of PTP1B

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix α7), which forms an allosteric binding site 20 ? away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

The Mechanism of Allosteric Inhibition of Protein Tyrosine Phosphatase 1B

As the prototypical member of the PTP family, protein tyrosine phosphatase 1B (PTP1B) is an attractive target for therapeutic interventions in type 2 diabetes. The extremely conserved catalytic site of PTP1B renders the design of selective PTP1B inhibitors intractable. Although discovered allosteric inhibitors containing a benzofuran sulfonamide scaffold offer fascinating opportunities to overcome selectivity issues, the allosteric inhibitory mechanism of PTP1B has remained elusive. Here, molecular dynamics (MD) simulations, coupled with a dynamic weighted community analysis, were performed to unveil the potential allosteric signal propagation pathway from the allosteric site to the catalytic site in PTP1B. This result revealed that the allosteric inhibitor compound-3 induces a conformational rearrangement in helix α7, disrupting the triangular interaction among helix α7, helix α3, and loop11. Helix α7 then produces a force, pulling helix α3 outward, and promotes Ser190 to interact with Tyr176. As a result, the deviation of Tyr176 abrogates the hydrophobic interactions with Trp179 and leads to the downward movement of the WPD loop, which forms an H-bond between Asp181 and Glu115. The formation of this H-bond constrains the WPD loop to its open conformation and thus inactivates PTP1B. The discovery of this allosteric mechanism provides an overall view of the regulation of PTP1B, which is an important insight for the design of potent allosteric PTP1B inhibitors.     

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 (WPD_closed) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPD_open) 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 WPD_closed 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 WPD_open conformation. When the active site water molecule network that is a part of the free WPD_closed 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. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

A targeted molecular dynamics study of WPD loop movement in PTP1B

Targeted molecular dynamics was used to examine the mechanism of WPD loop closure in PTP1B, which is essential for the activity of the enzyme. Two important regions are identified: the R-loop (residues 113-118), which assists in substrate binding, and the S-loop (residues 198-209), which undergoes a conformational change that appears to be vital for the movement of the WPD loop. The S-loop is adjacent to the alpha3-helix, and its conformational change is coupled with a change of interactions between the alpha3- and alpha7-helices. This latter observation is of particular interest in connection with a novel class of allosteric inhibitors of PTP1B [Wiesmann et al., Nat. Struc. Mol. Biol. 11 (2004) 730-737]. These compounds prevent the closure of the WPD loop, forcing the enzyme to remain in a catalytically inactive conformation, by blocking the rearrangement of the alpha3-helix relative to the alpha7-helix.     

Differential oxidation of protein-tyrosine phosphatases

Oxidation is emerging as an important regulatory mechanism of protein-tyrosine phosphatases (PTPs). Here we report that PTPs are differentially oxidized, and we provide evidence for the underlying mechanism. The membrane-proximal RPTPalpha-D1 was catalytically active but not readily oxidized as assessed by immunoprobing with an antibody that recognized oxidized catalytic site cysteines in PTPs (oxPTPs). In contrast, the membrane-distal RPTPalpha-D2, a poor PTP, was readily oxidized. Oxidized catalytic site cysteines in PTP immunoprobing and mass spectrometry demonstrated that mutation of two residues in the Tyr(P) loop and the WPD loop that reverse catalytic activity of RPTPalpha-D1 and RPTPalpha-D2 also reversed oxidizability, suggesting that oxidizability and catalytic activity are coupled. However, catalytically active PTP1B and LAR-D1 were readily oxidized. Oxidizability was strongly dependent on pH, indicating that the microenvironment of the catalytic cysteine has an important role. Crystal structures of PTP domains demonstrated that the orientation of the absolutely conserved PTP loop arginine correlates with oxidizability of PTPs, and consistently, RPTPmu-D1, with a similar conformation as RPTPalpha-D1, was not readily oxidized. In conclusion, PTPs are differentially oxidized at physiological pH and H(2)O(2) concentrations, and the PTP loop arginine is an important determinant for susceptibility to oxidation.     

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. I. ligand-induced changes in the protein motions.

Activity of enzymes, such as protein tyrosine phosphatases (PTPs), is often associated with structural changes in the enzyme, resulting in selective and stereospecific reactions with the substrate. To investigate the effect of a substrate on the motions occurring in PTPs, we have performed molecular dynamics simulations of PTP1B and PTP1B complexed with a high-affinity peptide DADEpYL, where pY stands for phosphorylated tyrosine. The peptide sequence is derived from the epidermal growth factor receptor (EGFR988-993). Simulations were performed in water for 1 ns, and the concerted motions in the protein were analyzed using the essential dynamics technique. Our results indicate that the predominately internal motions in PTP1B occur in a subspace of only a few degrees of freedom. Upon substrate binding, the flexibility of the protein is reduced by approximately 10%. The largest effect is found in the protein region, where the N-terminal of the substrate is located, and in the loop region Val198-Gly209. Displacements in the latter loop are associated with the motions in the WPD loop, which contains a catalytically important aspartic acid. Estimation of the pKa of the active-site cysteine along the trajectory indicates that structural inhomogeneity causes the pKa to vary by approximately +/-1 pKa unit. In agreement with experimental observations, the active-site cysteine is negatively charged at physiological pH.     

Single replacement constructs of all hydroxyl, basic, and acidic amino acids identify new function and structure-sensitive regions of the mitochondrial phosphate transport protein

The phosphate transport protein (PTP) catalyzes the proton cotransport of phosphate into the mitochondrial matrix. It functions as a homodimer, and thus residues of the phosphate and proton pores are somewhat scattered throughout the primary sequence. With 71 new single mutation per subunit PTPs, all its hydroxyl, basic, and acidic residues have now been replaced to identify these essential residues. We assayed the initial rate of pH gradient-dependent unidirectional phosphate transport activity and the liposome incorporation efficiency (LIE) of these mutants. Single mutations of Thr79, Tyr83, Lys90, Tyr94, and Lys98 inactivate transport. The spacings between these residues imply that they are located along the same face of transmembrane (TM) helix B, requiring an extension of its current model C-terminal domain by 10 residues. This extension superimposes very well onto the shorter bovine PTP helix B, leaving a 15-residue hydrophobic extension of the yeast helix B N-terminus. This is similar to the helix D and F regions of the yeast PTP. Only one transport-inhibiting mutation is located within loops: Ser158Thr in the matrix loop between helices C and D. All other transport-inhibiting mutations are located within the TM helices. Mutations that yield LIEs of <6% are all, except for four, within helices. The four exceptions are Tyr12Ala near the PTP N-terminus and Arg159Ala, Glu163Gln, and Glu164Gln in the loop between helices C and D. The PTP C-terminal segment beyond Thr214 at the N-terminus of helix E has 11 mutations with LIEs >20% and none with LIE <6%. Mutations with LIEs >20% are located near the ends of all the TM helices except TM helix D. Only a few mutations alter PTP structure (LIE) and also affect PTP transport activity. A novel observation is that Ser4Ala blocks the formation of PTP bacterial inclusion bodies.     

Visualizing active-site dynamics in single crystals of HePTP: opening of the WPD loop involves coordinated movement of the E loop

Phosphotyrosine hydrolysis by protein tyrosine phosphatases (PTPs) involves substrate binding by the PTP loop and closure over the active site by the WPD loop. The E loop, located immediately adjacent to the PTP and WPD loops, is conserved among human PTPs in both sequence and structure, yet the role of this loop in substrate binding and catalysis is comparatively unexplored. Hematopoietic PTP (HePTP) is a member of the kinase interaction motif (KIM) PTP family. Compared to other PTPs, KIM-PTPs have E loops that are unique in both sequence and structure. In order to understand the role of the E loop in the transition between the closed state and the open state of HePTP, we identified a novel crystal form of HePTP that allowed the closed-state-to-open-state transition to be observed within a single crystal form. These structures, which include the first structure of the HePTP open state, show that the WPD loop adopts an ‘atypically open’ conformation and, importantly, that ligands can be exchanged at the active site, which is critical for HePTP inhibitor development. These structures also show that tetrahedral oxyanions bind at a novel secondary site and function to coordinate the PTP, WPD, and E loops. Finally, using both structural and kinetic data, we reveal a novel role for E-loop residue Lys182 in enhancing HePTP catalytic activity through its interaction with Asp236 of the WPD loop, providing the first evidence for the coordinated dynamics of the WPD and E loops in the catalytic cycle, which, as we show, is relevant to multiple PTP families.     

Identification of flavonolignans from Silybum marianum seeds as allosteric protein tyrosine phosphatase 1B inhibitors

Protein tyrosine phosphatase 1B (PTP1B) is an attractive molecular target for anti-diabetes, anti-obesity, and anti-cancer drug development. From the seeds of Silybum marianum, nine flavonolignans, namely, silybins A, B (1, 2), isosilybins A, B (3, 4), silychristins A, B (5, 6), isosilychristin A (7), dehydrosilychristin A (8), and silydianin (11) were identified as a novel class of natural PTP1B inhibitors (IC₅₀ 1.3 7–23.87?µM). Analysis of structure–activity relationship suggested that the absolute configurations at C-7" and C-8" greatly affected the PTP1B inhibitory activity. Compounds 1–5 were demonstrated to be non-competitive inhibitors of PTP1B based on kinetic analyses. Molecular docking simulations resulted that 1–5 docked into the allosteric site, including α3, α6, and α7 helix of PTP1B. At a concentration inhibiting PTP1B completely, compounds 1–5 moderately inhibited VHR and SHP-2, and weakly inhibited TCPTP and SHP-1. These results suggested the potentiality of these PTP1B inhibitors as lead compounds for further drug developments.     

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.     

Dynamics of the WPD loop of the Yersinia protein tyrosine phosphatase

The bacterial protein tyrosine phosphatase YopH is an essential virulence determinant in Yersinia spp., causing gastrointestinal diseases and the plague. Like eukaryotic PTPases, YopH catalyzes the hydrolysis of the phosphate moiety of phosphotyrosine within a highly conserved binding pocket, which is also characterized by the closure of the so-called "WPD loop" upon ligand binding. In this study, we investigate the conformational changes and dynamics of the WPD loop by molecular dynamics simulations. Consistent with experimental observations, our simulations show that the WPD loop of YopH is intrinsically flexible and fluctuates between the open and closed conformation with a frequency of approximately 4 ns for the apo, native protein. The region of helix alpha4 spanning loop 384-392, which has been revealed experimentally as a second substrate-binding site in YopH, is found to be highly associated with the WPD loop, stabilizing it in the closed, active conformation, and providing a structural basis for the cooperation of the second-substrate binding site in substrate recognition. Loop L4 (residues 323-327) is shown to be involved in a parallel, correlated motion mode with the WPD loop that contributes the stabilization of a more extended open conformation. In addition, we have simulated the loop reopening in the ligand-bound protein complex by applying the locally enhanced sampling method. Finally, the dynamic behavior of the WPD loop for the C403S mutant differs from the wild-type YopH remarkably. These results shed light on the role of the WPD loop in PTPase-mediated catalysis, and are useful in structure-based design for novel, selective YopH inhibitors as antibacterial drugs.     

The molecular details of WPD-loop movement differ in the protein-tyrosine phosphatases YopH and PTP1B

The movement of a conserved protein loop (the WPD-loop) is important in catalysis by protein tyrosine phosphatases (PTPs). Using kinetics, isotope effects, and X-ray crystallography, the different effects arising from mutation of the conserved tryptophan in the WPD-loop were compared in two PTPs, the human PTP1B, and the bacterial YopH from Yersinia. Mutation of the conserved tryptophan in the WPD-loop to phenylalanine has a negligible effect on k(cat) in PTP1B and full loop movement is maintained. In contrast, the corresponding mutation in YopH reduces k(cat) by two orders of magnitude and the WPD loop locks in an intermediate position, disabling general acid catalysis. During loop movement the indole moiety of the WPD-loop tryptophan moves in opposite directions in the two enzymes. Comparisons of mammalian and bacterial PTPs reveal differences in the residues forming the hydrophobic pocket surrounding the conserved tryptophan. Thus, although WPD-loop movement is a conserved feature in PTPs, differences exist in the molecular details, and in the tolerance to mutation, in PTP1B compared to YopH. Despite high structural similarity of the active sites in both WPD-loop open and closed conformations, differences are identified in the molecular details associated with loop movement in PTPs from different organisms.     

Rab-subfamily-specific regions of Ypt7p are structurally different from other RabGTPases

The GTPase Ypt7p from S. cerevisiae is involved in late endosome-to-vacuole transport and homotypic vacuole fusion. We present crystal structures of the GDP- and GppNHp-bound conformation of Ypt7p solved at 1.35 and 1.6 A resolution, respectively. Despite the similarity of the overall structure to other Ypt/Rab proteins, Ypt7p displays small but significant differences. The Ypt7p-specific residues Tyr33 and Tyr37 cause a difference in the main chain trace of the RabSF2 region and form a characteristic surface epitope. Ypt7p*GppNHp does not display the helix alpha2, characteristic of the Ras-superfamily, but instead possess an extended loop L4/L5. Due to insertions in loops L3 and L7, the neighboring RabSF1 and RabSF4 regions are different in their conformations to those of other Ypt/Rab proteins.     

The structure of apo protein-tyrosine phosphatase 1B C215S mutant: more than just an S --> O change

Protein-tyrosine phosphatases catalyze the hydrolysis of phosphate monoesters via a two-step mechanism involving a covalent phospho-enzyme intermediate. Biochemical and site-directed mutagenesis experiments show that the invariant Cys residue present in the PTPase signature motif (H/V)CX(5)R(S/T) (i.e., C215 in PTP1B) is absolutely required for activity. Mutation of the invariant Cys to Ser results in a catalytically inactive enzyme, which still is capable of binding substrates and inhibitors. Although it often is assumed that substrate-trapping mutants such as the C215S retain, in solution, the structural and binding properties of wild-type PTPases, significant differences have been found in the few studies that have addressed this issue, suggesting that the mutation may lead to structural/conformational alterations in or near the PTP1B binding site. Several crystal structures of apo-WT PTP1B, and of WT- and C215S-mutant PTP1B in complex with different ligands are available, but no structure of the apo-PTP1B C215S has ever been reported. In all previously reported structures, residues of the PTPase signature motif have an identical conformation, while residues of the WPD loop (a surface loop which includes the catalytic Asp) assume a different conformation in the presence or absence of ligand. These observations led to the hypothesis that the different spectroscopic and thermodynamic properties of the mutant protein may be the result of a different conformation for the WPD loop. We report here the structure of the apo-PTP1B C215S mutant, which reveals that, while the WPD loop is in the open conformation observed in the apo WT enzyme crystal structure, the residues of the PTPases signature motif are in a dramatically different conformation. These results provide a structural basis for the differences in spectroscopic properties and thermodynamic parameters in inhibitor binding observed for the wild-type and mutant enzymes.     

Small molecule peptidomimetics containing a novel phosphotyrosine bioisostere inhibit protein tyrosine phosphatase 1B and augment insulin action

Protein tyrosine phosphatase 1B (PTP1B) attenuates insulin signaling by catalyzing dephosphorylation of insulin receptors (IR) and is an attractive target of potential new drugs for treating the insulin resistance that is central to type II diabetes. Several analogues of cholecystokinin(26)(-)(33) (CCK-8) were found to be surprisingly potent inhibitors of PTP1B, and a common N-terminal tripeptide, N-acetyl-Asp-Tyr(SO(3)H)-Nle-, was shown to be necessary and sufficient for inhibition. This tripeptide was modified to reduce size and peptide character, and to replace the metabolically unstable sulfotyrosyl group. This led to the discovery of a novel phosphotyrosine bioisostere, 2-carboxymethoxybenzoic acid, and to analogues that were >100-fold more potent than the CCK-8 analogues and >10-fold selective for PTP1B over two other PTP enzymes (LAR and SHP-2), a dual specificity phosphatase (cdc25b), and a serine/threonine phosphatase (calcineurin). These inhibitors disrupted the binding of PTP1B to activated IR in vitro and prevented the loss of tyrosine kinase (IRTK) activity that accompanied PTP1B-catalyzed dephosphorylation of IR. Introduction of these poorly cell permeant inhibitors into insulin-treated cells by microinjection (oocytes) or by esterification to more lipophilic proinhibitors (3T3-L1 adipocytes and L6 myocytes) resulted in increased potency, but not efficacy, of insulin. In some instances, PTP1B inhibitors were insulin-mimetic, suggesting that in unstimulated cells PTP1B may suppress basal IRTK activity. X-ray crystallography of PTP1B-inhibitor complexes revealed that binding of an inhibitor incorporating phenyl-O-malonic acid as a phosphotyrosine bioisostere occurred with the mobile WPD loop in the open conformation, while a closely related inhibitor with a 2-carboxymethoxybenzoic acid bioisostere bound with the WPD loop closed, perhaps accounting for its superior potency. These CCK-derived peptidomimetic inhibitors of PTP1B represent a novel template for further development of potent, selective inhibitors, and their cell activity further justifies the selection of PTP1B as a therapeutic target.     

Novel inter- and intrasubunit contacts between transport-relevant residues of the homodimeric mitochondrial phosphate transport protein

Ser158 is located near the middle of the matrix loop connecting transmembrane helices C and D of the mitochondrial phosphate transport protein (PTP). The mutant Ser158Thr PTP is transport-inactive. His32 is located near the middle of transmembrane helix A and Thr79 is located 5 residues away from transmembrane helix B and its N-terminal (matrix end). Single site mutant PTPs that have either residue replaced with Ala are transport-inactive. Based on the high resolution structure of a subunit of the bovine ADP/ATP translocase, on sequence similarities between members of the mitochondrial transport protein family, and on the PTP subunit/subunit contact site between transmembrane A helices, it is now suggested that the Ser158 site is at the PTP subunit/subunit contact site. This contact site is essential for keeping the transport cycles catalyzed by the two PTP subunits 180 degrees out of phase. The data also suggest that His32 and Thr79 of the same subunit interact and couple the phosphate and the proton transport paths.     

Synthesis,biological evaluation,and molecular docking study of novel allyl-retrochalcones as a new class of protein tyrosine phosphatase 1B inhibitors

We describe herein the design, synthesis, and biological evaluation of a series of novel protein tyrosine phosphatase 1B (PTP1B) inhibitor retrochalcones having an allyl chain at the C-5 position of their B ring. Biological screening results showed that the majority of these compounds exhibited an inhibitory activity against PTP1B. Thus, preliminary structure-activity relationship (SAR) and quantitative SAR analyses were conducted. Among the compounds, 23 was the most potent inhibitor, exhibiting the highest in vitro inhibitory activity against PTP1B with an IC₅₀ of 0.57?µM. Moreover, it displayed a significant hepatoprotective property via activation of the IR pathway in type 2 diabetic db/db mice. In addition, the results of our docking study showed that 23, as a specific inhibitor of PTP1B, effectively transformed the WPD loop from “close” to “open” in the active site. These results may reveal suitable compounds for the development of PTP1B inhibitors.     

The Distinct Conformational Dynamics of K-Ras and H-Ras A59G

Ras proteins regulate signaling cascades crucial for cell proliferation and differentiation by switching between GTP- and GDP-bound conformations. Distinct Ras isoforms have unique physiological functions with individual isoforms associated with different cancers and developmental diseases. Given the small structural differences among isoforms and mutants, it is currently unclear how these functional differences and aberrant properties arise. Here we investigate whether the subtle differences among isoforms and mutants are associated with detectable dynamical differences. Extensive molecular dynamics simulations reveal that wild-type K-Ras and mutant H-Ras A59G are intrinsically more dynamic than wild-type H-Ras. The crucial switch 1 and switch 2 regions along with loop 3, helix 3, and loop 7 contribute to this enhanced flexibility. Removing the gamma-phosphate of the bound GTP from the structure of A59G led to a spontaneous GTP-to-GDP conformational transition in a 20-ns unbiased simulation. The switch 1 and 2 regions exhibit enhanced flexibility and correlated motion when compared to non-transitioning wild-type H-Ras over a similar timeframe. Correlated motions between loop 3 and helix 5 of wild-type H-Ras are absent in the mutant A59G reflecting the enhanced dynamics of the loop 3 region. Taken together with earlier findings, these results suggest the existence of a lower energetic barrier between GTP and GDP states of the mutant. Molecular dynamics simulations combined with principal component analysis of available Ras crystallographic structures can be used to discriminate ligand- and sequence-based dynamic perturbations with potential functional implications. Furthermore, the identification of specific conformations associated with distinct Ras isoforms and mutants provides useful information for efforts that attempt to selectively interfere with the aberrant functions of these species.     

Solution NMR and computer simulation studies of active site loop motion in triosephosphate isomerase

Solution NMR spin relaxation experiments and classical MD simulations are used to study the dynamics of triosephosphate isomerase (TIM) in complex with glycerol 3-phosphate (G3P). Three regions in TIM exhibit conformational transitions on the micros-ms time scale as detected by chemical exchange broadening effects in NMR spectroscopy: residue Lys 84 on helix C, located at the dimeric interface; active site loop 6; and helix G. The results indicate that the conformational exchange process affecting the residues of loop 6 is the correlated opening and closing of the loop. Distinct processes are responsible for the chemical exchange linebroadening observed in the other regions of TIM. MD simulations confirm that motions of individual residues within the active site loop are correlated and suggest that the chemical exchange processes observed for residues in helix G arise from transitions between 3(10)- and alpha-helical structures. The results of the joint NMR and MD study provide global insight into the role of conformational dynamic processes in the function of TIM.     

Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics

Recent functional studies reported on human adult hemoglobin (HbA) show that heterotropic effector-linked tertiary structural changes are primarily responsible for modulating the oxygen affinity of hemoglobin. We present the results of 6-ns molecular dynamics simulations performed to gain insights into the dynamical and structural details of these effector-linked tertiary changes. All-atom simulations were carried out on a series of models generated for T- and R-state HbA, and for 2,3-diphosphoglycerate-bound models. Cross-correlation analyses identify both intra- and intersubunit correlated motions that are perturbed by the presence of the effector. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. It is found that HbA in both quaternary states exists as ensembles of tertiary conformations that introduce dynamic heterogeneity in the protein. 2,3-Diphosphoglycerate induces significant perturbations in the fluctuations of both HbA states that translate into the protein visiting different tertiary conformations within each quaternary state. The analysis reveals that the presence of the effector affects the most important components of HbA motions and that heterotropic effectors modify the overall dynamics of the quaternary equilibrium via tertiary changes occurring in regions where conserved functionally significant residues are located, namely in the loop regions between helices C and E, E and F, and F and G, and in concerted helix motions. The changes are not apparent when comparing the available x-ray crystal structures in the presence and absence of effector, but are striking when comparing the respective dynamic tertiary conformations of the R and T tetramers.