Search

Detection of homologous relationships among proteins and understanding their mechanisms of diversification are major topics in the fields of protein science, bioinformatics, and phylogenetics. Recent developments in sequence/profile-based and structural similarity-based methods have greatly facilitated the unification and classification of many protein families into superfamilies or folds, yet many proteins remain unclassified in current protein databases. As one of the three earliest identified RNases in biology, ribonuclease T2, also known as RNase I in Escherichia coli, RNase Rh in fungi, or S-RNase in plant, is thought to be an ancient RNase family due to its widespread distribution and distinct structure. In this study, we present evidence that RNase T2 represents a circularly permutated version of the BECR (Barnase-EndoU-Colicin E5/D-RelE) fold RNases. This subtle relationship cannot be detected by traditional methods such as sequence/profile-based comparisons, structure-similarity searches, and circular permutation detections. However, we were able to identify the structural similarity using rational reconstruction of a theoretical RNase T2 ancestor via a reverse circular permutation process, followed by structural modeling using AlphaFold2, and structural comparisons. This relationship is further supported by the fact that RNase T2 and other typical BECR RNases, namely Colicin D, RNase A, and BrnT, share similar catalytic site configurations, all involving an analogous set of conserved residues on the α0 helix and the β4 strand of the BECR fold. This study revealed a hidden root of RNase T2 in bacterial toxin systems and demonstrated that reconstruction and modeling of ancestral topology is an effective strategy to identify remote relationship between proteins., (© 2022 The Protein Society.)

Regulated degradation of mature, cytoplasmic mRNA is a key step in eukaryotic gene regulation. This process is typically initiated by the recruitment of deadenylase enzymes by cis-acting elements in the 3' untranslated region resulting in the shortening and removal of the 3' poly(A) tail of the target mRNA. The Ccr4-Not complex, a major eukaryotic deadenylase, contains two exoribonuclease subunits with selectivity toward poly(A): Caf1 and Ccr4. The Caf1 deadenylase subunit binds the MIF4G domain of the large subunit CNOT1 (Not1) that is the scaffold of the complex. The Ccr4 nuclease is connected to the complex via its leucine-rich repeat (LRR) domain, which binds Caf1, whereas the catalytic activity of Ccr4 is provided by its EEP domain. While the relative positions of the MIF4G domain of CNOT1, the Caf1 subunit, and the LRR domain of Ccr4 are clearly defined in current models, the position of the EEP nuclease domain of Ccr4 is ambiguous. Here, we use X-ray crystallography, the AlphaFold resource of predicted protein structures, and pulse electron paramagnetic resonance spectroscopy to determine and validate the position of the EEP nuclease domain of Ccr4 resulting in an improved model of the human Ccr4-Not nuclease module., (© 2021 The Protein Society.)

PhoH2 proteins are found in a very diverse range of microorganisms that span bacteria and archaea. These proteins are composed of two domains: an N-terminal PIN-domain fused with a C-terminal PhoH domain. Collectively this fusion functions as an RNA helicase and ribonuclease. In other genomic contexts, PINdomains and PhoHdomains are separate but adjacent suggesting association to achieve similar function. Exclusively among the mycobacteria, PhoH2 proteins are encoded in the genome with an upstream gene, phoAT, which is thought to play the role of an antitoxin (in place of the traditional VapB antitoxin that lies upstream of the 47 other PINdomains in the mycobacterial genome). This review examines PhoH2 proteins as a whole and describes the bioinformatics, biochemical, structural, and biological properties of the two domains that make up PhoH2: PIN and PhoH. We review the transcriptional regulators of phoH2 from two mycobacterial species and speculate on the function of PhoH2 proteins in the context of a Type II toxin-antitoxin system which are thought to play a role in the stress response in bacteria., (© 2019 The Protein Society.)

Antibiotic resistance and microorganism virulence have been consistently exhibited by bacteria and archaea, which survive in conditions of environmental stress through toxin-antitoxin (TA) systems. The HP0892-HP0893 TA system is one of the two known TA systems belonging to Helicobacter pylori. The antitoxin, HP0893, binds and inhibits the HP0892 toxin and regulates the transcription of the TA operon. Here, we present the crystal structure of the zinc-bound HP0892 toxin at 1.8 Å resolution. Reorientation of residues at the mRNase active site was shown. The involved residues, namely E58A, H86A, and H58A/ H60A, were mutated and the binding affinity was monitored by ITC studies. Through the structural difference between the apo and the metal-bound state, and using a homology modeling tool, the involvement of the metal ion in mRNase active site could be identified. The most catalytically important residue, His86, reorients itself to exhibit RNase activity. His47, Glu58, and His60 are involved in metal binding where Glu58 acts as a general base and His47 and His60 may also act as a general acid in enzymatic activity. Glu58 and Asp64 are involved in substrate binding and specific sequence recognition. Arg83 is involved in phosphate binding and stabilization of the transition state, and Phe90 is involved in base packing and substrate orientation., (© 2014 The Protein Society.)

Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein., (© 2014 The Protein Society.)

