Mutagenic evidence for the optimal control of evolutionary dynamics.

Elucidating the fitness measures optimized during the evolution of complex biological systems is a major challenge in evolutionary theory. We present experimental evidence and an analytical framework demonstrating how biochemical networks exploit optimal control strategies in their evolutionary dynamics. Optimal control theory explains a striking pattern of extremization in the redox potentials of electron transport proteins, assuming only that their fitness measure is a control objective functional with bounded controls.

Biochemical basis for enhanced binding of peptide dimers to X-linked inhibitor of apoptosis protein.

XIAP (X-linked inhibitor of apoptosis protein) is involved in the mediation of programmed cell death and, therefore, is a target for the development of cancer therapeutics. Peptide mimetics based upon Smac, the natural binding partner of XIAP, and specifically, dimeric peptides, have shown great promise in drug development. In the present work, the basis for enhanced dimer efficacy has been explored. Comparisons are made between the peptide binding site on the BIR3 domain of XIAP alone (residues 238-358) and a less truncated construct that includes both BIR2 and BIR3 domains (residues 151-350). This contingency differentially enhances the binding of dimeric tetrapeptides, potentially by providing additional hydrophobic binding surface. The effect of BIR2 on the BIR3 binding site is sustained, even if the BIR2 binding site is disrupted by mutagenesis, as shown by both a fluorescent competition assay and a polarity sensitive dye, badan. FRET measurements reveal an observed separation of >or=45 A between the BIR2 and BIR3 peptide binding pockets, thereby precluding a direct simultaneous interaction of the dimer molecules with both binding domains. Furthermore, variations in the linker length between dimeric tetrapeptides did not show a predictable trend in binding affinities, suggesting that local concentration effects were also an unlikely explanation for the enhanced dimeric affinities. Taken together, the results suggest that enhanced binding of dimeric peptides likely reflects the increased hydrophobic surface area on or near the BIR3 site and have significant ramifications for the design of therapeutics that target this class of proteins.

Structure-activity based study of the Smac-binding pocket within the BIR3 domain of XIAP.

A small series of peptide mimics was designed and synthesized to contain a heterocyclic ring in place of the potentially labile N-terminal peptide bond of the tetrapeptide containing the Smac-XIAP-binding motif. Two Smac mimics were shown to bind to the BIR3 domain of XIAP with moderate affinity and one displayed increased activity in cells relative to the Smac peptides. The structures of BIR3-XIAP in complex with a Smac peptide and a peptide mimic were solved and analyzed to elucidate the structure-activity relationship surrounding the Smac-binding domain within BIR3-XIAP.

Tryptophan and tyrosine to terbium fluorescence resonance energy transfer as a method to "map" aromatic residues and monitor docking.

Fluorescent lanthanide ions, with large Stokes shifts and narrow emission bands, are excellent tools for the development of FRET-based assays. In this work, a terbium ion is tethered to a peptide which binds to the BIR3 domain of XIAP, an anti-apoptotic protein. Excitation of tryptophan and tyrosine residues in the BIR3 domain causes the peptide bound terbium ion to fluoresce relative to its distance from these aromatic residues. By developing ligands with terbium ions tethered at different residues, the relative terbium emission can be used to "map" the aromatic residues within the ligand binding pocket.

Endonuclease-like activity of heme proteins.

Heme proteins, metmyoglobin, methemoglobin, and metcytochrome c showed unusual affinity for double-stranded DNA. Calorimetric studies show that binding of methemoglobin to calf thymus DNA (CTDNA) is weakly endothermic, and the binding constant is 4.9+/-0.7x10(5) M(-1). The Soret absorption bands of the heme proteins remained unchanged, in the presence of excess CTDNA, but a new circular dichroic band appeared at 210 nm. Helix melting studies indicated that the protein-DNA mixture denatures at a lower temperature than the individual components. Thermograms obtained by differential scanning calorimetry of the mixture indicated two distinct transitions, which are comparable to the thermograms obtained for individual components, but there was a reduction in the excess heat capacity. Activation of heme proteins by hydrogen peroxide resulted in the formation of high valent Fe(IV) oxo intermediates, and CTDNA reacted rapidly under these conditions. The rate was first-order in DNA concentration, and this reactivity resulted in DNA strand cleavage. Upon activation with hydrogen peroxide, for example, the heme proteins converted the supercoiled pUC18 DNA into nicked circular and linear DNA. No reaction occurred in the absence of the heme protein, or hydrogen peroxide. These data clearly indicate a novel property of several heme proteins, and this is first report of the endonuclease-like activity of the heme proteins.

