Structural mimicry of a native protein by a minimized binding domain

The affinity between molecules depends both on the nature and presentation of the contacts. Here, we observe coupling of functional and structural elements when a protein binding domain is evolved to a smaller functional mimic. Previously, a 38-residue form of the 59-residue B-domain of protein A, termed Z38, was selected by phage display. Z38 contains 13 mutations and binds IgG only 10-fold weaker than the native B-domain. We present the solution structure of Z38 and show that it adopts a tertiary structure remarkably similar to that observed for the first two helices of B-domain in the B-domain/Fc complex [Deisenhofer, J. (1981) Biochemistry 20, 2361–2370], although it is significantly less stable. Based on this structure, we have improved on Z38 by designing a 34-residue disulfide-bonded variant (Z34C) that has dramatically enhanced stability and binds IgG with 9-fold higher affinity. The improved stability of Z34C led to NMR spectra with much greater chemical shift dispersion, resulting in a more precisely determined structure. Z34C, like Z38, has a structure virtually identical to the equivalent region from native protein A domains. The well-defined hydrophobic core of Z34C reveals key structural features that have evolved in this small, functional domain. Thus, the stabilized two-helix peptide, about half the size and having one-third of the remaining residues altered, accurately mimics both the structure and function of the native domain.

Structure and function in rhodopsin: Packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain are coupled*

A previous study of the retinitis pigmentosa mutation L125R and two designed mutations at this site, L125A and L125F, showed that these mutations cause partial or total misfolding of the opsins expressed in COS cells from the corresponding mutant opsin genes. We now report on expression and characterization of the opsins from the following retinitis pigmentosa mutants in the transmembrane domain of rhodopsin that correspond to six of the seven helices: G51A and G51V (helix A), G89D (helix B), A164V (helix D), H211P (helix E), P267L and P267R (helix F), and T297R (helix G). All the mutations caused partial misfolding of the opsins as observed by the UV/visible absorption characteristics and by separation of the expressed opsins into fractions that bound 11-cis-retinal to form the corresponding mutant rhodopsins and those that did not bind 11-cis-retinal. Further, all the mutant rhodopsins prepared from the above mutants, except for G51A, showed strikingly abnormal bleaching behavior with abnormal metarhodopsin II photointermediates. The results show that retinitis pigmentosa mutations in every one of the transmembrane helices can cause misfolding of the opsin. Therefore, on the basis of these and previous results, we conclude that defects in the packing of the transmembrane helices resulting from these mutations are relayed to the intradiscal domain, where they cause misfolding of the opsin by inducing the formation of a disulfide bond other than the native Cys-110—Cys-187 disulfide bond. Thus, there is coupling between packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain.

A link between protein structure and enzyme catalyzed hydrogen tunneling

We present evidence that the size of an active site side chain may modulate the degree of hydrogen tunneling in an enzyme-catalyzed reaction. Primary and secondary kH/kT and kD/kT kinetic isotope effects have been measured for the oxidation of benzyl alcohol catalyzed by horse liver alcohol dehydrogenase at 25°C. As reported in earlier studies, the relationship between secondary kH/kT and kD/kT isotope effects provides a sensitive probe for deviations from classical behavior. In the present work, catalytic efficiency and the extent of hydrogen tunneling have been correlated for the alcohol dehydrogenase-catalyzed hydride transfer among a group of site-directed mutants at position 203. Val-203 interacts with the opposite face of the cofactor NAD+ from the alcohol substrate. The reduction in size of this residue is correlated with diminished tunneling and a two orders of magnitude decrease in catalytic efficiency. Comparison of the x-ray crystal structures of a ternary complex of a high-tunneling (Phe-93 → Trp) and a low-tunneling (Val-203 → Ala) mutant provides a structural basis for the observed effects, demonstrating an increase in the hydrogen transfer distance for the low-tunneling mutant. The Val-203 → Ala ternary complex crystal structure also shows a hyperclosed interdomain geometry relative to the wild-type and the Phe-93 → Trp mutant ternary complex structures. This demonstrates a flexibility in interdomain movement that could potentially narrow the distance between the donor and acceptor carbons in the native enzyme and may enhance the role of tunneling in the hydride transfer reaction.

A physical basis for protein secondary structure

A physical theory of protein secondary structure is proposed and tested by performing exceedingly simple Monte Carlo simulations. In essence, secondary structure propensities are predominantly a consequence of two competing local effects, one favoring hydrogen bond formation in helices and turns, the other opposing the attendant reduction in sidechain conformational entropy on helix and turn formation. These sequence specific biases are densely dispersed throughout the unfolded polypeptide chain, where they serve to preorganize the folding process and largely, but imperfectly, anticipate the native secondary structure.

Structure and function in rhodopsin: Topology of the C-terminal polypeptide chain in relation to the cytoplasmic loops*

