Extent of hydrogen-bond protection in folded proteins: a constraint on packing architectures

Progressive structuring and ultimately exclusion of water by hydrophobes surrounding backbone hydrogen bonds turn the latter into guiding factors of protein folding. Here we demonstrate that an arrangement of five hydrophobes yields an optimal hydrogen-bond stabilization. This motif is shown to be nearly ubiquitous in native folds.     

The nature of folded states of globular proteins.

We suggest, using dynamical simulations of a simple heteropolymer modelling the alpha-carbon sequence in a protein, that generically the folded states of globular proteins correspond to statistically well-defined metastable states. This hypothesis, called the metastability hypothesis, states that there are several free energy minima separated by barriers of various heights such that the folded conformations of a polypeptide chain in each of the minima have similar structural characteristics but have different energies from one another. The calculated structural characteristics, such as bond angle and dihedral angle distribution functions, are assumed to arise from only those configurations belonging to a given minimum. The validity of this hypothesis is illustrated by simulations of a continuum model of a heteropolymer whose low temperature state is a well-defined beta-barrel structure. The simulations were done using a molecular dynamics algorithm (referred to as the "noisy" molecular dynamics method) containing both friction and noise terms. It is shown that for this model there are several distinct metastable minima in which the structural features are similar. Several new methods of analyzing fluctuations in structures belonging to two distinct minima are introduced. The most notable one is a dynamic measure of compactness that can in principle provide the time required for maximal compactness to be achieved. The analysis shows that for a given metastable state in which the protein has a well-defined folded structure the transition to a state of higher compactness occurs very slowly, lending credence to the notion that the system encounters a late barrier in the process of folding to the most compact structure. The examination of the fluctuations in the structures near the unfolding----folding transition temperature indicates that the transition state for the unfolding to folding process occurs closer to the folded state.     

Helix-helix packing and interfacial pairwise interactions of residues in membrane proteins

Helix-helix packing plays a critical role in maintaining the tertiary structures of helical membrane proteins. By examining the overall distribution of voids and pockets in the transmembrane (TM) regions of helical membrane proteins, we found that bacteriorhodopsin and halorhodopsin are the most tightly packed, whereas mechanosensitive channel is the least tightly packed. Large residues F, W, and H have the highest propensity to be in a TM void or a pocket, whereas small residues such as S, G, A, and T are least likely to be found in a void or a pocket. The coordination number for non-bonded interactions for each of the residue types is found to correlate with the size of the residue. To assess specific interhelical interactions between residues, we have developed a new computational method to characterize nearest neighboring atoms that are in physical contact. Using an atom-based probabilistic model, we estimate the membrane helical interfacial pairwise (MHIP) propensity. We found that there are many residue pairs that have high propensity for interhelical interactions, but disulfide bonds are rarely found in the TM regions. The high propensity pairs include residue pairs between an aromatic residue and a basic residue (W-R, W-H, and Y-K). In addition, many residue pairs have high propensity to form interhelical polar-polar atomic contacts, for example, residue pairs between two ionizable residues, between one ionizable residue and one N or Q. Soluble proteins do not share this pattern of diverse polar-polar interhelical interaction. Exploratory analysis by clustering of the MHIP values suggests that residues similar in side-chain branchness, cyclic structures, and size tend to have correlated behavior in participating interhelical interactions. A chi-square test rejects the null hypothesis that membrane protein and soluble protein have the same distribution of interhelical pairwise propensity. This observation may help us to understand the folding mechanism of membrane proteins.     

The solvent dependence of hydrogen exchange kinetics of folded proteins.

The effects of ethanol, ethylene glycol, dioxane, and other organic co-solvents upon the hydrogen exchange rates of randomly coiled oxidized RNase, native RNase, and native trypsin have been measured. The exchange rate of oxidized RNase, the model compound for the proton transfer step in hydrogen exchange, is decreased by all of the co-solvents studied at temperatures in the range 3-20 degrees. This has been ascribed to the combined effects of the disruption of peptide bond solvation due to a reduction in the concentration of water, and of changes in [OH-] ion concentration due to changes in the acid dissociation constant of water, Kw. The solvent dependence for both native RNase and native trypsin is similar in all of the solvents studied. At a low temperature (3-20 degrees), the exchange rates go through a minimum as the solvent concentration is increased. At higher temperatures (20-35 degrees) the exchange rates are increased at all concentrations of the co-solvent. The apparent rate minimum at lower temperatures is due to two opposing effects. Co-solvents decrease the rate of exchange that occurs directly from the folded molecule. At higher concentrations and higer temperature. The decrease in rates for exchange directly from folded protein is primarily due to the effects on the proton transfer step, and not to binding or the solvent effects on protein structure. The solvents used in this study have no apparent effect on conformational processes contributing to the hydrogen exchange process in folded proteins.     

