Combining thermostable mutations increases the stability of lambda repressor

We have combined three mutations previously shown to stabilize lambda repressor against thermal denaturation. Two of these mutations are in helix 3, where Gly-46 and Gly-48 have been replaced by alanines [Hecht, M. H., et al. (1986) Proteins: Struct., Funct., Genet. 1, 43-46]. The other mutation, which replaces Tyr-88 with cysteine, allows the protein to form an intersubunit disulfide bond [Sauer, R. T., et al. (1986) Biochemistry 25, 5992-5998]. Calorimetric measurements show that the two alanine substitutions stabilize repressor by about 8 degrees C, that the disulfide bond stabilizes repressor by about 8 degrees C, and that the triple mutant is 16 degrees C more stable than wild-type repressor.     

The conformational and dynamic basis for ligand binding reactivity in hemoglobin Ypsilanti (beta 99 asp-->Tyr): origin of the quaternary enhancement effect

Hemoglobin Ypsilanti (HbY) is a stable tetrameric hemoglobin that binds oxygen with little or no cooperativity and with high affinity [Doyle, M. L., et al. (1992) Proteins: Struct., Funct., Genet. 14, 351-362]. It displays an especially large quaternary enhancement effect. An X-ray crystallographic study [Smith, F. R., et al. (1991) Proteins: Struct., Funct., Genet. 10, 81-91] of the carboxy derivative of this hemoglobin (COHbY) revealed a new quaternary structure that partially resembles the recently described R2 structure [Silva, M. M., et al. (1992) J. Biol. Chem. 267, 17248-17256]. Very little is known about either the solution phase conformations of the liganded and deoxy forms of HbY or the molecular basis for the large quaternary enhancement effect (Doyle et al., 1992). In this study, near-IR absorption, Soret-enhanced Raman, and UV (229 nm) resonance Raman spectroscopies are used to probe the liganded and deoxy derivatives of HbY in solution. Nanosecond time-resolved near-IR absorption measurements are used to expose the relaxation properties of the photoproduct of COHbY. Time-resolved (Soret band) absorption is used to generate the geminate and solvent phase ligand rebinding curves for photodissociated COHbY. The spectroscopic results indicate that COHbY has an R-like conformation with respect to both the proximal heme pocket and the hinge region of the alpha 1 beta 2 interface. The deoxy derivative of HbY has spectroscopic features that are very similar to those observed for species assigned to the deoxy R or half-liganded R conformations of human adult hemoglobin (HbA). The 10 ns to 100 micros relaxation properties of the photoproduct of COHbY are distinctly different from those of HbA in that for HbY, little if any tertiary or quaternary relaxation is observed. The near-absence of relaxation in the HbY photoproduct explains the differences in the geminate and solvent phase CO recombination between HbA and HbY. The impact of the conformational and relaxation properties of HbY on the geminate rebinding process forms the basis of a model that accounts for the large quaternary enhancement effect reported for HbY (Doyle et al., 1992). In addition, the spectroscopic data and the X-ray crystallographic results explain the slow relaxation for HbY and the near-absence of cooperative ligand binding for this protein based on the behavior of the penultimate tyrosines.     

The heme environment of mouse neuroglobin: histidine imidazole plane orientations obtained from solution NMR and EPR spectroscopy as compared with X-ray crystallography

The ¹H NMR chemical shifts of the heme methyl groups of the ferriheme complex of metneuroglobin (Du et al. in J. Am. Chem. Soc. 125:8080–8081, 2003) predict orientations of the axial histidine ligands (Shokhirev and Walker in J. Biol. Inorg. Chem. 3:581–594, 1998) that are not consistent with the X-ray data (Vallone et al. in Proteins Struct. Funct. Bioinf. 56:85–94, 2004), and the EPR spectrum (Vinck et al. in J. Am. Chem. Soc. 126:4516–4517, 2004) is only marginally consistent with these data. The reasons for these inconsistencies appear to be rooted in the high degree of aqueous solution exposure of the heme group and the fact that there are no strong hydrogen-bond acceptors for the histidine imidazole N–H protons provided by the protein. Similar inconsistencies may exist for other water-soluble heme proteins, and ¹H NMR spectroscopy provides a simple means to verify whether the solution structure of the heme center is the same as or different from that in the crystalline state.     

