Mutational analysis of allosteric activation and inhibition of glucokinase.

GK (glucokinase) is activated by glucose binding to its substrate site, is inhibited by GKRP (GK regulatory protein) and stimulated by GKAs (GK activator drugs). To explore further the mechanisms of these processes we studied pure recombinant human GK (normal enzyme and a selection of 31 mutants) using steady-state kinetics of the enzyme and TF (tryptophan fluorescence). TF studies of the normal binary GK-glucose complex corroborate recent crystallography studies showing that it exists in a closed conformation greatly different from the open conformation of the ligand-free structure, but indistinguishable from the ternary GK-glucose-GKA complex. GKAs did activate and GKRP did inhibit normal GK, whereas its TF was doubled by glucose saturation. However, the enzyme kinetics, GKRP inhibition, TF enhancement by glucose and responsiveness to GKA of the selected mutants varied greatly. Two predominant response patterns were identified accounting for nearly all mutants: (i) GK mutants with a normal or close to normal response to GKA, normally low basal TF (indicating an open conformation), some variability of kinetic parameters (k(cat), glucose S(0.5), h and ATP K(m)), but usually strong GKRP inhibition (13/31); and (ii) GK mutants that are refractory to GKAs, exhibit relatively high basal TF (indicating structural compaction and partial closure), usually show strongly enhanced catalytic activity primarily due to lowering of the glucose S(0.5), but with reduced or no GKRP inhibition in most cases (14/31). These results and those of previous studies are best explained by envisioning a common allosteric regulator region with spatially non-overlapping GKRP- and GKA-binding sites.

Hypermethylation of the ABCB1 downstream gene promoter accompanies ABCB1 gene amplification and increased expression in docetaxel-resistant MCF-7 breast tumor cells.

Drug transporters have been implicated in resistance of solid and non-solid tumors to a variety of chemotherapeutic agents. Higher expression of the ABCB1 drug transporter is often observed in drug-resistant tumor cells, although the precise mechanism remains unclear. During selection of MCF-7 cells for survival in increasing concentrations of docetaxel (MCF-7TXT cells), we observed in this study a temporal correlation between the acquisition of docetaxel resistance at selection dose 9 and the increased expression of ABCB1. Both the magnitude of docetaxel resistance and the level of ABCB1 expression then rose as the selection dose was further elevated. We also observed through bisulfite sequencing experiments that the ABCB1 downstream promoter became increasingly methylated following the acquisition of drug resistance (selection doses 10-12). Transcription was solely attributed to the upstream ABCB1 promoter within MCF-7TXT cells at the highest selection dose suggesting that hypermethylation caused a shift in promoter usage. The hypermethylation was also accompanied by regional amplification of chromosome 7 containing the ABCB1 gene and its neighbor ABCB4 but not DBF-4. The amplification of the ABCB1 gene correlated positively both with the hypermethylation of the ABCB1 downstream promoter (r=0.90) and the increased expression of ABCB1 (r=0.78). Moreover demethylation of the ABCB1 downstream promoter induced by 5-aza-2A'deoxycytidine treatment decreased the expression of ABCB1 mRNA in MCF-7TXT cells. Taken together, our findings suggest that the increased expression of ABCB1 upon acquisition of docetaxel resistance in breast tumor cells can be multifactorial, involving both epigenetic changes in promoter usage and regional chromosome amplification.

Role of drug transporters and drug accumulation in the temporal acquisition of drug resistance.

