Refolding simulations of an isolated fragment of barnase into a native-like β hairpin: Evidence for compactness and hydrogen bonding as concurrent stabilizing factors

Experimental evidence and theoretical models both suggest that protein folding is initiated within specific fragments intermittently adopting conformations close to that found in the protein native structure. These folding initiation sites encompassing short portions of the protein are ideally suited for study in isolation by computational methods aimed at peering into the very early events of folding. We have used Molecular Dynamics (MD) technique to investigate the behavior of an isolated protein fragment formed by residues 85 to 102 of barnase that folds into a β hairpin in the protein native structure. Three independent MD simulations of 1.3 to 1.8 ns starting from unfolded conformations of the peptide portrayed with an all-atom model in water were carried out at gradually decreasing temperature. A detailed analysis of the conformational preferences adopted by this peptide in the course of the simulations is presented. Two of the unfolded peptide conformations fold into a hairpin characterized by native and a larger bulk of nonnative interactions. Both refolding simulations substantiate the close relationship between interstrand compactness and hydrogen bonding network involving backbone atoms. Persistent compactness witnessed by side-chain interactions always occurs concomitantly with the formation of backbone hydrogen bonds. No highly populated conformations generated in a third simulation starting from the remotest unfolded conformer relative to the native structure are observed. However, nonnative long-range and medium-range contacts with the aromatic moiety of Trp94 are spotted, which are in fair agreement with a former nuclear magnetic resonance study of a denaturing solution of an isolated barnase fragment encompassing the β hairpin. All this lends reason to believe that the 85–102 barnase fragment is a strong initiation site for folding. Proteins 29:212–227, 1997. © 1997 Wiley-Liss, Inc.     

All-atom computer simulations of amyloid fibrils disaggregation

Amyloidlike fibrils are found in many fatal diseases, including Alzheimer's disease, type II diabetes mellitus, transmissible spongiform encephalopathies, and prion diseases. These diseases are linked to proteins that have partially unfolded, misfolded, and aggregated into amyloidlike fibrils. The kinetics of amyloidlike fibrils aggregation is still hotly debated and remains an important open question. We have utilized the GNNQQNY crystal structure and high-temperature molecular dynamics simulation in explicit solvent to study the disaggregation mechanism of the GNNQQNY fibrils and to infer its likely aggregation pathways. A hexamer model and a 12-mer model both with two parallel β-sheets separated by a dry side-chain interface were adopted in our computational analysis. A cumulative time of 1 μs was simulated for the hexamer model at five different temperatures (298 K, 348 K, 398 K, 448 K, and 498 K), and a cumulative time of 2.1 μs was simulated for the 12-mer model at four temperatures (298 K, 398 K, 448 K, and 498 K). Our disaggregation landscape and kinetics analyses indicate that tetramers probably act as the transition state in both the hexamer and the 12-mer simulations. In addition, the 12-mer simulations show that the initial aggregation nucleus is with eight peptides. Furthermore, the landscape is rather flat from 8-mers to 12-mers, indicating the absence of major barriers once the initial aggregation nucleus forms. Thus, the likely aggregation pathway is from monomers to the initial nucleus of 8-mers with tetramers as the transition state. Transition state structure analysis shows that the two dominant transition state conformations are tetramers in the 3-1 and 2-2 arrangements. The predominant nucleus conformations are in peptide arrangements maximizing dry side-chain contacts. Landscape and kinetics analyses also indicate that the parallel β-sheets form earlier than the dry side-chain contacts during aggregation. These results provide further insights in understanding the early fibrils aggregation.     

Free energy landscapes of two model peptides: α-helical and β-hairpin peptides explored with Brownian dynamics simulation

We applied an atomistic Brownian dynamics (BD) simulation with multiple time step method for the folding simulation of a 13-mer α-helical peptide and a 12-mer β-hairpin peptide, giving successful folding simulations. In this model, the driving energy contribution towards folding came from both electrostatic and van der Waals interactions for the α-helical peptide and from van der Waals interactions for the β-hairpin peptide. Although, many non-native structures having the same or lower energy than that of native structure were observed, the folded states formed the most populated cluster when the structures obtained by the BD simulations were subjected to the cluster analysis based on distance-based root mean square deviation of side-chains between different structures. This result indicates that we can predict the native structures from conformations sampled by BD simulation.     

