Dynamic Monte Carlo simulations of a new lattice model of globular protein folding, structure and dynamics

A long-standing problem of molecular biology is the prediction of globular protein tertiary structure from the primary sequence. In the context of a new, 24-nearest-neighbor lattice model of proteins that includes both alpha and beta-carbon atoms, the requirements for folding to a unique four-member beta-barrel, four-helix bundles and a model alpha/beta-bundle have been explored. A number of distinct situations are examined, but the common requirements for the formation of a unique native conformation are tertiary interactions plus the presence of relatively small (but not irrelevant) intrinsic turn preferences that select out the native conformer from a manifold of compact states. When side-chains are explicitly included, there are many conformations having the same or a slightly greater number of side-chain contacts as in the native conformation, and it is the local intrinsic turn preferences that produce the conformational selectivity on collapse. The local preference for helix or beta-sheet secondary structure may be at odds with the secondary structure ultimately found in the native conformation. The requisite intrinsic turn populations are about 0.3% for beta-proteins, 2% for mixed alpha/beta-proteins and 6% for helix bundles. In addition, an idealized model of an allosteric conformational transition has been examined. Folding occurs predominantly by a sequential on-site assembly mechanism with folding initiating either at a turn or from an isolated helix or beta-strand (where appropriate). For helical and beta-protein models, similar folding pathways were obtained in diamond lattice simulations, using an entirely different set of local Monte Carlo moves. This argues strongly that the results are universal; that is, they are independent of lattice, protein model or the particular realization of Monte Carlo dynamics. Overall, these simulations demonstrate that the folding of all known protein motifs can be achieved in the context of a single class of lattice models that includes realistic backbone structures and idealized side-chains.     

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.     

Organization of regions with amphiphilic alpha-helical potential within the three-dimensional structure of beta-sheet proteins.

The observation that strong amphiphilic alpha-helical potential exists in all proteins, including beta-sheet proteins, has given rise to the idea that alpha-helical intermediates may be critical to the folding paths of all proteins. Here we report that regions with amphiphilic alpha-helical potential in beta-sheet proteins are regularly spaced within the native structure of the proteins at an average interval of about 13 A. This regular spacing did not occur when the location of amphiphilic regions was randomly assigned (p = 0.0056), suggesting some degree of organization with respect to the native fold. However, in the native structure of various non-homologous proteins that contain the same fold, the location of the regions with amphiphilic alpha-helical potential was not conserved. Further, there was no apparent association of amphiphilic alpha-helical potential with any particular type of secondary structure, confirming that this potential is not involved in maintenance of native structure and suggesting that it may be associated with a highly adaptable process.     

Microsecond hydrophobic collapse in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein

Using small-angle X-ray scattering combined with a continuous-flow mixing device, we monitored the microsecond compaction dynamics in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein. A significant collapse of the radius of gyration from 30 A to 23.2 A occurs within 300 micros after the initiation of refolding by a urea dilution jump. The subsequent folding after the major chain collapse occurs on a considerably longer time-scale. The protein folding trajectories constructed by comparing the development of the compactness and the secondary structure suggest that the specific hydrophobic collapse model rather than the framework model better explains the experimental observations. The folding trajectory of this alpha/beta-type protein is located between those of alpha-helical and beta-sheet proteins, suggesting that native structure determines the folding landscape.     

Conformational propagation with prion-like characteristics in a simple model of protein folding

Protein refolding/misfolding to an alternative form plays an aetiologic role in many diseases in humans, including Alzheimer's disease, the systemic amyloidoses, and the prion diseases. Here we have discovered that such refolding can occur readily for a simple lattice model of proteins in a propagatable manner without designing for any particular alternative native state. The model uses a simple contact energy function for interactions between residues and does not consider the peculiarities of polypeptide geometry. In this model, under conditions where the normal (N) native state is marginally stable or unstable, two chains refold from the N native state to an alternative multimeric energetic minimum comprising a single refolded conformation that can then propagate itself to other protein chains. The only requirement for efficient propagation is that a two-faced mode of packing must be in the ground state as a dimer (a higher-energy state for this packing leads to less efficient propagation). For random sequences, these ground-state dimeric configurations tend to have more beta-sheet-like extended structure than almost any other sort of dimeric ground-state assembly. This implies that propagating states (such as for prions) are beta-sheet rich because the only likely propagating forms are beta-sheet rich. We examine the details of our simulations to see to what extent the observed properties of prion propagation can be predicted by a simple protein folding model. The formation of the alternative state in the present model shows several distinct features of amyloidogenesis and of prion propagation. For example, an analog of the phenomenon of conformationally distinct strains in prions is observed. We find a parallel between 'glassy' behavior in liquids and the formation of a propagatable state in proteins. This is the first report of simulation of conformational propagation using any heteropolymer model. The results imply that some (but not most) small protein sequences must maintain a sequence signal that resists refolding to propagatable alternative native states and that the ability to form such states is not limited to polypeptides (or reliant on regular hydrogen bonding per se) but can occur for other protein-like heteropolymers.     

