Evolution of the SARS-CoV-2 proteome in three dimensions (3D) during the first 6 months of the COVID-19 pandemic.

Understanding the molecular evolution of the SARS-CoV-2 virus as it continues to spread in communities around the globe is important for mitigation and future pandemic preparedness. Three-dimensional structures of SARS-CoV-2 proteins and those of other coronavirusess archived in the Protein Data Bank were used to analyze viral proteome evolution during the first 6 months of the COVID-19 pandemic. Analyses of spatial locations, chemical properties, and structural and energetic impacts of the observed amino acid changes in >48 000 viral isolates revealed how each one of 29 viral proteins have undergone amino acid changes. Catalytic residues in active sites and binding residues in protein-protein interfaces showed modest, but significant, numbers of substitutions, highlighting the mutational robustness of the viral proteome. Energetics calculations showed that the impact of substitutions on the thermodynamic stability of the proteome follows a universal bi-Gaussian distribution. Detailed results are presented for potential drug discovery targets and the four structural proteins that comprise the virion, highlighting substitutions with the potential to impact protein structure, enzyme activity, and protein-protein and protein-nucleic acid interfaces. Characterizing the evolution of the virus in three dimensions provides testable insights into viral protein function and should aid in structure-based drug discovery efforts as well as the prospective identification of amino acid substitutions with potential for drug resistance.

Undergraduate structural biology education: A shift from users to developers of computation and simulation tools.

The use of theory and simulation in undergraduate education in biochemistry, molecular biology, and structural biology is now common, but the skills students need and the curriculum instructors have to train their students are evolving. The global pandemic and the immediate switch to remote instruction forced instructors to reconsider how they can use computation to teach concepts previously approached with other instructional methods. In this review, we survey some of the curricula, materials, and resources for instructors who want to include theory, simulation, and computation in the undergraduate curriculum. There has been a notable progression from teaching students to use discipline-specific computational tools to developing interactive computational tools that promote active learning to having students write code themselves, such that they view computation as another tool for solving problems. We are moving toward a future where computational skills, including programming, data analysis, visualization, and simulation, will no longer be considered an optional bonus for students but a required skill for the 21st century STEM (Science, Technology, Engineering, and Mathematics) workforce; therefore, all physical and life science students should learn to program in the undergraduate curriculum.

Evolution of the SARS-CoV-2 proteome in three dimensions (3D) during the first six months of the COVID-19 pandemic.

Three-dimensional structures of SARS-CoV-2 and other coronaviral proteins archived in the Protein Data Bank were used to analyze viral proteome evolution during the first six months of the COVID-19 pandemic. Analyses of spatial locations, chemical properties, and structural and energetic impacts of the observed amino acid changes in >48,000 viral proteome sequences showed how each one of the 29 viral study proteins have undergone amino acid changes. Structural models computed for every unique sequence variant revealed that most substitutions map to protein surfaces and boundary layers with a minority affecting hydrophobic cores. Conservative changes were observed more frequently in cores versus boundary layers/surfaces. Active sites and protein-protein interfaces showed modest numbers of substitutions. Energetics calculations showed that the impact of substitutions on the thermodynamic stability of the proteome follows a universal bi-Gaussian distribution. Detailed results are presented for six drug discovery targets and four structural proteins comprising the virion, highlighting substitutions with the potential to impact protein structure, enzyme activity, and functional interfaces. Characterizing the evolution of the virus in three dimensions provides testable insights into viral protein function and should aid in structure-based drug discovery efforts as well as the prospective identification of amino acid substitutions with potential for drug resistance.

Virus-like glycodendrinanoparticles displaying quasi-equivalent nested polyvalency upon glycoprotein platforms potently block viral infection.

Ligand polyvalency is a powerful modulator of protein-receptor interactions. Host-pathogen infection interactions are often mediated by glycan ligand-protein interactions, yet its interrogation with very high copy number ligands has been limited to heterogenous systems. Here we report that through the use of nested layers of multivalency we are able to assemble the most highly valent glycodendrimeric constructs yet seen (bearing up to 1,620 glycans). These constructs are pure and well-defined single entities that at diameters of up to 32 nm are capable of mimicking pathogens both in size and in their highly glycosylated surfaces. Through this mimicry these glyco-dendri-protein-nano-particles are capable of blocking (at picomolar concentrations) a model of the infection of T-lymphocytes and human dendritic cells by Ebola virus. The high associated polyvalency effects (ß>10(6), ß/N ~10(2)-10(3)) displayed on an unprecedented surface area by precise clusters suggest a general strategy for modulation of such interactions.

