Structural principles of the globular organization of protein chains. A stereochemical theory of globular protein secondary structure

The paper reveals the types of amino acid sequences of polypeptide chain regions of globular protein which form a regular (α or β) or irregular conformation in the native globule. The study was made taking into account general “architectural” principles of packing of polypeptide chains in globular proteins and considering the interactions of proteins with water molecules. An a priori theory is developed which permits the identification, in good agreement with experiment, of α-helical and β-structural regions in globular proteins from their primary structure.     

A simple and novel interpretation of the three-dimensional structure of globular proteins based on quantum-mechanical computations on small model molecules. I

It is suggested that the three-dimensional structure of globular proteins is partly determined by a framework of strengthened hydrogen bonds that involves both ionic side chains and water molecules in addition to the polypeptide backbone. This conclusion follows from a combination of the results of ab initio molecular-orbital computations on small model molecules and high-accuracy x-ray data on the rubredoxin molecule. The computations yield the idea of hydrogen-bonded bridges that are built from tens of atoms, and the experimental information yields the idea that the bridges are assembled into clusters, each of which is built from hundreds of atoms. Some 10 such clusters then form a globular protein.     

Theory for the folding and stability of globular proteins

Using lattice statistical mechanics, we develop theory to account for the folding of a heteropolymer molecule such as a protein to the globular and soluble state. Folding is assumed to be driven by the association of solvophobic monomers to avoid solvent and opposed by the chain configurational entropy. Theory predicts a phase transition as a function of temperature or solvent character. Molecules that are too short or too long or that have too few solvophobic residues are predicted not to fold. Globular molecules should have a largely solvophobic core, but there is an entropic tendency for some residues to be "out of place", particularly in small molecules. For long chains, molecules comprised of globular domains are predicted to be thermodynamically more stable than spherical molecules. The number of accessible conformations in the globular state is calculated to be an exceedingly small fraction of the number available to the random coil. Previous estimates of this number, which have motivated kinetic theories of folding, err by many tens of orders of magnitude.     

Polymeric aspects of protein folding: a brief overview

Regardless of the differences in primary amino acid sequences, protein molecules in a number of conformational states behave as polymer homologues, allowing speculations as to the volume interactions being a driving force in formation of equilibrium structures. For instance, both native and molten globules exhibit key features of polymer globules, where the fluctuations of the molecular density are expected to be much less than the molecular density itself. Protein molecules in the compact denatured (pre-molten globule) states possess properties of squeezed coils. In fact, even high concentrations of strong denaturants (e.g., urea and GdmCl) more likely constitute bad solvents for protein chains. Thus, globular proteins are probably never random coils without positional correlations and biological polypeptide chains represent the macromolecular coils below a critical point even under harsh denaturing conditions. Several implications of these findings to protein folding are discussed.     

Thermodynamics and kinetics of protein folding: an evolutionary perspective

This article appeals to an evolutionary model which postulates that primordial proteins were described by small polypeptide chains which (i) lack disulfide bridges, and (ii) display slow folding rates with multi-state kinetics, to determine relations between structural properties of proteins and their folding kinetics. We parameterize the energy landscape of proteins in terms of thermodynamic activation variables. The model studies evolutionary changes in these thermodynamic parameters, and we invoke relations between these activation variables and structural properties of the protein to predict the following correspondence between protein structure and folding kinetics. 1. Proteins with inter- and intra-chain disulfide bridges: large variability in both folding rates and stability of intermediates, multi-state kinetics. 2. Proteins which lack inter and intra-chain disulfide bridges. 2.1 Single-domain chains: fast folding rates; unstable intermediates; two-state kinetics. 2.2 Multi-domain monomers: intermediate rates; metastable intermediates; multi-state kinetics. 2.3 Multi-domain oligomers: slow rates; metastable intermediates; multi-state kinetics. The evolutionary model thus provides a kinetic characterization of one important subfamily of proteins which we describe by the following properties: Folding dynamics of single-domain proteins which lack disulfide bridges are described by two-state kinetics. Folding rate of this class of proteins is positively correlated with the thermodynamic stability of the folded state.     

