Molecular dynamic study of SARS-CoV-2 with various S protein mutations and their effect on thermodynamic properties.
Ontology highlight
ABSTRACT:
SUBMITTER: Abdalla M
PROVIDER: S-EPMC8576119 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Structural proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are potential drug targets due to their role in the virus life cycle. The envelope (E) protein is one of the structural proteins; plays a critical role in virulency. However, the emergence of mutations oftenly leads to drug resistance and may also play a vital role in virus stabilization and evolution. In this study, we aimed to identify mutations in E proteins that affect the protein stability. About 0.3 million complete whole genome sequences were analyzed to screen mutations in E protein. All these mutations were subjected to stability prediction using the DynaMut server. The most common mutations that were detected at the C-terminal domain, Ser68Phe, Pro71Ser, and Leu73Phe, were examined through molecular dynamics (MD) simulations for a 100ns period. The sequence analysis shows the existence of 259 mutations in E protein. Interestingly, 16 of them were detected in the DFLV amino acid (aa) motif (aa72-aa75) that binds the host PALS1 protein. The results of root mean square deviation, fluctuations, radius of gyration, and free energy landscape show that Ser68Phe, Pro71Ser, and Leu73Phe are exhibiting a more stabilizing effect. However, a more comprehensive experimental study may be required to see the effect on virus pathogenicity. Potential antiviral drugs, and vaccines may be developed used after screening the genomic variations for better management of SARS-CoV-2 infections.
| S-EPMC8314890 | biostudies-literature
Project description:The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 (ACE2) binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work, and performed all-atom molecular dynamics simulations to get insight into the structural dynamics and interactions between S protein and antibodies. Investigation of the residues critical for S-antibody binding allows us to predict the potential influence of mutations in SARS-CoV-2 variants. Comparison of the glycan conformations between S-only and S-antibody systems reveals the roles of glycans in S-antibody binding. In addition, we explored the antibody binding modes, and the influences of antibody on the motion of S protein receptor binding domains. Overall, our analyses provide a better understanding of S-antibody interactions, and the simulation-based S-antibody interaction maps could be used to predict the influences of S mutation on S-antibody interactions, which will be useful for the development of vaccine and antibody-based therapy.
| S-EPMC8132224 | biostudies-literature
Project description:The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work and performed all-atom molecular dynamics simulations to gain insight into the structural dynamics and interactions between S protein and antibodies. Investigation of the residues critical for S-antibody binding allows us to predict the potential influence of mutations in SARS-CoV-2 variants. Comparison of the glycan conformations between S-only and S-antibody systems reveals the roles of glycans in S-antibody binding. In addition, we explored the antibody binding modes and the influences of antibody on the motion of S protein receptor binding domains. Overall, our analyses provide a better understanding of S-antibody interactions, and the simulation-based S-antibody interaction maps could be used to predict the influences of S mutation on S-antibody interactions, which will be useful for the development of vaccine and antibody-based therapy.
| S-EPMC8457324 | biostudies-literature
Project description:The current situation with the SARS-CoV-2 pandemic indicates the importance of new approaches in vaccine design. In order to design new attenuated vaccines, to decrease virulence of virus wild types, it is important to understand what allows a virus to hijack its host cell's metabolism, a property of all viruses. RNA and protein sequences obtained from databases were used to count the number of atoms of each element in the virions of SARS, MERS and SARS-CoV-2. The number of protein copies and carbohydrate composition were taken from the literature. The number of lipid molecules was estimated from the envelope surface area. Based on elemental composition, growth equations were balanced, and thermodynamic properties of the viruses were determined using Patel-Erickson and Battley equations. Elemental and molecular compositions of SARS, MERS and SARS-CoV-2 were found, as well as their standard thermodynamic properties of formation and growth. Standard Gibbs energy of growth of virus nucleocapsids was found to be significantly more negative than that of their host tissue. The ratio of Gibbs energies of growth of virus nucleocapsids and host cell is greater than unity. The more negative Gibbs energy of growth of viruses implies that virus multiplication has a greater driving force than synthesis of host cell components, giving a physical explanation of why viruses are able to hijack their host cell's metabolism. Knowing the mechanism of viral metabolism hijacking can open new paths for vaccine design. By manipulating chemical composition of viruses, virulence can be decreased by making the Gibbs energy of their growth less negative, resulting in decreased multiplication rate, while preserving antigenic properties.
| S-EPMC7489929 | biostudies-literature
Project description:The emergence of novel coronavirus mutants is a main factor behind the deterioration of the epidemic situation. Further studies into the pathogenicity of these mutants are thus urgently needed. Binding of the spinous protein receptor binding domain (RBD) of SARS-CoV-2 to the angiotensin-converting enzyme 2 (ACE2) receptor was shown to initiate coronavirus entry into host cells and lead to their infection. The receptor-binding motif (RBM, 438-506) is a region that directly interacts with ACE2 receptor in the RBD and plays a crucial role in determining affinity. To unravel how mutations in the non-RBM regions impact the interaction between RBD and ACE2, we selected three non-RBM mutant systems (N354D, D364Y, and V367F) from the documented clinical cases, and the Q498A mutant system located in the RBM region served as the control. Molecular dynamics simulation was conducted on the mutant systems and the wild-type (WT) system, and verified experiments also performed. Non-RBM mutations have been shown not only to change conformation of the RBM region but also to significantly influence its hydrogen bonding and hydrophobic interactions. In particular, the D364Y and V367F systems showed a higher affinity for ACE2 owing to their electrostatic interactions and polar solvation energy changes. In addition, although the binding free energy at this point increased after the mutation of N354D, the conformation of the random coil (Pro384-Asp389) was looser than that of other systems, and the combined effect weakened the binding free energy between RBD and ACE2. Interestingly, we also found a random coil (Ala475-Gly485). This random coil is very sensitive to mutations, and both types of mutations increase the binding free energy of residues in this region. We found that the binding loop (Tyr495-Tyr505) in the RBD domain strongly binds to Lys353, an important residue of the ACE2 domain previously identified. The binding free energy of the non-RBM mutant group at the binding loop had positive and negative changes, and these changes were more obvious than that of the Q498A system. The results of this study elucidate the effect of non-RBM mutation on ACE2-RBD binding, and provide new insights for SARS-CoV-2 mutation research.
| S-EPMC8353372 | biostudies-literature
Project description:pncB1 and pncB2 are two putative nicotinic acid phosphoribosyltransferases, playing a role in cofactor salvage and drug resistance in Mycobacterium tuberculosis. Mutations have been reported in first- and second-line drug targets, causing resistance. However, pncB1 and pncB2 mutational data are not available, and neither of their mutation effects have been investigated in protein structures. The current study has been designed to investigate mutations and also their effects on pncB1 and pncB2 structures. A total of 287 whole-genome sequenced data of drug-resistant Mycobacterium tuberculosis isolates from Khyber Pakhtunkhwa of Pakistan were retrieved (BioSample PRJEB32684, ERR2510337-ERR2510445, ERR2510546-ERR2510645) from NCBI. The genomic data were analyzed for pncB1 and pncB2 mutations using PhyResSE. All the samples harbored numerous synonymous and non-synonymous mutations in pncB1 and pncB2 except one. Mutations Pro447Ser, Arg286Arg, Gly127Ser, and delTCAGGCCG1499213>1499220 in pncB1 are novel and have not been reported in literature and TB databases. The most common non-synonymous mutations exhibited stabilizing effects on the pncB1 structure. Moreover, 36 out of 287 samples harbored two non-synonymous and 34 synonymous mutations in pncB2 among which the most common was Phe204Phe (TTT/TTC), present in 8 samples, which may have an important effect on the usage of specific codons that may increase the gene expression level or protein folding effect. Mutations Ser120Leu and Pro447Ser, which are present in the loop region, exhibited a gain in flexibility in the surrounding residues while Gly429Ala and Gly127Ser also demonstrated stabilizing effects on the protein structure. Inhibitors designed based on the most common pncB1 and pncB2 mutants may be a more useful strategy in high-burden countries. More studies are needed to elucidate the effect of synonymous mutations on organism phenotype.
| S-EPMC9031469 | biostudies-literature
Project description:The human spike protein sequences from Asia, Africa, Europe, North America, South America, and Oceania were analyzed by comparing with the reference severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) protein sequence from Wuhan-Hu-1, China. Out of 10333 spike protein sequences analyzed, 8155 proteins comprised one or more mutations. A total of 9654 mutations were observed that correspond to 400 distinct mutation sites. The receptor binding domain (RBD) which is involved in the interactions with human angiotensin-converting enzyme-2 (ACE-2) receptor and causes infection leading to the COVID-19 disease comprised 44 mutations that included residues within 3.2 Å interacting distance from the ACE-2 receptor. The mutations observed in the spike proteins are discussed in the context of their distribution according to the geographical locations, mutation sites, mutation types, distribution of the number of mutations at the mutation sites and mutations at the glycosylation sites. The density of mutations in different regions of the spike protein sequence and location of the mutations in protein three-dimensional structure corresponding to the RBD are discussed. The mutations identified in the present work are important considerations for antibody, vaccine, and drug development.
| S-EPMC8014176 | biostudies-literature
Project description:The SARS-CoV-2 virus and its mutations have affected human health globally and created significant danger for the health of people all around the world. To cure this virus, the human Angiotensin Converting Enzyme-2 (ACE2) receptor, the SARS-CoV-2 main protease (Mpro), and spike proteins were found to be likely candidates for the synthesis of novel therapeutic drug. In the past, proteins were capable of engaging in interaction with a wide variety of ligands, including both manmade and plant-derived small molecules. Pyrus communis L., Ginko bibola, Carica papaya, Syrian rue, and Pimenta dioica were some of the plant species that were studied for their tendency to interact with SARS-CoV-2 main protease (Mpro) in this research project (6LU7). This scenario investigates the geometry, electronic, and thermodynamic properties computationally. Assessing the intermolecular forces of phytochemicals with the targets of the SARS-CoV-2 Mpro spike protein (SP) resulted in the recognition of a compound, kaempferol, as the most potent binding ligand, -7.7 kcal mol-1. Kaempferol interacted with ASP-187, CYS-145, SER-144, LEU 141, MET-165, and GLU-166 residues. Through additional molecular dynamic simulations, the stability of ligand-protein interactions was assessed for 100 ns. GLU-166 remained intact with 33% contact strength with phenolic OH group. We noted a change in torsional conformation, and the molecular dynamics simulation showed a potential variation in the range from 3.36 to 7.44 against a 45-50-degree angle rotation. SAR, pharmacokinetics, and drug-likeness characteristic investigations showed that kaempferol may be the suitable candidate to serve as a model for designing and developing new anti-COVID-19 medicines.
| S-EPMC10934104 | biostudies-literature
Project description:Vaccination against SARS-CoV-2 just started in most of the countries. However, the development of specific vaccines against SARS-CoV-2 is not the only approach to control the virus and monoclonal antibodies (mAbs) start to merit special attention as a therapeutic option to treat COVID-19 disease. Here, the main conformations and interactions between the receptor-binding domain (RBD) of spike glycoprotein of SARS-CoV-2 (S protein) with two mAbs (CR3022 and S309) and the ACE2 cell receptor are studied as the main representatives of three different epitopes on the RBD of S protein. The combined approach of 1 μs accelerated molecular dynamics (aMD) and ab-initio hybrid molecular dynamics is used to identify the most predominant interactions under physiological conditions. Results allow to determine the main receptor-binding mapping, hydrogen bonding network and salt bridges in the most populated antigen-antibody interface conformations. The deep knowledge on the protein-protein interactions involving mAbs and ACE2 receptor with the spike glycoprotein of SARS-CoV-2 increases background knowledge to speed up the development of new vaccines and therapeutic drugs.
| S-EPMC8545699 | biostudies-literature
Project description: Not available
| EMPIAR-11168 | biostudies-other