The method of diagonalization in a mixed basis (DIMB) that was published previously (Mouawad.), L... more The method of diagonalization in a mixed basis (DIMB) that was published previously (Mouawad.), L. and Perahia D., Biopolymers 33, 599, 1993), allows the computation of the low-frequency vibrational modes for large macromolecules. Improvements to this method are presented here, namely the single and double truncation window techniques. The best convergence rate is obtained with the double truncation windows, which couple most efficiently the parts of the macromolecule which are far in sequence but close in space. Both methods were applied to the T-state of hemoglobin, to compare their efficiency. The resulting modes are analyzed in order to study the pathways between T- and R-states of this protein. They show that the quaternary conformational are mainly due to one mode at 2 cm-1.
Iron (II) basket-handle porphyrins (BHP) are a series of encumbered heme models designed several ... more Iron (II) basket-handle porphyrins (BHP) are a series of encumbered heme models designed several years ago to mimic the ligand binding site of hemoproteins. Contrary to expectations, kinetic investigations have revealed that the k(on) rates for CO and/or O2 binding were only marginally affected by the assumed central steric hindrance of the iron atom. Thus, it was hypothesized that the internal dynamics of the molecule might be at the origin of the poor steric protection. To address this issue, measurements of nuclear magnetic resonance relaxation rates, fluorescence anisotropy experiments, and molecular dynamics simulations were undertaken. The size of BHP is small enough to allow the simulation in explicit chloroform with an almost complete sampling of the conformational space. The order parameters calculated from the MD trajectory compare well with the NMR experimental data and the predicted rotational correlation time corresponding to the Brownian motion of the molecule is in good agreement with the fluorescence measurements. Moreover, combining the results obtained using the three techniques allows the attribution of each internal NMR correlation time to a particular internal motion, revealing that even such medium-sized molecules are able to display quite complex internal dynamics. In particular, the handle phenyls that were assumed to sandwich the porphyrin have in fact a vanishing probability to be found in the proximity of the iron atom. They are therefore unable to reduce ligand accessibility significantly, which may explain the behavior of the k(on) rates.
Horse muscle phosphoglycerate kinase, like other mammalian phosphoglycerate kinases, contains sev... more Horse muscle phosphoglycerate kinase, like other mammalian phosphoglycerate kinases, contains seven cysteine residues of which two react rapidly with 5,5'-dithio-bis(2-nitrobenzoate) (Nbs2) following second-order kinetics (k = 640 M-1.s-1). Selective cyanylation of the fast-reacting cysteines, followed by chemical cleavage and subsequent sodium dodecyl sulfate/polyacrylamide gel electrophoresis analysis of the resulting polypeptides, suggested that these cysteines are at positions 378 and 379. Cysteine residues were introduced into yeast phosphoglycerate kinase by site-directed mutagenesis. Mutant enzymes, each containing only one cysteine residue at position 364, 376, or 377, were constructed from a mutant devoid of cysteine (Cys97----Ala). In the last two mutants, the cysteines were at positions corresponding to Cys378 and Cys379, respectively, in the horse muscle enzyme. The chemical reactivity of the cysteine groups in these latter two yeast mutant enzymes was similar to that of the fast-reacting cysteines in the horse muscle enzyme. Furthermore, they were similarly modified upon substrate binding. All these data demonstrate unambiguously that the fast-reacting cysteines in the horse muscle enzyme are Cys378 and Cys379.
Protein kinase casein kinase 2 (CK2) is a serine/threonine kinase with evidence of implication in... more Protein kinase casein kinase 2 (CK2) is a serine/threonine kinase with evidence of implication in growth dysregulation and apoptosis resistance, making it a relevant target for cancer therapy. Several CK2 inhibitors have been developed showing variable efficiency, emphasizing the need to expand the chemical diversity of those inhibitors. We report the identification and characterization of 2,8-difurandicarboxylic acid derivatives as a new class of nanomolar ATP-competitive inhibitors. Selectivity profiling pointed out proviral insertion Moloney virus kinases (Pim kinases) as the only other kinases that are significantly inhibited. By combining structure-activity relationship analysis with structural determination, we were able to determine the binding mode of these inhibitors for both kinases and to explain their strong inhibitory potency. Essential chemical features necessary for activity on both kinases were then identified. The described compounds are not cell permeable: however, they could provide a lead for developing novel inhibitors usable also in vivo. Given the similar but not redundant pathophysiological functions of CK2 and Pim family members, such inhibitors would provide new attractive leads for targeted cancer therapy. This work highlights that 2 functionally related kinases from different kinome branches display exquisite sensitivity to a common inhibitor.
A microtubule is a dynamic system formed of αβ-tubulins. The presence of nonhydrolyzable guanosin... more A microtubule is a dynamic system formed of αβ-tubulins. The presence of nonhydrolyzable guanosine-5'-triphosphate (GTP)/guanosine diphosphate (GDP) on the β-tubulins provokes microtubule polymerization/depolymerization. Despite the large number of experimental studies of this dynamical process, its mechanism is still unclear. To provide insights into this mechanism we studied the first depolymerization steps of GDP/GTP-bound microtubules by normal-mode analysis with the all-atom model. We also constructed a depolymerizing microtubule and compared it to cryo-electron microscopy tomograms (cyro-ET). The results show that during depolymerization, the protofilaments not only curve but twist to weaken their lateral interactions. These interactions are stabilized by GTP, but not evenly. Not all of the interface residues are of equal importance: five of them, belonging to the H2-S3 loop, play a special role; acting as a lock whose key is the γ-phosphate of GTP. Sequence alignments of several tubulins confirm the importance of these residues.
Several mutants of yeast phosphoglycerate kinase, each containing only one internal cysteine resi... more Several mutants of yeast phosphoglycerate kinase, each containing only one internal cysteine residue, were constructed from a single mutant devoid of cysteine. These cysteines were introduced as local conformational probes in selected buried positions. The enzyme activity, conformational characteristics and stability indicated that the mutations introduced only small perturbations in the molecule. The folding-unfolding process mediated by guanidine hydrochloride under equilibrium conditions was studied by following the variations in ellipticity and the reactivity of the cysteine residue towards 5,5'-dithiobis(nitrobenzoate). The process was found to be reversible except for mutant C97A, V49C, suggesting that this region located in helix I might be crucial in determining an intermediate on the folding pathway. The transitions obtained by the two signals did not coincide, indicating that the local structures, in several parts inside the molecule, are more sensitive to the denaturant than the overall conformation.
This paper presents the comparative comprehensive analysis of NMR structural parameters (NOEs, sc... more This paper presents the comparative comprehensive analysis of NMR structural parameters (NOEs, scalar coupling, chemical shifts) of toxin gamma, a cardiotoxin isolated from the venom of Naja nigricollis, and three chemical derivatives, i.e. the 2-nitrophenylsulphonyl (NPS)-Trp11, 3-nitro-Tyr22 and 3-nitro-Tyr51 derivatives. In previous work, the chemical modifications of single side chains have suggested that these aromatic residues, in association with several lysine residues, contributed to the cytotoxicity of toxin gamma. Analysis of these results based on the refined solution structure of the toxin has resulted in the proposal of a conserved phospholipid binding site through which cardiotoxins are likely to interact with the membrane of target cells. The present work shows that modifications of either the tryptophan residue or the tyrosine residues, which are within or near the proposed binding site, have no influence on the three-dimensional structure of the protein. On the other hand, the proton exchange study of the backbone amides indicates that the structural core of the protein is destabilized in the three derivatives. This corresponds to a decrease of the overall stability of the protein as indicated by the comparative solvent denaturation study of the unmodified toxin gamma and the Trp11 derivative. More specifically, the dynamics of the three-stranded beta sheet, a part of the structural core, are highly perturbed by the chemical modifications. This sheet was previously proposed as a part of the phospholipid binding site of cardiotoxins. The dynamical perturbation of this site appears to be correlated with the decrease in toxicity of the chemical derivatives.
Aspartate transcarbamylase (ATCase) is an important control enzyme in the pyrimidine biosynthetic... more Aspartate transcarbamylase (ATCase) is an important control enzyme in the pyrimidine biosynthetic pathway in Escherichia coli. It is a classic example of an allosteric protein and has been extensively studied biochemically, kinetically and structurally. As yet, however, a detailed model for the cooperative transition between the tensed (T) and relaxed (R) forms of the protein does not exist. In this work we have calculated the low frequency normal modes of the CTP-ligated T-state of ATCase with the aim of identifying some of the motions that could be important in initiating the transition. The calculated modes, of frequencies lower than 5 per cm, produce root-mean-square coordinate deviations for the atoms which are a substantial fraction of those derived from the crystallographic B-factors. Some of the modes result in displacements which change the quaternary structure of the protein (in particular the elongation of the protein and the relative rotation of the subunits) in such a way that the R-state structure is approached. The implication of these mode motions for the overall T-->R transition process is discussed.
The normal mode analysis of human hemoglobin showed the presence in the deoxy T-state of one main... more The normal mode analysis of human hemoglobin showed the presence in the deoxy T-state of one main preferential direction that brings the structure close to the R-state, with a low-energy variation, while in the oxy R-state there are several modes that point towards the T-state, but with higher energy variations and less contribution to the transition. The displacement along a combination of normal modes, followed by energy minimization, starting from the R-state, did not allow one to obtain a structure significantly different from that of R, showing that the fully oxygenated hemoglobin is trapped in a deep and narrow potential energy minimum. On the contrary, starting from the deoxy T-state, the displacement along a combination of normal modes, followed by energy minimization, yielded an intermediate structure, that we designate Tmin(d1), which is closer to R; the normal modes of Tmin(d1) indicated that the potential energy minimum in the vicinity of this structure is as narrow as that of R but less deep. The procedure of displacement along the modes, followed by energy minimization, was applied to Tmin(d1), yielding Tmin(d2); then the procedure was repeated, yielding the intermediate structures Tmin(d3) and Tmin(d4). The structures Tmin(d2), Tmin(d3) and Tmin(d4) are not significantly different from each other, indicating that they are trapped in a narrow, deep energy minimum. This procedure revealed the existence of at least two intermediate sets of structures between T and R: the first one, Tmin(d1), is different from the T and R structures, while the second set, Tmin(d2), Tmin(d3) and Tmin(d4), is quaternary R-like and tertiary T-like, where the contacts at the interfaces alpha1 beta1 and alpha1 beta2 are R-like, and the alpha and beta heme environments are still T-like.
The proteins Arno and Gea2 of the Sec7 family can promote GDP-GTP exchange on Arf1, a small GTP-b... more The proteins Arno and Gea2 of the Sec7 family can promote GDP-GTP exchange on Arf1, a small GTP-binding protein, which coordinates coated vesicle formation for protein transport within the cell. Crystal structures of the essential Sec7 domain (Sec7d) of Gea2 in the free and Arf1-bound forms suggest that conformational dynamics of the Sec7d as well as those of the G-protein play a role in nucleotide exchange. Starting from a set of complementary crystal structures, we compared the collective movements of unbound Gea2 and Arno Sec7 domains, Arf1-GDP, and the Arf1-Gea2(Sec7d) nucleotide-free complex using normal modes analyses. In all unbound Sec7d analyses, significant low-energy movements were found to lead to closure of the hydrophobic groove towards the form seen in the Arf1-Gea2(Sec7d) complex, suggesting that groove closure is a general feature of the Sec7 family. Low-energy movements in Arf1-GDP implicate critical switch 1 and 2 residues which are coupled to modifications in the myristoylated N-terminal-helix binding site at the other end of the "interswitch" beta hairpin. It is suggested that Sec7d groove closure upon docking of the two molecules may permit extraction of switch 1 from Arf1-GDP and prepare the complex for movement of the interswitch, which is central to the membrane-linked exchange activity. Large-scale collective movements in the Arf1-Sec7d complex appear to participate in the insertion of the Sec7d Glu finger into the GDP binding site to promote actual nucleotide release.
The method of diagonalization in a mixed basis (DIMB) that was published previously (Mouawad.), L... more The method of diagonalization in a mixed basis (DIMB) that was published previously (Mouawad.), L. and Perahia D., Biopolymers 33, 599, 1993), allows the computation of the low-frequency vibrational modes for large macromolecules. Improvements to this method are presented here, namely the single and double truncation window techniques. The best convergence rate is obtained with the double truncation windows, which couple most efficiently the parts of the macromolecule which are far in sequence but close in space. Both methods were applied to the T-state of hemoglobin, to compare their efficiency. The resulting modes are analyzed in order to study the pathways between T- and R-states of this protein. They show that the quaternary conformational are mainly due to one mode at 2 cm-1.
Iron (II) basket-handle porphyrins (BHP) are a series of encumbered heme models designed several ... more Iron (II) basket-handle porphyrins (BHP) are a series of encumbered heme models designed several years ago to mimic the ligand binding site of hemoproteins. Contrary to expectations, kinetic investigations have revealed that the k(on) rates for CO and/or O2 binding were only marginally affected by the assumed central steric hindrance of the iron atom. Thus, it was hypothesized that the internal dynamics of the molecule might be at the origin of the poor steric protection. To address this issue, measurements of nuclear magnetic resonance relaxation rates, fluorescence anisotropy experiments, and molecular dynamics simulations were undertaken. The size of BHP is small enough to allow the simulation in explicit chloroform with an almost complete sampling of the conformational space. The order parameters calculated from the MD trajectory compare well with the NMR experimental data and the predicted rotational correlation time corresponding to the Brownian motion of the molecule is in good agreement with the fluorescence measurements. Moreover, combining the results obtained using the three techniques allows the attribution of each internal NMR correlation time to a particular internal motion, revealing that even such medium-sized molecules are able to display quite complex internal dynamics. In particular, the handle phenyls that were assumed to sandwich the porphyrin have in fact a vanishing probability to be found in the proximity of the iron atom. They are therefore unable to reduce ligand accessibility significantly, which may explain the behavior of the k(on) rates.
Horse muscle phosphoglycerate kinase, like other mammalian phosphoglycerate kinases, contains sev... more Horse muscle phosphoglycerate kinase, like other mammalian phosphoglycerate kinases, contains seven cysteine residues of which two react rapidly with 5,5'-dithio-bis(2-nitrobenzoate) (Nbs2) following second-order kinetics (k = 640 M-1.s-1). Selective cyanylation of the fast-reacting cysteines, followed by chemical cleavage and subsequent sodium dodecyl sulfate/polyacrylamide gel electrophoresis analysis of the resulting polypeptides, suggested that these cysteines are at positions 378 and 379. Cysteine residues were introduced into yeast phosphoglycerate kinase by site-directed mutagenesis. Mutant enzymes, each containing only one cysteine residue at position 364, 376, or 377, were constructed from a mutant devoid of cysteine (Cys97----Ala). In the last two mutants, the cysteines were at positions corresponding to Cys378 and Cys379, respectively, in the horse muscle enzyme. The chemical reactivity of the cysteine groups in these latter two yeast mutant enzymes was similar to that of the fast-reacting cysteines in the horse muscle enzyme. Furthermore, they were similarly modified upon substrate binding. All these data demonstrate unambiguously that the fast-reacting cysteines in the horse muscle enzyme are Cys378 and Cys379.
Protein kinase casein kinase 2 (CK2) is a serine/threonine kinase with evidence of implication in... more Protein kinase casein kinase 2 (CK2) is a serine/threonine kinase with evidence of implication in growth dysregulation and apoptosis resistance, making it a relevant target for cancer therapy. Several CK2 inhibitors have been developed showing variable efficiency, emphasizing the need to expand the chemical diversity of those inhibitors. We report the identification and characterization of 2,8-difurandicarboxylic acid derivatives as a new class of nanomolar ATP-competitive inhibitors. Selectivity profiling pointed out proviral insertion Moloney virus kinases (Pim kinases) as the only other kinases that are significantly inhibited. By combining structure-activity relationship analysis with structural determination, we were able to determine the binding mode of these inhibitors for both kinases and to explain their strong inhibitory potency. Essential chemical features necessary for activity on both kinases were then identified. The described compounds are not cell permeable: however, they could provide a lead for developing novel inhibitors usable also in vivo. Given the similar but not redundant pathophysiological functions of CK2 and Pim family members, such inhibitors would provide new attractive leads for targeted cancer therapy. This work highlights that 2 functionally related kinases from different kinome branches display exquisite sensitivity to a common inhibitor.
A microtubule is a dynamic system formed of αβ-tubulins. The presence of nonhydrolyzable guanosin... more A microtubule is a dynamic system formed of αβ-tubulins. The presence of nonhydrolyzable guanosine-5'-triphosphate (GTP)/guanosine diphosphate (GDP) on the β-tubulins provokes microtubule polymerization/depolymerization. Despite the large number of experimental studies of this dynamical process, its mechanism is still unclear. To provide insights into this mechanism we studied the first depolymerization steps of GDP/GTP-bound microtubules by normal-mode analysis with the all-atom model. We also constructed a depolymerizing microtubule and compared it to cryo-electron microscopy tomograms (cyro-ET). The results show that during depolymerization, the protofilaments not only curve but twist to weaken their lateral interactions. These interactions are stabilized by GTP, but not evenly. Not all of the interface residues are of equal importance: five of them, belonging to the H2-S3 loop, play a special role; acting as a lock whose key is the γ-phosphate of GTP. Sequence alignments of several tubulins confirm the importance of these residues.
Several mutants of yeast phosphoglycerate kinase, each containing only one internal cysteine resi... more Several mutants of yeast phosphoglycerate kinase, each containing only one internal cysteine residue, were constructed from a single mutant devoid of cysteine. These cysteines were introduced as local conformational probes in selected buried positions. The enzyme activity, conformational characteristics and stability indicated that the mutations introduced only small perturbations in the molecule. The folding-unfolding process mediated by guanidine hydrochloride under equilibrium conditions was studied by following the variations in ellipticity and the reactivity of the cysteine residue towards 5,5'-dithiobis(nitrobenzoate). The process was found to be reversible except for mutant C97A, V49C, suggesting that this region located in helix I might be crucial in determining an intermediate on the folding pathway. The transitions obtained by the two signals did not coincide, indicating that the local structures, in several parts inside the molecule, are more sensitive to the denaturant than the overall conformation.
This paper presents the comparative comprehensive analysis of NMR structural parameters (NOEs, sc... more This paper presents the comparative comprehensive analysis of NMR structural parameters (NOEs, scalar coupling, chemical shifts) of toxin gamma, a cardiotoxin isolated from the venom of Naja nigricollis, and three chemical derivatives, i.e. the 2-nitrophenylsulphonyl (NPS)-Trp11, 3-nitro-Tyr22 and 3-nitro-Tyr51 derivatives. In previous work, the chemical modifications of single side chains have suggested that these aromatic residues, in association with several lysine residues, contributed to the cytotoxicity of toxin gamma. Analysis of these results based on the refined solution structure of the toxin has resulted in the proposal of a conserved phospholipid binding site through which cardiotoxins are likely to interact with the membrane of target cells. The present work shows that modifications of either the tryptophan residue or the tyrosine residues, which are within or near the proposed binding site, have no influence on the three-dimensional structure of the protein. On the other hand, the proton exchange study of the backbone amides indicates that the structural core of the protein is destabilized in the three derivatives. This corresponds to a decrease of the overall stability of the protein as indicated by the comparative solvent denaturation study of the unmodified toxin gamma and the Trp11 derivative. More specifically, the dynamics of the three-stranded beta sheet, a part of the structural core, are highly perturbed by the chemical modifications. This sheet was previously proposed as a part of the phospholipid binding site of cardiotoxins. The dynamical perturbation of this site appears to be correlated with the decrease in toxicity of the chemical derivatives.
Aspartate transcarbamylase (ATCase) is an important control enzyme in the pyrimidine biosynthetic... more Aspartate transcarbamylase (ATCase) is an important control enzyme in the pyrimidine biosynthetic pathway in Escherichia coli. It is a classic example of an allosteric protein and has been extensively studied biochemically, kinetically and structurally. As yet, however, a detailed model for the cooperative transition between the tensed (T) and relaxed (R) forms of the protein does not exist. In this work we have calculated the low frequency normal modes of the CTP-ligated T-state of ATCase with the aim of identifying some of the motions that could be important in initiating the transition. The calculated modes, of frequencies lower than 5 per cm, produce root-mean-square coordinate deviations for the atoms which are a substantial fraction of those derived from the crystallographic B-factors. Some of the modes result in displacements which change the quaternary structure of the protein (in particular the elongation of the protein and the relative rotation of the subunits) in such a way that the R-state structure is approached. The implication of these mode motions for the overall T-->R transition process is discussed.
The normal mode analysis of human hemoglobin showed the presence in the deoxy T-state of one main... more The normal mode analysis of human hemoglobin showed the presence in the deoxy T-state of one main preferential direction that brings the structure close to the R-state, with a low-energy variation, while in the oxy R-state there are several modes that point towards the T-state, but with higher energy variations and less contribution to the transition. The displacement along a combination of normal modes, followed by energy minimization, starting from the R-state, did not allow one to obtain a structure significantly different from that of R, showing that the fully oxygenated hemoglobin is trapped in a deep and narrow potential energy minimum. On the contrary, starting from the deoxy T-state, the displacement along a combination of normal modes, followed by energy minimization, yielded an intermediate structure, that we designate Tmin(d1), which is closer to R; the normal modes of Tmin(d1) indicated that the potential energy minimum in the vicinity of this structure is as narrow as that of R but less deep. The procedure of displacement along the modes, followed by energy minimization, was applied to Tmin(d1), yielding Tmin(d2); then the procedure was repeated, yielding the intermediate structures Tmin(d3) and Tmin(d4). The structures Tmin(d2), Tmin(d3) and Tmin(d4) are not significantly different from each other, indicating that they are trapped in a narrow, deep energy minimum. This procedure revealed the existence of at least two intermediate sets of structures between T and R: the first one, Tmin(d1), is different from the T and R structures, while the second set, Tmin(d2), Tmin(d3) and Tmin(d4), is quaternary R-like and tertiary T-like, where the contacts at the interfaces alpha1 beta1 and alpha1 beta2 are R-like, and the alpha and beta heme environments are still T-like.
The proteins Arno and Gea2 of the Sec7 family can promote GDP-GTP exchange on Arf1, a small GTP-b... more The proteins Arno and Gea2 of the Sec7 family can promote GDP-GTP exchange on Arf1, a small GTP-binding protein, which coordinates coated vesicle formation for protein transport within the cell. Crystal structures of the essential Sec7 domain (Sec7d) of Gea2 in the free and Arf1-bound forms suggest that conformational dynamics of the Sec7d as well as those of the G-protein play a role in nucleotide exchange. Starting from a set of complementary crystal structures, we compared the collective movements of unbound Gea2 and Arno Sec7 domains, Arf1-GDP, and the Arf1-Gea2(Sec7d) nucleotide-free complex using normal modes analyses. In all unbound Sec7d analyses, significant low-energy movements were found to lead to closure of the hydrophobic groove towards the form seen in the Arf1-Gea2(Sec7d) complex, suggesting that groove closure is a general feature of the Sec7 family. Low-energy movements in Arf1-GDP implicate critical switch 1 and 2 residues which are coupled to modifications in the myristoylated N-terminal-helix binding site at the other end of the "interswitch" beta hairpin. It is suggested that Sec7d groove closure upon docking of the two molecules may permit extraction of switch 1 from Arf1-GDP and prepare the complex for movement of the interswitch, which is central to the membrane-linked exchange activity. Large-scale collective movements in the Arf1-Sec7d complex appear to participate in the insertion of the Sec7d Glu finger into the GDP binding site to promote actual nucleotide release.
Uploads