Dimensionality of Carbon Nanomaterials Determines the Binding and Dynamics of Amyloidogenic Peptides: Multiscale Theoretical Simulations

Experimental studies have demonstrated that nanoparticles can affect the rate of protein self-assembly, possibly interfering with the development of protein misfolding diseases such as Alzheimer''s, Parkinson''s and prion disease caused by aggregation and fibril formation of amyloid-prone proteins. We employ classical molecular dynamics simulations and large-scale density functional theory calculations to investigate the effects of nanomaterials on the structure, dynamics and binding of an amyloidogenic peptide apoC-II(60-70). We show that the binding affinity of this peptide to carbonaceous nanomaterials such as C60, nanotubes and graphene decreases with increasing nanoparticle curvature. Strong binding is facilitated by the large contact area available for π-stacking between the aromatic residues of the peptide and the extended surfaces of graphene and the nanotube. The highly curved fullerene surface exhibits reduced efficiency for π-stacking but promotes increased peptide dynamics. We postulate that the increase in conformational dynamics of the amyloid peptide can be unfavorable for the formation of fibril competent structures. In contrast, extended fibril forming peptide conformations are promoted by the nanotube and graphene surfaces which can provide a template for fibril-growth.     

Effects of Macromolecular Crowding on Protein Conformational Changes

Many protein functions can be directly linked to conformational changes. Inside cells, the equilibria and transition rates between different conformations may be affected by macromolecular crowding. We have recently developed a new approach for modeling crowding effects, which enables an atomistic representation of “test” proteins. Here this approach is applied to study how crowding affects the equilibria and transition rates between open and closed conformations of seven proteins: yeast protein disulfide isomerase (yPDI), adenylate kinase (AdK), orotidine phosphate decarboxylase (ODCase), Trp repressor (TrpR), hemoglobin, DNA β-glucosyltransferase, and Ap₄A hydrolase. For each protein, molecular dynamics simulations of the open and closed states are separately run. Representative open and closed conformations are then used to calculate the crowding-induced changes in chemical potential for the two states. The difference in chemical-potential change between the two states finally predicts the effects of crowding on the population ratio of the two states. Crowding is found to reduce the open population to various extents. In the presence of crowders with a 15 Å radius and occupying 35% of volume, the open-to-closed population ratios of yPDI, AdK, ODCase and TrpR are reduced by 79%, 78%, 62% and 55%, respectively. The reductions for the remaining three proteins are 20–44%. As expected, the four proteins experiencing the stronger crowding effects are those with larger conformational changes between open and closed states (e.g., as measured by the change in radius of gyration). Larger proteins also tend to experience stronger crowding effects than smaller ones [e.g., comparing yPDI (480 residues) and TrpR (98 residues)]. The potentials of mean force along the open-closed reaction coordinate of apo and ligand-bound ODCase are altered by crowding, suggesting that transition rates are also affected. These quantitative results and qualitative trends will serve as valuable guides for expected crowding effects on protein conformation changes inside cells.     

β-Hairpin-Mediated Formation of Structurally Distinct Multimers of Neurotoxic Prion Peptides

Protein misfolding disorders are associated with conformational changes in specific proteins, leading to the formation of potentially neurotoxic amyloid fibrils. During pathogenesis of prion disease, the prion protein misfolds into β-sheet rich, protease-resistant isoforms. A key, hydrophobic domain within the prion protein, comprising residues 109–122, recapitulates many properties of the full protein, such as helix-to-sheet structural transition, formation of fibrils and cytotoxicity of the misfolded isoform. Using all-atom, molecular simulations, it is demonstrated that the monomeric 109–122 peptide has a preference for α-helical conformations, but that this peptide can also form β-hairpin structures resulting from turns around specific glycine residues of the peptide. Altering a single amino acid within the 109–122 peptide (A₁₁₇V, associated with familial prion disease) increases the prevalence of β-hairpin formation and these observations are replicated in a longer peptide, comprising residues 106–126. Multi-molecule simulations of aggregation yield different assemblies of peptide molecules composed of conformationally-distinct monomer units. Small molecular assemblies, consistent with oligomers, comprise peptide monomers in a β-hairpin-like conformation and in many simulations appear to exist only transiently. Conversely, larger assemblies are comprised of extended peptides in predominately antiparallel β-sheets and are stable relative to the length of the simulations. These larger assemblies are consistent with amyloid fibrils, show cross-β structure and can form through elongation of monomer units within pre-existing oligomers. In some simulations, assemblies containing both β-hairpin and linear peptides are evident. Thus, in this work oligomers are on pathway to fibril formation and a preference for β-hairpin structure should enhance oligomer formation whilst inhibiting maturation into fibrils. These simulations provide an important new atomic-level model for the formation of oligomers and fibrils of the prion protein and suggest that stabilization of β-hairpin structure may enhance cellular toxicity by altering the balance between oligomeric and fibrillar protein assemblies.     

