Characterizing the conformational ensemble of monomeric polyglutamine

Studies of synthetic polyglutamine peptides in vitro have established that polyglutamine peptides aggregate via a classic nucleation and growth mechanism. Chen and colleagues [Proc Natl Acad Sci U S A 2002;99:11884-11889] have found that monomeric polyglutamine, which is a disordered statistical coil in solution, is the critical nucleus for aggregation. Therefore, nucleation of beta-sheet-rich aggregates requires an initial disorder to order transition, which is a highly unfavorable thermodynamic reaction. The questions of interest to us are as follows: What are the statistical fluctuations that drive beta-sheet formation in monomeric polyglutamine? How do these fluctuations vary with chain length? And why is this process thermodynamically unfavorable, that is, why is monomeric polyglutamine disordered? To answer these questions we use multiple molecular dynamics simulations to provide quantitative characterization of conformational ensembles for two short polyglutamine peptides. We find that the ensemble for polyglutamine is indeed disordered. However, the disorder is inherently different from that of denatured proteins and the average compactness and magnitude of conformational fluctuations increase with chain length. Most importantly, the effective concentration of sidechain primary amides around backbone units is inherently high and peptide units are solvated either by hydrogen bonds to sidechains or surrounding water molecules. Due to the multiplicity of backbone solvation modes the probability associated with any specific backbone conformation is small, resulting in a conformational entropy bottleneck which makes beta-sheet formation in monomeric polyglutamine thermodynamically unfavorable.     

Crystal architecture and conformational properties of the inclusion complex, neohesperidin dihydrochalcone-cyclomaltoheptaose (beta-cyclodextrin), by X-ray diffraction