Cysteine mutagenesis and site-directed spin labeling in the C-terminal region of rhodopsin have been used to probe the local structure and proximity of that region to the cytoplasmic loops. Each of the native amino acids in the sequence T335–T340 was replaced with Cys, one at a time. The sulfhydryl groups of all mutants reacted rapidly with the sulfhydryl reagent 4,4′-dithiodipyridine, which indicated a high degree of solvent accessibility. Furthermore, to probe the proximity relationships, a series of double Cys mutants was constructed. One Cys in all sets was at position 338 and the other was at a position in the sequence S240–V250 in the EF interhelical loop, at position 65 in the AB interhelical loop, or at position 140 in the CD interhelical loop. In the dark state, no significant disulfide formation was observed between C338 and C65 or C140 under the conditions used, whereas a relatively rapid disulfide formation was observed between C338 and C242 or C245. Spin labels in the double Cys mutants showed the strongest magnetic interactions between the nitroxides attached to C338 and C245 or C246. Light activation of the double mutant T242C/S338C resulted in slower disulfide formation, whereas interactions between nitroxides at C338 and C245 or C246 decreased. These results suggest the proximity of the C-terminal residue C338 to residues located on the outer face of a cytoplasmic helical extension of the F helix with an apparent increase of distance upon photoactivation.

Influenza hemagglutinin is spring-loaded by a metastable native conformation

Enveloped viruses enter cells by protein-mediated membrane fusion. For influenza virus, membrane fusion is regulated by the conformational state of the hemagglutinin (HA) protein, which switches from a native (nonfusogenic) structure to a fusion-active (fusogenic) conformation when exposed to the acidic environment of the cellular endosome. Here we demonstrate that destabilization of HA at neutral pH, with either heat or the denaturant urea, triggers a conformational change that is biochemically indistinguishable from the change triggered by low pH. In each case, the conformational change is coincident with induction of membrane-fusion activity, providing strong evidence that the fusogenic structure is formed. These results indicate that the native structure of HA is trapped in a metastable state and that the fusogenic conformation is released by destabilization of native structure. This strategy may be shared by other enveloped viruses, including those that enter the cell at neutral pH, and could have implications for understanding the membrane-fusion step of HIV infection.

Chaperone activity and structure of monomeric polypeptide binding domains of GroEL

The chaperonin GroEL is a large complex composed of 14 identical 57-kDa subunits that requires ATP and GroES for some of its activities. We find that a monomeric polypeptide corresponding to residues 191 to 345 has the activity of the tetradecamer both in facilitating the refolding of rhodanese and cyclophilin A in the absence of ATP and in catalyzing the unfolding of native barnase. Its crystal structure, solved at 2.5 Å resolution, shows a well-ordered domain with the same fold as in intact GroEL. We have thus isolated the active site of the complex allosteric molecular chaperone, which functions as a “minichaperone.” This has mechanistic implications: the presence of a central cavity in the GroEL complex is not essential for those representative activities in vitro, and neither are the allosteric properties. The function of the allosteric behavior on the binding of GroES and ATP must be to regulate the affinity of the protein for its various substrates in vivo, where the cavity may also be required for special functions.

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.

Substitutions in the aspartate transcarbamoylase domain of hamster CAD disrupt oligomeric structure

Aspartate transcarbamoylase (ATCase; EC 2.1.3.2) is one of three enzymatic domains of CAD, a protein whose native structure is usually a hexamer of identical subunits. Alanine substitutions for the ATCase residues Asp-90 and Arg-269 were generated in a bicistronic vector that encodes a 6-histidine-tagged hamster CAD. Stably transfected mammalian cells expressing high levels of CAD were easily isolated and CAD purification was simplified over previous procedures. The substitutions reduce the ATCase Vmax of the altered CADs by 11-fold and 46-fold, respectively, as well as affect the enzyme's affinity for aspartate. At 25 mM Mg2+, these substitutions cause the oligomeric CAD to dissociate into monomers. Under the same dissociating conditions, incubating the altered CAD with the ATCase substrate carbamoyl phosphate or the bisubstrate analogue N-phosphonacetyl-l-aspartate unexpectedly leads to the reformation of hexamers. Incubation with the other ATCase substrate, aspartate, has no effect. These results demonstrate that the ATCase domain is central to hexamer formation in CAD and suggest that the ATCase reaction mechanism is ordered in the same manner as the Escherichia coli ATCase. Finally, the data indicate that the binding of carbamoyl phosphate induces conformational changes that enhance the interaction of CAD subunits.

Pulling a single chromatin fiber reveals the forces that maintain its higher-order structure