Structural mechanisms of constitutive activation in the C5a receptors with mutations in the extracellular loops: molecular modeling study

Previously we demonstrated by random saturation mutagenesis a set of mutations in the extracellular (EC) loops that constitutively activate the C5a receptor (C5aR) (Klco et al., Nat Struct Mol Biol 2005;12:320-326; Klco et al., J Biol Chem 2006;281:12010-12019). In this study, molecular modeling revealed possible conformations for the extracellular loops of the C5a receptors with mutations in the EC2 loop or in the EC3 loop. Comparison of low-energy conformations of the EC loops defined two distinct clusters of conformations typical either for strongly constitutively active mutants of C5aR (the CAM cluster) or for nonconstitutively active mutants (the non-CAM cluster). In the CAM cluster, the EC3 loop was turned towards the transmembrane (TM) helical bundle and more closely interacted with EC2 than in the non-CAM cluster. This suggested a structural mechanism of constitutive activity where EC3 contacts EC2 leading to EC2 interactions with helix TM3, thus triggering movement of TM7 towards TM2 and TM3. The movement initiates rearrangement of the system of hydrogen bonds between TM2, TM3 and TM7 including formation of the hydrogen bond between the side chains of D82(2.50) in TM2 and N296(7.49) in TM7, which is crucial for formation of the activated states of the C5a receptors (Nikiforovich et al., Proteins: Struct Funct Gene 2011;79:787-802). Since the relative large length of EC3 in C5aR (13 residues) is comparable with those in many other members of rhodopsin family of GPCRs (13-19 residues), our findings might reflect general mechanisms of receptor constitutive activation. The very recent X-ray structure of the agonist-induced constitutively active mutant of rhodopsin (Standfuss et al., Nature 2011;471:656-660) is discussed in view of our modeling results.     

Free energy landscape of a biomolecule in dihedral principal component space: sampling convergence and correspondence between structures and minima

Dihedral principal component analysis (dPCA) has recently been developed and shown to display complex features of the free energy landscape of a biomolecule that may be absent in the free energy landscape plotted in principal component space due to mixing of internal and overall rotational motion that can occur in principal component analysis (PCA) [Mu et al., Proteins: Struct Funct Bioinfo 2005;58:45-52]. Another difficulty in the implementation of PCA is sampling convergence, which we address here for both dPCA and PCA using a tetrapeptide as an example. We find that for both methods the sampling convergence can be reached over a similar time. Minima in the free energy landscape in the space of the two largest dihedral principal components often correspond to unique structures, though we also find some distinct minima to correspond to the same structure.     

Nearest neighbour algorithm for predicting protein subcellular location by combining functional domain composition and pseudo-amino acid composition

In this paper, based on the approach by combining the "functional domain composition" [K.C. Chou, Y. D. Cai, J. Biol. Chem. 277 (2002) 45765] and the pseudo-amino acid composition [K.C. Chou, Proteins Struct. Funct. Genet. 43 (2001) 246; Correction Proteins Struct. Funct. Genet. 2044 (2001) 2060], the Nearest Neighbour Algorithm (NNA) was developed for predicting the protein subcellular location. Very high success rates were observed, suggesting that such a hybrid approach may become a useful high-throughput tool in the area of bioinformatics and proteomics.     

Transcellular transport of fluorescein dextran through an arterial endothelial cell monolayer

Transcellular transport of fluorescein dextran (FD) of various molecular weights (4K, 10K, 20K, 70K and 150K daltons) through porcine arterial endothelial cells cultured on a type I collagen gel supported by a dacron sheet was studied and compared with the transport of low density lipoprotein labeled with rhodamine B (RB-LDL) described previously (Hashida et al., Cell Struct. Funct. 11, 31-42, 1986). The rate of FD transport through the monolayer depended on the size of the FD. FD transport was not temperature-dependent and was not a saturable process. Our findings show that FD transport differs from RB-LDL transport which is temperature- and dose-dependent. The mechanism of the transport of FD is compared with that of RB-LDL.     

A monoclonal antibody against human urokinase: the epitope structure and sequence homology with a human tissue-type plasminogen activator