The crystal structure of the host-guest noncovalent complex of cyclomaltoheptaose (beta-cyclodextrin, betaCD) with the O-diglycosyl flavonoid neohesperidin dihydrochalcone [(3,5-dihydroxy-4-(3-hydroxy-4-methoxyhydrocinnamoyl)phenyl-2-O-(alpha-L-rhamnopyranosyl)-beta-D-glucopyranoside, NDC] has been determined from single-crystal X-ray diffraction data collected at low temperature (130 K), using synchrotron radiation. The crystal data are as follows: a =15.125(5), b =30.523(5), c =41.332(5) Angstroms, orthorhombic, space group C222(1). The structure contains 19 molecules of water, of which 11 appeared well positioned, whereas 9 are disordered over 23-positions. The betaCD-NDC complex is characterized by one aromatic part of NDC deeply inserted into the hydrophobic cavity of the betaCD through the primary OH rim, and it is present in the crystal as a dimer. The dimeric units, formed by head-to-head assemblies of CD molecules, each with its guest, are self-assembled in columns. The stability of the columns is provided by host-guest and guest-guest attractive interactions, thus showing a key role of the guest molecules in the crystal architecture. The guest conformation in the complex is different from that reported in the literature for uncomplexed NDC. The host-induced conformational changes on NDC provide the optimum geometry requirements for the assembly of the dimeric units.     

Coarse-grained simulation studies of peptide-induced pore formation

We investigate the interactions between lipid bilayers and amphiphilic peptides using a solvent-free coarse-grained simulation technique. In our model, each lipid is represented by one hydrophilic and three hydrophobic beads. The amphiphilic peptide is modeled as a hydrophobic-hydrophilic cylinder with hydrophilic caps. We find that with increasing peptide-lipid attraction the preferred state of the peptide changes from desorbed, to adsorbed, to inserted. A single peptide with weak attraction binds on the bilayer surface, while one with strong attraction spontaneously inserts into the bilayer. We show how several peptides, which individually bind only to the bilayer surface, cooperatively insert. Furthermore, hydrophilic strips along the peptide cylinder induce the formation of multipeptide pores, whose size and morphology depend on the peptides’ overall hydrophilicity, the distribution of hydrophilic residues, and the peptide-peptide interactions. Strongly hydrophilic peptides insert less readily, but prove to be more destructive to bilayer integrity.     

Membrane perturbation induced by interfacially adsorbed peptides

The structural and energetic characteristics of the interaction between interfacially adsorbed (partially inserted) alpha-helical, amphipathic peptides and the lipid bilayer substrate are studied using a molecular level theory of lipid chain packing in membranes. The peptides are modeled as "amphipathic cylinders" characterized by a well-defined polar angle. Assuming two-dimensional nematic order of the adsorbed peptides, the membrane perturbation free energy is evaluated using a cell-like model; the peptide axes are parallel to the membrane plane. The elastic and interfacial contributions to the perturbation free energy of the "peptide-dressed" membrane are evaluated as a function of: the peptide penetration depth into the bilayer's hydrophobic core, the membrane thickness, the polar angle, and the lipid/peptide ratio. The structural properties calculated include the shape and extent of the distorted (stretched and bent) lipid chains surrounding the adsorbed peptide, and their orientational (C-H) bond order parameter profiles. The changes in bond order parameters attendant upon peptide adsorption are in good agreement with magnetic resonance measurements. Also consistent with experiment, our model predicts that peptide adsorption results in membrane thinning. Our calculations reveal pronounced, membrane-mediated, attractive interactions between the adsorbed peptides, suggesting a possible mechanism for lateral aggregation of membrane-bound peptides. As a special case of interest, we have also investigated completely hydrophobic peptides, for which we find a strong energetic preference for the transmembrane (inserted) orientation over the horizontal (adsorbed) orientation.     

An accurate model of polyglutamine

Polyglutamine repeats in proteins are highly correlated with amyloid formation and neurological disease. To better understand the molecular basis of glutamine repeat diseases, structural analysis of polyglutamine peptides as soluble monomers, oligomers, and insoluble amyloid fibrils is necessary. In this study, fluorescence resonance energy transfer (FRET) experiments and molecular dynamics simulations using different theoretical models of polyglutamine were conducted. This study demonstrates that a previously proposed simple C(α)C(β) model of polyglutamine, denoted as FCO, accurately reproduced the present FRET results and the results of previously published FRET, triplet-state quenching, and fluorescence correlation studies. Other simple C(α)C(β) models with random coil and extended β-strand parameters, and all-atom models with parm96 and parm99SB force fields, did not match the FRET result well. The FCO is an intrinsically disordered model with a high-effective persistence length producing extended peptides at short lengths (Q(N) < 10). Because of an increasing number of attractive Q-Q interactions at longer lengths, the FCO model becomes increasingly more compact at lengths between Q(N) ～ 10-16 and is as compact as many folded proteins at Q(N) > 16.     

Tryptophan stabilizes His-heme loops in the denatured state only when it is near a loop end

We use a host-guest approach to evaluate the effect of Trp guest residues relative to Ala on the kinetics and thermodynamics of formation of His-heme loops in the denatured state of iso-1-cytochrome c at 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). Trp guest residues are inserted into an alanine-rich segment placed after a unique His near the N-terminus of iso-1-cytochrome c. Trp guest residues are either 4 or 10 residues from the His end of the 28-residue loop. We find the guest Trp stabilizes the His-heme loop at all GdnHCl concentrations when it is the 4th, but not the 10th, residue from the His end of the loop. Thus, residues near loop ends are most important in developing topological constraints in the denatured state that affect protein folding. In 1.5 M GdnHCl, the loop stabilization is ~0.7 kcal/mol, providing a thermodynamic rationale for the observation that Trp often mediates residual structure in the denatured state. Measurement of loop breakage rate constants, k(b,His), indicates that loop stabilization by the Trp guest residues occurs completely after the transition state for loop formation in 6.0 M GdnHCl. Under poorer solvent conditions, approximately half of the stabilization of the loop develops in the transition state, consistent with contacts in the denatured state being energetically downhill and providing evidence for funneling even near the rim of the folding funnel.     

An isolated helix persists in a sparsely populated form of KIX under native conditions

NMR relaxation dispersion techniques were used to investigate conformational exchange of the three-helix bundle protein KIX under native conditions. These experiments provide site-resolved kinetic information about microsecond-to-millisecond time scale motions along with structural (chemical shift) information without requiring a perturbation of the equilibrium. All kinetic data are consistent with an apparent two-state transition between natively folded KIX and a partially unfolded high-energy state that is populated to 3.0 +/- 0.2% at 27 degrees C. By combining (13)C- and (15)N-based experiments that probe specific structural aspects, we show that the sparsely populated high-energy state displays a strong conformational preference. An isolated secondary structural element, C-terminal helix alpha3, is highly populated, while the hydrophobic core of the domain and the remainder of the protein backbone, including helices alpha1 and alpha2, are disordered and devoid of specific interactions. This high-energy state presumably represents the equilibrium analogue of a folding intermediate that is transiently populated in stopped-flow kinetic experiments [Horng, J. C., Tracz, S. M., Lumb, K. J., and Raleigh, D. P. (2002) Biochemistry 44, 627-634].     

