Refolding of Myoglobins of Nonhomologous Sequences by the Island Model

The folding process of sea hare myoglobin was simulated by the island model, which does not rely on sequence homologies or statistical inference from database of known structure. Sea hare myoglobin has low sequence homology (28%), but high structural similarity, with sperm whale myoglobin, which was already simulated by the island model. Their structural similarity is shown physiochemically from the distribution of hydrophobic-residue pairs, that is, the key pairs for packing of the secondary structures. Irrelevant to the sequence homology, the secondary structures can be packed into the tertiary structure through the hydrophobic interactions among the amino acid pairs responsible for the local structure formation. The results on the two species of myoglobins indicate that, in contrast to other prediction methods, the island model is applicable to any type of protein without extra information other than the distribution of hydrophobic-residue pairs and the positions of the secondary structures. Consequently the present results provide another verification of the validity of the island model for elucidating the mechanisms of protein folding and predicting protein structures.     

Assembly of protein structure from sparse experimental data: An efficient Monte Carlo model

A new, efficient method for the assembly of protein tertiary structure from known, loosely encoded secondary structure restraints and sparse information about exact side chain contacts is proposed and evaluated. The method is based on a new, very simple method for the reduced modeling of protein structure and dynamics, where the protein is described as a lattice chain connecting side chain centers of mass rather than Cαs. The model has implicit built-in multibody correlations that simulate short- and long-range packing preferences, hydrogen bonding cooperativity and a mean force potential describing hydrophobic interactions. Due to the simplicity of the protein representation and definition of the model force field, the Monte Carlo algorithm is at least an order of magnitude faster than previously published Monte Carlo algorithms for structure assembly. In contrast to existing algorithms, the new method requires a smaller number of tertiary restraints for successful fold assembly; on average, one for every seven residues as compared to one for every four residues. For example, for smaller proteins such as the B domain of protein G, the resulting structures have a coordinate root mean square deviation (cRMSD), which is about 3 Å from the experimental structure; for myoglobin, structures whose backbone cRMSD is 4.3 Å are produced, and for a 247-residue TIM barrel, the cRMSD of the resulting folds is about 6 Å. As would be expected, increasing the number of tertiary restraints improves the accuracy of the assembled structures. The reliability and robustness of the new method should enable its routine application in model building protocols based on various (very sparse) experimentally derived structural restraints. Proteins 32:475–494, 1998. © 1998 Wiley-Liss, Inc.     

S-nitrosylation-induced conformational change in blackfin tuna myoglobin

S-nitrosylation is a post-translational protein modification that can alter the function of a variety of proteins. Despite the growing wealth of information that this modification may have important functional consequences, little is known about the structure of the moiety or its effect on protein tertiary structure. Here we report high-resolution x-ray crystal structures of S-nitrosylated and unmodified blackfin tuna myoglobin, which demonstrate that in vitro S-nitrosylation of this protein at the surface-exposed Cys-10 directly causes a reversible conformational change by "wedging" apart a helix and loop. Furthermore, we have demonstrated in solution and in a single crystal that reduction of the S-nitrosylated myoglobin with dithionite results in NO cleavage from the sulfur of Cys-10 and rebinding to the reduced heme iron, showing the reversibility of both the modification and the conformational changes. Finally, we report the 0.95-A structure of ferrous nitrosyl myoglobin, which provides an accurate structural view of the NO coordination geometry in the context of a globin heme pocket.     

Mechanism of protein folding: I. General considerations and refolding of myoglobin

To explain the rapidity of the process of protein folding, we cite two aspects of hydrophobic interaction: its long-range nature and the specificity of pairing after the formation of secondary structures. These two factors, when incorporated with the growth-type mechanism, can determine the folding pathway of proteins. This mechanism is applied to myoglobin. Appropriate introduction of side chains of amino acid residues and the heme group attached to His 93 yield a refolded tertiary structure that is in good agreement with the native structure.     