Our previous study showed that an epitope defined by a monoclonal antibody against human urokinase is located on the 33-Kdalton catalytic domain of the enzyme (Nakamura, M. et al., Cell Struct Funct., 9, 167-179, 1984). The epitope structure was further determined and characterized on one-dimensional SDS-polyacrylamide slab gel maps of CNBr-cleaved polypeptide fragments as well as on their Western blots. A single homogeneous polypeptide with an approximate molecular weight of 3.4-Kdaltons was found to be antigenic. The monoclonal antibody exhibited a stronger inhibition of the enzyme activity than the polyclonal antibodies tested, and cross-reacted with a 65-Kdalton tissue-type plasminogen activator present in Detroit 562 cells. From these results and data made up with the help of a computer comparison of known sequences of urokinase and a tissue-type plasminogen activator, we concluded that the epitope is Cys-Gln-Gly-Asp-Ser-Gly-Gly-Pro-Leu-Val-Cys and contains a catalytically active residue, serine.     

Quenching of fluorescein-conjugated lipids by antibodies. Quantitative recognition and binding of lipid-bound haptens in biomembrane models, formation of two-dimensional protein domains and molecular dynamics simulations.

Three model biomembrane systems, monolayers, micelles, and vesicles, have been used to study the influence of chemical and physical variables of hapten presentation at membrane interfaces on antibody binding. Hapten recognition and binding were monitored for the anti-fluorescein monoclonal antibody 4-4-20 generated against the hapten, fluorescein, in these membrane models as a function of fluorescein-conjugated lipid architecture. Specific recognition and binding in this system are conveniently monitored by quenching of fluorescein emission upon penetration of fluorescein into the antibody's active site. Lipid structure was shown to play a large role in affecting antibody quenching. Interestingly, the observed degrees of quenching were nearly independent of the lipid membrane model studied, but directly correlated with the chemical structure of the lipids. In all cases, the antibody recognized and quenched most efficiently a lipid based on dioctadecylamine where fluorescein is attached to the headgroup via a long, flexible hydrophilic spacer. Dipalmitoyl phosphatidylethanolamine containing a fluorescein headgroup demonstrated only partial binding/quenching. Egg phosphatidylethanolamine with a fluorescein headgroup showed no susceptibility to antibody recognition, binding, or quenching. Formation of two-dimensional protein domains upon antibody binding to the fluorescein-lipids in monolayers is also presented. Chemical and physical requirements for these antibody-hapten complexes at membrane surfaces have been discussed in terms of molecular dynamics simulations based on recent crystallographic models for this antibody-hapten complex (Herron et al., 1989. Proteins Struct. Funct. Genet. 5:271-280).     

Structure from NMR and molecular dynamics: Distance restraining inhibits motion in the essential subspace

Summary We address the question how well proteins can be modelled on the basis of NMR data, when these data are incorporated into the protein model using distance restraints in a molecular dynamics simulation. We found, using HPr as a model protein, that distance restraining freezes the essential motion of proteins, as defined by Amadei et al. [Amadei, A., Linssen, A.B.M. and Berendsen, H.J.C. (1993) Protein Struct. Funct. Genet., 17, 412–425]. We discuss how modelling protocols can be improved in order to solve this problem.     

Deletion of the omega-loop in the active site of staphylococcal nuclease. 1. Effect on catalysis and stability

The high-resolution X-ray structure of wild-type staphylococcal nuclease (E43 SNase) suggests that Glu 43 acts a general basic catalyst to assist the attack of water on a phosphodiester substrate [Loll, P., & Lattman, E. E. (1989) Proteins: Struct., Funct., Genet. 5, 183]. Glu 43 is located at the base of the solvent-exposed and conformationally mobile omega-loop in the active site of E43 SNase having the sequence Glu43-Thr44-Lys45-His46-Pro47-Lys48- Lys49-Gly50-Val51-Glu52, where the gamma-carboxylate of Glu 52 is hydrogen bonded to the amide hydrogen of Glu 43. With a metabolic selection for SNase activity produced in an Escherichia coli host, we detected an unexpected deletion of residues 44-49 of the omega-loop of E43 SNase in cassette mutagenesis experiments designed to randomize codons 44 and 45 in the omega-loop and increase the activity of the previously described E43D mutation (D43 SNase). A high-resolution X-ray structure of D43 SNase has revealed that the E43D substitution significantly changes the structure of the omega-loop, reduces the interaction of the essential Ca2+ ion with its active-site ligands, and diminishes the network of hydrogen-bonded water molecules in the active site [Loll, P., & Lattman, E. E. (1990) Biochemistry 29, 6866]. This deletion of six amino acids from the omega-loop generates a protein (E43 delta SNase) having a partially solvent-exposed, surface beta-turn with the sequence Glu43-Gly50-Val51-Glu52; the structure of this beta-turn is addressed in the following article [Baldisseri et al. (1991) Biochemistry (following paper in this issue)].(ABSTRACT TRUNCATED AT 250 WORDS)     