Design, folding, and activities of metal-assembled coiled coil proteins.

Metal ions serve many purposes in natural proteins, from the stabilization of tertiary structure to the direction of protein folding to crucial roles in electron transfer and catalysis. There is considerable interest in creating metal binding sites in designed proteins to understand the structural role of metal ions and to design new metalloproteins with useful functions. The de novo design of metalloproteins and the role of metals in the folding of designed proteins are reviewed here, with particular focus on the design, folding, and activities of the [M(bpy-peptide)(3)](2+) structure. This maquette is constructed by the covalent attachment of 2,2'-bipyridine to the N-termini of amphiphilic peptides, and it is assembled into a folded trimeric coiled coil by the addition of a six-coordinate transition metal ion and the resulting hydrophobic collapse of the peptides. The [M(bpy-peptide)(3)](2+) structure has been employed in diverse applications, ranging from electron transfer pathway studies to the study of optimal hydrophobic packing in a virtual library to the construction of receptors and biosensors.

Metal-assembled modular proteins: toward functional protein design.

Metal-assembled parallel helix-bundle proteins have been used to investigate electron transfer through alpha-helical structures. Fermi Golden Rule distance dependence of electron transfer rates was established in a family of designed metalloproteins, and the contribution of intrahelical hydrogen bonding to the matrix tunneling element was explored. The first steps toward the design of functional proteins using dynamic combinatorial assembly of alpha-helical structural elements are described.

Time-resolved absorption and UV resonance Raman spectra reveal stepwise formation of T quaternary contacts in the allosteric pathway of hemoglobin.

Hemoglobin undergoes a series of molecular changes on the nanosecond and microsecond time-scale following photodissociation of CO ligands. We have monitored these processes with a combination of transient absorption and resonance Raman (RR) spectroscopy. The latter have been acquired at higher data rates than previously available, thanks to kilohertz Ti:sapphire laser technology, with frequency-quadrupling into the ultraviolet. As a result of improved resolution of the UVRR time-course, a new intermediate has been identified in the pathway from the R (HbCO) to the T (deoxyHb) state. This intermediate is not detected via absorption transients, since the change in heme absorption is insignificant, but its lifetime agrees with a reported magnetic circular dichroism transient, which has been attributed to a quaternary tryptophan interaction. The new UVRR data allow elaboration of the allosteric pathway by establishing that the T-state quaternary contacts are formed in two well-separated steps, with time constants of 2.9 micros and 21 micros, instead of a single 20 micros process. The first step involves the "hinge" region contacts, as monitored by the Trp beta 37...Asp alpha 94 H-bond, while the second involves the "switch" region, as monitored by the Tyr alpha 42...Asp beta 99 H-bond. A working model for the allosteric pathway is presented.

Hemoglobin site-mutants reveal dynamical role of interhelical H-bonds in the allosteric pathway: time-resolved UV resonance Raman evidence for intra-dimer coupling.