Contribution of protein conformation to stereochemistry and reactivity of the active center of heme proteins and enzymes. The existence of horseradish peroxidase conformations and their possible role in the catalysis mechanism

In the spectral region 350-800 nm at 4.2 K we measured magnetic circular dichroism (MCD) spectra of the pentacoordinated complex of protcheme with 2-methylimidazole, deoxyleghemoglobin, neutral and alkaline forms of reduced horseradish peroxidase in the equilibrium states, as well as in non-equilibrium states produced by low-temperature photolysis of their carbon monoxide derivatives. Earlier the corresponding results have been obtained for myoglobin, hemoglobin and cytochromes P-450 and P-420. The energies of Fe-N (proximal His) and Fe-N(pyrroles) bonds and their changes upon ligand binding in heme proteins and enzymes were compared with those in the model heme complex thus providing conformational contribution into stereochemistry of the active site. The examples of weak and strong conformational "pressure" on stereochemistry were analysed and observed. If conformational energy contribution into stereochemistry prevails the electronic one the heme stereochemistry remains unchanged on ligand binding as it was observed for leghemoglobin and alkaline horseradish peroxidase. The change of bond energies in myoglobin and hemoglobin on ligand binding are comparable with those in protein free pentacoordinated protoheme, giving an example of weak conformational contribution to heme stereochemistry. The role of protein conformation energy in the modulation of ligand binding properties of heme in leghemoglobin relative to those in myoglobins is discussed. The most striking result were obtained in the study of reduced horseradish peroxidase in the pH region of 6.0-10.2. It was found that such different perturbations as ligand binding and heme-linked ionization of the distal amino acid residue induce identical changes in heme stereochemistry. Neither heme-linked ionization in the carbon monoxide complex nor the geometry of Fe-Co bond affect the heme local structure of photoproducts. These and other findings suggest a very low conformation mobility of horseradish peroxidase whose protein constraints appear to allow only two preferable geometries of specific amino acid residues that form the heme pocket. The role of the two tertiary structure constraints on the heme in the mechanism of horseradish peroxidase function is discussed. It is supposed that one conformation produces a heme environment suitable for two-electron oxidation of the native enzyme to compound I by hydrogen peroxide while another conformation changes the heme stereochemistry in the direction favourable for back reduction of compound I by the substrate to the resting enzyme through two one-electron steps. The switch from one tertiary structure to another is expected to be induced by substrate bind     

Effect of gamma-irradiation on the molecular properties of myoglobin

To elucidate the effect of gamma-irradiation on the molecular properties of myoglobin, the secondary and tertiary structures, as well as the molecular weight size of the protein, were examined after irradiation at various irradiation doses. Gamma-irradiation of myoglobin solutions caused the disruption of the ordered structure of the protein molecules, as well as degradation, crosslinking, and aggregation of the polypeptide chains. A SDSPAGE study indicated that irradiation caused initial fragmentation of the proteins and subsequent aggregation, due to cross-linking of the protein molecules. The effect of irradiation on the protein was more significant at lower protein concentrations. Ascorbic acid protected against the degradation and aggregation of proteins by scavenging oxygen radicals that are produced by irradiation. A circular dichroism study showed that an increase of the irradiation decreased the alpha-helical content of myoglobin with a concurrent increase of the aperiodic structure content. Fluorescence spectroscopy indicated that irradiation increased the emission intensity that was excited at 280 nm.     

Stereochemical theory of the 3-dimensional structure of globular proteins. I. Highly helical intermediate structures

