Molecular dynamics simulation revealed binding of nucleotide inhibitors to ZIKV polymerase over 444 nanoseconds.
Ontology highlight
ABSTRACT:
SUBMITTER: Elfiky AA
PROVIDER: S-EPMC7166584 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Influenza A virus is the main cause of worldwide epidemics and annual influenza outbreaks in humans. In this study, a virtual screen was performed to identify compounds that interact with the PB2 cap-binding domain (CBD) of influenza A polymerase. A virtual screening workflow based on Glide docking was used to screen an internal database containing 8417 molecules, and then the output compounds were selected based on solubility, absorbance, and structural fingerprints. Of the 16 compounds selected for biological evaluation, six compounds were identified that rescued cells from H1N1 virus-mediated death at non-cytotoxic concentrations, with EC50 values ranging from 2.5-55.43 μM, and that could bind to the PB2 CBD of H1N1, with Kd values ranging from 0.081-1.53 μM. Molecular dynamics (MD) simulations of the docking complexes of our active compounds revealed that each compound had its own binding characteristics that differed from those of VX-787. Our active compounds have novel structures and unique binding modes with PB2 proteins, and are suitable to serve as lead compounds for the development of PB2 inhibitors. An analysis of the MD simulation also helped us to identify the dominant amino acid residues that play a key role in binding the ligand to PB2, suggesting that we should focus on increasing and enhancing the interaction between inhibitors and these major amino acids during lead compound optimization to obtain more active PB2 inhibitors.
| S-EPMC8623395 | biostudies-literature
Project description:The RNA-dependent RNA polymerase (RdRp) of norovirus is an attractive target of antiviral agents aimed at providing protection against norovirus-associated gastroenteritis. Here, we perform molecular dynamics simulations of the crystal structure of norovirus RdRp in complex with several known binders, as well as free-energy simulations by free-energy perturbation (FEP) to determine binding free energies of these molecules relative to the natural nucleotide substrates. We determine experimental EC50 values and nucleotide incorporation efficiencies for several of these compounds. Moreover, we investigate the mechanism of inhibition of some of these ligands. Using FEP, we screened a virtual nucleotide library with 121 elements for binding to the polymerase and successfully identified two novel chain terminators.
| S-EPMC7869559 | biostudies-literature
Project description:SARS-CoV-2, a member of the coronavirus family, has caused a global public health emergency. Based on our analysis of hepatitis C virus and coronavirus replication, and the molecular structures and activities of viral inhibitors, we previously reasoned that the FDA-approved hepatitis C drug EPCLUSA (Sofosbuvir/Velpatasvir) should inhibit coronaviruses, including SARS-CoV-2. Here, using model polymerase extension experiments, we demonstrate that the active triphosphate form of Sofosbuvir is incorporated by low-fidelity polymerases and SARS-CoV RNA-dependent RNA polymerase (RdRp), and blocks further incorporation by these polymerases; the active triphosphate form of Sofosbuvir is not incorporated by a host-like high-fidelity DNA polymerase. Using the same molecular insight, we selected 3'-fluoro-3'-deoxythymidine triphosphate and 3'-azido-3'-deoxythymidine triphosphate, which are the active forms of two other anti-viral agents, Alovudine and AZT (an FDA-approved HIV/AIDS drug) for evaluation as inhibitors of SARS-CoV RdRp. We demonstrate the ability of two of these HIV reverse transcriptase inhibitors to be incorporated by SARS-CoV RdRp where they also terminate further polymerase extension. Given the 98% amino acid similarity of the SARS-CoV and SARS-CoV-2 RdRps, we expect these nucleotide analogues would also inhibit the SARS-CoV-2 polymerase. These results offer guidance to further modify these nucleotide analogues to generate more potent broad-spectrum anti-coronavirus agents.
| S-EPMC7596664 | biostudies-literature
Project description:Zika virus (ZIKV) was identified in 1947 in the Zika forest of Uganda and it has emerged recently as a global health threat, with recurring outbreaks and its associations with congenital microcephaly through maternal fetal transmission and Guillain-Barré syndrome. Currently, there are no United States (US) Food and Drug Administration (FDA)-approved vaccines or antivirals to treat ZIKV infections, which underscores an urgent medical need for the development of disease intervention strategies to treat ZIKV infection and associated disease. Drug repurposing offers various advantages over developing an entirely new drug by significantly reducing the timeline and resources required to advance a candidate antiviral into the clinic. Screening the ReFRAME library, we identified ten compounds with antiviral activity against the prototypic mammarenavirus lymphocytic choriomeningitis virus (LCMV). Moreover, we showed the ability of these ten compounds to inhibit influenza A and B virus infections, supporting their broad-spectrum antiviral activity. In this study, we further evaluated the broad-spectrum antiviral activity of the ten identified compounds by testing their activity against ZIKV. Among the ten compounds, Azaribine (SI-MTT = 146.29), AVN-944 (SI-MTT = 278.16), and Brequinar (SI-MTT = 157.42) showed potent anti-ZIKV activity in post-treatment therapeutic conditions. We also observed potent anti-ZIKV activity for Mycophenolate mofetil (SI-MTT = 20.51), Mycophenolic acid (SI-MTT = 36.33), and AVN-944 (SI-MTT = 24.51) in pre-treatment prophylactic conditions and potent co-treatment inhibitory activity for Obatoclax (SI-MTT = 60.58), Azaribine (SI-MTT = 91.51), and Mycophenolate mofetil (SI-MTT = 73.26) in co-treatment conditions. Importantly, the inhibitory effect of these compounds was strain independent, as they similarly inhibited ZIKV strains from both African and Asian/American lineages. Our results support the broad-spectrum antiviral activity of these ten compounds and suggest their use for the development of antiviral treatment options of ZIKV infection.
| S-EPMC7551609 | biostudies-literature
Project description:Microtubules (MTs) are the largest type of cellular filament, essential in processes ranging from mitosis and meiosis to flagellar motility. Many of the processes depend critically on the mechanical properties of the MT, but the elastic moduli, notably the Young's modulus, are not directly revealed in experiment, which instead measures either flexural rigidity or response to radial deformation. Molecular dynamics (MD) is a method that allows the mechanical properties of single biomolecules to be investigated through computation. Typically, MD requires an atomic resolution structure of the molecule, which is unavailable for many systems, including MTs. By combining structural information from cryo-electron microscopy and electron crystallography, we have constructed an all-atom model of a complete MT and used MD to determine its mechanical properties. The simulations revealed nonlinear axial stress-strain behavior featuring a pronounced softening under extension, a possible plastic deformation transition under radial compression, and a distinct asymmetry in response to the two senses of twist. This work demonstrates the possibility of combining different levels of structural information to produce all-atom models suitable for quantitative MD simulations, which extends the range of systems amenable to the MD method and should enable exciting advances in our microscopic knowledge of biology.
| S-EPMC2905083 | biostudies-literature
Project description:N-Myristoyltransferase (NMT) is a cytosolic monomeric enzyme involved in the allocation of the myristoyl group to the aminoterminal of glycine in several viral and eukaryotic cellular proteins. NMT has been validated as a potential drug target against kinetoplastid for parasitic protozoa. A multistep virtual screening protocol based on the pharmacophore modeling, molecular docking, and molecular dynamics simulation was carried out. Initially, Maybridge database was virtually screened via a validated pharmacophore model. The effective pharmacophore models were accompanied with exclusion volumes to improve their receiver operating characteristic curve to identify potential NMT inhibitors. The hits identified as actives based on the 3D-pharmacophore model were evaluated by molecular docking studies. In stepwise screening, six compounds were shortlisted for the dynamic simulation to get insights into their binding mode. In conclusion, this study provides fundamental information about the architecture of the binding site and some crucial residues that may provide insights into the development of new antiparasitic agents.
| S-EPMC6714517 | biostudies-literature
Project description:Amyloid peptides are known to self-assemble into larger aggregates that are linked to the pathogenesis of many neurodegenerative disorders. In contrast to primary nucleation, recent experimental and theoretical studies have shown that many toxic oligomeric species are generated through secondary processes on a pre-existing fibrillar surface. Nucleation, for example, can also occur along the surface of a pre-existing fibril-secondary nucleation-as opposed to the primary one. However, explicit pathways are still not clear. In this study, we use molecular dynamics simulation to explore the free energy landscape of a free Abeta monomer binding to an existing fibrillar surface. We specifically look into several potential Abeta structural precursors that might precede some secondary events, including elongation and secondary nucleation. We find that the overall process of surface-dependent events can be described at least by the following three stages: 1. Free diffusion 2. Downhill guiding 3. Dock and lock. And we show that the outcome of adding a new monomer onto a pre-existing fibril is pathway-dependent, which leads to different secondary processes. To understand structural details, we have identified several monomeric amyloid precursors over the fibrillar surfaces and characterize their heterogeneity using a probability contact map analysis. Using the frustration analysis (a bioinformatics tool), we show that surface heterogeneity correlates with the energy frustration of specific local residues that form binding sites on the fibrillar structure. We further investigate the helical twisting of protofilaments of different sizes and observe a length dependence on the filament twisting. This work presents a comprehensive survey over the properties of fibril growth using a combination of several openMM-based platforms, including the GPU-enabled openAWSEM package for coarse-grained modeling, MDTraj for trajectory analysis, and pyEMMA for free energy calculation. This combined approach makes long-timescale simulation for aggregation systems as well as all-in-one analysis feasible. We show that this protocol allows us to explore fibril stability, surface binding affinity/heterogeneity, as well as fibrillar twisting. All these properties are important for understanding the molecular mechanism of surface-catalyzed secondary processes of fibril growth.
| S-EPMC8378332 | biostudies-literature
Project description:A multiscale simulation of a complex between the lac repressor protein (LacI) and a 107-bp-long DNA segment is reported. The complex between the repressor and two operator DNA segments is described by all-atom molecular dynamics; the size of the simulated system comprises either 226,000 or 314,000 atoms. The DNA loop connecting the operators is modeled as a continuous elastic ribbon, described mathematically by the nonlinear Kirchhoff differential equations with boundary conditions obtained from the coordinates of the terminal base pairs of each operator. The forces stemming from the looped DNA are included in the molecular dynamics simulations; the loop structure and the forces are continuously recomputed because the protein motions during the simulations shift the operators and the presumed termini of the loop. The simulations reveal the structural dynamics of the LacI-DNA complex in unprecedented detail. The multiple domains of LacI exhibit remarkable structural stability during the simulation, moving much like rigid bodies. LacI is shown to absorb the strain from the looped DNA mainly through its mobile DNA-binding head groups. Even with large fluctuating forces applied, the head groups tilt strongly and keep their grip on the operator DNA, while the remainder of the protein retains its V-shaped structure. A simulated opening of the cleft of LacI by 500-pN forces revealed the interactions responsible for locking LacI in the V-conformation.
| S-EPMC1100768 | biostudies-literature
Project description:[This corrects the article DOI: 10.3389/fmolb.2021.719320.].
| S-EPMC9252589 | biostudies-literature
Project description:Bisnaphthalimide intercalators are anti-tumour agents composed of two planar rings linked by a flexible diazanonylene chain. The intercalated rings of three bisnaphthalimide analogues complexed to DNA are found here to undergo 180 degrees rotating motions that do not affect the diazanonylene linker atoms bound to the major groove. These ring rotations are detected by NMR spectroscopy in a broad range of sequence contexts and duplex lengths. A comparative analysis of the frequency and activation energies of such excited states in different complexes and conditions indicates that these motions (i) are unrelated to drug dissociation; (ii) are a consequence of concerted, sequence-dependent nucleotide movements taking place on the millisecond time scale; and (iii) may occur inside the DNA duplexes. The rotation frequencies range from 2 to 25 s(-1) at 25 degrees C, depending on DNA composition and the size of the rotating rings. The detected nucleotide dynamics are likely to play an important role in the binding kinetics of the numerous proteins and drugs that require base unstacking when interacting with DNA.
| S-EPMC484180 | biostudies-literature