Analysis of the dispersity in carbohydrate loading of synthetic glycoproteins using MALDI-TOF mass spectrometry.

Statistical correlation of mass spectrum peak broadening with product dispersity in protein conjugation reactions allows more detailed characterization of putative therapeutic conjugates.

Fluoroglycoproteins: ready chemical site-selective incorporation of fluorosugars into proteins.

A tag-and-modify strategy allows the practical synthesis of homogenous fluorinated glyco-amino acids, peptides and proteins carrying a fluorine label in the sugar and allows access to first examples of directly radiolabelled ([(18)F]-glyco)proteins.

Mutations in the fourth EGF-like domain affect thrombomodulin-induced changes in the active site of thrombin.

A number of alanine and more conservative mutants of residues in the fourth domain of thrombomodulin (TM) were prepared and assayed for protein C activation and for thrombin binding. Several of the alanine mutations appeared to cause misfolding or structural defects as assessed by poor expression and/or NMR HSQC experiments, while more conservative mutations at the same site appeared to allow correct folding and preserved activity. Several of the conservative mutants bound more weakly to thrombin despite the fact that the fourth domain does not directly contact thrombin in the crystal structure of the thrombin-TM complex. A few of the mutant TM fragments bound thrombin with an affinity similar to that of the wild type but exhibited decreases in k cat for protein C activation. These mutants were also less able to cause a change in the steady state fluorescence of fluorescein-EGR-chloromethylketone bound to the active site of thrombin. These results suggest that some residues within the fourth domain of TM may primarily interact with protein C but others are functionally important for altering the way TM interacts with thrombin. Residues in the fourth domain that primarily affect k cat for protein C activation may do this by changing the active site of thrombin.

Solvent accessibility of protein surfaces by amide H/2H exchange MALDI-TOF mass spectrometry.

One advantage of detecting amide H/2H exchange by mass spectrometry instead of NMR is that the more rapidly exchanging surface amides are still detectable. In this study, we present quench-flow amide H/2H exchange experiments to probe how rapidly the surfaces of two different proteins exchange. We compared the amide H/2H exchange behavior of thrombin, a globular protein, and IkappaBalpha, a nonglobular protein, to explore any differences in the determinants of amide H/2H exchange rates for each class of protein. The rates of exchange of only a few of the surface amides were as rapid as the "intrinsic" exchange rates measured for amides in unstructured peptides. Most of the surface amides exchanged at a slower rate, despite the fact that they were not seen to be hydrogen bonded to another protein group in the crystal structure. To elucidate the influence of the surface environment on amide H/2H exchange, we compared exchange data with the number of amides participating in hydrogen bonds with other protein groups and with the solvent accessible surface area. The best correlation with amide H/2H exchange was found with the total solvent accessible surface area, including side chains. In the case of the globular protein, the correlation was modest, whereas it was well correlated for the nonglobular protein. The nonglobular protein also showed a correlation between amide exchange and hydrogen bonding. These data suggest that other factors, such as complex dynamic behavior and surface burial, may alter the expected exchange rates in globular proteins more than in nonglobular proteins where all of the residues are near the surface.

Amide H/2H exchange reveals a mechanism of thrombin activation.

Thrombin is a dual action serine protease in the blood clotting cascade. Similar to other clotting factors, thrombin is mainly present in the blood in a zymogen form, prothrombin. Although the two cleavage events required to activate thrombin are well-known, little is known about why the thrombin precursors are inactive proteases. Although prothrombin is much larger than thrombin, prethrombin-2, which contains all of the same amino acids as thrombin, but has not yet been cleaved between Arg320 and Ile321, remains inactive. Crystal structures of both prethrombin-2 and thrombin are available and show almost no differences in the active site conformations. Slight differences were, however, seen in the loops surrounding the active site, which are larger in thrombin than in most other trypsin-like proteases, and have been shown to be important for substrate specificity. To explore whether the dynamics of the active site loops were different in the various zymogen forms of thrombin, we employed amide H/(2)H exchange experiments to compare the exchange rates of regions of thrombin with the same regions of prothrombin, prethrombin-2, and meizothrombin. Many of the surface loops showed less exchange in the zymogen forms, including the large loop corresponding to anion binding exosite 1. Conversely, the autolysis loop and sodium-binding site exchanged more readily in the zymogen forms. Prothrombin and prethrombin-2 gave nearly identical results while meizothrombin in some regions more closely resembled active thrombin. Thus, cleavage of the Arg320-Ile321 peptide bond is the key to formation of the active enzyme, which involves increased dynamics of the substrate-binding loops and decreased dynamics of the catalytic site.

