Thermodynamic origin of cooperativity in actomyosin interactions: The coupling of short-range interactions with actin bending stiffness in an Ising-like model

We present Monte Carlo simulations for a molecular motor system found in virtually all eukaryotic cells, the acto-myosin motor system, composed of a group of organic macromolecules. Cell motors were mapped to an Ising-like model, where the interaction field is transmitted through a tropomyosin polymer chain. The presence of Ca(2+) induces tropomyosin to block or unblock binding sites of the myosin motor leading to its activation or deactivation. We used the Metropolis algorithm to find the transient and the equilibrium states of the acto-myosin system composed of solvent, actin, tropomyosin, troponin, Ca(2+), and myosin-S1 at a given temperature, including the spatial configuration of tropomyosin on the actin filament surface. Our model describes the short- and long-range cooperativity during actin-myosin binding which emerges from the bending stiffness of the tropomyosin complex. We found all transition rates between the states only using the interaction energy of the constituents. The agreement between our model and experimental data also supports the recent theory of flexible tropomyosin.

Catabolic ornithine carbamoyltransferase of Pseudomonas aeruginosa. Importance of the N-terminal region for dodecameric structure and homotropic carbamoylphosphate cooperativity.

Pseudomonas aeruginosa has an anabolic (ArgF) and a catabolic (ArcB) ornithine carbamoyltransferase (OTCase). Despite extensive sequence similarities, these enzymes function unidirectionally in vivo. In the dodecameric catabolic OTCase, homotropic cooperativity for carbamoylphosphate strongly depresses the anabolic reaction; the residue Glu1O5 and the C-terminus are known to be essential for this cooperativity. When Glu1O5 and nine C-terminal amino acids of the catabolic OTCase were introduced, by in vitro genetic manipulation, into the closely related, trimeric, anabolic (ArgF) OTCase of Escherichia coli, the enzyme displayed Michaelis-Menten kinetics and no cooperativity was observed. This indicates that additional amino acid residues are required to produce homotropic cooperativity and a dodecameric assembly. To localize these residues, we constructed several hybrid enzymes by fusing, in vivo or in vitro, the E. coli argF gene to the P. aeruginosa arcB gene. A hybrid enzyme consisting of 101 N-terminal ArgF amino acids fused to 233 C-terminal ArcB residues and the reciprocal ArcB-ArgF hybrid were both trimers with little or no cooperativity. Replacing the seven N-terminal residues of the ArcB enzyme by the corresponding six residues of E. coli ArgF enzyme produced a dodecameric enzyme which showed a reduced affinity for carbamoylphosphate and an increase in homotropic cooperativity. Thus, the N-terminal amino acids of catabolic OTCase are important for interaction with carbamoylphosphate, but do not alone determine dodecameric assembly. Hybrid enzymes consisting of either 26 or 42 N-terminal ArgF amino acids and the corresponding C-terminal ArcB residues were both trimeric, yet they retained some homotropic cooperativity. Within the N-terminal ArcB region, a replacement of motif 28-33 by the corresponding ArgF segment destabilized the dodecameric structure and the enzyme existed in trimeric and dodecameric states, indicating that this region is important for dodecameric assembly. These findings were interpreted in the light of the three-dimensional structure of catabolic OTCase, which allows predictions about trimer-trimer interactions. Dodecameric assembly appears to require at least three regions: the N- and C-termini (which are close to each other in a monomer), residues 28-33 and residues 147-154. Dodecameric structure correlates with high carbamoylphosphate cooperativity and thermal stability, but some trimeric hybrid enzymes retain cooperativity, and the dodecameric Glu1O5-->Ala mutant gives hyperbolic carbamoylphosphate saturation, indicating that dodecameric structure is neither necessary nor sufficient to ensure cooperativity.

Topological comparison of methods for predicting transcriptional cooperativity in yeast

