Download (7Mb)

Total Page:16

File Type:pdf, Size:1020Kb

Download (7Mb) University of Warwick institutional repository: http://go.warwick.ac.uk/wrap A Thesis Submitted for the Degree of PhD at the University of Warwick http://go.warwick.ac.uk/wrap/74165 This thesis is made available online and is protected by original copyright. Please scroll down to view the document itself. Please refer to the repository record for this item for information to help you to cite it. Our policy information is available from the repository home page. MENS T A T A G I MOLEM U N IS IV S E EN RS IC ITAS WARW Coarse-Grained Simulations of Intrinsically Disordered Peptides by Gil Rutter Thesis Submitted to the University of Warwick for the degree of Doctor of Philosophy Department of Physics October 2015 Contents List of Tables iv List of Figures vi Acknowledgments xvii Declarations xviii Abstract xix Chapter 1 Introduction 1 1.1 Proteinstructure ............................. 1 1.1.1 Dependence on ambient conditions . 3 1.2 Proteindisorder.............................. 6 1.2.1 Determinants of protein disorder . 6 1.2.2 Functions of IDPs . 7 1.3 Experimentaltechniques. 8 1.3.1 Intrinsicdisorder ......................... 8 1.3.2 Biomineralisation systems . 10 1.4 Molecular dynamics for proteins . 12 1.4.1 Approaches to molecular simulation . 12 1.4.2 Statistical ensembles . 13 1.5 Biomineralisation . 15 1.6 n16N.................................... 16 1.7 Summary ................................. 20 Chapter 2 Accelerated simulation techniques 21 2.1 Coarse-graining . 21 2.1.1 DefiningtheCGmapping . 22 2.1.2 Explicit and implicit solvation . 24 i 2.1.3 Validating a CG model . 26 2.1.4 Choices of coarse-grained models . 26 2.2 Acceleratedsamplingtechniques . 29 2.2.1 Replica exchange molecular dynamics . 30 2.2.2 Statistical temperature molecular dynamics . 31 2.3 Summary ................................. 34 Chapter 3 Software modifications and validation 35 3.1 PRIME20 initial parametrisation . 35 3.1.1 PRIME model validation . 36 3.1.2 Parametrising a PRIME20-like model . 40 3.2 Implementation work in LAMMPS . 42 3.2.1 The PLUM model . 42 3.2.2 Statistical temperature molecular dynamics . 46 3.3 Summary ................................. 52 Chapter 4 Single-chain simulations 53 4.1 PLUM model simulations . 53 4.1.1 Over-stabilisation of the ↵-helix . 53 4.1.2 PLUM* validation work . 54 4.1.3 The n16N-1 system in PLUM* . 58 4.1.4 The n16NN-1 system . 61 4.1.5 TheS1system .......................... 62 4.2 PRIME20-like model simulations . 64 4.2.1 TheS1system .......................... 64 4.2.2 The n16N-1 system . 65 4.2.3 The n16NN-1 system . 70 4.3 Summary ................................. 71 Chapter 5 Multiple-chain simulations 73 5.1 The n16N-2 system . 73 5.1.1 PLUM* . 73 5.1.2 PLUM . 74 5.1.3 PRIME20-like . 78 5.2 The n16NN-2 system . 81 5.2.1 PLUM* . 81 5.2.2 PRIME20-like . 82 5.3 Discussionofdimersystems . 83 ii 5.4 The n16N-3 system . 84 5.4.1 PLUM* . 84 5.4.2 PRIME20-like . 85 5.5 The n16NN-3 system . 87 5.6 Discussionoftrimersystems. 90 5.7 The n16N-6 system . 91 5.7.1 PLUM* . 91 5.7.2 PRIME20-like . 92 5.8 The n16NN-6 system . 97 5.9 Discussion of hexamer systems . 99 5.10Summary ................................. 101 Chapter 6 Conclusions 103 6.1 PLUM*model .............................. 103 6.1.1 Simulation results . 103 6.1.2 Furtherwork ........................... 105 6.2 ThePRIME20-likemodel ........................ 106 6.2.1 Simulation results . 106 6.2.2 Furtherwork ........................... 107 6.3 Outlook for coarse-grained models to study IDPs . 108 iii List of Tables 1.1 Relative frequencies of residues in common secondary structure mo- tifs [Creighton, 1992] [Berg et al., 2002, page 67] and disorder propen- sity according to the TOP-IDP scale [Campen et al., 2008], in which higher values are more disordered. This is one of many scales on which residue disorder propensity has been ranked [Galzitskaya et al., 2006; Oldfield et al., 2005; Vihinen et al., 1994; Garbuzynskiy et al., 2004]. The highest values are highlighted. .............. 7 1.2 Summary of experimental research on n16N peptide, with a focus on experimental techniques used. All assays contain n16N in water solution with calcium and carbonate ions for crystal growth. AFM - Atomic force microscopy CD - Circular dichroism spectroscopy DLS - Dynamic light scattering EDX - Energy dispersive X-ray spectroscopy FM - Fluorescence microscopy FTIR - Fourier-transformed infrared spectroscopy Raman - Raman spectroscopy SEM - Scanning electron microscopy TEM - Tunnelling electron microscopy XANES - X-ray absorption near edge spectromicroscopy XRD - X-ray di↵raction X-PEEM - X-ray photoelectron emission spectromicroscopy .... 11 1.3 Suggested roles of the subdomains of n16N [Brown et al., 2014]. .. 20 iv 4.1 Top interactions between residues ranked by frequency of occurrence. To qualify, any atom on residue A has to be in the square well of any atom on residue B, with A and B separated with at least three residues in-between. Tyrosine dominates the rankings, and is clearly a cornerstone of intrapeptide stabilisation. .............. 69 v List of Figures 1.1 Wikimedia user Dcrjsr, available under the Creative Commons Attri- bution 3.0 Unported license. The backbone atoms of a central amino acid residue is shown, with two neighbours labelled (-1) and (+1) partially shown. Each residue except glycine also has a side-chain, denoted Cβ. The dihedral angles φ, and ! which determine the secondary structure are labelled. 2 1.2 [Hollingsworth and Karplus, 2010] adapted from [Ramachandran and Sasisekharan, 1968]. A Ramachandran plot showing both example data and outlines of commonly accepted regions. The data is comprised of 63,149 residues from crystal structures and each point indicates the (φ, ) values rep- resenting the secondary structure of a single residue. The outlines are divided into core allowed regions (solid lines) and allowed regions (dashed lines). Several forms of secondary structure have their location indicated; these are ↵-helix (↵), 310-helix (3), ⇡- helix (⇡), left-handed ↵-helix (↵l), 2.27 ribbon (2), polyproline-II (II), collagen (C), parallel β-sheet ( ) and anti-parallel β-sheet ( ). 4 "" "# vi 1.3 A steric map, in which steric clashes leading to forbidden regions are shown in (φ, ) phase space and labelled according to the clashing pair of atoms [Ho et al., 2003]. Dashed blue lines and blue labels denote clashes. Favourable dipole-dipole backbone interactions are also plot- ted, with red dashed lines and red labels. Areas of the map are white if forbidden, and otherwise coloured according to the legend below. The atom labels match the illustration in fig. 1.1. Boundaries ---Attractive dipole-dipole interactions ---Steric clashes Accessible regions ⌅⌅⌅ Sterically permitted ⌅⌅⌅ Single steric clash (outlier region) ⌅⌅⌅ Left- and right-handed ↵-helix regions ⌅⌅⌅ β-strand region .......................... 5 1.4 [Levi-Kalisman et al., 2001] A putative scheme for the organic biomineralisation matrix which grows nacre. Interlamellar matrix sheets are composed mainly of aligned β-chitin fibres in several layers, with acidic glycoproteins at their surface which lead to electron-dense patches (not labelled). Arag- onite biomineralisation occurs in the disordered silk fibroin gel region, most likely nucleating epitaxially on the chitin framework [Weiner et al., 1984]. .............................. 17 1.5 Amino acid sequence of the 30AA N-terminal region of n16, called n16N. An ellipsis indicates where the full n16 sequence continues, and braces indicate suggested subdomains [Brown et al., 2014], sum- marised in table 1.3. Cationic amino acid residues shown in bold blue, anionic residues shown in bold red. The last 14 residues, labelled SD3, represent a highly charged region which may be the mineral assembly subdomain. ............................... 18 vii 1.6 [Amos et al., 2011] Two-stage mineralisation experiment to test whether pre-formed n16N assemblies could nucleate calcium carbonate polymorphs. In stage one, n16N deposits are allowed to form in typical mineralisation con- ditions (50 µMn16N,16h,16◦C). In stage two, washed supports with n16N deposits are transferred to a mineralisation solution with or without n16N content. (A) Calcium carbonate mineralisation without n16N IDP. Calcite polymorph dominates. (B) n16N present during first stage of mineralisation, yielding more vaterite (v). (C,D) n16N present during both stages of mineralisation. The ar- rows indicate fibril-spheroidal deposits of n16N on the exposed surface of a crystal. ............................... 19 2.1 The PRIME20 and PLUM models are identical in their mappings of atoms to coarse-grained sites. All residues feature a single side-chain site, except for glycine which has no side-chain site. The PLUM model backbone maintains its spatial arrangement through continuous- potential bonds, angles and dihedrals as in equation (2.2),whilethe PRIME20 model relies upon square bonds and square pseudo-bonds, seen in fig. 2.1b. ............................. 27 3.1 Ramachandran plots showing the exploration of (φ, ) phase space for PRIME2001 chains, according to the reference [Voegler Smith and Hall, 2001], our independent reproduction, and a reproduction with a modified parameter set. Note that T =0.15 is a high simulation temperature, certainly above that used by the original authors. By modifying the parameter set, it is possible to reproduce the reference data’s accessible regions with moderate accuracy. ........... 37 3.2 Behaviour of a chain of A20 in PRIME2001, PRIME2001-like and PRIME20-like models. After altering the 2001 model to reproduce the authors’ accessible areas (see fig. 3.1), we see here that the hydrogen bond interaction strength, ✏HB,needsraisingto1.26 to match the reference data best. In the PRIME20-like model, ↵-helices are far more stable. Unfortunately, no data for the behaviour of A20 in the canonical PRIME20 model is available. ................ 39 viii 3.3 Ramachandran plots in the custom PRIME2001-like model at T = 0.125,with✏HB =1.26. These plots illustrate that the PRIME2001 data in fig.
Recommended publications
  • Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development
    processes Review Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development Outi M. H. Salo-Ahen 1,2,* , Ida Alanko 1,2, Rajendra Bhadane 1,2 , Alexandre M. J. J. Bonvin 3,* , Rodrigo Vargas Honorato 3, Shakhawath Hossain 4 , André H. Juffer 5 , Aleksei Kabedev 4, Maija Lahtela-Kakkonen 6, Anders Støttrup Larsen 7, Eveline Lescrinier 8 , Parthiban Marimuthu 1,2 , Muhammad Usman Mirza 8 , Ghulam Mustafa 9, Ariane Nunes-Alves 10,11,* , Tatu Pantsar 6,12, Atefeh Saadabadi 1,2 , Kalaimathy Singaravelu 13 and Michiel Vanmeert 8 1 Pharmaceutical Sciences Laboratory (Pharmacy), Åbo Akademi University, Tykistökatu 6 A, Biocity, FI-20520 Turku, Finland; ida.alanko@abo.fi (I.A.); rajendra.bhadane@abo.fi (R.B.); parthiban.marimuthu@abo.fi (P.M.); atefeh.saadabadi@abo.fi (A.S.) 2 Structural Bioinformatics Laboratory (Biochemistry), Åbo Akademi University, Tykistökatu 6 A, Biocity, FI-20520 Turku, Finland 3 Faculty of Science-Chemistry, Bijvoet Center for Biomolecular Research, Utrecht University, 3584 CH Utrecht, The Netherlands; [email protected] 4 Swedish Drug Delivery Forum (SDDF), Department of Pharmacy, Uppsala Biomedical Center, Uppsala University, 751 23 Uppsala, Sweden; [email protected] (S.H.); [email protected] (A.K.) 5 Biocenter Oulu & Faculty of Biochemistry and Molecular Medicine, University of Oulu, Aapistie 7 A, FI-90014 Oulu, Finland; andre.juffer@oulu.fi 6 School of Pharmacy, University of Eastern Finland, FI-70210 Kuopio, Finland; maija.lahtela-kakkonen@uef.fi (M.L.-K.); tatu.pantsar@uef.fi
    [Show full text]
  • Density Functional Theory for Protein Transfer Free Energy
    Density Functional Theory for Protein Transfer Free Energy Eric A Mills, and Steven S Plotkin∗ Department of Physics & Astronomy, University of British Columbia, Vancouver, British Columbia V6T1Z4 Canada E-mail: [email protected] Phone: 1 604-822-8813. Fax: 1 604-822-5324 arXiv:1310.1126v1 [q-bio.BM] 3 Oct 2013 ∗To whom correspondence should be addressed 1 Abstract We cast the problem of protein transfer free energy within the formalism of density func- tional theory (DFT), treating the protein as a source of external potential that acts upon the sol- vent. Solvent excluded volume, solvent-accessible surface area, and temperature-dependence of the transfer free energy all emerge naturally within this formalism, and may be compared with simplified “back of the envelope” models, which are also developed here. Depletion contributions to osmolyte induced stability range from 5-10kBT for typical protein lengths. The general DFT transfer theory developed here may be simplified to reproduce a Langmuir isotherm condensation mechanism on the protein surface in the limits of short-ranged inter- actions, and dilute solute. Extending the equation of state to higher solute densities results in non-monotonic behavior of the free energy driving protein or polymer collapse. Effective inter- action potentials between protein backbone or sidechains and TMAO are obtained, assuming a simple backbone/sidechain 2-bead model for the protein with an effective 6-12 potential with the osmolyte. The transfer free energy dg shows significant entropy: d(dg)=dT ≈ 20kB for a 100 residue protein. The application of DFT to effective solvent forces for use in implicit- solvent molecular dynamics is also developed.
    [Show full text]
  • Computational Redox Potential Predictions: Applications to Inorganic and Organic Aqueous Complexes, and Complexes Adsorbed to Mineral Surfaces
    Minerals 2014, 4, 345-387; doi:10.3390/min4020345 OPEN ACCESS minerals ISSN 2075-163X www.mdpi.com/journal/minerals Review Computational Redox Potential Predictions: Applications to Inorganic and Organic Aqueous Complexes, and Complexes Adsorbed to Mineral Surfaces Krishnamoorthy Arumugam and Udo Becker * Department of Earth and Environmental Sciences, University of Michigan, 1100 North University Avenue, 2534 C.C. Little, Ann Arbor, MI 48109-1005, USA; E-Mail: [email protected] * Author to whom correspondence should be addressed; E-Mail: [email protected]; Tel.: +1-734-615-6894; Fax: +1-734-763-4690. Received: 11 February 2014; in revised form: 3 April 2014 / Accepted: 13 April 2014 / Published: 24 April 2014 Abstract: Applications of redox processes range over a number of scientific fields. This review article summarizes the theory behind the calculation of redox potentials in solution for species such as organic compounds, inorganic complexes, actinides, battery materials, and mineral surface-bound-species. Different computational approaches to predict and determine redox potentials of electron transitions are discussed along with their respective pros and cons for the prediction of redox potentials. Subsequently, recommendations are made for certain necessary computational settings required for accurate calculation of redox potentials. This article reviews the importance of computational parameters, such as basis sets, density functional theory (DFT) functionals, and relativistic approaches and the role that physicochemical processes play on the shift of redox potentials, such as hydration or spin orbit coupling, and will aid in finding suitable combinations of approaches for different chemical and geochemical applications. Identifying cost-effective and credible computational approaches is essential to benchmark redox potential calculations against experiments.
    [Show full text]
  • Comparative Study of Implicit and Explicit Solvation Models for Probing Tryptophan Side Chain Packing in Proteins†
    828 Bull. Korean Chem. Soc. 2012, Vol. 33, No. 3 Changwon Yang and Youngshang Pak http://dx.doi.org/10.5012/bkcs.2012.33.3.828 Comparative Study of Implicit and Explicit Solvation Models for Probing Tryptophan Side Chain Packing in Proteins† Changwon Yang and Youngshang Pak* Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 609-735, Korea *E-mail: [email protected] Received November 10, 2011, Accepted December 29, 2011 We performed replica exchange molecular dynamics (REMD) simulations of the tripzip2 peptide (beta- hairpin) using the GB implicit and TI3P explicit solvation models. By comparing the resulting free energy surfaces of these two solvation model, we found that the GB solvation model produced a distorted free energy map, but the explicit solvation model yielded a reasonable free energy landscape with a precise location of the native structure in its global free energy minimum state. Our result showed that in particular, the GB solvation model failed to describe the tryptophan packing of trpzip2, leading to a distorted free energy landscape. When the GB solvation model is replaced with the explicit solvation model, the distortion of free energy shape disappears with the native-like structure in the lowest free energy minimum state and the experimentally observed tryptophan packing is precisely recovered. This finding indicates that the main source of this problem is due to artifact of the GB solvation model. Therefore, further efforts to refine this model are needed for better predictions of various aromatic side chain packing forms in proteins. Key Words : Replica exchange molecular dynamics simulation, Implicit solvation model, Tryptophan side chain packing, Protein folding Introduction develop a transferable all-atom protein force field with the GBSA model, so that free energy based native structure Developments of efficient simulation strategies in con- predictions of diverse structural motifs become feasible at junction with all-atom force fields have provided an the all-atom resolution.
    [Show full text]
  • Atomistic Modeling of DNA and Protein Structures
    Atomistic Modeling of DNA and Protein Structures George C. Schatz Introduction In this lecture we introduce one of the most commonly used computational tools for studying protein and DNA structure, the use of empirical potential energy functions, and we describe their use in molecular mechanics and molecular dynamics calculations. As an application, we present results of a molecular mechanics study of the mechanical properties (mechanical pulling) of a small protein. We also describe modeling of DNA hairpin structures. References: Computer Simulation of Liquids, M. P. Allen and D. J. Tildesley, Clarendon Press, Oxford, 1987; Introduction to Modern Statistical Mechanics, D. Chandler, Oxford, New York, 1987; Understanding Molecular Simulations: from Algorithms to Applications, D. Frenkel and B. Smit, Academic Press; San Diego, 1996. Molecular Modeling: Principles and Applications, (2nd ed) Andrew Leach, Prentice Hall, Englewood Cliffs, 2001. What is a force field? This is the potential energy function V(R) that determines the interactions between the atoms in a molecule or solid. The derivatives of this function (or more technically the gradients) give the forces that these atoms exert against each other. Thus Newton’s equations say that F=ma, where the force F is related to V via F=-dV/dR, m is the mass and a is the acceleration. A simple example of a force field is a Lennard-Jones, or 6-12 potential. This describes the interaction of two argon atoms pretty well, and it is often used to describe the interaction of two methane molecules or even, to a lesser degree of rigor, two water molecules. Lennard-Jones potential 7 6 5 12 6 ⎛⎞⎡⎤σσ ⎡⎤ 4 V4=ε − 3 ⎜⎟⎢⎥ ⎢⎥ RR energy/well depth ⎝⎠⎣⎦ ⎣⎦ 2 1 0 -1 0.80 1.00 1.20 1.40 1.60 1.80 2.00 r/σ For a liquid composed of argon, or methane, one can approximate that the overall potential is simply the sum of all pair-wise interactions: ⎛⎞12 6 ⎡ σσij⎤⎡⎤ ij V4=ε⎜⎟⎢ ⎥⎢⎥ − ∑ ij ⎜⎟RR ij ⎝⎠⎣⎢ ij⎦⎣⎦⎥⎢⎥ ij This neglects three-body and higher interactions, which works pretty well for argon, but is a problem for more complex materials.
    [Show full text]
  • Simulations of the Role of Water in the Protein- Folding Mechanism
    Simulations of the role of water in the protein- folding mechanism Young Min Rhee*, Eric J. Sorin*, Guha Jayachandran†, Erik Lindahl‡, and Vijay S. Pande*‡§ Departments of *Chemistry, †Computer Science, and ‡Structural Biology, Stanford University, Stanford, CA 94305 Edited by Harold A. Scheraga, Cornell University, Ithaca, NY, and approved March 15, 2004 (received for review November 26, 2003) There are many unresolved questions regarding the role of water ments. To address this question, it is natural to look to computer in protein folding. Does water merely induce hydrophobic forces, simulations of protein folding by using explicit solvation models or does the discrete nature of water play a structural role in for direct comparisons to both experimental and implicit solva- folding? Are the nonadditive aspects of water important in deter- tion simulation results. mining the folding mechanism? To help to address these questions, However, even with advances in computational methodologies we have performed simulations of the folding of a model protein in recent years, a folding kinetics simulation in atomistic detail (BBA5) in explicit solvent. Starting 10,000 independent trajectories remains a demanding task (4), and there have been only a few from a fully unfolded conformation, we have observed numerous reports of such studies (5, 6). Moreover, these studies have been folding events, making this work a comprehensive study of the confined to the use of implicit and indirect models of solvation kinetics of protein folding starting from the unfolded state and to reduce the required central processing unit (CPU) time and reaching the folded state and with an explicit solvation model and storage requirements (5).
    [Show full text]
  • Simulations of the Role of Water in the Protein- Folding Mechanism
    Simulations of the role of water in the protein- folding mechanism Young Min Rhee*, Eric J. Sorin*, Guha Jayachandran†, Erik Lindahl‡, and Vijay S. Pande*‡§ Departments of *Chemistry, †Computer Science, and ‡Structural Biology, Stanford University, Stanford, CA 94305 Edited by Harold A. Scheraga, Cornell University, Ithaca, NY, and approved March 15, 2004 (received for review November 26, 2003) There are many unresolved questions regarding the role of water ments. To address this question, it is natural to look to computer in protein folding. Does water merely induce hydrophobic forces, simulations of protein folding by using explicit solvation models or does the discrete nature of water play a structural role in for direct comparisons to both experimental and implicit solva- folding? Are the nonadditive aspects of water important in deter- tion simulation results. mining the folding mechanism? To help to address these questions, However, even with advances in computational methodologies we have performed simulations of the folding of a model protein in recent years, a folding kinetics simulation in atomistic detail (BBA5) in explicit solvent. Starting 10,000 independent trajectories remains a demanding task (4), and there have been only a few from a fully unfolded conformation, we have observed numerous reports of such studies (5, 6). Moreover, these studies have been folding events, making this work a comprehensive study of the confined to the use of implicit and indirect models of solvation kinetics of protein folding starting from the unfolded state and to reduce the required central processing unit (CPU) time and reaching the folded state and with an explicit solvation model and storage requirements (5).
    [Show full text]
  • Solvation in the Cramer-Truhlar Groups
    Solvation in the Cramer-Truhlar Groups Developing solvation models that permit the accurate study of dynamical processes in Adam Chamberlin Natalie Elmasry aqueous and nonaqueous solution Casey Kelly Alek Marenich Ryan Olson Jonathan Smith* February 2007 Chris Cramer Don Truhlar *Gustavus Adolphus College TheThe process process of of moving moving a a moleculemolecule from from the the gas gas phase phase to to What is solvation? aa condensed condensed phase. phase. The solute charge density is AA solute solute shaped shaped cavity cavity of of vacuum vacuum The solute charge density is placed in the solute cavity. isis introduced introduced into into the the solvent. solvent. placed in the solute cavity. TheThe solvent solvent molecules The solute molecules The solute reorientreorient and and polarizespolarizes in in polarize in response to polarize in response to responseresponse to to the the thethe solvent solvent solute charge polarization solute charge polarization densitydensity Objective To develop a simple and general model to predict the solubility of high-energy compounds in binary supercritical solvents Milestones 1. Develop a simple and general solvent model to predict solvent properties at 298 K for a broad variety of compounds including ions in a broad variety of solvents 2. Extend the model to apply to a broad range of temperatures Work In Progress 1. Extend our temperature dependent solvation model (SM6T) to handle N-containing compounds and multiple solvents 2. Modify the model to work with binary solvents 3. Extend the model to apply to a range of pressures 4. Improve the accuracy of the model for ionic compounds in non-aqueous solvents tructure ispersion CDSavitation (SM6) G 6 + olarization uclear odel lectronic , 1133 ENP term is a parameterized 1 M G , Non-bulk electrostaticCDS contributions to the free energyG of hydration: first 2005 solvation shell effects =Δ The o term intended to minimize the S deviation between the predictions G and experiment.
    [Show full text]
  • Advances in Molecular Quantum Chemistry Contained in the Q-Chem 4 Program Package Yihan Shao Q-Chem Inc
    Chemistry Publications Chemistry 2015 Advances in Molecular Quantum Chemistry Contained in the Q-Chem 4 Program Package Yihan Shao Q-Chem Inc. Zhengting Gan Q-Chem Inc. Evgeny Epifanovsky Q-Chem Inc. Andrew T. B. Gilbert Australian National University Michael Wormit Ruprecht-Karls University See next page for additional authors Follow this and additional works at: http://lib.dr.iastate.edu/chem_pubs Part of the Chemistry Commons The ompc lete bibliographic information for this item can be found at http://lib.dr.iastate.edu/ chem_pubs/640. For information on how to cite this item, please visit http://lib.dr.iastate.edu/ howtocite.html. This Article is brought to you for free and open access by the Chemistry at Iowa State University Digital Repository. It has been accepted for inclusion in Chemistry Publications by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Advances in Molecular Quantum Chemistry Contained in the Q-Chem 4 Program Package Abstract A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open- shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given.
    [Show full text]
  • Sampling of the Conformational Landscape of Small Proteins With
    www.nature.com/scientificreports OPEN Sampling of the conformational landscape of small proteins with Monte Carlo methods Nana Heilmann, Moritz Wolf, Mariana Kozlowska, Elaheh Sedghamiz, Julia Setzler, Martin Brieg & Wolfgang Wenzel* Computer simulation provides an increasingly realistic picture of large-scale conformational change of proteins, but investigations remain fundamentally constrained by the femtosecond timestep of molecular dynamics simulations. For this reason, many biologically interesting questions cannot be addressed using accessible state-of-the-art computational resources. Here, we report the development of an all-atom Monte Carlo approach that permits the modelling of the large- scale conformational change of proteins using standard of-the-shelf computational hardware and standard all-atom force felds. We demonstrate extensive thermodynamic characterization of the folding process of the α-helical Trp-cage, the Villin headpiece and the β-sheet WW-domain. We fully characterize the free energy landscape, transition states, energy barriers between diferent states, and the per-residue stability of individual amino acids over a wide temperature range. We demonstrate that a state-of-the-art intramolecular force feld can be combined with an implicit solvent model to obtain a high quality of the folded structures and also discuss limitations that still remain. Conformational change is one of the most prominent mechanisms of protein function and regulation1–4. Te stability of the protein as a whole, as well as the stability of its particular conformational sub-ensembles is essential to understand and regulate protein function. Due to its dynamic nature, it is difcult to observe conformational changes directly at the single molecule level in experiment 5,6, while computational methods, generally, struggle with the timescales involved5,7–10.
    [Show full text]
  • Mathematical Methods for Implicit Solvation Models in Quantum Chemistry Chaoyu Quan
    Mathematical methods for implicit solvation models in quantum chemistry Chaoyu Quan To cite this version: Chaoyu Quan. Mathematical methods for implicit solvation models in quantum chemistry. Mathemat- ical Physics [math-ph]. Université Pierre et Marie Curie - Paris VI; Rheinisch-westfälische technische Hochschule (Aix-la-Chapelle, Allemagne), 2017. English. NNT : 2017PA066587. tel-01919793 HAL Id: tel-01919793 https://tel.archives-ouvertes.fr/tel-01919793 Submitted on 12 Nov 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Université Pierre et Marie Curie Laboratoire Jacques-Louis Lions Paris 6 UMR 7598 Thèse de Doctorat de l’Université Pierre et Marie Curie Présentée et soutenue publiquement le 21 Novembre 2017 pour l’obtention du grade de Docteur de l’Université Pierre et Marie Curie Spécialité : Mathématiques Appliquées par Chaoyu QUAN sous la direction de Yvon MADAY et Benjamin STAMM Mathematical Methods for Implicit Solvation Models in Quantum Chemistry après avis des rapporteurs M. Martin J. GANDER & M. Aihui ZHOU devant le jury composé de M. Eric CANCES Examinateur M. Pascal FREY Examinateur M. Martin J. GANDER Rapporteur Mme. Laura GRIGORI Examinateur M. Yvon MADAY Directeur de thèse M.
    [Show full text]
  • Atomistic Protein Folding Simulations on the Submillisecond Time Scale
    Vijay S. Pande1–4 Ian Baker1 Atomistic Protein Folding Jarrod Chapman4,* Sidney P. Elmer1 Simulations on the Siraj Khaliq5 Submillisecond Time Scale Stefan M. Larson2 Using Worldwide Distributed Young Min Rhee1 Michael R. Shirts1 Computing Christopher D. Snow2 Eric J. Sorin1 Bojan Zagrovic2 1 Department of Chemistry, 4 Stanford Synchrotron Stanford University, Radiation Laboratory, Stanford, CA 94305-5080 Stanford University, Stanford, CA 94305-5080 2 Biophysics Program, Stanford University, 5 Department of Computer Stanford, CA 94305-5080 Science, Stanford University, 3 Department of Structural Stanford, CA 94305-5080 Biology, Received 27 March 2002; Stanford University, accepted 22 May 2002 Stanford, CA 94305-5080 Abstract: Atomistic simulations of protein folding have the potential to be a great complement to experimental studies, but have been severely limited by the time scales accessible with current computer hardware and algorithms. By employing a worldwide distributed computing network of tens of thousands of PCs and algorithms designed to efficiently utilize this new many-processor, highly heterogeneous, loosely coupled distributed computing paradigm, we have been able to simulate hundreds of microseconds of atomistic molecular dynamics. This has allowed us to directly Correspondence to: Vijay S. Pande; email: pande@stanford. edu *Present address: Physics Department, University of California, Berkeley, CA Contact grant sponsor: ACS PRF, NSF MRSEC CPIMA, NIH BISTI, ARO, and Stanford University Contact grant numbers: 36028-AC4 (ACS PRF), DMR-9808677 (NSF MRSEC CPIMA), IP20 6M64782-01 (NIH BISTI), 41778- LS-RIP (ARO) Biopolymers, Vol. 68, 91–109 (2003) © 2002 Wiley Periodicals, Inc. 91 92 Pande et al. simulate the folding mechanism and to accurately predict the folding rate of several fast-folding proteins and polymers, including a nonbiological helix, polypeptide ␣-helices, a ␤-hairpin, and a three-helix bundle protein from the villin headpiece.
    [Show full text]