BACKGROUND: Anthracyclines and taxanes are commonly used in the treatment of breast cancer. However, tumor resistance to these drugs often develops, possibly due to overexpression of drug transporters. It remains unclear whether drug resistance in vitro occurs at clinically relevant doses of chemotherapy drugs and whether both the onset and magnitude of drug resistance can be temporally and causally correlated with the enhanced expression and activity of specific drug transporters. To address these issues, MCF-7 cells were selected for survival in increasing concentrations of doxorubicin (MCF-7DOX-2), epirubicin (MCF-7EPI), paclitaxel (MCF-7TAX-2), or docetaxel (MCF-7TXT). During selection cells were assessed for drug sensitivity, drug uptake, and the expression of various drug transporters. RESULTS: In all cases, resistance was only achieved when selection reached a specific threshold dose, which was well within the clinical range. A reduction in drug uptake was temporally correlated with the acquisition of drug resistance for all cell lines, but further increases in drug resistance at doses above threshold were unrelated to changes in cellular drug uptake. Elevated expression of one or more drug transporters was seen at or above the threshold dose, but the identity, number, and temporal pattern of drug transporter induction varied with the drug used as selection agent. The pan drug transporter inhibitor cyclosporin A was able to partially or completely restore drug accumulation in the drug-resistant cell lines, but had only partial to no effect on drug sensitivity. The inability of cyclosporin A to restore drug sensitivity suggests the presence of additional mechanisms of drug resistance. CONCLUSION: This study indicates that drug resistance is achieved in breast tumour cells only upon exposure to concentrations of drug at or above a specific selection dose. While changes in drug accumulation and the expression of drug transporters does occur at the threshold dose, the magnitude of resistance cannot be attributed solely to changes in drug accumulation or the activity of drug transporters. The identities of these additional drug-transporter-independent mechanisms are discussed, including their likely clinical relevance.

Protein dynamics explain the allosteric behaviors of hemoglobin.

Bohr, Hasselbalch, and Krogh discovered homotropic and heterotropic allosteric behaviors of hemoglobin (Hb) in 1903/1904. A chronological description since then of selected principal models of the allosteric mechanism of Hb, such as the Adair scheme, the MWC two-state concerted model, the KNF induced-fit sequential model, the Perutz stereochemical model, the tertiary two-state model, and the global allostery model (an expanded MWC models), is concisely presented, followed by analysis and discussion of their limitations and deficiencies. The determination of X-ray crystallographic structures of deoxy- and ligated-Hb and the structure-based stereochemical model by Perutz are an epoch-making event in this history. However, his assignment of low-affinity deoxy- and high-affinity oxy-quaternary structures of Hb to the T- and R-states, respectively, though apparently reasonable, and as well as his hypothesis that the T-/R-quaternary structural transition regulates the oxygen-affinity, have created confusions and side-tracked studies of Hb on the structure-function relationship. The differences in static molecular structures of Hb between T(deoxy)- and R(oxy)-quaternary states reported in detail by Perutz and others are ligation-linked structural changes, but not related to the control/regulation of the oxygen-affinity. The oxygen-affinity (K(T) and K(R)) of Hb has been shown to be regulated by the heterotropic effector-linked tertiary structural changes without involving the T/R-quaternary changes. However, a recent high-resolution crystallographic analysis of Hb with different oxygen-affinities shows that static molecular structures of Hb determined by crystallography can neither identify the nature of the T(low-affinity) functional state nor decipher the mechanism by which Hb stores free energy in the T(low-affinity) functional state. Molecular dynamics simulations show that fluctuations of helices of oxy-Hb are increased upon de-oxygenation and/or binding 2,3-biphosphoglycerate. These are known to lower the oxygen-affinity of Hb. It is proposed that the coordination mode of the heme Fe with proximal and distal His is modulated by these helical fluctuations, resulting in the modulation of the oxygen-affinity of Hb. Therefore, it is proposed that the oxygen-affinity of Hb is regulated by pentanary (the 5th-order time-dependent or dynamic) tertiary structural changes rather than the T-/R-quaternary structural transitions in Hb. Homotropic and heterotropic allosteric effects of Hb are oxygen- and effector-linked, conformational entropy-driven entropy-enthalpy compensation phenomena and not much to do with static structural changes. The dynamic allostery model, which integrates these observations, provides the structural basis for the global allostery model (an expanded MWC model).