The dynamical effect of eliminating specific tertiary H-bonds in the hemoglobin (Hb) tetramer has been investigated by site-directed mutagenesis and time-resolved absorption and ultraviolet resonance Raman (UVRR) spectroscopy. The Trp alpha 14...Thr alpha 67 and Trp beta 15...Ser beta 72 H-bonds connect the A and E helices in the alpha and beta chains, and are proposed to break in the earliest protein intermediate (Rdeoxy) following photo-deligation of HbCO, along with a second pair of H-bonds involving tyrosine residues. Mutation of the acceptor residues Thr alpha 67 and Ser beta 72 to Val and Ala eliminates the A-E H-bonds, but has been shown to have no significant effect on ligand-binding affinity or cooperativity, or on spectroscopic markers of the T-state quaternary interactions. However, the mutations have profound and unexpected effects on the character of the Rdeoxy intermediate, and on the dynamics of the subsequent steps leading to the T state. Formation of the initial quaternary contact (RT intermediate) is accelerated, by an order of magnitude, but the locking-in of the T state is delayed by a factor of 2. These rate effects are essentially the same for either mutation, or for the double mutation, suggesting that the alpha beta dimer behaves as a mechanically coupled dynamical unit. Further evidence for intra-dimer coupling is provided by the Rdeoxy UVRR spectrum, in which either or both mutations eliminate the tyrosine difference intensity, although only tryptophan H-bonds are directly affected. A possible mechanism for mechanical coupling is outlined, involving transmission of forces through the alpha(1)beta(1) (and alpha(2)beta(2)) interface. The present observations establish that quaternary motions can occur on the approximately 100 ns time-scale. They show also that a full complement of interhelical H-bonds actually slows the initial quaternary motion in Hb, but accelerates the locking in of the T-contacts.

Design of a functional protein for molecular recognition: specificity of ligand binding in a metal-assembled protein cavity probed by 19f NMR.

A metal-assembled homotrimeric coiled coil based on the GCN4-p1 sequence has been designed that noncovalently binds hexafluorobenzene and other similar ligands in a hydrophobic cavity, created by making the core substitution Asn16Ala ([Fe(bpyGCN4-N16A)3]2+). The KD of binding of hexafluorobenzene with [Fe(bpyGCN4-N16A)3]2+ was observed to be 1.1(9) x 10(-4) M by diffusion NMR experiments. A control coiled coil with the core substitution Asn16Val ([Fe(bpyGCN4-N16V)3]2+) exhibited a significantly weaker association with hexafluorobenzene, providing evidence that even in the absence of structural data, benzene-like ligands bind in the cavity created by the Asn16Ala substitution. 19F NMR was employed to observe hexafluorobenzene binding and to monitor titrations with competing hydrophobic and polar ligands similar in size and shape to hexafluorobenzene. All hydrophobic ligands bound with greater affinity than the polar ligands in the hydrophobic core, although the cavity seems to be somewhat flexible in terms of the sizes of molecules it can accommodate. Thus 19F NMR has proved to be a useful spectral tool to probe molecular recognition in a hydrophobic cavity of a metal-assembled coiled coil.

Dynamics of carbon monoxide binding to CooA.

CooA is a dimeric CO-sensing heme protein from Rhodospirillum rubrum. The heme iron in reduced CooA is six-coordinate; the axial ligands are His-77 and Pro-2. CO displaces Pro-2 and induces a conformation change that allows CooA to bind DNA and activate transcription of coo genes. Equilibrium CO binding is cooperative, with a Hill coefficient of n = 1.4, P(50) = 2.2 microm, and estimated Adair constants K(1) = 0.16 and K(2) = 1.3 microm(-1). The rates of CO binding and release are both strongly biphasic, with roughly equal amplitudes for the fast and slow phases. The association rates show a hyperbolic dependence on [CO], consistent with Pro-2 dissociation being rate-limiting. The kinetic characteristics of the transiently formed five-coordinate heme are probed via flash photolysis. These observations are integrated into a kinetic model, in which CO binding to one subunit decreases the rate of Pro-2 rebinding in the second, leading to a net increase in affinity for the second CO. The CO adduct exists in slowly interconverting "open" and "closed" forms. This interconversion probably involves the large-scale motions required to bring the DNA-binding domains into proper orientation. The combination of low CO affinity, slow CO binding, and slow conformational transitions ensures that activation of CooA only occurs at high (micromolar) and sustained (> or =1 min) levels of CO. When micromolar levels do occur, positive cooperativity allows efficient activation over a narrow range of CO concentrations.

Effect of redox potential of the heme on the peroxidase activity of cytochrome b562.

Measurements of peroxidase activities of two site-specific mutants and wild type cytochrome b562 suggest that the enzymatic activity correlates with the redox potential of the metal center. A lower value of the Fe(3+)/Fe(2+) redox potential seems to be important for promoting peroxidase activity of the hemeprotein possibly by stabilization of the high-valent redox intermediate involved in the catalytic function. The results provide an approach towards rational tuning of enzyme function when 'grafted' into a new protein environment.