Order through disorder: hyper-mobile C-terminal residues stabilize the folded state of a helical peptide. a molecular dynamics study

Conventional wisdom has it that the presence of disordered regions in the three-dimensional structures of polypeptides not only does not contribute significantly to the thermodynamic stability of their folded state, but, on the contrary, that the presence of disorder leads to a decrease of the corresponding proteins' stability. We have performed extensive 3.4 μs long folding simulations (in explicit solvent and with full electrostatics) of an undecamer peptide of experimentally known helical structure, both with and without its disordered (four residue long) C-terminal tail. Our simulations clearly indicate that the presence of the apparently disordered (in structural terms) C-terminal tail, increases the thermodynamic stability of the peptide's folded (helical) state. These results show that at least for the case of relatively short peptides, the interplay between thermodynamic stability and the apparent structural stability can be rather subtle, with even disordered regions contributing significantly to the stability of the folded state. Our results have clear implications for the understanding of peptide energetics and the design of foldable peptides.     

Left-handed Z form in superhelical DNA: a theoretical study

This is a comprehensive statistical mechanical treatment of the Z form formation in purinepyrimidine stretches of different length inserted into superhelical DNA. The B-Z transition for short inserts is shown to follow the "all-or-none" principle. Over some critical value of the insert length n, the B-Z transition in the insert proceeds in two stages. The flipping of m base pairs into the Z form is followed by a gradual growth of the Z-form stretch until it occupies the whole insert. By fitting the theoretical transition curves to experimental ones the fundamental thermodynamic parameters of the B-Z transition have been determined: the B-Z junction energy Fj = 4-5kcal.mol-1 and the free energy change delta FBZ = 0.5-7.0 kcal.mol-1 under standard salt conditions. Calculations show that the B-Z transition in short purinepyrimidine inserts may be seriously affected by cruciform formation in the carrier DNA.     

Conformational properties of a peptide model for unfolded alpha-helices

Models of protein folding often hypothesize that the first step is local secondary structure formation. The assumption is that unfolded polypeptide chains possess an intrinsic propensity to form these local secondary structures. On the basis of this idea, it is tempting to model the local conformational properties of unfolded proteins using well-established residue secondary structure propensities, in particular, alpha-helix forming propensities. We have used spectroscopic methods to investigate the conformational behavior of a host-guest series of peptides designed to model unfolded alpha-helices. A suitable peptide model for unfolded alpha-helices was determined from studies of the length dependence of the conformational properties of alanine-based peptides. The chosen host peptide possessed a small, detectable, alpha-helix content. Substituting various representative guest residues into the central position of the host peptide at times changed the conformational behavior dramatically, and often in ways that could not be predicted from known alpha-helix forming propensities. The data presented can be used to rationalize some of these propensities. However, it is clear that secondary structure propensities cannot be used to predict the local conformational properties of unfolded proteins.     

Energetics and self-assembly of amphipathic peptide pores in lipid membranes

