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].

Corticotropin-releasing factor-binding protein ligand inhibitor blunts excessive weight gain in genetically obese Zucker rats and rats during nicotine withdrawal

Elevation of the neuropeptide corticotropin-releasing factor (CRF) in the brain is associated with a reduction of food intake and body weight gain in normal and obese animals. A protein that binds CRF and the related peptide, urocortin, with high affinity, CRF-binding protein (CRF-BP), may play a role in energy homeostasis by inactivating members of this peptide family in ingestive and metabolic regulatory brain regions. Intracerebroventricular administration in rats of the high-affinity CRF-BP ligand inhibitor, rat/human CRF (6-33), which dissociates CRF or urocortin from CRF-BP and increases endogenous brain levels of “free” CRF or urocortin significantly blunted exaggerated weight gain in Zucker obese subjects and in animals withdrawn from chronic nicotine. Chronic administration of CRF suppressed weight gain nonselectively by 60% in both Zucker obese and lean control rats, whereas CRF-BP ligand inhibitor treatment significantly reduced weight gain in obese subjects, without altering weight gain in lean control subjects. Nicotine abstinent subjects, but not nicotine-naive controls, experienced a 35% appetite suppression and a 25% weight gain reduction following acute and chronic administration, respectively, of CRF-BP ligand inhibitor. In marked contrast to the effects of a CRF-receptor agonist, the CRF-BP ligand inhibitor did not stimulate adrenocorticotropic hormone secretion or elevate heart rate and blood pressure. These results provide support for the hypothesis that the CRF-BP may function within the brain to limit selected actions of CRF and/or urocortin. Furthermore, CRF-BP may represent a novel and functionally selective target for the symptomatic treatment of excessive weight gain associated with obesity of multiple etiology.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Role of metallothionein in nitric oxide signaling as revealed by a green fluorescent fusion protein

Although the function of metallothionein (MT), a 6- to 7-kDa cysteine-rich metal binding protein, remains unclear, it has been suggested from in vitro studies that MT is an important component of intracellular redox signaling, including being a target for nitric oxide (NO). To directly study the interaction between MT and NO in live cells, we generated a fusion protein consisting of MT sandwiched between two mutant green fluorescent proteins (GFPs). In vitro studies with this chimera (FRET-MT) demonstrate that fluorescent resonance energy transfer (FRET) can be used to follow conformational changes indicative of metal release from MT. Imaging experiments with live endothelial cells show that agents that increase cytoplasmic Ca2+ act via endogenously generated NO to rapidly and persistently release metal from MT. A role for this interaction in intact tissue is supported by the finding that the myogenic reflex of mesenteric arteries is absent in MT knockout mice (MT−/−) unless endogenous NO synthesis is blocked. These results are the first application of intramolecular green fluorescent protein (GFP)-based FRET in a native protein and demonstrate the utility of FRET-MT as an intracellular surrogate indicator of NO production. In addition, an important role of metal thiolate clusters of MT in NO signaling in vascular tissue is revealed.

Native protein sequences are close to optimal for their structures

How large is the volume of sequence space that is compatible with a given protein structure? Starting from random sequences, low free energy sequences were generated for 108 protein backbone structures by using a Monte Carlo optimization procedure and a free energy function based primarily on Lennard–Jones packing interactions and the Lazaridis–Karplus implicit solvation model. Remarkably, in the designed sequences 51% of the core residues and 27% of all residues were identical to the amino acids in the corresponding positions in the native sequences. The lowest free energy sequences obtained for ensembles of native-like backbone structures were also similar to the native sequence. Furthermore, both the individual residue frequencies and the covariances between pairs of positions observed in the very large SH3 domain family were recapitulated in core sequences designed for SH3 domain structures. Taken together, these results suggest that the volume of sequence space optimal for a protein structure is surprisingly restricted to a region around the native sequence.

A meanfield approach to the thermodynamics of a protein-solvent system with application to the oligomerization of the tumor suppressor p53

The thermodynamic stability and oligomerization status of the tumor suppressor p53 tetramerization domain have been studied experimentally and theoretically. A series of hydrophilic mutations at Met-340 and Leu-344 of human p53 were designed to disrupt the hydrophobic dimer–dimer interface of the tetrameric oligomerization domain of p53 (residues 325–355). Meanfield calculations of the free energy of the solvated mutants as a function of interdimer distance were compared with experimental data on the thermal stability and oligomeric state (tetramer, dimer, or equilibrium mixture of both) of each mutant. The calculations predicted a decreasing stability and oligomeric state for the following amino acids at residue 340: Met (tetramer) > Ser Asp, His, Gln, > Glu, Lys (dimer), whereas the experimental results showed the following order: Met (tetramer) > Ser > Gln > His, Lys > Asp, Glu (dimers). For residue 344, the calculated trend was Leu (tetramer) > Ala > Arg, Gln, Lys (dimer), and the experimental trend was Leu (tetramer) > Ala, Arg, Gln, Lys (dimer). The discrepancy for the lysine side chain at residue 340 is attributed to the dual nature of lysine, both hydrophobic and charged. The incorrect prediction of stability of the mutant with Asp at residue 340 is attributed to the fact that within the meanfield approach, we use the wild-type backbone configuration for all mutants, but low melting temperatures suggest a softening of the α-helices at the dimer–dimer interface. Overall, this initial application of meanfield theory toward a protein-solvent system is encouraging for the application of the theoretical model to more complex systems.