Recently, we showed that leucine zipper (LZ) motifs of basic leucine zipper (bZIP) transcription factors GCN4 and c-Jun are capable of catalyzing degradation of RNA (Nikolaev et al., PLoS ONE 2010; 5:e10765). This observation is intriguing given the tight regulation of RNA turnover control and the antiquity of bZIP transcription factors. To support further mechanistic studies, herein, we elucidated RNA binding interface of the GCN4 leucine zipper motif from yeast. Solution NMR experiments showed that the LZ-RNA interaction interface is located in the first two heptads of LZ moiety, and that only the dimeric (coiled coil) LZ conformation is capable of binding RNA. Site-directed mutagenesis of the LZ-GCN4 RNA binding interface showed that substrate binding is facilitated by lysine and arginine side chains, and that at least one nucleophilic residue is located in proximity to the RNA phosphate backbone. Further studies in the context of full-length bZIP factors are envisaged to address the biological relevance of LZ RNase activity., (Copyright © 2012 The Protein Society.)

Protein topology defined by the matrix of residue contacts has proved to be a fruitful basis for the study of protein dynamics. The widely implemented coarse-grained elastic network model of backbone fluctuations has been used to describe crystallographic temperature factors, allosteric couplings, and some aspects of the folding pathway. In the present study, we develop a model of protein dynamics based on the classical equations of motion of a damped network model (DNM) that describes the folding path from a completely unfolded state to the native conformation through a single-well potential derived purely from the native conformation. The kinetic energy gained through the collapse of the protein chain is dissipated through a friction term in the equations of motion that models the water bath. This approach is completely general and sufficiently fast that it can be applied to large proteins. Folding pathways for various proteins of different classes are described and shown to correlate with experimental observations and molecular dynamics and Monte Carlo simulations. Allosteric transitions between alternative protein structures are also modeled within the DNM through an asymmetric double-well potential., (Copyright © 2010 The Protein Society.)

Increasing the conformational stability of proteins is an important goal for both basic research and industrial applications. In vitro selection has been used successfully to increase protein stability, but more often site-directed mutagenesis is used to optimize the various forces that contribute to protein stability. In previous studies, we showed that improving electrostatic interactions on the protein surface and improving the beta-turn sequences were good general strategies for increasing protein stability, and used them to increase the stability of RNase Sa. By incorporating seven of these mutations in RNase Sa, we increased the stability by 5.3 kcal/mol. Adding one more mutation, D79F, gave a total increase in stability of 7.7 kcal/mol, and a melting temperature 28 degrees C higher than the wild-type enzyme. Surprisingly, the D79F mutation lowers the change in heat capacity for folding, DeltaC(p), by 0.6 kcal/mol/K. This suggests that this mutation stabilizes structure in the denatured state ensemble. We made other mutants that give some insight into the structure present in the denatured state. Finally, the thermodynamics of folding of these stabilized variants of RNase Sa are compared with those observed for proteins from thermophiles.

The advantages of electrospray ionization mass spectrometry (ESIMS) to measure relative solution-phase affinities of tightly bound protein-protein complexes are demonstrated with selected variants of the Bacillus amyloliquefaciens protein barstar (b*) and the RNAase barnase (bn), which form protein-protein complexes with a range of picomolar to nanomolar dissociation constants. A novel chemical annealing procedure rapidly establishes equilibrium in solutions containing competing b* variants with limiting bn. The relative ion abundances of the complexes and those of the competing unbound monomers are shown to reflect the relative solution-phase concentrations of those respective species. No measurable dissociation of the complexes occurs either during ESI or mass detection, nor is there any evidence for nonspecific binding at protein concentrations < 25 microM. Differences in DeltaDeltaG of dissociation between variants were determined with precisions < 0.1 kcal/mol. The DeltaDeltaG values obtained deviate on average by 0.26 kcal/mol from those measured with a solution-phase enzyme assay. It is demonstrated that information about the protein conformation and covalent modifications can be obtained from differences in mass and charge state distributions. This method serves as a rapid and precise means to interrogate protein-protein-binding surfaces for complexes that have affinities in the picomolar to nanomolar range.

Nuclease type colicins and related bacteriocins possess the unprecedented ability to translocate an enzymatic polypeptide chain across the Gram-negative cell envelope. Here we use the rRNase domain of the cytotoxic ribonuclease colicin E3 to examine the structural changes on its interaction with the membrane. Using phospholipid vesicles as model membranes we show that anionic membranes destabilize the nuclease domain of the rRNase type colicin E3. Intrinsic tryptophan fluorescence and circular dichroism show that vesicles consisting of pure DOPA act as a powerful protein denaturant toward the rRNase domain, although this interaction can be entirely prevented by the addition of salt. Binding of E3 rRNase to DOPA vesicles is an endothermic process (DeltaH=24 kcal mol-1), reflecting unfolding of the protein. Consistent with this, binding of a highly destabilized mutant of the E3 rRNase to DOPA vesicles is exothermic. With mixed vesicles containing anionic and neutral phospholipids at a ratio of 1:3, set to mimic the charge of the Escherichia coli inner membrane, destabilization of E3 rRNase is lessened, although the melting temperature of the protein at pH 7.0 is greatly reduced from 50 degrees C to 30 degrees C. The interaction of E3 rRNase with 1:3 DOPA:DOPC vesicles is also highly dependent on both ionic strength and temperature. We discuss these results in terms of the likely interaction of the E3 rRNase and the related E9 DNase domains with the E. coli inner membrane and their subsequent translocation to the cell cytoplasm.

Many proteins that accumulate in the form of insoluble aggregates when they are overproduced in Escherichia coli can be rendered soluble by fusing them to E. coli maltose binding protein (MBP), and this will often enable them to fold in to their biologically active conformations. Yet, although it is an excellent solubility enhancer, MBP is not a particularly good affinity tag for protein purification. To compensate for this shortcoming, we have engineered and successfully tested Gateway destination vectors for the production of dual His6MBP-tagged fusion proteins in the cytoplasm and periplasm of E. coli. The MBP moiety improves the yield and solubility of its fusion partners while the hexahistidine tag (His-tag) serves to facilitate their purification. The availability of a vector that targets His6MBP fusion proteins to the periplasm expands the utility of this dual tagging approach to include proteins that contain disulfide bonds or are toxic in the bacterial cytoplasm.

