scholarly journals Molecular Dynamics Simulation of Cholera Toxin A-1 Polypeptide

2016 ◽  
Vol 14 (1) ◽  
pp. 188-196 ◽  
Author(s):  
Syed Lal Badshah ◽  
Abdul Naeem Khan ◽  
Yahia Nasser Mabkhot

AbstractA molecular dynamics (MD) simulation study of the enzymatic portion of cholera toxin; cholera toxin A-1 polypeptide (CTA1) was performed at 283, 310 and 323 K. From total energy analysis it was observed that this toxin is stable thermodynamically and these outcomes were likewise confirmed by root mean square deviations (RMSD) investigations. The Cα root mean square fluctuation (RMSF) examinations revealed that there are a number of residues inside CTA1, which can be used as target for designing and synthesizing inhibitory drugs, in order to inactivate cholera toxin inside the human body. The fluctuations in the radius of gyration and hydrogen bonding in CTA1 proved that protein unfolding and refolding were normal routine phenomena in its structure at all temperatures. Solvent accessible surface area study identified the hydrophilic nature of the CTA1, and due to this property it can be a potential biological weapon. The structural identification (STRIDE) algorithm for proteins was successfully used to determine the partially disordered secondary structure of CTA1. On account of this partially disordered secondary structure, it can easily deceive the proteolytic enzymes of the endoplasmic reticulum of host cells.

Human insulin, a small protein hormone consisting of A-chain (21 residues) and B-chain (30 residues) linked by three disulfide bonds, is crucial for controlling the hyperglycemia in type I diabetes. In the present work molecular dynamics simulation (MD) with human insulin and its mutants was used to assess the influence of 10 point mutations (HisA8, ValA10, AspB10, GlnB17, AlaB17, GlnB18, AspB25, ThrB26, GluB27, AspB28), 6 double mutations (GluA13+GluB10, SerA13+GluB27, GluB1+GluB27, SerB2+AspB10, AspB9+GluB27, GluB16+GluB27) and one triple mutation (GluA15+AspA18+AspB3) in the protein sequence on the structure and dynamics of human insulin. A series of thermal unfolding MD simulations with wild type (WT) human insulin and its mutants was performed at 400 K with GROMACS software (version 5.1) using the CHARMM36m force field. The MD results have been analyzed in terms of the parameters characterizing both the global and local protein structure, such as the backbone root mean-square deviation, gyration radius, solvent accessible surface area, the root mean-square fluctuations and the secondary structure content. The MD simulation data showed that depending on time evolution of integral characteristics, the examined mutants can be tentatively divided into three groups: 1) the mutants HisA8, ValA10, AlaB17, AspB25, ThrB26, GluB27, GluA13+GluB10, GluB1+GluB27 and GluB16+GluB27, which exert stabilizing effect on the protein structure in comparison with wild type insulin; 2) the mutants GlnB17, AspB10, SerB2+AspB10 and GluA15+AspA18+AspB3 that did not significantly affect the dynamical properties of human insulin with a minimal stabilizing impact; 3) the mutants AspB28, AspB9+GluB27 and SerA13+GluB27, GlnB18, destabilizing the protein structure. Analysis of the secondary structure content provided evidence for the influence of AspB28, AspB9+GluB27 and SerA13+GluB27, GlnB18 on the insulin unfolding. Our MD results indicate that the replacement of superficial nonpolar residues in the insulin structure by hydrophilic ones gives rise to the increase in protein stability in comparison with the wild type protein.


2021 ◽  
Vol 1 ◽  
Author(s):  
Shafi Mahmud ◽  
Md. Robiul Hasan ◽  
Suvro Biswas ◽  
Gobindo Kumar Paul ◽  
Shamima Afrose ◽  
...  

Coronavirus disease 2019 (COVID-19) is a potentially lethal and devastating disease that has quickly become a public health threat worldwide. Due to its high transmission rate, many countries were forced to implement lockdown protocols, wreaking havoc on the global economy and the medical crisis. The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative virus for COVID-19, represent an effective target for the development of a new drug/vaccine because it is well-conserved and plays a vital role in viral replication. Mpro inhibition can stop the replication, transcription as well as recombination of SARS-CoV-2 after the infection and thus can halt the formation of virus particles, making Mpro a viable therapeutic target. Here, we constructed a phytochemical dataset based on a rigorous literature review and explored the probability that various phytochemicals will bind with the main protease using a molecular docking approach. The top three hit compounds, medicagol, faradiol, and flavanthrin, had binding scores of −8.3, −8.6, and −8.8 kcal/mol, respectively, in the docking analysis. These three compounds bind to the active groove, consisting of His41, Cys45, Met165, Met49, Gln189, Thr24, and Thr190, resulting in main protease inhibition. Moreover, the multiple descriptors from the molecular dynamics simulation, including the root-mean-square deviation, root-mean-square fluctuation, solvent-accessible surface area, radius of gyration, and hydrogen bond analysis, confirmed the stable nature of the docked complexes. In addition, absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis confirmed a lack of toxicity or carcinogenicity for the screened compounds. Our computational analysis may contribute toward the design of an effective drug against the main protease of SARS-CoV-2.