Ligand-induced conformational changes in the acetylcholine-binding protein analyzed by hydrogen-deuterium exchange mass spectrometry.

Recent x-ray crystallographic studies of the acetylcholine-binding protein (AChBP) suggest that loop C, found at the circumference of the pentameric molecule, shows distinctive conformational changes upon antagonist and agonist occupation. We have employed hydrogen-deuterium exchange mass spectrometry to examine the influence of bound ligands on solvent exposure of AChBP. Quantitative measurements of deuterium incorporation are possible for approximately 56% of the Lymnaea AChBP sequence, covering primarily the outer surface of AChBP. In the apoprotein, two regions flanking the ligand occupation site at the subunit interface, loop C (residues 175-193) and loop F (residues 164-171), show greater extents of solvent exchange than other regions of the protein including the N- and C-terminal regions. Occupation by nicotinic agonists, epibatidine and lobeline, and nicotinic antagonists, methyllycaconitine, alpha-bungarotoxin, and alpha-cobratoxin, markedly restricts the exchange of loop C amide protons, influencing both the rates and degrees of exchange. Solvent exposure of loop C and its protection by ligand suggest that in the apoprotein, loop C exhibits rapid fluctuations in an open conformation. Bound agonists restrict solvent exposure through loop closure, whereas the larger antagonists restrict solvent exposure largely through occlusion of solvent. Loop F, found on the complementary subunit surface at the interface, also reveals ligand selective changes in amide proton exchange rates. Agonists do not affect solvent accessibility of loop F, whereas certain antagonists cause subtle accessibility changes. These results reveal dynamic states and fluctuating movements in the vicinity of the binding site for unligated AChBP that can be influenced selectively by ligands.

Thrombomodulin tightens the thrombin active site loops to promote protein C activation.

Thrombomodulin (TM) forms a 1:1 complex with thrombin. Whereas thrombin alone cleaves fibrinogen to make the fibrin clot, the thrombin-TM complex cleaves protein C to initiate the anticoagulant pathway. Crystallographic investigations of the complex between thrombin and TMEGF456 did not show any changes in the thrombin active site. Therefore, research has focused recently on how TM may provide a docking site for the protein C substrate. Previous work, however, showed that when the thrombin active site was occupied with substrate analogues labeled with fluorophores, the fluorophores responded differently to active (TMEGF1-6) versus inactive (TMEGF56) fragments of TM. To investigate this further, we have carried out amide H/(2)H exchange experiments on thrombin in the presence of active (TMEGF45) and inactive (TMEGF56) fragments of TM. Both on-exchange and off-exchange experiments show changes in the thrombin active site loops, some of which are observed only when the active TM fragment is bound. These results are consistent with the previously observed fluorescence changes and point to a mechanism by which TM changes the thrombin substrate specificity in favor of protein C cleavage.

Allosteric changes in solvent accessibility observed in thrombin upon active site occupation.

The solvent accessibility of thrombin in its substrate-free and substrate-bound forms has been compared by amide hydrogen/deuterium (H/(2)H) exchange. The optimized inhibitor peptide dPhe-Pro-Arg chloromethyl ketone (PPACK) was used to simulate the substrate-bound form of thrombin. These studies were motivated by the lack of observed changes in the active site of thrombin in the crystal structure of the thrombin-thrombomodulin complex. This result appeared to contradict amide exchange studies on the thrombin-thrombomodulin complex that suggested subtle changes occur in the active site loops upon thrombomodulin binding. Our results show that two active site loops, residues 214-222 and residues 126-132, undergo decreases in solvent accessibility due to steric contacts with PPACK substrate. However, we also observe two regions outside the active site undergoing solvent protection upon substrate binding. The first region corresponds to anion binding exosite 1, and the second is a beta-strand-containing loop which runs through the core of the molecule and contains Trp141 which makes critical contacts with anion binding exosite 1. These results indicate two pathways of allosteric change that connect the active site to the distal anion binding exosite 1.