Influence of secondary structure of polyglycine on some conformational properties

Values of four conformational properties, namely unperturbed dimension [r2]0, dipole moment [mu 2], mean squared optical anisotropy [gamma 2], and molar Kerr constant [mK], have been calculated for polyglycine chains allowing several combinations of the secondary structure with the aim of studying the dependence of these magnitudes on the secondary structure of the chain. Two different approaches to the secondary structure have been used. In the first, chains with all their units in a given conformation (random coil, alpha-helix or beta-sheet) are interrupted at several positions by one unit in a different conformation. In the second, chains with varying composition of two conformations alpha-helix/beta-sheet and beta-sheet/random coil were allowed and the results obtained compared with previous work for alpha-helix/random coil chains.     

The solution structure of concanavalin A probed by FT-IR spectroscopy

The secondary structural properties of various forms of concanavalin A in solution were investigated by Fourier-transform infrared spectroscopy in the Amide I region. As in the crystal, the solution structure of the native protein consists mainly of antiparallel beta-sheet. Carbohydrate binding does not produce major changes in the overall secondary structure of concanavalin A, but affects infrared bands due to loops and beta-turns. Upon demetallization, the spectrum of concanavalin A shows only a small change in the Amide I band, indicating that whereas the beta-sheet structure is conserved, the tertiary properties may be altered. There are also changes in the bands from the tyrosine residues which are compatible with local changes in structure. Confirming tertiary structural differences, the cation-depleted apoprotein is much less stable, denaturing around 63 degrees C, while the native protein denatures only at temperatures around 85 degrees C. Tetramerization proceeds without significant secondary structural change. However, aggregation of the tetramers leads to a significant decrease of the bands corresponding to beta-sheet structure, and changes in the tyrosine bands.     

A new protein folding algorithm based on hydrophobic compactness: Rigid Unconnected Secondary Structure Iterative Assembly (RUSSIA). II: Applications

The RUSSIA procedure (Rigid Unconnected Secondary Structure Iterative Assembly) produces structural models of cores of small- and medium-sized proteins. Loops are omitted from this treatment and regular secondary structures are reduced to points, the centers of their hydrophobic faces. This methodology relies on the maximum compactness of the hydrophobic residues, as described in detail in Part I. Starting data are the sequence and the predicted limits and natures of regular secondary structures (alpha or beta). Helices are treated as rigid cylinders, whereas beta-strands are collectively taken into account within beta-sheets modeled by helicoid surfaces. Strands are allowed to shift along their mean axis to allow some flexibility and the alpha-helices can be placed on either side of beta-sheets. Numerous initial conformations are produced by discrete rotations of the helices and sheets around the direction going from the center of their hydrophobic face to the global center of the protein. Selection of proposed models is based upon a criterion lying on the minimization of distances separating hydrophobic residues belonging to different regular secondary structures. The procedure is rapid and appears to be robust relative to the quality of starting data (nature and length of regular secondary structures). This dependence of the quality of the model on secondary structure prediction and in particular the beta-sheet topology, is one of the limits of the present algorithm. We present here some results for a set of 12 proteins (alpha, beta and alpha/beta classes) of lengths 40-166 amino acids. The r.m.s. deviations for core models with respect to the native proteins are in the range 1.4-3.7 A.     

Vibrational circular dichroism studies of epidermal growth factor and basic fibroblast growth factor.

