scholarly journals Static all-atom energetic mappings of the SARS-Cov-2 spike protein and dynamic stability analysis of “Up” versus “Down” protomer states

PLoS ONE ◽  
2020 ◽  
Vol 15 (11) ◽  
pp. e0241168 ◽  
Author(s):  
Michael H. Peters ◽  
Oscar Bastidas ◽  
Daniel S. Kokron ◽  
Christopher E. Henze

The SARS-CoV-2 virion responsible for the current world-wide pandemic COVID-19 has a characteristic Spike protein (S) on its surface that embellishes both a prefusion state and fusion state. The prefusion Spike protein (S) is a large trimeric protein where each protomer may be in a so-called Up state or Down state, depending on the configuration of its receptor binding domain (RBD) within its distal, prefusion S1 domain. The Up state is believed to allow binding of the virion to ACE-2 receptors on human epithelial cells, whereas the Down state is believed to be relatively inactive or reduced in its binding behavior. We have performed detailed all-atom, dominant energy landscape mappings for noncovalent interactions (charge, partial charge, and van der Waals) of the SARS-CoV-2 Spike protein in its static prefusion state based on two recent and independent experimental structure publications. We included both interchain interactions and intrachain (domain) interactions in our mappings in order to determine any telling differences (different so-called “glue” points) between residues in the Up and Down state protomers. The S2 proximal, fusion domain demonstrated no appreciable energetic differences between Up and Down protomers, including interchain as well as each protomer’s intrachain, S1-S2 interactions. However, the S1 domain interactions across neighboring protomers, which include the RBD-NTD cross chain interactions, showed significant energetic differences between Up-Down and Down-Down neighboring protomers. This included, for example, a key RBD residue ARG357 in the Up-Down interaction and a three residue sequence ALA520-PRO521-ALA522, associated with a turn structure in the RBD of the Up state protomer, acting as a stabilizing interaction with the NTD of its neighbor protomer. Additionally, our intra chain dominant energy mappings within each protomer, identified a significant “glue” point or possible “latch” for the Down state protomer between the S1 subdomain, SD1, and the RBD domain of the same protomer that was completely missing in the Up state protomer analysis. Ironically, this dominant energetic interaction in the Down state protomer involved the backbone atoms of the same three residue sequence ALA520-PRO521-ALA522 of the RBD with the amino acid R-group of GLN564 in the SD1 domain. Thus, this same three residue sequence acts as a stabilizer of the RBD in the Up conformation through its interactions with its neighboring NTD chain and a kind of latch in the Down state conformation through its interactions with its own SD1 domain. The dominant interaction energy residues identified here are also conserved across reported variations of SARS-CoV-2, as well as the closely related virions SARS-Cov and the bat corona virus RatG13. We conducted preliminary molecular dynamics simulations across 0.1 μ seconds to see if this latch provided structural stability and indeed found that a single point mutation (Q564G) resulted in the latch releasing transforming the protomer from the Down to the Up state conformation. Full trimeric Spike protein studies of the same mutation across all protomers, however, did not exhibit latch release demonstrating the critical importance of interchain interactions across the S1 domain, including RBD-NTD neighboring chain interactions. Therapies aimed at disrupting these noncovalent interactions could be a viable route for the physico-chemical mitigation of this deadly virion.

2020 ◽  
Author(s):  
Michael H. Peters ◽  
Oscar Bastidas ◽  
Daniel S. Kokron ◽  
Christopher Henze

AbstractThe SARS-Cov-2 virion responsible for the current world-wide pandemic Covid-19 has a characteristic Spike protein (S) on its surface that embellishes both a prefusion state and fusion state. The prefusion Spike protein (S) is a large trimeric protein where each protomer may be in a so-called Up state or Down state, depending on the configuration of its receptor binding domain (RBD). The Up state is believed to allow binding of the virion to ACE-2 receptors on human epithelial cells, whereas the Down state is believed to be relatively inactive or reduced in its binding behavior. We have performed detailed all-atom, dominant energy landscape mappings for noncovalent interactions (charge, partial charge, and van der Waals) of the SARS-Cov-2 Spike protein in its static prefusion state based on recent structural information. We included both interchain interactions and intrachain (domain) interactions in our mappings in order to determine any telling differences (different so-called “glue” points) between residues in the Up and Down state protomers. In general, the S2 or fusion machinery domain of S is relatively rigid with strong noncovalent interactions facilitated by helical secondary structures, whereas the S1 domain, which contains the RBD and N-terminal domain (NTD), is relatively more flexible and characterized by beta strand structural motifs. The S2 domain demonstrated no appreciable energetic differences between Up and Down protomers, including interchain as well as each protomer’s intrachain, S1-S2 interactions. However, the S1 domain interactions across neighboring protomers, which include the RBD-NTD cross chain interactions, showed significant energetic differences between Up-Down and Down-Down neighboring protomers. Surprisingly, the Up-Down, RBD-NTD interactions were overall stronger and more numerous than the Down-Down cross chain interactions, including the appearance of the three residue sequence ALA520-PRO521-ALA522 associated with a turn structure in the RBD of the Up state protomer. Additionally, our intrachain dominant energy mappings within each protomer, identified a significant “glue” point or possible “latch” for the Down state protomer between the S1 subdomain, SD1, and the RBD domain of the same protomer that was completely missing in the Up state protomer analysis. Ironically, this dominant energetic interaction in the Down state protomer involved the backbone atoms of the same three residue sequence ALA520-PRO521-ALA522 of the RBD with the R-group of GLN564 in the SD1 domain. Thus, this same three residue sequence acts as a stabilizer of the RBD in the Up conformation through its interactions with its neighboring NTD chain and a kind of latch in the Down state conformation through its interactions with its own SD1 domain. The dominant interaction energy residues identified here are also conserved across reported variations of SARS-Cov-2, as well as the closely related virions SARS-Cov and the bat corona virus RatG13. To help verify the potential latch for the Down state protomer, we conducted some preliminary molecular dynamic simulations that effectively turn off this specific latch glue point via a single point mutation of GLN564. Interestingly, the single point mutation lead to the latch releasing in less than a few nanoseconds, but the latch remained fixed in the wild state protomer for up to 0.1 microseconds that were simulated. Many more detailed studies are needed to understand the dynamics of the Up and Down states of the Spike protein, including the stabilizing chain-chain interactions and the mechanisms of transition from Down to Up state protomers. Nonetheless, static dominant energy landscape mappings and preliminary molecular dynamic studies given here may represent a useful starting point for more detailed dynamic analyses and hopefully an improved understanding of the structure-function relationship of this highly complex protein associated with COVID-19.