A high-throughput screen for identification of molecular mimics of Smac/DIABLO utilizing a fluorescence polarization assay.

Resistance to apoptosis is afforded by inhibitor of apoptosis proteins (IAPs) which bind to and inhibit the caspases responsible for cleavage of substrates leading to apoptotic cell death. Smac (or DIABLO), a proapoptotic protein released from the mitochondrial intermembrane space into the cytosol, promotes apoptosis by binding to IAPs, thus reversing their inhibitory effects on caspases. We have developed a high-throughput fluorescence polarization assay utilizing a fluorescein-labeled peptide similar to the "IAP binding" domain of Smac N terminus complexed with the BIR3 domain of X-linked IAP (XIAP) to identify small-molecule mimics of the action of Smac. The IC(50)s of peptides and a tetrapeptidomimetic homologous to the N terminus of Smac demonstrated the specificity and utility of this assay. We have screened the National Cancer Institute "Training Set" of 230 compounds, with well-defined biological actions, and the "Diversity Set" of 2000 chemically diverse structures for compounds which significantly reduced fluorescence polarization. Highly fluorescing or fluorescence-quenching compounds (false positives) were distinguished from those which interfered with Smac peptide binding to the XIAP-BIR3 in a dose-dependent manner (true positives). This robust assay offers potential for high-throughput screening discovery of novel compounds simulating the action of Smac/DIABLO.

Electrospray ionization Fourier transform ion cyclotron resonance mass spectrometric analysis of metal-ion selected dynamic protein libraries.

The application of electrospray ionization (ESI) Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry to the investigation of the relative stabilities (and thus packing efficiencies) of Fe-bound trihelix peptide bundles is demonstrated. Small dynamic protein libraries are created by metal-ion assisted assembly of peptide subunits. Control of the trimeric aggregation state is coupled to stability selection by exploiting the coordination requirements of Fe(2+) in the presence of bidentate 2,2'-bipyridyl ligands covalently appended to the peptide monomers. At limiting metal-ion concentration, the most thermodynamically stable, optimally packed peptide trimers dominate the mass spectrum. The identities of optimally stable candidate trimers observed in the ESI FT-ICR mass spectra are confirmed by resynthesis of exchange-inert analogues and measurement of their folding free energies. The peptide composition of the trimers may be determined by infrared multiphoton dissociation (IRMPD) MS(3) experiments. Additional sequence information for the peptide subunits is obtained from electron capture dissociation (ECD) of peptides and metal-bound trimers. The experiments also suggest the presence of secondary structure in the gas phase, possibly due to partial retention of the solution-phase coiled coil structure.

Kinetics of proton-coupled electron-transfer reactions to the manganese-oxo "cubane" complexes containing the Mn4O4(6+) and Mn4O4(7+) core types.

The kinetics of proton-coupled electron-transfer (pcet) reactions are reported for Mn(4)O(4)(O(2)PPh(2))(6), 1, and [Mn(4)O(4)(O(2)PPh(2))(6)](+), 1(+), with phenothiazine (pzH). Both pcet reactions form 1H, by H transfer to 1 and by hydride transfer to 1(+). Surprisingly, the rate constants differ by only 25% despite large differences in the formal charges and driving force. The driving force is proportional to the difference in the bond-dissociation energies (BDE >94 kcalmol for homolytic, 1H --> H + 1, vs. approximately 127 kcalmol for heterolytic, 1H --> H(-) + 1(+), dissociation of the OH bond in 1H). The enthalpy and entropy of activation for the homolytic reaction (deltaH = -1.2 kcalmol and deltaS= -32 calmol.K; 25-6.7 degrees C) reveal a low activation barrier and an appreciable entropic penalty in the transition state. The rate-limiting step exhibits no HD kinetic isotope effect (k(H)k(D) = 0.96) for the first H atom-transfer step and a small kinetic isotope effect (1.4) for the second step (1H + pzH --> 1H(2) + pz(*)). These lines of evidence indicate that formation of a reactive precursor complex before atom transfer is rate-limiting (conformational gating), and that little or no NH bond cleavage occurs in the transition state. H-atom transfer from pzH to alkyl, alkoxyl, and peroxyl radicals reveals that BDEs are not a good predictor of the rates of this reaction. Hydride transfer to 1(+) provides a concrete example of two-electron pcet that is hypothesized for the OH bond cleavage step during catalysis of photosynthetic water oxidation.