Vibrational circular dichroism (VCD) studies are reported for two unrelated recombinant growth factor proteins: epidermal growth factor and basic fibroblast growth factor (bFGF). NMR, electronic CD, and bFGF X-ray studies indicate that these two proteins are primarily composed of beta-sheet and loop secondary structure elements with no detectable alpha-helices. Two reports on solution conformation of these proteins using FTIR absorption spectroscopy with subsequent resolution enhancement confirmed the presence of a large fraction of a beta-sheet conformation but in addition indicated the presence of large absorption bands in the 1650-1656 cm-1 region, which are typically assigned to alpha-helices. The VCD spectra of both proteins have band shapes that strongly resemble those of other high beta-sheet fraction proteins, such as the trypsin family of proteins. Quantitative analysis of the VCD spectra also indicates that these proteins are predominantly in beta-sheet and extended ("other") conformations with very little alpha-helix fraction. These results agree with the CD interpretation and affirm that the FTIR peaks in the region 1650-1656 cm-1 can be assigned to loops. This study provides an example of the limitations of using FTIR frequencies alone for examination of protein secondary structure.     

Structural composition of betaI- and betaII-proteins

Circular dichroism spectra of proteins are sensitive to protein secondary structure. The CD spectra of alpha-rich proteins are similar to those of model alpha-helices, but beta-rich proteins exhibit CD spectra that are reminiscent of CD spectra of either model beta-sheets or unordered polypeptides. The existence of these two types of CD spectra for beta-rich proteins form the basis for their classification as betaI- and betaII-proteins. Although the conformation of beta-sheets is largely responsible for the CD spectra of betaI-proteins, the source of betaII-protein CD, which resembles that of unordered polypeptides, is not completely understood. The CD spectra of unordered polypeptides are similar to that of the poly(Pro)II helix, and the poly(Pro)II-type (P2) structure forms a significant fraction of the unordered conformation in globular proteins. We have compared the beta-sheet and P2 structure contents in beta-rich proteins to understand the origin of betaII-protein CD. We find that betaII-proteins have a ratio of P2 to beta-sheet content greater than 0.4, whereas for betaI-proteins this ratio is less than 0.4. The beta-sheet content in betaI-proteins is generally higher than that in betaII-proteins. The origin of two classes of CD spectra for beta-rich proteins appears to lie in their relative beta-sheet and P2 structure contents.     

Statistical analyses of the vibrational circular dichroism of selected proteins and relationship to secondary structures

The vibrational circular dichroism (VCD) spectra of 20 proteins dissolved in D2O are presented in the amide I' region. These data are decomposed into a linear combination of orthogonal subspectra generated by the principal component method of factor analysis, and the results for 13 of them are compared to their secondary structures as determined from X-ray crystallography. Factor analysis of the VCD yields six statistically significant subspectra that can be used to reproduce the spectra. Their coefficients can then be used to characterize a given protein. Comparison of cluster analyses of these VCD coefficients and of the secondary structure fractional coefficients from X-ray crystallography showed that proteins clustered in the VCD analysis were also clustered in the X-ray analysis. The relative fractions of alpha-helix and beta-sheet in the protein dominate the clustering in both data sets. Qualitative characterization of the secondary structure of a given protein is obtained from its clustering on the basis of spectral characteristics. A strong linear correlation was found between the coefficient of the second subspectrum and the alpha-helical fraction for the proteins studied. The second coefficient also correlated to the beta-sheet fraction, and the first coefficient weakly correlated to the fraction for "other". Subsequent multiple-parameter regression analyses of the VCD factor analysis coefficients, constrained to include only significant dependencies, yielded reliable determination of the alpha-helix fraction and somewhat less confident determination of beta-sheet, bend, and "other" components. Predictive capability for proteins not in the regression was good. Varimax rotation of the coefficients transformed the subspectra and gave simple correlations to secondary structure components but had less reliability and more restrictions than the multiple regression on the original coefficients. The partial least-squares analysis method was also used to predict fractional secondary structures for the training set proteins but resulted in somewhat higher average error, particularly for beta-sheet, than the multiple regression. The turn fraction was effectively undetermined in both the regression and partial least-squares analyses. These statistical analyses represent the first determination of a quantitative relationship between VCD spectra and secondary structure in proteins.     

Lattice simulations of cotranslational folding of single domain proteins