The authors analyze the physical prerequisites on which the proposed stereochemical theory of the three-dimensional structure of globular proteins is based. The theory represents a stereochemical modelling of the mechanism of protein self-organization suggested earlier by one of the authors. According to this mechanism, a highly helical intermediate structure(s) is formed at first and then it passes into the native one. In the highly-helical intermediate structure the arrangement of the polypeptide chain in space is the same as in the native structure. These two structures differ mainly by the secondary structure of the chain. The transition into the native structure proceeds under the effect of long-range interactions which transform the excess alpha-helices into beta-structural and irregular conformations. The so-called s-helices are considered (the alpha-helix, whose hydrophobic groups form a separate cluster on its surface). s-Helices can be obtained on the greater part of the polypeptide chain of any globular protein. In the unfolded protein chain they are the most stable and rapidly formed structures. It has been shown that namely s-helices are the initial blocks for the formation of the highly-helical intermediate structure. Stereochemical principles of the s-helix packing that permit to predict the three-dimensional structure of highly helical proteins have been found. According to these principles the highly helical structure represents the packing of hydrophobic surfaces and s-helices. In their turn, hydrophobic surfaces are formed as a result of complementary interaction of borders of hydrophobic clusters of two s-helices according to the "knob-hole" principle.     

Folding myoglobin within a sol-gel glass: protein folding constrained to a small volume

The unfolding and refolding reaction of myoglobin was examined in solution and within a porous silica sol-gel glass. The sol-gel pores constrain the protein to a volume that is the same size and shape as the folded native state accompanied by a few layers of water solvation. Denaturants such as low pH buffers can be diffused through the gel pores to the protein to initiate unfolding and refolding. Acid-induced unfolding was hindered by the steric constraints imposed by the gel pores such that more denaturing conditions were required within the gel than in solution to create the unfolded state. No new folding intermediates were observed. Refolding of myoglobin was not complete in millimolar pH 7 buffer alone. Addition of 25% glycerol to the pH 7 buffer resulted in nearly complete refolding, and the use of 1 M phosphate buffer resulted in complete refolding. The role of this cosolvent and salt in disrupting the ordered water surrounding the protein within the gel is discussed in light of the Hofmeister series and entropic trapping via a diminished hydrophobic effect within the gel. These results are consistent with the premises of folding models in which secondary and tertiary structures are considered to form within a compact conformation of the protein backbone.     

RNA kink turns to the left and to the right

A helix-loop-helix within the group I intron has most of the canonical sequence elements of a kink turn (K-turn), yet it bends in the opposite direction. The reverse K-turn kinks toward the major rather than the minor grooves of the flanking helices. This suggests that there are two distinct subclasses of tertiary structures that a K-turn secondary structure can adopt. The final structure may be specified by external factors, such as protein binding or the tertiary structural context, rather than the intrinsic conformation of the RNA.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

Structural and functional relationships in DnaK and DnaK756 heat-shock proteins from Escherichia coli.

The secondary structures of DnaK and the mutant DnaK756 heat-shock proteins from Escherichia coli have been investigated by Fourier transform infrared spectroscopy. The analysis of infrared data showed that DnaK and DnaK756 proteins have different secondary structures that are not affected by the presence of ATP or beta, gamma-methyleneadenosine 5'-triphosphate. The infrared data indicate also that the tertiary structures of DnaK and DnaK756 proteins are different and that DnaK protein undergoes conformational changes in its tertiary structure not only during binding of ATP but also during ATP hydrolysis. Using fluorescence spectroscopy of a single tryptophan located in the N-terminal domain of DnaK protein and fluorescence of 1,1'-bis(4-anilino)naphthalene-5,5'-disulfonic acid, which interacts with hydrophobic domains of DnaK protein, we were able to distinguish between two conformational states of DnaK protein. After binding of triphosphonucleotides, the C-terminal domain of DnaK protein changes in tertiary structure in such a way that fewer hydrophobic segments are exposed on the surface of the protein. After ATP hydrolysis, the number of hydrophobic segments on the surface of the protein is further reduced, and moreover the tertiary structure of the N-terminal domain of the protein changes. These data are discussed in terms of structural and functional relationships of both DnaK and DnaK756 proteins.     

Hydrogen exchange and the proteolytic degradation of ribonuclease A. The position of the centers of the high-temperature splitting by proteases and the structure of denatured molecules