The role of ribosomal proteins in in vitro ribosome-membrane interactions]

The in vitro binding of total ribosomal proteins with rough endoplasmic membranes, from which 70% of ribosomes are eliminated by EDTA (ME) is studied. It is found that in conditions of specific interaction of ribosomes with membranes about 75% of total ribosomal proteins are bound with ME. Membranes, heterogenous in their content (different protein/lipid ratio), became homogenous in their buyoant density after the binding with proteins. The ability of membrane-ribosomal protein complex to bind ribosomes is not decreased, as it can be expected, but is considerablly increased, thus indicating on a non-specific character of ribosome binding. Ribosomal subunits lacking about half of structural protein are capable to bind with ribosome-binding membrane receptors and with some additional sites. This binding is also non-specific, because the binding efficiency of large and small subunits is the same.     

Probing origins of molecular interactions stabilizing the membrane proteins halorhodopsin and bacteriorhodopsin

Single-molecule atomic force microscopy and spectroscopy were applied to detect molecular interactions stabilizing the structure of halorhodopsin (HR), a light-driven chloride pump from Halobacterium salinarum. Because of the high structural and sequence similarities between HR and bacteriorhodopsin, we compared their unfolding pathways and polypeptide regions that established structurally stable segments against unfolding. Unfolding pathways and structural segments stabilizing the proteins both exhibited a remarkably high similarity. This suggests that different amino acid compositions can establish structurally indistinguishable energetic barriers. These stabilizing domains rather result from comprehensive interactions of all amino acids within a structural region than from specific interactions. However, one additional unfolding barrier located within a short segment of helix E was detected for HR. This barrier correlated with a Pi-bulk interaction, which locally disrupts helix E and divides a structural stabilizing segment.     

Helix packing in polytopic membrane proteins: role of glycine in transmembrane helix association.

The nature and distribution of amino acids in the helix interfaces of four polytopic membrane proteins (cytochrome c oxidase, bacteriorhodopsin, the photosynthetic reaction center of Rhodobacter sphaeroides, and the potassium channel of Streptomyces lividans) are studied to address the role of glycine in transmembrane helix packing. In contrast to soluble proteins where glycine is a noted helix breaker, the backbone dihedral angles of glycine in transmembrane helices largely fall in the standard alpha-helical region of a Ramachandran plot. An analysis of helix packing reveals that glycine residues in the transmembrane region of these proteins are predominantly oriented toward helix-helix interfaces and have a high occurrence at helix crossing points. Moreover, packing voids are generally not formed at the position of glycine in folded protein structures. This suggests that transmembrane glycine residues mediate helix-helix interactions in polytopic membrane proteins in a fashion similar to that seen in oligomers of membrane proteins with single membrane-spanning helices. The picture that emerges is one where glycine residues serve as molecular notches for orienting multiple helices in a folded protein complex.     

The role of non-native interactions in the folding of knotted proteins

Stochastic simulations of coarse-grained protein models are used to investigate the propensity to form knots in early stages of protein folding. The study is carried out comparatively for two homologous carbamoyltransferases, a natively-knotted N-acetylornithine carbamoyltransferase (AOTCase) and an unknotted ornithine carbamoyltransferase (OTCase). In addition, two different sets of pairwise amino acid interactions are considered: one promoting exclusively native interactions, and the other additionally including non-native quasi-chemical and electrostatic interactions. With the former model neither protein shows a propensity to form knots. With the additional non-native interactions, knotting propensity remains negligible for the natively-unknotted OTCase while for AOTCase it is much enhanced. Analysis of the trajectories suggests that the different entanglement of the two transcarbamylases follows from the tendency of the C-terminal to point away from (for OTCase) or approach and eventually thread (for AOTCase) other regions of partly-folded protein. The analysis of the OTCase/AOTCase pair clarifies that natively-knotted proteins can spontaneously knot during early folding stages and that non-native sequence-dependent interactions are important for promoting and disfavouring early knotting events.     