We present a theoretical study of the energetics, equilibrium size, and size distribution of membrane pores composed of electrically charged amphipathic peptides. The peptides are modeled as cylinders (mimicking alpha-helices) carrying different amounts of charge, with the charge being uniformly distributed over a hydrophilic face, defined by the angle subtended by polar amino acid residues. The free energy of a pore of a given radius, R, and a given number of peptides, s, is expressed as a sum of the peptides' electrostatic charging energy (calculated using Poisson-Boltzmann theory), and the lipid-perturbation energy associated with the formation of a membrane rim (which we model as being semitoroidal) in the gap between neighboring peptides. A simple phenomenological model is used to calculate the membrane perturbation energy. The balance between the opposing forces (namely, the radial free energy derivatives) associated with the electrostatic free energy that favors large R, and the membrane perturbation term that favors small R, dictates the equilibrium properties of the pore. Systematic calculations are reported for circular pores composed of various numbers of peptides, carrying different amounts of charge (1-6 elementary, positive charges) and characterized by different polar angles. We find that the optimal R's, for all (except, possibly, very weakly) charged peptides conform to the "toroidal" pore model, whereby a membrane rim larger than approximately 1 nm intervenes between neighboring peptides. Only weakly charged peptides are likely to form "barrel-stave" pores where the peptides essentially touch one another. Treating pore formation as a two-dimensional self-assembly phenomenon, a simple statistical thermodynamic model is formulated and used to calculate pore size distributions. We find that the average pore size and size polydispersity increase with peptide charge and with the amphipathic polar angle. We also argue that the transition of peptides from the adsorbed to the inserted (membrane pore) state is cooperative and thus occurs rather abruptly upon a change in ambient conditions.     

Elucidating quantitative stability/flexibility relationships within thioredoxin and its fragments using a distance constraint model

Numerous quantitative stability/flexibility relationships, within Escherichia coli thioredoxin (Trx) and its fragments are determined using a minimal distance constraint model (DCM). A one-dimensional free energy landscape as a function of global flexibility reveals Trx to fold in a low-barrier two-state process, with a voluminous transition state. Near the folding transition temperature, the native free energy basin is markedly skewed to allow partial unfolded forms. Under native conditions the skewed shape is lost, and the protein forms a compact structure with some flexibility. Predictions on ten Trx fragments are generally consistent with experimental observations that they are disordered, and that complementary fragments reconstitute. A hierarchical unfolding pathway is uncovered using an exhaustive computational procedure of breaking interfacial cross-linking hydrogen bonds that span over a series of fragment dissociations. The unfolding pathway leads to a stable core structure (residues 22-90), predicted to act as a kinetic trap. Direct connection between degree of rigidity within molecular structure and non-additivity of free energy is demonstrated using a thermodynamic cycle involving fragments and their hierarchical unfolding pathway. Additionally, the model provides insight about molecular cooperativity within Trx in its native state, and about intermediate states populating the folding/unfolding pathways. Native state cooperativity correlation plots highlight several flexibly correlated regions, giving insight into the catalytic mechanism that facilitates access to the active site disulfide bond. Residual native cooperativity correlations are present in the core substructure, suggesting that Trx can function when it is partly unfolded. This natively disordered kinetic trap, interpreted as a molten globule, has a wide temperature range of metastability, and it is identified as the "slow intermediate state" observed in kinetic experiments. These computational results are found to be in overall agreement with a large array of experimental data.     

Molecular origin of polyglutamine aggregation in neurodegenerative diseases

Expansion of polyglutamine (polyQ) tracts in proteins results in protein aggregation and is associated with cell death in at least nine neurodegenerative diseases. Disease age of onset is correlated with the polyQ insert length above a critical value of 35-40 glutamines. The aggregation kinetics of isolated polyQ peptides in vitro also shows a similar critical-length dependence. While recent experimental work has provided considerable insights into polyQ aggregation, the molecular mechanism of aggregation is not well understood. Here, using computer simulations of isolated polyQ peptides, we show that a mechanism of aggregation is the conformational transition in a single polyQ peptide chain from random coil to a parallel beta-helix. This transition occurs selectively in peptides longer than 37 glutamines. In the beta-helices observed in simulations, all residues adopt beta-strand backbone dihedral angles, and the polypeptide chain coils around a central helical axis with 18.5 +/- 2 residues per turn. We also find that mutant polyQ peptides with proline-glycine inserts show formation of antiparallel beta-hairpins in their ground state, in agreement with experiments. The lower stability of mutant beta-helices explains their lower aggregation rates compared to wild type. Our results provide a molecular mechanism for polyQ-mediated aggregation.     

Intrinsically disordered protein from a pathogenic mesophile Mycobacterium tuberculosis adopts structured conformation at high temperature