The information on the high-temperature proteolytic degradation of RNase A has been analyzed. It has been shown that a few peptide bonds primarily splitted by trypsin, chymotrypsin and thermolysin are localized only in the N-terminal part of structural domain II of the native molecule. The same peptide bonds are splitted by proteases with the highest rate upon the denaturation in the presense of trifluoroethanol or the renaturation from concentrated urea solutions and after the desorganization of the native structure by the reduction of all four S-S bonds of RNase A. According to the data on hydrogen exchange in the native RNase A molecule, the dynamic stability of the tertiary structure of domain II is lower than that of domain I because of the lesser amount of the internal bulky nonpolar residues Val, Ile, and Phe. For the same reason, this part of the molecule in different nonnative forms of RNase A is less compact and more flexible and is splitted with the highest rate in the segment 31-39 enriched by long cationic residues Lys and Arg. A common feature of the conformation of the flexible disordered backbone of all RNase A nonnative structures considered is the predominance of short PPII helices, which provides a high rate of the restoration of the native secondary and tertiary structures upon renaturation or self-organization and global fluctuations of the native structure revealed by the hydrogen exchange and proteolytic degradation.     

Looking at proteins: representations, folding, packing, and design. Biophysical Society National Lecture, 1992.

Looking at proteins is an active process of interpretation and selection, emphasizing some features and deleting others. Multiple representations are needed, for such purposes as showing motions or conveying both the chain connectivity and the three-dimensional shape simultaneously. In studying and comparing protein structures, ideas are suggested about the determinants of tertiary structure and of folding (e.g., that Greek key beta barrels may fold up two strands at a time). The design and synthesis of new proteins "from scratch" provides a route toward the experimental testing of such ideas. It has also been a fruitful new perspective from which to look at structures, requiring such things as statistics on very narrowly defined structural categories and explicit attention to "negative design" criteria that actively block unwanted alternatives (e.g., reverse topology of a helix bundle, or edge-to-edge aggregation of beta sheets). Recently, the field of protein design has produced a rather unexpected general result: apparently we do indeed know enough to successfully design proteins that fold into approximately correct structures, but not enough to design unique, native-like structures. The degree of order varies considerably, but even the best designed material shows multiple conformations by NMR, more similar to a "molten globule" folding intermediate than to a well ordered native tertiary structure. In response to this conclusion, we are now working on systems that test useful questions with approximate structures (such as determining which factors most influence the choice of helix-bundle topology) and also analyzing how natural proteins achieve unique core conformations (e.g., for side chains on the interior side of a beta sheet, illustrated in the kinemages).     

Amino acid code of protein secondary structure

The calculation of protein three-dimensional structure from the amino acid sequence is a fundamental problem to be solved. This paper presents principles of the code theory of protein secondary structure, and their consequence--the amino acid code of protein secondary structure. The doublet code model of protein secondary structure, developed earlier by the author (Shestopalov, 1990), is part of this theory. The theory basis are: 1) the name secondary structure is assigned to the conformation, stabilized only by the nearest (intraresidual) and middle-range (at a distance no more than that between residues i and i + 5) interactions; 2) the secondary structure consists of regular (alpha-helical and beta-structural) and irregular (coil) segments; 3) the alpha-helices, beta-strands and coil segments are encoded, respectively, by residue pairs (i, i + 4), (i, i + 2), (i, i = 1), according to the numbers of residues per period, 3.6, 2, 1; 4) all such pairs in the amino acid sequence are codons for elementary structural elements, or structurons; 5) the codons are divided into 21 types depending on their strength, i.e. their encoding capability; 6) overlappings of structurons of one and the same structure generate the longer segments of this structure; 7) overlapping of structurons of different structures is forbidden, and therefore selection of codons is required, the codon selection is hierarchic; 8) the code theory of protein secondary structure generates six variants of the amino acid code of protein secondary structure. There are two possible kinds of model construction based on the theory: the physical one using physical properties of amino acid residues, and the statistical one using results of statistical analysis of a great body of structural data. Some evident consequences of the theory are: a) the theory can be used for calculating the secondary structure from the amino acid sequence as a partial solution of the problem of calculation of protein three-dimensional structure from the amino acid sequence, and the calculated secondary structure and codon strength distribution can be used for simulating the next step of protein folding; b) one can propose that the same secondary structures can be folded into different tertiary structures and, vice versa, different secondary structures can be folded into the same tertiary structures, provided codon distributions are considered also; c) codons can be considered as first elements of protein three-dimensional structure language.     