Gaining a better understanding of the denatured state ensemble of proteins is important for understanding protein stability and the mechanism of protein folding. We studied the folding kinetics of ribonuclease Sa (RNase Sa) and a charge-reversal variant (D17R). The refolding kinetics are similar, but the unfolding rate constant is 10-fold greater for the variant. This suggests that charge-charge interactions in the denatured state and the transition state ensembles are more favorable in the variant than in RNase Sa, and shows that charge-charge interactions can influence the kinetics and mechanism of protein folding.

The near-UV absorption spectra of barnase double-point mutants are calculated using a combination of molecular dynamics and ab initio techniques. The atoms of the fluorescent probes are placed in a cloud of point charges, generated by molecular dynamics simulations. Ab initio calculations (CASPT2) are performed on these systems. Three molecular dynamics packages are compared-Amber5.0, CHARMM-c27b1, and GROMOS96-using indole as the fluorescent probe. It was found that calculated absorption spectra reproduce experimental values very well, provided detailed charge cloud descriptions are included. These calculations further sustain the hypothesis that different tryptophan rotamers can be present in proteins. Molecular dynamics calculations of the double-point mutants also point to the structural effect of counter ions.

Colicin E3 is a cytotoxic ribonuclease that specifically cleaves 16S rRNA at the ribosomal A-site to abolish protein synthesis in sensitive Escherichia coli cells. We have performed extensive mutagenesis of the 96-residue colicin E3 cytotoxic domain (E3 rRNase), assayed mutant colicins for in vivo cytotoxicity, and tested the corresponding E3 rRNase domains for their ability to inactivate ribosome function in vitro. From 21 alanine mutants, we identified five positions where mutation resulted in a colicin with no measurable cytotoxicity (Y52, D55, H58, E62, and Y64) and four positions (R40, R42, E60, and R90) where mutation caused a significant reduction in cytotoxicity. Mutations that were found to have large in vivo and in vitro effects were tested for structural integrity through circular dichroism and fluorescence spectroscopy using purified rRNase domains. Our data indicate that H58 and E62 likely act as the acid-base pair during catalysis with other residues likely involved in transition state stabilization. Both the Y52 and Y64 mutants were found to be highly destabilized and this is the likely origin of the loss of their cytotoxicity. The identification of important active site residues and sequence alignments of known rRNase homologs has allowed us to identify other proteins containing the putative rRNase active site motif. Proteins that contained this active site motif included three hemagglutinin-type adhesins and we speculate that these have evolved to deliver a cytotoxic rRNase into eukaryotic cells during pathogenesis.

The crystal structures of bovine pancreatic ribonuclease A (RNase A) in complex with 3',5'-ADP, 2',5'-ADP, 5'-ADP, U-2'-p and U-3'-p have been determined at high resolution. The structures reveal that each inhibitor binds differently in the RNase A active site by anchoring a phosphate group in subsite P1. The most potent inhibitor of all five, 5'-ADP (Ki = 1.2 microM), adopts a syn conformation (in contrast to 3',5'-ADP and 2',5'-ADP, which adopt an anti), and it is the beta- rather than the alpha-phosphate group that binds to P1. 3',5'-ADP binds with the 5'-phosphate group in P1 and the adenosine in the B2 pocket. Two different binding modes are observed in the two RNase A molecules of the asymmetric unit for 2',5'-ADP. This inhibitor binds with either the 3' or the 5' phosphate groups in subsite P1, and in each case, the adenosine binds in two different positions within the B2 subsite. The two uridilyl inhibitors bind similarly with the uridine moiety in the B1 subsite but the placement of a different phosphate group in P1 (2' versus 3') has significant implications on their potency against RNase A. Comparative structural analysis of the RNase A, eosinophil-derived neurotoxin (EDN), eosinophil cationic protein (ECP), and human angiogenin (Ang) complexes with these and other phosphonucleotide inhibitors provides a wealth of information for structure-based design of inhibitors specific for each RNase. These inhibitors could be developed to therapeutic agents that could control the biological activities of EDN, ECP, and ANG, which play key roles in human pathologies.

Ribonuclease Sa and two charge-reversal variants can be converted into amyloid in vitro by the addition of 2,2,2-triflouroethanol (TFE). We report here amyloid fibril formation for these proteins as a function of pH. The pH at maximal fibril formation correlates with the pH dependence of protein solubility, but not with stability, for these variants. Additionally, we show that the pH at maximal fibril formation for a number of well-characterized proteins is near the pI, where the protein is expected to be the least soluble. This suggests that protein solubility is an important determinant of fibril formation.