Molecules ◽  
2021 ◽  
Vol 26 (21) ◽  
pp. 6538
Author(s):  
Md. Mominur Rahman ◽  
Md. Junaid ◽  
S. M. Zahid Hosen ◽  
Mohammad Mostafa ◽  
Lei Liu ◽  
...  

Inflammation plays an important role in different chronic diseases. Brominated indoles derived from the Australian marine mollusk Dicathais orbita (D. orbita) are of interest for their anti-inflammatory properties. This study evaluates the binding mechanism and potentiality of several brominated indoles (tyrindoxyl sulfate, tyrindoleninone, 6-bromoisatin, and 6,6′-dibromoindirubin) against inflammatory mediators cyclooxygenases-1/2 (COX-1/2) using molecular docking, followed by molecular dynamics simulation, along with physicochemical, drug-likeness, pharmacokinetic (pk), and toxicokinetic (tk) properties. Molecular docking identified that these indole compounds are anchored, with the main amino acid residues, positioned in the binding pocket of the COX-1/2, required for selective inhibition. Moreover, the molecular dynamics simulation based on root mean square deviation (RMSD), radius of gyration (Rg), solvent accessible surface area (SASA), and root mean square fluctuation (RMSF) analyses showed that these natural brominated molecules transit rapidly to a progressive constant configuration during binding with COX-1/2 and seem to accomplish a consistent dynamic behavior by maintaining conformational stability and compactness. The results were comparable to the Food and Drug Administration (FDA)-approved selective COX inhibitor, aspirin. Furthermore, the free energy of binding for the compounds assessed by molecular mechanics–Poisson–Boltzmann surface area (MM–PBSA) confirmed the binding capacity of indoles towards COX-1/2, with suitable binding energy values except for the polar precursor tyrindoxyl sulfate (with COX-1). The physicochemical and drug-likeness analysis showed zero violations of Lipinski’s rule, and the compounds are predicted to have excellent pharmacokinetic profiles. These indoles are projected to be non-mutagenic and free from hepatotoxicity, with no inhibition of human ether-a-go–go gene (hERG) I inhibitors, and the oral acute toxicity LD50 in rats is predicted to be similar or lower than aspirin. Overall, this work has identified a plausible mechanism for selective COX inhibition by natural marine indoles as potential therapeutic candidates for the mitigation of inflammation.


2020 ◽  
Author(s):  
Ajita Pandey

Abstract Novel coronavirus disease 2019 (COVID-19) is a highly infectious disease that is caused by the recently discovered severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Because there are no specific vaccines or drugs for SARS-CoV-2, drug repurposing may be a promising approach. SARS-CoV-2 has a positive-sense RNA genome that encodes non-structural proteins (Nsps), which are essential for viral replication in the host cell. Non-structural protein 3 (Nsp3) is a multidomain protein and is the largest protein of the replicase complex. Nsp3 contains an ADP-ribose phosphatase (ADRP) domain, also called the macrodomain, which interferes with the host immune response. In the present study, we used computational regression methods to target the ADRP domain of Nsp3, using FDA-approved drugs. We virtually screened 2,892 FDA-approved drugs, using a combination of molecular docking and scoring functions. Saquinavir and trifluperidol were identified as potential leads and were further investigated using molecular dynamics simulation (MDS) to predict the stability and behavior of the ADRP-drug complexes. Analysis of root mean square deviation, root mean square fluctuation, radius of gyration, solvent accessible surface area and number of hydrogen bonds showed that the ADRP-trifluperidol complex is more stable than the ADRP-saquinavir complex. The screening and the MDS results suggest that trifluperidol is a novel inhibitor of the ADRP domain of Nsp3. Trifluperidol could, therefore, potentially be used to help control the spread of COVID-19, either alone or in combination with antiviral agents. Further in-vitro and in-vivo experiments are necessary to confirm our in silico results.


Molecules ◽  
2021 ◽  
Vol 26 (8) ◽  
pp. 2210
Author(s):  
Shafi Mahmud ◽  
Gobindo Kumar Paul ◽  
Mirola Afroze ◽  
Shirmin Islam ◽  
Swagota Briti Ray Gupt ◽  
...  

