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.

Magnesium-dependent folding of self-splicing RNA: Exploring the link between cooperativity, thermodynamics, and kinetics

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

Membrane-bound state of the colicin E1 channel domain as an extended two-dimensional helical array

Atomic level structures have been determined for the soluble forms of several colicins and toxins, but the structural changes that occur after membrane binding have not been well characterized. Changes occurring in the transition from the soluble to membrane-bound state of the C-terminal 190-residue channel polypeptide of colicin E1 (P190) bound to anionic membranes are described. In the membrane-bound state, the α-helical content increases from 60–64% to 80–90%, with a concomitant increase in the average length of the helical segments from 12 to 16 or 17 residues, close to the length required to span the membrane bilayer in the open channel state. The average distance between helical segments is increased and interhelix interactions are weakened, as shown by a major loss of tertiary structure interactions, decreased efficiency of fluorescence resonance energy transfer from an energy donor on helix V of P190 to an acceptor on helix IX, and decreased resonance energy transfer at higher temperatures, not observed in soluble P190, implying freedom of motion of helical segments. Weaker interactions are also shown by a calorimetric thermal transition of low cooperativity, and the extended nature of the helical array is shown by a 3- to 4-fold increase in the average area subtended per molecule to 4,200 Å2 on the membrane surface. The latter, with analysis of the heat capacity changes, implies the absence of a developed hydrophobic core in the membrane-bound P190. The membrane interfacial layer thus serves to promote formation of a highly helical extended two-dimensional flexible net. The properties of the membrane-bound state of the colicin channel domain (i.e., hydrophobic anchor, lengthened and loosely coupled α-helices, and close association with the membrane interfacial layer) are plausible structural features for the state that is a prerequisite for voltage gating, formation of transmembrane helices, and channel opening.

Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

The inv(16) encodes an acute myeloid leukemia 1 transcriptional corepressor

The inv(16) is one of the most frequent chromosomal translocations associated with acute myeloid leukemia (AML). The inv(16) fusion protein acts by dominantly interfering with AML-1/core binding factor β-dependent transcriptional regulation. Here we demonstrate that the inv(16) fusion protein cooperates with AML-1B to repress transcription. This cooperativity requires the ability of the translocation fusion protein to bind to AML-1B. Mutational analysis and cell fractionation experiments indicated that the inv(16) fusion protein acts in the nucleus and that repression occurs when the complex is bound to DNA. We also found that the inv(16) fusion protein binds to AML-1B when it is associated with the mSin3A corepressor. An AML-1B mutant that fails to bind mSin3A was impaired in cooperative repression, suggesting that the inv(16) fusion protein acts through mSin3 and possibly other corepressors. Finally, we demonstrate that the C-terminal portion of the inv(16) fusion protein contains a repression domain, suggesting a molecular mechanism for AML-1-mediated repression.

In vitro properties of the first ORF protein from mouse LINE-1 support its role in ribonucleoprotein particle formation during retrotransposition

LINEs are transposable elements, widely distributed among eukaryotes, that move via reverse transcription of an RNA intermediate. Mammalian LINEs have two ORFs (ORF1 and ORF2). The proteins encoded by these ORFs play important roles in the retrotransposition process. Although the predicted amino acid sequence of ORF1 is not closely related to any known proteins, it is highly basic; thus, it has long been hypothesized that ORF1 protein functions to bind LINE-1 (L1) RNA during retrotransposition. Cofractionation of ORF1 protein and L1 RNA in extracts from both mouse and human embryonal carcinoma cells indicated that ORF1 protein binds L1 RNA, forming a ribonucleoprotein particle. Based on UV crosslinking and electrophoretic mobility-shift assays using purified components, we demonstrate here that the ORF1 protein encoded by mouse L1 binds nucleic acids with a strong preference for RNA and other single-stranded nucleic acids. Furthermore, multiple copies of ORF1 protein appear to bind single-stranded nucleic acid in a manner suggesting positive cooperativity; such binding characteristics are likely to be facilitated by the protein–protein interactions detected among molecules of ORF1 polypeptide by coimmunoprecipitation. These observations are consistent with the formation of ribonucleoprotein particles containing L1 RNA and ORF1 protein and provide additional evidence for the role of ORF1 protein during retrotransposition of L1.

A model of excitation and adaptation in bacterial chemotaxis

Bacterial chemotaxis is widely studied because of its accessibility and because it incorporates processes that are important in a number of sensory systems: signal transduction, excitation, adaptation, and a change in behavior, all in response to stimuli. Quantitative data on the change in behavior are available for this system, and the major biochemical steps in the signal transduction/processing pathway have been identified. We have incorporated recent biochemical data into a mathematical model that can reproduce many of the major features of the intracellular response, including the change in the level of chemotactic proteins to step and ramp stimuli such as those used in experimental protocols. The interaction of the chemotactic proteins with the motor is not modeled, but we can estimate the degree of cooperativity needed to produce the observed gain under the assumption that the chemotactic proteins interact directly with the motor proteins.