Why do proteins divide into domains? Insights from lattice model simulations

It is known that larger globular proteins are built from domains, relatively independent structural units. A domain size seems to be limited, and a single domain consists of from few tens to a couple of hundred amino acids. Based on Monte Carlo simulations of a reduced protein model restricted to the face centered simple cubic lattice, with a minimal set of short-range and long-range interactions, we have shown that some model sequences upon the folding transition spontaneously divide into separate domains. The observed domain sizes closely correspond to the sizes of real protein domains. Short chains with a proper sequence pattern of the hydrophobic and polar residues undergo a two-state folding transition to the structurally ordered globular state, while similar longer sequences follow a multistate transition. Homopolymeric (uniformly hydrophobic) chains and random heteropolymers undergo a continuous collapse transition into a single globule, and the globular state is much less ordered. Thus, the factors responsible for the multidomain structure of proteins are sufficiently long polypeptide chain and characteristic, protein-like, sequence patterns. These findings provide some hints for the analysis of real sequences aimed at prediction of the domain structure of large proteins.     

Protein sequences yield a proteomic code

Analysis of crystallized protein structures suggests that globular proteins are organized as consecutively connected units of 25-35 residues. These units are closed loops, that is returns of the polypeptide chain trajectory to a close contact with itself. This universal feature of apparently polymer-statistical nature is a basis for a principally novel view on the globular proteins as loop fold structures. The same unit size has been detected in protein sequences translated from complete prokaryotic genomes by positional autocorrelation analysis, which strongly indicates the evolutionary connection of the units. The units are further characterized by prototype sequences matching to their numerous derivatives in the translated genomes. The matches to five strongest prokaryotic prototypes and three prototypes of C. elegans are identified in the sequences of crystallized proteins, and their structures analyzed. Corresponding segments of the polypeptide chains in majority of cases form closed loops, though evolutionary fate of every prototype element is shown to be rather diverse. Then loop ends can be separated by a sequence-wise distant segments and stabilized by the spatial interactions in the context of the overall globular structure. The units belong to a presumably limited spectrum of the sequence prototypes, full repertoire of which would constitute a proteomic code.     

Structure and hydrodynamic properties of plectin molecules

Plectin is a cytoskeletal, high molecular weight protein of widespread and abundant occurrence in cultured cells and tissues. To study its molecular structure, the protein was purified from rat glioma C6 cells and subjected to chemical and biophysical analyses. Plectin's polypeptide chains have an apparent molecular weight of 300,000, as shown by one-dimensional sodium dodecyl sulfate/polyacrylamide electrophoresis. Cross-linking of non-denatured plectin in solution with dimethyl suberimidate and electrophoretic analyses on sodium dodecyl sulfate/agarose gels revealed that the predominant soluble plectin species was a molecule of 1200 X 10(3) Mr consisting of four 300 X 10(3) Mr polypeptide chains. Hydrodynamic properties of plectin in solution were obtained by sedimentation velocity centrifugation and high-pressure liquid chromatography analysis yielding a sedimentation coefficient of 10 S and a Stokes radius of 27 nm. The high f/fmin ratio of 4.0 indicated a very elongated shape of plectin molecules and an axial ratio of about 50. Shadowing and negative staining electron microscopy of plectin molecules revealed multiple domains: a rigid rod of 184 nm in length and 2 nm in diameter, and two globular heads of 9 nm diameter at each end of the rod. Circular dichroism spectra suggested a composition of 30% alpha-helix, 9% beta-structure and 61% random coil or aperiodic structure. The rod-like shape, the alpha-helix content as well as the thermal transition within a midpoint of 45 degrees C and the transition enthalpy (168 kJ/mol) of secondary structure suggested a double-stranded, alpha-helical coiled coil rod domain. Based on the available data, we favor a model of native plectin as a dumb-bell-like association of four 300 X 10(3) Mr polypeptide chains. Electron microscopy and turbidity measurements showed that plectin molecules self-associate into various oligomeric states in solutions of nearly physiological ionic strength. These interactions apparently involved the globular end domains of the molecule. Given its rigidity and elongated shape, and its tendency towards self-association, plectin may well be an interlinking element of the cytoskeleton that may also form a network of its own.     

