Molecular dynamics simulations of the flexibility and inhibition of SARS-CoV-2 NSP 13 helicase.
Ontology highlight
ABSTRACT:
SUBMITTER: Raubenolt BA
PROVIDER: S-EPMC8730789 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:SARS-CoV-2 is accountable for severe social and economic disruption around the world causing COVID-19. Non-structural protein-15 (NSP15) possesses a domain that is vital to the viral life cycle and is known as uridylate-specific endoribonuclease (EndoU). This domain binds to the uridine 5'-monophosphate (U5P) so that the protein may carry out its native activity. It is considered a vital drug target to inhibit the growth of the virus. Thus, in this current study, ML-based QSAR and virtual screening of U5P analogues targeting Nsp15 were performed to identify potential molecules against SARS-CoV-2. Screening of 816 unique U5P analogues using ML-based QSAR identified 397 compounds ranked on their predicted bioactivity (pIC50). Further, molecular docking and hydrogen bond interaction analysis resulted in the selection of the top three compounds (53309102, 57398422, and 76314921). Molecular dynamics simulation of the most promising compounds showed that two molecules 53309102 and 57398422 acted as potential binders of Nsp15. The compound was able to inhibit nsp15 activity as it was successfully bound to the active site of the nsp15 protein. This was achieved by the formation of relevant contacts with enzymatically critical amino acid residues (His235, His250, and Lys290). Principal component analysis and free energy landscape studies showed stable complex formation while MM/GBSA calculation showed lower binding energies for 53309102 (ΔGTOTAL = -29.4 kcal/mol) and 57398422 (ΔGTOTAL = -39.4 kcal/mol) compared to the control U5P (ΔGTOTAL = -18.8 kcal/mol). This study aimed to identify analogues of U5P inhibiting the NSP15 function that potentially could be used for treating COVID-19.
| S-EPMC10727945 | biostudies-literature
Project description:Helicases are ubiquitous enzymes involved in nucleic acid metabolism. The PcrA DNA helicase is an essential bacterial protein involved in rolling circle plasmid replication and DNA repair. Recent crystal structures of PcrA bound to DNA indicate that a flexible loop mediates a functionally important rigid-body-domain rotation. In this study, we report stochastic boundary molecular dynamics simulations focused on this region for wild-type and mutants designed to increase the rigidity of the region. Residues in the region that were helix-disfavoring, such as glycine, threonine, and others, were mutated to alanine. The simulated dynamics, analyzed with a variety of measures of structure and mobility, indicate that a few point mutations will substantially increase helix formation in this region. Subnanosecond stochastic boundary molecular dynamics simulations at several temperatures offer a rapid protocol for assessing large numbers of mutants and provides a novel strategy for the design of experiments to test the role of this flexible loop region in the function of PcrA.
| S-EPMC3034217 | biostudies-literature
Project description:The SARS coronavirus main proteinase (M(pro)) is a key enzyme in the processing of the viral polyproteins and thus an attractive target for the discovery of drugs directed against SARS. The enzyme has been shown by X-ray crystallography to undergo significant pH-dependent conformational changes. Here, we assess the conformational flexibility of the M(pro) by analysis of multiple crystal structures (including two new crystal forms) and by molecular dynamics (MD) calculations. The MD simulations take into account the different protonation states of two histidine residues in the substrate-binding site and explain the pH-activity profile of the enzyme. The low enzymatic activity of the M(pro) monomer and the need for dimerization are also discussed.
| S-EPMC7094468 | biostudies-literature
Project description:The outbreak of COVID-19 caused by SARS-CoV-2 virus continually led to infect a large population worldwide. Currently, there is no specific viral protein-targeted therapeutics. The Nucleocapsid (N) protein of the SARS-CoV-2 virus is necessary for viral RNA replication and transcription. The C-terminal domain of N protein (CTD) involves in the self-assembly of N protein into a filament that is packaged into new virions. In this study, the CTD (PDB ID: 6WJI) was targeted for the identification of possible inhibitors of oligomerization of N protein. Herein, multiple computational approaches were employed to explore the potential mechanisms of binding and inhibitor activity of five antiviral drugs toward CTD. The five anti-N drugs studied in this work are 4E1RCat, Silmitasertib, TMCB, Sapanisertib, and Rapamycin. Among the five drugs, 4E1RCat displayed highest binding affinity (-10.95 kcal/mol), followed by rapamycin (-8.91 kcal/mol), silmitasertib (-7.89 kcal/mol), TMCB (-7.05 kcal/mol), and sapanisertib (-6.14 kcal/mol). Subsequently, stability and dynamics of the protein-drug complex were examined with molecular dynamics (MD) simulations. Overall, drug binding increases the stability of the complex with maximum stability observed in the case of 4E1RCat. The CTD-drug complex systems behave differently in terms of the free energy landscape and showed differences in population distribution. Overall, the MD simulation parameters like RMSD, RMSF, Rg, hydrogen bonds analysis, PCA, FEL, and DCCM analysis indicated that 4E1RCat and TMCB complexes were more stable as compared to silmitasertib and sapanisertib and thus could act as effective drug compounds against CTD.Communicated by Ramaswamy H. Sarma.
| S-EPMC7663461 | biostudies-literature
Project description:The functional role of the highly conserved stem-loop II motif (s2m) in SARS-CoV and SARS-CoV-2 in the viral lifecycle remains enigmatic and an intense area of research. Structure and dynamics of the s2m are key to establishing a structure-function connection, yet a full set of atomistic resolution coordinates is not available for SARS-CoV-2 s2m. Our work constructs three-dimensional coordinates consistent with NMR solution phase data for SARS-CoV-2 s2m and provides a comparative analysis with its counterpart SARS-CoV s2m. We employed initial coordinates based on PDB ID 1XJR for SARS-CoV s2m and two models for SARS-CoV-2 s2m: one based on 1XJR in which we introduced the mutations present in SARS-CoV-2 s2m and the second based on the available SARS-CoV-2 NMR NOE data supplemented with knowledge-based methods. For each of the three systems, 3.5 μs molecular dynamics simulations were used to sample the structure and dynamics, and principal component analysis (PCA) reduced the ensembles to hierarchal conformational substates for detailed analysis. Dilute solution simulations of SARS-CoV s2m demonstrate that the GNRA-like terminal pentaloop is rigidly defined by base stacking uniquely positioned for possible kissing dimer formation. However, the SARS-CoV-2 s2m simulation did not retain the reported crystallographic SARS-CoV motifs and the terminal loop expands to a highly dynamic "nonaloop." Increased flexibility and structural disorganization are observed for the larger terminal loop, where an entropic penalty is computed to explain the experimentally observed reduction in kissing complex formation. Overall, both SARS-CoV and SARS-CoV-2 s2m elements have a similarly pronounced L-shape due to different motif interactions. Our study establishes the atomistic three-dimensional structure and uncovers dynamic differences that arise from s2m sequence changes, which sets the stage for the interrogation of different mechanistic pathways of suspected biological function.
| S-EPMC9578647 | biostudies-literature
Project description:Novel coronavirus SARS-CoV-2 has infected millions of people with thousands of mortalities globally. The main protease (Mpro) is vital in processing replicase polyproteins. Both the CoV's Mpro shares 97% identity, with 12 mutations, but none are present in the active site. Although many therapeutics and vaccines are available to combat SARS-CoV-2, these treatments may not be practical due to their high mutational rate. On the other hand, Mpro has a high degree of conservation throughout variants, making Mpro a stout drug target. Here, we report a detailed comparison of both the monomeric Mpro and the biologically active dimeric Mpro using MD simulation to understand the impact of the 12 divergent residues (T35V, A46S, S65N, L86V, R88K, S94A, H134F, K180N, L202V, A267S, T285A and I286L) on the molecular microenvironment and the interaction between crucial residues. The present study concluded that the change in the microenvironment of residues at the entrance (T25, T26, M49 and Q189), near the catalytic site (F140, H163, H164, M165 and H172) and in the substrate-binding site (V35, N65, K88 and N180) is due to 12 mutations in the SARS-CoV-2 Mpro. Furthermore, the involvement of F140, E166 and H172 residues in dimerization stabilizes the Mpro dimer, which should be considered. We anticipate that networks and microenvironment changes identified here might guide repurposing attempts and optimization of new Mpro inhibitors.Supplementary informationThe online version contains supplementary material available at 10.1007/s11224-022-02089-6.
| S-EPMC9686461 | biostudies-literature
Project description:SARS-CoV-2 has led to a global pandemic of new crown pneumonia, which has had a tremendous impact on human society. Antibody drug therapy is one of the most effective way of combating SARS-CoV-2. In order to design potential antibody drugs with high affinity, we used antibody S309 from patients with SARS-CoV as the target antibody and RBD of S protein as the target antigen. Systems with RBD glycosylated and non-glycosylated were constructed to study the influence of glycosylation. From the results of molecular dynamics simulations, the steric effects of glycans on the surface of RBD plays a role of "wedge", which makes the L335-E340 region of RBD close to the CDR3 region of the heavy chain of antibody and increases the contact area between antigen and antibody. By mutating the key residues of antibody at the interaction interface, we found that the binding affinities of antibody mutants G103A, P28W and Y100W were all stronger than that of the wild-type, especially for the G103A mutant. G103A significantly reduces the distance between the binding region of L335-K356 in the antigen and P28-Y32 of heavy chain in the antibody through structural transition. Taken together, the antibody design method described in this work can provide theoretical guidance and a time-saving method for antibody drug design.
| S-EPMC9820416 | biostudies-literature
Project description:ORF8a protein is 39 residues long and contains a single transmembrane domain. The protein is synthesized using solid phase peptide synthesis and reconstituted into artificial lipid bilayers that forms cation-selective ion channels with a main conductance level of 8.9±0.8pS at elevated temperature (38.5°C). Computational modeling studies including multi nanosecond molecular dynamics simulations in a hydrated POPC lipid bilayer are done with a 22 amino acid transmembrane helix to predict a putative homooligomeric helical bundle model. A structural model of a pentameric bundle is proposed with cysteines, serines and threonines facing the pore.
| S-EPMC7094593 | biostudies-literature
Project description:At present, the most powerful new drugs for COVID-19 are antibody proteins. In addition, there are some star small molecule drugs. However, there are few studies on nanomaterials. Here, we study the intact graphene (IG), defective graphene (DG), and graphene oxide (GO) interacting with COVID-19 protein. We find that they show progressive inhibition of COVID-19 protein. By using molecular dynamics simulations, we study the interactions between SARS-CoV-2 3CL Mpro and graphene-related materials (GRMs): IG, DG, and GO. The results show that Mpro can be absorbed onto the surfaces of investigated materials. DG and GO interacted with Mpro more intensely, causing the decisive part of Mpro to become more flexible. Further analysis shows that compared to IG and GO, DG can inactivate Mpro and inhibit its expression effectively by destroying the active pocket of Mpro. Our work not only provides detailed and reliable theoretical guidance for the application of GRMs in treating with SARS-CoV-2 but also helps in developing new graphene-based anti-COVID-19 materials.
| S-EPMC8713398 | biostudies-literature
Project description:Flexibility and dynamics are protein characteristics that are essential for the process of molecular recognition. Conformational changes in the protein that are coupled to ligand binding are described by the biophysical models of induced fit and conformational selection. Different concepts that incorporate protein flexibility into protein-ligand docking within the context of these two models are reviewed. Several computational studies that discuss the validity and possible limitations of such approaches will be presented. Finally, different approaches that incorporate protein dynamics, e.g., configurational entropy, and solvation effects into docking will be highlighted.
| S-EPMC3172316 | biostudies-literature