Interaction of antagonists with calmodulin: insights from molecular dynamics simulations.

We report results of 12 ns, all-atom molecular dynamics simulation (MDS) and Poisson-Boltzmann free energy calculations (PBFE) on calmodulin (CaM) bound to two molecules of trifluoperazine (TFP) and of N-(3,3, diphenylpropyl)- N'-[1- R-(3,4-bis-butoxyphenyl)-ethyl]-propylenediamine (DPD). X-ray data show very similar structures for the two complexes, yet the antagonists significantly differ with respect to their CaM binding affinities, the neutral DPD is much more potent. The goal of the study was to unravel the reason why TFP is less potent although its positive charge should facilitate binding. The electrostatic energy terms in CHARMM and binding free energy terms of the PBFE approach showed TFP a better antagonist, while inspection of hydrophobic contacts supports DPD binding. Detailed inspection of the amino acid contributions of PBFE calculations unravel that steric reasons oppose the favorable binding of TFP. Structural conditions are given for a successful drug design strategy, which may benefit also from charge-charge interactions.

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.

Common dynamics of globin family proteins.

The recently discovered new members of the globin family, neurogobin and cytoglobin, are the object of sustained structural and functional studies aimed at understanding their physiological role and elucidating the impact of their bis-his heme hexacoordination. However, no studies have yet considered the dynamics of this protein family, an essential link between structure and function. In this communication, we present normal mode analysis results for neuroglobin, cytoglobin, hemoglobin and myoglobin to provide exploratory insights into globin characteristic motions. Our results show a clear correlation in the protein dynamics of this family. All four globins exhibit a high degree of correlated displacements involving residues in the C, E and F helices and link regions. They suggest that these motions play an important role in the reversible oxygen binding function of these proteins. Further, our results may help rationalize some functional features of the 6c-globins in that they alone exhibit correlated displacements of the G-helix region.

Static normal coordinate deformations of the heme group in mutants of ferrocytochrome c from Saccharomyces cerevisiae probed by resonance Raman spectroscopy.

The function of heme proteins is, to a significant extent, influenced by the ligand field probed by the heme iron, which itself can be affected by deformations of the heme macrocycle. The exploration of this field is difficult because the heme structure obtained from X-ray crystallography is not resolved enough to unambiguously identify structural changes on the scale of 10(-2) A. However, asymmetric deformations in this order of magnitude affect the depolarization ratio of the resonance Raman lines assignable to normal vibrations of the heme group. We have measured the dispersion of the depolarization ratios of four structure sensitive Raman bands (i.e., nu4, nu11, nu21, and nu28) in yeast iso-1-ferrocytochrome c and its mutants N52V, Y67F, and N52VY67F with B- and Q-band excitation. The DPR dispersion of all bands indicates the presence of asymmetric in-plane and out-of-plane deformations. The replacement of the polar tyrosine residue at position 67 by phenylalanine significantly increases the triclinic B2g deformation, which involves a distortion of the pyrrole symmetry. We relate this deformation to changes of the electronic structure of pyrrole A, which modulates the interaction between its propionate substituents and the protein environment. This specific heme deformation is eliminated in the double mutant N52VY67F. The additional substitution of N52 by valine induces a tetragonal B1g deformation which involves asymmetric changes of the Fe-N distances and increases the rhombicity of the ligand field probed by the heme iron. This heme deformation might be caused by the elimination of the water-protein hydrogen-bonding network in the heme cavity. The single mutation N52V does not significantly perturb the heme symmetry, but a small B1g deformation is consistent with our data and the heme structure obtained from a 1 ns molecular dynamics simulation of the protein.

The influence of interdomain interactions on the intradomain motions in yeast phosphoglycerate kinase: a molecular dynamics study.