Fluorescence study of conformational transitions in the structure of myoglobin

Fluorescence studies of myoglobin and Mb-like structures, apomyoglobin and the complex of apo-Mb with protoporphyrin IX, reveal both the similarity between them, which is due to a common type of polypeptide chain folding, and the distinctions imposed by the influence of the prosthetic group. Close resemblance of structures of holomyoglobin and its metal-free analog, PPIX--apo-Mb, points to a key role of specific interactions between the protein and the protoporphyrin macrocycle rather than the Fe-protein bond in the formation of Mb-like structures. In PPIX--apo-Mb, both the hydrophobic core and the important ionic bonds between different structural elements () stabilizing the Mb structure are almost completely retained. The bond between Fe and proximal His-F8 allows additional integration of the structures of the heme cavity and the myoglobin molecule as a whole, providing its functional activity and highly cooperative conformational transitions. In all the myoglobin-like structures studied, a certain relationship is found between conformational states of the , the heme cavity, and the N-terminal part of the molecule. This is probably due to variations in the mutual orientation of the ABCDE and FGH helical domains, depending on the interactions between the protein, the prosthetic group, and the ligand in the heme crevice. The correlation between conformations of the N-terminal and heme regions found at a level of the globin tertiary structure is very important for understanding the mechanisms of homo- and heterotropic regulation in tetrameric hemoglobins.     

Defining the interacting regions between apomyoglobin and lipid membrane by hydrogen/deuterium exchange coupled to mass spectrometry

Sperm whale myoglobin can be considered as the model protein of the globin family. The pH-dependence of the interactions of apomyoglobin with lipid bilayers shares some similarities with the behavior of pore-forming domains of bacterial toxins belonging also to the globin family. Two different states of apomyoglobin bound to a lipid bilayer have been characterized by using hydrogen/deuterium exchange experiments and mass spectrometry. When bound to the membrane at pH 5.5, apomyoglobin remains mostly native-like and interacts through alpha-helix A. At pH 4, the binding is related to the stabilization of a partially folded state. In that case, alpha-helices A and G are involved in the interaction. At this pH, alpha-helix G, which is the most hydrophobic region of apomyoglobin, is available for interaction with the lipid bilayer because of the loss of the tertiary structure. Our results show the feasibility of such experiments and their potential for the characterization of various membrane-bound states of amphitropic proteins such as pore-forming domains of bacterial toxins. This is not possible with other high-resolution methods, because these proteins are usually in partially folded states when interacting with membranes.     

Multi-wavelength anomalous diffraction using medium-angle X-ray solution scattering (MADMAX)

Proteins are dynamic molecules whose function in virtually all biological processes requires conformational motion. Direct experimental probes of protein structure in solution are needed to characterize these motions. Anomalous scattering from proteins in solution has the potential to act as a precise molecular ruler to determine the positions of specific chemical groups or atoms within proteins under conditions in which structural changes can take place free from the constraints of crystal contacts. In solution, anomalous diffraction has two components: a set of cross-terms that depend on the relative location of the anomalous centers and the rest of the protein, and a set of pure anomalous terms that depend on the distances between the anomalous centers. The cross-terms are demonstrated here to be observable and to provide direct information about the distance between the anomalous center and the center of mass of the protein. The second set of terms appears immeasurably small in the context of current experimental capabilities. Here, we outline the theory underlying anomalous scattering from proteins in solution, predict the anomalous differences expected on the basis of atomic coordinate sets, and demonstrate the measurement of anomalous differences at the iron edge for solutions of myoglobin and hemoglobin.     

