A Multiscale Approach to Characterize the Early Aggregation Steps of the Amyloid-Forming Peptide GNNQQNY from the Yeast Prion Sup-35
Identifieur interne : 000F92 ( Pmc/Corpus ); précédent : 000F91; suivant : 000F93A Multiscale Approach to Characterize the Early Aggregation Steps of the Amyloid-Forming Peptide GNNQQNY from the Yeast Prion Sup-35
Auteurs : Jessica Nasica-Labouze ; Massimiliano Meli ; Philippe Derreumaux ; Giorgio Colombo ; Normand MousseauSource :
- PLoS Computational Biology [ 1553-734X ] ; 2011.
Abstract
The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.
Url:
DOI: 10.1371/journal.pcbi.1002051
PubMed: 21625573
PubMed Central: 3098217
Links to Exploration step
PMC:3098217Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">A Multiscale Approach to Characterize the Early Aggregation Steps of the Amyloid-Forming Peptide GNNQQNY from the Yeast Prion Sup-35</title><author><name sortKey="Nasica Labouze, Jessica" sort="Nasica Labouze, Jessica" uniqKey="Nasica Labouze J" first="Jessica" last="Nasica-Labouze">Jessica Nasica-Labouze</name><affiliation><nlm:aff id="aff1"><addr-line>Département de Physique and GEPROM, Université de Montréal, Montréal, Québec, Canada</addr-line></nlm:aff></affiliation></author><author><name sortKey="Meli, Massimiliano" sort="Meli, Massimiliano" uniqKey="Meli M" first="Massimiliano" last="Meli">Massimiliano Meli</name><affiliation><nlm:aff id="aff2"><addr-line>Istituto di Chimica del Riconoscimento Molecolare, CNR, Milano, Italy</addr-line></nlm:aff></affiliation></author><author><name sortKey="Derreumaux, Philippe" sort="Derreumaux, Philippe" uniqKey="Derreumaux P" first="Philippe" last="Derreumaux">Philippe Derreumaux</name><affiliation><nlm:aff id="aff3"><addr-line>Laboratoire de Biochimie Théorique, UPR9080 CNRS, Institut de Biologie Physico-Chimique, Université Paris 7, and Institut Universitaire de France, Paris, France</addr-line></nlm:aff></affiliation></author><author><name sortKey="Colombo, Giorgio" sort="Colombo, Giorgio" uniqKey="Colombo G" first="Giorgio" last="Colombo">Giorgio Colombo</name><affiliation><nlm:aff id="aff2"><addr-line>Istituto di Chimica del Riconoscimento Molecolare, CNR, Milano, Italy</addr-line></nlm:aff></affiliation></author><author><name sortKey="Mousseau, Normand" sort="Mousseau, Normand" uniqKey="Mousseau N" first="Normand" last="Mousseau">Normand Mousseau</name><affiliation><nlm:aff id="aff1"><addr-line>Département de Physique and GEPROM, Université de Montréal, Montréal, Québec, Canada</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">21625573</idno><idno type="pmc">3098217</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3098217</idno><idno type="RBID">PMC:3098217</idno><idno type="doi">10.1371/journal.pcbi.1002051</idno><date when="2011">2011</date><idno type="wicri:Area/Pmc/Corpus">000F92</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000F92</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">A Multiscale Approach to Characterize the Early Aggregation Steps of the Amyloid-Forming Peptide GNNQQNY from the Yeast Prion Sup-35</title><author><name sortKey="Nasica Labouze, Jessica" sort="Nasica Labouze, Jessica" uniqKey="Nasica Labouze J" first="Jessica" last="Nasica-Labouze">Jessica Nasica-Labouze</name><affiliation><nlm:aff id="aff1"><addr-line>Département de Physique and GEPROM, Université de Montréal, Montréal, Québec, Canada</addr-line></nlm:aff></affiliation></author><author><name sortKey="Meli, Massimiliano" sort="Meli, Massimiliano" uniqKey="Meli M" first="Massimiliano" last="Meli">Massimiliano Meli</name><affiliation><nlm:aff id="aff2"><addr-line>Istituto di Chimica del Riconoscimento Molecolare, CNR, Milano, Italy</addr-line></nlm:aff></affiliation></author><author><name sortKey="Derreumaux, Philippe" sort="Derreumaux, Philippe" uniqKey="Derreumaux P" first="Philippe" last="Derreumaux">Philippe Derreumaux</name><affiliation><nlm:aff id="aff3"><addr-line>Laboratoire de Biochimie Théorique, UPR9080 CNRS, Institut de Biologie Physico-Chimique, Université Paris 7, and Institut Universitaire de France, Paris, France</addr-line></nlm:aff></affiliation></author><author><name sortKey="Colombo, Giorgio" sort="Colombo, Giorgio" uniqKey="Colombo G" first="Giorgio" last="Colombo">Giorgio Colombo</name><affiliation><nlm:aff id="aff2"><addr-line>Istituto di Chimica del Riconoscimento Molecolare, CNR, Milano, Italy</addr-line></nlm:aff></affiliation></author><author><name sortKey="Mousseau, Normand" sort="Mousseau, Normand" uniqKey="Mousseau N" first="Normand" last="Mousseau">Normand Mousseau</name><affiliation><nlm:aff id="aff1"><addr-line>Département de Physique and GEPROM, Université de Montréal, Montréal, Québec, Canada</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS Computational Biology</title><idno type="ISSN">1553-734X</idno><idno type="eISSN">1553-7358</idno><imprint><date when="2011">2011</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Chiti, F" uniqKey="Chiti F">F Chiti</name></author><author><name sortKey="Dobson, Cm" uniqKey="Dobson C">CM Dobson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ross, Ca" uniqKey="Ross C">CA Ross</name></author><author><name sortKey="Poirier, Ma" uniqKey="Poirier M">MA Poirier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Miranker, Ad" uniqKey="Miranker A">AD Miranker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zerovnik, E" uniqKey="Zerovnik E">E Zerovnik</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dobson, Cm" uniqKey="Dobson C">CM Dobson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sawaya, Mr" uniqKey="Sawaya M">MR Sawaya</name></author><author><name sortKey="Sambashivan, S" uniqKey="Sambashivan S">S Sambashivan</name></author><author><name sortKey="Nelson, R" uniqKey="Nelson R">R Nelson</name></author><author><name sortKey="Ivanova, Mi" uniqKey="Ivanova M">MI Ivanova</name></author><author><name sortKey="Sievers, Sa" uniqKey="Sievers S">SA Sievers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nelson, R" uniqKey="Nelson R">R Nelson</name></author><author><name sortKey="Sawaya, Mr" uniqKey="Sawaya M">MR Sawaya</name></author><author><name sortKey="Balbirnie, M" uniqKey="Balbirnie M">M Balbirnie</name></author><author><name sortKey="Madsen, A" uniqKey="Madsen A">AØ Madsen</name></author><author><name sortKey="Riekel, C" uniqKey="Riekel C">C Riekel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kayed, R" uniqKey="Kayed R">R Kayed</name></author><author><name sortKey="Head, E" uniqKey="Head E">E Head</name></author><author><name sortKey="Thompson, Jl" uniqKey="Thompson J">JL Thompson</name></author><author><name sortKey="Mcintire, Tm" uniqKey="Mcintire T">TM McIntire</name></author><author><name sortKey="Milton, Sc" uniqKey="Milton S">SC Milton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Thirumalai, D" uniqKey="Thirumalai D">D Thirumalai</name></author><author><name sortKey="Klimov, Dk" uniqKey="Klimov D">DK Klimov</name></author><author><name sortKey="Dima, Ri" uniqKey="Dima R">RI Dima</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sunde, M" uniqKey="Sunde M">M Sunde</name></author><author><name sortKey="Serpell, Lc" uniqKey="Serpell L">LC Serpell</name></author><author><name sortKey="Bartlam, M" uniqKey="Bartlam M">M Bartlam</name></author><author><name sortKey="Fraser, Pe" uniqKey="Fraser P">PE Fraser</name></author><author><name sortKey="Pepys, Mb" uniqKey="Pepys M">MB Pepys</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="F Ndrich, M" uniqKey="F Ndrich M">M Fändrich</name></author><author><name sortKey="Meinhardt, J" uniqKey="Meinhardt J">J Meinhardt</name></author><author><name sortKey="Grigorieff, N" uniqKey="Grigorieff N">N Grigorieff</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Celej, Ms" uniqKey="Celej M">MS Celej</name></author><author><name sortKey="Caarls, W" uniqKey="Caarls W">W Caarls</name></author><author><name sortKey="Demchenko, Ap" uniqKey="Demchenko A">AP Demchenko</name></author><author><name sortKey="Jovin, Tm" uniqKey="Jovin T">TM Jovin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Meinhardt, J" uniqKey="Meinhardt J">J Meinhardt</name></author><author><name sortKey="Sachse, C" uniqKey="Sachse C">C Sachse</name></author><author><name sortKey="Hortschansky, P" uniqKey="Hortschansky P">P Hortschansky</name></author><author><name sortKey="Grigorieff, N" uniqKey="Grigorieff N">N Grigorieff</name></author><author><name sortKey="F Ndrich, M" uniqKey="F Ndrich M">M Fändrich</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Paravastu, Ak" uniqKey="Paravastu A">AK Paravastu</name></author><author><name sortKey="Leapman, Rd" uniqKey="Leapman R">RD Leapman</name></author><author><name sortKey="Yau, Wm" uniqKey="Yau W">WM Yau</name></author><author><name sortKey="Tycko, R" uniqKey="Tycko R">R Tycko</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Madine, J" uniqKey="Madine J">J Madine</name></author><author><name sortKey="Jack, E" uniqKey="Jack E">E Jack</name></author><author><name sortKey="Stockley, Pg" uniqKey="Stockley P">PG Stockley</name></author><author><name sortKey="Radford, Se" uniqKey="Radford S">SE Radford</name></author><author><name sortKey="Serpell, Lc" uniqKey="Serpell L">LC Serpell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Verel, R" uniqKey="Verel R">R Verel</name></author><author><name sortKey="Tomka, It" uniqKey="Tomka I">IT Tomka</name></author><author><name sortKey="Bertozzi, C" uniqKey="Bertozzi C">C Bertozzi</name></author><author><name sortKey="Cadalbert, R" uniqKey="Cadalbert R">R Cadalbert</name></author><author><name sortKey="Kammerer, Ra" uniqKey="Kammerer R">RA Kammerer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kodali, R" uniqKey="Kodali R">R Kodali</name></author><author><name sortKey="Wetzel, R" uniqKey="Wetzel R">R Wetzel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Petkova, At" uniqKey="Petkova A">AT Petkova</name></author><author><name sortKey="Yau, Wm" uniqKey="Yau W">WM Yau</name></author><author><name sortKey="Tycko, R" uniqKey="Tycko R">R Tycko</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kreplak, L" uniqKey="Kreplak L">L Kreplak</name></author><author><name sortKey="Aebi, U" uniqKey="Aebi U">U Aebi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Anderson, M" uniqKey="Anderson M">M Anderson</name></author><author><name sortKey="Bocharova, Ov" uniqKey="Bocharova O">OV Bocharova</name></author><author><name sortKey="Makarava, N" uniqKey="Makarava N">N Makarava</name></author><author><name sortKey="Breydo, L" uniqKey="Breydo L">L Breydo</name></author><author><name sortKey="Salnikov, Vv" uniqKey="Salnikov V">VV Salnikov</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Petkova, At" uniqKey="Petkova A">AT Petkova</name></author><author><name sortKey="Leapman, Rd" uniqKey="Leapman R">RD Leapman</name></author><author><name sortKey="Guo, Wm" uniqKey="Guo W">WM Guo</name></author><author><name sortKey="Mattson, Mp" uniqKey="Mattson M">MP Mattson</name></author><author><name sortKey="Tycko, R" uniqKey="Tycko R">R Tycko</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mo, Y" uniqKey="Mo Y">Y Mo</name></author><author><name sortKey="Lu, Y" uniqKey="Lu Y">Y Lu</name></author><author><name sortKey="Wei, G" uniqKey="Wei G">G Wei</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Park, J" uniqKey="Park J">J Park</name></author><author><name sortKey="Kahng, B" uniqKey="Kahng B">B Kahng</name></author><author><name sortKey="Hwang, W" uniqKey="Hwang W">W Hwang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Berryman, Jt" uniqKey="Berryman J">JT Berryman</name></author><author><name sortKey="Radford, Se" uniqKey="Radford S">SE Radford</name></author><author><name sortKey="Harris, Sa" uniqKey="Harris S">SA Harris</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Klein, Wl" uniqKey="Klein W">WL Klein</name></author><author><name sortKey="Stine, Wb" uniqKey="Stine W">WB Stine</name></author><author><name sortKey="Teplow, Db" uniqKey="Teplow D">DB Teplow</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bucciantini, M" uniqKey="Bucciantini M">M Bucciantini</name></author><author><name sortKey="Giannoni, E" uniqKey="Giannoni E">E Giannoni</name></author><author><name sortKey="Chiti, F" uniqKey="Chiti F">F Chiti</name></author><author><name sortKey="Baroni, F" uniqKey="Baroni F">F Baroni</name></author><author><name sortKey="Formigli, L" uniqKey="Formigli L">L Formigli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shankar, Gm" uniqKey="Shankar G">GM Shankar</name></author><author><name sortKey="Li, S" uniqKey="Li S">S Li</name></author><author><name sortKey="Mehta, Th" uniqKey="Mehta T">TH Mehta</name></author><author><name sortKey="Garcia Munoz, A" uniqKey="Garcia Munoz A">A Garcia-Munoz</name></author><author><name sortKey="Shepardson, Ne" uniqKey="Shepardson N">NE Shepardson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ono, K" uniqKey="Ono K">K Ono</name></author><author><name sortKey="Condron, Mm" uniqKey="Condron M">MM Condron</name></author><author><name sortKey="Teplow, Db" uniqKey="Teplow D">DB Teplow</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lashuel, Ha" uniqKey="Lashuel H">HA Lashuel</name></author><author><name sortKey="Hartley, D" uniqKey="Hartley D">D Hartley</name></author><author><name sortKey="Petre, Bm" uniqKey="Petre B">BM Petre</name></author><author><name sortKey="Walz, T" uniqKey="Walz T">T Walz</name></author><author><name sortKey="Lansbury, Pt" uniqKey="Lansbury P">PT Lansbury</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Goldberg, Ms" uniqKey="Goldberg M">MS Goldberg</name></author><author><name sortKey="Lansbury, Pt" uniqKey="Lansbury P">PT Lansbury</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lansbury, Pt" uniqKey="Lansbury P">PT Lansbury</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zheng, J" uniqKey="Zheng J">J Zheng</name></author><author><name sortKey="Jang, H" uniqKey="Jang H">H Jang</name></author><author><name sortKey="Ma, B" uniqKey="Ma B">B Ma</name></author><author><name sortKey="Nussinov, R" uniqKey="Nussinov R">R Nussinov</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xue, Wf" uniqKey="Xue W">WF Xue</name></author><author><name sortKey="Hellewell, Al" uniqKey="Hellewell A">AL Hellewell</name></author><author><name sortKey="Gosal, Ws" uniqKey="Gosal W">WS Gosal</name></author><author><name sortKey="Homans, Sw" uniqKey="Homans S">SW Homans</name></author><author><name sortKey="Hewitt, Ew" uniqKey="Hewitt E">EW Hewitt</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Knowles, Tp" uniqKey="Knowles T">TP Knowles</name></author><author><name sortKey="Waudby, Ca" uniqKey="Waudby C">CA Waudby</name></author><author><name sortKey="Devlin, Gl" uniqKey="Devlin G">GL Devlin</name></author><author><name sortKey="Cohen, Si" uniqKey="Cohen S">SI Cohen</name></author><author><name sortKey="Aguzzi, A" uniqKey="Aguzzi A">A Aguzzi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Balbirnie, M" uniqKey="Balbirnie M">M Balbirnie</name></author><author><name sortKey="Grothe, R" uniqKey="Grothe R">R Grothe</name></author><author><name sortKey="Eisenberg, Ds" uniqKey="Eisenberg D">DS Eisenberg</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Diaz Avalos, R" uniqKey="Diaz Avalos R">R Diaz-Avalos</name></author><author><name sortKey="Long, C" uniqKey="Long C">C Long</name></author><author><name sortKey="Fontano, E" uniqKey="Fontano E">E Fontano</name></author><author><name sortKey="Balbirnie, M" uniqKey="Balbirnie M">M Balbirnie</name></author><author><name sortKey="Grothe, R" uniqKey="Grothe R">R Grothe</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Der Wel, Pca" uniqKey="Van Der Wel P">PCA Van der Wel</name></author><author><name sortKey="Lewandowski, Jr" uniqKey="Lewandowski J">JR Lewandowski</name></author><author><name sortKey="Griffin R, G" uniqKey="Griffin R G">G Griffin R</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wei, G" uniqKey="Wei G">G Wei</name></author><author><name sortKey="Song, W" uniqKey="Song W">W Song</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Song, W" uniqKey="Song W">W Song</name></author><author><name sortKey="Wei, G" uniqKey="Wei G">G Wei</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Irb Ck, A" uniqKey="Irb Ck A">A Irbäck</name></author><author><name sortKey="Mitternacht, S" uniqKey="Mitternacht S">S Mitternacht</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De Simone, A" uniqKey="De Simone A">A De Simone</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Periole, X" uniqKey="Periole X">X Periole</name></author><author><name sortKey="Rampioni, A" uniqKey="Rampioni A">A Rampioni</name></author><author><name sortKey="Vendruscolo, M" uniqKey="Vendruscolo M">M Vendruscolo</name></author><author><name sortKey="Mark, Ae" uniqKey="Mark A">AE Mark</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, J" uniqKey="Wang J">J Wang</name></author><author><name sortKey="Tan, C" uniqKey="Tan C">C Tan</name></author><author><name sortKey="Chen, Hf" uniqKey="Chen H">HF Chen</name></author><author><name sortKey="Luo, R" uniqKey="Luo R">R Luo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Meli, M" uniqKey="Meli M">M Meli</name></author><author><name sortKey="Morra, G" uniqKey="Morra G">G Morra</name></author><author><name sortKey="Colombo, G" uniqKey="Colombo G">G Colombo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name></author><author><name sortKey="Chen, H" uniqKey="Chen H">H Chen</name></author><author><name sortKey="Bai, H" uniqKey="Bai H">H Bai</name></author><author><name sortKey="Lai, L" uniqKey="Lai L">L Lai</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Strodel, B" uniqKey="Strodel B">B Strodel</name></author><author><name sortKey="Whittleston, Cs" uniqKey="Whittleston C">CS Whittleston</name></author><author><name sortKey="Wales, Dj" uniqKey="Wales D">DJ Wales</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Reddy, As" uniqKey="Reddy A">AS Reddy</name></author><author><name sortKey="Chopra, M" uniqKey="Chopra M">M Chopra</name></author><author><name sortKey="De Pablo, Jj" uniqKey="De Pablo J">JJ de Pablo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vitagliano, L" uniqKey="Vitagliano L">L Vitagliano</name></author><author><name sortKey="Esposito, L" uniqKey="Esposito L">L Esposito</name></author><author><name sortKey="Pedone, C" uniqKey="Pedone C">C Pedone</name></author><author><name sortKey="De Simone, A" uniqKey="De Simone A">A De Simone</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zheng, J" uniqKey="Zheng J">J Zheng</name></author><author><name sortKey="Ma, B" uniqKey="Ma B">B Ma</name></author><author><name sortKey="Tsai, Cj" uniqKey="Tsai C">CJ Tsai</name></author><author><name sortKey="Nussinov, R" uniqKey="Nussinov R">R Nussinov</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Esposito, L" uniqKey="Esposito L">L Esposito</name></author><author><name sortKey="Pedone, C" uniqKey="Pedone C">C Pedone</name></author><author><name sortKey="Vitagliano, L" uniqKey="Vitagliano L">L Vitagliano</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gsponer, J" uniqKey="Gsponer J">J Gsponer</name></author><author><name sortKey="Haberthur, U" uniqKey="Haberthur U">U Haberthür</name></author><author><name sortKey="Caflisch, A" uniqKey="Caflisch A">A Caflisch</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sugita, Y" uniqKey="Sugita Y">Y Sugita</name></author><author><name sortKey="Okamoto, Y" uniqKey="Okamoto Y">Y Okamoto</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maupetit, J" uniqKey="Maupetit J">J Maupetit</name></author><author><name sortKey="Tuffery, P" uniqKey="Tuffery P">P Tufféry</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hess, B" uniqKey="Hess B">B Hess</name></author><author><name sortKey="Kutzner, C" uniqKey="Kutzner C">C Kutzner</name></author><author><name sortKey="Van Der Spoel, D" uniqKey="Van Der Spoel D">D Van Der Spoel</name></author><author><name sortKey="Lindahl, E" uniqKey="Lindahl E">E Lindahl</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, W" uniqKey="Chen W">W Chen</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Melquiond, A" uniqKey="Melquiond A">A Melquiond</name></author><author><name sortKey="Dong, X" uniqKey="Dong X">X Dong</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maupetit, J" uniqKey="Maupetit J">J Maupetit</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author><author><name sortKey="Tuffery, P" uniqKey="Tuffery P">P Tufféry</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maupetit, J" uniqKey="Maupetit J">J Maupetit</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author><author><name sortKey="Tuffery, P" uniqKey="Tuffery P">P Tufféry</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Berendsen, Hjc" uniqKey="Berendsen H">HJC Berendsen</name></author><author><name sortKey="Postma, Jpm" uniqKey="Postma J">JPM Postma</name></author><author><name sortKey="Van Gunsteren, Wf" uniqKey="Van Gunsteren W">WF van Gunsteren</name></author><author><name sortKey="Dinola, A" uniqKey="Dinola A">A Dinola</name></author><author><name sortKey="Haak, Jr" uniqKey="Haak J">JR Haak</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chebaro, Y" uniqKey="Chebaro Y">Y Chebaro</name></author><author><name sortKey="Dong, X" uniqKey="Dong X">X Dong</name></author><author><name sortKey="Laghaei, R" uniqKey="Laghaei R">R Laghaei</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Carter, P" uniqKey="Carter P">P Carter</name></author><author><name sortKey="Andersen, Ca" uniqKey="Andersen C">CA Andersen</name></author><author><name sortKey="Rost, B" uniqKey="Rost B">B Rost</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, L" uniqKey="Yang L">L Yang</name></author><author><name sortKey="Shao, Q" uniqKey="Shao Q">Q Shao</name></author><author><name sortKey="Goa, Yq" uniqKey="Goa Y">YQ Goa</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Holm, L" uniqKey="Holm L">L Holm</name></author><author><name sortKey="Sander, C" uniqKey="Sander C">C Sander</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hartmann, C" uniqKey="Hartmann C">C Hartmann</name></author><author><name sortKey="Antes, I" uniqKey="Antes I">I Antes</name></author><author><name sortKey="Lengauer, T" uniqKey="Lengauer T">T Lengauer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dunbrack, Rl" uniqKey="Dunbrack R">RL Dunbrack</name></author><author><name sortKey="Cohen, Fe" uniqKey="Cohen F">FE Cohen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Berendsen, Hjc" uniqKey="Berendsen H">HJC Berendsen</name></author><author><name sortKey="Grigera, Jr" uniqKey="Grigera J">JR Grigera</name></author><author><name sortKey="Straatsma, Tp" uniqKey="Straatsma T">TP Straatsma</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hess, B" uniqKey="Hess B">B Hess</name></author><author><name sortKey="Bekker, H" uniqKey="Bekker H">H Bekker</name></author><author><name sortKey="Berendsen, Hjc" uniqKey="Berendsen H">HJC Berendsen</name></author><author><name sortKey="Fraaije, Jgem" uniqKey="Fraaije J">JGEM Fraaije</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Miyamoto, S" uniqKey="Miyamoto S">S Miyamoto</name></author><author><name sortKey="Kollman, Pa" uniqKey="Kollman P">PA Kollman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Berendsen, Hjc" uniqKey="Berendsen H">HJC Berendsen</name></author><author><name sortKey="Van Der Spoel, D" uniqKey="Van Der Spoel D">D van der Spoel</name></author><author><name sortKey="Van Drunen, R" uniqKey="Van Drunen R">R van Drunen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lindahl, E" uniqKey="Lindahl E">E Lindahl</name></author><author><name sortKey="Hess, B" uniqKey="Hess B">B Hess</name></author><author><name sortKey="Van Der Spoel, D" uniqKey="Van Der Spoel D">D van der Spoel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Der Spoel, D" uniqKey="Van Der Spoel D">D van der Spoel</name></author><author><name sortKey="Lindahl, E" uniqKey="Lindahl E">E Lindahl</name></author><author><name sortKey="Hess, B" uniqKey="Hess B">B Hess</name></author><author><name sortKey="Van Buuren, Ar" uniqKey="Van Buuren A">AR van Buuren</name></author><author><name sortKey="Apol, E" uniqKey="Apol E">E Apol</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Gunsteren, Wf" uniqKey="Van Gunsteren W">WF van Gunsteren</name></author><author><name sortKey="Berendsen, Hjc" uniqKey="Berendsen H">HJC Berendsen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Gunsteren, Wf" uniqKey="Van Gunsteren W">WF van Gunsteren</name></author><author><name sortKey="Daura, Wfx" uniqKey="Daura W">WFX Daura</name></author><author><name sortKey="Mark, Ae" uniqKey="Mark A">AE Mark</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Scott, Wrp" uniqKey="Scott W">WRP Scott</name></author><author><name sortKey="Hunenberger, Ph" uniqKey="Hunenberger P">PH Hunenberger</name></author><author><name sortKey="Tironi, Ig" uniqKey="Tironi I">IG Tironi</name></author><author><name sortKey="Mark, Ae" uniqKey="Mark A">AE Mark</name></author><author><name sortKey="Billeter, Sr" uniqKey="Billeter S">SR Billeter</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Gunsteren, Wf" uniqKey="Van Gunsteren W">WF van Gunsteren</name></author><author><name sortKey="Billeter, Sr" uniqKey="Billeter S">SR Billeter</name></author><author><name sortKey="Eising, Aa" uniqKey="Eising A">AA Eising</name></author><author><name sortKey="Huenenberger, Ph" uniqKey="Huenenberger P">PH Huenenberger</name></author><author><name sortKey="Krueger, P" uniqKey="Krueger P">P Krueger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, P" uniqKey="Liu P">P Liu</name></author><author><name sortKey="Friesner, Ra" uniqKey="Friesner R">RA Friesner</name></author><author><name sortKey="Berne, Bj" uniqKey="Berne B">BJ Berne</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kubitzki, Mb" uniqKey="Kubitzki M">MB Kubitzki</name></author><author><name sortKey="De Groot, Bl" uniqKey="De Groot B">BL de Groot</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Der Spoel, D" uniqKey="Van Der Spoel D">D van der Spoel</name></author><author><name sortKey="Seibert, Mm" uniqKey="Seibert M">MM Seibert</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wei, G" uniqKey="Wei G">G Wei</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Santini, S" uniqKey="Santini S">S Santini</name></author><author><name sortKey="Wei, G" uniqKey="Wei G">G Wei</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pellarin, R" uniqKey="Pellarin R">R Pellarin</name></author><author><name sortKey="Schuetz, P" uniqKey="Schuetz P">P Schuetz</name></author><author><name sortKey="Guarnera, E" uniqKey="Guarnera E">E Guarnera</name></author><author><name sortKey="Caflisch, A" uniqKey="Caflisch A">A Caflisch</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Santini, S" uniqKey="Santini S">S Santini</name></author><author><name sortKey="Mousseau, N" uniqKey="Mousseau N">N Mousseau</name></author><author><name sortKey="Derreumaux, P" uniqKey="Derreumaux P">P Derreumaux</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, Dw" uniqKey="Li D">DW Li</name></author><author><name sortKey="Mohanty, S" uniqKey="Mohanty S">S Mohanty</name></author><author><name sortKey="Irb Ck, A" uniqKey="Irb Ck A">A Irbäck</name></author><author><name sortKey="Huo, S" uniqKey="Huo S">S Huo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tycko, R" uniqKey="Tycko R">R Tycko</name></author><author><name sortKey="Sciarretta, Kl" uniqKey="Sciarretta K">KL Sciarretta</name></author><author><name sortKey="Orgel, Jp" uniqKey="Orgel J">JP Orgel</name></author><author><name sortKey="Meredith, Sc" uniqKey="Meredith S">SC Meredith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Salemme, Fr" uniqKey="Salemme F">FR Salemme</name></author><author><name sortKey="Weatherford, Dw" uniqKey="Weatherford D">DW Weatherford</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Salemme, Fr" uniqKey="Salemme F">FR Salemme</name></author><author><name sortKey="Weatherford, Dw" uniqKey="Weatherford D">DW Weatherford</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Antzutkin, On" uniqKey="Antzutkin O">ON Antzutkin</name></author><author><name sortKey="Balbach, Jj" uniqKey="Balbach J">JJ Balbach</name></author><author><name sortKey="Leapman, Rd" uniqKey="Leapman R">RD Leapman</name></author><author><name sortKey="Rizzo, Nw" uniqKey="Rizzo N">NW Rizzo</name></author><author><name sortKey="Reed, J" uniqKey="Reed J">J Reed</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shorter, J" uniqKey="Shorter J">J Shorter</name></author><author><name sortKey="Lindquist, S" uniqKey="Lindquist S">S Lindquist</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Behrends, C" uniqKey="Behrends C">C Behrends</name></author><author><name sortKey="Langer, Ca" uniqKey="Langer C">CA Langer</name></author><author><name sortKey="Boteva, R" uniqKey="Boteva R">R Boteva</name></author><author><name sortKey="Bottcher, Um" uniqKey="Bottcher U">UM Böttcher</name></author><author><name sortKey="Stemp, Mj" uniqKey="Stemp M">MJ Stemp</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ehrnhoefer, De" uniqKey="Ehrnhoefer D">DE Ehrnhoefer</name></author><author><name sortKey="Bieschke, J" uniqKey="Bieschke J">J Bieschke</name></author><author><name sortKey="Boeddrich, A" uniqKey="Boeddrich A">A Boeddrich</name></author><author><name sortKey="Herbst, M" uniqKey="Herbst M">M Herbst</name></author><author><name sortKey="Masino, L" uniqKey="Masino L">L Masino</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS Comput Biol</journal-id><journal-id journal-id-type="iso-abbrev">PLoS Comput. Biol</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">ploscomp</journal-id><journal-title-group><journal-title>PLoS Computational Biology</journal-title></journal-title-group><issn pub-type="ppub">1553-734X</issn><issn pub-type="epub">1553-7358</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">21625573</article-id><article-id pub-id-type="pmc">3098217</article-id><article-id pub-id-type="publisher-id">PCOMPBIOL-D-10-00055</article-id><article-id pub-id-type="doi">10.1371/journal.pcbi.1002051</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology</subject><subj-group><subject>Biochemistry</subject><subj-group><subject>Biochemistry Simulations</subject><subject>Proteins</subject></subj-group></subj-group><subj-group><subject>Biophysics</subject><subj-group><subject>Biophysics Simulations</subject><subject>Macromolecular Assemblies</subject></subj-group></subj-group><subj-group><subject>Computational Biology</subject><subj-group><subject>Biochemical Simulations</subject><subject>Biophysic Al Simulations</subject><subject>Macromolecular Structure Analysis</subject></subj-group></subj-group></subj-group></article-categories><title-group><article-title>A Multiscale Approach to Characterize the Early Aggregation Steps of the Amyloid-Forming Peptide GNNQQNY from the Yeast Prion Sup-35</article-title><alt-title alt-title-type="running-head">Study of GNNQQNY Peptide by Numerical Methods</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Nasica-Labouze</surname><given-names>Jessica</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Meli</surname><given-names>Massimiliano</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Derreumaux</surname><given-names>Philippe</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Colombo</surname><given-names>Giorgio</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author"><name><surname>Mousseau</surname><given-names>Normand</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Département de Physique and GEPROM, Université de Montréal, Montréal, Québec, Canada</addr-line></aff><aff id="aff2"><label>2</label><addr-line>Istituto di Chimica del Riconoscimento Molecolare, CNR, Milano, Italy</addr-line></aff><aff id="aff3"><label>3</label><addr-line>Laboratoire de Biochimie Théorique, UPR9080 CNRS, Institut de Biologie Physico-Chimique, Université Paris 7, and Institut Universitaire de France, Paris, France</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Pande</surname><given-names>Vijay S.</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1">Stanford University, United States of America</aff><author-notes><corresp id="cor1">* E-mail: <email>g.colombo@icrm.cnr.it</email> (GC); <email>normand.mousseau@umontreal.ca</email> (NM)</corresp><fn fn-type="con"><p>Conceived and designed the experiments: JNL MM PD GC NM. Performed the experiments: JNL MM. Analyzed the data: JNL MM PD GC NM. Contributed reagents/materials/analysis tools: JNL MM PD GC NM. Wrote the paper: JNL MM PD GC NM.</p></fn></author-notes><pub-date pub-type="collection"><month>5</month><year>2011</year></pub-date><pub-date pub-type="epub"><day>19</day><month>5</month><year>2011</year></pub-date><volume>7</volume><issue>5</issue><elocation-id>e1002051</elocation-id><history><date date-type="received"><day>5</day><month>10</month><year>2010</year></date><date date-type="accepted"><day>28</day><month>3</month><year>2011</year></date></history><permissions><copyright-statement>Nasica-Labouze et al.</copyright-statement><copyright-year>2011</copyright-year><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.</license-p></license></permissions><abstract><p>The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.</p></abstract><abstract abstract-type="summary"><title>Author Summary</title><p>The formation of amyloid fibrils is associated with many neurodegenerative diseases such as Alzheimer's, Creutzfeld-Jakob, Parkinson's, the Prion disease and diabetes mellitus. In all cases, proteins misfold to form highly ordered insoluble aggregates called amyloid fibrils that deposit intra- and extracellularly and are resistant to proteases. All proteins are believed to have the instrinsic capability of forming amyloid fibrils that share common specific structural properties that have been observed by X-ray crystallography and by NMR. However, little is known about the aggregation dynamics of amyloid assemblies, and their toxicity mechanism is therefore poorly understood. It is believed that small amyloid oligomers, formed on the aggregation pathway of full amyloid fibrils, are the toxic species. A detailed atomic characterization of the oligomerization process is thus necessary to further our understanding of the amyloid oligomer's toxicity. Our approach here is to study the aggregation dynamics of a 7-residue amyloid peptide GNNQQNY through a combination of numerical techniques. Our results suggest that this amyloid sequence can form fibril-like structures and is polymorphic, which agrees with recent experimental observations. The ability to fully characterize and describe the aggregation pathway of amyloid sequences numerically is key to the development of future drugs to target amyloid oligomers.</p></abstract><counts><page-count count="18"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>The aggregation of soluble peptides and proteins first into soluble oligomeric assemblies and then into insoluble amyloid fibrils is associated with the onset of misfolding diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes and transmissible spongiform encephalopathies <xref rid="pcbi.1002051-Chiti1" ref-type="bibr">[1]</xref>–<xref rid="pcbi.1002051-Dobson1" ref-type="bibr">[5]</xref>. Though there is no sequence similarity, the final products of all amyloidogenic proteins display a similar cross-β structure <xref rid="pcbi.1002051-Sawaya1" ref-type="bibr">[6]</xref>, <xref rid="pcbi.1002051-Nelson1" ref-type="bibr">[7]</xref> and the soluble oligomers of several proteins appear to share similar structural properties <xref rid="pcbi.1002051-Kayed1" ref-type="bibr">[8]</xref>, suggesting common pathways for amyloid formation <xref rid="pcbi.1002051-Kayed1" ref-type="bibr">[8]</xref>–<xref rid="pcbi.1002051-Sunde1" ref-type="bibr">[10]</xref>. Structural similarity does not, however, exclude diversity or polymorphism in the intermediates and products of amyloid assembly <xref rid="pcbi.1002051-Fndrich1" ref-type="bibr">[11]</xref>–<xref rid="pcbi.1002051-Berryman1" ref-type="bibr">[24]</xref>.</p><p>Many studies have shown that soluble oligomeric intermediates are more toxic than the full fibrils themselves <xref rid="pcbi.1002051-Klein1" ref-type="bibr">[25]</xref>, <xref rid="pcbi.1002051-Bucciantini1" ref-type="bibr">[26]</xref>. These transient oligomers include low molecular weight aggregates (e.g. dimers <xref rid="pcbi.1002051-Shankar1" ref-type="bibr">[27]</xref> and tetramers <xref rid="pcbi.1002051-Ono1" ref-type="bibr">[28]</xref>) and high molecular weight species (e.g., β-sheet rich annular protofibrils similar to pore-forming toxins <xref rid="pcbi.1002051-Lashuel1" ref-type="bibr">[29]</xref>–<xref rid="pcbi.1002051-Zheng1" ref-type="bibr">[32]</xref>). While oligomers are considered as primary toxic species for most neurodegenerative diseases, there is recent experimental evidence that fragmentation or breakage of fibrils can contribute to the kinetics of aggregation and the amyloid cytotoxicity itself <xref rid="pcbi.1002051-Xue1" ref-type="bibr">[33]</xref>, <xref rid="pcbi.1002051-Knowles1" ref-type="bibr">[34]</xref>.</p><p>One important way for investigating amyloid fibril formation, polymorphism and cytotoxicity is offered by short protein fragments. Among them, GNNQQNY, from the N-terminal prion-determining domain of the yeast protein Sup35, is a paradigmatic example of a short sequence with the same properties as its corresponding full-length protein <xref rid="pcbi.1002051-Balbirnie1" ref-type="bibr">[35]</xref>, <xref rid="pcbi.1002051-DiazAvalos1" ref-type="bibr">[36]</xref>. These properties include an amyloid fibril with a core cross-β spine, Congo-red binding and a nucleated-growth aggregation process <xref rid="pcbi.1002051-Balbirnie1" ref-type="bibr">[35]</xref>. In particular, X-ray diffraction of several micro-crystals provides a detailed atomic structure for different GNNQQNY fibrillar morphologies where the side-chains form self-complementing steric zippers <xref rid="pcbi.1002051-Sawaya1" ref-type="bibr">[6]</xref>, <xref rid="pcbi.1002051-Nelson1" ref-type="bibr">[7]</xref>, <xref rid="pcbi.1002051-Balbirnie1" ref-type="bibr">[35]</xref>–<xref rid="pcbi.1002051-VanderWel1" ref-type="bibr">[37]</xref>. As for all amyloid sequences, structural characterization of the intermediate GNNQQNY oligomers has been however precluded experimentally due to the high complexity of the aggregation process, and the short-lived and meta-stable character of the early aggregates.</p><p>Computer simulations have proved useful complements to experiments for looking at the initial aggregation steps providing information, for example, about the presence of amorphous states in dynamic equilibrium with fibrillar and annular states <xref rid="pcbi.1002051-Wei1" ref-type="bibr">[38]</xref>–<xref rid="pcbi.1002051-DeSimone1" ref-type="bibr">[41]</xref> and the final steps of the polymerization-nucleation process <xref rid="pcbi.1002051-Park1" ref-type="bibr">[23]</xref>, <xref rid="pcbi.1002051-Periole1" ref-type="bibr">[42]</xref>. They can provide atomic-resolution insights into several factors, ranging from the effect of sequence variations on aggregation tendencies to information on the stability of aggregates and the kinetics of aggregation. Due to lighter computational costs, short peptides are more amenable to simulations of the aggregation process than full-length proteins. For example, a number of numerical works have characterized the structures and free energy of small GNNQQNY aggregates ranging from 2-mers to 8-mers starting from disordered states or studied the stability of pre-formed GNNQQNY assemblies with cross-β or annular morphologies <xref rid="pcbi.1002051-Park1" ref-type="bibr">[23]</xref>, <xref rid="pcbi.1002051-Berryman1" ref-type="bibr">[24]</xref>, <xref rid="pcbi.1002051-Zheng1" ref-type="bibr">[32]</xref>, <xref rid="pcbi.1002051-Periole1" ref-type="bibr">[42]</xref>–<xref rid="pcbi.1002051-Gsponer1" ref-type="bibr">[51]</xref>.</p><p>In this paper, we push the boundaries of the GNNQQNY oligomer size and investigate, through a multi-scale simulation approach, the aggregation and polymorphism of three GNNQQNY oligomer sizes: 3-mers, 12-mers and 20-mers. Our approach takes advantage of the accelerated sampling properties of replica exchange molecular dynamics (REMD) simulations <xref rid="pcbi.1002051-Sugita1" ref-type="bibr">[52]</xref> based on coarse-grained models and of the accurate description of the physico-chemical interactions between the peptides and the solvent by using an all-atom model. More precisely, we first use REMD simulations <xref rid="pcbi.1002051-Sugita1" ref-type="bibr">[52]</xref> with the coarse-grained potential OPEP <xref rid="pcbi.1002051-Maupetit1" ref-type="bibr">[53]</xref>, <xref rid="pcbi.1002051-Derreumaux1" ref-type="bibr">[54]</xref>, and then analyze the stability and conformational properties of selected aggregates by room temperature MD as well as REMD simulations using the GROMOS force-field <xref rid="pcbi.1002051-Hess1" ref-type="bibr">[55]</xref>. In total, we accumulated more than 23.60 µs and 2.66 µs of coarse-grained and all-atom simulations, respectively, allowing relevant statistical analysis. To our knowledge, the present study reports the largest simulations of spontaneous self-organization carried out at the atomic resolution on an amyloid peptide without any pre-formed seed. Overall, the results of our integrated simulations and analysis show the existence of a high degree of polymorphism for the GNNQQNY sequence, even for oligomeric assemblies containing as many as 20 monomers.</p></sec><sec sec-type="materials|methods" id="s2"><title>Materials and Methods</title><p>Simulations and analyses presented here couple a number of approaches, which are described briefly in this section. The first set of simulations uses the coarse-grained OPEP potential with replica-exchange molecular dynamics (REMD). These are followed by all-atom simulations using GROMACS with MD and REMD. All simulations are labeled as follows: a number, which indicates the number of monomers, two letters indicating the force field (OP for OPEP and GR for GROMACS), a letter or number indicating the simulation and a label for the specific conformation studied (when appropriate) giving, for example: 01OP2-A1.</p><sec id="s2a"><title>Replica-Exchange Molecular Dynamics (REMD)</title><p>REMD is a thermodynamical sampling method that requires the running of N MD trajectories (or replica) in parallel at N different temperatures selected in order to optimize thermodynamical sampling <xref rid="pcbi.1002051-Sugita1" ref-type="bibr">[52]</xref>. At regular time intervals, conformational exchanges are attempted between adjacent simulation pairs according to the Metropolis criterion with accept-reject probability:<disp-formula><graphic xlink:href="pcbi.1002051.e001.jpg" mimetype="image" position="float"></graphic></disp-formula>where, before the exchange, trajectory <italic>i</italic> at temperature <italic>T<sub>i</sub></italic> has an energy <italic>E<sub>i</sub></italic> and trajectory <italic>j</italic> has an energy <italic>E<sub>j</sub></italic> at temperature <italic>T<sub>j</sub></italic>.</p><p>This broadly used method allows for conformations in a deep local minimum to explore other regions of the energy landscape by migrating to higher temperatures. While thermodynamical properties converge faster than with single temperature standard MD, dynamical information is lost due to temperature exchanges. It is still possible, however, to derive thermodynamically putative aggregation pathways by following the continuous trajectories through temperature space.</p></sec><sec id="s2b"><title>The Optimal Potential For Efficient Peptide-Structure Prediction (OPEP) Force-Field</title><p>OPEP is a coarse-grained protein model that uses a detailed representation of all backbone atoms (N, H, Cα, C and O) and reduces each side-chain to one single bead with appropriate geometrical parameters and van der Waals radius. The OPEP energy function, which includes implicit effects of aqueous solution, is expressed as a sum of local potentials (taking into account the changes in bond lengths, bond angles, improper torsions of the side-chains and backbone torsions), non-bonded potentials (taking into account the hydrophobic and hydrophilic properties of each amino acid) and hydrogen-bonding potentials (taking into account two- and four- body interactions) <xref rid="pcbi.1002051-Maupetit1" ref-type="bibr">[53]</xref>. OPEP has been extensively tested on peptides using multiple approaches such as the activation-relaxation technique <xref rid="pcbi.1002051-Chen1" ref-type="bibr">[56]</xref>, Monte Carlo <xref rid="pcbi.1002051-Derreumaux1" ref-type="bibr">[54]</xref>, MD <xref rid="pcbi.1002051-Song1" ref-type="bibr">[39]</xref> and REMD simulations <xref rid="pcbi.1002051-Melquiond1" ref-type="bibr">[57]</xref>, and greedy-based algorithms <xref rid="pcbi.1002051-Maupetit2" ref-type="bibr">[58]</xref>, <xref rid="pcbi.1002051-Maupetit3" ref-type="bibr">[59]</xref>. OPEP is also appropriate for simulations of GNNQQNY. Preliminary test simulations on this peptide's dimer indicate that, at 300 K, the GNNQQNY relative orientation is a 60 to 40 probability in favor of the antiparallel dimer with a least two hydrogen bonds. This result is in general agreement with what was found by Strodel <italic>et al.</italic> with CHARMM19 and the implicit solvation potential EEFI <xref rid="pcbi.1002051-Strodel1" ref-type="bibr">[46]</xref> where both orientations of the strands are visited with similar probabilities.</p></sec><sec id="s2c"><title>OPEP Simulation Details</title><p>REMD were carried out using a 1.5 fs time-step, periodic boundary conditions with box sizes depending on the systems and a weak coupling to an external temperature bath <xref rid="pcbi.1002051-Berendsen1" ref-type="bibr">[60]</xref>, <xref rid="pcbi.1002051-Derreumaux2" ref-type="bibr">[61]</xref>, <xref rid="pcbi.1002051-Chebaro1" ref-type="bibr">[62]</xref>. Replica exchanges were attempted every 5000 steps and configurations saved every 5000 steps. Initial structures for 3-mer and the 20-mer simulations were constructed by placing random coil monomers between 12 Å to 50 Å apart (<xref ref-type="fig" rid="pcbi-1002051-g001"><bold>Figure 1a</bold></xref><bold>and</bold><xref ref-type="fig" rid="pcbi-1002051-g002"><bold>Figure 2</bold></xref>). For the 12-mer, the initial chains occupied four rows, with each peptide separated from the others by 15 Å (<xref ref-type="fig" rid="pcbi-1002051-g001"><bold>Figure 1b</bold></xref>). Because of the extensive sampling of REMD, all results are independent from this initial setup. For the 20-mer system, three OPEP-REMD simulations were launched. A preliminary REMD simulation (20OPp) was used to obtain a first estimate of the melting temperature (Tm = 283 K) and generate some representative conformations for all-atom MD refinement. The 20 initial temperatures were logarithmically distributed between 230 K and 450 K. Despite 200-ns simulation per replica, we found that the configuration space was not optimally sampled because of the existence of a large discontinuity in the potential energy when the system orders. Thus, a Gaussian distribution of temperatures around 283 K was deemed preferable to allow a better sampling of the phase space. The other two REMD simulations, running for 400 ns at each temperature, were started from the same random configuration (<xref ref-type="fig" rid="pcbi-1002051-g002"><bold>Figure 2</bold></xref>), but with an optimized Gaussian temperature distribution centered around 283 K: 20OP1 uses 20 temperatures (in Kelvins: 223.8, 249.2, 260.1, 266.0, 270.3, 273.8, 277.1, 280.1, 283.0, 285.9, 288.9, 292.2, 295.7, 300.1, 305.9, 316.8, 342.2, 370.1, 398.0, 426.0) and 20OP2 uses 22 temperatures, with two more temperatures below the transition, at 236.5 and 254.7 K, to increase exchanges between low-energy structures. All REMD simulations are summarized in <xref ref-type="table" rid="pcbi-1002051-t001"><bold>Table 1</bold></xref>.</p><fig id="pcbi-1002051-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g001</object-id><label>Figure 1</label><caption><title>Starting structures for a) the trimer and b) dodecamer.</title><p>The concentration for both systems is set at 4.15 mM. The random coil monomers are placed 15 Å apart.</p></caption><graphic xlink:href="pcbi.1002051.g001"></graphic></fig><fig id="pcbi-1002051-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g002</object-id><label>Figure 2</label><caption><title>Starting structure with random coils and no seed for the 20-mer simulations.</title><p>The concentration is also 4.15 mM. The monomers are randomly placed 12 to 50 Å apart.</p></caption><graphic xlink:href="pcbi.1002051.g002"></graphic></fig><table-wrap id="pcbi-1002051-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.t001</object-id><label>Table 1</label><caption><title>Details of all simulations run for the trimer, dodecamer and 20-mer systems.</title></caption><alternatives><graphic id="pcbi-1002051-t001-1" xlink:href="pcbi.1002051.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Length of OPEP simulations (ns)<xref ref-type="table-fn" rid="nt102">(a)</xref></td><td align="left" rowspan="1" colspan="1">OPEP - LABEL of structures extracted<xref ref-type="table-fn" rid="nt103">(b)</xref></td><td align="left" rowspan="1" colspan="1">Temperatures min-max (K) & number of temperatures<xref ref-type="table-fn" rid="nt104">(c)</xref></td><td align="left" rowspan="1" colspan="1">GROMACS - LABEL of reconstructed OPEP extracted structures<xref ref-type="table-fn" rid="nt103">(b)</xref></td><td align="left" rowspan="1" colspan="1">Length of GROMACS simulations (ns)<xref ref-type="table-fn" rid="nt105">(d)</xref></td><td align="left" rowspan="1" colspan="1">Total number of atoms<xref ref-type="table-fn" rid="nt106">(e)</xref></td><td align="left" rowspan="1" colspan="1">Temperature (K)<xref ref-type="table-fn" rid="nt107">(f)</xref></td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">3-mer</td><td align="left" rowspan="1" colspan="1">50×16</td><td align="left" rowspan="1" colspan="1">03OP1-A</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">03-GR1-A</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">6504</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">03OP1-B</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">03-GR1-B</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">5997</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">03OP1-C</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">03-GR1-C</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">5619</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">03OP1-D</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">03-GR1-D</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">5949</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">03OP1-E</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">03-GR1-E</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">6621</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1">12-mer</td><td align="left" rowspan="1" colspan="1">125×16</td><td align="left" rowspan="1" colspan="1">12OP1-A</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">12-GR1-A</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">56271</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">12OP1-B</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">12-GR1-B</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">40476</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">12OP1-C</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">12-GR1-C</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">15381</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">12OP1-D</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">12-GR1-D</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">16134</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">12OP1-E</td><td align="left" rowspan="1" colspan="1">222.5–525 16</td><td align="left" rowspan="1" colspan="1">12-GR1-E</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">17157</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1">20-mer OPp</td><td align="left" rowspan="1" colspan="1">200×20</td><td align="left" rowspan="1" colspan="1">20OPp-A</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-A1</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">25620</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OPp-A</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-A2</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">25620</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OPp-B</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-B1</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">27816</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OPp-B</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-B2</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">27816</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OPp-B</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-B3 (REMD)</td><td align="left" rowspan="1" colspan="1">10×12</td><td align="left" rowspan="1" colspan="1">27816</td><td align="left" rowspan="1" colspan="1">see text</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OPp-C</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-C1</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">49674</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OPp-C</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-C2</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">49674</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OPp-D</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-D1</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">23746</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OPp-D</td><td align="left" rowspan="1" colspan="1">234.6–447.6 20</td><td align="left" rowspan="1" colspan="1">20GRp-D2</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">23746</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1">20-mer OP2</td><td align="left" rowspan="1" colspan="1">400×22</td><td align="left" rowspan="1" colspan="1">20OP2-A</td><td align="left" rowspan="1" colspan="1">223.8–425.9 22</td><td align="left" rowspan="1" colspan="1">20GR2-A</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">21586</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OP2-B</td><td align="left" rowspan="1" colspan="1">223.8–425.9 22</td><td align="left" rowspan="1" colspan="1">20GR2-B</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">24640</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OP2-C</td><td align="left" rowspan="1" colspan="1">223.8–425.9 22</td><td align="left" rowspan="1" colspan="1">20GR2-C</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">23554</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OP2-E</td><td align="left" rowspan="1" colspan="1">223.8–425.9 22</td><td align="left" rowspan="1" colspan="1">20GR2-E</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">26065</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OP2-N</td><td align="left" rowspan="1" colspan="1">223.8–425.9 22</td><td align="left" rowspan="1" colspan="1">20GR2-N1</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">32206</td><td align="left" rowspan="1" colspan="1">300</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20OP2-N</td><td align="left" rowspan="1" colspan="1">223.8–425.9 22</td><td align="left" rowspan="1" colspan="1">20GR2-N2 (REMD)</td><td align="left" rowspan="1" colspan="1">10×12</td><td align="left" rowspan="1" colspan="1">32206</td><td align="left" rowspan="1" colspan="1">see text</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt101"><label></label><p>This table presents simulations done with OPEP (coarse-grained potential) and GROMACS (all-atom potential).</p></fn><fn id="nt102"><label>(a)</label><p>The total simulation time for OPEP REMD simulations in the format time_per_replica x number_of_replicas.</p></fn><fn id="nt103"><label>(b)</label><p>The label of the OPEP/GROMACS structures extracted. The label indicates the number of monomers, the potential used (OP for OPEP and GR for GROMACS), the simulation index (1,2 or p (preliminary)) and the letter ID of the structure.</p></fn><fn id="nt104"><label>(c)</label><p>The range of temperatures (in K) used for OPEP REMD simulations.</p></fn><fn id="nt105"><label>(d)</label><p>The total simulation time for GROMACS simulations. MD simulations are indicated by only one number while, for REMD simulations, the total simulation time is given in the format time_per_replica x number_of_replicas.</p></fn><fn id="nt106"><label>(e)</label><p>The total number of atoms in the system including protein and solvatation water atoms.</p></fn><fn id="nt107"><label>(f)</label><p>The temperature used in GROMACS simulations (in K).</p></fn></table-wrap-foot></table-wrap><p>Determining whether equilibrium has been reached, even for the trimer, is difficult. It is always possible that a system is stuck in a minimum and thermodynamical properties will then appear as though they are converged. Here, we use the specific heat to track convergence. This quantity, the second derivative of the free energy, is very sensitive to convergence at all temperatures, and provides a very stringent test even near transitions. Because we are mostly interested in the qualitative properties of the systems under study here, we consider that a system is converged when the overall shape of the specific heat near the transition is converged. This ensures that the dominant structures are found with the proper weight, within the limits of our simulations.</p></sec><sec id="s2d"><title>OPEP Analysis and Structure Selection</title><p>Analysis for these simulations was performed, in part, using a new clustering code that enables us to identify the dominant configuration types in terms of clusters formed in β-sheet structures based on strand attachment. The criterion set to define a hydrogen bond between two given strands is similar to the one used in the DSSP algorithm. <xref rid="pcbi.1002051-Carter1" ref-type="bibr">[63]</xref>. A cutoff of one hydrogen bond is used for distinguishing random from β-strands since we are dealing with a very short sequence and not considering the hydrogen bonds with the N-terminal glycines. The configuration types are defined here in terms of the number of sheets and the number of strands per sheet in the structure. For instance, a configuration type 8 7 5 for the 20-mer describes a structure with 3 β-sheets containing 8, 7 and 5 strands, respectively. The clustering code also provides information about the orientation of the strands in a β-sheet (i.e., parallel or anti-parallel), alignment of the β-strands within a β-sheet (i.e., in register or out-of register) and nature of the β-sheets (i.e. fully parallel, full anti-parallel or mixed orientations within a sheet). In addition to the clustering analysis, a PTWHAM analysis <xref rid="pcbi.1002051-Yang1" ref-type="bibr">[64]</xref> was also performed on all of our data to compute thermodynamical properties.</p><p>In all cases, structures for all-atom simulations were taken among those of lower-energy OPEP that resisted most efficiently to a temperature increase during replica exchanges. For one preliminary simulation (20OPp) however, the structures were selected based on their frequency of occurrence.</p></sec><sec id="s2e"><title>All-Atom MD Analysis of the Conformational and Stability Properties of Opep-Generated, Selected Oligomeric Structures</title><p>The initial structures for all-atom, explicit solvent Molecular Dynamics (MD) simulations were built by reconstructing the atomic detail of selected conformations from the OPEP coarse-grained runs. Reconstruction was carried out using the MAXSPROUT server <xref rid="pcbi.1002051-Holm1" ref-type="bibr">[65]</xref>. Refinement of side-chain rotameric states was performed using the program IRECS <xref rid="pcbi.1002051-Hartmann1" ref-type="bibr">[66]</xref>, <xref rid="pcbi.1002051-Dunbrack1" ref-type="bibr">[67]</xref>, where the prediction is guided by a combination of potential interaction and rotamer scores calculated with probabilities from the backbone dependent rotamer library. Resulting all-atom structures obtained with this procedure were first minimized using the Macromodel package (Schrodinger Incorporated, USA) for 5000 steps with Polak-Ribier Conjugate Gradient method and an energy gradient criterion for convergence set to 0.05 kJ/mol. This minimization protocol was intended to initially remove unphysical contacts between atoms resulting from the reconstruction procedure, and not to optimize structures. At this stage, the Cα atoms were constrained to their positions with the default force constant (25 kcal/mol Å<sup>2</sup>).</p><p>The resulting minimized systems were then solvated in a cubic-shaped box large enough to contain 1.0<italic>nm</italic> of solvent around each initial aggregate. The simple point charge (SPC) water model was used <xref rid="pcbi.1002051-Berendsen2" ref-type="bibr">[68]</xref> to solvate each oligomer in the simulation box. Each system was subsequently energy minimized with a steepest descent method for 5000 steps. The minimization was considered to be converged when the maximum force was smaller than 0.0001 kJ mol<sup>−1</sup> nm<sup>−1</sup>. The initial step size for minimization was 0.01 nm. The calculation of electrostatic forces was done with the PME implementation of the Ewald summation method. The LINCS <xref rid="pcbi.1002051-Hess2" ref-type="bibr">[69]</xref> algorithm was used to constrain all bond lengths and the SETTLE algorithm <xref rid="pcbi.1002051-Miyamoto1" ref-type="bibr">[70]</xref> for the water molecules. Simulations were performed with a dielectric permittivity, <inline-formula><inline-graphic xlink:href="pcbi.1002051.e002.jpg" mimetype="image"></inline-graphic></inline-formula> = 1, and a time step of 2 fs. Initial velocities were taken from a Maxwellian distribution at the desired initial temperature of 300 K. The density of the system was adjusted performing the first equilibration runs at NPT condition by weak coupling to a bath of constant pressure (P<sub>0</sub> = 1 bar, coupling time τ<italic><sub>P</sub></italic> = 0.5 ps) <xref rid="pcbi.1002051-Berendsen1" ref-type="bibr">[60]</xref> and the system was weakly coupled to an external temperature bath <xref rid="pcbi.1002051-Berendsen1" ref-type="bibr">[60]</xref> with a coupling constant of 0.1 ps. The proteins and the rest of the system were coupled separately to the temperature bath. <xref ref-type="table" rid="pcbi-1002051-t001">Table 1</xref> summarizes the simulation conditions and number of peptides for each simulation. All simulations and analysis were carried out using the GROMACS package (version 3.3) <xref rid="pcbi.1002051-Berendsen3" ref-type="bibr">[71]</xref>–<xref rid="pcbi.1002051-vanderSpoel1" ref-type="bibr">[73]</xref> and the GROMOS96 43A1 force field <xref rid="pcbi.1002051-vanGunsteren1" ref-type="bibr">[74]</xref>–<xref rid="pcbi.1002051-vanGunsteren3" ref-type="bibr">[77]</xref>.</p><p>For all MD simulations, aggregates were simulated at 300 K for 100 ns. REMD simulations were also used to investigate the stability and the conformational preferences of two 20-mer aggregates. The replica exchange simulations were carried out using the Solute Tempering REMD <xref rid="pcbi.1002051-Liu1" ref-type="bibr">[78]</xref> protocol using the version implemented in GROMACS by de Groot and coworkers <xref rid="pcbi.1002051-Kubitzki1" ref-type="bibr">[79]</xref>. Twelve temperatures between 308 K and 419 K were selected according to <xref rid="pcbi.1002051-vanderSpoel2" ref-type="bibr">[80]</xref> for an exchange probability of around 40%.</p></sec></sec><sec id="s3"><title>Results/Discussion</title><p>The aggregation process for the three types of GNNQQNY oligomers – containing 3, 12 and 20 chains, respectively – was studied by a multi-scale approach consisting in a preliminary, thorough exploration of the phase space through REMD with the OPEP coarse-grained potential, followed by the refinement of the most representative aggregate structures obtained via all-atom MD or REMD simulations in explicit solvent. The initial concentration for the OPEP runs was around 4.15 mM. This concentration is 10 times higher than the concentration at which amyloid GNNQQNY fibrils form in a few hours according to Nelson et al. <xref rid="pcbi.1002051-Thirumalai1" ref-type="bibr">[9]</xref> allowing for the formation of ordered structures within our simulation time frame. The diversity in the number of chains allows us to examine possible intermediates and analyze molecular mechanisms of polymorphism in amyloid aggregates.</p><p>For clarity, we first present and discuss results for the trimeric and dodecameric systems as they will serve as basis for understanding the results observed for the 20-mer presented in the last part of this section.</p><sec id="s3a"><title>Simulations of Trimeric Systems</title><sec id="s3a1"><title>Coarse-grained simulations</title><p>Coarse-Grained REMD simulations were performed with 16 replicas for 50 ns at temperatures discussed in the materials and methods section. Although the system is not fully converged for the very low-temperature replicas, the PTWHAM-generated specific heat computed over two different time intervals shows that the melting temperature, T<sub>m</sub>, is well-established at ∼294 K (<xref ref-type="fig" rid="pcbi-1002051-g003"><bold>Figure 3</bold></xref>). Below this temperature, a clustering analysis shows that GNNQQNY monomers are assembled into ordered structures with high β-sheet content, while above T<sub>m</sub>, the system visits mostly disordered structures with very low secondary structure composition. The alignment of individual strands within oligomers, the secondary structures and the configuration types of the aggregates are summarized in <xref ref-type="table" rid="pcbi-1002051-t002"><bold>Table 2</bold></xref>.</p><fig id="pcbi-1002051-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g003</object-id><label>Figure 3</label><caption><title>Specific heat as a function of temperature for the trimer and dodecamer systems.</title><p>The specific heat is calculated over two time intervals for each system (trimer on the left panel and dodecamer on the right panel). Both systems have converged over the time windows displayed here.</p></caption><graphic xlink:href="pcbi.1002051.g003"></graphic></fig><table-wrap id="pcbi-1002051-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.t002</object-id><label>Table 2</label><caption><title>Structural characteristics for small aggregates as a function of temperature.</title></caption><alternatives><graphic id="pcbi-1002051-t002-2" xlink:href="pcbi.1002051.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td colspan="7" align="left" rowspan="1">Population</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">222.5 K</td><td align="left" rowspan="1" colspan="1">235.8 K</td><td align="left" rowspan="1" colspan="1">250.8 K</td><td align="left" rowspan="1" colspan="1">266.7 K</td><td align="left" rowspan="1" colspan="1">283.4 K</td><td align="left" rowspan="1" colspan="1">300 K</td><td align="left" rowspan="1" colspan="1">313.8 K</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Trimer</td><td align="left" rowspan="1" colspan="1">Configuration types (%)<xref ref-type="table-fn" rid="nt109">(a)</xref></td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">96.3</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">98.2</td><td align="left" rowspan="1" colspan="1">98.2</td><td align="left" rowspan="1" colspan="1">92.6</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">1.9</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">2 1</td><td align="left" rowspan="1" colspan="1">3.7</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1.8</td><td align="left" rowspan="1" colspan="1">1.8</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">22.2</td><td align="left" rowspan="1" colspan="1">9.2</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">1 1 1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">5.5</td><td align="left" rowspan="1" colspan="1">64.8</td><td align="left" rowspan="1" colspan="1">88.9</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% parallel<xref ref-type="table-fn" rid="nt110">(b)</xref></td><td align="left" rowspan="1" colspan="1">45.4</td><td align="left" rowspan="1" colspan="1">39.8</td><td align="left" rowspan="1" colspan="1">22.2</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">6.5</td><td align="left" rowspan="1" colspan="1">7.4</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% antiparallel<xref ref-type="table-fn" rid="nt110">(b)</xref></td><td align="left" rowspan="1" colspan="1">54.6</td><td align="left" rowspan="1" colspan="1">60.2</td><td align="left" rowspan="1" colspan="1">77.8</td><td align="left" rowspan="1" colspan="1">88</td><td align="left" rowspan="1" colspan="1">88</td><td align="left" rowspan="1" colspan="1">24.1</td><td align="left" rowspan="1" colspan="1">5.6</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% fully parallel sheets<xref ref-type="table-fn" rid="nt111">(c)</xref></td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">21.1</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% fully antiparallel sheets<xref ref-type="table-fn" rid="nt111">(c)</xref></td><td align="left" rowspan="1" colspan="1">9.3</td><td align="left" rowspan="1" colspan="1">20.4</td><td align="left" rowspan="1" colspan="1">57.4</td><td align="left" rowspan="1" colspan="1">75.9</td><td align="left" rowspan="1" colspan="1">86.3</td><td align="left" rowspan="1" colspan="1">68.4</td><td align="left" rowspan="1" colspan="1">50</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% mixed sheets<xref ref-type="table-fn" rid="nt111">(c)</xref></td><td align="left" rowspan="1" colspan="1">90.7</td><td align="left" rowspan="1" colspan="1">79.6</td><td align="left" rowspan="1" colspan="1">40.7</td><td align="left" rowspan="1" colspan="1">24.1</td><td align="left" rowspan="1" colspan="1">13.7</td><td align="left" rowspan="1" colspan="1">10.5</td><td align="left" rowspan="1" colspan="1">50</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">β-sheet content (%)<xref ref-type="table-fn" rid="nt112">(d)</xref></td><td align="left" rowspan="1" colspan="1">72.4</td><td align="left" rowspan="1" colspan="1">70.6</td><td align="left" rowspan="1" colspan="1">66.6</td><td align="left" rowspan="1" colspan="1">61.6</td><td align="left" rowspan="1" colspan="1">55.1</td><td align="left" rowspan="1" colspan="1">33.3</td><td align="left" rowspan="1" colspan="1">28.3</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% Strands in-register/out-of-register by 1 residue<xref ref-type="table-fn" rid="nt113">(e)</xref></td><td align="left" rowspan="1" colspan="1">51.9/46.3</td><td align="left" rowspan="1" colspan="1">47.2/48.2</td><td align="left" rowspan="1" colspan="1">47.2/38.9</td><td align="left" rowspan="1" colspan="1">38.0/43.5</td><td align="left" rowspan="1" colspan="1">17.7/45.1</td><td align="left" rowspan="1" colspan="1">21.1/42.1</td><td align="left" rowspan="1" colspan="1">33.3/41.7</td></tr><tr><td align="left" rowspan="1" colspan="1">Dodecamer</td><td align="left" rowspan="1" colspan="1">Configuration types (%)<xref ref-type="table-fn" rid="nt109">(a)</xref></td><td align="left" rowspan="1" colspan="1">7 5</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">9.8</td><td align="left" rowspan="1" colspan="1">54.1</td><td align="left" rowspan="1" colspan="1">75.2</td><td align="left" rowspan="1" colspan="1">34.6</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">7 4 1</td><td align="left" rowspan="1" colspan="1">30.8</td><td align="left" rowspan="1" colspan="1">33.1</td><td align="left" rowspan="1" colspan="1">12.8</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">8 4</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">11.3</td><td align="left" rowspan="1" colspan="1">38.4</td><td align="left" rowspan="1" colspan="1">36.8</td><td align="left" rowspan="1" colspan="1">5.3</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">66.9</td><td align="left" rowspan="1" colspan="1">7.5</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">11 1</td><td align="left" rowspan="1" colspan="1">2.3</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">10 1 1</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% parallel<xref ref-type="table-fn" rid="nt110">(b)</xref></td><td align="left" rowspan="1" colspan="1">32.3</td><td align="left" rowspan="1" colspan="1">31.3</td><td align="left" rowspan="1" colspan="1">40.9</td><td align="left" rowspan="1" colspan="1">47.9</td><td align="left" rowspan="1" colspan="1">41.9</td><td align="left" rowspan="1" colspan="1">37.9</td><td align="left" rowspan="1" colspan="1">17.1</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% antiparallel<xref ref-type="table-fn" rid="nt110">(b)</xref></td><td align="left" rowspan="1" colspan="1">60.2</td><td align="left" rowspan="1" colspan="1">48.4</td><td align="left" rowspan="1" colspan="1">45.8</td><td align="left" rowspan="1" colspan="1">46.7</td><td align="left" rowspan="1" colspan="1">49.5</td><td align="left" rowspan="1" colspan="1">41.1</td><td align="left" rowspan="1" colspan="1">16</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% fully parallel sheets<xref ref-type="table-fn" rid="nt111">(c)</xref></td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">1.2</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">15.3</td><td align="left" rowspan="1" colspan="1">22.6</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% fully antiparallel sheets<xref ref-type="table-fn" rid="nt111">(c)</xref></td><td align="left" rowspan="1" colspan="1">6.3</td><td align="left" rowspan="1" colspan="1">10.8</td><td align="left" rowspan="1" colspan="1">6.2</td><td align="left" rowspan="1" colspan="1">2.3</td><td align="left" rowspan="1" colspan="1">2.7</td><td align="left" rowspan="1" colspan="1">11.4</td><td align="left" rowspan="1" colspan="1">21.8</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% mixed sheets<xref ref-type="table-fn" rid="nt111">(c)</xref></td><td align="left" rowspan="1" colspan="1">93.7</td><td align="left" rowspan="1" colspan="1">88.8</td><td align="left" rowspan="1" colspan="1">92.6</td><td align="left" rowspan="1" colspan="1">97.3</td><td align="left" rowspan="1" colspan="1">96.5</td><td align="left" rowspan="1" colspan="1">73.3</td><td align="left" rowspan="1" colspan="1">55.6</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">β-sheet content (%)<xref ref-type="table-fn" rid="nt112">(d)</xref></td><td align="left" rowspan="1" colspan="1">75.8</td><td align="left" rowspan="1" colspan="1">59.1</td><td align="left" rowspan="1" colspan="1">61.7</td><td align="left" rowspan="1" colspan="1">63.5</td><td align="left" rowspan="1" colspan="1">61.8</td><td align="left" rowspan="1" colspan="1">43.8</td><td align="left" rowspan="1" colspan="1">11.4</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% Strands in-register/out-of-register by 1 residue<xref ref-type="table-fn" rid="nt113">(e)</xref></td><td align="left" rowspan="1" colspan="1">52.5/27.1</td><td align="left" rowspan="1" colspan="1">33.5/46.3</td><td align="left" rowspan="1" colspan="1">39.2/41.7</td><td align="left" rowspan="1" colspan="1">44.7/43.4</td><td align="left" rowspan="1" colspan="1">38.8/54.0</td><td align="left" rowspan="1" colspan="1">30.6/45.6</td><td align="left" rowspan="1" colspan="1">29.5/39.7</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt108"><label></label><p>Temperatures above 313.8 K are not displayed here since they are populated essentially by conformations with random coil monomers with no secondary structure. The percentages are calculated over all the structures obtained in the last 40 ns (trimer) and in the last 100 ns (dodecamer) of the OPEP REMD simulations, where the systems have converged.</p></fn><fn id="nt109"><label>(a)</label><p>The dominant configuration types (as described in the OPEP Analysis and Structure Selection section).</p></fn><fn id="nt110"><label>(b)</label><p>The average amount of parallel and anti-parallel strands in the β-sheets formed. The sum of parallel and antiparallel strands in a structure does not always total 100% if the structure sees strands in an undefined orientation, i.e. attached by only one hydrogen bond.</p></fn><fn id="nt111"><label>(c)</label><p>The average amount of fully parallel, fully antiparallel and mixed sheets.</p></fn><fn id="nt112"><label>(d)</label><p>The average amount of residues in a β conformation.</p></fn><fn id="nt113"><label>(e)</label><p>The average amount of strands in-register and out-of-register (by one residue) in β-sheets.</p></fn></table-wrap-foot></table-wrap><p>Structurally, the trimer displays a strong tendency to form ordered planar β-sheets below T<sub>m</sub> (<xref ref-type="fig" rid="pcbi-1002051-g004"><bold>Figure 4</bold></xref><bold>, left part of the panel</bold>). These appear rapidly, within 1 to 8 ns, in a mostly anti-parallel organization. Following trajectories leading to ordered structures, we see that the three-stranded β-sheet is always preceded by the formation of a mostly anti-parallel dimer seed. Averaging over all structures below T<sub>m</sub>, we find that only a very small proportion of structures just below T<sub>m</sub> consist of a two-stranded β-sheet interacting with one chain in coil conformation (1.9%) or three random coil chains (1.1%). The peptides at a temperature just below T<sub>m</sub> prefer an anti-parallel β-strand order (87% at 267 K) over a parallel arrangement (13% at 267 K), while this proportion falls to 55–60% at the lowest temperatures. As seen in <xref ref-type="table" rid="pcbi-1002051-t002"><bold>Table 2</bold></xref>, the three β-strands prefer to be perfectly aligned or in-registered at the lowest temperatures and are typically shifted by one residue, i.e. out-of-registered, at temperatures close to the melting point. As the temperature increases, the population of two-stranded and three-stranded β-sheets becomes very low, amounting to 8% and 0% at 333 K and 352 K, respectively. Except for the lowest temperatures, where mixed parallel/antiparallel sheets are most common, there is a clear dominance of fully antiparallel sheets for three-stranded structures while fully parallel sheets are rare, even among the few three-stranded sheets found above T<sub>m</sub>, where they reach 21%, to 68% for fully anti-parallel.</p><fig id="pcbi-1002051-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g004</object-id><label>Figure 4</label><caption><title>Structures obtained for the trimeric simulations.</title><p>We show, on the left-hand side panel, representative structures obtained from the OPEP simulations and, on the right-hand side panel, the representative structures obtained after all-atom MD refinements. 03OP1-A,-B,-C,-D and –E were extracted respectively at 222.5 K (probability of occurrence for this β-strand organization: 91%), 235.7 K (80%), 250.8 K (41%), 266.7 K (76%) and 283.4 K (86%). 03OP1-A to -C are mixed β-sheets while 03OP1-D and –E are fully antiparallel β-sheets. The all-atom structures are represented in secondary structure cartoon and only the tyrosines (most hydrophobic residues in the sequence) are shown in blue sticks (hydrogen atoms are omitted).</p></caption><graphic xlink:href="pcbi.1002051.g004"></graphic></fig></sec><sec id="s3a2"><title>All-atom MD simulations</title><p>Five representative OPEP-generated structures, labeled 03OP1-A, 03OP1-B, 03OP1-C, 03OP1-D and 03OP1-E (<xref ref-type="fig" rid="pcbi-1002051-g004"><bold>Figure 4</bold></xref><bold>, left side of the panel</bold>), were then subjected to all-atom MD simulations as described in materials and methods. These structures can be divided in two sets: 03OP1-A, 03OP1-B, 03OP1-C are characterized by three-stranded β-sheets with mixed parallel/anti-parallel β-strands, while 03OP1-D and 03OP1-E display a fully anti-parallel three-stranded β-sheet.</p><p>As seen in the final structures of the all-atom simulations displayed in <xref ref-type="fig" rid="pcbi-1002051-g004">Figure 4</xref> (right side of the panel), the five structures show different evolutions after the 100 ns all-atom MD. The three structures 03OP1-A, -B and -D tend towards configuration types 2-1, i.e. with one chain converted from β-strand to random coil and the two other chains enhancing their β-sheet contents. This inter-conversion is independent on the initial orientation of the strands. In contrast, the other two structures 03OP1-C and –E preserve their three-stranded β-sheet configurations and enhance their β-sheet contents. Simulation 03GR1-C keeps its starting mixed parallel/anti-parallel configuration of the strands; in the simulation 03GR1-E, one of the peptide flips orientation leading to a perfectly aligned mixed β-sheet from an initial fully anti-parallel sheet.</p><p>Even though all-atom simulations cannot capture fully disordered chains within 100 ns at 300 K, the coarse-grained and all-atom simulations indicate that both parallel and anti-parallel arrangements can be found in multiple meta-stable two-stranded and three-stranded structures, with various registers of hydrogen bonds contributing to the structural richness and conformational variability of the trimeric aggregates.</p><p>Our trimeric results point to the existence of three minima associated with parallel, antiparallel and mixed parallel/antiparallel β-sheet structures, and are consistent with previous computational studies at the all-atom level on the GNNQQNY trimer <xref rid="pcbi.1002051-Meli1" ref-type="bibr">[44]</xref>, <xref rid="pcbi.1002051-Zhang1" ref-type="bibr">[45]</xref>, <xref rid="pcbi.1002051-Gsponer1" ref-type="bibr">[51]</xref>. Our conformational distribution for the trimer is not biased, therefore, from the use of the OPEP coarse-grained potential. We emphasize that the population of the fully parallel and antiparallel β-structures in small aggregates vary substantially with the selected force field. Using CHARMM force field and the EEF1 implicit water model, Wales <italic>et al.</italic> predicted equal populations for both states from free energy calculations <xref rid="pcbi.1002051-Strodel1" ref-type="bibr">[46]</xref>. Lai <italic>et al.</italic> using multiple MD simulations with the Gromos force field and the SPC explicit water models found many transitions between both states <xref rid="pcbi.1002051-Zhang1" ref-type="bibr">[45]</xref>, while Reddy <italic>et al.</italic> using the same Gromos force field and the SPC explicit water model predicted a much higher population for the parallel geometry <xref rid="pcbi.1002051-Reddy1" ref-type="bibr">[47]</xref>.</p></sec></sec><sec id="s3b"><title>Simulations of Dodecameric Systems</title><sec id="s3b1"><title>Coarse-grained simulations</title><p>OPEP-REMD was performed with the 16 replicas as in the case of the trimer, but each for 125 ns. Within the first 25 ns, the system converges at low temperature to β-sheet rich structures where the strands prefer an antiparallel orientation, as for the trimer, but with a lower melting temperature of 283 K (see <xref ref-type="fig" rid="pcbi-1002051-g003"><bold>Figure 3</bold></xref>) even though the potential energy per monomer in the ordered phase is much lower, reaching −37.0 kcal/mol/monomer for the 12-mer compared to −18.4 kcal/mol/monomer for the trimer, indicating a clear bias toward aggregation and resulting in a much more marked peak in the specific heat.</p><p>Kinetically, the aggregation tendency for the dodecamer is to first form one or two stable four-stranded β-sheets that show little dissociation and that trigger the transient formation of one or two longer β-sheets. The formation of a trimer that precedes the four-stranded β-sheet shows, however, a higher dissociation/association rate. Interestingly, the tendency of the GNNQQNY sequence to form stable tetrameric aggregation nuclei had already been noticed in a previous investigation on the system <xref rid="pcbi.1002051-Meli1" ref-type="bibr">[44]</xref> and was proposed by the Eisenberg group on the basis of entropic and energetic arguments <xref rid="pcbi.1002051-Nelson1" ref-type="bibr">[7]</xref>. The final stable ordered structures are shown in <xref ref-type="fig" rid="pcbi-1002051-g005"><bold>Figure 5</bold></xref> (left side of the panel).</p><fig id="pcbi-1002051-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g005</object-id><label>Figure 5</label><caption><title>Structures obtained for the dodecameric simulations.</title><p>We show, on the left-hand side panel, representative structures obtained from the OPEP simulations and, on the right-hand side panel, representative structures obtained after all-atom MD refinements. 12OP1-A,-B,-C,-D and –E were extracted respectively at 222.5 K, 235.7 K, 250.8 K, 266.7 K and 283.4 K. 12OP1-A (top left structure) is a long flat beta-sheet. 12OP1-B to -E (second left to bottom left structures) are made of 2 beta-sheets facing each other. Monomers forming β-sheets in the initial state are colored red or green. These colors are kept in the final structure. The tyrosines are shown in blue sticks for the all-atom structures. During the all-atom MD simulation the structures tend to be more globular but the strands see no exchange between the β-sheets, i.e. the red and green β-sheets do not dissociate for the 12-mer system.</p></caption><graphic xlink:href="pcbi.1002051.g005"></graphic></fig><p>As would be expected, a rich set of ordered configurations is visited for the 12-mer (<xref ref-type="table" rid="pcbi-1002051-t002"><bold>Table 2</bold></xref>). Regrouping all structures below melting, the dominant conformation, visited 63% of the time, is a two β-sheet structure with a 7 or 8-strand sheet stabilized by a smaller, 4–5 strand sheet positioned on top (<xref ref-type="fig" rid="pcbi-1002051-g005"><bold>Figure 5</bold></xref><bold>, structures 12OP1-B to -E</bold>). Single sheets, with 11 or 12 strands also appear with a frequency of 23.3% below melting (<xref ref-type="fig" rid="pcbi-1002051-g005"><bold>Figure 5</bold></xref><bold>, structure 12OP1-A</bold>). Surprisingly, strand orientation probabilities vary significantly going from the 3-peptide to the 12-peptide system. As for the 3-peptide system, the anti-parallel orientation is favored below melting for the 12-peptide system especially at the lowest two temperatures where the probability of forming anti-parallel is between 60% and 45% compared to 30% for the parallel. Then, as the temperature is increased, the amount of parallel and anti-parallel orientation becomes almost the same, suggesting that while anti-parallel orientation is energetically favored, it is rapidly overcome by the entropic gain of mixing orientations. The alignment of the β-strands is a mix of perfectly aligned strands and strands misaligned by one residue at all temperatures below the melting point. Because sheets are longer than for the trimer, the 12-mer comprises mostly β-sheets with strands in mixed orientations at low temperatures below T<sub>m</sub> with a low probability of forming fully parallel or fully antiparallel sheets (<xref ref-type="table" rid="pcbi-1002051-t002"><bold>Table 2</bold></xref>). Interestingly in the few and much smaller sheets observed just above T<sub>m</sub>, fully parallel and antiparallel β sheets form with almost identical probability (data not shown), suggesting that with slower growth, structures visited below T<sub>m</sub> could be more ordered.</p></sec><sec id="s3b2"><title>All-atom MD simulations</title><p>The 5 most representative structures obtained from OPEP REMD (labeled 12OP1-A to 12OP1-E) were further studied by all-atom MD. Representative structures obtained from the latter simulations are shown in <xref ref-type="fig" rid="pcbi-1002051-g005">Figure 5</xref> right panel. The 12OP1-A OPEP structure is characterized by the presence of a flat arrangement of β-sheets. It undergoes significant rearrangements during the all-atom evolution in explicit solvent (12GR1-A), as shown by the time evolution of the radius of gyration (<bold><xref ref-type="supplementary-material" rid="pcbi.1002051.s001">Figure S1</xref></bold>), with the planar β-sheet breaking into four fragments of two to four stranded β-sheets that assemble on top of each other, with two central parallel β-sheets covered on both sides by a perpendicular β-sheet. The overall amount of β-sheet structure is conserved during the all-atom simulation (<xref ref-type="table" rid="pcbi-1002051-t003"><bold>Table 3</bold></xref>).</p><table-wrap id="pcbi-1002051-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.t003</object-id><label>Table 3</label><caption><title>Comparison of the structural properties between OPEP and GROMACS structures for selected stable structures.</title></caption><alternatives><graphic id="pcbi-1002051-t003-3" xlink:href="pcbi.1002051.t003"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td colspan="4" align="left" rowspan="1">β-Sheet Contenent (based on the dssp program)<xref ref-type="table-fn" rid="nt116">(b)</xref></td><td colspan="4" align="left" rowspan="1">% Parallel - Antiparallel<xref ref-type="table-fn" rid="nt117">(c)</xref></td><td colspan="4" align="left" rowspan="1">% Res Align (% 100 Align - % Shift by 1 Res - % Shift by 2 Res…)<xref ref-type="table-fn" rid="nt118">(d)</xref></td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="left" rowspan="1">GROMACS - GR</td><td colspan="2" align="left" rowspan="1">OPEP - OP</td><td colspan="2" align="left" rowspan="1">GROMACS - GR</td><td colspan="2" align="left" rowspan="1">OPEP - OP</td><td colspan="2" align="left" rowspan="1">GROMACS - GR</td><td colspan="2" align="left" rowspan="1">OPEP - OP</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Configuration Type<xref ref-type="table-fn" rid="nt115">(a)</xref></td><td align="left" rowspan="1" colspan="1">Structure</td><td align="left" rowspan="1" colspan="1">First Cluster</td><td align="left" rowspan="1" colspan="1">Final Structure</td><td align="left" rowspan="1" colspan="1">CG</td><td align="left" rowspan="1" colspan="1">Min</td><td align="left" rowspan="1" colspan="1">First Cluster</td><td align="left" rowspan="1" colspan="1">Final Structure</td><td align="left" rowspan="1" colspan="1">CG</td><td align="left" rowspan="1" colspan="1">Min</td><td align="left" rowspan="1" colspan="1">First Cluster</td><td align="left" rowspan="1" colspan="1">Final Structure</td><td align="left" rowspan="1" colspan="1">CG</td><td align="left" rowspan="1" colspan="1">Min</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Trimers</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">03_1 - A</td><td align="left" rowspan="1" colspan="1">38%</td><td align="left" rowspan="1" colspan="1">19%</td><td align="left" rowspan="1" colspan="1">71% (100%)</td><td align="left" rowspan="1" colspan="1">67%</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 0</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">03_1 - B</td><td align="left" rowspan="1" colspan="1">38%</td><td align="left" rowspan="1" colspan="1">24%</td><td align="left" rowspan="1" colspan="1">19% (27%)</td><td align="left" rowspan="1" colspan="1">67%</td><td align="left" rowspan="1" colspan="1">0 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 0 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">03_1 - C</td><td align="left" rowspan="1" colspan="1">57%</td><td align="left" rowspan="1" colspan="1">62%</td><td align="left" rowspan="1" colspan="1">62% (87%)</td><td align="left" rowspan="1" colspan="1">67%</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">03_1 - D</td><td align="left" rowspan="1" colspan="1">48%</td><td align="left" rowspan="1" colspan="1">57%</td><td align="left" rowspan="1" colspan="1">52% (73%)</td><td align="left" rowspan="1" colspan="1">67%</td><td align="left" rowspan="1" colspan="1">0 - 50</td><td align="left" rowspan="1" colspan="1">0 - 50</td><td align="left" rowspan="1" colspan="1">0 - 50</td><td align="left" rowspan="1" colspan="1">0 - 100</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">03_1 - E</td><td align="left" rowspan="1" colspan="1">48%</td><td align="left" rowspan="1" colspan="1">52%</td><td align="left" rowspan="1" colspan="1">57% (80%)</td><td align="left" rowspan="1" colspan="1">53%</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">0 - 50</td><td align="left" rowspan="1" colspan="1">0 - 50</td><td align="left" rowspan="1" colspan="1">0 - 100</td><td align="left" rowspan="1" colspan="1">50 - 0 - 50</td><td align="left" rowspan="1" colspan="1">50 - 0 - 50</td><td align="left" rowspan="1" colspan="1">0 - 50 - 50</td><td align="left" rowspan="1" colspan="1">0 - 50 - 50</td></tr><tr><td align="left" rowspan="1" colspan="1">12-mers</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">12_1 - A</td><td align="left" rowspan="1" colspan="1">49%</td><td align="left" rowspan="1" colspan="1">35%</td><td align="left" rowspan="1" colspan="1">64% (90%)</td><td align="left" rowspan="1" colspan="1">60%</td><td align="left" rowspan="1" colspan="1">60 - 20</td><td align="left" rowspan="1" colspan="1">40 - 10</td><td align="left" rowspan="1" colspan="1">36 - 64</td><td align="left" rowspan="1" colspan="1">36 - 64</td><td align="left" rowspan="1" colspan="1">80 - 20</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">64 - 18 - 18</td><td align="left" rowspan="1" colspan="1">64 - 18 - 18</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">6 4 1</td><td align="left" rowspan="1" colspan="1">12_1 - B</td><td align="left" rowspan="1" colspan="1">14%</td><td align="left" rowspan="1" colspan="1">14%</td><td align="left" rowspan="1" colspan="1">45% (63%)</td><td align="left" rowspan="1" colspan="1">44%</td><td align="left" rowspan="1" colspan="1">13 - 25</td><td align="left" rowspan="1" colspan="1">20 - 40</td><td align="left" rowspan="1" colspan="1">22 - 22</td><td align="left" rowspan="1" colspan="1">33 - 33</td><td align="left" rowspan="1" colspan="1">35 - 25 - 37</td><td align="left" rowspan="1" colspan="1">0 - 40 - 60</td><td align="left" rowspan="1" colspan="1">33 - 56 - 11</td><td align="left" rowspan="1" colspan="1">33 - 56 - 11</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">8 4</td><td align="left" rowspan="1" colspan="1">12_1 - C</td><td align="left" rowspan="1" colspan="1">35%</td><td align="left" rowspan="1" colspan="1">29%</td><td align="left" rowspan="1" colspan="1">49% (68%)</td><td align="left" rowspan="1" colspan="1">42%</td><td align="left" rowspan="1" colspan="1">10 - 20</td><td align="left" rowspan="1" colspan="1">40 - 20</td><td align="left" rowspan="1" colspan="1">60 - 40</td><td align="left" rowspan="1" colspan="1">40 - 30</td><td align="left" rowspan="1" colspan="1">40 - 40 - 20</td><td align="left" rowspan="1" colspan="1">50 - 20 - 30</td><td align="left" rowspan="1" colspan="1">60 - 20 - 20</td><td align="left" rowspan="1" colspan="1">60 - 20 - 20</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">7 5</td><td align="left" rowspan="1" colspan="1">12_1 - D</td><td align="left" rowspan="1" colspan="1">33%</td><td align="left" rowspan="1" colspan="1">24%</td><td align="left" rowspan="1" colspan="1">54% (75%)</td><td align="left" rowspan="1" colspan="1">55%</td><td align="left" rowspan="1" colspan="1">11 - 22</td><td align="left" rowspan="1" colspan="1">0 - 30</td><td align="left" rowspan="1" colspan="1">40 - 60</td><td align="left" rowspan="1" colspan="1">10 - 60</td><td align="left" rowspan="1" colspan="1">45 - 44 - 11</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">50 - 50</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">7 5</td><td align="left" rowspan="1" colspan="1">12_1 - E</td><td align="left" rowspan="1" colspan="1">56%</td><td align="left" rowspan="1" colspan="1">51%</td><td align="left" rowspan="1" colspan="1">46%</td><td align="left" rowspan="1" colspan="1">51%</td><td align="left" rowspan="1" colspan="1">60 - 40</td><td align="left" rowspan="1" colspan="1">50 - 20</td><td align="left" rowspan="1" colspan="1">60 - 40</td><td align="left" rowspan="1" colspan="1">60 - 40</td><td align="left" rowspan="1" colspan="1">60 - 30 - 10</td><td align="left" rowspan="1" colspan="1">50 - 50</td><td align="left" rowspan="1" colspan="1">60 - 30 - 10</td><td align="left" rowspan="1" colspan="1">50 - 40 -10</td></tr><tr><td align="left" rowspan="1" colspan="1">20-mers</td><td align="left" rowspan="1" colspan="1">11 7 1 1</td><td align="left" rowspan="1" colspan="1">20_p - A1</td><td align="left" rowspan="1" colspan="1">43%</td><td align="left" rowspan="1" colspan="1">38%</td><td align="left" rowspan="1" colspan="1">38% (53%)</td><td align="left" rowspan="1" colspan="1">37%</td><td align="left" rowspan="1" colspan="1">44 - 6</td><td align="left" rowspan="1" colspan="1">56 - 13</td><td align="left" rowspan="1" colspan="1">76 - 24</td><td align="left" rowspan="1" colspan="1">44 - 25</td><td align="left" rowspan="1" colspan="1">56 - 31 - 13</td><td align="left" rowspan="1" colspan="1">63 - 19 - 18</td><td align="left" rowspan="1" colspan="1">34 - 33 - 33</td><td align="left" rowspan="1" colspan="1">44 - 37 - 19</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20_p - A2</td><td align="left" rowspan="1" colspan="1">31%</td><td align="left" rowspan="1" colspan="1">34%</td><td align="left" rowspan="1" colspan="1">38% (53%)</td><td align="left" rowspan="1" colspan="1">37%</td><td align="left" rowspan="1" colspan="1">33 - 0</td><td align="left" rowspan="1" colspan="1">50 - 0</td><td align="left" rowspan="1" colspan="1">76 - 24</td><td align="left" rowspan="1" colspan="1">44 - 25</td><td align="left" rowspan="1" colspan="1">40 - 40 - 20</td><td align="left" rowspan="1" colspan="1">63 - 31 - 6</td><td align="left" rowspan="1" colspan="1">34 - 33 - 33</td><td align="left" rowspan="1" colspan="1">44 - 37 - 19</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">11 7 1 1</td><td align="left" rowspan="1" colspan="1">20_p - B1</td><td align="left" rowspan="1" colspan="1">52%</td><td align="left" rowspan="1" colspan="1">49%</td><td align="left" rowspan="1" colspan="1">36% (50%)</td><td align="left" rowspan="1" colspan="1">66%</td><td align="left" rowspan="1" colspan="1">69 - 19</td><td align="left" rowspan="1" colspan="1">59 - 12</td><td align="left" rowspan="1" colspan="1">71 -24</td><td align="left" rowspan="1" colspan="1">75 - 25</td><td align="left" rowspan="1" colspan="1">31 - 44 - 25</td><td align="left" rowspan="1" colspan="1">35 - 47 -18</td><td align="left" rowspan="1" colspan="1">33 - 39 - 28</td><td align="left" rowspan="1" colspan="1">31 - 44 - 25</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20_p - B2</td><td align="left" rowspan="1" colspan="1">47%</td><td align="left" rowspan="1" colspan="1">38%</td><td align="left" rowspan="1" colspan="1">36% (50%)</td><td align="left" rowspan="1" colspan="1">66%</td><td align="left" rowspan="1" colspan="1">40 - 13</td><td align="left" rowspan="1" colspan="1">38 - 13</td><td align="left" rowspan="1" colspan="1">71 -24</td><td align="left" rowspan="1" colspan="1">75 - 25</td><td align="left" rowspan="1" colspan="1">34 - 53 - 13</td><td align="left" rowspan="1" colspan="1">31 - 56 - 13</td><td align="left" rowspan="1" colspan="1">33 - 39 - 28</td><td align="left" rowspan="1" colspan="1">31 - 44 - 25</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">11 7 1 1</td><td align="left" rowspan="1" colspan="1">20_p - C1</td><td align="left" rowspan="1" colspan="1">39%</td><td align="left" rowspan="1" colspan="1">35%</td><td align="left" rowspan="1" colspan="1">44% (44%)</td><td align="left" rowspan="1" colspan="1">59%</td><td align="left" rowspan="1" colspan="1">44 - 19</td><td align="left" rowspan="1" colspan="1">47 - 20</td><td align="left" rowspan="1" colspan="1">72 - 22</td><td align="left" rowspan="1" colspan="1">69 - 25</td><td align="left" rowspan="1" colspan="1">31 - 50 - 19</td><td align="left" rowspan="1" colspan="1">54 - 33 -13</td><td align="left" rowspan="1" colspan="1">33 - 39 - 28</td><td align="left" rowspan="1" colspan="1">31 - 44 - 25</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20_p - C2</td><td align="left" rowspan="1" colspan="1">59%</td><td align="left" rowspan="1" colspan="1">36%</td><td align="left" rowspan="1" colspan="1">44% (44%)</td><td align="left" rowspan="1" colspan="1">59%</td><td align="left" rowspan="1" colspan="1">69 - 25</td><td align="left" rowspan="1" colspan="1">69 - 8</td><td align="left" rowspan="1" colspan="1">72 - 22</td><td align="left" rowspan="1" colspan="1">69 - 25</td><td align="left" rowspan="1" colspan="1">27 - 47 - 26</td><td align="left" rowspan="1" colspan="1">46 - 54</td><td align="left" rowspan="1" colspan="1">33 - 39 - 28</td><td align="left" rowspan="1" colspan="1">31 - 44 - 25</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">12 7 1</td><td align="left" rowspan="1" colspan="1">20_p - D1</td><td align="left" rowspan="1" colspan="1">44%</td><td align="left" rowspan="1" colspan="1">46%</td><td align="left" rowspan="1" colspan="1">44% (62%)</td><td align="left" rowspan="1" colspan="1">64%</td><td align="left" rowspan="1" colspan="1">24 - 35</td><td align="left" rowspan="1" colspan="1">24 - 35</td><td align="left" rowspan="1" colspan="1">35 - 47</td><td align="left" rowspan="1" colspan="1">41 -47</td><td align="left" rowspan="1" colspan="1">59 - 35 - 6</td><td align="left" rowspan="1" colspan="1">65 - 33 -6</td><td align="left" rowspan="1" colspan="1">59 - 35 - 6</td><td align="left" rowspan="1" colspan="1">59 - 35 - 6</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20_p - D2</td><td align="left" rowspan="1" colspan="1">27%</td><td align="left" rowspan="1" colspan="1">27%</td><td align="left" rowspan="1" colspan="1">44% (62%)</td><td align="left" rowspan="1" colspan="1">64%</td><td align="left" rowspan="1" colspan="1">29 - 24</td><td align="left" rowspan="1" colspan="1">22 - 22</td><td align="left" rowspan="1" colspan="1">35 - 47</td><td align="left" rowspan="1" colspan="1">41 - 47</td><td align="left" rowspan="1" colspan="1">63 - 31 - 6</td><td align="left" rowspan="1" colspan="1">61 - 33 -6</td><td align="left" rowspan="1" colspan="1">59 - 35 - 6</td><td align="left" rowspan="1" colspan="1">59 - 35 - 6</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">9 6 5</td><td align="left" rowspan="1" colspan="1">20_2 - A</td><td align="left" rowspan="1" colspan="1">44%</td><td align="left" rowspan="1" colspan="1">44%</td><td align="left" rowspan="1" colspan="1">34% (47%)</td><td align="left" rowspan="1" colspan="1">28%</td><td align="left" rowspan="1" colspan="1">69 - 13</td><td align="left" rowspan="1" colspan="1">61 - 11</td><td align="left" rowspan="1" colspan="1">72 - 22</td><td align="left" rowspan="1" colspan="1">73 - 13</td><td align="left" rowspan="1" colspan="1">31 - 44 -25</td><td align="left" rowspan="1" colspan="1">39 - 39 - 22</td><td align="left" rowspan="1" colspan="1">56 - 39 - 5</td><td align="left" rowspan="1" colspan="1">20 - 53 - 27</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">8 7 5</td><td align="left" rowspan="1" colspan="1">20_2 - B</td><td align="left" rowspan="1" colspan="1">45%</td><td align="left" rowspan="1" colspan="1">39%</td><td align="left" rowspan="1" colspan="1">39% (55%)</td><td align="left" rowspan="1" colspan="1">71%</td><td align="left" rowspan="1" colspan="1">41 - 41</td><td align="left" rowspan="1" colspan="1">17 - 22</td><td align="left" rowspan="1" colspan="1">71 - 18</td><td align="left" rowspan="1" colspan="1">41 - 59</td><td align="left" rowspan="1" colspan="1">41 - 18 - 41</td><td align="left" rowspan="1" colspan="1">56 - 28 - 16</td><td align="left" rowspan="1" colspan="1">35 - 53 - 12</td><td align="left" rowspan="1" colspan="1">41 - 24 - 35</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">7 7 6</td><td align="left" rowspan="1" colspan="1">20_2 - C</td><td align="left" rowspan="1" colspan="1">36%</td><td align="left" rowspan="1" colspan="1">37%</td><td align="left" rowspan="1" colspan="1">52% (73%)</td><td align="left" rowspan="1" colspan="1">65%</td><td align="left" rowspan="1" colspan="1">44 - 17</td><td align="left" rowspan="1" colspan="1">41 - 35</td><td align="left" rowspan="1" colspan="1">76 - 24</td><td align="left" rowspan="1" colspan="1">59 - 41</td><td align="left" rowspan="1" colspan="1">50 - 39 - 11</td><td align="left" rowspan="1" colspan="1">35 - 47 - 18</td><td align="left" rowspan="1" colspan="1">30 - 41 - 29</td><td align="left" rowspan="1" colspan="1">41 - 47 - 12</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">8 8 4</td><td align="left" rowspan="1" colspan="1">20_2 - E</td><td align="left" rowspan="1" colspan="1">32%</td><td align="left" rowspan="1" colspan="1">33%</td><td align="left" rowspan="1" colspan="1">51% (72%)</td><td align="left" rowspan="1" colspan="1">63%</td><td align="left" rowspan="1" colspan="1">28 - 11</td><td align="left" rowspan="1" colspan="1">42 - 16</td><td align="left" rowspan="1" colspan="1">59 - 41</td><td align="left" rowspan="1" colspan="1">50 - 39</td><td align="left" rowspan="1" colspan="1">44 - 39 - 17</td><td align="left" rowspan="1" colspan="1">42 - 37 - 21</td><td align="left" rowspan="1" colspan="1">54 - 23 -23</td><td align="left" rowspan="1" colspan="1">39 - 33 - 28</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">10 8 2</td><td align="left" rowspan="1" colspan="1">20_2 - N</td><td align="left" rowspan="1" colspan="1">41%</td><td align="left" rowspan="1" colspan="1">39%</td><td align="left" rowspan="1" colspan="1">42% (59%)</td><td align="left" rowspan="1" colspan="1">57%</td><td align="left" rowspan="1" colspan="1">39 - 22</td><td align="left" rowspan="1" colspan="1">41 - 24</td><td align="left" rowspan="1" colspan="1">31 - 56</td><td align="left" rowspan="1" colspan="1">31 - 56</td><td align="left" rowspan="1" colspan="1">44 - 28 - 28</td><td align="left" rowspan="1" colspan="1">29 - 47 - 24</td><td align="left" rowspan="1" colspan="1">44 - 25 -31</td><td align="left" rowspan="1" colspan="1">44 - 25 - 31</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt114"><label></label><p>“First Cluster” means the most representative structure of the GROMACS simulations. “Final structure” is the final conformation obtained at the end of the GROMACS simulations. “CG” is the structure extracted at the end of the OPEP simulations before the reconstruction of the side chains. “Min” indicates the structure resulting from the reconstruction of the side chains after a minimization step.</p></fn><fn id="nt115"><label>(a)</label><p>The configuration type (as described in the OPEP Analysis and Structure Selection section).</p></fn><fn id="nt116"><label>(b)</label><p>The average amount (percentage) of residues in a β conformation. For OPEP, the percentage in brackets has been calculated without taking the Glycines into account.</p></fn><fn id="nt117"><label>(c)</label><p>The average amount (percentage) of parallel and anti-parallel strands in a structure. The sum of parallel and antiparallel strands in a structure does not always total 100% if the structure sees strands in an undefined orientation, i.e. attached by only one hydrogen bond.</p></fn><fn id="nt118"><label>(d)</label><p>The average amount of strands in-register and out-of-register (by one residue).</p></fn></table-wrap-foot></table-wrap><p>Structure 12OP1-B is characterized by a mainly parallel twisted β-sheet, with four strands packed on top. This structure is not stable in the all-atom MD setting, simulation 12GR1-B, and evolves towards a compact globular structure as shown by the evolution of the radius of gyration in time (<bold><xref ref-type="supplementary-material" rid="pcbi.1002051.s001">Figure S1</xref></bold>). Interestingly, the external side of the final aggregate is lined with hydrophilic Asn and Gln side chains that provide favorable contacts with the solvent. No specific order is observed for contacts among these side chains, although some cases of interdigitation as seen in the final steric zipper are noticed. The interior of the final aggregate is lined with Tyr aromatic side chains.</p><p>Such a supramolecular organization of the peptides may be representative of one of the soluble intermediates on the pathway to fibril formation. Solubility is favored by the presence of hydrophilic side chains on the external surface of the aggregate. At the same time, the packing of the interior is not optimal, so that the resulting structure may not be in the most favorable arrangement to ensure lasting stability. Water can also access the interior of the globular aggregate, disrupting inter-strand hydrogen bonds, eventually favoring conformational changes.</p><p>Structures 12OP1-C and 12OP1-D are similar to 12OP1-B: the main difference is that four strand pack with their long axis almost perpendicular to the long axis of the extended β-sheet. The main difference between 12OP1-C and 12OP1-D is that the planes defined by the four strands have different inclinations with respect to the plane of the long extended β-sheet. In the all-atom MD setting — simulations 12GR1-C and 12GR1-D — these structures evolve to less globular, but more compact final arrangements than that observed above, with most of the Tyr side-chains in contact with the solvent (<bold><xref ref-type="supplementary-material" rid="pcbi.1002051.s001">Figure S1</xref></bold>). The exterior of the aggregates is lined with Asn, while the interior is more compact than for 12GR1-A and 12GR1-B and packed with the side-chains of Gln, that form a network of van der Waals and hydrogen bonding contacts.</p><p>Finally, structure 12OP1-E is characterized by two orthogonal twisted β-sheets. The OPEP structure is very stable: it does not undergo significant rearrangement during the all-atom MD, contrary to the previous cases, and the β-sheet content remains constant (<xref ref-type="table" rid="pcbi-1002051-t003"><bold>Table 3</bold></xref>). The oligomer is trapped in this conformation by the extensive contacts packing determined by the Tyr side chains in the two sheets. Moreover, the inter-sheet space is filled by Asn and Gln side chains. However no specific packing into the ordered steric zipper is evident.</p><p><xref ref-type="table" rid="pcbi-1002051-t003"><bold>Table 3</bold></xref> recapitulates the conformational heterogeneity and plasticity of the 12-mer aggregates. As a general case, the presence of explicit solvent tend to condense OPEP-generated structures, at the expense of structured β-sheets and the associated parallel-antiparallel structure, strand alignment and register. It must be kept in mind, though, that MD simulations may be affected by sampling limitations associated with the short runs and the presence of solvent.</p><p>Overall, the combined results indicate that the configurational richness increases from the trimer to the 12-mer and that the critical nucleus has not yet been found. Though, the strands do not see much exchange between sheets as seen in <xref ref-type="fig" rid="pcbi-1002051-g005"><bold>Figure 5</bold></xref>. While ordered 12-mers are energetically much more favorable than the trimers, entropic factors may be considered prevalent, favoring a wide variety of metastable structures. The presence of explicit solvent decreases significantly the stability of elongated β-sheets either by increasing the effective hydrophobic interactions or decreasing entropic gains, favoring rather more compact structures. Different molecular mechanisms may be responsible for the stabilization of different conformations, endowed with different solubility properties. Indeed, we have observed globular-like structures with an external region decorated with hydrophilic groups that may determine the oligomers to be soluble in aqueous solution. In contrast, more ordered structures with higher β-sheet content appear to expose more hydrophobic area to the contact with the solvent. In turn, the latter may recruit more monomers or preformed oligomers that can aggregate by the juxtaposition of hydrophobic surfaces. The observations on the 12-mer systems also underline the enormous structural diversity that characterizes the aggregation of amyloidogenic peptides, which is reflected at the macroscopic level in a high degree of polymorphism.</p></sec></sec><sec id="s3c"><title>Simulations of 20-mer Systems</title><p>Next, we turned to the study of 20-mers in order to assess the importance of the number of chains on the final supra-molecular organization and determine whether new structural motifs can emerge.</p><sec id="s3c1"><title>Coarse-grained simulations</title><p>Three REMD simulations with OPEP were thus generated for the GNNQQNY 20-mer systems: 20OPp, 20OP1 and 20OP2. A preliminary run 20OPp was run to identify the four most common low-energy clusters, from which we extract the central structure for each: 20OPp-A, 20OPp-B, 20OPp-C and 20OPp-D (<xref ref-type="fig" rid="pcbi-1002051-g006"><bold>Figure 6</bold></xref><bold>, left panel</bold>). These were used as starting points for MD simulations with GROMACS. The first three are two-sheet structures while the fourth is a three-sheet configuration. What is particularly interesting here is that we obtain a protofibril-like structure (20OPp-B) among the most dominant clusters after only 200 ns starting from a random coil configuration. Interestingly, the protofibril-like structure is possible but not dominant in this preliminary simulation.</p><fig id="pcbi-1002051-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g006</object-id><label>Figure 6</label><caption><title>Structures obtained for the 20-mer preliminary simulations.</title><p>The stable 20-mer structures obtained from OPEP's preliminary simulation 20OPp are shown on the left-hand side panel. The final primary clusters obtained from the OPEP structures with all-atom MD or all-atom REMD are displayed on the right-hand side panel. 20OPp-A,-B,-C and -D were extracted at 283.4 K. The color code is the same as in <xref ref-type="fig" rid="pcbi-1002051-g005">Figure 5</xref>. 20OPp-A is composed of 2 perpendicular β-sheets. 20OPp-B is a twisted β-sandwich fibril-like structure. 20OPp-C is made of 2 sheets on top of one another. 20OPp-D consists of a folded sheet (green) facing another shorter sheet (red). During the all-atom MD simulation the structures tend to be more globular with the strands seeing some exchange between the β-sheets, i.e. the red and green β-sheets from the OPEP structures dissociate and re-associate during the all-atom MD simulations except for structures 20GRp-D1 and -D2.</p></caption><graphic xlink:href="pcbi.1002051.g006"></graphic></fig><p>Following this preliminary run, we have performed two additional simulations 400 ns-long 20OP1 and 20OP2 (<xref ref-type="fig" rid="pcbi-1002051-g007"><bold>Figure 7</bold></xref>) to attempt to better sample the phase space to determine the degree of preference and the importance of the protofibril-like structure among the morphologies accessible to that sequence for twenty peptides. Even after 400 ns, however, neither simulation is fully converged and the melting temperature is evaluated, from specific heat, to be at 280 K or higher, with ordered structures forming successfully below this temperature: the melting temperature is likely to continue to increase with the simulation length as the average nucleation time for the density used here appears to be around 1 µs based on the fact that slightly more than half the trajectories have not yet visited ordered structures during the 400 ns simulation. In spite of this limitation, we observe significant exchange among the trajectories below melting, suggesting that these achieve some degree of thermodynamic equilibrium.</p><fig id="pcbi-1002051-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g007</object-id><label>Figure 7</label><caption><title>Specific heat as a function of temperature for the two 20-mer simulations sets.</title><p>The specific heat is calculated over two time intervals for the systems 20OP1 (left panel) and 20OP2 (right panel) during the last 200 ns.</p></caption><graphic xlink:href="pcbi.1002051.g007"></graphic></fig><p>As for the 12-mer, aggregation is extremely favorable energetically. The melting temperature for 20OP1 varies between 280.4 K and 289.2 K during the last 200 ns of simulation and the energy of ordered structures at the lowest temperature, 223.8 K, is on average −27.8 kcal/mol/monomer for 20OP1, as calculated from the PTWHAM analysis. For 20OP2, the transition is happening between 260.2 K and 290.5 K and the potential energy of aggregated structures at the lowest temperature, 223.8 K, is on average −28.1 kcal/mol, which is comparable to the energies of aggregated structures for 20OP1. Those energies are about 10 Kcal/mol/monomer above the dodecamer structures' energies at 222.5 K: clearly, the structures generated for the 20-mer are not as ordered as those found for the 12-mer due to the much longer time needed to sample these energetically-favorable conformations, but also because the entropic loss associated with full-ordering is larger for the 20-mer. For both the 20OP1 and the 20OP2 simulation sets, random coil structures dominate at simulations whose temperature is above 280 K.</p><p>Following specific trajectories, as they move through temperatures, it is possible to identify sequences of steps leading to low-energy ordered structures. In the more than 25 such events observed in 20OP1 and 20OP2, the aggregation process is systematically triggered by the formation of a few dimers, trimers and/or tetramers seeds. The conformations obtained from both 20OP1 and 20OP2 are structurally similar in the sense that they are almost always composed of three sheets composed of 5 to 9 strands each either facing each other in a triangle-like or organized in a propeller-like or β-sandwich conformation (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref>). Irrespective of the final shape, the system displays a strong tendency to form β-sheets. The five final ordered structures selected from 20OP2 and shown in <xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref> are representative of all three REMD simulation sets: below melting, the 20-chain system mostly forms three β-sheets, but can also form two-sheet structures. Looking at the statistics collected for 20OP1 and 20OP2 (<xref ref-type="table" rid="pcbi-1002051-t004"><bold>Table 4</bold></xref>), we observe that various three-β-sheet configurations with juxtaposed β-sheets containing 8-7-5 monomers or 9-6-5 monomers are frequent below T<sub>m</sub>. Two-β-sheet systems appear to be less frequent but are populated close to the melting temperature as well. Although, sheet-lengths differ slightly between simulations 20OP1 and 20OP2 for the dominant structures, the overall results are consistent. We also observe a high number of possible β-sandwich morphologies for the three-sheet configuration among the final structures obtained where some of them see their three sheets facing each other in a triangle or twisted-around-each-other arrangements (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref>, structures 20OP2-A, -B, -C and -E). These topologies run from a rather well defined 3-fold symmetry (20OP2-E) to more disordered conformations with little symmetry (20OP2-A, -B and –C). The minimal β-sheet unit contains four strands in 20OP2-E, five strands in 20OP2-A,-B and six strands in 20OP2-C. We note that the structure 20OP2-E is reminiscent of the recently proposed structure of Aβ1–40 fibrils with a three-fold symmetry <xref rid="pcbi.1002051-Paravastu1" ref-type="bibr">[14]</xref>. Interestingly, ordered two-sheet conformations with one sheet slightly longer than the other, such as 11 9 and 12 8 (not shown, but close to 20OP2-N, a 8 10 2 configuration in <xref ref-type="fig" rid="pcbi-1002051-g008">Figure 8</xref>), represent a significant fraction of the accessible states either as a β-sandwich –with an occasional insertion of a short β-sheet – or as two perpendicular sheets.</p><fig id="pcbi-1002051-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g008</object-id><label>Figure 8</label><caption><title>Structures obtained for the 20-mer 20OP2 and 20GR2 simulations.</title><p>The final primary clusters obtained from the OPEP structures with all-atom MD or all-atom REMD are displayed on the right-hand side panel. 20OP2-A,-B,-C,-E and -N were extracted respectively at 260.1 K, 249.2 K, 254.7 K, 265.9 K and 292.2 K. The different sheets are distinguished by either a green, red or yellow color and the tyrosines are shown in blue sticks for the all-atom structures. Structures 20OP2-A,-B,-C and -E are composed of 3 sheets twisted around each other while structure 20OP2-N is a 2-sheet fibril-like conformation. During the all-atom MD simulation the structures tend to be more globular with the strands seeing some exchange between the β-sheets, i.e. the red, green and yellow β-sheets from the OPEP structures dissociate and re-associate during the all-atom MD simulations except for the fibril-like structures 20GR2-N1 and -N2.</p></caption><graphic xlink:href="pcbi.1002051.g008"></graphic></fig><table-wrap id="pcbi-1002051-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.t004</object-id><label>Table 4</label><caption><title>Structural characteristics for the 20-mer aggregates as a function of temperature below the melting point.</title></caption><alternatives><graphic id="pcbi-1002051-t004-4" xlink:href="pcbi.1002051.t004"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td colspan="8" align="left" rowspan="1">Population</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">223.8 K</td><td align="left" rowspan="1" colspan="1">249.3 K</td><td align="left" rowspan="1" colspan="1">260.1 K</td><td align="left" rowspan="1" colspan="1">265.9 K</td><td align="left" rowspan="1" colspan="1">270.3 K</td><td align="left" rowspan="1" colspan="1">273.8 K</td><td align="left" rowspan="1" colspan="1">277.1 K</td><td align="left" rowspan="1" colspan="1">280.1 K</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">20OP1</td><td align="left" rowspan="1" colspan="1">Configuration types (%)<sup>(a)</sup></td><td align="left" rowspan="1" colspan="1">8 7 5</td><td align="left" rowspan="1" colspan="1">93.2</td><td align="left" rowspan="1" colspan="1">12.4</td><td align="left" rowspan="1" colspan="1">7.5</td><td align="left" rowspan="1" colspan="1">7.1</td><td align="left" rowspan="1" colspan="1">9.8</td><td align="left" rowspan="1" colspan="1">6.8</td><td align="left" rowspan="1" colspan="1">7.1</td><td align="left" rowspan="1" colspan="1">9.4</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">8 7 3 2</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">44.7</td><td align="left" rowspan="1" colspan="1">2.6</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">7 7 6</td><td align="left" rowspan="1" colspan="1">2.6</td><td align="left" rowspan="1" colspan="1">12.4</td><td align="left" rowspan="1" colspan="1">13.9</td><td align="left" rowspan="1" colspan="1">9.8</td><td align="left" rowspan="1" colspan="1">12.8</td><td align="left" rowspan="1" colspan="1">10.5</td><td align="left" rowspan="1" colspan="1">12.4</td><td align="left" rowspan="1" colspan="1">7.9</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">10 6 4</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">7.5</td><td align="left" rowspan="1" colspan="1">8.7</td><td align="left" rowspan="1" colspan="1">5.3</td><td align="left" rowspan="1" colspan="1">11.3</td><td align="left" rowspan="1" colspan="1">5.3</td><td align="left" rowspan="1" colspan="1">3.4</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">11 5 4</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">7.5</td><td align="left" rowspan="1" colspan="1">6.4</td><td align="left" rowspan="1" colspan="1">6.0</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">11 5 2 2</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">10.5</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">14 6</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">5.3</td><td align="left" rowspan="1" colspan="1">6.8</td><td align="left" rowspan="1" colspan="1">7.1</td><td align="left" rowspan="1" colspan="1">7.5</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">4.9</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">15 5</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">7.1</td><td align="left" rowspan="1" colspan="1">3.8</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">6.8</td><td align="left" rowspan="1" colspan="1">8.3</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">10 10</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">8.7</td><td align="left" rowspan="1" colspan="1">6.0</td><td align="left" rowspan="1" colspan="1">5.3</td><td align="left" rowspan="1" colspan="1">7.9</td><td align="left" rowspan="1" colspan="1">2.3</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">10 5 5</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">3.8</td><td align="left" rowspan="1" colspan="1">8.3</td><td align="left" rowspan="1" colspan="1">7.1</td><td align="left" rowspan="1" colspan="1">3.8</td><td align="left" rowspan="1" colspan="1">2.6</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">11 9</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">3.8</td><td align="left" rowspan="1" colspan="1">9.0</td><td align="left" rowspan="1" colspan="1">11.7</td><td align="left" rowspan="1" colspan="1">6.8</td><td align="left" rowspan="1" colspan="1">7.1</td><td align="left" rowspan="1" colspan="1">4.5</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">1.1</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">4.9</td><td align="left" rowspan="1" colspan="1">5.3</td><td align="left" rowspan="1" colspan="1">5.3</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">0.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% parallel<sup>(b)</sup></td><td align="left" rowspan="1" colspan="1">73.9</td><td align="left" rowspan="1" colspan="1">48.8</td><td align="left" rowspan="1" colspan="1">53.4</td><td align="left" rowspan="1" colspan="1">55.1</td><td align="left" rowspan="1" colspan="1">56.2</td><td align="left" rowspan="1" colspan="1">56.8</td><td align="left" rowspan="1" colspan="1">58.4</td><td align="left" rowspan="1" colspan="1">46.8</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% antiparallel<sup>(b)</sup></td><td align="left" rowspan="1" colspan="1">25.1</td><td align="left" rowspan="1" colspan="1">48.2</td><td align="left" rowspan="1" colspan="1">41.3</td><td align="left" rowspan="1" colspan="1">37.6</td><td align="left" rowspan="1" colspan="1">36.7</td><td align="left" rowspan="1" colspan="1">35.6</td><td align="left" rowspan="1" colspan="1">34.5</td><td align="left" rowspan="1" colspan="1">29.4</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% fully parallel sheets<sup>(c)</sup></td><td align="left" rowspan="1" colspan="1">26.2</td><td align="left" rowspan="1" colspan="1">22.4</td><td align="left" rowspan="1" colspan="1">19.7</td><td align="left" rowspan="1" colspan="1">15.1</td><td align="left" rowspan="1" colspan="1">11.8</td><td align="left" rowspan="1" colspan="1">14.5</td><td align="left" rowspan="1" colspan="1">11.1</td><td align="left" rowspan="1" colspan="1">16.1</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% fully antiparallel sheets<sup>(c)</sup></td><td align="left" rowspan="1" colspan="1">0.9</td><td align="left" rowspan="1" colspan="1">7.9</td><td align="left" rowspan="1" colspan="1">3.1</td><td align="left" rowspan="1" colspan="1">4.6</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">2.5</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">6.5</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% mixed sheets<sup>(c)</sup></td><td align="left" rowspan="1" colspan="1">72.9</td><td align="left" rowspan="1" colspan="1">69.7</td><td align="left" rowspan="1" colspan="1">77.2</td><td align="left" rowspan="1" colspan="1">80.3</td><td align="left" rowspan="1" colspan="1">86.3</td><td align="left" rowspan="1" colspan="1">83</td><td align="left" rowspan="1" colspan="1">87.9</td><td align="left" rowspan="1" colspan="1">77.4</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">β-sheet content (%)<sup>(d)</sup></td><td align="left" rowspan="1" colspan="1">59.1</td><td align="left" rowspan="1" colspan="1">61.3</td><td align="left" rowspan="1" colspan="1">60.5</td><td align="left" rowspan="1" colspan="1">60.5</td><td align="left" rowspan="1" colspan="1">60.5</td><td align="left" rowspan="1" colspan="1">59.8</td><td align="left" rowspan="1" colspan="1">57.1</td><td align="left" rowspan="1" colspan="1">45.7</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% Strands in-register/out-of-register by 1 residue<sup>(e)</sup></td><td align="left" rowspan="1" colspan="1">32.1/38.0</td><td align="left" rowspan="1" colspan="1">43.1/31.2</td><td align="left" rowspan="1" colspan="1">42.6/33.8</td><td align="left" rowspan="1" colspan="1">41.2/36.9</td><td align="left" rowspan="1" colspan="1">40.1/38.4</td><td align="left" rowspan="1" colspan="1">38.8/38.6</td><td align="left" rowspan="1" colspan="1">39.7/38.6</td><td align="left" rowspan="1" colspan="1">38.7/39.8</td></tr><tr><td align="left" rowspan="1" colspan="1">20OP2</td><td align="left" rowspan="1" colspan="1">Configuration types (%)<sup>(a)</sup></td><td align="left" rowspan="1" colspan="1">9 6 5</td><td align="left" rowspan="1" colspan="1">29.3</td><td align="left" rowspan="1" colspan="1">29.3</td><td align="left" rowspan="1" colspan="1">24.4</td><td align="left" rowspan="1" colspan="1">24.1</td><td align="left" rowspan="1" colspan="1">18.4</td><td align="left" rowspan="1" colspan="1">15.4</td><td align="left" rowspan="1" colspan="1">8.7</td><td align="left" rowspan="1" colspan="1">8.3</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">9 7 4</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">7.5</td><td align="left" rowspan="1" colspan="1">4.9</td><td align="left" rowspan="1" colspan="1">6.0</td><td align="left" rowspan="1" colspan="1">3.8</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">8 8 4</td><td align="left" rowspan="1" colspan="1">1.1</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">1.1</td><td align="left" rowspan="1" colspan="1">3.8</td><td align="left" rowspan="1" colspan="1">6.0</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">4.1</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">7 7 6</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">6.0</td><td align="left" rowspan="1" colspan="1">8.7</td><td align="left" rowspan="1" colspan="1">12.0</td><td align="left" rowspan="1" colspan="1">12.8</td><td align="left" rowspan="1" colspan="1">16.2</td><td align="left" rowspan="1" colspan="1">15.0</td><td align="left" rowspan="1" colspan="1">12.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">15 5</td><td align="left" rowspan="1" colspan="1">12.4</td><td align="left" rowspan="1" colspan="1">9.4</td><td align="left" rowspan="1" colspan="1">7.5</td><td align="left" rowspan="1" colspan="1">2.3</td><td align="left" rowspan="1" colspan="1">6.8</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">1.5</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">13 7</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">2.3</td><td align="left" rowspan="1" colspan="1">4.9</td><td align="left" rowspan="1" colspan="1">7.1</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">14 6</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">3.8</td><td align="left" rowspan="1" colspan="1">3.0</td><td align="left" rowspan="1" colspan="1">4.9</td><td align="left" rowspan="1" colspan="1">10.9</td><td align="left" rowspan="1" colspan="1">6.4</td><td align="left" rowspan="1" colspan="1">8.7</td><td align="left" rowspan="1" colspan="1">6.4</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">11 9</td><td align="left" rowspan="1" colspan="1">9.4</td><td align="left" rowspan="1" colspan="1">14.7</td><td align="left" rowspan="1" colspan="1">6.8</td><td align="left" rowspan="1" colspan="1">3.8</td><td align="left" rowspan="1" colspan="1">2.3</td><td align="left" rowspan="1" colspan="1">4.5</td><td align="left" rowspan="1" colspan="1">5.6</td><td align="left" rowspan="1" colspan="1">3.8</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">11 5 4</td><td align="left" rowspan="1" colspan="1">16.5</td><td align="left" rowspan="1" colspan="1">12.4</td><td align="left" rowspan="1" colspan="1">3.8</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">8 7 5</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">1.5</td><td align="left" rowspan="1" colspan="1">1.9</td><td align="left" rowspan="1" colspan="1">5.6</td><td align="left" rowspan="1" colspan="1">6.8</td><td align="left" rowspan="1" colspan="1">6.4</td><td align="left" rowspan="1" colspan="1">6.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">10 10</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.4</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.4</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% parallel<sup>(b)</sup></td><td align="left" rowspan="1" colspan="1">53.1</td><td align="left" rowspan="1" colspan="1">56.6</td><td align="left" rowspan="1" colspan="1">57.1</td><td align="left" rowspan="1" colspan="1">56.8</td><td align="left" rowspan="1" colspan="1">56.7</td><td align="left" rowspan="1" colspan="1">55.7</td><td align="left" rowspan="1" colspan="1">56.7</td><td align="left" rowspan="1" colspan="1">55.1</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% antiparallel<sup>(b)</sup></td><td align="left" rowspan="1" colspan="1">38.1</td><td align="left" rowspan="1" colspan="1">35.7</td><td align="left" rowspan="1" colspan="1">36.5</td><td align="left" rowspan="1" colspan="1">37.5</td><td align="left" rowspan="1" colspan="1">37.9</td><td align="left" rowspan="1" colspan="1">39.3</td><td align="left" rowspan="1" colspan="1">38.4</td><td align="left" rowspan="1" colspan="1">38.5</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% fully parallel sheets<sup>(c)</sup></td><td align="left" rowspan="1" colspan="1">7.4</td><td align="left" rowspan="1" colspan="1">16.9</td><td align="left" rowspan="1" colspan="1">15.0</td><td align="left" rowspan="1" colspan="1">14.2</td><td align="left" rowspan="1" colspan="1">14.1</td><td align="left" rowspan="1" colspan="1">14.1</td><td align="left" rowspan="1" colspan="1">18.7</td><td align="left" rowspan="1" colspan="1">17.1</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% fully antiparallel sheets<sup>(c)</sup></td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">4.7</td><td align="left" rowspan="1" colspan="1">3.2</td><td align="left" rowspan="1" colspan="1">4.2</td><td align="left" rowspan="1" colspan="1">3.9</td><td align="left" rowspan="1" colspan="1">4.4</td><td align="left" rowspan="1" colspan="1">5.1</td><td align="left" rowspan="1" colspan="1">8.1</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% mixed sheets<sup>(c)</sup></td><td align="left" rowspan="1" colspan="1">92.0</td><td align="left" rowspan="1" colspan="1">78.5</td><td align="left" rowspan="1" colspan="1">81.8</td><td align="left" rowspan="1" colspan="1">81.6</td><td align="left" rowspan="1" colspan="1">82.0</td><td align="left" rowspan="1" colspan="1">81.5</td><td align="left" rowspan="1" colspan="1">76.2</td><td align="left" rowspan="1" colspan="1">74.8</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">β-sheet content (%)<sup>(d)</sup></td><td align="left" rowspan="1" colspan="1">66.5</td><td align="left" rowspan="1" colspan="1">63.4</td><td align="left" rowspan="1" colspan="1">63.0</td><td align="left" rowspan="1" colspan="1">62.5</td><td align="left" rowspan="1" colspan="1">61.5</td><td align="left" rowspan="1" colspan="1">60.6</td><td align="left" rowspan="1" colspan="1">58.7</td><td align="left" rowspan="1" colspan="1">56.0</td></tr><tr><td align="left" rowspan="1" colspan="1"></td><td colspan="2" align="center" rowspan="1">% Strands in-register/out-of-register by 1 residue<sup>(e)</sup></td><td align="left" rowspan="1" colspan="1">51.0/36.3</td><td align="left" rowspan="1" colspan="1">50.0/34.0</td><td align="left" rowspan="1" colspan="1">48.1/35.1</td><td align="left" rowspan="1" colspan="1">44.6/37.0</td><td align="left" rowspan="1" colspan="1">42.4/37.2</td><td align="left" rowspan="1" colspan="1">42.3/37.2</td><td align="left" rowspan="1" colspan="1">39.9/39.1</td><td align="left" rowspan="1" colspan="1">39.7/39.4</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt119"><label></label><p>Temperatures above 280.1 K are not displayed here since they are populated essentially by conformations with random coil monomers with no secondary structure. The percentages are calculated over all the structures obtained in the last 200 ns of both OPEP REMD simulations. For details on <sup>(a)–(e)</sup>, see <xref ref-type="table" rid="pcbi-1002051-t002">Table 2</xref>.</p></fn></table-wrap-foot></table-wrap><p>In both REMDs, the β-sheets have a high probability of being in a mixed anti-parallel/parallel orientation state due to their length (over 70% below Tm) (<xref ref-type="table" rid="pcbi-1002051-t004"><bold>Table 4</bold></xref>). As for the dodecamer, as T rises, sheets become shorter and we notice a rise of fully parallel sheets at the five temperatures above T<sub>m</sub> (data not shown). Contrary to the 3-mer and 12-mer where the β-strands prefer an anti-parallel orientation below 300 K, the parallel orientation of the strands is preferred in the 20-mer at all temperatures. The extended chains are also dominantly perfectly aligned at the lowest temperatures and the structures see a mix of perfectly aligned strands and strands misaligned by one residue close to the melting point. The dominance of perfectly aligned strands and parallel orientations of the chains is consistent with the experimental observations that GNNQQNY fibrils display parallel β-strands <xref rid="pcbi.1002051-Thirumalai1" ref-type="bibr">[9]</xref>.</p></sec><sec id="s3c2"><title>All-atom MD simulations of dominant clusters generated by OPEP-REMD 20OPp</title><p>Consistent with the protocols described above for the 3-mers and 12-mers, the dynamical properties of 9 selected oligomeric conformations generated by the OPEP simulations were refined by all-atom MD simulations in explicit solvent (see <xref ref-type="table" rid="pcbi-1002051-t001"><bold>Table 1</bold></xref>).</p><p>The first set of all-atom MD simulations was run on the structures selected from 20OPp calculations. OPEP runs identified two main types of 3D organization for the 20mer: extended β-sheet and globular-like structures. The former are characterized by the presence of two parallel sheets, while the latter are characterized by a circular organization of the strands, in a mostly parallel arrangement. The major representatives of the extended β-sheet like structures obtained from OPEP are labeled 20OPp-B and 20OPp-C; the globular structures are recapitulated by 20OPp-A (which shows a compact part packed by a more extended sheet) and 20OPp-D, see <xref ref-type="fig" rid="pcbi-1002051-g006"><bold>Figure 6</bold></xref>. Two sets of 100 ns MD simulations, starting with different initial velocities, were run for each structure.</p><p>A<sc>ggregates with facing</sc><sc>B-sheets</sc>. The structure of 20OPp-B is characterized by the presence of two perpendicular β-sheets. Each β-sheet consists mainly of twisted parallel strands.</p><p>After the first MD run (20GRp-B1) the two sheets are oriented anti-parallel to each other, forming a tight and elongated structure (<xref ref-type="fig" rid="pcbi-1002051-g006"><bold>Figure 6</bold></xref>). The Asn side chains from opposite sheets occupy the inter-sheet space with some intertwining of the amide side chains. Significant packing is also provided by the aromatic rings of Tyrosine belonging to adjacent strands in the same sheet. The 100-ns generated aggregate also displays a significant degree of twisting in the strands that make up the two facing antiparallel β-sheets. This final structure is similar to that observed by others <xref rid="pcbi.1002051-Esposito1" ref-type="bibr">[50]</xref>. It is important to underline that this twisted configuration forms spontaneously during the MD simulation time starting from a less compact structure.</p><p>In the second MD simulation (20GRp-B2), starting from the same initial structure with a different set of velocities, the β-sheet content decreases due a lower degree of packing of the Tyr side chains and interdigitation of the Gln and Asn side-chains (<xref ref-type="fig" rid="pcbi-1002051-g009"><bold>Figure 9</bold></xref>). Packing interactions, however, still appear to be important in stabilizing the compact structure.</p><fig id="pcbi-1002051-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g009</object-id><label>Figure 9</label><caption><title>Structure 20-mer showing interdigitation of the side-chains.</title><p>This structure is representative of the 20GRp-B2 simulation. It shows the inter-sheet space filled by Asn and Gln side chains, which in some case interdigitate.</p></caption><graphic xlink:href="pcbi.1002051.g009"></graphic></fig><p>During both 20GRp-B1 and 20GRp-B2 MD simulations, the strands are dynamically interchanged between the two β-sheets.</p><p>Two elongated sheets may, however, evolve towards very different supramolecular organizations. In 20OPp-C two elongated β-sheets, with mainly parallel β-strands, are in contact through the terminal Tyr aromatic chains in an extended and non-compact structure (<xref ref-type="fig" rid="pcbi-1002051-g006"><bold>Figure 6</bold></xref>). In both all-atom MD simulations, 20GRp-C1 and 20GRp-C2, the two β-sheets break up and reorganize in more compact conformations. Though, for 20GRp-C2, the strands do not interchange between sheets while 20GRp-C1 sees some mixing between the strands of the 2 β-sheets. The compact final structures still show some of the features we have described in the previous case: anti-parallel orientations of facing sheets, compaction of the side chains. One part of the long parallel β-sheet, represented by a four-stranded unit, detaches from the initial complex and packs onto the compact structure described above, minimizing the exposed surface area. The resulting structures correspond to a β-sandwich composed of three sheets. The time evolution of the radius of gyration is reported in (<bold><xref ref-type="supplementary-material" rid="pcbi.1002051.s002">Figures S2</xref>, <xref ref-type="supplementary-material" rid="pcbi.1002051.s003">S3</xref></bold>).</p><p>E<sc>volution of globular structures</sc>. All-atom MD simulations were also started from more densely packed structures obtained from OPEP runs, represented by 20OPp-D and 20OPp-A. Each starting structure was subjected to two 100 ns MD simulations, resulting in simulations 20GRp-A1, 20GRp-A2, 20GRp-D1 and 20GRp-D2 (<xref ref-type="fig" rid="pcbi-1002051-g006"><bold>Figure 6</bold></xref>). In the case of 20OPp-D both all-atom MD runs show that the starting conformation retains its overall structure with mostly parallel β-strands. Most of the Tyr side chains form a compact hydrophobic core at the interior of the structure. During all-atom MD evolution they reorganize and form small compact clusters that are optimally packed through pi-stacking face-to-face interactions, minimizing unfavorable contacts with the solvent. In general, the interior of the aggregate is compact and hydrophobic while the external surface is decorated with the more hydrophilic amide side chains. From the biochemical point of view, this structure could be described as a “soluble” oligomeric state reminiscent of the ones observed for the 12-mer run.</p><p>Starting from 20OPp-A, simulation 20GRp-A1 evolves towards a compact globular structure, in which parts of the ordered β-sheets are lost and strands are dynamically interchanged between sheets. The aromatic Tyr side chains cluster in the hydrophobic core of the structure and Asn side chains align on the surface. Most of the conserved sheets are still in parallel orientation.</p><p>Strikingly, in simulation 2 (20GRp-A2) the initial structure evolves to form two twisted antiparallel sheets in which the constitutive strands are parallel to each other. This structure resembles the twisted conformation observed for 20OPp-B and was observed in previous simulations <xref rid="pcbi.1002051-Esposito1" ref-type="bibr">[50]</xref>. The inter-sheet space is filled by Asn and Gln side chains, which in some case interdigitate (<xref ref-type="fig" rid="pcbi-1002051-g009"><bold>Figure 9</bold></xref>). The conformation of this intermediate still shows some Tyr aromatic side chains in the interior, disrupting the optimal interdigitation of amidic side chains and causing the structure to be non-perfectly ordered (<xref ref-type="fig" rid="pcbi-1002051-g006"><bold>Figure 6</bold></xref>).</p><p>These results suggest that the sheet organization in the twisted antiparallel conformation(s) may be accessible on the aggregation pathway, once two sheets are formed and docked upon each other. Interestingly, we have observed the formation of elongated, twisted antiparallel structures in MD only in the 20-mer system. The latter appear to evolve preferentially towards globular structures, suggesting that elongated, fibril-like conformations of the oligomers may be accessible only in the presence of a higher number of monomers. At the atomic level, sheet-locking is favored by the packing of Asn and Tyr side-chains. The Tyr aromatic packing and the initial formation of steric-zipper-like structures also provide important contributions in determining the ordering and stabilization of the growing aggregate and, possibly, its evolution to a stable fibril.</p></sec><sec id="s3c3"><title>All-atom MD simulations of the dominant clusters generated by OPEP-REMD 20OP2</title><p>The detailed role of side-chains in determining the conformational characteristics of compact aggregate structures was further evaluated by analyzing at atomic resolution a set of diverse OPEP structures: 20OP2-A, 20OP2-B, 20OP2-C, 20OP2-E and 20OP2-N (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref>).</p><p>The starting structures of the MD simulations from 20OP2-A, 20OP2-C and 20OP2-E of the aggregates all consist of three extended β-sheets, organized in different tertiary arrangements (see <xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref>). The all-atom evolution of structure 20OP2-A leads to a more compact and globular-like conformation with a mixing of the β-strands and a partial loss of ordered β-structure involving external strands (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref><bold>, structure 20GR2-A</bold>). Tyr side-chains mainly pack in the interior of the globule in the representative structures of MD simulations while most of the hydrophilic side chains (Asn and Gln) are located on the exterior of the oligomer pointing towards the solvent. This arrangement would confer water solubility to the oligomer by sequestering hydrophobic-aromatic side-chains to the interior of the aggregate and decorating the external surface with hydrophilic groups.</p><p>In the case of 20OP2-C, the starting structure constituted by three β-sheets, which are lined and twisted along a common axis, is not stable in the all-atom MD setting, and immediately evolves to a more compact globular structure that however does not display specific supramolecular properties or preferential orientations of the strands within the aggregate (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref><bold>, structure 20GR2-C</bold>).</p><p>In the case of 20OP2-E, the evolution of all-atom MD simulations at 300 K determines a large decrease in the degree of ordered β-structure leading to the formation of disordered, amorphous conformations (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref><bold>, structure 20GR2-E</bold>).</p><p>The remaining two representative clusters obtained from OPEP simulations display different three-dimensional arrangements. In the case of 20OP2-B, the structure is characterized by parallel β-sheet motifs that form a less compact conformation than the one observed above. All-atom MD evolution leads to a globular structure with a global reorganization of the β-strands (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref><bold>, structure 20GR2-B</bold>). A large number of Asn and Gln amidic side-chains point towards the interior of the globular structure, and a number of hydrophobic Tyr side-chains are aligned in contact with the solvent. The reorganization of the starting structure to this globular intermediate determines a partial loss of ordered secondary structure in some of the constituent peptides. This structure retains a large hydrophobic surface in contact with the solvent. The large hydrophobic area exposed to water may be one of the causes of the insoluble character of these intermediates.</p><p>Finally, we simulated the structure of cluster 20OP2-N at all-atom resolution as this aggregate forms an elongated structure with two facing β-sheets. MD evolution at 300 K for this system shows no reorganization of the β-strands (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref>, structure 20GR2-N1). In contrast to what we observed for 20OP2-B, the Asn and Gln amidic side-chains fill the space between the two sheets, establishing H-bonding interactions, and showing interdigitation of side-chains reminiscent of the experimentally observed dry steric zipper. The sheets are not all perfectly parallel, and this might oppose the formation of a perfectly packed steric zipper.</p><p>Summarizing, as shown in <xref ref-type="table" rid="pcbi-1002051-t003"><bold>Table 3</bold></xref>, the structural properties of the aggregates cannot be easily rationalized into specific classes: great variability is observed in terms of β-sheet content, percentage of parallel vs. antiparallel arrangements, register and relative orientations of the strands. Interestingly, it appears that the fibril-like structure 20GR2-N1 is stable during the 100-ns all-atom MD simulation, while the other arrangements undergo large reorganization. In general, the consistency of the results in terms of conformational plasticity for all constructs clearly indicates that a wide range of different structures is accessible to relatively large sized oligomers.</p></sec><sec id="s3c4"><title>All-atom REMD simulations of two elongated 20mer structures generated by OPEP-REMD</title><p>The all-atom simulations of 100 ns starting from 20OPp-B and 20OP2-N showed the possibility for the aggregates to remain elongated partially ordered oligomers whose structures are reminiscent of the arrangements observed by X-rays of micro-crystals. In order to gain more insights into the stability and conformational evolution properties of these structures, we set out to run all-atom REMD simulations starting from OPEP structures 20OPp-B (<xref ref-type="fig" rid="pcbi-1002051-g006"><bold>Figure 6</bold></xref>) and 20OP2-N (<xref ref-type="fig" rid="pcbi-1002051-g008"><bold>Figure 8</bold></xref>). The resulting all-atom REMD simulations are labeled 20GRp-B3 and 20GR2-N2, respectively.</p><p>In the REMD simulation labeled 20GRp-B3, and similarly to what is seen in OPEP REMDs, we observe that structures interconvert between compact and elongated conformations with a pair of sheets facing each other. The main representative structures for simulation 20GRp-B3, and their relative stabilities, are reported in <xref ref-type="fig" rid="pcbi-1002051-g010"><bold>Figure 10a</bold></xref>. In the elongated conformation, the relative orientation of the strands within the sheets tends to be parallel. The two facing strands are oriented antiparallel to each other. The Tyr aromatic side chains form clusters of packed rings that are reminiscent of the arrangements observed in the crystals from the Eisenberg group (<xref ref-type="fig" rid="pcbi-1002051-g010"><bold>Figure 10a</bold></xref>). Partial ordering of the Gln and Asn side chains into the zipper spine arrangement is also observed. This conformation, labeled as conformation 3 in <xref ref-type="fig" rid="pcbi-1002051-g010"><bold>Figure 10a</bold></xref>, is however not stable enough to be the most populated structure at the lowest temperatures. Representative structures of the most representative clusters are extracted from the trajectories and their relative free energies evaluated with the GB/SA approach implemented in the program MacroModel, and according to what was already reported <xref rid="pcbi.1002051-Meli1" ref-type="bibr">[44]</xref>. These calculations provide an approximated energetic value for the stability of the aggregates in solution, and show that the more elongated structures tend to undergo transitions to more globular like conformations.</p><fig id="pcbi-1002051-g010" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002051.g010</object-id><label>Figure 10</label><caption><title>Representative structures of the most populated clusters from all-atom REMD simulations of the 20mer: (a) structure 20GRp-B3 and (b) structure 20OP2-N.</title><p>The number identifying each structure represents the cluster rank (1. being the most populated cluster). The value of the GB/SA energy in water of the complex is reported. Arrows represent transitions between clusters, indicating possible paths between cluster structures.</p></caption><graphic xlink:href="pcbi.1002051.g010"></graphic></fig><p>Several conformational transitions among different structural clusters are observed, and highlighted by direction arrows in <xref ref-type="fig" rid="pcbi-1002051-g010"><bold>Figure 10a</bold></xref>. Summarizing, the starting conformation for 20GRp-B3 with two facing sheets antiparallel to each other shows a high tendency to interconvert rapidly into globular, compact structures. The high conformational plasticity and limited stability of elongated conformations had already been noticed in the present OPEP-based simulations. Analysis of the time evolution of REMD trajectories shows that globular conformations form by detachment and re-docking of β-sheet motifs. During the formation of the more globular structures a slight decrease of the amount of ordered β-sheet could also be noticed.</p><p>All-atom REMD simulations were also used to analyze the structural evolution of cluster 20OP2-N. The representative structures obtained from the all-atom simulation 20GR2-N2 are reported in <xref ref-type="fig" rid="pcbi-1002051-g010"><bold>Figure 10b</bold></xref>, with their relative stabilities. The packing of the interior of the initial aggregate is not optimal. Consistently with the previous case, the overall aggregate still shows a high degree of structural plasticity undergoing transitions to conformations with a lower degree of ordered β sheet and more globular-like shapes. Elongated structures featuring two facing β-strands still appear and show partial interdigitation of Asn and Gln side chains in the inter-strand space, with formation of a steric zipper motif. This does not seem sufficient, however, to provide the necessary stabilization for the aggregate to evolve to a protofibril structure. Once more, a number of transitions among different structural families (<xref ref-type="fig" rid="pcbi-1002051-g010"><bold>Figure 10b</bold></xref>) are observed.</p><p>Overall, we observe that the elongated structures reminiscent of the one observed by the Eisenberg group in the fibril microcrystals is accessible at room temperature but is not stable and can easily interconvert into globular or more disordered conformations, even in the presence of 20 monomers (see <xref ref-type="table" rid="pcbi-1002051-t003"><bold>Table 3</bold></xref>). According to our observations, this happens both by detachment-reattachment of strands or small oligomers from existing structures, or through sliding/reptation moves without detaching from the aggregates <xref rid="pcbi.1002051-Wei2" ref-type="bibr">[81]</xref>, <xref rid="pcbi.1002051-Santini1" ref-type="bibr">[82]</xref>. This behavior is widely observed across our 20-mer simulations while the 12-mer system does not undergo much reorganization during all-atom MD simulations. Interestingly, when ordered supramolecular structure forms, with either an elongated or a compact globular shape, there is no space for water in the interior of the aggregate. Water exclusion from the core of the aggregate may represent the first step leading to the formation of a dry interface.</p></sec></sec><sec id="s3d"><title>Conclusions</title><p>The self-organization process of peptides and proteins into oligomeric soluble and insoluble aggregates that eventually evolve to fibrils is still difficult, if not impossible, to study at atomic resolution using experimental approaches. In this paper, we have carried out an extensive and comprehensive study of the formation of oligomers of the model peptide GNNQQNY in different conditions combining coarse-grained and all-atom simulation approaches. Different numbers of peptides were used in several simulations. In the smallest systems, composed of three strands, a diverse set of structural motifs is accessible at room temperatures. When bigger systems consisting of 12 chains are analyzed, compact and globular structures begin to appear. Interestingly, in some cases, globular oligomers expose hydrophilic side chains to the contact with the water solvent, providing a viable model for soluble intermediates that have been observed on pathway to the formation of the final fibril. In parallel, at 300 K, globular structures exposing a large amount of hydrophobic surface also appear. These may represent possible nuclei for the growth of bigger supramolecular structures.</p><p>In simulations conducted using 20 monomers, we have noticed the appearance of elongated structures characterized by the juxtaposition of two mainly parallel β-sheets with partial interdigitation of amidic side chains reminiscent of the zipper-spine observed in fibril microcrystals. It is important to notice, however, that these structures are not stable in water solvent and evolve towards more globular conformations. This observation suggests that while ordered fibril-like structures are accessible on the energy landscape, they need further stabilization by establishing contacts with multiple copies of similar structures in order to evolve to a fully fibrillar geometry. In this context, the formation of this geometry would require the constructive interplay of many factors and the entropic expense of such process would be clearly very high, explaining the long lag phase times and very slow kinetics of amyloid fibril formation. Moreover, the rich variety of structures and conformational changes observed for the aggregates may also reverberate into the fibril polymorphism observed at the experimental scale.</p><p>In summary, our data and structural models represent valid complements to experimental approaches in the attempt to shed light on the supramolecular arrangements of amyloidogenic oligomers, and lead the following conclusions.</p><p>First, the 20-mers of GNNQQNY are polymorphic and endowed with a high degree of structural plasticity. Polymorphism of the fibrillar products of amyloid aggregation has been observed for many sequences by X-ray diffraction and solid-state NMR experiments <xref rid="pcbi.1002051-Sawaya1" ref-type="bibr">[6]</xref>, <xref rid="pcbi.1002051-Paravastu1" ref-type="bibr">[14]</xref>, <xref rid="pcbi.1002051-Petkova2" ref-type="bibr">[21]</xref> and a recent computational study using a simplified description, where the peptide has a single degree of freedom with two minima <xref rid="pcbi.1002051-Pellarin1" ref-type="bibr">[83]</xref>. Polymorphism of the low-n oligomers with different registers of inter-peptide H-bonds and orientations of the chains has also been reported using various force fields and computational procedures <xref rid="pcbi.1002051-Mo1" ref-type="bibr">[22]</xref>, <xref rid="pcbi.1002051-DeSimone1" ref-type="bibr">[41]</xref>, <xref rid="pcbi.1002051-Meli1" ref-type="bibr">[44]</xref>, <xref rid="pcbi.1002051-Strodel1" ref-type="bibr">[46]</xref>, <xref rid="pcbi.1002051-Santini2" ref-type="bibr">[84]</xref>.</p><p>Second, the 20-mers of GNNQQNY in explicit water are in dynamic equilibrium, within at least 100 ns, between amorphous structures (high probability) and configurations with three β-sheets in various orientations (medium probability) and two β-sheets (low probability). These two-β sheets, reminiscent of the cross-β structures and the dry steric zipper observed experimentally for mature fibrils, are not parallel, however, suggesting the existence of a free energy barrier preventing the formation of a perfectly packed steric zipper.</p><p>Third, there is a reorientation of the β-strands between the GNNQQNY oligomers and fibrils. We find that an anti-parallel β-strand alignment dominates over the parallel one in the 3 and 12 peptide systems. This contradiction with the fibrillar parallel β-strand orientation <xref rid="pcbi.1002051-Nelson1" ref-type="bibr">[7]</xref> is however reconciled by the 20-peptide systems, where a significant increase in the amount of parallel strands and fraction of fully parallel β-sheets is observed.</p><p>Fourth, a common observation is that short amyloid peptide fragments assume antiparallel β-strand geometries whereas longer peptides, and proteins, often assume parallel geometries. Our simulations along with other recent studies show this geometrical property is more complex and depends strongly on the amino acid composition. The dependence of β-strand orientation with oligomer size occurs in the GNNQQNY (Sup35) peptide and the VQIVYK (PH6) peptide, as reported by another computational study <xref rid="pcbi.1002051-Li1" ref-type="bibr">[85]</xref>, because both short peptide fragments do not contain any opposite charged amino acids at the extremities as opposed to numerous fragments, such as the KLVFFAE (Aβ(16–22)) and the KFFE peptides. <xref ref-type="table" rid="pcbi-1002051-t002">Table 2</xref>, for example, shows that the free energy difference in favor of the antiparallel sheet decreases from 0.70 kcal/mol for the trimer to only 0.05 kcal/mol for the dodecamer at 300 K. This competition between the two β-strand orientations during polymerization is also supported experimentally by the co-existence of either antiparallel or parallel strands in seeded hIAPP20–29 (SNNFGAILSS) fibrils <xref rid="pcbi.1002051-Madine1" ref-type="bibr">[15]</xref>, and the fact that D23N-Aβ40 forms fibrils with a majority having antiparallel β-sheet structures and a minority having parallel β-sheet structures <xref rid="pcbi.1002051-Tycko1" ref-type="bibr">[86]</xref>. We hypothesize that both short polar peptides and longer peptides could display this β-strand orientation transition during aggregation.</p><p>In addition, antiparallel β-sheets allow a higher potential variability of the inter-chain H-bond geometry <xref rid="pcbi.1002051-Salemme1" ref-type="bibr">[87]</xref>, <xref rid="pcbi.1002051-Salemme2" ref-type="bibr">[88]</xref> accommodating a higher number of possible strand conformations, which is the favored situation in smaller aggregates. A parallel organization could be more favorable for larger aggregates, which can for instance create double-layered sheets that can pack with a parallel arrangement of the strands as indicated by SS-NMR <xref rid="pcbi.1002051-Antzutkin1" ref-type="bibr">[89]</xref>. Moreover, in larger aggregates, the parallel arrangement might be stabilized by the definition and stabilization of large numbers of side-chain side-chain contacts. Therefore, our results are in agreement with experiments but also suggest a more complex relationship between the monomers during the assembly process, something that cannot be measured experimentally.</p><p>Fifth, in terms of experimental relevance, it is important to note that evidence exists showing that aggregation pathways can be manipulated by the use of molecular chaperones. In the case of the Sup35 prion protein, the chaperone Hsp104 catalyzes the polymerization of seeds that are crucial for efficient amyloid formation <xref rid="pcbi.1002051-Shorter1" ref-type="bibr">[90]</xref>. Different chaperones such as Hsp70 and Hsp40, on the contrary, prevent the assembly of aggregating species when added during the polymerization reactions. In these cases, soluble aggregates are formed, showing that chaperones can redirect amyloidogenic polypeptides into non-amyloidogenic species. Finally, the chaperonin TRiC has also been shown to be a modulator of amyloid formation <xref rid="pcbi.1002051-Behrends1" ref-type="bibr">[91]</xref>. TRiC, in combination with Hsp70 and 40, stimulates the reassembly of huntingtin oligomers into soluble species, which are non-toxic. Clearly, the type of binding and structural remodeling determined by different chaperones is dictated by the details of the molecular recognition between the oligomers and the interaction surfaces of the chaperones. In this context, it is worth noting that generating plausible models of possible oligomeric substrates for the chaperones may be of great help in the design of optimized systems aimed at modulating aggregate properties. Based on the structural and surface properties of the oligomers, models of their complexes with chaperones may be generated.</p><p>On this basis, the hydrophobic-hydrophilic profile of the chaperone interaction surfaces could, for instance, be changed by means of site-directed mutagenesis, affecting their activity and ultimately the properties of the remodeled oligomers. This would allow a rational manipulation of the amyloidogenic pathways, helping to shed light on a very complex biological phenomenon.</p><p>A final consideration helpful to put our results in a biological perspective is related to the importance of the knowledge of oligomeric structures in the design of amyloidogenic inhibitors. In this context, we are currently exploring the characterization of the solvent accessible hydrophobic surface area of the 20-mers to guide docking-experiments of small-molecule compounds (Congo Red and EGCG in particular), in order to derive possible rules for the rational selection of aggregation inhibitors. Preliminary data and results show that this could be helpful in alleviating the difficulties associated to drug-design when dealing with amyloid-targets. Indeed, compared to classical drug-design efforts where the target is an active site, with well-defined structure and cavities, the variety of structures, mechanisms and conformational plasticity of oligomers shown here confirm that rational design of aggregation inhibitors is a daunting challenge. However, careful characterization of oligomeric structures provides useful suggestions for the design of possible inhibitors. Selective compounds or peptidomimetics could be designed/selected to target the oligomer conformations characterized by the presence of aromatic groups on their external surface. These compounds would actually target intermediates that are more prone to be insoluble or to favor the addition of monomers through hydrophobic interactions. Interestingly, most of the existing inhibitors of amyloidogenic pathways are small molecules rich in aromatic functionalities, which can target more than one single aggregating species, showing a general mechanism of action <xref rid="pcbi.1002051-Ehrnhoefer1" ref-type="bibr">[92]</xref>.</p><p>Alternatively, one could design peptidomimetic-based or small molecule chaperones that can stabilize soluble species, subtracting them from the amyloidogenic pathway. This would lead to a redirection of otherwise amyloidogenic peptides into non-amyloidogenic species.</p></sec></sec><sec sec-type="supplementary-material" id="s4"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pcbi.1002051.s001"><label>Figure S1</label><caption><p>Time evolution of the radius of gyration of the 12-mer oligomers. From top to bottom: structures 12GR1-A, 12GR1-B, 12GR1-C, 12GR1-D and 12GR1-E. The structures shown are the final structures of the all-atom MD simulations with GROMACS.</p><p>(TIF)</p></caption><media xlink:href="pcbi.1002051.s001.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1002051.s002"><label>Figure S2</label><caption><p>Time evolution of the radius of gyration of the 20-mer oligomers for the preliminary simulation. Structures 20GRp-B1, 20GRp-B2, 20GRp-C1 and 20GRp-C2. The structures shown are the final structures of the all-atom MD simulations with GROMACS.</p><p>(TIF)</p></caption><media xlink:href="pcbi.1002051.s002.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pcbi.1002051.s003"><label>Figure S3</label><caption><p>Time evolution of the radius of gyration of the 20-mer oligomers for the preliminary simulation. Structures 20GRp-A1, 20GRp-A2, 20GRp-D1 and 20GRp-D2. The structures shown are the final structures of the all-atom MD simulations with GROMACS.</p><p>(TIF)</p></caption><media xlink:href="pcbi.1002051.s003.tif"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work was partially supported by grants from the Natural Sciences and Engineering Research Council of Canada (<ext-link ext-link-type="uri" xlink:href="http://www.nserc-crsng.gc.ca/">http://www.nserc-crsng.gc.ca/</ext-link>), the Canada Research Chair Foundation (<ext-link ext-link-type="uri" xlink:href="http://www.chairs-chaires.gc.ca/">http://www.chairs-chaires.gc.ca/</ext-link>), the Fonds de la recherche en santé du Québec (<ext-link ext-link-type="uri" xlink:href="http://www.frsq.gouv.qc.ca/">http://www.frsq.gouv.qc.ca/</ext-link>) and the Réseau québécois de calcul de haute performance (<ext-link ext-link-type="uri" xlink:href="http://rqchp.ca">http://rqchp.ca</ext-link>). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group><ref-list><title>References</title><ref id="pcbi.1002051-Chiti1"><label>1</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chiti</surname><given-names>F</given-names></name><name><surname>Dobson</surname><given-names>CM</given-names></name></person-group><year>2006</year><article-title>Protein misfolding, functional amyloid and human disease.</article-title><source>Annu Rev Biochem</source><volume>75</volume><fpage>333</fpage><lpage>366</lpage><pub-id pub-id-type="pmid">16756495</pub-id></element-citation></ref><ref id="pcbi.1002051-Ross1"><label>2</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ross</surname><given-names>CA</given-names></name><name><surname>Poirier</surname><given-names>MA</given-names></name></person-group><year>2004</year><article-title>Protein aggregation and neurodegenerative disease.</article-title><source>Nat Med</source><volume>10</volume><fpage>S10</fpage><lpage>S17</lpage><pub-id pub-id-type="pmid">15272267</pub-id></element-citation></ref><ref id="pcbi.1002051-Miranker1"><label>3</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miranker</surname><given-names>AD</given-names></name></person-group><year>2004</year><article-title>Unzipping the mysteries of amyloid fibril formation.</article-title><source>Proc Natl Acad Sci U S A</source><volume>101</volume><fpage>4335</fpage><lpage>4336</lpage><pub-id pub-id-type="pmid">15070716</pub-id></element-citation></ref><ref id="pcbi.1002051-Zerovnik1"><label>4</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zerovnik</surname><given-names>E</given-names></name></person-group><year>2002</year><article-title>Amyloid-fibril formation. Proposed mechanisms and relevance to conformational disease.</article-title><source>Eur J Biochem</source><volume>269</volume><fpage>3362</fpage><lpage>3371</lpage><pub-id pub-id-type="pmid">12135474</pub-id></element-citation></ref><ref id="pcbi.1002051-Dobson1"><label>5</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dobson</surname><given-names>CM</given-names></name></person-group><year>1999</year><article-title>Protein misfolding, evolution and disease.</article-title><source>Trends Biochem Sci</source><volume>24</volume><fpage>329</fpage><lpage>332</lpage><pub-id pub-id-type="pmid">10470028</pub-id></element-citation></ref><ref id="pcbi.1002051-Sawaya1"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sawaya</surname><given-names>MR</given-names></name><name><surname>Sambashivan</surname><given-names>S</given-names></name><name><surname>Nelson</surname><given-names>R</given-names></name><name><surname>Ivanova</surname><given-names>MI</given-names></name><name><surname>Sievers</surname><given-names>SA</given-names></name><etal></etal></person-group><year>2007</year><article-title>Atomic structures of amyloid cross-β spines reveal varied steric zippers.</article-title><source>Nature</source><volume>447</volume><fpage>453</fpage><lpage>457</lpage><pub-id pub-id-type="pmid">17468747</pub-id></element-citation></ref><ref id="pcbi.1002051-Nelson1"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nelson</surname><given-names>R</given-names></name><name><surname>Sawaya</surname><given-names>MR</given-names></name><name><surname>Balbirnie</surname><given-names>M</given-names></name><name><surname>Madsen</surname><given-names>AØ</given-names></name><name><surname>Riekel</surname><given-names>C</given-names></name><etal></etal></person-group><year>2005</year><article-title>Structure of the cross-β spine of amyloid-like fibrils.</article-title><source>Nature</source><volume>435</volume><fpage>773</fpage><lpage>778</lpage><pub-id pub-id-type="pmid">15944695</pub-id></element-citation></ref><ref id="pcbi.1002051-Kayed1"><label>8</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kayed</surname><given-names>R</given-names></name><name><surname>Head</surname><given-names>E</given-names></name><name><surname>Thompson</surname><given-names>JL</given-names></name><name><surname>McIntire</surname><given-names>TM</given-names></name><name><surname>Milton</surname><given-names>SC</given-names></name><etal></etal></person-group><year>2003</year><article-title>Common structure of soluble amyloid oligomers implies common mechanism of pathogenesis.</article-title><source>Science</source><volume>300</volume><fpage>486</fpage><lpage>489</lpage><pub-id pub-id-type="pmid">12702875</pub-id></element-citation></ref><ref id="pcbi.1002051-Thirumalai1"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thirumalai</surname><given-names>D</given-names></name><name><surname>Klimov</surname><given-names>DK</given-names></name><name><surname>Dima</surname><given-names>RI</given-names></name></person-group><year>2003</year><article-title>Emerging ideas on the molecular basis of protein and peptide aggregation.</article-title><source>Curr Opin Struct Biol</source><volume>13</volume><fpage>146</fpage><lpage>159</lpage><pub-id pub-id-type="pmid">12727507</pub-id></element-citation></ref><ref id="pcbi.1002051-Sunde1"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sunde</surname><given-names>M</given-names></name><name><surname>Serpell</surname><given-names>LC</given-names></name><name><surname>Bartlam</surname><given-names>M</given-names></name><name><surname>Fraser</surname><given-names>PE</given-names></name><name><surname>Pepys</surname><given-names>MB</given-names></name><etal></etal></person-group><year>1997</year><article-title>Common structure of amyloid fibrils by synchrotron X-ray diffraction.</article-title><source>J Mol Biol</source><volume>273</volume><fpage>729</fpage><lpage>739</lpage><pub-id pub-id-type="pmid">9356260</pub-id></element-citation></ref><ref id="pcbi.1002051-Fndrich1"><label>11</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fändrich</surname><given-names>M</given-names></name><name><surname>Meinhardt</surname><given-names>J</given-names></name><name><surname>Grigorieff</surname><given-names>N</given-names></name></person-group><year>2009</year><article-title>Structural polymorphism of Alzheimer Aβ and other amyloid fibrils.</article-title><source>Prion</source><volume>3</volume><fpage>89</fpage><lpage>93</lpage><pub-id pub-id-type="pmid">19597329</pub-id></element-citation></ref><ref id="pcbi.1002051-Celej1"><label>12</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Celej</surname><given-names>MS</given-names></name><name><surname>Caarls</surname><given-names>W</given-names></name><name><surname>Demchenko</surname><given-names>AP</given-names></name><name><surname>Jovin</surname><given-names>TM</given-names></name></person-group><year>2009</year><article-title>A triple emission fluorescent probe reveals distinctive amyloid fibrillar polymorphism of wild-type α-synuclein and its familial Parkinsons disease-mutants.</article-title><source>Biochem</source><volume>48</volume><fpage>7465</fpage><lpage>7472</lpage><pub-id pub-id-type="pmid">19586054</pub-id></element-citation></ref><ref id="pcbi.1002051-Meinhardt1"><label>13</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meinhardt</surname><given-names>J</given-names></name><name><surname>Sachse</surname><given-names>C</given-names></name><name><surname>Hortschansky</surname><given-names>P</given-names></name><name><surname>Grigorieff</surname><given-names>N</given-names></name><name><surname>Fändrich</surname><given-names>M</given-names></name></person-group><year>2009</year><article-title>Aβ(1–40) fibril polymorphism implies diverse interaction patterns in amyloid fibrils.</article-title><source>J Mol Biol</source><volume>386</volume><fpage>869</fpage><lpage>877</lpage><pub-id pub-id-type="pmid">19038266</pub-id></element-citation></ref><ref id="pcbi.1002051-Paravastu1"><label>14</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Paravastu</surname><given-names>AK</given-names></name><name><surname>Leapman</surname><given-names>RD</given-names></name><name><surname>Yau</surname><given-names>WM</given-names></name><name><surname>Tycko</surname><given-names>R</given-names></name></person-group><year>2008</year><article-title>Molecular structural basis for polymorphism in Alzheimer's beta-amyloid fibrils.</article-title><source>Proc Natl Acad Sci U S A</source><volume>105</volume><fpage>18349</fpage><lpage>54</lpage><pub-id pub-id-type="pmid">19015532</pub-id></element-citation></ref><ref id="pcbi.1002051-Madine1"><label>15</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Madine</surname><given-names>J</given-names></name><name><surname>Jack</surname><given-names>E</given-names></name><name><surname>Stockley</surname><given-names>PG</given-names></name><name><surname>Radford</surname><given-names>SE</given-names></name><name><surname>Serpell</surname><given-names>LC</given-names></name><etal></etal></person-group><year>2008</year><article-title>Structural insights into the polymorphism of amyloid-like fibrils formed by region 20–29 of amylin revealed by solid-state NMR and X-ray fiber diffraction.</article-title><source>J Am Chem Soc</source><volume>130</volume><fpage>14990</fpage><lpage>5001</lpage><pub-id pub-id-type="pmid">18937465</pub-id></element-citation></ref><ref id="pcbi.1002051-Verel1"><label>16</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Verel</surname><given-names>R</given-names></name><name><surname>Tomka</surname><given-names>IT</given-names></name><name><surname>Bertozzi</surname><given-names>C</given-names></name><name><surname>Cadalbert</surname><given-names>R</given-names></name><name><surname>Kammerer</surname><given-names>RA</given-names></name><etal></etal></person-group><year>2008</year><article-title>Polymorphism in an amyloid-like fibril-forming model peptide.</article-title><source>Angew Chem Int Ed Engl</source><volume>47</volume><fpage>5842</fpage><lpage>5</lpage><pub-id pub-id-type="pmid">18528917</pub-id></element-citation></ref><ref id="pcbi.1002051-Kodali1"><label>17</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kodali</surname><given-names>R</given-names></name><name><surname>Wetzel</surname><given-names>R</given-names></name></person-group><year>2007</year><article-title>Polymorphism in the intermediates and products of amyloid assembly.</article-title><source>Curr Opin Struct Biol</source><volume>17</volume><fpage>48</fpage><lpage>57</lpage><pub-id pub-id-type="pmid">17251001</pub-id></element-citation></ref><ref id="pcbi.1002051-Petkova1"><label>18</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Petkova</surname><given-names>AT</given-names></name><name><surname>Yau</surname><given-names>WM</given-names></name><name><surname>Tycko</surname><given-names>R</given-names></name></person-group><year>2006</year><article-title>Experimental constraints on quaternary structure in Alzheimer's beta-amyloid fibrils.</article-title><source>Biochem</source><volume>45</volume><fpage>498</fpage><lpage>512</lpage><pub-id pub-id-type="pmid">16401079</pub-id></element-citation></ref><ref id="pcbi.1002051-Kreplak1"><label>19</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kreplak</surname><given-names>L</given-names></name><name><surname>Aebi</surname><given-names>U</given-names></name></person-group><year>2006</year><article-title>From the polymorphism of amyloid fibrils to their assembly mechanism and cytotoxicity.</article-title><source>Adv Protein Chem</source><volume>73</volume><fpage>217</fpage><lpage>33</lpage><pub-id pub-id-type="pmid">17190615</pub-id></element-citation></ref><ref id="pcbi.1002051-Anderson1"><label>20</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Anderson</surname><given-names>M</given-names></name><name><surname>Bocharova</surname><given-names>OV</given-names></name><name><surname>Makarava</surname><given-names>N</given-names></name><name><surname>Breydo</surname><given-names>L</given-names></name><name><surname>Salnikov</surname><given-names>VV</given-names></name><etal></etal></person-group><year>2006</year><article-title>Polymorphism and ultrastructural organizations of prion protein amyloid fibrils: an insight from high-resolution atomic force microscopy.</article-title><source>J Mol Biol</source><volume>358</volume><fpage>580</fpage><lpage>96</lpage><pub-id pub-id-type="pmid">16519898</pub-id></element-citation></ref><ref id="pcbi.1002051-Petkova2"><label>21</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Petkova</surname><given-names>AT</given-names></name><name><surname>Leapman</surname><given-names>RD</given-names></name><name><surname>Guo</surname><given-names>WM</given-names></name><name><surname>Mattson</surname><given-names>MP</given-names></name><name><surname>Tycko</surname><given-names>R</given-names></name></person-group><year>2005</year><article-title>Self-propagating, molecular-level polymorphism in Alzheimer's beta-amyloid fibrils.</article-title><source>Science</source><volume>307</volume><fpage>262</fpage><lpage>5</lpage><pub-id pub-id-type="pmid">15653506</pub-id></element-citation></ref><ref id="pcbi.1002051-Mo1"><label>22</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mo</surname><given-names>Y</given-names></name><name><surname>Lu</surname><given-names>Y</given-names></name><name><surname>Wei</surname><given-names>G</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2009</year><article-title>Structural diversity of the soluble trimers of the human amylin(20–29) peptide revealed by molecular dynamics simulations.</article-title><source>J Chem Phys</source><volume>130</volume><fpage>125101</fpage><pub-id pub-id-type="pmid">19334894</pub-id></element-citation></ref><ref id="pcbi.1002051-Park1"><label>23</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Park</surname><given-names>J</given-names></name><name><surname>Kahng</surname><given-names>B</given-names></name><name><surname>Hwang</surname><given-names>W</given-names></name></person-group><year>2009</year><article-title>Thermodynamic selection of steric zipper patterns in the amyloid cross-beta spine.</article-title><source>PLoS Comput Biol</source><volume>5</volume><fpage>e1000492</fpage><pub-id pub-id-type="pmid">19730673</pub-id></element-citation></ref><ref id="pcbi.1002051-Berryman1"><label>24</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berryman</surname><given-names>JT</given-names></name><name><surname>Radford</surname><given-names>SE</given-names></name><name><surname>Harris</surname><given-names>SA</given-names></name></person-group><year>2009</year><article-title>Thermodynamic description of polymorphism in Q- and N-rich peptide aggregates revealed by atomistic simulation.</article-title><source>Biophys J</source><volume>97</volume><fpage>1</fpage><lpage>11</lpage><pub-id pub-id-type="pmid">19580739</pub-id></element-citation></ref><ref id="pcbi.1002051-Klein1"><label>25</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Klein</surname><given-names>WL</given-names></name><name><surname>Stine</surname><given-names>WB</given-names></name><name><surname>Teplow</surname><given-names>DB</given-names></name></person-group><year>2004</year><article-title>Small assemblies of unmodified amyloid β-protein are the proximate neurotoxin in Alzheimer's disease.</article-title><source>Neurobiol Aging</source><volume>25</volume><fpage>569</fpage><lpage>580</lpage><pub-id pub-id-type="pmid">15172732</pub-id></element-citation></ref><ref id="pcbi.1002051-Bucciantini1"><label>26</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bucciantini</surname><given-names>M</given-names></name><name><surname>Giannoni</surname><given-names>E</given-names></name><name><surname>Chiti</surname><given-names>F</given-names></name><name><surname>Baroni</surname><given-names>F</given-names></name><name><surname>Formigli</surname><given-names>L</given-names></name><etal></etal></person-group><year>2002</year><article-title>Inherent toxicity of aggregates implies a common mechanism for protein misfolding diseases.</article-title><source>Nature</source><volume>416</volume><fpage>507</fpage><lpage>511</lpage><pub-id pub-id-type="pmid">11932737</pub-id></element-citation></ref><ref id="pcbi.1002051-Shankar1"><label>27</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shankar</surname><given-names>GM</given-names></name><name><surname>Li</surname><given-names>S</given-names></name><name><surname>Mehta</surname><given-names>TH</given-names></name><name><surname>Garcia-Munoz</surname><given-names>A</given-names></name><name><surname>Shepardson</surname><given-names>NE</given-names></name><etal></etal></person-group><year>2008</year><article-title>Amyloid-β protein dimers isolated directly from Alzheimer's brains impair synaptic plasticity and memory.</article-title><source>Nat Med</source><volume>14</volume><fpage>837</fpage><lpage>842</lpage><pub-id pub-id-type="pmid">18568035</pub-id></element-citation></ref><ref id="pcbi.1002051-Ono1"><label>28</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ono</surname><given-names>K</given-names></name><name><surname>Condron</surname><given-names>MM</given-names></name><name><surname>Teplow</surname><given-names>DB</given-names></name></person-group><year>2009</year><article-title>Structure-neurotoxicity relationships of amyloid beta-protein oligomers.</article-title><source>Proc Natl Acad Sci U S A</source><volume>106</volume><fpage>14745</fpage><lpage>50</lpage><pub-id pub-id-type="pmid">19706468</pub-id></element-citation></ref><ref id="pcbi.1002051-Lashuel1"><label>29</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lashuel</surname><given-names>HA</given-names></name><name><surname>Hartley</surname><given-names>D</given-names></name><name><surname>Petre</surname><given-names>BM</given-names></name><name><surname>Walz</surname><given-names>T</given-names></name><name><surname>Lansbury</surname><given-names>PT</given-names></name></person-group><year>2002</year><article-title>Amyloid pores from pathogenic mutations.</article-title><source>Nature</source><volume>418</volume><fpage>291</fpage></element-citation></ref><ref id="pcbi.1002051-Goldberg1"><label>30</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goldberg</surname><given-names>MS</given-names></name><name><surname>Lansbury</surname><given-names>PT</given-names></name></person-group><year>2000</year><article-title>Is there a cause-and-effect relationship between alpha-synuclein fibrillization and Parkinson's disease?</article-title><source>Nature Cell Biol</source><volume>2</volume><fpage>E115</fpage><lpage>E119</lpage><pub-id pub-id-type="pmid">10878819</pub-id></element-citation></ref><ref id="pcbi.1002051-Lansbury1"><label>31</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lansbury</surname><given-names>PT</given-names></name></person-group><year>1999</year><article-title>Evolution of amyloid: what normal protein folding may tell us about fibrillogenesis and disease.</article-title><source>Proc Natl Acad Sci U S A</source><volume>96</volume><fpage>3342</fpage><lpage>3344</lpage><pub-id pub-id-type="pmid">10097040</pub-id></element-citation></ref><ref id="pcbi.1002051-Zheng1"><label>32</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zheng</surname><given-names>J</given-names></name><name><surname>Jang</surname><given-names>H</given-names></name><name><surname>Ma</surname><given-names>B</given-names></name><name><surname>Nussinov</surname><given-names>R</given-names></name></person-group><year>2008</year><article-title>Annular structures as intermediates in fibril formation of Alzheimer's Aβ-17–42.</article-title><source>J Phys Chem B</source><volume>112</volume><fpage>6856</fpage><lpage>6865</lpage><pub-id pub-id-type="pmid">18457440</pub-id></element-citation></ref><ref id="pcbi.1002051-Xue1"><label>33</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xue</surname><given-names>WF</given-names></name><name><surname>Hellewell</surname><given-names>AL</given-names></name><name><surname>Gosal</surname><given-names>WS</given-names></name><name><surname>Homans</surname><given-names>SW</given-names></name><name><surname>Hewitt</surname><given-names>EW</given-names></name><etal></etal></person-group><year>2009</year><article-title>Fibril fragmentation enhances amyloid cytotoxicity.</article-title><source>J Biol Chem</source><volume>284</volume><fpage>34272</fpage><lpage>34282</lpage><pub-id pub-id-type="pmid">19808677</pub-id></element-citation></ref><ref id="pcbi.1002051-Knowles1"><label>34</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Knowles</surname><given-names>TP</given-names></name><name><surname>Waudby</surname><given-names>CA</given-names></name><name><surname>Devlin</surname><given-names>GL</given-names></name><name><surname>Cohen</surname><given-names>SI</given-names></name><name><surname>Aguzzi</surname><given-names>A</given-names></name><etal></etal></person-group><year>2009</year><article-title>An analytical solution to the kinetics of breakable filament assembly.</article-title><source>Science</source><volume>326</volume><fpage>1533</fpage><lpage>1537</lpage><pub-id pub-id-type="pmid">20007899</pub-id></element-citation></ref><ref id="pcbi.1002051-Balbirnie1"><label>35</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Balbirnie</surname><given-names>M</given-names></name><name><surname>Grothe</surname><given-names>R</given-names></name><name><surname>Eisenberg</surname><given-names>DS</given-names></name></person-group><year>2001</year><article-title>An amyloid-forming peptide from the yeast prion Sup35 reveals a dehydrated β-sheet structure for amyloid.</article-title><source>Proc Natl Acad Sci U S A</source><volume>98</volume><fpage>2375</fpage><lpage>2380</lpage><pub-id pub-id-type="pmid">11226247</pub-id></element-citation></ref><ref id="pcbi.1002051-DiazAvalos1"><label>36</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Diaz-Avalos</surname><given-names>R</given-names></name><name><surname>Long</surname><given-names>C</given-names></name><name><surname>Fontano</surname><given-names>E</given-names></name><name><surname>Balbirnie</surname><given-names>M</given-names></name><name><surname>Grothe</surname><given-names>R</given-names></name><etal></etal></person-group><year>2003</year><article-title>Cross-beta order and diversity in nanocrystals of an amyloid-forming peptide.</article-title><source>J Mol Biol</source><volume>330</volume><fpage>1165</fpage><lpage>1175</lpage><pub-id pub-id-type="pmid">12860136</pub-id></element-citation></ref><ref id="pcbi.1002051-VanderWel1"><label>37</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Van der Wel</surname><given-names>PCA</given-names></name><name><surname>Lewandowski</surname><given-names>JR</given-names></name><name><surname>Griffin R</surname><given-names>G</given-names></name></person-group><year>2006</year><article-title>Solid-state NMR study of amyloid nanocrystals and fibrils formed by the peptide GNNQQNY from yeast prion protein Sup35p.</article-title><source>J Am Chem Soc</source><volume>129</volume><fpage>5117</fpage><lpage>5130</lpage></element-citation></ref><ref id="pcbi.1002051-Wei1"><label>38</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wei</surname><given-names>G</given-names></name><name><surname>Song</surname><given-names>W</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name></person-group><year>2008</year><article-title>Self-assembly of amyloid-forming peptides using molecular dynamics simulations and the OPEP coarse-grained force field.</article-title><source>Frontiers in Biosciences</source><volume>13</volume><fpage>5681</fpage><lpage>5692</lpage></element-citation></ref><ref id="pcbi.1002051-Song1"><label>39</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Song</surname><given-names>W</given-names></name><name><surname>Wei</surname><given-names>G</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2008</year><article-title>Self-assembly of the beta<sub>2</sub>-microglobulin NHVTLSQ peptide using coarse-grained protein model reveals β-barrel species.</article-title><source>J Phys Chem B</source><volume>112</volume><fpage>4410</fpage><lpage>4418</lpage><pub-id pub-id-type="pmid">18341325</pub-id></element-citation></ref><ref id="pcbi.1002051-Irbck1"><label>40</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Irbäck</surname><given-names>A</given-names></name><name><surname>Mitternacht</surname><given-names>S</given-names></name></person-group><year>2007</year><article-title>Spontaneous β-barrel formation: an all-atom Monte-Carlo study of Aβ-16–22 oligomerization.</article-title><source>Proteins</source><volume>71</volume><fpage>207</fpage><lpage>214</lpage></element-citation></ref><ref id="pcbi.1002051-DeSimone1"><label>41</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>De Simone</surname><given-names>A</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2010</year><article-title>Low molecular weight oligomers of amyloid peptides display beta-barrel conformations: a replica exchange molecular dynamics study in explicit solvent.</article-title><source>J Chem Phys</source><volume>132</volume><fpage>165103</fpage><pub-id pub-id-type="pmid">20441311</pub-id></element-citation></ref><ref id="pcbi.1002051-Periole1"><label>42</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Periole</surname><given-names>X</given-names></name><name><surname>Rampioni</surname><given-names>A</given-names></name><name><surname>Vendruscolo</surname><given-names>M</given-names></name><name><surname>Mark</surname><given-names>AE</given-names></name></person-group><year>2009</year><article-title>Factors that affect the degree of twist in beta-sheet structures: a molecular dynamics simulation study of a cross-beta filament of the GNNQQNY peptide.</article-title><source>J Phys Chem B</source><volume>113</volume><fpage>1728</fpage><lpage>1737</lpage><pub-id pub-id-type="pmid">19154133</pub-id></element-citation></ref><ref id="pcbi.1002051-Wang1"><label>43</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>J</given-names></name><name><surname>Tan</surname><given-names>C</given-names></name><name><surname>Chen</surname><given-names>HF</given-names></name><name><surname>Luo</surname><given-names>R</given-names></name></person-group><year>2008</year><article-title>All-atom simulations of amyloid fibrils disaggregation.</article-title><source>Biophys J</source><volume>95</volume><fpage>5037</fpage><lpage>5047</lpage><pub-id pub-id-type="pmid">18757563</pub-id></element-citation></ref><ref id="pcbi.1002051-Meli1"><label>44</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meli</surname><given-names>M</given-names></name><name><surname>Morra</surname><given-names>G</given-names></name><name><surname>Colombo</surname><given-names>G</given-names></name></person-group><year>2008</year><article-title>Investigating the mechanism of peptide aggregation: Insights from mixed Monte-Carlo molecular dynamics simulations.</article-title><source>Biophys J</source><volume>94</volume><fpage>4414</fpage><lpage>4426</lpage><pub-id pub-id-type="pmid">18263661</pub-id></element-citation></ref><ref id="pcbi.1002051-Zhang1"><label>45</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>Z</given-names></name><name><surname>Chen</surname><given-names>H</given-names></name><name><surname>Bai</surname><given-names>H</given-names></name><name><surname>Lai</surname><given-names>L</given-names></name></person-group><year>2007</year><article-title>Molecular dynamics simulations on the oligomer-formation process of the GNNQQNY peptide from yeast prion protein Sup35.</article-title><source>Biophys J</source><volume>93</volume><fpage>1484</fpage><lpage>1492</lpage><pub-id pub-id-type="pmid">17483185</pub-id></element-citation></ref><ref id="pcbi.1002051-Strodel1"><label>46</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Strodel</surname><given-names>B</given-names></name><name><surname>Whittleston</surname><given-names>CS</given-names></name><name><surname>Wales</surname><given-names>DJ</given-names></name></person-group><year>2007</year><article-title>Thermodynamics and kinetics of aggregation for the GNNQQNY peptide.</article-title><source>J Am Chem Soc</source><volume>129</volume><fpage>16005</fpage><lpage>16014</lpage><pub-id pub-id-type="pmid">18052168</pub-id></element-citation></ref><ref id="pcbi.1002051-Reddy1"><label>47</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reddy</surname><given-names>AS</given-names></name><name><surname>Chopra</surname><given-names>M</given-names></name><name><surname>de Pablo</surname><given-names>JJ</given-names></name></person-group><year>2010</year><article-title>GNNQQNY—investigation of early steps during amyloid formation.</article-title><source>Biophys J</source><volume>98</volume><fpage>1038</fpage><lpage>1045</lpage><pub-id pub-id-type="pmid">20303861</pub-id></element-citation></ref><ref id="pcbi.1002051-Vitagliano1"><label>48</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vitagliano</surname><given-names>L</given-names></name><name><surname>Esposito</surname><given-names>L</given-names></name><name><surname>Pedone</surname><given-names>C</given-names></name><name><surname>De Simone</surname><given-names>A</given-names></name></person-group><year>2008</year><article-title>Stability of single sheet GNNQQNY aggregates analyzed by replica exchange molecular dynamics: Antiparallel versus parallel association.</article-title><source>Biochem Biophys Res Commun</source><volume>377</volume><fpage>1036</fpage><lpage>1041</lpage><pub-id pub-id-type="pmid">18938138</pub-id></element-citation></ref><ref id="pcbi.1002051-Zheng2"><label>49</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zheng</surname><given-names>J</given-names></name><name><surname>Ma</surname><given-names>B</given-names></name><name><surname>Tsai</surname><given-names>CJ</given-names></name><name><surname>Nussinov</surname><given-names>R</given-names></name></person-group><year>2006</year><article-title>Structural stability and dynamics of an amyloid-forming peptide GNNQQNY from the yeast prion Sup35.</article-title><source>Biophys J</source><volume>91</volume><fpage>824</fpage><lpage>833</lpage><pub-id pub-id-type="pmid">16679374</pub-id></element-citation></ref><ref id="pcbi.1002051-Esposito1"><label>50</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Esposito</surname><given-names>L</given-names></name><name><surname>Pedone</surname><given-names>C</given-names></name><name><surname>Vitagliano</surname><given-names>L</given-names></name></person-group><year>2006</year><article-title>Molecular dynamics analyses of cross-β-spine steric zipper models: β-sheet twisting and aggregation.</article-title><source>Proc Natl Acad Sci U S A</source><volume>103</volume><fpage>11533</fpage><lpage>11538</lpage><pub-id pub-id-type="pmid">16864786</pub-id></element-citation></ref><ref id="pcbi.1002051-Gsponer1"><label>51</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gsponer</surname><given-names>J</given-names></name><name><surname>Haberthür</surname><given-names>U</given-names></name><name><surname>Caflisch</surname><given-names>A</given-names></name></person-group><year>2003</year><article-title>The role of side-chain interactions in the early steps of aggregation: Molecular dynamics simulations of an amyloid-forming peptide from the yeast prion Sup35.</article-title><source>Proc Natl Acad Sci U S A</source><volume>100</volume><fpage>5154</fpage><lpage>5159</lpage><pub-id pub-id-type="pmid">12700355</pub-id></element-citation></ref><ref id="pcbi.1002051-Sugita1"><label>52</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sugita</surname><given-names>Y</given-names></name><name><surname>Okamoto</surname><given-names>Y</given-names></name></person-group><year>1999</year><article-title>Replica-exchange molecular dynamics method for protein folding.</article-title><source>Chem Phys Lett</source><volume>314</volume><fpage>141</fpage><lpage>151</lpage></element-citation></ref><ref id="pcbi.1002051-Maupetit1"><label>53</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maupetit</surname><given-names>J</given-names></name><name><surname>Tufféry</surname><given-names>P</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2007</year><article-title>A coarse-grained protein force field for folding and structure prediction.</article-title><source>Proteins</source><volume>69</volume><fpage>394</fpage><lpage>408</lpage><pub-id pub-id-type="pmid">17600832</pub-id></element-citation></ref><ref id="pcbi.1002051-Derreumaux1"><label>54</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>1999</year><article-title>From polypeptide sequences to structures using Monte-Carlo simulations and an optimized potential.</article-title><source>J Chem Phys</source><volume>111</volume><fpage>2301</fpage><lpage>2310</lpage></element-citation></ref><ref id="pcbi.1002051-Hess1"><label>55</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hess</surname><given-names>B</given-names></name><name><surname>Kutzner</surname><given-names>C</given-names></name><name><surname>Van Der Spoel</surname><given-names>D</given-names></name><name><surname>Lindahl</surname><given-names>E</given-names></name></person-group><year>2008</year><article-title>GROMACS 4: Algorithms for highly efficient load-balanced, scalable molecular simulation.</article-title><source>J Chem Theory Comput</source><volume>4</volume><fpage>435</fpage><lpage>447</lpage><pub-id pub-id-type="pmid">26620784</pub-id></element-citation></ref><ref id="pcbi.1002051-Chen1"><label>56</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>W</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2006</year><article-title>The conformations of the amyloid-beta(21–30) fragment can be described by three families in solution.</article-title><source>J Chem Phys</source><volume>125</volume><fpage>084911</fpage><pub-id pub-id-type="pmid">16965061</pub-id></element-citation></ref><ref id="pcbi.1002051-Melquiond1"><label>57</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Melquiond</surname><given-names>A</given-names></name><name><surname>Dong</surname><given-names>X</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2008</year><article-title>Role of the region 23–28 in Abeta fibril formation: insights from simulations of the monomers and dimmers of Alzheimer's peptides Abeta40 and Abeta42.</article-title><source>Curr Alzheimer Res</source><volume>5</volume><fpage>244</fpage><lpage>250</lpage><pub-id pub-id-type="pmid">18537541</pub-id></element-citation></ref><ref id="pcbi.1002051-Maupetit2"><label>58</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maupetit</surname><given-names>J</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name><name><surname>Tufféry</surname><given-names>P</given-names></name></person-group><year>2009</year><article-title>PEP-FOLD: an online resource for de novo peptide structure prediction.</article-title><source>Nucleic Acids Res</source><volume>37</volume><fpage>W498</fpage><lpage>503</lpage><pub-id pub-id-type="pmid">19433514</pub-id></element-citation></ref><ref id="pcbi.1002051-Maupetit3"><label>59</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maupetit</surname><given-names>J</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name><name><surname>Tufféry</surname><given-names>P</given-names></name></person-group><year>2010</year><article-title>A fast method for large-scale de novo peptide and miniprotein structure prediction.</article-title><source>J Comput Chem</source><volume>31</volume><fpage>726</fpage><lpage>738</lpage><pub-id pub-id-type="pmid">19569182</pub-id></element-citation></ref><ref id="pcbi.1002051-Berendsen1"><label>60</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berendsen</surname><given-names>HJC</given-names></name><name><surname>Postma</surname><given-names>JPM</given-names></name><name><surname>van Gunsteren</surname><given-names>WF</given-names></name><name><surname>Dinola</surname><given-names>A</given-names></name><name><surname>Haak</surname><given-names>JR</given-names></name></person-group><year>1984</year><article-title>Molecular dynamics with coupling to an external bath.</article-title><source>J Chem Phys</source><volume>81</volume><fpage>3684</fpage><lpage>3690</lpage></element-citation></ref><ref id="pcbi.1002051-Derreumaux2"><label>61</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Derreumaux</surname><given-names>P</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name></person-group><year>2007</year><article-title>Coarse-grained protein molecular dynamics simulations.</article-title><source>J Chem Phys</source><volume>126</volume><fpage>025101</fpage><pub-id pub-id-type="pmid">17228975</pub-id></element-citation></ref><ref id="pcbi.1002051-Chebaro1"><label>62</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chebaro</surname><given-names>Y</given-names></name><name><surname>Dong</surname><given-names>X</given-names></name><name><surname>Laghaei</surname><given-names>R</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name></person-group><year>2009</year><article-title>Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent.</article-title><source>J Phys Chem B</source><volume>113</volume><fpage>267</fpage><lpage>274</lpage><pub-id pub-id-type="pmid">19067549</pub-id></element-citation></ref><ref id="pcbi.1002051-Carter1"><label>63</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Carter</surname><given-names>P</given-names></name><name><surname>Andersen</surname><given-names>CA</given-names></name><name><surname>Rost</surname><given-names>B</given-names></name></person-group><year>2003</year><article-title>DSSPcont: Continuous secondary structure assignments for proteins.</article-title><source>Nucleic Acids Res</source><volume>31</volume><fpage>3293</fpage><lpage>3295</lpage><pub-id pub-id-type="pmid">12824310</pub-id></element-citation></ref><ref id="pcbi.1002051-Yang1"><label>64</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>L</given-names></name><name><surname>Shao</surname><given-names>Q</given-names></name><name><surname>Goa</surname><given-names>YQ</given-names></name></person-group><year>2009</year><article-title>Comparison between integrated and parallel tempering methods in enhanced sampling simulations.</article-title><source>J Chem Phys</source><volume>130</volume><fpage>124111</fpage><pub-id pub-id-type="pmid">19334812</pub-id></element-citation></ref><ref id="pcbi.1002051-Holm1"><label>65</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Holm</surname><given-names>L</given-names></name><name><surname>Sander</surname><given-names>C</given-names></name></person-group><year>1991</year><article-title>Database algorithm for generating protein backbone and side-chain co-ordinates from a C-alpha trace application to model building and detection of co-ordinate errors.</article-title><source>J Mol Biol</source><volume>218</volume><fpage>183</fpage><lpage>194</lpage><pub-id pub-id-type="pmid">2002501</pub-id></element-citation></ref><ref id="pcbi.1002051-Hartmann1"><label>66</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hartmann</surname><given-names>C</given-names></name><name><surname>Antes</surname><given-names>I</given-names></name><name><surname>Lengauer</surname><given-names>T</given-names></name></person-group><year>2007</year><article-title>IRECS: A new algorithm for the selection of most probable ensembles of side-chain conformations in protein models.</article-title><source>Protein Sci</source><volume>16</volume><fpage>294</fpage><lpage>1307</lpage></element-citation></ref><ref id="pcbi.1002051-Dunbrack1"><label>67</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dunbrack</surname><given-names>RL</given-names><suffix>Jr</suffix></name><name><surname>Cohen</surname><given-names>FE</given-names></name></person-group><year>1997</year><article-title>Bayesian statistical analysis of protein side-chain rotamer preferences.</article-title><source>Protein Sci</source><volume>6</volume><fpage>1661</fpage><lpage>1681</lpage><pub-id pub-id-type="pmid">9260279</pub-id></element-citation></ref><ref id="pcbi.1002051-Berendsen2"><label>68</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berendsen</surname><given-names>HJC</given-names></name><name><surname>Grigera</surname><given-names>JR</given-names></name><name><surname>Straatsma</surname><given-names>TP</given-names></name></person-group><year>1987</year><article-title>The missing term in effective pair potentials.</article-title><source>J Phys Chem</source><volume>91</volume><fpage>6269</fpage><lpage>6271</lpage></element-citation></ref><ref id="pcbi.1002051-Hess2"><label>69</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hess</surname><given-names>B</given-names></name><name><surname>Bekker</surname><given-names>H</given-names></name><name><surname>Berendsen</surname><given-names>HJC</given-names></name><name><surname>Fraaije</surname><given-names>JGEM</given-names></name></person-group><year>1997</year><article-title>LINCS: A linear constraint solver for molecular simulations.</article-title><source>J Comp Chem</source><volume>18</volume><fpage>1463</fpage><lpage>1472</lpage></element-citation></ref><ref id="pcbi.1002051-Miyamoto1"><label>70</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miyamoto</surname><given-names>S</given-names></name><name><surname>Kollman</surname><given-names>PA</given-names></name></person-group><year>1992</year><article-title>SETTLE: An analytical version of the SHAKE and RATTLE algorithms for rigid water models.</article-title><source>J Comp Chem</source><volume>13</volume><fpage>952</fpage><lpage>962</lpage></element-citation></ref><ref id="pcbi.1002051-Berendsen3"><label>71</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berendsen</surname><given-names>HJC</given-names></name><name><surname>van der Spoel</surname><given-names>D</given-names></name><name><surname>van Drunen</surname><given-names>R</given-names></name></person-group><year>1995</year><article-title>GROMACS: A message passing parallel molecular dynamics implementation.</article-title><source>Comp Phys Comm</source><volume>91</volume><fpage>43</fpage><lpage>56</lpage></element-citation></ref><ref id="pcbi.1002051-Lindahl1"><label>72</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lindahl</surname><given-names>E</given-names></name><name><surname>Hess</surname><given-names>B</given-names></name><name><surname>van der Spoel</surname><given-names>D</given-names></name></person-group><year>2001</year><article-title>Gromacs 3.0: A package for molecular simulation and trajectory analysis.</article-title><source>J Mol Model</source><volume>7</volume><fpage>306</fpage><lpage>317</lpage></element-citation></ref><ref id="pcbi.1002051-vanderSpoel1"><label>73</label><element-citation publication-type="other"><person-group person-group-type="author"><name><surname>van der Spoel</surname><given-names>D</given-names></name><name><surname>Lindahl</surname><given-names>E</given-names></name><name><surname>Hess</surname><given-names>B</given-names></name><name><surname>van Buuren</surname><given-names>AR</given-names></name><name><surname>Apol</surname><given-names>E</given-names></name><etal></etal></person-group><year>2004</year><article-title>Gromacs User Manual version 3.2.</article-title><comment>Nijenborgh 4, 9747 AG Groningen, The Netherlands</comment><comment>Available: <ext-link ext-link-type="uri" xlink:href="http://www.gromacs.org">www.gromacs.org</ext-link></comment><size units="page">335</size></element-citation></ref><ref id="pcbi.1002051-vanGunsteren1"><label>74</label><element-citation publication-type="other"><person-group person-group-type="author"><name><surname>van Gunsteren</surname><given-names>WF</given-names></name><name><surname>Berendsen</surname><given-names>HJC</given-names></name></person-group><year>1987</year><article-title>GROMOS-87 manual.</article-title><comment>Nijenborgh 4, 9747 AG Groningen, The Netherlands, Biomos BV</comment></element-citation></ref><ref id="pcbi.1002051-vanGunsteren2"><label>75</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>van Gunsteren</surname><given-names>WF</given-names></name><name><surname>Daura</surname><given-names>WFX</given-names></name><name><surname>Mark</surname><given-names>AE</given-names></name></person-group><year>1998</year><article-title>GROMOS Force Field.</article-title><source>Encyclopedia of Computational Chemistry</source><volume>2</volume><fpage>1211</fpage><lpage>1216</lpage></element-citation></ref><ref id="pcbi.1002051-Scott1"><label>76</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Scott</surname><given-names>WRP</given-names></name><name><surname>Hunenberger</surname><given-names>PH</given-names></name><name><surname>Tironi</surname><given-names>IG</given-names></name><name><surname>Mark</surname><given-names>AE</given-names></name><name><surname>Billeter</surname><given-names>SR</given-names></name><etal></etal></person-group><year>1999</year><article-title>The GROMOS biomolecular simulation program package.</article-title><source>J Phys Chem A</source><volume>103</volume><fpage>3596</fpage><lpage>3607</lpage></element-citation></ref><ref id="pcbi.1002051-vanGunsteren3"><label>77</label><element-citation publication-type="other"><person-group person-group-type="author"><name><surname>van Gunsteren</surname><given-names>WF</given-names></name><name><surname>Billeter</surname><given-names>SR</given-names></name><name><surname>Eising</surname><given-names>AA</given-names></name><name><surname>Huenenberger</surname><given-names>PH</given-names></name><name><surname>Krueger</surname><given-names>P</given-names></name><etal></etal></person-group><year>1996</year><article-title>Biomolecular Simulation: The GROMOS96 Manual and User Guide.</article-title><comment>Switzerland, vdf Hochschulverlag, ETH Zurich</comment></element-citation></ref><ref id="pcbi.1002051-Liu1"><label>78</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>P</given-names></name><name><surname>Friesner</surname><given-names>RA</given-names></name><name><surname>Berne</surname><given-names>BJ</given-names></name></person-group><year>2005</year><article-title>Replica exchange with solute tempering: a method for sampling biological systems in explicit water.</article-title><source>Proc Natl Acad Sci U S A</source><volume>102</volume><fpage>13749</fpage><lpage>13754</lpage><pub-id pub-id-type="pmid">16172406</pub-id></element-citation></ref><ref id="pcbi.1002051-Kubitzki1"><label>79</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kubitzki</surname><given-names>MB</given-names></name><name><surname>de Groot</surname><given-names>BL</given-names></name></person-group><year>2007</year><article-title>Molecular dynamics simulations using temperature-enhanced essential dynamics replica exchange.</article-title><source>Biophys J</source><volume>92</volume><fpage>4262</fpage><lpage>4270</lpage><pub-id pub-id-type="pmid">17384062</pub-id></element-citation></ref><ref id="pcbi.1002051-vanderSpoel2"><label>80</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>van der Spoel</surname><given-names>D</given-names></name><name><surname>Seibert</surname><given-names>MM</given-names></name></person-group><year>2006</year><article-title>Protein folding kinetics and thermodynamics from atomistic simulations.</article-title><source>Phys Rev E</source><volume>96</volume><fpage>238102</fpage></element-citation></ref><ref id="pcbi.1002051-Wei2"><label>81</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wei</surname><given-names>G</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2007</year><article-title>Computational simulations of the early steps of protein aggregation.</article-title><source>Prion</source><volume>1</volume><fpage>3</fpage><lpage>8</lpage><pub-id pub-id-type="pmid">19164927</pub-id></element-citation></ref><ref id="pcbi.1002051-Santini1"><label>82</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Santini</surname><given-names>S</given-names></name><name><surname>Wei</surname><given-names>G</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2004</year><article-title>Pathway complexity of Alzheimer's beta-amyloid Abeta16–22 peptide assembly.</article-title><source>Structure</source><volume>12</volume><fpage>1245</fpage><lpage>1255</lpage><pub-id pub-id-type="pmid">15242601</pub-id></element-citation></ref><ref id="pcbi.1002051-Pellarin1"><label>83</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pellarin</surname><given-names>R</given-names></name><name><surname>Schuetz</surname><given-names>P</given-names></name><name><surname>Guarnera</surname><given-names>E</given-names></name><name><surname>Caflisch</surname><given-names>A</given-names></name></person-group><year>2010</year><article-title>Amyloid fibril polymorphism is under kinetic control.</article-title><source>J Am Chem Soc</source><volume>132</volume><fpage>14960</fpage><lpage>14970</lpage><pub-id pub-id-type="pmid">20923147</pub-id></element-citation></ref><ref id="pcbi.1002051-Santini2"><label>84</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Santini</surname><given-names>S</given-names></name><name><surname>Mousseau</surname><given-names>N</given-names></name><name><surname>Derreumaux</surname><given-names>P</given-names></name></person-group><year>2004</year><article-title>In silico assembly of Alzheimer's Abeta16–22 peptide into beta-sheets.</article-title><source>J Am Chem Soc</source><volume>126</volume><fpage>11509</fpage><lpage>11516</lpage><pub-id pub-id-type="pmid">15366896</pub-id></element-citation></ref><ref id="pcbi.1002051-Li1"><label>85</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>DW</given-names></name><name><surname>Mohanty</surname><given-names>S</given-names></name><name><surname>Irbäck</surname><given-names>A</given-names></name><name><surname>Huo</surname><given-names>S</given-names></name></person-group><year>2008</year><article-title>Formation and growth of oligomers: a Monte Carlo study of an amyloid tau fragment.</article-title><source>PLoS Comput Biol</source><volume>4</volume><fpage>e1000238</fpage><pub-id pub-id-type="pmid">19057640</pub-id></element-citation></ref><ref id="pcbi.1002051-Tycko1"><label>86</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tycko</surname><given-names>R</given-names></name><name><surname>Sciarretta</surname><given-names>KL</given-names></name><name><surname>Orgel</surname><given-names>JP</given-names></name><name><surname>Meredith</surname><given-names>SC</given-names></name></person-group><year>2009</year><article-title>Evidence for novel beta-sheet structures in Iowa mutant beta-amyloid fibrils.</article-title><source>Biochemistry</source><volume>48</volume><fpage>6072</fpage><lpage>6084</lpage><pub-id pub-id-type="pmid">19358576</pub-id></element-citation></ref><ref id="pcbi.1002051-Salemme1"><label>87</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Salemme</surname><given-names>FR</given-names></name><name><surname>Weatherford</surname><given-names>DW</given-names></name></person-group><year>1981</year><article-title>Conformational and geometrical properties of beta-sheets in proteins. I. Parallel beta-sheets.</article-title><source>J Mol Biol</source><volume>146</volume><fpage>101</fpage><lpage>117</lpage><pub-id pub-id-type="pmid">7265225</pub-id></element-citation></ref><ref id="pcbi.1002051-Salemme2"><label>88</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Salemme</surname><given-names>FR</given-names></name><name><surname>Weatherford</surname><given-names>DW</given-names></name></person-group><year>1981</year><article-title>Conformational and geometrical properties of beta-sheets in proteins. II. Antiparallel and mixed beta-sheets.</article-title><source>J Mol Biol</source><volume>146</volume><fpage>119</fpage><lpage>141</lpage><pub-id pub-id-type="pmid">7265226</pub-id></element-citation></ref><ref id="pcbi.1002051-Antzutkin1"><label>89</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Antzutkin</surname><given-names>ON</given-names></name><name><surname>Balbach</surname><given-names>JJ</given-names></name><name><surname>Leapman</surname><given-names>RD</given-names></name><name><surname>Rizzo</surname><given-names>NW</given-names></name><name><surname>Reed</surname><given-names>J</given-names></name><etal></etal></person-group><year>2000</year><article-title>Multiple quantum solid-state NMW indicates a parallel, not antiparallel, organization of beta-sheets in Alzheimer's beta-amyloid fibrils.</article-title><source>Proc Natl Acad Sci U S A</source><volume>97</volume><fpage>13045</fpage><lpage>13050</lpage><pub-id pub-id-type="pmid">11069287</pub-id></element-citation></ref><ref id="pcbi.1002051-Shorter1"><label>90</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shorter</surname><given-names>J</given-names></name><name><surname>Lindquist</surname><given-names>S</given-names></name></person-group><year>2006</year><article-title>Destruction or potentiation of different prions catalyzed by similar Hsp104 remodeling activities.</article-title><source>Mol Cell</source><volume>23</volume><fpage>425</fpage><lpage>438</lpage><pub-id pub-id-type="pmid">16885031</pub-id></element-citation></ref><ref id="pcbi.1002051-Behrends1"><label>91</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Behrends</surname><given-names>C</given-names></name><name><surname>Langer</surname><given-names>CA</given-names></name><name><surname>Boteva</surname><given-names>R</given-names></name><name><surname>Böttcher</surname><given-names>UM</given-names></name><name><surname>Stemp</surname><given-names>MJ</given-names></name><etal></etal></person-group><year>2006</year><article-title>Chaperonin TRiC promotes the assembly of polyQ expansion proteins into nontoxic oligomers.</article-title><source>Mol Cell</source><volume>23</volume><fpage>887</fpage><lpage>97</lpage><pub-id pub-id-type="pmid">16973440</pub-id></element-citation></ref><ref id="pcbi.1002051-Ehrnhoefer1"><label>92</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ehrnhoefer</surname><given-names>DE</given-names></name><name><surname>Bieschke</surname><given-names>J</given-names></name><name><surname>Boeddrich</surname><given-names>A</given-names></name><name><surname>Herbst</surname><given-names>M</given-names></name><name><surname>Masino</surname><given-names>L</given-names></name><etal></etal></person-group><year>2008</year><article-title>EGCG redirects amyloidogenic polypeptides into unstructured, off-pathway oligomers.</article-title><source>Nat Struct Mol Biol</source><volume>15</volume><fpage>558</fpage><lpage>566</lpage><pub-id pub-id-type="pmid">18511942</pub-id></element-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000F92 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000F92 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Sante
|area= MersV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:3098217
|texte= A Multiscale Approach to Characterize the Early Aggregation Steps of the Amyloid-Forming Peptide GNNQQNY from the Yeast Prion Sup-35
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:21625573" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a MersV1
| This area was generated with Dilib version V0.6.33. |  |