A 3-ns molecular dynamics simulation in explicit solvent was performed to examine the inter- and intradomain motions of the two-domain enzyme yeast phosphoglycerate kinase without the presence of substrates. To elucidate contributions from individual domains, simulations were carried out on the complete enzyme as well as on each isolated domain. The enzyme is known to undergo a hinge-bending type of motion as it cycles from an open to a closed conformation to allow the phosphoryl transfer occur. Analysis of the correlation of atomic movements during the simulations confirms hinge bending in the nanosecond timescale: the two domains of the complete enzyme exhibit rigid body motions anticorrelated with respect to each other. The correlation of the intradomain motions of both domains converges, yielding a distinct correlation map in the enzyme. In the isolated domain simulations-in which interdomain interactions cannot occur-the correlation of domain motions no longer converges and shows a very small correlation during the same simulation time. This result points to the importance of interdomain contacts in the overall dynamics of the protein. The secondary structure elements responsible for interdomain contacts are also discussed.

Functionally relevant electric-field induced perturbations of the prosthetic group of yeast ferrocytochrome c mutants obtained from a vibronic analysis of low-temperature absorption spectra.

We have measured the low temperature (T = 20 K) absorption spectra of the N52A, N52V, N52I, Y67F, and N52AY67F mutants of ferrous Saccharomyces cerevisiae (baker's yeast) cytochrome c. All the bands in the Q0- and Q(v)-band region are split, and the intensity distributions among the split bands are highly asymmetric. The spectra were analyzed by a decomposition into Voigtian profiles. The spectral parameters thus obtained were further analyzed in terms of the vibronic coupling model of Schweitzer-Stenner and Bigman (Schweitzer-Stenner, R.; Bigman, D. J. Phys. Chem. B 2001, 7064-7073) to identify parameters related to electronic and vibronic perturbations of the heme macrocycle. We report that the electronic perturbation is of B(1g) symmetry and reflects the heterogeneity of the electric field at the heme, that is, the difference between the gradients along the perpendicular N-Fe-N axis of the heme core. We found that all the investigated mutations substantially increase this electronic perturbation, so that the spectral properties become similar to those of horse heart cytochrome c. Moreover, the electronic perturbation was found to correlate nonlinearly with the enthalpy changes associated with the reduction of the heme iron. Group theoretical arguments are invoked to propose a simple model which explains how a perturbation of the obtained symmetry can stabilize the reduced state of the heme iron. Finally, vibronic coupling parameters obtained from the analysis of the Q(v)-band region suggest that the investigated mutations decrease the nonplanar deformations of the heme group. This finding was reproduced by a normal mode structural decomposition (NSD) analysis of the N52V and N52VY67F heme conformations obtained from a 1 ns molecular dynamics simulation. We argue that the reduced nonplanarity contributes to the stabilization of the reduced state.

Normal mode analysis of the horseradish peroxidase collective motions: correlation with spectroscopically observed heme distortions.

Horseradish peroxidase C is a class III peroxidase whose structure is stabilized by the presence of two endogenous calcium atoms. Calcium removal has been shown to decrease the enzymatic activity of the enzyme and significantly affect the spectroscopically detectable properties of the heme, such as the spin state of the iron, heme normal modes, and distortions from planarity. In this work, we report on normal mode analysis (NMA) performed on models subjected to 2 ns of molecular dynamics simulations to describe the effect of calcium removal on protein collective motions and to investigate the correlation between active site (heme) and protein matrix fluctuations. We show that in the native peroxidase model, heme fluctuations are correlated to matrix fluctuations while they are not in the calcium-depleted model.

The importance of vibronic perturbations in ferrocytochrome c spectra: a reevaluation of spectral properties based on low-temperature optical absorption, resonance Raman, and molecular-dynamics simulations.

We have measured and analyzed the low-temperature (T=10 K) absorption spectrum of reduced horse heart and yeast cytochrome c. Both spectra show split and asymmetric Q(0) and Q(upsilon) bands. The spectra were first decomposed into the individual split vibronic sidebands assignable to B(1g) (nu15) and A(2g) (nu19, nu21, and nu22) Herzberg-Teller active modes due to their strong intensity in resonance Raman spectra acquired with Q(0) and Q(upsilon) excitations. The measured band splittings and asymmetries cannot be rationalized solely in terms of electronic perturbations of the heme macrocycle. On the contrary, they clearly point to the importance of considering not only electronic perturbations but vibronic perturbations as well. The former are most likely due to the heterogeneity of the electric field produced by charged side chains in the protein environment, whereas the latter reflect a perturbation potential due to multiple heme-protein interactions, which deform the heme structure in the ground and excited states. Additional information about vibronic perturbations and the associated ground-state deformations are inferred from the depolarization ratios of resonance Raman bands. The results of our analysis indicate that the heme group in yeast cytochrome c is more nonplanar and more distorted along a B(2g) coordinate than in horse heart cytochrome c. This conclusion is supported by normal structural decomposition calculations performed on the heme extracted from molecular-dynamic simulations of the two investigated proteins. Interestingly, the latter are somewhat different from the respective deformations obtained from the x-ray structures.

R-state hemoglobin bound to heterotropic effectors: models of the DPG, IHP and RSR13 binding sites.

We performed a docking study followed by a 500-ps molecular dynamics simulation of R-state human adult hemoglobin (HbA) complexed to different heterotropic effectors [2,3-diphosphoglycerate (DPG), inositol hexaphosphate (IHP), and 2-[4-[(3,5-dichlorophenylcarbamoyl)-]methyl]-phenoxy]-2-methylpropionic acid (RSR13)) to propose a molecular basis for recently reported interactions of effectors with oxygenated hemoglobin. The simulations were carried out with counterions and explicit solvation. As reported for T-state HbA, the effector binding sites are also located in the central cavity of the R-state and differ depending on effector anionic character. DPG and IHP bind between the alpha-subunits and the RSR13 site spans the alpha1-, alpha2- and beta2-subunits. The generated models provide the first report of the molecular details of R-state HbA bound to heterotropic effectors.