We have used site-specific mutagenesis to study the contribution of Glu 74 and the active site residues Gln 38, Glu 41, Glu 54, Arg 65, and His 85 to the catalytic activity and thermal stability of ribonuclease Sa. The activity of Gln38Ala is lowered by one order of magnitude, which confirms the involvement of this residue in substrate binding. In contrast, Glu41Lys had no effect on the ribonuclease Sa activity. This is surprising, because the hydrogen bond between the guanosine N1 atom and the side chain of Glu 41 is thought to be important for the guanine specificity in related ribonucleases. The activities of Glu54Gln and Arg65Ala are both lowered about 1000-fold, and His85Gln is totally inactive, confirming the importance of these residues to the catalytic function of ribonuclease Sa. In Glu74Lys, k(cat) is reduced sixfold despite the fact that Glu 74 is over 15 A from the active site. The pH dependence of k(cat)/K(M) is very similar for Glu74Lys and wild-type RNase Sa, suggesting that this is not due to a change in the pK values of the groups involved in catalysis. Compared to wild-type RNase Sa, the stabilities of Gln38Ala and Glu74Lys are increased, the stabilities of Glu41Lys, Glu54Gln, and Arg65Ala are decreased and the stability of His85Gln is unchanged. Thus, the active site residues in the ribonuclease Sa make different contributions to the stability.

The rapid binding of cytotoxic colicin E3 by its cognate immunity protein Im3 is essential in safeguarding the producing cell. The X-ray structure of the E3/Im3 complex shows that the Im3 molecule interfaces with both the C-terminal ribonuclease (RNase) domain and the N-terminal translocation domain of E3. The association and dissociation rates of the RNase domain and Im3 show drastically different sensitivities to ionic strength, as previously rationalized for electrostatically enhanced diffusion-limited protein-protein associations. Relative to binding to the RNase domain, binding to full-length E3 shows a comparable association rate but a significantly lower dissociation rate. This outcome is just what was anticipated by a theory for the binding of two linked domains to a protein. The E3/Im3 system thus provides a powerful paradigm for the interplay of theory and experiment.

Osmolytes increase the thermodynamic conformational stability of proteins, shifting the equilibrium between native and denatured states to favor the native state. However, their effects on conformational equilibria within native-state ensembles of proteins remain controversial. We investigated the effects of sucrose, a model osmolyte, on conformational equilibria and fluctuations within the native-state ensembles of bovine pancreatic ribonuclease A and S and horse heart cytochrome c. In the presence of sucrose, the far- and near-UV circular dichroism spectra of all three native proteins were slightly altered and indicated that the sugar shifted the native-state ensemble toward species with more ordered, compact conformations, without detectable changes in secondary structural contents. Thermodynamic stability of the proteins, as measured by guanidine HCl-induced unfolding, increased in proportion to sucrose concentration. Native-state hydrogen exchange (HX) studies monitored by infrared spectroscopy showed that addition of 1 M sucrose reduced average HX rate constants at all degrees of exchange of the proteins, for which comparison could be made in the presence and absence of sucrose. Sucrose also increased the exchange-resistant core regions of the proteins. A coupling factor analysis relating the free energy of HX to the free energy of unfolding showed that sucrose had greater effects on large-scale than on small-scale fluctuations. These results indicate that the presence of sucrose shifts the conformational equilibria toward the most compact protein species within native-state ensembles, which can be explained by preferential exclusion of sucrose from the protein surface.

Ribonuclease Sa (pI = 3.5) from Streptomyces aureofaciens and its 3K (D1K, D17K, E41K) (pI = 6.4) and 5K (3K + D25K, E74K) (pI = 10.2) mutants were tested for cytotoxicity. The 5K mutant was cytotoxic to normal and v-ras-transformed NIH3T3 mouse fibroblasts, but RNase Sa and 3K were not. The structure, stability, and activity of the three proteins are comparable, but the net charge at pH 7 increases from -7 for RNase Sa to -1 for 3K and to +3 for 5K. These results suggest that a net positive charge is a key determinant of ribonuclease cytotoxicity. The cytotoxic 5K mutant preferentially attacks v-ras-NIH3T3 fibroblasts, suggesting that mammalian cells expressing the ras-oncogene are potential targets for ribonuclease-based drugs.

Van der Waals (vdW) interaction energies between different atom types, energies of hydrogen bonds (H-bonds), and atomic solvation parameters (ASPs) have been derived from the published thermodynamic stabilities of 106 mutants with available crystal structures by use of an originally designed model for the calculation of free-energy differences. The set of mutants included substitutions of uncharged, inflexible, water-inaccessible residues in alpha-helices and beta-sheets of T4, human, and hen lysozymes and HI ribonuclease. The determined energies of vdW interactions and H-bonds were smaller than in molecular mechanics and followed the "like dissolves like" rule, as expected in condensed media but not in vacuum. The depths of modified Lennard-Jones potentials were -0.34, -0.12, and -0.06 kcal/mole for similar atom types (polar-polar, aromatic-aromatic, and aliphatic-aliphatic interactions, respectively) and -0.10, -0.08, -0.06, -0.02, and nearly 0 kcal/mole for different types (sulfur-polar, sulfur-aromatic, sulfur-aliphatic, aliphatic-aromatic, and carbon-polar, respectively), whereas the depths of H-bond potentials were -1.5 to -1.8 kcal/mole. The obtained solvation parameters, that is, transfer energies from water to the protein interior, were 19, 7, -1, -21, and -66 cal/moleA(2) for aliphatic carbon, aromatic carbon, sulfur, nitrogen, and oxygen, respectively, which is close to the cyclohexane scale for aliphatic and aromatic groups but intermediate between octanol and cyclohexane for others. An analysis of additional replacements at the water-protein interface indicates that vdW interactions between protein atoms are reduced when they occur across water.