The recent coronavirus disease 2019 (COVID-19) pandemic is a global threat for healthcare management and the economic system, and effective treatments against the pathogenic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus responsible for this disease have not yet progressed beyond the developmental phases. As drug refinement and vaccine progression require enormously broad investments of time, alternative strategies are urgently needed. In this study, we examined phytochemicals extracted from Avicennia officinalis and evaluated their potential effects against the main protease of SARS-CoV-2. The antioxidant activities of A. officinalis leaf and fruit extracts at 150 µg/mL were 95.97% and 92.48%, respectively. Furthermore, both extracts displayed low cytotoxicity levels against Artemia salina. The gas chromatography–mass spectroscopy analysis confirmed the identifies of 75 phytochemicals from both extracts, and four potent compounds, triacontane, hexacosane, methyl linoleate, and methyl palminoleate, had binding free energy values of −6.75, −6.7, −6.3, and −6.3 Kcal/mol, respectively, in complexes with the SARS-CoV-2 main protease. The active residues Cys145, Met165, Glu166, Gln189, and Arg188 in the main protease formed non-bonded interactions with the screened compounds. The root-mean-square difference (RMSD), root-mean-square fluctuations (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA), and hydrogen bond data from a molecular dynamics simulation study confirmed the docked complexes′ binding rigidity in the atomistic simulated environment. However, this study′s findings require in vitro and in vivo validation to ensure the possible inhibitory effects and pharmacological efficacy of the identified compounds.


The mutants of human lysozyme are capable of fibril formation implicated in the etiology of familial systemic or renal amyloidosis. A series of 100 ns thermal unfolding molecular dynamics (MD) simulations with WT human lysozyme and its seven amyloidogenic variants (I56T, D67H, F57I, W64R, Y54N, F57I/T70N and T70N/W112R) have been performed at 500 K. The molecular dynamics simulations were performed with GROMACS software (version 5.1) using the CHARMM36m force field. The MD results have been analysed in terms of the parameters characterizing both the global and local protein structure, such as the backbone root mean-square deviation, gyration radius, solvent accessible surface area, the root mean-square fluctuations and the secondary structure content. Depending on the observed effects, the examined variants of human lysozyme have been roughly divided into three groups comprising of mutants with faster (Y54N and F57I/T70N), similar (D67H and I56T) or slower (W64, F57I and T70N/W112R) unfolding rate compared to the wild-type counterpart. The analysis of the protein fluctuational behavior revealed that in most mutants the β-domain displays stronger fluctuations (except the W64R and F57I) and higher flexibility of the C- and D-helices relative to the native lysozyme with the exception of W64R and Y54N which show marked decrease (W64R) or increase (Y54N) in mobility of almost all residues. The analysis of secondary structure evolution provided evidence for higher stability of α-domain compared to β-domain. The results obtained reinforce the idea that mutation-induced global structural destabilization is not the only factor contributing to protein misfolding, the modifications in conformation and dynamics of selected protein regions may also play significant role in amyloid fibril formation.


2021 ◽  
Vol 12 (6) ◽  
pp. 7239-7248

The novel coronavirus, recognized as COVID-19, is the cause of an infection outbreak in December 2019. The effect of temperature and pH changes on the main protease of SARS-CoV-2 were investigated using all-atom molecular dynamics simulation. The obtained results from the root mean square deviation (RMSD) and root mean square fluctuations (RMSF) analyses showed that at a constant temperature of 25℃ and pH=5, the conformational change of the main protease is more significant than that of pH=6 and 7. Also, by increasing temperature from 25℃ to 55℃ at constant pH=7, a remarkable change in protein structure was observed. The radial probability of water molecules around the main protease was decreased by increasing temperature and decreasing pH. The weakening of the binding energy between the main protease and water molecules due to the increasing temperature and decreasing pH has reduced the number of hydrogen bonds between the main protease and water molecules. Finding conditions that alter the conformation of the main protease could be fundamental because this change could affect the virus’s functionality and its ability to impose illness.


Molecules ◽  
2020 ◽  
Vol 25 (20) ◽  
pp. 4829
Author(s):  
Sajjad Haider ◽  
Assem Barakat ◽  
Zaheer Ul-Haq