The charge transfer band in horseradish peroxidase correlates with heme in-plane distortions induced by calcium removal.

Horseradish peroxidase C (HRPC) is a class III peroxidase whose structure is stabilized by the presence of two endogenous calcium atoms. Calcium removal has been shown to decrease the enzymatic activity of the enzyme. The spin state of the iron, a mixture of high spin (HS) and mixed quantum spin state (QS) consisting of intermediate spin (IS) 3/2 + (HS) 5/2, is also significantly affected by calcium removal, going from a predominant QS component to a predominant HS component upon removal of one calcium. Removal of both calcium ions, however, results in the appearance of a significant LS contribution, easily monitored in the charge transfer (CT) band region by low-T absorption. Normal structural decomposition (NSD) calculations of the in-plane (ip) modes of the heme extracted from HRPC native and Ca(2+)-depleted models show that removal of the proximal calcium is associated with perturbed E(u) and increased A(1g) ip distortions of the heme. The effect of complete or distal calcium removal on the heme also results in increased A(1g) ip distortions, but in significantly decreased E(u) distortions. The overall effect is to decrease the nonplanarity of the heme: the total ip distortion of the native HRPC heme is 0.200 and 0.134 A for the Ca(2+)-depleted species. Our NSD results corroborate the role proposed for the protein matrix, namely to fine-tune the active site by inducing subtle changes in heme planarity and spin state of the iron.

Resonance Raman spectroscopy study of change of iron spin state in horseradish peroxidase C induced by removal of calcium.