Single chicken erythrocyte chromatin fibers were stretched and released at room temperature with force-measuring laser tweezers. In low ionic strength, the stretch-release curves reveal a process of continuous deformation with little or no internucleosomal attraction. A persistence length of 30 nm and a stretch modulus of ≈5 pN is determined for the fibers. At forces of 20 pN and higher, the fibers are modified irreversibly, probably through the mechanical removal of the histone cores from native chromatin. In 40–150 mM NaCl, a distinctive condensation-decondensation transition appears between 5 and 6 pN, corresponding to an internucleosomal attraction energy of ≈2.0 kcal/mol per nucleosome. Thus, in physiological ionic strength the fibers possess a dynamic structure in which the fiber locally interconverting between “open” and “closed” states because of thermal fluctuations.

Native display of complete foreign protein domains on the surface of hepatitis B virus capsids

The nucleocapsid of hepatitis B virus (HBV), or HBcAg, is a highly symmetric structure formed by multiple dimers of a single core protein that contains potent T helper epitopes in its 183-aa sequence. Both factors make HBcAg an unusually strong immunogen and an attractive candidate as a carrier for foreign epitopes. The immunodominant c/e1 epitope on the capsid has been suggested as a superior location to convey high immunogenicity to a heterologous sequence. Because of its central position, however, any c/e1 insert disrupts the core protein’s primary sequence; hence, only peptides, or rather small protein fragments seemed to be compatible with particle formation. According to recent structural data, the epitope is located at the tips of prominent surface spikes formed by the very stable dimer interfaces. We therefore reasoned that much larger inserts might be tolerated, provided the individual parts of a corresponding fusion protein could fold independently. Using the green fluorescent protein (GFP) as a model insert, we show that the chimeric protein efficiently forms fluorescent particles; hence, all of its structurally important parts must be properly folded. We also demonstrate that the GFP domains are surface-exposed and that the chimeric particles elicit a potent humoral response against native GFP. Hence, proteins of at least up to 238 aa can be natively displayed on the surface of HBV core particles. Such chimeras may not only be useful as vaccines but may also open the way for high resolution structural analyses of nonassembling proteins by electron microscopy.

Structure and function in rhodopsin: Further elucidation of the role of the intradiscal cysteines, Cys-110, -185, and -187, in rhodopsin folding and function

The disulfide bond between Cys-110 and Cys-187 in the intradiscal domain is required for correct folding in vivo and function of mammalian rhodopsin. Misfolding in rhodopsin, characterized by the loss of ability to bind 11-cis-retinal, has been shown to be caused by an intradiscal disulfide bond different from the above native disulfide bond. Further, naturally occurring single mutations of the intradiscal cysteines (C110F, C110Y, and C187Y) are associated with retinitis pigmentosa (RP). To elucidate further the role of every one of the three intradiscal cysteines, mutants containing single-cysteine replacements by alanine residues and the above three RP mutants have been studied. We find that C110A, C110F, and C110Y all form a disulfide bond between C185 and C187 and cause loss of retinal binding. C185A allows the formation of a C110–C187 disulfide bond, with wild-type-like rhodopsin phenotype. C187A forms a disulfide bond between C110 and C185 and binds retinal, and the pigment formed has markedly altered bleaching behavior. However, the opsin from the RP mutant C187Y forms no rhodopsin chromophore.

NMR investigation of ferricytochrome c unfolding: Detection of an equilibrium unfolding intermediate and residual structure in the denatured state

Horse ferricytochrome c (cyt c) undergoes exchange of one of its axial heme ligands (Met-80) for one or more non-native ligands under denaturing conditions. We have used 1H NMR spectroscopy to detect two conformations of paramagnetic cyt c with non-native heme ligation through a range of urea concentrations. One non-native form is an equilibrium unfolding intermediate observed under partially denaturing conditions and is attributed to replacement of Met-80 with one or more Lys side chains. The second non-native form, in which the native Met ligand is replaced by a His, is observed under strongly denaturing conditions. Thermodynamic analysis of these data indicates a relatively small ΔG (17 kJ/mol) for the transition from native to the Lys-ligated intermediate and a significantly larger ΔG (47 kJ/mol) for the transition from native to the His-ligated species. Although CD and fluorescence data indicate that the equilibrium unfolding of cyt c is a two-state process, these NMR results implicate an intermediate with His-Lys ligation.

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.

Statistical significance of protein structure prediction by threading

In this study, we estimate the statistical significance of structure prediction by threading. We introduce a single parameter ɛ that serves as a universal measure determining the probability that the best alignment is indeed a native-like analog. Parameter ɛ takes into account both length and composition of the query sequence and the number of decoys in threading simulation. It can be computed directly from the query sequence and potential of interactions, eliminating the need for sequence reshuffling and realignment. Although our theoretical analysis is general, here we compare its predictions with the results of gapless threading. Finally we estimate the number of decoys from which the native structure can be found by existing potentials of interactions. We discuss how this analysis can be extended to determine the optimal gap penalties for any sequence-structure alignment (threading) method, thus optimizing it to maximum possible performance.