The role of internal packing interactions in determining the structure and stability of a protein

Cassette mutagenesis has been used to investigate how internal packing interactions help to specify a protein's three-dimensional structure and stability. Three interacting residues in the hydrophobic core of the N-terminal domain of lambda repressor were randomized combinatorially. The randomization was restricted to the five amino acids Val, Leu, Ile, Met and Phe, thereby generating a sterically diverse set of core sequences composed solely of hydrophobic residues. We have isolated 78 of the 125 possible sequences generated by this randomization. Approximately 70% of the isolated sequences show some level of biological activity, and thus still carry sufficient information to encode the basic structure of lambda repressor. An assay based on the temperature dependence of activity in vivo has been used to estimate the relative activities and thermal stabilities of the set of mutants. In addition, nine mutants have been purified and their stabilities and DNA binding activities characterized in vitro. Of the 56 active sequences, only two, in addition to the wild-type, maintain the wild-type level of stability and activity. All three of these proteins satisfy stringent requirements for specifically shaped residues at each position. All of the remaining active sequences have reduced stabilities and/or reduced DNA binding affinities. These and previous results suggest that there are two levels of structural information encoded in core residues. At the first level, the basic structural information appears to reside largely in the hydrophobic character of these residues. The majority of sequences that simply maintain hydrophobicity at core positions are able to adopt the overall lambda repressor fold and maintain moderate stability. At the second, more detailed level, specific steric features of these residues and their packing interactions clearly act as important determinants of the protein's precise structure and stability. These results imply that many of the basic structural features of a protein could be predicted from relatively simple, degenerate sequence patterns.     

The role of ion-lipid interactions and lipid packing in transient defects caused by phenolic compounds

The transient disruption of membranes for the passive permeation of ions or small molecules is a complex process relevant to understanding physiological processes and biotechnology applications. Phenolic compounds are widely studied for their antioxidant and antimicrobial properties, and some of these activities are based on the interactions of the phenolic compound with membranes. Ions are ubiquitous in cells and are known to alter the structure of phospholipid bilayers. Yet, ion-lipid interactions are usually ignored when studying the membrane-altering properties of phenolic compounds. This study aims to assess the role of Ca²⁺ ions on the membrane-disrupting activity of two phenolic acids and to highlight the role of local changes in lipid packing in forming transient defects or pores. Results from tethered bilayer lipid membrane electrical impedance spectroscopy experiments showed that Ca²⁺ significantly reduces membrane disruption by caffeic acid methyl ester and caffeic acid. As phenolic acids are known metal chelators, we used UV-vis and fluorescence spectroscopy to exclude the possibility that Ca²⁺ interferes with membrane disruption by binding to the phenolic compound and subsequently preventing membrane binding. Molecular dynamics simulations showed that Ca²⁺ but not caffeic acid methyl ester or caffeic acid increases lipid packing in POPC bilayers. The combined data confirm that Ca²⁺ reduces the membrane-disrupting activity of the phenolic compounds, and that Ca²⁺-induced changes to lipid packing govern this effect. We discuss our data in the context of ion-induced pores and transient defects and how lipid packing affects membrane disruption by small molecules.     

The role of neisserial Opa proteins in interactions with host cells

Pathogenic Neisseria spp. possess a repertoire of phase-variable Opa proteins that mediate various pathogen–host cell interactions, including bacterial engulfment by epithelial cells and opsonin-independent phagocytosis by professional phagocytes. Recent studies have identified cellular targets recognized by defined Opa proteins and have begun to reveal host signalling events involved in mediating these Opa-dependent cellular processes.     

Arginine residues as stabilizing elements in proteins.

Site-specific substitutions of arginine for lysine in the thermostable D-xylose isomerase (XI) from Actinoplanes missouriensis are shown to impart significant heat stability enhancement in the presence of sugar substrates most probably by interfering with nonenzymatic glycation. The same substitutions are also found to increase heat stability in the absence of any sugar derivatives, where a mechanism based on prevention of glycation can no longer be invoked. This rather conservative substitution is moreover shown to improve thermostability in two other structurally unrelated proteins, human copper, zinc-superoxide dismutase (CuZnSOD) and D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from Bacillus subtilis. The stabilizing effect of Lys----Arg substitutions is rationalized on the basis of a detailed analysis of the crystal structures of wild-type XI and of engineered variants with Lys----Arg substitution at four distinct locations, residues 253, 309, 319, and 323. Molecular model building analysis of the structures of wild-type and mutant CuZnSOD (K9R) and GAPDH (G281K and G281R) is used to explain the observed stability enhancement in these proteins. In addition to demonstrating that even thermostable proteins can lend themselves to further stability improvement, our findings provide direct evidence that arginine residues are important stabilizing elements in proteins. Moreover, the stabilizing role of electrostatic interactions, particularly between subunits in oligomeric proteins, is documented.     