All-Atom Monte Carlo Approach to Protein-Peptide Binding

We develop a procedure for exploring the free energy landscape of protein-peptide binding at atomic detail and apply it to PDZ domain-peptide interactions. The procedure involves soft constraints on receptor proteins providing limited chain flexibility, including backbone motions. Peptide chains are left fully flexible and kept in spatial proximity of the protein through periodic boundary conditions. By extensive Monte Carlo simulations, full representative conformational ensembles at temperatures where bound and unbound states coexist are obtained. To make this approach computationally feasible, we develop an effective all-atom energy function centering on hydrophobicity, hydrogen bonding, and electrostatic interactions. Our initial focus is a set of 11 PDZ domain-peptide pairs with experimentally determined complex structures. Minimum-energy conformations are found to be highly similar to the respective native structures in eight of the cases (all-atom peptide RMSDs < 6 Å). Having achieved that, we turn to a more complete characterization of the bound peptide state through a clustering scheme applied on the full ensembles of peptide structures. We find a significant diversity among bound peptide conformations for several PDZ domains, in particular involving the N terminal side of the peptide chains. Our computational model is then tested further on a set of nine PDZ domain-peptide pairs where the peptides are not originally present in the experimentally determined structures. We find a similar success rate in terms of the nativeness of minimum-energy conformations. Finally, we investigate the ability of our approach to capture variations in binding affinities for different peptide sequences. This is done in particular for a set of related sequences binding to the third PDZ domain of PSD-95 with encouraging results.     

Ab initio construction of polypeptide fragments: Accuracy of loop decoy discrimination by an all-atom statistical potential and the AMBER force field with the Generalized Born solvation model

The accuracy of model selection from decoy ensembles of protein loop conformations was explored by comparing the performance of the Samudrala-Moult all-atom statistical potential (RAPDF) and the AMBER molecular mechanics force field, including the Generalized Born/surface area solvation model. Large ensembles of consistent loop conformations, represented at atomic detail with idealized geometry, were generated for a large test set of protein loops of 2 to 12 residues long by a novel ab initio method called RAPPER that relies on fine-grained residue-specific phi/psi propensity tables for conformational sampling. Ranking the conformers on the basis of RAPDF scores resulted in selected conformers that had an average global, non-superimposed RMSD for all heavy mainchain atoms ranging from 1.2 A for 4-mers to 2.9 A for 8-mers to 6.2 A for 12-mers. After filtering on the basis of anchor geometry and RAPDF scores, ranking by energy minimization of the AMBER/GBSA potential energy function selected conformers that had global RMSD values of 0.5 A for 4-mers, 2.3 A for 8-mers, and 5.0 A for 12-mers. Minimized fragments had, on average, consistently lower RMSD values (by 0.1 A) than their initial conformations. The importance of the Generalized Born solvation energy term is reflected by the observation that the average RMSD accuracy for all loop lengths was worse when this term is omitted. There are, however, still many cases where the AMBER gas-phase minimization selected conformers of lower RMSD than the AMBER/GBSA minimization. The AMBER/GBSA energy function had better correlation with RMSD to native than the RAPDF. When the ensembles were supplemented with conformations extracted from experimental structures, a dramatic improvement in selection accuracy was observed at longer lengths (average RMSD of 1.3 A for 8-mers) when scoring with the AMBER/GBSA force field. This work provides the basis for a promising hybrid approach of ab initio and knowledge-based methods for loop modeling.     