Exon structure conservation despite low sequence similarity: a relic of dramatic events in evolution?

The evolutionary significance of introns remains a mystery. The current availability of several complete eukaryotic genomes permits new studies to probe the possible function of these peculiar genomic features. Here we investigate the degree to which gene structure (intron position, phase and length) is conserved between homologous protein domains. We find that for certain extracellular-signalling and nuclear domains, gene structures are similar even when protein sequence similarity is low or not significant and sequences can only be aligned with a knowledge of protein tertiary structure. In contrast, other domains, including most intracellular signalling modules, show little gene structure conservation. Intriguingly, many domains with conserved gene structures, such as cytokines, are involved in similar biological processes, such as the immune response. This suggests that gene structure conservation may be a record of key events in evolution, such as the origin of the vertebrate immune system or the duplication of nuclear receptors in nematodes. The results suggest ways to detect new and potentially very remote homologues, and to construct phylogenies for proteins with limited sequence similarity.     

Finding the needle in the haystack: towards solving the protein-folding problem computationally

Prediction of protein tertiary structures from amino acid sequence and understanding the mechanisms of how proteins fold, collectively known as “the protein folding problem,” has been a grand challenge in molecular biology for over half a century. Theories have been developed that provide us with an unprecedented understanding of protein folding mechanisms. However, computational simulation of protein folding is still difficult, and prediction of protein tertiary structure from amino acid sequence is an unsolved problem. Progress toward a satisfying solution has been slow due to challenges in sampling the vast conformational space and deriving sufficiently accurate energy functions. Nevertheless, several techniques and algorithms have been adopted to overcome these challenges, and the last two decades have seen exciting advances in enhanced sampling algorithms, computational power and tertiary structure prediction methodologies. This review aims at summarizing these computational techniques, specifically conformational sampling algorithms and energy approximations that have been frequently used to study protein-folding mechanisms or to de novo predict protein tertiary structures. We hope that this review can serve as an overview on how the protein-folding problem can be studied computationally and, in cases where experimental approaches are prohibitive, help the researcher choose the most relevant computational approach for the problem at hand. We conclude with a summary of current challenges faced and an outlook on potential future directions.     

Simulating the folding of small proteins by use of the local minimum energy and the free solvation energy yields native-like structures

Assuming that the protein primary sequence contains all information required to fold a protein into its native tertiary structure, we propose a new computational approach to protein folding by distributing the total energy of the macromolecular system along the torsional axes.We further derive a new semiempirical equation to calculate the total energy of a macromolecular system including its free energy of solvation. The energy of solvation makes an important contribution to the stability of biological structures. The segregation of hydrophilic and hydrophobic domains is essential for the formation of micelles, lipid bilayers, and biological membranes, and it is also important for protein folding. The free energy of solvation consists of two components: one derived from interactions between the atoms of the protein, and the second resulting from interactions between the protein and the solvent. The latter component is expressed as a function of the fractional area of protein atoms accessible to the solvent.The protein-folding procedure described in this article consists of two successive steps: a theoretical transition from an ideal α helix to an ideal β sheet is first imposed on the protein conformation, in order to calculate an initial secondary structure. The most stable secondary structure is built from a combination of the lowest energy structures calculated for each amino acid during this transition. An angular molecular dynamics step is then applied to this secondary structure. In this computational step, the total energy of the system consisting of the sum of the torsional energy, the van der Waals energy, the electrostatic energy, and the solvation energy is minimized. This process yields 3-D structures of minimal total energy that are considered to be the most probable native-like structures for the protein.This method therefore requires no prior hypothesis about either the secondary or the tertiary structure of the protein and restricts the input of data to its sequence. The validity of the results is tested by comparing the crystalline and computed structures of four proteins, i.e., the avian and bovine pancreatic polypeptide (36 residues each), uteroglobin (70 residues), and the calcium-binding protein (75 residues); the C_α-C_α maps show significant homologies and the position of secondary structure domains; that of the α helices is particularly close.