Novel protein science enabled by total chemical synthesis

Chemical synthesis is a well‐established method for the preparation in the research laboratory of multiple‐tens‐of‐milligram amounts of correctly folded, high purity protein molecules. Chemically synthesized proteins enable a broad spectrum of novel protein science. Racemic mixtures consisting of d ‐protein and l ‐protein enantiomers facilitate crystallization and determination of protein structures by X‐ray diffraction. d ‐Proteins enable the systematic development of unnatural mirror image protein molecules that bind with high affinity to natural protein targets. The d ‐protein form of a therapeutic target can also be used to screen natural product libraries to identify novel small molecule leads for drug development. Proteins with novel polypeptide chain topologies including branched, circular, linear‐loop, and interpenetrating polypeptide chains can be constructed by chemical synthesis. Medicinal chemistry can be applied to optimize the properties of therapeutic protein molecules. Chemical synthesis has been used to redesign glycoproteins and for the a priori design and construction of covalently constrained novel protein scaffolds not found in nature. Versatile and precise labeling of protein molecules by chemical synthesis facilitates effective application of advanced physical methods including multidimensional nuclear magnetic resonance and time‐resolved FTIR for the elucidation of protein structure–activity relationships. The chemistries used for total synthesis of proteins have been adapted to making artificial molecular devices and protein‐inspired nanomolecular constructs. Research to develop mirror image life in the laboratory is in its very earliest stages, based on the total chemical synthesis of d ‐protein forms of polymerase enzymes.     

Theoretical study of a landscape of protein folding-unfolding pathways. Folding rates at midtransition

This paper presents a new method for calculating the folding-unfolding rates of globular proteins. The method is based on solution of kinetic equations for a network of folding-unfolding pathways of the proteins. The rates are calculated in the point of thermodynamic equilibrium between the native and completely unfolded states. The method has been applied to all the proteins listed by Jackson [Jackson, S. E. (1998) Folding Des. 3, R81-R91] and some peptides. Although the studied protein chains differ by more than 1 order of magnitude in size and exhibit two- as well as three-state kinetics in water, and their folding rates cover more than 11 orders of magnitude, the theoretical estimates are reasonable close to the experimentally measured folding rates in midtransition (the correlation coefficient being as high as 0.78). This means that the presented theory (having no adjustable parameters at all) is consistent with the experimental observations.     

Measuring “Unmeasurable” Folding Kinetics of Proteins by Single-Molecule Force Spectroscopy

Folding and unfolding are fundamental biological processes in cell and are important for the biological functions of proteins. Characterizing the folding and unfolding kinetics of proteins is important for understanding the energetic landscape leading to the active native conformations of these molecules. However, the thermal or chemical-induced unfolding of many proteins is irreversible in vitro, precluding characterization of the folding kinetics of such proteins, just as it is impossible to “un-boil” an egg. Irreversible unfolding often manifests as irreversible aggregation of unfolded polypeptide chains, which typically occurs between denatured protein molecules in response to the exposure of hydrophobic residues to solvent. An example of such a protein where thermal denaturation results in irreversible aggregation is the β-1,4 endoxylanase from Bacillus circulans (BCX). Here, we report the use of single-molecule atomic force microscopy to directly measure the folding kinetics of BCX in vitro. By mechanically unfolding BCX, we essentially allowed only one unfolded molecule to exist in solution at a given time, effectively eliminating the possibility for aggregation. We found that BCX can readily refold back to the native state, allowing us to measure its folding kinetics for the first time. Our results demonstrate that single-molecule force-spectroscopy-based methods can adequately tackle the challenge of “un-boiling eggs”, providing a general methodology to characterize the folding kinetics of many proteins that suffer from irreversible denaturation and thus cannot be characterized using traditional equilibrium methodologies.     