mBio ◽  
2016 ◽  
Vol 7 (5) ◽  
Author(s):  
Ana R. Pereira ◽  
Jen Hsin ◽  
Ewa Król ◽  
Andreia C. Tavares ◽  
Pierre Flores ◽  
...  

ABSTRACT A mechanistic understanding of the determination and maintenance of the simplest bacterial cell shape, a sphere, remains elusive compared with that of more complex shapes. Cocci seem to lack a dedicated elongation machinery, and a spherical shape has been considered an evolutionary dead-end morphology, as a transition from a spherical to a rod-like shape has never been observed in bacteria. Here we show that a Staphylococcus aureus mutant (M5) expressing the ftsZ G193D allele exhibits elongated cells. Molecular dynamics simulations and in vitro studies indicate that FtsZ G193D filaments are more twisted and shorter than wild-type filaments. In vivo , M5 cell wall deposition is initiated asymmetrically, only on one side of the cell, and progresses into a helical pattern rather than into a constricting ring as in wild-type cells. This helical pattern of wall insertion leads to elongation, as in rod-shaped cells. Thus, structural flexibility of FtsZ filaments can result in an FtsZ-dependent mechanism for generating elongated cells from cocci. IMPORTANCE The mechanisms by which bacteria generate and maintain even the simplest cell shape remain an elusive but fundamental question in microbiology. In the absence of examples of coccus-to-rod transitions, the spherical shape has been suggested to be an evolutionary dead end in morphogenesis. We describe the first observation of the generation of elongated cells from truly spherical cocci, occurring in a Staphylococcus aureus mutant containing a single point mutation in its genome, in the gene encoding the bacterial tubulin homologue FtsZ. We demonstrate that FtsZ-dependent cell elongation is possible, even in the absence of dedicated elongation machinery.


2018 ◽  
Vol 20 (40) ◽  
pp. 25901-25909 ◽  
Author(s):  
Rodrigo Ochoa ◽  
Miguel A. Soler ◽  
Alessandro Laio ◽  
Pilar Cossio

Single-point mutation protocols based on backbone-dependent rotamer libraries show the best performance in predicting equilibrium configurations from molecular dynamics simulations.


1998 ◽  
Vol 140 (1) ◽  
pp. 91-99 ◽  
Author(s):  
Malini Vashishtha ◽  
Thomas Phalen ◽  
Marianne T. Marquardt ◽  
Jae S. Ryu ◽  
Alice C. Ng ◽  
...  

Membrane fusion and budding are key steps in the life cycle of all enveloped viruses. Semliki Forest virus (SFV) is an enveloped alphavirus that requires cellular membrane cholesterol for both membrane fusion and efficient exit of progeny virus from infected cells. We selected an SFV mutant, srf-3, that was strikingly independent of cholesterol for growth. This phenotype was conferred by a single amino acid change in the E1 spike protein subunit, proline 226 to serine, that increased the cholesterol independence of both srf-3 fusion and exit. The srf-3 mutant emphasizes the relationship between the role of cholesterol in membrane fusion and virus exit, and most significantly, identifies a novel spike protein region involved in the virus cholesterol requirement.


Pathogens ◽  
2021 ◽  
Vol 10 (7) ◽  
pp. 869
Author(s):  
Chee-Hing Yang ◽  
Hui-Chun Li ◽  
Wen-Han Lee ◽  
Shih-Yen Lo

