Coarse-grained pulling simulations snapshots and contact maps analysis scripts for SARS-CoV-2 RBD variants in complex with H11-H4 nanobody

The dataset consists of 600 snapshots for each replica, totaling 50 pulling trajectories per system, across the four SARS-CoV-2 RBD/H11-H4 complexes.

SARS-CoV-2 RBD WT H11-H4.tar.xz

SARS-CoV-2 RBD Alpha H11-H4.tar.xz

SARS-CoV-2 RBD Delta H11-H4.tar.xz

SARS-CoV-2 RBD XBB.1.5 H11-H4.tar.xz

The CG topology files for each protein complex were created with martinize2 (Kroon et al. 2022) and the Martini 3 force field (Souza et al. 2021). The secondary structure was identified using the DSSP v3.0.0 program (Touw et al. 2015). The GōMartini 3 approach (Poma, Cieplak, and Theodorakis 2017) was employed, substituting conventional harmonic bonds with Lennard-Jones (LJ) interactions based on contact maps obtained from AA-MD simulations. This method utilizes LJ potentials for virtual sites, which allows for the exploration of a broader conformational space, including critical unfolded states, enhancing our understanding of protein dynamics and functionality. A parameter study from standard values of 9.14 kJ/mol to 20 kJ/mol was conducted to determine the optimal value for the LJ potential depth (data not shown). The interaction energy of the contacts in GōMartini was set at 15.0 kJ/mol. This effective value allowed us to recover agreement with all-atom results of previously reported studies (Golcuk et al. 2022). The CG structures were initially minimized in vacuum during 5,000 steps using the steepest descent algorithm. Subsequently, the complexes were solvated in a 10x10x60 nm³ box using the Martini water model. This solvation involved about 38,511 coarse-grained water beads, corresponding to around 154,044 water molecules, with the addition of Na+ and Cl- ions to form a 0.15 M NaCl solution. Systems were then minimized using the same parameters as above. Positional restraints were imposed on the BB beads of each protein to prevent drifting during the equilibration phases. Both the NVT and NPT equilibrations, along with the production phase, utilized the V-rescale thermostat (Bussi, Donadio, and Parrinello 2007). The temperature coupling time constant was set at 1.0 ps for both protein and non-protein parts of the system, keeping the temperature at 300 K. The NVT equilibration was run for 2 ns, with an integration time of 20 fs. During the NPT equilibration and production, an isotropic pressure coupling was used with a compressibility set at 10-4 bar-1 and 1 bar pressure. The NPT equilibration was run for 5 ns, using the C-rescale barostat (Bernetti and Bussi 2020) with a pressure coupling time constant of 18 ps and an integration time of 10 fs. For the production phase, the Parrinello-Rahman barostat (Parrinello and Rahman 1981) was used, with a pressure coupling time constant of 15 ps. The cutoff distances for Coulomb and Van der Waals interactions were set at 1.2 nm across all equilibration and production phases. The pulling simulations (see Figure 1C) were conducted over 1.2 microseconds, with a time step of 20 fs. For the RBD/H11-H4 complexes, specific constraints were applied: the positions of the heavy atoms of the final three residues from the C-terminus of RBD were frozen along the z-axis, and similarly, the coordinates of heavy atoms of residues S126, S127, and K128 in H11-H4 were fixed along the x- and y-axes. The center of mass (COM) of these coordinates was targeted for steered molecular dynamics simulation at a constant speed of 1x10-4 nm/ps and a spring constant of 37.6 kJ/mol nm². A total of 50 independent replicas were conducted for each system using GROMACS 2023 (Abraham et al. 2015).

For further details on the trajectories, please contact Luis F. Cofas-Vargas ([e-mail]).

Gromacs input files

WT input.tar.xz

Alpha input.tar.xz

Delta input.tar.xz

XBB.1.5 input.tar.xz

6ZH9_WT_pull.pdb - Starting structure for contact map calculation (http://pomalab.ippt.pan.pl/GoContactMap/)

6ZH9_WT_pull.map - Contact map

go_system.gro - Protein complex + Gō contacts molecular structure 

go_martini.itp - GōMartini 3 contacts topology

go_molecule1.itp - Topology file for RDB

go_molecule1.itp - Topology file for H11-H4

*.mdp -  Molecular dynamic parameters files required for running simulations

*.ndx - Index file for pulling simulations.

Gromacs output files

WT output.tar.xz

Alpha output.tar.xz

Delta output.tar.xz

XBB.1.5 output.tar.xz

*.gro - Final snapshots for minimization, NVT and NPT equilibrations

*.tpr - Topology and coordinate information for minimization, NVT and NPT equilibrations

pullf*.xvg - Pull-force files for each replica

pullx*.xvg - Pull-coordinate files for each replica

Contact map analysis scripts 

Native_contact_analysis.tar.xz

Native contacts

Distances_native_WT.ipynb - A script that calculates the distances and persistence of native contact pairs along their respective trajectories.

npt_dry.pdb - Topology for calculations

distances.py - Script that is designed to identify and calculate the critical breaking points for each contact pair in a molecular structure. This file must be in the "native_distances" folder

average_native_distance.py - Script for statistics calculation of breaking distance of contact pairs. This file must be in the  "distances" folder.

Graph.ipynb - Script to graph the contact persistence along trajectories. This file must be in the "distances" folder

Nonnative (NON) contacts

NON_contact_analysis.tar.xz

Distances_nonnative_WT.ipynb - A script that calculates the distances and persistence of nonnative contact pairs along their respective trajectories.

sort_contacts.py - Script that sorts native and nonnative contacts.

distances.py -  Script that is designed to identify and calculate the critical breaking points for each contact pair in a molecular structure. This file must be in the "nonnative" folder

count_nonnative_contacts.py -  Script to count the number of nonnative contacts along the trajectories. This file must be in the "nonnative" folder

average_distances.py - Script for statistics calculation of breaking distance of contact pairs. This file must be in the  "distances" folder.

Graph_nonnative.pynb - Script to graph the contact persistence along trajectories. This file must be in the "distances" folder

RMSD analysis

RMSD.tar.xz

RMSD.ipynb - Script to calculate and graph the RMSD of RBD and H11-H4