Correlation of folding kinetics with the number and isomerization states of prolines in three homologous proteins of the RNase family

Defining Electron Bifurcation in the Electron-Transferring Flavoprotein Family

Electron bifurcation is the coupling of exergonic and endergonic redox reactions to simultaneously generate (or utilize) low- and high-potential electrons. It is the third recognized form of energy conservation in biology and was recently described for select electron-transferring flavoproteins (Etfs). Etfs are flavin-containing heterodimers best known for donating electrons derived from fatty acid and amino acid oxidation to an electron transfer respiratory chain via Etf-quinone oxidoreductase. Canonical examples contain a flavin adenine dinucleotide (FAD) that is involved in electron transfer, as well as a non-redox-active AMP. However, Etfs demonstrated to bifurcate electrons contain a second FAD in place of the AMP. To expand our understanding of the functional variety and metabolic significance of Etfs and to identify amino acid sequence motifs that potentially enable electron bifurcation, we compiled 1,314 Etf protein sequences from genome sequence databases and subjected them to informatic and structural analyses. Etfs were identified in diverse archaea and bacteria, and they clustered into five distinct well-supported groups, based on their amino acid sequences. Gene neighborhood analyses indicated that these Etf group designations largely correspond to putative differences in functionality. Etfs with the demonstrated ability to bifurcate were found to form one group, suggesting that distinct conserved amino acid sequence motifs enable this capability. Indeed, structural modeling and sequence alignments revealed that identifying residues occur in the NADH- and FAD-binding regions of bifurcating Etfs. Collectively, a new classification scheme for Etf proteins that delineates putative bifurcating versus nonbifurcating members is presented and suggests that Etf-mediated bifurcation is associated with surprisingly diverse enzymes.IMPORTANCE Electron bifurcation has recently been recognized as an electron transfer mechanism used by microorganisms to maximize energy conservation. Bifurcating enzymes couple thermodynamically unfavorable reactions with thermodynamically favorable reactions in an overall spontaneous process. Here we show that the electron-transferring flavoprotein (Etf) enzyme family exhibits far greater diversity than previously recognized, and we provide a phylogenetic analysis that clearly delineates bifurcating versus nonbifurcating members of this family. Structural modeling of proteins within these groups reveals key differences between the bifurcating and nonbifurcating Etfs.

Rigid Single Carbon-Carbon Bond That Does Not Rotate in Water

Carbon-carbon bond is one of the most ubiquitous molecular building blocks for natural and man-made materials. Rotational isomerization is fundamentally important for understanding the structure and reactivity of chemical and biological molecules. Reported herein is the first demonstration that a single C-C bond does not rotate in water. The two distal C-S bonds in both 1,2-ethanedithiolate ((-)S-CH2-CH2-S(-), 1,2-EDT(2-)) and 2,3-butanedithiolate (2,3-BuDT(2-)) are exclusively in the trans conformer with reference to their respective center single C-C bond. In contrast, both trans and gauche conformers are observed in neutral 1,2-ethanedithiol (1,2-EDT) and 2, 3-butanedithiol (2,3-BuDT). The insight from this work should be important for understanding the charge effect on the molecular conformation in aqueous solutions.

Dynamic Formation and Breaking of Disulfide Bonds in Molecular Dynamics Simulations with the UNRES Force Field

Many proteins contain disulfide bonds that are usually essential for maintaining function and a stable structure. Several algorithms attempt to predict the arrangement of disulfide bonds in the context of protein structure prediction, but none can simulate the entire process of oxidative folding, including dynamic formation and breaking of disulfide bonds. In this work, a potential function developed to model disulfide bonds is coupled with the united-residue (UNRES) force field, and used in both canonical and replica exchange molecular dynamics simulations to produce complete oxidative folding pathways. The potential function is obtained by introducing a transition barrier that separates the bonded and nonbonded states of the half-cystine residues. Tests on several helical proteins show that improved predictions are obtained when dynamic disulfide-bond formation and breaking are considered. The effect of the disulfide bonds on the folding kinetics is also investigated, particularly their role in stabilizing folding intermediates, resulting in slower folding.

Stability and folding of amphibian ribonuclease A superfamily members in comparison with mammalian homologues

Comparative studies on homologous proteins can provide knowledge on how limited changes in the primary structure find their expression in large effects on catalytic activity, stability or the folding behavior. For more than half a century, members of the ribonuclease A superfamily have been the subject of a myriad of studies on protein folding and stability. Both the unfolding and refolding kinetics as well as the structure of several folding intermediates of ribonuclease A have been characterized in detail. Moreover, the RNA-degrading activity of these enzymes provides a basis for their cytotoxicity, which renders them potential tumor therapeutics. Because amphibian ribonuclease A homologues evade the human ribonuclease inhibitor, they emerged as particularly promising candidates. Interestingly, the amphibian ribonuclease A homologues investigated to date are more stable than the mammalian homologues. Nevertheless, despite the generation of numerous genetically engineered variants, knowledge of the folding of amphibian ribonuclease A homologues remains rather limited. An exception is onconase, a ribonuclease A homologue from Rana pipiens, which has been characterized in detail. This review summarizes the data on the unfolding and refolding kinetics and pathways, as well on the stability of amphibian ribonuclease A homologues compared with those of ribonuclease A, the best known member of this superfamily.