The COVID-19 pandemic was caused by SARS-CoV-2 infection. To prevent the spread of SARS-CoV-2, an effective vaccine is required. Two linear peptides from potential B-cell epitopes in the spike protein of SARS-CoV-2 (a.a. 440–460; a.a. 494–506) were synthesized and used to immunize rabbits. High-titer antibodies of IgG were produced, purified, and verified by Western blot analysis. Antibodies against these two epitopes could effectively neutralize SARS-CoV-2 pseudoviral particles with the spike proteins from not only the original strain (basal; wild-type), but also a strain with a single point mutation (D614G), and two other emerging variants (the Alpha and Beta variants) prevalent around the world, but not from SARS-CoV. In conclusion, antibodies against these two epitopes are protective. This information is important for the development of vaccines against SARS-CoV-2.


Author(s):  
G. Zhang ◽  
S. Pomplun ◽  
A. R. Loftis ◽  
X. Tan ◽  
A. Loas ◽  
...  

AbstractCoronavirus disease 19 (COVID-19) is an emerging global health crisis. With over 7 million confirmed cases to date, this pandemic continues to expand, spurring research to discover vaccines and therapies. SARS-CoV-2 is the novel coronavirus responsible for this disease. It initiates entry into human cells by binding to angiotensin-converting enzyme 2 (ACE2) via the receptor binding domain (RBD) of its spike protein (S). Disrupting the SARS-CoV-2-RBD binding to ACE2 with designer drugs has the potential to inhibit the virus from entering human cells, presenting a new modality for therapeutic intervention. Peptide-based binders are an attractive solution to inhibit the RBD-ACE2 interaction by adequately covering the extended protein contact interface. Using molecular dynamics simulations based on the recently solved cryo-EM structure of ACE2 in complex with SARS-CoV-2-RBD, we observed that the ACE2 peptidase domain (PD) α1 helix is important for binding SARS-CoV-2-RBD. Using automated fast-flow peptide synthesis, we chemically synthesized a 23-mer peptide fragment of the ACE2 PD α1 helix (SBP1) composed entirely of proteinogenic amino acids. Chemical synthesis of SBP1 was complete in 1.5 hours, and after work up and isolation >20 milligrams of pure material was obtained. Bio-layer interferometry (BLI) revealed that SBP1 associates with micromolar affinity to insect-derived SARS-CoV-2-RBD protein obtained from Sino Biological. Association of SBP1 was not observed to an appreciable extent to HEK cell-expressed SARS-CoV-2-RBD proteins and insect-derived variants acquired from other vendors. Moreover, competitive BLI assays showed SBP1 does not outcompete ACE2 binding to Sino Biological insect-derived SARS-CoV-2-RBD. Further investigations are ongoing to gain insight into the molecular and structural determinants of the variable binding behavior to different SARS-CoV-2-RBD protein variants.


1998 ◽  
Vol 18 (01) ◽  
pp. 1-10
Author(s):  
A. Hillarp ◽  
S. Rosen ◽  
B. Zöller ◽  
B. Dahlbäck

SummaryVitamin K-dependent protein C is an important regulator of blood coagulation. After its activation on the endothelial cell surface by thrombin bound to thrombomodulin, it cleaves and inactivates procoagulant cofactors Va and Villa, protein S and intact factor V working as cofactors. Until recently, genetic defects of protein C or protein S were, together with antithrombin III deficiency, the established major causes of familial venous thromboembolism, but they were found in fewer than 5-10% of patients with thrombosis. In 1993, inherited resistance to activated protein C (APC) was described as a major risk factor for venous thrombosis. It is found in up to 60% of patients with venous thrombosis. In more than 90% of cases, the molecular background for the APC resistance is a single point mutation in the factor V gene, which predicts substitution of an arginine (R) at position 506 by a glutamine (Q). Mutated factor V (FV: Q506) is activated by thrombin or factor Xa in normal way, but impaired inactivation of mutated factor Va by APC results in life-long hypercoagulability. The prevalence of the FV:Q506 allele in the general population of Western countries varies between 2 and 15%, whereas it is not found in several other populations with different ethnic backgrounds. Owing to the high prevalence of FV:Q506 in Western populations, it occasionally occurs in patients with deficiency of protein S, protein C, or antithrombin III. Individuals with combined defects suffer more severely from thrombosis, and often at a younger age, than those with single defects, suggesting severe thrombophilia to be a multigenetic disease.


2019 ◽  
Vol 92 (2) ◽  
pp. 241-247
Author(s):  
Padmaja D. Wakchaure ◽  
Bishwajit Ganguly

Riboswitches are the type of regulatory elements present in the untranslated region of mRNA and specifically bind to the natural ligand to regulate gene expression. This binding specificity can be affected by even single point mutation incorporated in the core of the riboswitch. In this work, we have examined the mutations at the binding site residue in Flavin Mononucleotide (FMN) riboswitch structure with 30ns molecular dynamics simulations. The interaction of ligand (FMN) with riboswitch has been characterized using root mean square deviation, hydrogen bonding analysis, and the calculated binding affinities. Mutation at A48G and G62U show the enhanced binding energy however, the mutation at A85G, are energetically unfavorable compared to the wild type. This work gives valuable insight into the structures and energetics of the mutated FMN riboswitch to design new hits for biological applications.


Sign in / Sign up

Export Citation Format

Share Document