Using functional domain composition and support vector machines for prediction of protein subcellular location

Proteins are generally classified into the following 12 subcellular locations: 1) chloroplast, 2) cytoplasm, 3) cytoskeleton, 4) endoplasmic reticulum, 5) extracellular, 6) Golgi apparatus, 7) lysosome, 8) mitochondria, 9) nucleus, 10) peroxisome, 11) plasma membrane, and 12) vacuole. Because the function of a protein is closely correlated with its subcellular location, with the rapid increase in new protein sequences entering into databanks, it is vitally important for both basic research and pharmaceutical industry to establish a high throughput tool for predicting protein subcellular location. In this paper, a new concept, the so-called "functional domain composition" is introduced. Based on the novel concept, the representation for a protein can be defined as a vector in a high-dimensional space, where each of the clustered functional domains derived from the protein universe serves as a vector base. With such a novel representation for a protein, the support vector machine (SVM) algorithm is introduced for predicting protein subcellular location. High success rates are obtained by the self-consistency test, jackknife test, and independent dataset test, respectively. The current approach not only can play an important complementary role to the powerful covariant discriminant algorithm based on the pseudo amino acid composition representation (Chou, K. C. (2001) Proteins Struct. Funct. Genet. 43, 246-255; Correction (2001) Proteins Struct. Funct. Genet. 44, 60), but also may greatly stimulate the development of this area.     

Inhibitory effect of cholesterol on interaction between cytoplasmic actin and liposomes, and restorative effect of high osmotic pressure

Although cholesterol is one of the major components of plasma membranes in eukaryotic cells, very little is known about its role in biological membranes. We reported previously (Okimasu et al., Cell Struct. Funct. 11, 273-283, 1986) that introduction of cholesterol into the liposomal membrane caused a decrease in membrane permeability, especially by the binding of cytoplasmic proteins to the liposomal membrane. The present study was carried out to further clarify the biochemical function of cholesterol in the membrane-protein interactions, especially under high osmotic pressure. The association of membranes with cytoplasmic proteins and their permeability were decreased by the introduction of cholesterol, but its effects were diminished in a hypertonic medium. The protein species associated with cholesterol-containing liposomes vary depending on the sort of hypertonic condition. It was suggested that since the degree of lipid packing by the cholesterol was reduced by the locally increased curvature in the lipid bilayer under high osmotic pressure, some cytoplasmic proteins can penetrate into the liposomal membrane.     

A common evolutionary origin of two elementary enzyme folds

The (beta alpha)(8)-barrel is the most frequent and most versatile fold among enzymes [H?cker et al., Curr. Opin. Biotechnol. 12 (2001) 376-381; Wierenga, FEBS Lett. 492 (2001) 193-198]. Structural and functional evidence suggests that (beta alpha)(8)-barrels evolved from an ancestral half-barrel, which consisted of four (beta alpha) units stabilized by dimerization [Lang et al., Science 289 (2000) 1546-550; H?cker et al., Nat. Struct. Biol. 8 (2001) 32-36; Gerlt and Babbitt, Nat. Struct. Biol. 8 (2001) 5-7]. Here, by performing a comprehensive database search, we detect a striking and unexpected structural and amino acid sequence similarity between (beta alpha)(4) half-barrels and members of the (beta alpha)(5) flavodoxin-like fold. These findings provoke the hypothesis that a large fraction of the modern-day enzymes evolved from a basic structural building block, which can be identified by a combination of sequence and structural analyses.     

Predicting protein structural classes with pseudo amino acid composition: an approach using geometric moments of cellular automaton image