Dynamic polyhedral models of globular proteins

We have devised several mechanical models of globular proteins by approximating them to various polyhedra (dodecahedron, truncated octahedron, icosahedron, truncated icosahedron). The models comprise hollow blocks linked together in a flexible chain. Between blocks there is a set of several reversible, weak magnetic interactions such that when the chain is agitated, it will fold into a stable polyhedral structure about the size of a hand. Folding may be followed in real time with a video camera. Key to the success of the folding process is the lightness of the chain. Several side chains may also be added to the blocks such that they come together to create a polyhedral core when the chain folds. The models have a number of similarities to globular proteins: each chain folds into a unique, but dynamic, three-dimensional structure; the instructions that determine this structure are built into the configuration of blocks; and it is difficult to predict this structure given the unfolded block configuration. Furthermore, the chains fold quickly, generally in less than a minute, several pathways are involved, and these pathways progress through elements of "native" structure. In particular, the models emphasize the importance of restricted conformational mobility in assisting the chain to fold, and also in eliminating undesirable interactions. Because of these similarities to globular proteins, we believe that the polyhedral models will, with continued development, be helpful in understanding the protein folding process, while at the same time acting as valuable educational visual aids. They might also inspire the construction of new types of microscopic, self-assembling devices.     

Insight into the structure of amyloid fibrils from the analysis of globular proteins

The conversion from soluble states into cross-β fibrillar aggregates is a property shared by many different proteins and peptides and was hence conjectured to be a generic feature of polypeptide chains. Increasing evidence is now accumulating that such fibrillar assemblies are generally characterized by a parallel in-register alignment of β-strands contributed by distinct protein molecules. Here we assume a universal mechanism is responsible for β-structure formation and deduce sequence-specific interaction energies between pairs of protein fragments from a statistical analysis of the native folds of globular proteins. The derived fragment–fragment interaction was implemented within a novel algorithm, prediction of amyloid structure aggregation (PASTA), to investigate the role of sequence heterogeneity in driving specific aggregation into ordered self-propagating cross-β structures. The algorithm predicts that the parallel in-register arrangement of sequence portions that participate in the fibril cross-β core is favoured in most cases. However, the antiparallel arrangement is correctly discriminated when present in fibrils formed by short peptides. The predictions of the most aggregation-prone portions of initially unfolded polypeptide chains are also in excellent agreement with available experimental observations. These results corroborate the recent hypothesis that the amyloid structure is stabilised by the same physicochemical determinants as those operating in folded proteins. They also suggest that side chain–side chain interaction across neighbouring β-strands is a key determinant of amyloid fibril formation and of their self-propagating ability.     

Observations on the gel chromatography of collagen-like peptides and use of a simple calibration method for molecular weight determination.

Elution characteristics of collagen-derived polypeptides and of globular proteins were compared under identical experimental conditions with agarose gels. This comparison permitted calculation of the hydrodynamic radii of collagen polypeptide chains of different molecular weight, and these radii were shown to be in reasonable agreement with estimates made from intrinsic viscosity data. Two distinet linear relationships were observed for collagen polypeptide chains, relating logarithm of molecular weight to elution parameters. Peptide chains of MW 3300 and lower fell on a line of a steeper slope than did larger polypeptide chains.A simple procedure for molecular weight estimation of an unknown polypeptide chain of the collagen class is described, using only three commercially available standards for calibration: reduced, carboxymethylated Ascaris cuticle collagen, and the synthetic peptides (l-Pro-l-Pro-Gly)₁₀ and (l-Pro-L-Pro-Gly)₅.     

Correlations between internal mobility and stability of globular proteins.