Background: The cooperative interaction between transcription factors has a decisive role in the control of the fate of the eukaryotic cell. Computational approaches for characterizing cooperative transcription factors in yeast, however, are based on different rationales and provide a low overlap between their results. Because the wealth of information contained in protein interaction networks and regulatory networks has proven highly effective in elucidating functional relationships between proteins, we compared different sets of cooperative transcription factor pairs (predicted by four different computational methods) within the frame of those networks. Results: Our results show that the overlap between the sets of cooperative transcription factors predicted by the different methods is low yet significant. Cooperative transcription factors predicted by all methods are closer and more clustered in the protein interaction network than expected by chance. On the other hand, members of a cooperative transcription factor pair neither seemed to regulate each other nor shared similar regulatory inputs, although they do regulate similar groups of target genes. Conclusion: Despite the different definitions of transcriptional cooperativity and the different computational approaches used to characterize cooperativity between transcription factors, the analysis of their roles in the framework of the protein interaction network and the regulatory network indicates a common denominator for the predictions under study. The knowledge of the shared topological properties of cooperative transcription factor pairs in both networks can be useful not only for designing better prediction methods but also for better understanding the complexities of transcriptional control in eukaryotes.

Molecular dynamics study of orientational cooperativity in water

Recent experiments on liquid water show collective dipole orientation fluctuations dramatically slower than expected (with relaxation time >tation, the self-dipole randomization time tr, which is an upper limit on ta; we find that tr5ta. Third, to check if there are correlated domains of dipoles in water which have large relaxation times compared to the individual dipoles, we calculate the randomization time tbox of the site-dipole field, the net dipole moment formed by a set of molecules belonging to a box of edge Lbox. We find that the site-dipole randomization time tbox2.5ta for Lbox3 , i.e., it is shorter than the same quantity calculated for the self-dipole. Finally, we find that the orientational correlation length is short even at low T.

Functional and topological characterization of transcriptional cooperativity in yeast

This work was supported by grants from Spanish Ministry of Science andInnovation (MICINN) BIO2011-22568 & BIO2008-205.

Methionine-321 in the C-terminal alpha-helix of catabolic ornithine carbamoyltransferase from Pseudomonas aeruginosa is important for positive homotropic cooperativity.

Pseudomonas aeruginosa has a pair of distinct ornithine carbamoyltransferases. The anabolic ornithine carbamoyltransferase encoded by the argF gene catalyzes the formation of citrulline from ornithine and carbamoylphosphate. The catabolic ornithine carbamoyltransferase encoded by the arcB gene promotes the reverse reaction in vivo; although this enzyme can be assayed in vitro for citrulline synthesis, its unidirectionality in vivo is determined by its high concentration at half maximum velocity for carbamoylphosphate ([S]0.5) and high cooperativity toward this substrate. We have isolated mutant forms of catabolic ornithine carbamoyltransferase catalyzing the anabolic reaction in vivo. The corresponding arcB mutant alleles on a multicopy plasmid specifically suppressed an argF mutation of P. aeruginosa. Two new mutant enzymes were obtained. When methionine 321 was replaced by isoleucine, the mutant enzyme showed loss of homotropic cooperativity at physiological carbamoylphosphate concentrations. Substitution of glutamate 105 by lysine resulted in a partial loss of the sigmoidal response to increasing carbamoylphosphate concentrations. However, both mutant enzymes were still sensitive to the allosteric activator AMP and to the inhibitor spermidine. These results indicate that at least two residues of catabolic ornithine carbamoyltransferase are critically involved in positive carbamoylphosphate cooperativity: glutamate 105 (previously known to be important) and methionine 321. Mutational changes in either amino acid will affect the geometry of helix H2, which contains several residues required for carbamoylphosphate binding.

Investigation of cooperativity in the binding of ligands to the D-2 dopamine receptor

The D 2 dopamine receptor exists as dimers or as higher-order oligomers, as determined from data from physical experiments. In this study, we sought evidence that this oligomerization leads to cooperativity by examining the binding of three radioligands ([H-3] nemonapride, [H-3] raclopride, and [H-3] spiperone) to D 2 dopamine receptors expressed in membranes of Sf9 cells. In saturation binding experiments, the three radioligands exhibited different B-max values, and the B-max values could be altered by the addition of sodium ions to assays. Despite labeling different numbers of sites, the different ligands were able to achieve full inhibition in competition experiments. Some ligand pairs also exhibited complex inhibition curves in these experiments. In radioligand dissociation experiments, the rate of dissociation of [H-3] nemonapride or [H-3] spiperone depended on the sodium ion concentration but was independent of the competing ligand. Although some of the data in this study are consistent with the behavior of a cooperative oligomeric receptor, not all of the data are in agreement with this model. It may, therefore, be necessary to consider more complex models for the behavior of this receptor.

Investigation of cooperativity in the binding of ligands to the D-2 dopamine receptor

The D 2 dopamine receptor exists as dimers or as higher-order oligomers, as determined from data from physical experiments. In this study, we sought evidence that this oligomerization leads to cooperativity by examining the binding of three radioligands ([H-3] nemonapride, [H-3] raclopride, and [H-3] spiperone) to D 2 dopamine receptors expressed in membranes of Sf9 cells. In saturation binding experiments, the three radioligands exhibited different B-max values, and the B-max values could be altered by the addition of sodium ions to assays. Despite labeling different numbers of sites, the different ligands were able to achieve full inhibition in competition experiments. Some ligand pairs also exhibited complex inhibition curves in these experiments. In radioligand dissociation experiments, the rate of dissociation of [H-3] nemonapride or [H-3] spiperone depended on the sodium ion concentration but was independent of the competing ligand. Although some of the data in this study are consistent with the behavior of a cooperative oligomeric receptor, not all of the data are in agreement with this model. It may, therefore, be necessary to consider more complex models for the behavior of this receptor.

Insights into Phosphate Cooperativity and Influence of Substrate Modifications on Binding and Catalysis of Hexameric Purine Nucleoside Phosphorylases

The hexameric purine nucleoside phosphorylase from Bacillus subtilis (BsPNP233) displays great potential to produce nucleoside analogues in industry and can be exploited in the development of new anti-tumor gene therapies. In order to provide structural basis for enzyme and substrates rational optimization, aiming at those applications, the present work shows a thorough and detailed structural description of the binding mode of substrates and nucleoside analogues to the active site of the hexameric BsPNP233. Here we report the crystal structure of BsPNP233 in the apo form and in complex with 11 ligands, including clinically relevant compounds. The crystal structure of six ligands (adenine, 2'deoxyguanosine, aciclovir, ganciclovir, 8-bromoguanosine, 6-chloroguanosine) in complex with a hexameric PNP are presented for the first time. Our data showed that free bases adopt alternative conformations in the BsPNP233 active site and indicated that binding of the co-substrate (2'deoxy) ribose 1-phosphate might contribute for stabilizing the bases in a favorable orientation for catalysis. The BsPNP233-adenosine complex revealed that a hydrogen bond between the 5' hydroxyl group of adenosine and Arg(43*) side chain contributes for the ribosyl radical to adopt an unusual C3'-endo conformation. The structures with 6-chloroguanosine and 8-bromoguanosine pointed out that the Cl-6 and Br-8 substrate modifications seem to be detrimental for catalysis and can be explored in the design of inhibitors for hexameric PNPs from pathogens. Our data also corroborated the competitive inhibition mechanism of hexameric PNPs by tubercidin and suggested that the acyclic nucleoside ganciclovir is a better inhibitor for hexameric PNPs than aciclovir. Furthermore, comparative structural analyses indicated that the replacement of Ser(90) by a threonine in the B. cereus hexameric adenosine phosphorylase (Thr(91)) is responsible for the lack of negative cooperativity of phosphate binding in this enzyme.

Order-Disorder Transitions Govern Kinetic Cooperativity and Allostery of Monomeric Human Glucokinase

Glucokinase (GCK) catalyzes the rate-limiting step of glucose catabolism in the pancreas, where it functions as the body's principal glucose sensor. GCK dysfunction leads to several potentially fatal diseases including maturity-onset diabetes of the young type II (MODY-II) and persistent hypoglycemic hyperinsulinemia of infancy (PHHI). GCK maintains glucose homeostasis by displaying a sigmoidal kinetic response to increasing blood glucose levels. This positive cooperativity is unique because the enzyme functions exclusively as a monomer and possesses only a single glucose binding site. Despite nearly a half century of research, the mechanistic basis for GCK's homotropic allostery remains unresolved. Here we explain GCK cooperativity in terms of large-scale, glucose-mediated disorder-order transitions using 17 isotopically labeled isoleucine methyl groups and three tryptophan side chains as sensitive nuclear magnetic resonance (NMR) probes. We find that the small domain of unliganded GCK is intrinsically disordered and samples a broad conformational ensemble. We also demonstrate that small-molecule diabetes therapeutic agents and hyperinsulinemia-associated GCK mutations share a strikingly similar activation mechanism, characterized by a population shift toward a more narrow, well-ordered ensemble resembling the glucose-bound conformation. Our results support a model in which GCK generates its cooperative kinetic response at low glucose concentrations by using a millisecond disorder-order cycle of the small domain as a "time-delay loop," which is bypassed at high glucose concentrations, providing a unique mechanism to allosterically regulate the activity of human GCK under physiological conditions.

Propagation of melting cooperativity along the phosphodiester backbone of DNA

The role of the DNA phosphodiester backbone in the transfer of melting cooperativity between two helical domains was experimentally addressed with a helix-bulge-helix DNA model, in which the bulge consisted of a varying number of either conformationally flexible propanediol or conformationally constrained bicyclic anucleosidic phosphodiester backbone units. We found that structural communication between two double helical domains is transferred along the DNA backbone over the equivalent of ca. 12-20 backbone units, depending on whether there is a symmetric or asymmetric distribution of the anucleosidic units on both strands. We observed that extension of anucleosidic units on one strand only suffices to disrupt cooperativity transfer in a similar way as if extension occurs on both strands, indicating that the length of the longest anucleosidic inset determines cooperativity transfer. Furthermore, conformational rigidity of the sugar unit increases the distance of coopertivity transfer along the phosphodiester backbone. This is especially the case when the units are asymmetrically distributed in both strands

Analysis of human cytochrome P450 3A4 cooperativity: Construction and characterization of a site-directed mutant that displays hyperbolic steroid hydroxylation kinetics

Cytochrome P450 3A4 is generally considered to be the most important human drug-metabolizing enzyme and is known to catalyze the oxidation of a number of substrates in a cooperative manner. An allosteric mechanism is usually invoked to explain the cooperativity. Based on a structure–activity study from another laboratory using various effector–substrate combinations and on our own studies using site-directed mutagenesis and computer modeling of P450 3A4, the most likely location of effector binding is in the active site along with the substrate. Our study was designed to test this hypothesis by replacing residues Leu-211 and Asp-214 with the larger Phe and Glu, respectively. These residues were predicted to constitute a portion of the effector binding site, and the substitutions were designed to mimic the action of the effector by reducing the size of the active site. The L211F/D214E double mutant displayed an increased rate of testosterone and progesterone 6β-hydroxylation at low substrate concentrations and a decreased level of heterotropic stimulation elicited by α-naphthoflavone. Kinetic analyses of the double mutant revealed the absence of homotropic cooperativity with either steroid substrate. At low substrate concentrations the steroid 6β-hydroxylase activity of the wild-type enzyme was stimulated by a second steroid, whereas L211F/D214E displayed simple substrate inhibition. To analyze L211F/D214E at a more mechanistic level, spectral binding studies were carried out. Testosterone binding by the wild-type enzyme displayed homotropic cooperativity, whereas substrate binding by L211F/D214E displayed hyperbolic behavior.

The unique hetero-oligomeric nature of the subunits in the catalytic cooperativity of the yeast Cct chaperonin complex

The structural and functional organization of the Cct complex was addressed by genetic analyses of subunit interactions and catalytic cooperativity among five of the eight different essential subunits, Cct1p–Cct8p, in the yeast Saccharomyces cerevisiae. The cct1–1, cct2–3, and cct3–1 alleles, containing mutations at the conserved putative ATP-binding motif, GDGTT, are cold-sensitive, whereas single and multiple replacements of the corresponding motif in Cct6p are well tolerated by the cell. We demonstrated herein that cct6–3 (L19S), but not the parolog cct1–5 (R26I), specifically suppresses the cct1–1, cct2–3, and cct3–1 alleles, and that this suppression can be modulated by mutations in a putative phosphorylation motif, RXS, and the putative ATP-binding pocket of Cct6p. Our results suggest that the Cct ring is comprised of a single hetero-oligomer containing eight subunits of differential functional hierarchy, in which catalytic cooperativity of ATP-binding/hydrolysis takes place in a sequential manner different from the concerted cooperativity proposed for GroEL.