Title: All atom simulations snapshots and contact maps analysis scripts for SARS-CoV-2002 and SARS-CoV-2 spike proteins with and without ACE2 enzyme

Type Dataset Rodrigo A. Moreira, Mateusz Chwastyk, Joseph L. Baker, Horacio V Guzman, Adolfo B. Poma (2020): All atom simulations snapshots and contact maps analysis scripts for SARS-CoV-2002 and SARS-CoV-2 spike proteins with and without ACE2 enzyme. Zenodo. Dataset. https://zenodo.org/record/3817447

Authors: Rodrigo A. Moreira (Institute of Fundamental Technological Research, Polish Academy of Sciences, Pawinskiego 5B, 02-106 Warsaw, Poland) ; Mateusz Chwastyk (Institute of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warasw, Poland) ; Joseph L. Baker (Department of Chemistry, The College of New Jersey, 2000 Pennington Road, Ewing, NJ 08628, United States) ; Horacio V Guzman (Jožef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia) ; Adolfo B. Poma (Institute of Fundamental Technological Research, Polish Academy of Sciences, Pawinskiego 5B, 02-106 Warsaw, Poland) ;

Links

• Item record in Zenodo
• Digital object URL

Summary

The dataset contains a total of 40 snapshots of the four trajectories (10 snapshots each system = two per replica x 5 replicas/system):

Molecular dynamics simulation trajectories (320ns each) have been performed using the Amber ff14SB force field running with the Amber18 package at the the NSF-funded (OAC-1826915, OAC-1828163) ELSA high performance computing cluster at The College of New Jersey. Under the following simulation methodology:

All-atom simulations were carried out with Amber18 (ambermd.org), and system components (protein, ions, water) were modeled with the included FF14SB and TIP3P parameter sets. Energy minimization used CPU pmemd, while later simulation stages used GPU pmemd. CoV2 and CoV1 systems with one RBD up (with/without ACE2) were solvated in 12 angstrom water shells. Cysteine residues identified in the initial models as having a disulfide bond (DB) were bonded using tLeap. All simulations used 0.150 M NaCl. Hydrogen mass repartitioning was applied only to the protein to enable a 4 fs timestep (https://pubs.acs.org/doi/abs/10.1021/ct5010406). The SHAKE algorithm was applied to hydrogens, and a real-space cutoff of 8 angstroms was used. Periodic boundary conditions were applied and PME was used for long-range electrostatics. Minimization was by steepest descent (2000 steps) followed by conjugate gradient (3000 steps). Heating used two stages: (1) NVT heating from 0 K to 100 K (50 ps), and (2) NPT heating from 100 K to 300 K (100 ps). Restraints of 10 kcal mol-1 angstrom-2 were applied during minimization and heating to C-alpha atoms. During 6 ns of equilibration at 300 K C-alpha restraints were gradually reduced from 10 kcal mol-1 angstrom-2 to 0.1 kcal mol-1 angstrom-2. Finally, restraints were released and 320 ns unrestrained production simulations were carried out for CoV2 and CoV1 systems. Production simulations began from the final equilibrated snapshots, and five copies of each system were simulated. As unrestrained systems can freely rotate we monitored simulations for any close contacts and found that in one copy of the CoV1 simulation without ACE2 and one RBD up that a few contacts close to 8 angstrom occur near the end of the 320 ns between the RBD and a different subdomain of the spike complex in a periodic image. However this did not influence analyzed structural properties which is verified by comparing results across simulations. The Monte Carlo barostat was used to maintain pressure (1 atm), and the Langevin thermostat was used to maintain 300 K temperature (collision frequency 1 ps-1), as implemented in Amber18. In aggregate, nearly 7 microseconds of simulation of systems ranging from 396,147 to 879,100 atoms was carried out for this work. For further details on the trajectories, please contact Joseph Baker (bakerj@tcnj.edu).

Regarding the contact map analysis scripts (contactMaps_Analysis.tar.gz), they contain the following workflow:

contactmap      --> source files from contact_map executable process_nc.sh   --> convert raw data from all-atom simulation to numbered PDB files and get the contact maps frequency.lua   --> read a set of PDB files and output the frequency count for each contact consensus.fasta --> align sequence of Covid19 and SARS from Chimera consensus.lua   --> read data previously generated and compute the frequency per residue, among other things. consensus.sh    --> input information to consensus.lua consensus.gp    --> gnuplot script to plot figures

This dataset and the code is part of tripartite collaboration between:

More information

• DOI: 10.5281/zenodo.3817447
• Language: en

Subjects

• COVID-19, SARS, Mechanical Stability, RBD, Molecular Dynamics, Coarse-Grain, coronavirus

Dates

• Publication date: 2020
• Issued: May 08, 2020

Rights

• https://creativecommons.org/licenses/by/4.0/legalcode Creative Commons Attribution 4.0 International
• info:eu-repo/semantics/openAccess Open Access

Format

electronic resource

Relateditems