Distributions in protein conformation space: implications for structure prediction and entropy

By considering how polymer structures are distributed in conformation space, we show that it is possible to quantify the difficulty of structural prediction and to provide a measure of progress for prediction calculations. The critical issue is the probability that a conformation is found within a specified distance of another conformer. We address this question by constructing a cumulative distribution function (CDF) for the average probability from observations about its limiting behavior at small displacements and numerical simulations of polyalanine chains. We can use the CDF to estimate the likelihood that a structure prediction is better than random chance. For example, the chance of randomly predicting the native backbone structure of a 150-amino-acid protein to low resolution, say within 6 A, is 10(-14). A high-resolution structural prediction, say to 2 A, is immensely more difficult (10(-57)). With additional assumptions, the CDF yields the conformational entropy of protein folding from native-state coordinate variance. Or, using values of the conformational entropy change on folding, we can estimate the native state's conformational span. For example, for a 150-mer protein, equilibrium alpha-carbon displacements in the native ensemble would be 0.3-0.5 A based on T Delta S of 1.42 kcal/(mol residue).     

Design and synthesis of novel nonpolar host peptides for the determination of the 310- and α-helix compatibilities of α-amino acid buildig blocks: An assessment of α,α-disubstituted glycines

The present work describes three novel nonpolar host peptide sequences that provide a ready assessment of the 3₁₀- and α-helix compatibilities of natural and unnatural amino acids at different positions of small- to medium-size peptides. The unpolar peptides containing Ala, Aib, and a C-terminal p-iodoanilide group were designed in such a way that the peptides could be rapidly assembled in a modular fashion, were highly soluble in solvent mixtures of triflouroethanol and H₂O for CD- and two-dimensional (2D) nmr spectroscopic analyses, and showed excellent crystallinity suited for x-ray structure analysis. To validate our approach we synthesized 9-mer peptides 79a–96 (Table IV), 12-mer peptides 99–110c (Table V), and 10-mer peptides 120a–125d and 129–133 (Table VI and Scheme 8) incorporating a series of optically pure cyclic and open-chain (R)- and (S)-α,α-disubstituted glycines 1–10 (Figure 2). These amino acids are known to significantly modulate the conformations of small peptides. Based on x-ray structures of 9-mers 79a, 80, and 87 (Figures 4–7), 10-mers 124c, 131, and 132 (Figures 9–12), and 12-mer peptide 102b (Figure 13), CD spectra of all peptides recorded in acidic, neutral, and basic media and detailed 2D-nmr analyses of 9-mer peptide 86 and 12-mer 102b, several interesting conformational observations were made. Especially interesting results were obtained using the convex constraint CD analysis proposed by Fasman on 9-mer peptides 79a–d, 80, 81, 86, and 87, which allowed us to determine the relative content of 3₁₀- and α-helical conformations. These results were fully supported by the corresponding x-ray and 2D-nmr analyses. As a striking example we found that the (S)- and (R)-β-tetralin derived amino acids (R)- and (S)-1 show excellent α-helix stabilisation, more pronounced than Aib and Ala. These novel reference peptide sequences should help establish a scale for natural and unnatural amino acids concerning their intrinsic 3₁₀- and α-helix compatibilities at different positions of medium-sized peptides and thus improve our understanding in the folding processes of peptides. © 1997 John Wiley & Sons, Inc. Biopoly 42: 575–626, 1997     

Packing and recognition of protein structural elements: A new approach applied to the 4-helix bundle of myohemerythrin

We present a novel search strategy for determining the optimal packing of protein secondary structure elements. The approach is based on conformational energy optimization using a predetermined set of side chain rotamers and appropriate methods for sampling the conformational space of peptide fragments having fixed backbone geometries. An application to the 4-helix bundle of myohemerythrin is presented. It is shown that the conformations of the amino acid side chains are largely determined at the level of helix pairs and that superposition of these results can be used to construct the full bundle. The final solution obtained, taking into account restrictions due to the lateral amphiphilicity of the helices, differs from the native structure by only a 20° rotation of a single helix. © 1993 Wiley-Liss, Inc.     