We present a novel and efficient approach for assessing protein-protein complex formation, which combines ab initio docking calculations performed with the protein docking algorithm BiGGER and chemical shift perturbation data collected with heteronuclear single quantum coherence (HSQC) or TROSY nuclear magnetic resonance (NMR) spectroscopy. This method, termed "restrained soft-docking," is validated for several known protein complexes. These data demonstrate that restrained soft-docking extends the size limitations of NMR spectroscopy and provides an alternative method for investigating macromolecular protein complexes that requires less experimental time, effort, and resources. The potential utility of this novel NMR and simulated docking approach in current structural genomic initiatives is discussed.

Ribotoxins are a family of potent cytotoxic proteins from Aspergillus whose members display a high sequence identity (85% for about 150 amino acid residues). The three-dimensional structures of two of these proteins, alpha-sarcin and restrictocin, are known. They interact with phospholipid bilayers, according to their ability to enter cells, and cleave a specific phosphodiester bond in the large subunit of ribosome thus inhibiting protein biosynthesis. Two nonconservative sequence changes between these proteins are located at the amino-terminal beta-hairpin of alpha-sarcin, a characteristic structure that is absent in other nontoxic structurally related microbial RNases. These two residues of alpha-sarcin, Lys 11 and Thr 20, have been substituted with the equivalent amino acids in restrictocin. The single mutants (K11L and T20D) and the corresponding K11L/T20D double mutant have been produced in Escherichia coli and purified to homogeneity. The spectroscopic characterization of the purified proteins reveals that the overall native structure is preserved. The ribonuclease and lipid-perturbing activities of the three mutants and restrictocin have been evaluated and compared with those of alpha-sarcin. These proteins exhibit the same ability to specifically inactivate ribosomes, although they show different activity against nonspecific substrate analogs such as poly(A). The mutant variant K11L and restrictocin display a lower phospholipid-interacting ability correlated with a decreased cytotoxicity. The results obtained are interpreted in terms of the involvement of the amino-terminal beta-hairpin in the interaction with both membranes and polyadenylic acid.

The net charge and isoelectric pH (pI) of a protein depend on the content of ionizable groups and their pK values. Ribonuclease Sa (RNase Sa) is an acidic protein with a pI = 3.5 that contains no Lys residues. By replacing Asp and Glu residues on the surface of RNase Sa with Lys residues, we have created a 3K variant (D1K, D17K, E41K) with a pI = 6.4 and a 5K variant (3K + D25K, E74K) with a pI = 10.2. We show that pI values estimated using pK values based on model compound data can be in error by >1 pH unit, and suggest how the estimation can be improved. For RNase Sa and the 3K and 5K variants, the solubility, activity, and stability have been measured as a function of pH. We find that the pH of minimum solubility varies with the pI of the protein, but that the pH of maximum activity and the pH of maximum stability do not.

The tRNA splicing endoribonuclease EndA from Methanococcus jannaschii is a homotetramer formed via heterologous interaction between the two pairs of homodimers. Each monomer consists of two alpha/beta domains, the N-terminal domain (NTD) and the C-terminal domain (CTD) containing the RNase A-like active site. Comparison of the EndA coordinates with the publicly available protein structure database revealed the similarity of both domains to site-specific deoxyribonucleases: the NTD to the LAGLIDADG family and the CTD to the PD-(D/E)XK family. Superposition of the NTD on the catalytic domain of LAGLIDADG homing endonucleases allowed a suggestion to be made about which amino acid residues of the tRNA splicing nuclease might participate in formation of a presumptive cryptic deoxyribonuclease active site. On the other hand, the CTD and PD-(D/E)XK endonucleases, represented by restriction enzymes and a phage lambda exonuclease, were shown to share extensive similarities of the structural framework, to which entirely different active sites might be attached in two alternative locations. These findings suggest that EndA evolved from a fusion protein with at least two distinct endonuclease activities: the ribonuclease, which made it an essential "antitoxin" for the cells whose RNA genes were interrupted by introns, and the deoxyribonuclease, which provided the means for homing-like mobility. The residues of the noncatalytic CTDs from the positions corresponding to the catalytic side chains in PD-(D/E)XK deoxyribonucleases map to the surface at the opposite side to the tRNA binding site, for which no function has been implicated. Many restriction enzymes from the PD-(D/E)XK superfamily might have the potential to maintain an additional active or binding site at the face opposite the deoxyribonuclease active site, a property that can be utilized in protein engineering.

Sequence similarity is the most common measure currently used to infer homology between proteins. Typically, homologous protein domains show sequence similarity over their entire lengths. Here we identify Asp box motifs, initially found as repeats in sialidases and neuraminidases, in new structural and sequence contexts. These motifs represent significantly similar sequences, localized to beta hairpins within proteins that are otherwise different in sequence and three-dimensional structure. By performing a combined sequence- and structure-based analysis we detect Asp boxes in more than nine protein families, including bacterial ribonucleases, sulfite oxidases, reelin, netrins, some lipoprotein receptors, and a variety of glycosyl hydrolases. Although the function common to each of these proteins, if any, remains unclear, we discuss possible functions of Asp boxes on the basis of previously determined experimental results and discuss different evolutionary scenarios for the origin of Asp-box containing proteins.