A probabilistic approach to the effect of water hydrogen bonds on the kinetics of protein folding and protein denaturation

Previously, we presented a review of our kinetic models for the nucleation mechanism of protein folding and for the protein thermal denaturation in a barrierless way. A protein was treated as a random heteropolymer consisting of hydrophobic, hydrophilic, and neutral beads. As a crucial idea of the model, an overall potential around the cluster of native residues wherein a residue performs a chaotic motion was considered as the combination of the average dihedral, effective pairwise, and confining potentials. The overall potential as a function of the distance from the cluster has a double well shape. This allowed one to develop kinetic models for the nucleation mechanism of protein folding (NMPF) and barrierless protein denaturation (BPD) by using the mean first passage time analysis. In the original models, however, hydrogen bonding effects were taken into account only indirectly which affected the accuracy of the models because hydrogen bonding does play a crucial role in the folding, stability, and denaturation of proteins. To improve the NMPF and BPD models and explicitly take into account the hydrogen bonding "water-water" and "water-protein residue", we have developed a probabilistic hydrogen bond (PHB) model for the effect of hydrogen bond networks of water molecules around two solute particles (immersed in water) on their interaction, and have then combined the PHB model with the NMPF and BPD models. In this paper, that can be regarded as sequel of our previous review, we analyze the modified NMPF and BPD models that explicitly take into account the effect of water-water hydrogen bonding on these processes. As expected, the application of the modified models to the folding/unfolding of two model proteins (one short, consisting of 124 residues and the other large, consisting of 2500 residues) demonstrate that the hydrogen bond networks play a very important role in the protein folding/unfolding phenomena.

Ribonucleases as potential modalities in anticancer therapy

Antitumor ribonucleases are small (10-28 kDa) basic proteins. They were found among members of both, ribonuclease A and T1 superfamilies. Their cytotoxic properties are conferred by enzymatic activity, i.e., the ability to catalyze cleavages of phosphodiester bonds in RNA. They bind to negatively charged cell membrane, enter cells by endocytosis and translocate to cytosol where they evade mammalian protein ribonuclease inhibitor and degrade RNA. Here, we discuss structures, functions and mechanisms of antitumor activity of several cytotoxic ribonucleases with particular emphasis to the amphibian Onconase, the only enzyme of this class that reached clinical trials. Onconase is the smallest, very stable, less catalytically efficient and more cytotoxic than most RNase A homologues. Its cytostatic, cytotoxic and anticancer effects were extensively studied. It targets tRNA, rRNA, mRNA as well as the non-coding RNA (microRNAs). Numerous cancer lines are sensitive to Onconase; their treatment with 10-100 nM enzyme leads to suppression of cell cycle progression, predominantly through G(1), followed by apoptosis or cell senescence. Onconase also has anticancer properties in animal models. Many effects of this enzyme are consistent with the microRNAs, one of its critical targets. Onconase sensitizes cells to a variety of anticancer modalities and this property is of particular interest, suggesting its application as an adjunct to chemotherapy or radiotherapy in treatment of different tumors. Cytotoxic RNases as exemplified by Onconase represent a new class of antitumor agents, with an entirely different mechanism of action than the drugs currently used in the clinic. Further studies on animal models including human tumors grafted on severe combined immunodefficient (SCID) mice and clinical trials are needed to explore clinical potential of cytotoxic RNases.

Effect of hydrogen bond networks on the nucleation mechanism of protein folding

We have recently developed a kinetic model for the nucleation mechanism of protein folding (NMPF) in terms of ternary nucleation by using the first passage time analysis. A protein was considered as a random heteropolymer consisting of hydrophobic, hydrophilic (some of which are negatively or positively ionizable), and neutral beads. The main idea of the NMPF model consisted of averaging the dihedral potential in which a selected residue is involved over all possible configurations of all neighboring residues along the protein chain. The combination of the average dihedral, effective pairwise (due to Lennard-Jones-type and electrostatic interactions), and confining (due to the polymer connectivity constraint) potentials gives rise to an overall potential around the cluster that, as a function of the distance from the cluster center, has a double-well shape. This allows one to evaluate the protein folding time. In the original NMPF model hydrogen bonding was not taken into account explicitly. To improve the NMPF model and make it more realistic, in this paper we modify our (previously developed) probabilistic hydrogen bond model and combine it with the former. Thus, a contribution due to the disruption of hydrogen bond networks around the interacting particles (cluster of native residues and residue in the protein unfolded part) appears in the overall potential field around a cluster. The modified model is applied to the folding of the same model proteins that were examined in the original model: a short protein consisting of 124 residues (roughly mimicking bovine pancreatic ribonuclease) and a long one consisting of 2500 residues (as a representative of large proteins with superlong polypeptide chains), at pH=8.3 , 7.3, and 6.3. The hydrogen bond contribution now plays a dominant role in the total potential field around the cluster (except for very short distances thereto where the repulsive energy tends to infinity). It is by an order of magnitude stronger for hydrophobic residues than for hydrophilic ones. The range of "residue-cluster" distances, at which the hydrogen bond effect exists, is twice as long for hydrophobic residues as for hydrophilic ones.