Solvent dependence of dynamic transitions in protein solutions

A transition as a function of increasing temperature from harmonic to anharmonic dynamics has been observed in globular proteins by using spectroscopic, scattering, and computer simulation techniques. We present here results of a dynamic neutron scattering analysis of the solvent dependence of the picosecond-time scale dynamic transition behavior of solutions of a simple single-subunit enzyme, xylanase. The protein is examined in powder form, in D2O, and in four two-component perdeuterated single-phase cryosolvents in which it is active and stable. The scattering profiles of the mixed solvent systems in the absence of protein are also determined. The general features of the dynamic transition behavior of the protein solutions follow those of the solvents. The dynamic transition in all of the mixed cryosolvent–protein systems is much more gradual than in pure D2O, consistent with a distribution of energy barriers. The differences between the dynamic behaviors of the various cryosolvent protein solutions themselves are remarkably small. The results are consistent with a picture in which the picosecond-time scale atomic dynamics respond strongly to melting of pure water solvent but are relatively invariant in cryosolvents of differing compositions and melting points.

Novel role of phosphorylation in Fe–S cluster stability revealed by phosphomimetic mutations at Ser-138 of iron regulatory protein 1

Animals regulate iron metabolism largely through the action of the iron regulatory proteins (IRPs). IRPs modulate mRNA utilization by binding to iron-responsive elements (IRE) in the 5′ or 3′ untranslated region of mRNAs encoding proteins involved in iron homeostasis or energy production. IRP1 is also the cytosolic isoform of aconitase. The activities of IRP1 are mutually exclusive and are modulated through the assembly/disassembly of its [4Fe–4S] cluster, reversibly converting it between an IRE-binding protein and cytosolic aconitase. IRP1 is also phosphoregulated by protein kinase C, but the mechanism by which phosphorylation posttranslationally increases IRE binding activity has not been fully defined. To investigate this, Ser-138 (S138), a PKC phosphorylation site, was mutated to phosphomimetic glutamate (S138E), aspartate (S138D), or nonphosphorylatable alanine (S138A). The S138E IRP1 mutant and, to a lesser extent, the S138D IRP1 mutant were impaired in aconitase function in yeast when grown aerobically but not when grown anaerobically. Purified wild-type and mutant IRP1s could be reconstituted to active aconitases anaerobically. However, when exposed to oxygen, the [4Fe–4S] cluster of the S138D and S138E mutants decayed 5-fold and 20-fold faster, respectively, than was observed for wild-type IRP1. Our findings suggest that stability of the Fe–S cluster of IRP1 can be regulated by phosphorylation and reveal a mechanism whereby the balance between the IRE binding and [4Fe–4S] forms of IRP1 can be modulated independently of cellular iron status. Furthermore, our results show that IRP1 can function as an oxygen-modulated posttranscriptional regulator of gene expression.

Protein fluctuations are sensed by stimulated infrared echoes of the vibrations of carbon monoxide and azide probes

The correlation functions of the fluctuations of vibrational frequencies of azide ions and carbon monoxide in proteins are determined directly from stimulated photon echoes generated with femtosecond infrared pulses. The asymmetric stretching vibration of azide bound to carbonic anhydrase II exhibits a pronounced evolution of its vibrational frequency distribution on the time scale of a few picoseconds, which is attributed to modifications of the ligand structure through interactions with the nearby Thr-199. When azide is bound in hemoglobin, a more complex evolution of the protein structure is required to interchange the different ligand configurations, as evidenced by the much slower relaxation of the frequency distribution in this case. The time evolution of the distribution of frequencies of carbon monoxide bound in hemoglobin occurs on the ≈10-ps time scale and is very nonexponential. The correlation functions of the frequency fluctuations determine the evolution of the protein structure local to the probe and the extent to which the probe can navigate those parts of the energy landscape where the structural configurations are able to modify the local potential energy function of the probe.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Proteolytic cleavage product of 30-kDa adipocyte complement-related protein increases fatty acid oxidation in muscle and causes weight loss in mice