Side-chain entropy and packing in proteins.

What role does side-chain packing play in protein stability and structure? To address this question, we compare a lattice model with side chains (SCM) to a linear lattice model without side chains (LCM). Self-avoiding configurations are enumerated in 2 and 3 dimensions exhaustively for short chains and by Monte Carlo sampling for chains up to 50 main-chain monomers long. This comparison shows that (1) side-chain degrees of freedom increase the entropy of open conformations, but side-chain steric exclusion decreases the entropy of compact conformations, thus producing a substantial entropy that opposes folding; (2) there is a side-chain “freezing” or ordering, i.e., a sharp decrease in entropy, near maximum compactness; and (3) the different types of contacts among side chains (s) and main-chain elements (m) have different frequencies, and the frequencies have different dependencies on compactness. mm contacts contribute significantly only at high densities, suggesting that main-chain hydrogen bonding in proteins may be promoted by compactness. The distributions of mm, ms, and ss contacts in compact SCM configurations are similar to the distributions in protein structures in the Brookhaven Protein Data Bank. We propose that packing in proteins is more like the packing of nuts and bolts in a jar than like the pairwise matching of jigsaw puzzle pieces.     

Import of stably folded proteins into peroxisomes.

By virtue of their synthesis in the cytoplasm, proteins destined for import into peroxisomes are obliged to traverse the single membrane of this organelle. Because the targeting signal for most peroxisomal matrix proteins is a carboxy-terminal tripeptide sequence (SKL or its variants), these proteins must remain import competent until their translation is complete. We sought to determine whether stably folded proteins were substrates for peroxisomal import. Prefolded proteins stabilized with disulfide bonds and chemical cross-linkers were shown to be substrates for peroxisomal import, as were mature folded and disulfide-bonded IgG molecules containing the peroxisomal targeting signal. In addition, colloidal gold particles conjugated to proteins bearing the peroxisomal targeting signal were translocated into the peroxisomal matrix. These results support the concept that proteins may fold in the mammalian cytosol, before their import into the peroxisome, and that protein unfolding is not a prerequisite for peroxisomal import.     

The role of hydrophobic interactions in positioning of peripheral proteins in membranes

Background  

Three-dimensional (3D) structures of numerous peripheral membrane proteins have been determined. Biological activity, stability, and conformations of these proteins depend on their spatial positions with respect to the lipid bilayer. However, these positions are usually undetermined.     

The number of charge-charge interactions stabilizing the ends of nucleosome DNA.

It has been shown by others that the melting of DNA in the nucleosome core particle is biphasic (ref. 1) and that the initial denaturation phase is due to melting of the DNA termini (refs. 1 & 2). We analyze the salt dependence of the melting temperature of this first transition and estimate that only 15% of the phosphates of the DNA termini are involved in intimate charge-charge interactions with histones. (The simplest model yields approximately 9%, whereas a calculated overestimate yields approximately 21% neutralization.) This is a surprisingly small number of interactions but we suggest that it may nonetheless be representative of all the core particle DNA.     

Membrane translocation of folded proteins

An ever-increasing number of proteins have been shown to translocate across various membranes of bacterial as well as eukaryotic cells in their folded states as a part of physiological and/or pathophysiological processes. Herein, we provide an overview of the systems/processes that are established or likely to involve the membrane translocation of folded proteins, such as protein export by the twin-arginine translocation system in bacteria and chloroplasts, unconventional protein secretion and protein import into the peroxisome in eukaryotes, and the cytosolic entry of proteins (e.g., bacterial toxins) and viruses into eukaryotes. We also discuss the various mechanistic models that have previously been proposed for the membrane translocation of folded proteins including pore/channel formation, local membrane disruption, membrane thinning, and transport by membrane vesicles. Finally, we introduce a newly discovered vesicular transport mechanism, vesicle budding and collapse, and present evidence that vesicle budding and collapse may represent a unifying mechanism that drives some (and potentially all) of folded protein translocation processes.