The recent work is surveyed which leads to the suggestions that the conformation of globular proteins in solution corresponds to a dynamic ensemble of rapidly interconverting spatial structures, that clusters of hydrophobic amino acid side chains have an important role in the architecture of protein molecules, and that mechanistic aspects of protein denaturation can be correlated with internal mobility seen in the native conformation. These conclusions resulted originally from high resolution 1H nuclear magnetic resonance (NMR) studies of aromatic ring mobility, exchange of interior amide protons and thermal denaturation of the basic pancreatic trypsin inhibitor and a group of related proteins. Various new approaches to further characterize proteins in solution have now been taken and preliminary data are presented. These include computer graphics to outline hydrophobic clusters in globular protein structures, high resolution 1H-NMR experiments at variable hydrostatic pressure and 13C-NMR relaxation measurements. At the present early stage of these new investigations it appears that the hydrophobic cluster model for globular proteins is compatible with the data obtained.     

Purification and properties of the membrane-bound form of acetylcholinesterase from chicken brain. Evidence for two distinct polypeptide chains

The major molecular form of acetylcholinesterase (AChE) from chicken brain is a membrane-bound glycoprotein with an apparent sedimentation coefficient of 11.4 S. Analysis of the purified protein by gel filtration, velocity sedimentation, and sodium dodecyl sulfate-gel electrophoresis shows that the solubilized enzyme is a globular tetramer with an apparent Mr = 420,000. This membrane-bound form of AChE is hydrophobic and readily aggregates in the absence of detergent. These aggregates are concentration-dependent, relatively stable in the presence of high salt concentrations, yet readily dissociate upon addition of detergent to the 11.4 S form, indicating that the interactions are hydrophobic. Polyclonal and monoclonal antibodies raised against chicken brain AChE purified by ion exchange chromatography, affinity chromatography, and preparative gel electrophoresis precipitate AChE enzyme activity. However, these antibodies do not cross-react with the enzyme from chicken muscle which preferentially hydrolyses butyrylcholine. Immunoprecipitation of isotopically labeled enzyme molecules from tissue cultured brain cells and analysis by sodium dodecyl sulfate-gel electrophoresis shows that AChE consists of two polypeptide chains with apparent Mr = 105,000 (alpha) and 100,000 (beta) in a 1:1 ratio. Immunoblotting of brain AChE with either the polyclonal or monoclonal antibodies indicates that the alpha and beta chains share antigenic determinants. Furthermore, both polypeptide chains can be labeled with [3H]diisopropyl fluorophosphate, indicating that they each contain a catalytic site. This is the first indication that globular forms of AChE may consist of multiple polypeptide chains.     

FoldUnfold: web server for the prediction of disordered regions in protein chain

Identification of disordered regions in polypeptide chains is very important because such regions are essential for protein function. A new parameter, namely mean packing density of residues has been introduced to detect disordered regions in a protein sequence. We have demonstrated that regions with weak expected packing density would be responsible for the appearance of disordered regions. Our method (FoldUnfold) has been tested on datasets of globular proteins (559 proteins) and long disordered protein segments (129 proteins) and showed improved performance over some other widely used methods, such as DISOPRED, PONDR VL3H, IUPred and GlobPlot. AVAILABILITY: The FoldUnfold server is available for users at http://skuld.protres.ru/~mlobanov/ogu/ogu.cgi. There is a link to our server through the web site of DisProt (http://www.disprot.org/predictors.php).     

Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

Novel forms of chemical protein diversity -- in nature and in the laboratory

Novel chemical variants of proteins have been found in nature, including potent 'microprotein' natural products and folded protein molecules that contain a cyclic polypeptide chain. Researchers have used chemical synthesis and genetic methods to make these proteins and more: protein catenanes, neoglycoproteins, and artificial protein molecules with novel architectures or made from novel building blocks. De novo design has taken a big step forward with the accurate design and construction of proteins with complex molecular structure. A variety of non-coded amino acids and other building blocks has been used to make increasingly sophisticated protein molecular devices for use as biosensors and for the study of signal transduction inside living cells.