A novel approach was developed for predicting the structural classes of proteins based on their sequences. It was assumed that proteins belonging to the same structural class must bear some sort of similar texture on the images generated by the cellular automaton evolving rule [Wolfram, S., 1984. Cellular automation as models of complexity. Nature 311, 419-424]. Based on this, two geometric invariant moment factors derived from the image functions were used as the pseudo amino acid components [Chou, K.C., 2001. Prediction of protein cellular attributes using pseudo amino acid composition. Proteins: Struct., Funct., Genet. (Erratum: ibid., 2001, vol. 44, 60) 43, 246-255] to formulate the protein samples for statistical prediction. The success rates thus obtained on a previously constructed benchmark dataset are quite promising, implying that the cellular automaton image can help to reveal some inherent and subtle features deeply hidden in a pile of long and complicated amino acid sequences.     

Using optimized evidence-theoretic K-nearest neighbor classifier and pseudo-amino acid composition to predict membrane protein types

Knowledge of membrane protein type often provides crucial hints toward determining the function of an uncharacterized membrane protein. With the avalanche of new protein sequences emerging during the post-genomic era, it is highly desirable to develop an automated method that can serve as a high throughput tool in identifying the types of newly found membrane proteins according to their primary sequences, so as to timely make the relevant annotations on them for the reference usage in both basic research and drug discovery. Based on the concept of pseudo-amino acid composition [K.C. Chou, Proteins: Struct. Funct. Genet. 43 (2001) 246-255; Erratum: Proteins: Struct. Funct. Genet. 44 (2001) 60] that has made it possible to incorporate a considerable amount of sequence-order effects by representing a protein sample in terms of a set of discrete numbers, a novel predictor, the so-called "optimized evidence-theoretic K-nearest neighbor" or "OET-KNN" classifier, was proposed. It was demonstrated via the self-consistency test, jackknife test, and independent dataset test that the new predictor, compared with many previous ones, yielded higher success rates in most cases. The new predictor can also be used to improve the prediction quality for, among many other protein attributes, structural class, subcellular localization, enzyme family class, and G-protein coupled receptor type. The OET-KNN classifier will be available as a web-server at http://www.pami.sjtu.edu.cn/kcchou.     

Analysis of the differences in the folding kinetics of structurally homologous proteins based on predictions of the gross features of residue contacts

It is a general notion that proteins with very similar three-dimensional structures would show very similar folding kinetics. However, recent studies reveal that the folding kinetic properties of some proteins contradict this thought (i.e., the members in a same protein family fold through different pathways). For example, it has been reported that some beta-proteins in the intracellular lipid-binding protein family fold through quite different pathways (Burns et al., Proteins 1998;33:107-118). Similar differences in folding kinetics are also observed in the members of the globin family (Nishimura et al., Nat Struct Biol 2000;7:679-686). In our study, we examine the possibility of predicting qualitative differences in folding kinetics of the intracellular lipid-binding proteins and two globin proteins (i.e., myoglobin and leghemoglobin). The problem is tackled by means of a contact map based on the average distance statistics between residues, the Average Distance Map (ADM), as constructed from sequence. The ADMs for the three proteins show overall similarity, but some local differences among maps are also observed. Our results demonstrate that some properties of the protein folding kinetics are consistent with local differences in the ADMs. We also discuss the general possibility of predicting folding kinetics from sequence information.     

Folding-unfolding equilibrium and kinetics of equine beta-lactoglobulin: equivalence between the equilibrium molten globule state and a burst-phase folding intermediate

The denaturant-induced equilibrium unfolding transition of equine beta-lactoglobulin was investigated by ultraviolet absorption, fluorescence, and circular dichroism (CD) spectra. An equilibrium intermediate populates at moderate denaturant concentrations, and its CD spectrum is similar to that of the molten globule state previously observed for this protein at acid pH [Ikeguchi, M., Kato, S., Shimizu, A., and Sugai, S. (1997) Proteins: Struct., Funct., Genet. 27, 567-575]. The unfolding and refolding kinetics were also investigated by the stopped-flow CD and fluorescence. A significant change in the CD intensity was observed within the dead time of measurements (25 ms) when the refolding reaction was initiated by diluting the urea-unfolded protein solution, indicating the transient accumulation of the folding intermediate. The CD spectrum of this burst-phase intermediate agrees well with that of the molten globule state at acid pH. The stability of the burst-phase intermediate was also estimated from the urea-concentration dependence of the burst-phase amplitude, and it shows a fair agreement with that of the equilibrium intermediate. These results indicate that the molten globule state of equine beta-lactoglobulin populates at moderate urea concentration as well as at acid pH and it is equivalent with the kinetic folding intermediate.