Theoretical and experimental studies have shown that the large desolvation penalty required for polar and charged groups frequently precludes their involvement in electrostatic interactions that contribute strongly to net stability in the folding or binding of proteins in aqueous solution near room temperature. We have previously developed a theoretical framework for computing optimized electrostatic interactions and illustrated use of the algorithm with simplified geometries. Given a receptor and model assumptions, the method computes the ligand-charge distribution that provides the most favorable balance of desolvation and interaction effects on binding. In this paper the method has been extended to treat complexes using actual molecular shapes. The barnase-barstar protein complex was investigated with barnase treated as a target receptor. The atomic point charges of barstar were varied to optimize the electrostatic binding free energy. Barnase and natural barstar form a tight complex (K(d) approximately 10(-14) M) with many charged and polar groups near the interface that make this a particularly relevant system for investigating the role of electrostatic effects on binding. The results show that sets of barstar charges (resulting from optimization with different constraints) can be found that give rise to relatively large predicted improvements in electrostatic binding free energy. Principles for enhancing the effect of electrostatic interactions in molecular binding in aqueous environments are discussed in light of the optima. Our findings suggest that, in general, the enhancements in electrostatic binding free energy resulting from modification of polar and charged groups can be substantial. Moreover, a recently proposed definition of electrostatic complementarity is shown to be a useful tool for examining binding interfaces. Finally, calculational results suggest that wild-type barstar is closer to being affinity optimized than is barnase for their mutual binding, consistent with the known roles of these proteins.

The effects of engineered disulfide bonds on protein stability are poorly understood because they can influence the structure, dynamics, and energetics of both the native and denatured states. To explore the effects of two engineered disulfide bonds on the stability of barnase, we have conducted a combined molecular dynamics and NMR study of the denatured state of the two mutants. As expected, the disulfide bonds constrain the denatured state. However, specific extended beta-sheet structure can also be detected in one of the mutant proteins. This mutant is also more stable than would be predicted. Our study suggests a possible cause of the very high stability conferred by this disulfide bond: the wild-type denatured ensemble is stabilized by a nonnative hydrophobic cluster, which is constrained from occurring in the mutant due to the formation of secondary structure.

Proteins often require cofactors to perform their biological functions and must fold in the presence of their cognate ligands. Using circular dichroism spectroscopy. we investigated the effects of divalent metal binding upon the folding pathway of Escherichia coli RNase HI. This enzyme binds divalent metal in its active site, which is proximal to the folding core of RNase HI as defined by hydrogen/deuterium exchange studies. Metal binding increases the apparent stability of native RNase HI chiefly by reducing the unfolding rate. As with the apo-form of the protein, refolding from high denaturant concentrations in the presence of Mg2+ follows three-state kinetics: formation of a rapid burst phase followed by measurable single exponential kinetics. Therefore, the overall folding pathway of RNase HI is minimally perturbed by the presence of metal ions. Our results indicate that the metal cofactor enters the active site pocket only after the enzyme reaches its native fold, and therefore, divalent metal binding stabilizes the protein by decreasing its unfolding rate. Furthermore, the binding of the cofactor is dependent upon a carboxylate critical for activity (Asp10). A mutation in this residue (D10A) alters the folding kinetics in the absence of metal ions such that they are similar to those observed for the unaltered enzyme in the presence of metal.

Several recent studies have shown that it is possible to increase protein stability by improving electrostatic interactions among charged groups on the surface of the folded protein. However, the stability increases are considerably smaller than predicted by a simple Coulomb's law calculation, and in some cases, a charge reversal on the surface leads to a decrease in stability when an increase was predicted. These results suggest that favorable charge-charge interactions are important in determining the denatured state ensemble, and that the free energy of the denatured state may be decreased more than that of the native state by reversing the charge of a side chain. We suggest that when the hydrophobic and hydrogen bonding interactions that stabilize the folded state are disrupted, the unfolded polypeptide chain rearranges to compact conformations with favorable long-range electrostatic interactions. These charge-charge interactions in the denatured state will reduce the net contribution of electrostatic interactions to protein stability and will help determine the denatured state ensemble. To support this idea, we show that the denatured state ensemble of ribonuclease Sa is considerably more compact at pH 7 where favorable charge-charge interactions are possible than at pH 3, where unfavorable electrostatic repulsion among the positive charges causes an expansion of the denatured state ensemble. Further support is provided by studies of the ionic strength dependence of the stability of charge-reversal mutants of ribonuclease Sa. These results may have important implications for the mechanism of protein folding.

We describe a simple experimental approach for the rapid determination of protein global folds. This strategy utilizes site-directed spin labeling (SDSL) in combination with isotope enrichment to determine long-range distance restraints between amide protons and the unpaired electron of a nitroxide spin label using the paramagnetic effect on relaxation rates. The precision and accuracy of calculating a protein global fold from only paramagnetic effects have been demonstrated on barnase, a well-characterized protein. Two monocysteine derivatives of barnase, (H102C) and (H102A/Q15C), were 15N enriched, and the paramagnetic nitroxide spin label, MTSSL, attached to the single Cys residue of each. Measurement of amide 1H longitudinal relaxation times, in both the oxidized and reduced states, allowed the determination of the paramagnetic contribution to the relaxation processes. Correlation times were obtained from the frequency dependence of these relaxation processes at 800, 600, and 500 MHz. Distances in the range of 8 to 35 A were calculated from the magnitude of the paramagnetic contribution to the relaxation processes and individual amide 1H correlation times. Distance restraints from the nitroxide spin to amide protons were used as restraints in structure calculations. Using nitroxide to amide 1H distances as long-range restraints and known secondary structure restraints, barnase global folds were calculated having backbone RMSDs <3 A from the crystal structure. This approach makes it possible to rapidly obtain the overall topology of a protein using a limited number of paramagnetic distance restraints.