CXCL12 are small pro-inflammatory chemo-attractant cytokines that bind to a specific receptor CXCR4 with a role in angiogenesis, tumor progression, metastasis, and cell survival. Globally, cancer metastasis is a major cause of morbidity and mortality. In this study, we targeted CXCL12 rather than the chemokine receptor (CXCR4) because most of the drugs failed in clinical trials due to unmanageable toxicities. Until now, no FDA approved medication has been available against CXCL12. Therefore, we aimed to find new inhibitors for CXCL12 through virtual screening followed by molecular dynamics simulation. For virtual screening, active compounds against CXCL12 were taken as potent inhibitors and utilized in the generation of a pharmacophore model, followed by validation against different datasets. Ligand based virtual screening was performed on the ChEMBL and in-house databases, which resulted in successive elimination through the steps of pharmacophore-based and score-based screenings, and finally, sixteen compounds of various interactions with significant crucial amino acid residues were selected as virtual hits. Furthermore, the binding mode of these compounds were refined through molecular dynamic simulations. Moreover, the stability of protein complexes, Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), and radius of gyration were analyzed, which led to the identification of three potent inhibitors of CXCL12 that may be pursued in the drug discovery process against cancer metastasis.


2013 ◽  
Vol 2013 ◽  
pp. 1-12 ◽  
Author(s):  
Zhiguo Chen ◽  
Yi Fu ◽  
Wenbo Xu ◽  
Ming Li

Bacillus amyloliquefaciensribonuclease Barnase (RNase Ba) is a 12 kD (kilodalton) small extracellular ribonuclease. It has broad application prospects in agriculture, clinical medicine, pharmaceutical, and so forth. In this work, the thermal stability of Barnase has been studied using molecular dynamics simulation at different temperatures. The present study focuses on the contribution of noncovalent intramolecular interaction to protein stability and how they affect the thermal stability of the enzyme. Profiles of root mean square deviation and root mean square fluctuation identify thermostable and thermosensitive regions of Barnase. Analyses of trajectories in terms of secondary structure content, intramolecular hydrogen bonds and salt bridge interactions indicate distinct differences in different temperature simulations. In the simulations, Four three-member salt bridge networks (Asp8-Arg110-Asp12, Arg83-Asp75-Arg87, Lys66-Asp93-Arg69, and Asp54-Lys27-Glu73) have been identified as critical salt bridges for thermostability which are maintained stably at higher temperature enhancing stability of three hydrophobic cores. The study may help enlighten our knowledge of protein structural properties, noncovalent interactions which can stabilize secondary peptide structures or promote folding, and also help understand their actions better. Such an understanding is required for designing efficient enzymes with characteristics for particular applications at desired working temperatures.


Biology ◽  
2021 ◽  
Vol 10 (8) ◽  
pp. 789
Author(s):  
Mycal Dutta ◽  
Abu Montakim Tareq ◽  
Ahmed Rakib ◽  
Shafi Mahmud ◽  
Saad Ahmed Sami ◽  
...  

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a contemporary coronavirus, has impacted global economic activity and has a high transmission rate. As a result of the virus’s severe medical effects, developing effective vaccinations is vital. Plant-derived metabolites have been discovered as potential SARS-CoV-2 inhibitors. The SARS-CoV-2 main protease (Mpro) is a target for therapeutic research because of its highly conserved protein sequence. Gas chromatography–mass spectrometry (GC-MS) and molecular docking were used to screen 34 compounds identified from Leucas zeylanica for potential inhibitory activity against the SARS-CoV-2 Mpro. In addition, prime molecular mechanics–generalized Born surface area (MM-GBSA) was used to screen the compound dataset using a molecular dynamics simulation. From molecular docking analysis, 26 compounds were capable of interaction with the SARS-CoV-2 Mpro, while three compounds, namely 11-oxa-dispiro[4.0.4.1]undecan-1-ol (−5.755 kcal/mol), azetidin-2-one 3,3-dimethyl-4-(1-aminoethyl) (−5.39 kcal/mol), and lorazepam, 2TMS derivative (−5.246 kcal/mol), exhibited the highest docking scores. These three ligands were assessed by MM-GBSA, which revealed that they bind with the necessary Mpro amino acids in the catalytic groove to cause protein inhibition, including Ser144, Cys145, and His41. The molecular dynamics simulation confirmed the complex rigidity and stability of the docked ligand–Mpro complexes based on the analysis of mean radical variations, root-mean-square fluctuations, solvent-accessible surface area, radius of gyration, and hydrogen bond formation. The study of the postmolecular dynamics confirmation also confirmed that lorazepam, 11-oxa-dispiro[4.0.4.1]undecan-1-ol, and azetidin-2-one-3, 3-dimethyl-4-(1-aminoethyl) interact with similar Mpro binding pockets. The results of our computerized drug design approach may assist in the fight against SARS-CoV-2.


Sign in / Sign up

Export Citation Format

Share Document