Mapping of the interleukin-10/interleukin-10 receptor combining site.

The discontinuous interleukin-10(IL-10)/interleukin-10 receptor (IL-10R) combining site was mapped using sets of overlapping peptides derived from both binding partners bound to continuous cellulose membranes. Low affinity binding of single regions of the discontinuous contact sites on IL-10 and IL-10R could be identified due to (1) high peptide density on the membrane support, (2) incubation with high protein concentrations, (3) indirect immunodetection of the ligates after electrotransfer onto polyvinylene difluoride membranes, and (4) use of highly overlapping peptide scans of different length (6-mers and 15-mers). The single binding regions identified for each protein species are separated in the protein sequences, but form continuous areas on the surface of IL-10 (X-ray structure) and IL-10R (computer model). Furthermore, four epitopes of neutralizing anti-IL-10 and anti-IL-10R antibodies were mapped and overlap with these binding regions. Soluble peptides (15- to 19-mers) each spanning one of the three identified IL-10-derived receptor binding regions displayed no significant affinity to IL-10R as expected, whereas a peptide (35-mer) comprising two of these regions had considerably higher binding activity. The data are consistent with a previously published computer model of the IL-10/IL-10R complex. This approach should be generally applicable for the mapping of non-linear protein-protein contact sites.     

Structural and dynamic determinants of protein-peptide recognition

Protein-peptide interactions play important roles in many cellular processes, including signal transduction, trafficking, and immune recognition. Protein conformational changes upon binding, an ill-defined peptide binding surface, and the large number of peptide degrees of freedom make the prediction of protein-peptide interactions particularly challenging. To address these challenges, we perform rapid molecular dynamics simulations in order to examine the energetic and dynamic aspects of protein-peptide binding. We find that, in most cases, we recapitulate the native binding sites and native-like poses of protein-peptide complexes. Inclusion of electrostatic interactions in simulations significantly improves the prediction accuracy. Our results also highlight the importance of protein conformational flexibility, especially side-chain movement, which allows the peptide to optimize its conformation. Our findings not only demonstrate the importance of sufficient sampling of the protein and peptide conformations, but also reveal the possible effects of electrostatics and conformational flexibility on peptide recognition.     

Interactive NMR and computer simulation studies of lanthionine-ring structures

We report progress in elucidating the structure of nisin, a naturally occurring peptide antibiotic. Nisin contains five rings constrained by lanthionine or methyllanthionine bridges, as well as alpha, beta-unsaturated amino acids. We have determined conformations for two model compounds of ring A and a derivative of ring B through interactive nmr and computer simulation studies. High-resolution nmr techniques provides structural information, which was further refined through molecular dynamics simulations. These methods are being applied to the remaining constrained fragments of the molecule. This conformational information will be employed in an aufbau approach to determining the structure of the entire molecule.     

Conformational adjustments of SNase R and its N-terminal fragments probed by monoclonal antibodies

Two monoclonal antibodies specific for staphylococcal nuclease R (SNase R) (McAb2C9 and McAb1B8) were prepared and used to probe protein folding during peptide elongation, by measuring antibody binding to seven N-terminal fragments (SNR141, SNR135, SNR121, SNR110, SNR102, SNR79 and SNR52) of SNase R. Comparative studies of the conformations of the N-terminal fragments have shown that all seven fragments of SNase R have a certain amount of residual structure, indicating that folding may occur during elongation of the nascent peptide chain. We show that the binding abilities of the intact enzyme and its seven fragments to the monoclonal antibodies are not simply proportional to the length of the peptide chain, suggesting that there may be continuous conformational adjustment in the nascent peptide chain as new C-terminal amino acids are added. A folding intermediate close in structure to the native state but with structural features in common with SNR121 is highly populated in 0.6 M GuHCl, and is also formed transiently during folding.     

The RecB protein of Escherichia coli translocates along single-stranded DNA in the 3′ to 5′ direction: a proposed ratchet mechanism