The folding pathway of onconase is directed by a conserved intermediate

A promising approach to unravel the relationship between sequence information, tertiary structure, and folding mechanism of proteins is the analysis of the folding behavior of proteins with low sequence identity but comparable tertiary structures. Ribonuclease A (RNase A) and its homologues, forming the RNase A superfamily, provide an excellent model system for respective studies. RNase A has been used extensively as a model protein for folding studies. However, little is known about the folding of homologous RNases. Here, we analyze the folding pathway of onconase, a homologous protein from the Northern leopard frog with great potential as a tumor therapeutic, by high-resolution techniques. Although onconase and RNase A significantly differ in the primary structure (28% sequence identity) and in thermodynamic stability (DeltaDeltaG = 20 kJ mol(-1)), both enzymes possess very similar tertiary structures. The present folding studies on onconase by rapid mixing techniques in combination with fluorescence and NMR spectroscopy allow the structural assignment of the three kinetic phases observed in stopped-flow fluorescence spectroscopy. After a slow peptidyl-prolyl cis-to-trans isomerization reaction in the unfolded state, ONC folds via an on-pathway intermediate to the native state. By quenched-flow hydrogen/deuterium exchange experiments coupled with 2D NMR spectroscopy, 31 amino acid residues were identified to be involved in the structure formation of the intermediate. Twelve of these residues are identical in the RNase A sequence, which is a significantly higher percentage (39%) than the overall 28% sequence identity. Moreover, the structure of this intermediate closely resembles two of the intermediates that occur early during the refolding of RNase A. Obviously, in spite of considerable differences in their amino acid sequence the initial folding events of both proteins are comparable, guided by a limited number of conserved residues.

Untangling the folding mechanism of the 5(2)-knotted protein UCH-L3

Proteins possessing deeply embedded topological knots in their structure add a stimulating new challenge to the already complex protein-folding problem. The most complicated knotted topology observed to date belongs to the human enzyme ubiquitin C-terminal hydrolase UCH-L3, which is an integral part of the ubiquitin-proteasome system. The structure of UCH-L3 contains five distinct crossings of its polypeptide chain, and it adopts a 5(2)-knotted topology, making it a fascinating target for folding studies. Here, we provide the first in depth characterization of the stability and folding of UCH-L3. We show that the protein can unfold and refold reversibly in vitro without the assistance of molecular chaperones, demonstrating that all the information necessary for the protein to find its knotted native structure is encoded in the amino acid sequence, just as with any other globular protein, and that the protein does not enter into any deep kinetic traps. Under equilibrium conditions, the unfolding of UCH-L3 appears to be two-state, however, multiphasic folding and unfolding kinetics are observed and the data are consistent with a folding pathway in which two hyperfluorescent intermediates are formed. In addition, a very slow phase in the folding kinetics is shown to be limited by proline-isomerization events. Overall, the data suggest that a knotted topology, even in its most complex form, does not necessarily limit folding in vitro, however, it does seem to require a complex folding mechanism which includes the formation of several distinct intermediate species.

Back to the future: Ribonuclease A

Pancreatic ribonuclease A (EC 3.1.27.5, RNase) is, perhaps, the best-studied enzyme of the 20th century. It was isolated by René Dubos, crystallized by Moses Kunitz, sequenced by Stanford Moore and William Stein, and synthesized in the laboratory of Bruce Merrifield, all at the Rockefeller Institute/University. It has proven to be an excellent model system for many different types of experiments, both as an enzyme and as a well-characterized protein for biophysical studies. Of major significance was the demonstration by Chris Anfinsen at NIH that the primary sequence of RNase encoded the three-dimensional structure of the enzyme. Many other prominent protein chemists/enzymologists have utilized RNase as a dominant theme in their research. In this review, the history of RNase and its offspring, RNase S (S-protein/S-peptide), will be considered, especially the work in the Merrifield group, as a preface to preliminary data and proposed experiments addressing topics of current interest. These include entropy-enthalpy compensation, entropy of ligand binding, the impact of protein modification on thermal stability, and the role of protein dynamics in enzyme action. In continuing to use RNase as a prototypical enzyme, we stand on the shoulders of the giants of protein chemistry to survey the future.