Compared to eukaryotes, the occurrence of "intrinsically disordered" or "natively unfolded" proteins in prokaryotes has not been explored extensively. Here, we report the occurrence of an intrinsically disordered protein from the mesophilic human pathogen Mycobacterium tuberculosis. The Histidine-tagged recombinant Rv3221c biotin-binding protein is intrinsically disordered at ambient and physiological growth temperatures as revealed by circular dichroism and Fourier transform infrared (FTIR) spectroscopic studies. However, an increase in temperature induces a transition from disordered to structured state with a folding temperature of approximately 53 degrees C. Addition of a structure inducing solvent trifluoroethanol (TFE) causes the protein to fold at lower temperatures suggesting that TFE fosters hydrophobic interactions, which drives protein folding. Differential Scanning Calorimetry studies revealed that folding is endothermic and the transition from a disordered to structured state is continuous (higher-order), implying existence of intermediates during folding process. Secondary structure analysis revealed that the protein has propensity to form beta-sheets. This is in conformity with FTIR spectrum that showed an absorption peak at wave number of 1636 cm(-1), indicative of disordered beta-sheet conformation in the native state. These data suggest that although Rv3221c may be disordered under ambient or optimal growth temperature conditions, it has the potential to fold into ordered structure at high temperature driven by increased hydrophobic interactions. In contrast to the generally known behavior of other intrinsically disordered proteins folding at high temperature, Rv3221c does not appear to oligomerize or aggregate as revealed through numerous experiments including Congo red binding, Thioflavin T-binding, turbidity measurements, and examining molar ellipticity as a function of protein concentration. The amino acid composition of Rv3221c reveals that it has 24% charged and 54.9% hydrophobic amino acid residues. In this respect, this protein, although belonging to the class of intrinsically disordered proteins, has distinct features. The intrinsically disordered state and the biotin-binding feature of this protein suggest that it may participate in many biochemical processes requiring biotin as a cofactor and adopt suitable conformations upon binding other folded targets.     

Conversion of two-state to multi-state folding kinetics on fusion of two protein foldons

Chymotrypsin inhibitor 2 (CI2) is the archetypal single-foldon protein that folds in simple two-state kinetics without the accumulation of a folding intermediate. To model the effects of fusion of single foldons to give a multi-foldon protein, we engineered a "double-CI2" protein, in which another CI2 polypeptide was inserted into the loop region of the parent CI2. CD and HSQC spectra demonstrated that while the double-CI2 protein adopted two kinds of native conformations, CI2-like structure was almost preserved in both the domains of double-CI2. In the folding kinetic studies, double-CI2 exhibited a remarkable rollover of the observed folding rates at low denaturant concentrations, indicating that double-CI2 accumulated a kinetic folding intermediate. The different folding mechanisms between WT-CI2 and double-CI2 support the present view that protein size or number of domains is an important determinant for formation of folding intermediates.     

Effects of chain length on the aggregation of model polyglutamine peptides: molecular dynamics simulations

Aggregation in the brain of polyglutamine-containing proteins is either a cause or an associated symptom of nine hereditary neurodegenerative disorders including Huntington's disease. The molecular level mechanisms by which these proteins aggregate are still unclear. In an effort to shed light on this important phenomenon, we are investigating the aggregation of model polyglutamine peptides using molecular-level computer simulation with a simplified model of polyglutamine that we have developed. This model accounts for the most important types of intra- and inter-molecular interactions-hydrogen bonding and hydrophobic interactions-while allowing the folding process to be simulated in a reasonable time frame. The model is used to examine the folding of isolated polyglutamine peptides 16, 32, and 48 residues long and the folding and aggregation of systems of 24 model polyglutamine peptides 16, 24, 32, 36, 40, and 48 residues long. Although the isolated polyglutamine peptides did form some alpha and beta backbone-backbone hydrogen bonds they did not have as many of these bonds as they would have if they had folded into a complete alpha helix or beta sheet. In one of the simulations on the isolated polyglutamine peptide 48 residues long, we observed a structure that resembles a beta helix. In the multi-chain simulations we observed amorphous aggregates at low temperatures, ordered aggregates with significant beta sheet character at intermediate temperatures, and random coils at high temperatures. We have found that the temperature at which the model peptides undergo the transition from amorphous aggregates to ordered aggregates and the temperature at which the model peptides undergo the transition from ordered aggregates to random coils increase with increasing chain length. Our finding that the stability of the ordered aggregates increases as the peptide chain length increases may help to explain the experimentally observed relation between polyglutamine tract length and aggregation in vitro and disease progression in vivo. We have also observed in our simulations that the optimal temperature for the formation of beta sheets increases with chain length up to 36 glutamine residues but not beyond. Equivalently, at fixed temperature we find a transition from a region dominated by random coils at chain lengths less than 36 to a region dominated by relatively ordered beta sheet structures at chain lengths greater than 36. Our finding of this critical chain length of 36 glutamine residues is interesting because a critical chain length of 37 glutamine residues has been observed experimentally.     