Thermal stability of mutant proteins has been investigated using temperature dependent molecular dynamics (MD) simulations in vacuo. The numerical modeling was aimed at mimicking protein expansion upon heating. After the conditions for an expanding protein accessible surface area were established for T4 lysozyme and barnase wild-type proteins, MD simulations were carried out under the same conditions using the crystal structures of several mutant proteins. The computed thermal expansion of the accessible surface area of mutant proteins was found to be strongly correlated with their experimentally measured stabilities. A similar, albeit weaker, correlation was observed for model mutant proteins. This opens the possibility of obtaining stability information directly from protein structure.

Dimeric seminal RNase (BS-RNase) is an equilibrium mixture of conformationally different quaternary structures, one characterized by the interchange between subunits of their N-terminal ends (the MXM form); the other with no interchange (the M=M form). Controlled tryptic digestion of each isolated quaternary form generates, as limit digest products, folded and enzymatically active molecules, very resistant to further tryptic degradation. Electrospray mass spectrometric analyses and N-terminal sequence determinations indicate that trypsin can discriminate between the conformationally different quaternary structures of seminal RNase, and exerts a differential and asymmetric action on the two dimeric forms, depending on the original quaternary conformation of each form. The two digestion products from the MXM and the M=M dimeric forms have different structures, which are reminiscent of the original quaternary conformation of the dimers: one with interchange, the other with no interchange, of the N-terminal ends. The surprising resistance of these tryptic products to further tryptic action is explained by the persistence in each digestion product of the original intersubunit interface.

Analysis of the heterogeneity of packing in proteins showed that different groups of the protein preferentially contribute to low- or high-density regions. Statistical distribution reveals the two preferable values for packing density in the form of two peaks. One peak occurs in the range of densities 0.55-0.65, the other occurs in the range 0.75-0.8. The high-density peak is originated primarily by high packing inside the hydrogen bonded backbone and to some extent by side chains. Polar/charged and apolar side chains both contribute to the low-density peak. The average packing density values of individual atomic groups significantly vary for backbone atoms as well as for side chain atoms. The carbonyl oxygen atoms of protein backbone and the end groups of side chains show lower packing density than the rest of the protein. The side-chain atomic groups of a secondary structure element when packed against the neighboring secondary structure element form stronger contacts with the side chains of this element than with its backbone. Analysis of the low-density regions around each buried peptide group was done for the set of proteins with different types of packing, including alpha-alpha, alpha-beta, and beta-beta packing. It was shown that cavities are regularly situated in the groove of secondary structure element packed against neighboring elements for all types of packing. Low density in the regions surrounding the peptide groups and the end groups of side chains can be explained by their positioning next to a cavity formed upon the association of secondary structure elements. The model proposed can be applied to the analysis of protein internal motions, mechanisms of cellular signal transduction, diffusion through protein matrix, and other events.

Recent evidence suggests that the net effect of electrostatics is generally to destabilize protein binding due to large desolvation penalties. A novel method for computing ligand-charge distributions that optimize the tradeoff between ligand desolvation penalty and favorable interactions with a binding site has been applied to a model for barnase. The result is a ligand-charge distribution with a favorable electrostatic contribution to binding due, in part, to ligand point charges whose direct interaction with the binding site is unfavorable, but which make strong intra-molecular interactions that are uncloaked on binding and thus act to lessen the ligand desolvation penalty.

RNase S is a complex consisting of two proteolytic fragments of RNase A: the S peptide (residues 1-20) and S protein (residues 21-124). RNase S and RNase A have very similar X-ray structures and enzymatic activities. Previous experiments have shown increased rates of hydrogen exchange and greater sensitivity to tryptic cleavage for RNase S relative to RNase A. It has therefore been asserted that the RNase S complex is considerably more dynamically flexible than RNase A. In the present study we examine the differences in the dynamics of RNase S and RNase A computationally, by MD simulations, and experimentally, using trypsin cleavage as a probe of dynamics. The fluctuations around the average solution structure during the simulation were analyzed by measuring the RMS deviation in coordinates. No significant differences between RNase S and RNase A dynamics were observed in the simulations. We were able to account for the apparent discrepancy between simulation and experiment by a simple model. According to this model, the experimentally observed differences in dynamics can be quantitatively explained by the small amounts of free S peptide and S protein that are present in equilibrium with the RNase S complex. Thus, folded RNase A and the RNase S complex have identical dynamic behavior, despite the presence of a break in polypeptide chain between residues 20 and 21 in the latter molecule. This is in contrast to what has been widely believed for over 30 years about this important fragment complementation system.

The failure to appreciate that the hydration of polar groups is a major contribution to the entropy of protein unfolding has led to considerable underestimates for the loss of configurational freedom when a protein chain folds.

Denaturant m values, the dependence of the free energy of unfolding on denaturant concentration, have been collected for a large set of proteins. The m value correlates very strongly with the amount of protein surface exposed to solvent upon unfolding, with linear correlation coefficients of R = 0.84 for urea and R = 0.87 for guanidine hydrochloride. These correlations improve to R = 0.90 when the effect of disulfide bonds on the accessible area of the unfolded protein is included. A similar dependence on accessible surface area has been found previously for the heat capacity change (delta Cp), which is confirmed here for our set of proteins. Denaturant m values and heat capacity changes also correlate well with each other. For proteins that undergo a simple two-state unfolding mechanism, the amount of surface exposed to solvent upon unfolding is a main structural determinant for both m values and delta Cp.