Adipocyte complement-related protein (30 kDa) (Acrp30), a secreted protein of unknown function, is exclusively expressed in differentiated adipocytes; its mRNA is decreased in obese humans and mice. Here we describe novel pharmacological properties of the protease-generated globular head domain of Acrp30 (gAcrp30). Acute treatment of mice with gAcrp30 significantly decreased the elevated levels of plasma free fatty acids caused either by administration of a high fat test meal or by i.v. injection of Intralipid. This effect of gAcrp30 was caused, at least in part, by an acute increase in fatty acid oxidation by muscle. As a result, daily administration of a very low dose of gAcrp30 to mice consuming a high-fat/sucrose diet caused profound and sustainable weight reduction without affecting food intake. Thus, gAcrp30 is a novel pharmacological compound that controls energy homeostasis and exerts its effect primarily at the peripheral level.

Protein-free parallel triple-stranded DNA complex formation

A 14 nt DNA sequence 5′-AGAATGTGGCAAAG-3′ from the zinc finger repeat of the human KRAB zinc finger protein gene ZNF91 bearing the intercalator 2-methoxy,6-chloro,9-amino acridine (Acr) attached to the sugar–phosphate backbone in various positions has been shown to form a specific triple helix (triplex) with a 16 bp hairpin (intramolecular) or a two-stranded (intermolecular) duplex having the identical sequence in the same (parallel) orientation. Intramolecular targets with the identical sequence in the antiparallel orientation and a non-specific target sequence were tested as controls. Apparent binding constants for formation of the triplex were determined by quantitating electrophoretic band shifts. Binding of the single-stranded oligonucleotide probe sequence to the target led to an increase in the fluorescence anisotropy of acridine. The parallel orientation of the two identical sequence segments was confirmed by measurement of fluorescence resonance energy transfer between the acridine on the 5′-end of the probe strand as donor and BODIPY-Texas Red on the 3′-amino group of either strand of the target duplex as acceptor. There was full protection from OsO4-bipyridine modification of thymines in the probe strand of the triplex, in accordance with the presumed triplex formation, which excluded displacement of the homologous duplex strand by the probe–intercalator conjugate. The implications of these results for the existence of protein-independent parallel triplexes are discussed.

Human MTH1 protein hydrolyzes the oxidized ribonucleotide, 2-hydroxy-ATP

The human nucleotide pool sanitization enzyme, MTH1, hydrolyzes 2-hydroxy-dATP and 8-hydroxy-dATP in addition to 8-hydroxy-dGTP. We report here that human MTH1 is highly specific for 2-hydroxy-ATP, among the cognate ribonucleoside triphosphates. The pyrophosphatase activities for 8-hydroxy-GTP, 2-hydroxy-ATP and 8-hydroxy-ATP were measured by high-performance liquid chromatography. The kinetic parameters thus obtained indicate that the catalytic efficiencies of MTH1 are in the order of 2-hydroxy-dATP > 2-hydroxy-ATP > 8-hydroxy-dGTP > 8-hydroxy-dATP >> dGTP > 8-hydroxy-GTP > 8-hydroxy-ATP. Notably, MTH1 had the highest affinity for 2-hydroxy-ATP among the known substrates. ATP is involved in energy metabolism and signal transduction, and is a precursor in RNA synthesis. We suggest that the 2-hydroxy-ATP hydrolyzing activity of MTH1 might prevent the perturbation of these ATP-related pathways by the oxidized ATP.

The role of structure, energy landscape, dynamics, and allostery in the enzymatic function of myoglobin

The grail of protein science is the connection between structure and function. For myoglobin (Mb) this goal is close. Described as only a passive dioxygen storage protein in texts, we argue here that Mb is actually an allosteric enzyme that can catalyze reactions among small molecules. Studies of the structural, spectroscopic, and kinetic properties of Mb lead to a model that relates structure, energy landscape, dynamics, and function. Mb functions as a miniature chemical reactor, concentrating and orienting diatomic molecules such as NO, CO, O2, and H2O2 in highly conserved internal cavities. Reactions can be controlled because Mb exists in distinct taxonomic substates with different catalytic properties and connectivities of internal cavities.

Kinetics of protein import into isolated Xenopus oocyte nuclei

An in vitro assay for nucleocytoplasmic transport was established in which signal-dependent protein import is reproduced faithfully by isolated purified nuclei. The assay permits the precise quantification of import kinetics and the discrimination between translocation through the nuclear envelope and intranuclear transport. Nuclei were manually isolated from Xenopus oocytes and after manual purification incubated with a medium containing a green fluorescent transport substrate, karyopherins α2 and β1, a red fluorescent control substrate, an energy mix and, for keeping an osmotic balance, 20% (wt/vol) BSA. Import of transport substrates into the nucleus and exclusion of the control substrate were monitored simultaneously by two-color confocal microscopy. Two widely differing import substrates were used: the recombinant protein P4K [480 kDa, four nuclear localization sequences (NLSs) per P4K tetramer], and NLS-BSA (90 kDa, 15 NLSs). The measurements suggested that import, at the specific conditions used in this study, consisted of two consecutive processes: (i) the rapid equilibration of the concentration difference across the nuclear envelope, a process involving binding and translocation of substrate by the nuclear pore complex , and (ii) the dissipation of the intranuclear concentration difference by diffusion.