Systematic derivation of partition functions for ligand binding to two-dimensional lattices

The Ising problem consists in finding the analytical solution of the partition function of a lattice once the interaction geometry among its elements is specified. No general analytical solution is available for this problem, except for the one-dimensional case. Using site-specific thermodynamics, it is shown that the partition function for ligand binding to a two-dimensional lattice can be obtained from those of one-dimensional lattices with known solution. The complexity of the lattice is reduced recursively by application of a contact transformation that involves a relatively small number of steps. The transformation implemented in a computer code solves the partition function of the lattice by operating on the connectivity matrix of the graph associated with it. This provides a powerful new approach to the Ising problem, and enables a systematic analysis of two-dimensional lattices that model many biologically relevant phenomena. Application of this approach to finite two-dimensional lattices with positive cooperativity indicates that the binding capacity per site diverges as Na (N = number of sites in the lattice) and experiences a phase-transition-like discontinuity in the thermodynamic limit N → ∞. The zeroes of the partition function tend to distribute on a slightly distorted unit circle in complex plane and approach the positive real axis already for a 5×5 square lattice. When the lattice has negative cooperativity, its properties mimic those of a system composed of two classes of independent sites with the apparent population of low-affinity binding sites increasing with the size of the lattice, thereby accounting for a phenomenon encountered in many ligand-receptor interactions.

Dual role of the nuclear factor of activated T cells insert region in DNA recognition and cooperative contacts to activator protein 1

The transcription factors nuclear factor of activated T cells (NFAT) and activator protein 1 (AP-1) coordinately regulate cytokine gene expression in activated T-cells by binding to closely juxtaposed sites in cytokine promoters. The structural basis for cooperative binding of NFAT and AP-1 to these sites, and indeed for the cooperative binding of transcription factors to composite regulatory elements in general, is not well understood. Mutagenesis studies have identified a segment of AP-1, which lies at the junction of its DNA-binding and dimerization domains (basic region and leucine zipper, respectively), as being essential for protein–protein interactions with NFAT in the ternary NFAT/AP-1/DNA complex. In a model of the ternary complex, the segment of NFAT nearest AP-1 is the Rel insert region (RIR), a feature that is notable for its hypervariability in size and in sequence amongst members of the Rel transcription factor family. Here we have used mutational analysis to study the role of the NFAT RIR in binding to DNA and AP-1. Parallel yeast one-hybrid screening assays in combination with alanine-scanning mutagenesis led to the identification of four amino acid residues in the RIR of NFAT2 (also known as NFATC1 or NFATc) that are essential for cooperativity with AP-1 (Ile-544, Glu-545, Thr-551, and Ile-553), and three residues that are involved in interactions with DNA (Lys-538, Arg-540, and Asn-541). These results were confirmed and extended through in vitro binding assays. We thus conclude that the NFAT RIR plays an essential dual role in DNA recognition and cooperative binding to AP-1 family transcription factors.

The dTTPase mechanism of T7 DNA helicase resembles the binding change mechanism of the F1-ATPase

Bacteriophage T7 DNA helicase is a ring-shaped hexamer that catalyzes duplex DNA unwinding using dTTP hydrolysis as an energy source. Of the six potential nucleotide binding sites on the hexamer, we have found that three are noncatalytic sites and three are catalytic sites. The noncatalytic sites bind nucleotides with a high affinity, but dTTPs bound to these sites do not dissociate or hydrolyze through many dTTPase turnovers at the catalytic sites. The catalytic sites show strong cooperativity which leads to sequential binding and hydrolysis of dTTP. The elucidated dTTPase mechanism of the catalytic sites of T7 helicase is remarkably similar to the binding change mechanism of the ATP synthase. Based on the similarity, a general mechanism for hexameric helicases is proposed. In this mechanism, an F1-ATPase-like rotational movement around the single-stranded DNA, which is bound through the central hole of the hexamer, is proposed to lead to unidirectional translocation along single-stranded DNA and duplex DNA unwinding.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Meis1 and pKnox1 bind DNA cooperatively with Pbx1 utilizing an interaction surface disrupted in oncoprotein E2a-Pbx1

