Normand Mousseau

EVENT-BASED RELAXATION OF CONTINUOUS DISORDERED SYSTEMS

A computational approach is presented to obtain energy-minimized structures in glassy materials. This approach, the activation-relaxation technique (ART), achieves its efficiency by focusing on significant changes in the microscopic structure (events). The application of ART is illustrated with two examples: the structure of amorphous silicon and the structure of Ni{sub 80}P{sub 20}, a metallic glass. {copyright} {ital 1996 The American Physical Society.}

Amyloid β Protein and Alzheimer’s Disease: When Computer Simulations Complement Experimental Studies

Simulations Complement Experimental Studies Jessica Nasica-Labouze,† Phuong H. Nguyen,† Fabio Sterpone,† Olivia Berthoumieu,‡ Nicolae-Viorel Buchete, Sebastien Cote, Alfonso De Simone, Andrew J. Doig, Peter Faller,‡ Angel Garcia, Alessandro Laio, Mai Suan Li, Simone Melchionna, Normand Mousseau, Yuguang Mu, Anant Paravastu, Samuela Pasquali,† David J. Rosenman, Birgit Strodel, Bogdan Tarus,† John H. Viles, Tong Zhang,†,▲ Chunyu Wang, and Philippe Derreumaux*,†,□ †Laboratoire de Biochimie Theorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Universite Paris Diderot, Sorbonne Paris Cite, 13 rue Pierre et Marie Curie, 75005 Paris, France ‡LCC (Laboratoire de Chimie de Coordination), CNRS, Universite de Toulouse, Universite Paul Sabatier (UPS), Institut National Polytechnique de Toulouse (INPT), 205 route de Narbonne, BP 44099, Toulouse F-31077 Cedex 4, France School of Physics & Complex and Adaptive Systems Laboratory, University College Dublin, Belfield, Dublin 4, Ireland Deṕartement de Physique and Groupe de recherche sur les proteines membranaires (GEPROM), Universite de Montreal, C.P. 6128, succursale Centre-ville, Montreal, Quebec H3C 3T5, Canada Department of Life Sciences, Imperial College London, London SW7 2AZ, United Kingdom Manchester Institute of Biotechnology, University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom Department of Physics, Applied Physics, & Astronomy, and Department of Biology, Rensselaer Polytechnic Institute, Troy, New York 12180, United States The International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland Institute for Computational Science and Technology, SBI Building, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam Instituto Processi Chimico-Fisici, CNR-IPCF, Consiglio Nazionale delle Ricerche, 00185 Roma, Italy School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, 637551 Singapore Department of Chemical and Biomedical Engineering, Florida A&M University-Florida State University (FAMU-FSU) College of Engineering, 2525 Pottsdamer Street, Tallahassee, Florida 32310, United States National High Magnetic Field Laboratory, 1800 East Paul Dirac Drive, Tallahassee, Florida 32310, United States Institute of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Julich GmbH, 52425 Julich, Germany School of Biological and Chemical Sciences, Queen Mary University of London, London E1 4NS, United Kingdom Institut Universitaire de France, 75005 Paris, France

Traveling through potential energy landscapes of disordered materials: The activation-relaxation technique

A detailed description of the activation-relaxation technique (ART) is presented. This method defines events in the configurational energy landscape of disordered materials such as amorphous semiconductors, glasses and polymers, in a two-step process: first, a configuration is activated from a local minimum to a nearby saddle point; next, the configuration is relaxed to a new minimum; this allows for jumps over energy barriers much higher than what can be reached with standard techniques. Such events can serve as basic steps in equilibrium and kinetic Monte Carlo schemes.

High-quality continuous random networks

The continuous random network (CRN) model is an idealized model for perfectly coordinated amorphous semiconductors. The quality of a CRN can be assessed in terms of topological and configurational properties, including coordination, bond-angle distributions, and deformation energy. Using a variation on the sillium approach proposed 14 years ago by Wooten, Winer, and Weaire, we present 1000-atom and 4096-atom configurations with a degree of strain significantly less than the best CRN available at the moment and comparable to experimental results. The low strain is also reflected in the electronic properties. The electronic density of state obtained from ab initio calculation shows a perfect band gap, without any defect, in agreement with experimental data.

Thermodynamics and dynamics of amyloid peptide oligomerization are sequence dependent

Aggregation of the full‐length amyloid‐β (Aβ) and β2‐microglobulin (β2m) proteins is associated with Alzheimers disease and dialysis‐related amyloidosis, respectively. This assembly process is not restricted to full‐length proteins, however, many short peptides also assemble into amyloid fibrils in vitro. Remarkably, the kinetics of amyloid‐fibril formation of all these molecules is generally described by a nucleation‐polymerization process characterized by a lag phase associated with the formation of a nucleus, after which fibril elongation occurs rapidly. In this study, we report using long molecular dynamics simulations with the OPEP coarse‐grained force field, the thermodynamics and dynamics of the octamerization for two amyloid 7‐residue peptides: the β2m83‐89 NHVTLSQ and Aβ16‐22 KLVFFAE fragments. Based on multiple trajectories run at 310 K, totaling 2.2 μs (β2m83‐89) and 4.8 μs (Aβ16‐22) and starting from random configurations and orientations of the chains, we find that the two peptides not only share common but also very different aggregation properties. Notably, an increase in the hydrophobic character of the peptide, as observed in Aβ16‐22 with respect to β2m83‐89 impacts the thermodynamics by reducing the population of bilayer β‐sheet assemblies. Higher hydrophobicity is also found to slow down the dynamics of β‐sheet formation by enhancing the averaged lifetime of all configuration types (CT) and by reducing the complexity of the CT transition probability matrix. Proteins 2009.