Resonance Raman spectroscopy is used to probe the effect of calcium depletion on the heme group of horseradish peroxidase C at pH 8. Polarized Raman spectra are recorded with an argon ion laser at eight different wavelengths to provide a sound database for a reliable spectral decomposition. Upon calcium depletion, the spectrum is indicative of a predominantly pentacoordinated high spin state of the heme iron coexisting with small fractions of hexacoordinated high and low spin states. The dominant quantum mixed spin state of native ferric horseradish peroxidase, which is characteristic for class III peroxidases, is not detectable in the spectrum of the enzyme with partial distal Ca(2+) depletion. The quenching of the quantum mixed spin state and the predominance of the pentacoordinated high spin state are likely to arise from distortions induced by distal calcium depletion, which translates into a weaker Fe-N(epsilon)(His) bond and a more tilted imidazole. A correlation is proposed between the lower enzyme activity and the elimination of the pentacoordinated quantum mixed state upon Ca(2+) depletion.

The endogenous calcium ions of horseradish peroxidase C are required to maintain the functional nonplanarity of the heme.

Horseradish peroxidase C (HRPC) binds 2 mol calcium per mol of enzyme with binding sites located distal and proximal to the heme group. The effect of calcium depletion on the conformation of the heme was investigated by combining polarized resonance Raman dispersion spectroscopy with normal coordinate structural decomposition analysis of the hemes extracted from models of Ca(2+)-bound and Ca(2+)-depleted HRPC generated and equilibrated using molecular dynamics simulations. Results show that calcium removal causes reorientation of heme pocket residues. We propose that these rearrangements significantly affect both the in-plane and out-of-plane deformations of the heme. Analysis of the experimental depolarization ratios are clearly consistent with increased B(1g)- and B(2g)-type distortions in the Ca(2+)-depleted species while the normal coordinate structural decomposition results are indicative of increased planarity for the heme of Ca(2+)-depleted HRPC and of significant changes in the relative contributions of three of the six lowest frequency deformations. Most noteworthy is the decrease of the strong saddling deformation that is typical of all peroxidases, and an increase in ruffling. Our results confirm previous work proposing that calcium is required to maintain the structural integrity of the heme in that we show that the preferred geometry for catalysis is lost upon calcium depletion.

Normal coordinate structural decomposition of the heme distortions of hemoglobin in various quaternary states and bound to allosteric effectors.

The distortions of the alpha1, alpha2, beta1, and beta2 hemes of human hemoglobin (HbA) in various quaternary states and as affected by the presence of allosteric effectors was investigated by subjecting CHARMM energy-minimized models to normal coordinate structural decomposition (NSD) analysis. NSD was applied to the individual hemes extracted from the R, T, and R2-state models of HbA and to HbA bound to DPG and to IHP. Overall, NSD results are indicative of characteristic distortions, not only for the hemes of the different HbA quaternary states, but also for the hemes of the HbA models bound to allosteric effectors. Comparing the distortions of the inequivalent alpha and beta hemes in T-state HbA, we show good correlation between NSD and the experimentally observed low-frequency nu52 (Eg) and gamma7 (A2u) modes reported in the literature for alpha and beta HbA hemes while noting substantial differences between these types for B2u and B1u distortions. For the R2 hemes, NSD yields heme distortions that are more comparable to those of the R-state, especially in magnitude. However, the R2 hemes do not exhibit inequivalence of alpha and beta heme distortions, a result that may contribute to an understanding of the functional importance of this state. Relative to T-state heme distortions, NSD results on the effector-bound hemes show that tertiary changes induced in T-state HbA as a result of binding DPG and IHP drastically affect heme distortions. In the alpha hemes extracted from the HbA-DPG model, most noteworthy are the increased wav(x) and wav(y) distortions and enhancement of ruf and dom deformations. In the beta hemes, the wav(y) is the most affected distortion with increase in sad. The NSD results are also different for the hemes of the HbA-IHP model, in that the beta sad and ruf deformations are more enhanced with increase of doming in the alpha hemes. Our results describe the effect of the subtle protein-induced changes on the nonplanarity of the HbA hemes that may play a role in the regulation of their oxygen affinity.