Residues P19, L28, C31, and C32 have been implicated (Di Donato A, Cafaro V, D'Alessio G, 1994, J Biol Chem 269:17394-17396; Mazzarella L, Vitagliano L, Zagari A, 1995, Proc Natl Acad Sci USA: forthcoming) with key roles in determining the dimeric structure and the N-terminal domain swapping of seminal RNase. In an attempt to have a clearer understanding of the structural and functional significance of these residues in seminal RNase, a series of mutants of pancreatic RNase A was constructed in which one or more of the four residues were introduced into RNase A. The RNase mutants were examined for: (1) the ability to form dimers; (2) the capacity to exchange their N-terminal domains; (3) resistance to selective cleavage by subtilisin; and (4) antitumor activity. The experiments demonstrated that: (1) the presence of intersubunit disulfides is both necessary and sufficient for engendering a stably dimeric RNase; (2) all four residues play a role in determining the exchange of N-terminal domains; (3) the exchange is the molecular basis for the RNase antitumor action; and (4) this exchange is not a prerequisite in an evolutionary mechanism for the generation of dimeric RNases.

A method is described to objectively identify hydrophobic clusters in proteins of known structure. Clusters are found by examining a protein for compact groupings of side chains. Compact clusters contain seven or more residues, have an average of 65% hydrophobic residues, and usually occur in protein interiors. Although smaller clusters contain only side-chain moieties, larger clusters enclose significant portions of the peptide backbone in regular secondary structure. These clusters agree well with hydrophobic regions assigned by more intuitive methods and many larger clusters correlate with protein domains. These results are in striking contrast with the clustering algorithm of J. Heringa and P. Argos (1991, J Mol Biol 220:151-171). That method finds that clusters located on a protein's surface are not especially hydrophobic and average only 3-4 residues in size. Hydrophobic clusters can be correlated with experimental evidence on early folding intermediates. This correlation is optimized when clusters with less than nine hydrophobic residues are removed from the data set. This suggests that hydrophobic clusters are important in the folding process only if they have enough hydrophobic residues.

The thermodynamics of barnase denaturation has been studied calorimetrically over a broad range of temperature and pH. It is shown that in acidic solutions the heat denaturation of barnase is well approximated by a 2-state transition. The heat denaturation of barnase proceeds with a significant increase of heat capacity, which determines the temperature dependencies of the enthalpy and entropy of its denaturation. The partial specific heat capacity of denatured barnase is very close to that expected for the completely unfolded protein. The specific denaturation enthalpy value extrapolated to 130 degrees C is also close to the value expected for the full unfolding. Therefore, the calorimetrically determined thermodynamic characteristics of barnase denaturation can be considered as characteristics of its complete unfolding and can be correlated with structural features--the number of hydrogen bonds, extent of van der Waals contacts, and the surface areas of polar and nonpolar groups. Using this information and thermodynamic information on transfer of protein groups into water, the contribution of various factors to the stabilization of the native structure of barnase has been estimated. The main contributors to the stabilization of the native state of barnase appear to be intramolecular hydrogen bonds. The contributions of van der Waals interactions between nonpolar groups and those of hydration effects of these groups are not as large if considered separately, but the combination of these 2 factors, known as hydrophobic interactions, is of the same order of magnitude as the contribution of hydrogen bonding.

The structures of two catalytically modified semisynthetic RNases obtained by replacing phenylalanine 120 with leucine and tyrosine have been determined and refined at a resolution of 2.0 A (R = 0.161 and 0.184, respectively). These structures have been compared with the refined 1.8-A structure (R = 0.204) of the fully active phenylalanine-containing enzyme (Martin PD, Doscher MS, Edwards BFP, 1987, J Biol Chem 262:15930-15938) and with the catalytically defective D121A (2.0 A, R = 0.172) and D121N (2.0 A, R = 0.186) analogs (deMel VSJ, Martin PD, Doscher MS, Edwards BFP, 1992, J Biol Chem 267:247-256). The movement away from the active site of the loop containing residues 65-72 is seen in all three catalytically defective analogs--F120L, D121A, and D121N--but not in the fully active (or hyperactive) F120Y. The insertion of the phenolic hydroxyl of Tyr 120 into a hydrogen-bonding network involving the hydroxyl group of Ser 123 and a water molecule in F120Y is the likely basis for the hyperactivity toward uridine 2',3'-cyclic phosphate previously found for this analog (Hodges RS, Merrifield RB, 1974, Int J Pept Protein Res 6:397-405) as well as the threefold increase in KM for cytidine 2',3'-cyclic phosphate found for this analog by ourselves.

S-peptide (residues 1-20) and S-protein (residues 21-124) are the enzymatically inactive products of the limited digestion of ribonuclease A by subtilisin. S-peptide binds S-protein with high affinity to form ribonuclease S, which has full enzymatic activity. Recombinant DNA technology was used to produce a fusion protein having three parts: carrier, spacer, and target. The two carriers used were the first 15 residues of S-peptide (S15) and a mutant S15 in which Asp 14 had been changed to Asn (D14N S15). The spacer consisted of three proline residues and a four-residue sequence recognized by factor Xa protease. The target was beta-galactosidase. The interaction between the S-peptide portion of the fusion protein and immobilized S-protein allowed for affinity purification of the fusion protein under denaturing (S15 as carrier) or nondenaturing (D14N S15 as carrier) conditions. A sensitive method was developed to detect the fusion protein after sodium dodecyl sulfate-polyacrylamide gel electrophoresis by its ribonuclease activity following activation with S-protein. S-peptide has distinct advantages over existing carriers in fusion proteins in that it combines a small size (> or = 15 residues), a tunable affinity for ligand (Kd > or = 10(-9) M), and a high sensitivity of detection (> or = 10(-16) mol in a gel).