To investigate the role that the individual subunits play in the ATP-dependent helicase activity of the RecBCD protein we have investigated the ability of the RecB, RecC and RecD proteins to displace various 20-mer oligonucleotides annealed to either end or to the centre of an oligonucleotide 60 bases long. The results show that the only subunit which can displace the 20-mers in the absence of the other subunits is the RecB protein. Moreover, the 20-mer is displaced only if it is annealed to the 60-mer at the 5′ end or the middle, suggesting that the RecB protein translocates along the 60-mer in the 3′ to 5′ direction, displacing annealed 20-mers as it proceeds. We have shown that reconstituted RecBC and RecBCD complexes displace the 20-mers but, unlike RecB, they do not require a 3′-ended single-stranded region for helicase action, but can displace the 20-mers from either end of the 60-mer. The level of helicase activity of the RecBC complex is considerably greater than that of RecB alone, and the activity of the RecBCD complex appears to be greater still. This hierarchy of activity is also shown by DNA binding studies, but is not reflected in the ATPase activities of the enzymes. We have also shown that the ability of trypsin to cleave various sites on the RecB molecule is modified by the presence of ATP or ATP-γ-S, suggesting that conformational changes may be induced in RecB upon ATP binding. We discuss a model for the ATP-driven, unidirectional motion of the RecB translocase along single-stranded DNA. In this model, the RecB molecule binds to single-stranded DNA and then translocates along it, one base at a time, in the 3′ to 5′ direction, by a `ratchet' mechanism in which repeated stretching and contraction of the protein is coupled to ATP hydrolysis. The RecC protein in the RecBC complex is proposed to act as a `sliding clamp' which increases processivity by preventing dissociation.     

Enhanced Regulatory Sequence Prediction Using Gapped k-mer Features

Oligomers of length k, or k-mers, are convenient and widely used features for modeling the properties and functions of DNA and protein sequences. However, k-mers suffer from the inherent limitation that if the parameter k is increased to resolve longer features, the probability of observing any specific k-mer becomes very small, and k-mer counts approach a binary variable, with most k-mers absent and a few present once. Thus, any statistical learning approach using k-mers as features becomes susceptible to noisy training set k-mer frequencies once k becomes large. To address this problem, we introduce alternative feature sets using gapped k-mers, a new classifier, gkm-SVM, and a general method for robust estimation of k-mer frequencies. To make the method applicable to large-scale genome wide applications, we develop an efficient tree data structure for computing the kernel matrix. We show that compared to our original kmer-SVM and alternative approaches, our gkm-SVM predicts functional genomic regulatory elements and tissue specific enhancers with significantly improved accuracy, increasing the precision by up to a factor of two. We then show that gkm-SVM consistently outperforms kmer-SVM on human ENCODE ChIP-seq datasets, and further demonstrate the general utility of our method using a Naïve-Bayes classifier. Although developed for regulatory sequence analysis, these methods can be applied to any sequence classification problem.     

A Nine-Residue Synthetic Propeptide Enhances Secretion Efficiency of Heterologous Proteins in Lactococcus lactis

Lactococcus lactis, a gram-positive organism widely used in the food industry, is a potential candidate for the secretion of biologically useful proteins. We examined the secretion efficiency and capacity of L. lactis by using the Staphylococcus aureus nuclease (Nuc) as a heterologous model protein. When expressed in L. lactis from an efficient lactococcal promoter and its native signal peptide, only ~60% of total Nuc was present in a secreted form at ~5 mg per liter. The remaining 40% was found in a cell-associated precursor form. The secretion efficiency was reduced further to ~30% by the deletion of 17 residues of the Nuc native propeptide (resulting in NucT). We identified a modification which improved secretion efficiency of both native Nuc and NucT. A 9-residue synthetic propeptide, LEISSTCDA, which adds two negative charges at the +2 and +8 positions, was fused immediately after the signal peptide cleavage site. In the case of Nuc, secretion efficiency was increased to ~80% by LEISSTCDA insertion without altering the signal peptide cleavage site, and the yield was increased two- to fourfold (up to ~20 mg per liter). The improvement of NucT secretion efficiency was even more marked and rose from 30 to 90%. Similarly, the secretion efficiency of a third protein, the α-amylase of Bacillus stearothermophilus, was also improved by LEISSTCDA. These data indicate that the LEISSTCDA synthetic propeptide improves secretion of different heterologous proteins in L. lactis.