Atomistic simulations of the effects of polyglutamine chain length and solvent quality on conformational equilibria and spontaneous homodimerization

Aggregation of expanded polyglutamine tracts is associated with nine different neurodegenerative diseases, including Huntington's disease. Experiments and computer simulations have demonstrated that monomeric forms of polyglutamine molecules sample heterogeneous sets of collapsed structures in water. The current work focuses on a mechanistic characterization of polyglutamine homodimerization as a function of chain length and temperature. These studies were carried out using molecular simulations based on a recently developed continuum solvation model that was designed for studying conformational and binding equilibria of intrinsically disordered molecules such as polyglutamine systems. The main results are as follows: Polyglutamine molecules form disordered, collapsed globules in aqueous solution. These molecules spontaneously associate at conditions approaching those of typical in vitro experiments for chains of length N ≥ 15. The spontaneity of these homotypic associations increases with increasing chain length and decreases with increasing temperature. Similar and generic driving forces govern both collapse and spontaneous homodimerization of polyglutamine in aqueous milieus. Collapse and dimerization maximize self-interactions and reduce the interface between polyglutamine molecules and the surrounding solvent. Other than these generic considerations, there do not appear to be any specific structural requirements for either chain collapse or chain dimerization; that is, both collapse and dimerization are nonspecific in that disordered globules form disordered dimers. In fact, it is shown that the driving force for intermolecular associations is governed by spontaneous conformational fluctuations within monomeric polyglutamine. These results suggest that polyglutamine aggregation is unlikely to follow a homogeneous nucleation mechanism with the monomer as the critical nucleus. Instead, the results support the formation of disordered, non-β-sheet-like soluble molten oligomers as early intermediates—a proposal that is congruent with recent experimental data.     

Does alpha-helix folding necessarily provide an energy source for the protein-lipid binding?

Lipid-induced alpha-helix folding, which occurs in many lipid surface-binding proteins and peptides such as apolipoproteins and synucleins, has been proposed to provide an energy source for protein-lipid interactions. We propose that in a system comprised of a phospholipid surface and a small polypeptide that is unfolded in solution and binds reversibly to lipid surface, helical folding involves expenditure of free energy as compared to a similar polypeptide that is alpha-helical in solution. This is a consequence of the entropic cost of helix folding that is illustrated in a simple thermodynamic model and exemplifies the general "key-into-lock" paradigm of protein-ligand binding. Even though this simple model does not explicitly address the protein-induced lipid re-arrangement and may not directly apply to large proteins that undergo significant tertiary structural changes upon lipid binding, it suggests that the notion of helix folding as an energy source for lipid binding should be treated with caution.     

Side-chain interactions determine amyloid formation by model polyglutamine peptides in molecular dynamics simulations

The pathological manifestation of nine hereditary neurodegenerative diseases is the presence within the brain of aggregates of disease-specific proteins that contain polyglutamine tracts longer than a critical length. To improve our understanding of the processes by which polyglutamine-containing proteins misfold and aggregate, we have conducted molecular dynamics simulations of the aggregation of model polyglutamine peptides. This work was accomplished by extending the PRIME model to polyglutamine. PRIME is an off-lattice, unbiased, intermediate-resolution protein model based on an amino acid representation of between three and seven united atoms, depending on the residue being modeled. The effects of hydrophobicity on the system are studied by varying the strength of the hydrophobic interaction from 12.5% to 5% of the hydrogen-bonding interaction strength. In our simulations, we observe the spontaneous formation of aggregates and annular structures that are made up of beta-sheets starting from random configurations of random coils. This result was interesting because tubular protofibrils were recently found in experiments on polyglutamine aggregation and because of Perutz's prediction that polyglutamine would form water-filled nanotubes.