Utilizing Computational Techniques to Accelerate Discovery in Peanut Allergenicity: A Case Study

Computational tools and techniques are attractive due to the low-cost and scope of application and are becoming more prevalent and accepted across a variety of fields. With the advancement of high performance computing resources and the recent utility of GPU architecture, molecular dynamics (MD) simulations have become a nearly ubiquitous tool to interrogate protein structure-function relationships. Understanding these structure-function relationships is fundamental to protein biochemistry and facilitates avenues to modulate activity in disease states. Herein, we studied and compared the ability of several GPU- and CPU-based clusters to accelerate mechanistic and therapeutic discovery with MD simulations. MD simulations on the peanut protein Ara h 2, the causative agent of peanut allergenicity, from 10 different species of the genus Arachis were conducted in triplicate for 1 μs to probe structural dynamics and epitope solvent accessibility relevant to allergenicity. These simulations provide molecular insight into the causes and potential potency of the allergenic response. Ultimately, these computational findings and insights can guide experimental design, increasing the impact of therapeutics and understanding of disease states. Further, integrating GPUs into these parallel computing operations accelerates computations and provides researchers with tools to perform MD simulations without large CPU centers. We have thus provided an increased understanding of protein structure-function relationships through the usage of high performance computing and have shown GPU scaling power with MD simulations.

[1]  M. Sippl Recognition of errors in three‐dimensional structures of proteins , 1993, Proteins.

[2]  Torsten Schwede,et al.  BIOINFORMATICS Bioinformatics Advance Access published November 12, 2005 The SWISS-MODEL Workspace: A web-based environment for protein structure homology modelling , 2022 .

[3]  R. Helm,et al.  Identification and mutational analysis of the immunodominant IgE binding epitopes of the major peanut allergen Ara h 2. , 1997, Archives of biochemistry and biophysics.

[4]  K. Lindorff-Larsen,et al.  Picosecond to Millisecond Structural Dynamics in Human Ubiquitin. , 2016, The journal of physical chemistry. B.

[5]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[6]  Ian W. Davis,et al.  Structure validation by Cα geometry: ϕ,ψ and Cβ deviation , 2003, Proteins.

[7]  F. Melo,et al.  Assessing protein structures with a non-local atomic interaction energy. , 1998, Journal of molecular biology.

[8]  X. Daura,et al.  Peptide Folding: When Simulation Meets Experiment , 1999 .

[9]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[10]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[11]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[12]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[13]  S. Vieths,et al.  Structure and stability of 2S albumin-type peanut allergens: implications for the severity of peanut allergic reactions. , 2006, The Biochemical journal.

[14]  B. L. de Groot,et al.  CHARMM36m: an improved force field for folded and intrinsically disordered proteins , 2016, Nature Methods.

[15]  Marco Biasini,et al.  Toward the estimation of the absolute quality of individual protein structure models , 2010, Bioinform..

[16]  J. Wal,et al.  Capacity of purified peanut allergens to induce degranulation in a functional in vitro assay: Ara h 2 and Ara h 6 are the most efficient elicitors , 2009, Clinical and experimental allergy : journal of the British Society for Allergy and Clinical Immunology.

[17]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[18]  Joseph L. Baumert,et al.  Conformational stability of digestion-resistant peptides of peanut conglutins reveals the molecular basis of their allergenicity , 2016, Scientific Reports.

[19]  J. Hourihane Peanut allergy. , 2011, Pediatric clinics of North America.

[20]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[21]  Diwakar Shukla,et al.  OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular Simulation. , 2013, Journal of chemical theory and computation.

[22]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.

[23]  Martin Fechner,et al.  Best bang for your buck: GPU nodes for GROMACS biomolecular simulations , 2015, J. Comput. Chem..

[24]  G. N. Ramachandran,et al.  Stereochemistry of polypeptide chain configurations. , 1963, Journal of molecular biology.

[25]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[26]  M. Kulis,et al.  The 2S albumin allergens of Arachis hypogaea, Ara h 2 and Ara h 6, are the major elicitors of anaphylaxis and can effectively desensitize peanut‐allergic mice , 2012, Clinical and experimental allergy : journal of the British Society for Allergy and Clinical Immunology.

[27]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[28]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[29]  Manfred J. Sippl,et al.  Thirty years of environmental health research--and growing. , 1996, Nucleic Acids Res..

[30]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..