We use lattice protein models and Monte Carlo simulations to study cotranslational folding of small single domain proteins. We show that the assembly of native structure begins during late extrusion stages, but final formation of native state occurs during de novo folding, when all residues are extruded. There are three main results in our study. First, for the sequences displaying two-state refolding mechanism de novo cotranslational folding pathway differs from that sampled in in vitro refolding. The change in folding pathways is due to partial assembly of native interactions during extrusion that results in different starting conditions for in vitro refolding and for de novo cotranslational folding. For small single domain proteins cotranslational folding is slower than in vitro refolding, but is generally fast enough to be completed before the release from a ribosome. Second, we found that until final stages of biosynthesis cotranslational folding is essentially equilibrium. This observation is explained by low stability of structured states for partially extruded chains. Finally, our data suggest that the proteins, which refold in vitro slowly via intermediates, complete their de novo folding after the release from a ribosome. Comparison of our lattice cotranslational simulations with recent experimental and computational studies is discussed.     

Probing local secondary structure by fluorescence: time-resolved and circular dichroism studies of highly purified neurotoxins.

The relationship between beta-sheet secondary structure and intrinsic tryptophan fluorescence parameters of erabutoxin b, alpha-cobratoxin, and alpha-bungarotoxin were examined. Nuclear magnetic resonance and x-ray crystallography have shown that these neurotoxins have comparable beta-sheet, beta-turn, and random coil secondary structures. Each toxin contains a single tryptophan (Trp) residue within its beta-sheet. The time-resolved fluorescence properties of native erabutoxin b and alpha-cobratoxin are best described by triple exponential decay kinetics, whereas native alpha-bungarotoxin exhibits more than four lifetimes. The disulphide bonds of each toxin were reduced to facilitate carboxymethylation and amidocarboxymethylation. The two different toxin derivatives of all three neurotoxins displayed triple exponential decay kinetics and were completely denatured as evidenced by circular dichroism (random coil). The concentration (c) values of the three fluorescence decay times (time-resolved fluorescence spectroscopy (TRFS)) were dramatically different from those of the native toxins. Each neurotoxin, treated with different concentrations of guanidinium hydrochloride (GuHCl), was studied both by circular dichroism and TRFS. Disappearance of the beta-sheet secondary structural features with increasing concentrations of GuHCl was accompanied by a shift in the relative contribution (c value) of each fluorescence decay time (TRFS). It was found that certain disulphide residues confer added stability to the beta-sheet secondary structure of these neurotoxins and that the center of the beta-sheet is last to unfold. These titrations show that Trp can be used as a very localized probe of secondary structure.     

Unraveling double stranded alpha-helical coiled coils: an x-ray diffraction study on hard alpha-keratin fibers

Transformations of proteins secondary and tertiary structures are generally studied in globular proteins in solution. In fibrous proteins, such as hard alpha-keratin, that contain long and well-defined double stranded alpha-helical coiled coil domains, such study can be directly done on the native fibrous tissue. In order to assess the structural behavior of the coiled coil domains under an axial mechanical stress, wide angle x-ray scattering and small angle x-ray scattering experiments have been carried out on stretched horse hair fibers at relative humidity around 30%. Our observations of the three major axial spacings as a function of the applied macroscopic strain have shown two rates. Up to 4% macroscopic strain the coiled coils were slightly distorted but retained their overall conformation. Above 4% the proportion of coiled coil domains progressively decreased. The main and new result of our study is the observation of the transition from alpha-helical coiled coils to disordered chains instead of the alpha-helical coiled coil to beta-sheet transition that occurs in wet fibers.     

Proteins can convert to beta-sheet in single crystals

Raman microscopy was used to follow conformational changes in single protein crystals. Crystals of native insulin and of the 5S and 12S subunits of the enzyme transcarboxylase showed a mixture of Raman marker bands signifying alpha-helix, beta-sheet and nonordered secondary structure. However, by reducing the S-S bonds in the insulin crystal, or by lowering the pH for the 5S crystal, or by soaking substrates into the 12S crystal, the secondary structure in each crystal became predominantly beta-sheet. The beta-form crystals could be dissolved only with difficulty and yielded high-molecular weight protein aggregates, indicating that the beta-sheet formation involves intermolecular contacts. Although their morphology appeared unchanged, the crystals no longer diffracted X-rays. Using crystals that had not been exposed to laser light, the dye thioflavin T formed highly fluorescent complexes with the "beta-transformed" crystals but not with the native crystals.     