E2a-Pbx1 is a chimeric transcription factor oncoprotein produced by the t(1;19) translocation in human pre-B cell leukemia. Class I Hox proteins bind DNA cooperatively with both Pbx proteins and oncoprotein E2a-Pbx1, suggesting that leukemogenesis by E2a-Pbx1 and Hox proteins may alter transcription of cellular genes regulated by Pbx–Hox motifs. Likewise, in murine myeloid leukemia, transcriptional coactivation of Meis1 with HoxA7/A9 suggests that Meis1–HoxA7/9 heterodimers may evoke aberrant gene transcription. Here, we demonstrate that both Meis1 and its relative, pKnox1, dimerize with Pbx1 on the same TGATTGAC motif selected by dimers of Pbx proteins and unidentified partner(s) in nuclear extracts, including those from t(1;19) pre-B cells. Outside their homeodomains, Meis1 and pKnox1 were highly conserved only in two motifs required for cooperativity with Pbx1. Like the unidentified endogenous partner(s), both Meis1 and pKnox1 failed to dimerize significantly with E2a-Pbx1. The Meis1/pKnox1-interaction domain in Pbx1 resided predominantly in a conserved N-terminal Pbx domain deleted in E2a-Pbx1. Thus, the leukemic potential of E2a-Pbx1 may require abrogation of its interaction with members of the Meis and pKnox families of transcription factors, permitting selective targeting of genes regulated by Pbx–Hox complexes. In addition, because most motifs bound by Pbx–Meis1/pKnox1 were not bound by Pbx1–Hox complexes, the leukemic potential of Meis1 in myeloid leukemias may involve shifting Pbx proteins from promoters containing Pbx–Hox motifs to those containing Pbx–Meis motifs.

The PH Domain and the Polybasic c Domain of Cytohesin-1 Cooperate specifically in Plasma Membrane Association and Cellular Function

Recruitment of intracellular proteins to the plasma membrane is a commonly found requirement for the initiation of signal transduction events. The recently discovered pleckstrin homology (PH) domain, a structurally conserved element found in ∼100 signaling proteins, has been implicated in this function, because some PH domains have been described to be involved in plasma membrane association. Furthermore, several PH domains bind to the phosphoinositides phosphatidylinositol-(4,5)-bisphosphate and phosphatidylinositol-(3,4,5)-trisphosphate in vitro, however, mostly with low affinity. It is unclear how such weak interactions can be responsible for observed membrane binding in vivo as well as the resulting biological phenomena. Here, we investigate the structural and functional requirements for membrane association of cytohesin-1, a recently discovered regulatory protein of T cell adhesion. We demonstrate that both the PH domain and the adjacent carboxyl-terminal polybasic sequence of cytohesin-1 (c domain) are necessary for plasma membrane association and biological function, namely interference with Jurkat cell adhesion to intercellular adhesion molecule 1. Biosensor measurements revealed that phosphatidylinositol-(3,4,5)-trisphosphate binds to the PH domain and c domain together with high affinity (100 nM), whereas the isolated PH domain has a substantially lower affinity (2–3 μM). The cooperativity of both elements appears specific, because a chimeric protein, consisting of the c domain of cytohesin-1 and the PH domain of the β-adrenergic receptor kinase does not associate with membranes, nor does it inhibit adhesion. Moreover, replacement of the c domain of cytohesin-1 with a palmitoylation–isoprenylation motif partially restored the biological function, but the specific targeting to the plasma membrane was not retained. Thus we conclude that two elements of cytohesin-1, the PH domain and the c domain, are required and sufficient for membrane association. This appears to be a common mechanism for plasma membrane targeting of PH domains, because we observed a similar functional cooperativity of the PH domain of Bruton’s tyrosine kinase with the adjacent Bruton’s tyrosine kinase motif, a novel zinc-containing fold.

Ordered water molecules as key allosteric mediators in a cooperative dimeric hemoglobin

One of the most remarkable structural aspects of Scapharca dimeric hemoglobin is the disruption of a very well-ordered water cluster at the subunit interface upon ligand binding. We have explored the role of these crystallographically observed water molecules by site-directed mutagenesis and osmotic stress techniques. The isosteric mutation of Thr-72 → Val in the interface increases oxygen affinity more than 40-fold with a surprising enhancement of cooperativity. The only significant structural effect of this mutation is to destabilize two ordered water molecules in the deoxy interface. Wild-type Scapharca hemoglobin is strongly sensitive to osmotic conditions. Upon addition of glycerol, striking changes in Raman spectrum of the deoxy form are observed that indicate a transition toward the liganded form. Increased osmotic pressure, which lowers the oxygen affinity in human hemoglobin, raises the oxygen affinity of Scapharca hemoglobin regardless of whether the solute is glycerol, glucose, or sucrose. Analysis of these results provides an estimate of six water molecules lost upon oxygen binding to the dimer, in good agreement with eight predicted from crystal structures. These experiments suggest that the observed cluster of interfacial water molecules plays a crucial role in communication between subunits.