Coronavirus Exercises Part 1

Comparing structures with Qres

Recall the following two toy protein structures from the part 1 exercises.

image-center Two hypothetical protein structures with vectorizations into eight points each.

The formula for Qres of the i-th sampled point in a protein with N amino acids is reproduced below:

\[Q_{res}^{(i)} = \dfrac{1}{N-k} \sum^{residues}_{j\neq i-1,i,i+1} \textrm{exp}[-\dfrac{[d(s_i,s_j)-d(t_i,t_j)]^2}{2\sigma^2_{i,j}}]\, .\]

Exercise: Assuming that the first sampled point in the two figures are, respectively, (1,2) and (1,3), compute the Qres of the fourth sampled point.

In a tutorial, we computed Qres between SARS and SARS-CoV-2 RBD to find local differences between the two structures. We will repeat this analysis using the hemoglobin subunit alpha protein. Using the PDB files of human hemoglobin (1SI4) and mako shark hemoglobin (3MKB), upload the protein structures into VMD. Align the structures using “Multiseq” and visualize the structures using Qres as the coloring parameter.

Exercise: Use the visualizations to determine where the protein structures vary the most. Does what you find make sense intuitively?

Calculating interaction energy

In this exercise, we will be performing a short analysis of the SARS-CoV-2 Spike protein interaction energy with the human receptor ACE2 using NAMD. Go to the protein databank entry for the complex (6VW1) and download 6VW1.pdb. Upload the structure into VMD and navigate to Extension > Analysis > NAMD Energy. Calculate the energy types “VDW” and “Elec” using the following selections:

  Selection 1 Selection 2
A protein and chain B protein and chain F and resid 482 to 486
B protein and chain B protein and chain F and resid 383 to 394

Exercise: How do the results of part A compare with that of part B? Can you justify the results? (Hint: Visualize the specific residues).

Note: For help, visit the module tutorials for VMD visualization and NAMD Energy. For additional troubleshooting, download NamdEnergyAssignment.zip and follow README.txt.

Visualizing glycans

In the conclusion of part 2 of this module, we mentioned that the surface of the surface of viruses and host cells are “fuzzy” because they are covered by structures called glycans. As you may expect, SARS-CoV and SARS-CoV-2 are also covered by this “glycan shield”, in which glycosylation of surface antigens allows the virus to hide from antibody detection.

We will use what we learned about VMD in a tutorial to visualize glycans on the spike protein of both viruses. First, we will need to download the protein structures of the SARS-CoV spike protein (PDB entry: 5xlr and the SARS-CoV-2 spike protein (PDB entry: 6vxx. Load each structure into VMD and navigate to Graphics > Representations. For VMD, there is no specific keyword to select glycans, and so we will use a workaround with the keywords: "not protein and not water".

image-center

Exercise: Create the representation. Assuming that all the non-proteins are glycans, which structure contains the most number of glycans? Do you think this supports SARS-CoV-2’s higher infectivity compared to SARS-CoV?

Contact maps

Exercise: In our GNM tutorial, we created the contact map using the threshold of 20Å. Try making the contact map of one of the chains of SARS-CoV-2 spike protein 6vxx with thresholds of 10Å and 40Å. How different are the resulting contact maps?

Comments