Coarse-grained protein molecular dynamics simulations

A limiting factor in biological science is the time-scale gap between experimental and computational trajectories. At this point, all-atom explicit solvent molecular dynamics (MD) are clearly too expensive to explore long-range protein motions and extract accurate thermodynamics of proteins in isolated or multimeric forms. To reach the appropriate time scale, we must then resort to coarse graining. Here we couple the coarse-grained OPEP model, which has already been used with activated methods, to MD simulations. Two test cases are studied: the stability of three proteins around their experimental structures and the aggregation mechanisms of the Alzheimers Abeta16-22 peptides. We find that coarse-grained isolated proteins are stable at room temperature within 50 ns time scale. Based on two 220 ns trajectories starting from disordered chains, we find that four Abeta16-22 peptides can form a three-stranded beta sheet. We also demonstrate that the reptation move of one chain over the others, first observed using the activation-relaxation technique, is a kinetically important mechanism during aggregation. These results show that MD-OPEP is a particularly appropriate tool to study qualitatively the dynamics of long biological processes and the thermodynamics of molecular assemblies.

Fitting the Stillinger–Weber potential to amorphous silicon

Abstract Modifications are proposed to the Stillinger–Weber (SW) potential, an empirical interaction potential for silicon. The modifications are specifically intended to improve the description of the amorphous phase and are obtained by a direct fit to the amorphous structure. The potential is adjusted to reproduce the location of the transverse optic (TO) and transverse acoustic (TA) peaks of the vibrational density of states (VDOS), properties insensitive to the details of experimental preparation. These modifications also lead to excellent agreement with structural properties. Comparison with other empirical potentials shows that amorphous silicon configurations generated with the modified potential have overall better vibrational and structural properties.

Role of the Region 23-28 in Aβ Fibril Formation: Insights from Simulations of the Monomers and Dimers of Alzheimers Peptides Aβ40 and Aβ42

Self-assembly of the 40/42 amino acid A! peptide is a key player in Alzheimers disease. A! 40 is the most prevalent species, while A! 42 is the most toxic. It has been suggested that the amino acids 21-30 could nucleate the fold- ing of A! monomer and a bent in this region could be the rate-limiting step in A! fibril formation. In this study, we review our current understanding of the computer-predicted conformations of amino acids 23-28 in the monomer of A! (21-30) and the monomers A! 40 and A! 42. On the basis of new simulations on dimers of full-length A! , we propose that the rate- limiting step involves the formation of a multimeric ! -sheet spanning the central hydrophobic core (residues 17-21).

Structures of Aβ17–42 trimers in isolation and with five small-molecule drugs using a hierarchical computational procedure

The amyloid-β protein (Aβ) oligomers are believed to be the main culprits in the cytoxicity of Alzheimers disease (AD) and p3 peptides (Aβ17-42 fragments) are present in AD amyloid plaques. Many small-molecule or peptide-based inhibitors are known to slow down Aβ aggregation and reduce the toxicity in vitro, but their exact modes of action remain to be determined since there has been no atomic level of Aβ(p3)-drug oligomers. In this study, we have determined the structure of Aβ17-42 trimers both in aqueous solution and in the presence of five small-molecule inhibitors using a multiscale computational study. These inhibitors include 2002-H20, curcumin, EGCG, Nqtrp, and resveratrol. First, we used replica exchange molecular dynamics simulations coupled to the coarse-grained (CG) OPEP force field. These CG simulations reveal that the conformational ensemble of Aβ17-42 trimer can be described by 14 clusters with each peptide essentially adopting turn/random coil configurations, although the most populated cluster is characterized by one peptide with a β-hairpin at Phe19-Leu31. Second, these 14 dominant clusters and the less-frequent fibril-like state with parallel register of the peptides were subjected to atomistic Autodock simulations. Our analysis reveals that the drugs have multiple binding modes with different binding affinities for trimeric Aβ17-42 although they interact preferentially with the CHC region (residues 17-21). The compounds 2002-H20 and Nqtrp are found to be the worst and best binders, respectively, suggesting that the drugs may interfere at different stages of Aβ oligomerization. Finally, explicit solvent molecular dynamics of two predicted Nqtrp-Aβ17-42 conformations describe at atomic level some possible modes of action for Nqtrp.

Complex folding pathways in a simple β‐hairpin

The determination of the folding mechanisms of proteins is critical to understand the topological change that can propagate Alzheimer and Creutzfeld‐Jakobs diseases, among others. The computational community has paid considerable attention to this problem; however, the associated time scale, typically on the order of milliseconds or more, represents a formidable challenge. Ab initio protein folding from long molecular dynamics simulations or ensemble dynamics is not feasible with ordinary computing facilities and new techniques must be introduced. Here we present a detailed study of the folding of a 16‐residue β‐hairpin, described by a generic energy model and using the activation‐relaxation technique. From a total of 90 trajectories at 300 K, three folding pathways emerge. All involve a simultaneous optimization of the complete hydrophobic and hydrogen bonding interactions. The first two pathways follow closely those observed by previous theoretical studies (folding starting at the turn or by interactions between the termini). The third pathway, never observed by previous all‐atom folding, unfolding, and equilibrium simulations, can be described as a reptation move of one strand of the β‐sheet with respect to the other. This reptation move indicates that non‐native interactions can play a dominant role in the folding of secondary structures. Furthermore, such a mechanism mediated by non‐native hydrogen bonds is not available for study by unfolding and Gō model simulations. The exact folding path followed by a given β‐hairpin is likely to be influenced by its sequence and the solvent conditions. Taken together, these results point to a more complex folding picture than expected for a simple β‐hairpin. Proteins 2004;9999:000–000.