The kinetics of translocation of Smac/DIABLO from the mitochondria to the cytosol in HeLa cells.

Smac (second mitochondrial activator of caspases) is released from the mitochondria during apoptosis to relieve inhibition of caspases by the inhibitor of apoptosis proteins (IAPs). The release of Smac antagonizes several IAPs and assists the initiator caspase-9 and effector caspases (caspase-3, caspase-6, and caspase-7) in becoming active, ultimately leading to death of the cell. Translocation of Smac along with cytochrome c and other mitochondrial pro-apoptotic proteins represent important regulatory checkpoints for mitochondria-mediated apoptosis. Whether Smac and cytochrome c translocate by the same mechanism is not known. Here, we show that the time required for Smac efflux from the mitochondria of cells subjected to staurosporine-induced apoptosis is approximately four times longer than the time required for cytochrome c efflux. These results suggest that Smac and cytochrome c may exit the mitochondria by different pathways.

Molecular targeting of inhibitor of apoptosis proteins based on small molecule mimics of natural binding partners.

An assay based on a solvent-sensitive fluorogenic dye molecule, badan, is used to test the binding affinity of a library of tetrapeptide molecules for the BIR3 (baculovirus IAP repeat) domain of XIAP (X-linked inhibitor of apoptosis protein). The fluorophore is attached to a tetrapeptide, Ala-Val-Pro-Cys-NH(2), through a thiol linkage and, upon binding to XIAP, undergoes a solvatochromic shift in fluorescence emission. When a molecule (e.g., a natural protein known to bind to XIAP or a tetrapeptide mimic) displaces the dye, the emission shifts back to the spectrum observed in water. As emission intensity is related to the binding of the tetrapeptide, the intensity can be used to determine the equilibrium constant, K, for the displacement of the dye by the tetrapeptide. The results permit residue-specific analysis of the interaction. Furthermore, we show that hydrophobic effects in the fourth position are general and can effectively increase overall affinity.

A multigeneration analysis of cytochrome b(562) redox variants: evolutionary strategies for modulating redox potential revealed using a library approach.

The redox potential of cytochromes sets the energy yield possible in metabolism and is also a key determinant of the rate at which redox reactions proceed. Here, the heme protein, cytochrome b(562), is used to study the in vitro evolution of redox potential within a library of variants containing the same structural archetype, the four-helix bundle. Multisite variations in the active site of cytochrome b(562) were introduced. A library of variants containing random mutations in place of R98 and R106 was created, and the redox potentials of a statistical sampling of this library were measured. This procedure was carried out for both the low- and high-potential variants of a previously studied F61X/F65X, first-generation library [Springs, S. L., Bass, S. E., and McLendon, G. L. (2000) Biochemistry 39, 6075]. The second-generation library reported here has a range of redox potentials which is greater than 40% (160 mV) of the known accessible potential among cytochromes with identical axial ligands (but different folds) and exceeds the range exhibited phylogenetically by the cytochrome c' family which internally maintains the same axial ligation and fold. A statistical analysis of the libraries examined reveals that the redox potential of WT cyt b(562) is found at the high-potential extremum of the distribution, indicating that this protein apparently evolved to differentially stabilize the reduced protein. The 2.7 A crystal structure of F61I/F65Y/R106L (low-potential variant of the second-generation library) was solved and is compared to the wild-type structure and the 2.2 A resolution structure of the F61I/F65Y variant (low-potential variant of the first-generation library). The structures indicate that charge-dipole effects are responsible for shifting the redox equilibrium toward the oxidized state in both the F61I/F65Y and F61I/F65Y/R106L variants. Specifically, a new protein dipole is introduced into the heme microenvironment as a result of the F65Y mutation, two new internal water molecules (one in hydrogen-bonding distance of Y65) are found, and in the case of F61I/F65Y/R106L (DeltaE(m) = 158 mV vs NHE), increased solvent exposure of the heme as a result of the R106L substitution is identified.