Conformational transition between native and reactive center cleaved forms of alpha 1-antitrypsin by Fourier transform infrared spectroscopy and small-angle neutron scattering

alpha 1-Antitrypsin (alpha 1-AT) is the best-characterized member of the serpin superfamily of plasma proteins. Protease inhibitor members of this family undergo a characteristic reactive-center cleavage during expression of their inhibitory activity. The physical basis of this transition in alpha 1-AT from the stressed native conformation to the more stable reactive center cleaved (split) form was studied by Fourier transform infrared (FT-IR) spectroscopy and neutron scattering. The FT-IR spectra show that, while split alpha 1-AT has three intense well-resolved components associated with the presence of antiparallel beta-sheet and alpha-helix conformations, the amide I band of native alpha 1-AT has only one intense component, associated with the presence of beta-sheet structure. 1H-2H exchange within the polypeptide backbone, studied by FT-IR and NMR spectroscopy, shows that the native form undergoes greater exchange than the split form. Under the same conditions, neutron scattering shows no differences in the radius of gyration RG of the native and the split forms. In contrast, in high concentrations of phosphate approaching those used for crystallization, the native form (unlike the split form) undergoes dimerization. These data indicate that the conformational transition largely involves localized secondary and tertiary structure rearrangements. We propose that the energetically stressed native alpha 1-AT structure is the consequence of a significantly reduced number of hydrogen bonds in secondary structure components and that reactive-site cleavage between Met358 and Ser359 is the key for the development of the fully hydrogen bonded more stable serpin structure.     

Study of the "molten globule" intermediate state in protein folding by a hydrophobic fluorescent probe

Binding of the hydrophobic fluorescent probe, 1-anilino-naphthalene-8-sulfonate (ANS), to synthetic polypeptides and proteins with a different structural organization has been studied. It has been shown that ANS has a much stronger affinity to the protein "molten globule" state, with a pronounced secondary structure and compactness, but without a tightly packed tertiary structure as compared with its affinity to the native and coil-like proteins, or to coil-like, alpha-helical, or beta-structural hydrophilic homopolypeptides. The possibility of using ANS for the study of equilibrium and kinetic molten globule intermediates is demonstrated, with carbonic anhydrase, beta-lactamase, and alpha-lactalbumin as examples.     

Prediction of secondary structure, spatial organization and distribution of antigenic determinants for hepatitis A virus proteins

On the basis of the secondary structure calculations from the known amino acid sequence we came to the conclusion that hepatitis A virus capsid proteins have the typical antiparallel beta-sheet bilayer structure. The predicted secondary structure of the HAV proteins can be well aligned with those of the poliovirus (type 1 Mahoney) and human rhinovirus (type 14). It enabled us to use the X-ray structure of the PV-1M and HRV-14 proteins as a template and then, firstly, to localize the positions of alpha and beta regions in the architecture of the HAV protein molecules and, secondly, to discover the amino acid homologies of the secondary structure regions aligned. The obtained model of the three-dimensional structure for HAV proteins helped us to indicate the exposed regions of the polypeptide chains and to pinpoint the potential neutralizing antigenic sites.     

Prediction of Secondary Structure,Spatial Organization and Distribution of Antigenic Determinants for Hepatitis A Virus Proteins

Abstract

On the basis of the secondary structure calculations from the known amino acid sequence we came to the conclusion that hepatitis A virus capsid proteins have the typical antiparallel beta-sheet bilayer structure.

The predicted secondary structure of the HAV proteins can be well aligned with those of the poliovirus (type 1 Mahoney) and human rhinovirus (type 14). It enabled us to use the X-ray structure of the PV-1M and HRV-14 proteins as a template and then, firstly, to localize the positions of alpha and beta regions in the architecture of the HAV protein molecules and, secondly, to discover the amino acid homologies of the secondary structure regions aligned. The obtained model of the three-dimensional structure for HAV proteins helped us to indicate the exposed regions of the polypeptide chains and to pinpoint the potential neutralizing antigenic sites.