The main protease of SARS-CoV-2 is one of the important targets to design and develop antiviral drugs. In this study, we have selected 40 antiviral phytochemicals to find out the best candidates ...which can act as potent inhibitors against the main protease. Molecular docking is performed using AutoDock Vina and GOLD suite to determine the binding affinities and interactions between the phytochemicals and the main protease. The selected candidates strongly interact with the key Cys145 and His41 residues. To validate the docking interactions, 100 ns molecular dynamics (MD) simulations on the five top-ranked inhibitors including hypericin, cyanidin 3-glucoside, baicalin, glabridin, and α-ketoamide-11r are performed. Principal component analysis (PCA) on the MD simulation discloses that baicalin, cyanidin 3-glucoside, and α-ketoamide-11r have structural similarity with the apo-form of the main protease. These findings are also strongly supported by root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), and solvent accessible surface area (SASA) investigations. PCA is also used to find out the quantitative structure-activity relationship (QSAR) for pattern recognition of the best ligands. Multiple linear regression (MLR) of QSAR reveals the R
2
value of 0.842 for the training set and 0.753 for the test set. Our proposed MLR model can predict the favorable binding energy compared with the binding energy detected from molecular docking. ADMET analysis demonstrates that these candidates appear to be safer inhibitors. Our comprehensive computational and statistical analysis show that these selected phytochemicals can be used as potential inhibitors against the SARS-CoV-2.
Communicated by Ramaswamy H. Sarma
Objective
Colorectal cancer (CRC) is the third leading cause of cancer-related deaths in both men and women. Toll-like receptor 5 (TLR5), an autoimmune signaling receptor that plays a role in cancer, ...can be exploited for the suppression of human colon cancer. Salmonella flagellin protein, a novel agonist of TLR5 activating downstream signaling, could be a basis for designing anticancer peptides.
Methods
The three-dimensional crystal structure of TLR5 (PDB ID: 3J0A, Resolution = 26.0 Å) was optimized using the AMBER force field in the YASARA suit. In silico enzymatic digestion tool, PeptideCutter, was used to identify peptides from Salmonella flagellin, an agonist against human TLR5. The 3D structure of the peptides was generated using PEP-FOLD3. These peptides were screened against human TLR5 using shape complementarity principles based on the binding affinity and interactions with the active residue of TLR5 monomer, and the selected peptides were further validated by molecular dynamic (MD) simulation.
Results
In this study, we generated 42 peptides from Salmonella flagellin protein by in silico protein digestion. Then, based on a new hidden Markov model sub-optimal conformation sampling approach as well as the size of the fragments, we select 38 effective peptides from these 42 cleavages. These peptides were screened against the monomeric Xray structure of human TLR5 using shape complementarity principles. Based on the binding affinity and interactions with the active residue of TLR5 monomer (residues 294 and 366 of TLR5), nine top-scored peptides were selected for the initial molecular dynamic (MD) simulation. Among these peptides, Clv10, Clv17, and Clv28 showed high stability and less flexibility during MD simulation. A 1 μs MD simulation was performed on TLR5-Clv10, TLR-Clv17, and TLR5-Clv28 complexes to further analyze the stability, conformational changes, and binding mode (Clv10, Clv17, and Clv28). During this MD study, the peptides showed high salt bridges and ionic interactions with residue ASP294 and residue ASP366 throughout the simulation and remained in the concave of the human TLR5 monomer. The RMSD and Rg values showed that the peptide-protein complexes become stable after 200 ns of contraction and extraction.
Conclusion
These findings can facilitate the rational design of selected peptides as an agonist of TLR5, which have antitumor activity, suppress colorectal cancer tumors, and can be used as promising candidates and novel agonists of TLR5.
Graphical Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmittable and pathogenic human coronavirus that caused a pandemic situation of acute respiratory syndrome, called ...COVID-19, which has posed a significant threat to global health security. The aim of the present study is to computationally design an effective peptide-based multi-epitope vaccine (MEV) against SARS-CoV-2. The overall model quality of the vaccine candidate, immunogenicity, allergenicity, and physiochemical analysis have been conducted and validated. Molecular dynamics studies confirmed the stability of the candidate vaccine. The docked complexes during the simulation revealed a strong and stable binding interactions of MEV with human and mice toll-like receptors (TLR), TLR3 and TLR4. Finally, candidate vaccine codons have been optimized for their in silico cloning in E. coli expression system, to confirm increased expression. The proposed MEV can be a potential candidate against SARS-CoV-2, but experimental validation is needed to ensure its safety and immunogenicity status.
Display omitted
•Potential T and B cell epitopes were screened from structural, non-structural and accessory proteins through immunoinformatics approach.•An immunogenic multi-subunit vaccine was hypothesized by combining the best immunodominant peptides with TLR agonist and universal PADRE sequence.•Molecular Dynamics revealed the stability of vaccine with the Toll-like-Receptors.•Epitopes from structural proteins were highly conserved with SARS CoV-1 epitopes.
Dengue, a mosquito-borne disease, has appeared as a major infectious disease globally. The virus requires its proteins to replicate and reproduce in the host cell. The NS3 protease converts the ...polyprotein to functional proteins with the help of the NS2B cofactor. Thus, NS3 protease is a promising target to develop antiviral inhibitors against the dengue virus. A systematic screening including ADMET properties, molecular docking, molecular dynamics (MD) simulation, binding free energy calculation, and QSAR studies is carried out to predict potent inhibitors against the NS3 protease. From the screening of 40 antiviral phytochemicals, ADMET properties analysis was used to screen out ligands that violate ADME rules and have probable toxicity. Cyanidin 3-Glucoside, Dithymoquinone, and Glabridin were predicted to be potent inhibitors against the NS3 protease according to their binding affinity. These ligands showed several noncovalent interactions, including hydrogen bond, hydrophobic interaction, electrostatic interaction, pi-sulfur interactions. The ligand-protein complexes were further scrutinized using 250 ns molecular dynamics simulation. The MM-PBSA binding free energy calculation was conducted to investigate their binding stability in dynamic conditions. The calculated pIC50(mM) value was predicted using the QSAR model with 89.91% goodness of fit. The predicted biologocal activity value for the ligands indicates they might have good potency.
Display omitted
•40 Antiviral phytochemicals screened against DENV-2 NS3-NS2B protease.•Initially, ADMET and Molecular docking were performed for screening.•250ns MD and MM/PBSA were implemented for the three selected candidates.•Phytochemicals Glabridin altered the protein-cofactor protein interactions interface.•QSAR model were studied for in-silico activity prediction.
Saprolegnia parasitica is an oomycete responsible for a fish disease called saprolegniosis, which poses an economic and environmental burden on aquaculture production. In Saprolegnia, CHS5 of S. ...parasitica (SpCHS5) contains an N-terminal domain, a catalytic domain of the glycosyltransferase −2 family containing a GT-A fold, and a C-terminal transmembrane domain. No three-dimensional structure of SpCHS5 is reported yet disclosing the structural details of this protein. We have developed a structural model of full-length SpCHS5 and validated it by molecular dynamics simulation technique. From the 1 microsecond simulations, we retrieved the stable RoseTTAFold model SpCHS5 protein to explain characteristics and structural features. Furthermore, from the analysis of the movement of chitin in the protein cavity, we assumed that ARG 482, GLN 527, PHE 529, PHE 530, LEU 540, SER 541, TYR 544, ASN 634, THR 641, TYR 645, THR 641, ASN 772 residues as a main cavity lining site. In SMD analysis, we investigated the opening of the transmembrane cavity required for chitin translocation. The pulling of chitin from the internal cavity to the extracellular region was observed through steered molecular dynamics simulations. A comparison of the initial and final structures of chitin complex showed that there's a transmembrane cavity opening in the simulations. Overall, this present work will help us understand the structural and functional basis of CHS5 and design inhibitors against SpCHS5.
Communicated by Ramaswamy H. Sarma
Conformational changes of chitin at the protein tunnel where 'a' is the orientation of chitin at 1000 ns, and 'b' is the initial orientation at 0 ns of the molecular dynamic simulation.
The COVID‐19 pandemic has been a public health emergency, with deadly forms constantly emerging around the world, highlighting the dire need for highly effective antiviral therapeutics. Peptide ...therapeutics show significant potential for this viral disease due to their efficiency, safety, and specificity. Here, two thousand seven hundred eight antibacterial peptides were screened computationally targeting the Main protease (Mpro) of SARS CoV‐2. Six top‐ranked peptides according to their binding scores, binding pose were investigated by molecular dynamics to explore the interaction and binding behavior of peptide‐Mpro complexes. The structural and energetic characteristics of Mpro‐DRAMP01760 and Mpro‐DRAMP01808 complexes fluctuated less during a 250 ns MD simulation. In addition, three peptides (DRAMP01760, DRAMP01808, and DRAMP01342) bind strongly to Mpro protein, according to the free energy landscape and principal component analysis. Peptide helicity and secondary structure analysis are in agreement with our findings. Interaction analysis of protein‐peptide complexes demonstrated that Mpro's residue CYS145, HIS41, PRO168, GLU166, GLN189, ASN142, MET49, and THR26 play significant contributions in peptide‐protein attachment. Binding free energy analysis (MM‐PBSA) demonstrated the energy profile of interacting residues of Mpro in peptide‐Mpro complexes. To summarize, the peptides DRAMP01808 and DRAMP01760 may be highly Mpro specific, resulting disruption in a viral replication and transcription. The results of this research are expected to assist future research toward the development of antiviral peptide‐based therapeutics for Covid‐19 treatment.
When SARS‐CoV‐2 enters inside the cell by fusion process it encodes few proteins including main protease with the help of cellular ribosomes. This main protease mediates the maturation cleavage of polyproteins at 11 different regions, which is a notable step in virus's life cycle. SARS‐CoV‐2 Mpro is mainly a homo‐dimeric enzyme consisting of three domains. The large substrate‐binding site with a catalytic dyad of His41 and Cys145 is located between chymotrypsin‐like domains I and picornavirus 3C protease‐like domain II. Inhibition of this substrate binding site disrupt the replication process. But being large substrate‐binding site, peptide is more suitable than small molecules for inhibiting this crucial site. The outcome of this computational study displays that the attachment of peptide with Mpro protease can be a good strategy to inhibit the substrate‐binding site which may lead to disruption of replication process of SARS CoV‐2.
The ORF3a is a large accessory protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which plays an important role in virulence and viral replication; especially in inflammasome ...activation and apoptosis. However,, the existing cryo-EM structure of SARS-CoV-2 ORF3a is incomplete, . making it challenging to understand its structural and functional features. The aim of this study is to investigate the dynamic behaviors of the full-sequence homology model of ORF3a and compare it with the cryo-EM structure using microsecond molecular dynamics simulations. The previous studies indicated that the unresolved residues of the cryo-EM structure are not only involved in the pathogenesis of the SARS-CoV-2 but also exhibit a significant antigenicity. The dynamics scenario of homology model revealed higher RMSD, Rg, and SASA values with stable pattern when compared to the cryo-EM structure. Moreover, the RMSF analysis demonstrated higher fluctuations at specific positions (1-43, 97-110, 172-180, 219-243) in the model structure, whereas the cryo-EM structure displayed lower overall drift (except 1-43) in comparison to the model structure.Secondary structural features indicated that a significant unfolding in the transmembrane domains and β-strand at positions 166 to 172, affecting the stability and compactness of the cryo-EM structure , whereas the model exhibited noticeable unfolding in transmembrane domains and small-coiled regions in the N-terminal. , The results from molecular docking and steered molecular dynamics investigations showed the model structure had a greater number of non-bonding interactions, leading to enhanced stability when compared to the cryo-EM structure. Consequently, higher forces were necessary for unbinding of the baricitinib and ruxolitinib inhibitors from the model structure.. Our findings can help better understanding of the significance of unresolved residues at the molecular level. Additionally, this information can guide researchers for experimental endeavors aimed at completing the full-sequence structure of the ORF3a.
Communicated by Ramaswamy H. Sarma
Nucleoprotein is a conserved structural protein of SARS-CoV-2, which is involved in several functions, including replication, packaging, and transcription. In this research, 21 antiviral peptides ...that are known to have inhibitory function against nucleoprotein in several other viruses, were screened computationally against the nucleoprotein of SARS-CoV-2. The complexes of five best performing peptides (AVP1142, AVP1145, AVP1148, AVP1150, AVP1155) with nucleoprotein were selected for subsequent screening via 5 ns molecular dynamics (MD) simulation. Two peptides, namely AVP1145 and AVP1155, came out as promising candidates and hence were selected for 200 ns MD simulation for further validation, incorporating a DMPC-based membrane environment. In the long MD simulation, both AVP1155 and AVP1145 utilized multiple residues—mainly aromatic, acidic, and nonpolar residues—as interacting points to remain in contact with the nucleoprotein and formed predominantly hydrogen bonds along with hydrophobic and electrostatic interactions. However, AVP1155 proved to be superior to AVP1145 when its complex with nucleoprotein was analyzed in terms of root-mean-square deviation, root-mean-square fluctuation, radius of gyration, solvent accessible surface area and free energy landscape. In a nutshell, the findings of this research may guide future studies in the development of selective peptide inhibitors of SARS-CoV-2 nucleoprotein.
Mixed surfactant systems are used in different applied fields like pharmaceutical formulation rather than single surfactant. Therefore, the determination of the clouding nature of the triton X-100 ...(TX-100) + Tween 80 (TW-80) mixture was carried out in H
O and polyvinyl alcohol (PVA). In the occurrence of PVA, the cloud point (CP) values of TX-100 initially enhance with enhancing the concentration of PVA and tend to decrease after a certain concentration. For different ratios of TX-100 and TW-80 mixture having the same concentration of both solutions, CP values increase through the decreasing ratios of TX-100 with/without PVA. In the presence of polymer, at higher ratios of TX-100 than TW-80, the CP values are higher in magnitudes in comparison to the aqueous medium but at lower ratios of TX-100, the value of CP are lower in magnitudes in comparison to the aqueous system. The CP values of the TX-100 + TW-80 mixture in the salt system are lower in magnitudes than the aqueous medium in both the absence/presence of PVA. However, a reduction of CP values was obtained to a large extent for Na
SO
over NaCl in the case of lower volume ratios of TX-100. Various thermodynamic variables (standard free energy (
), standard enthalpy (
), standard entropy (
) change, thermodynamic parameters of transfer (free energy of transfer (
), and transfer of enthalpies (
)) of phase transition) were also determined.
The multidrug transporter P-glycoprotein is an ATP binding cassette (ABC) exporter responsible for resistance to tumor cells during chemotherapy. This study was designed with computational approaches ...aimed at identifying the best potent inhibitors of P-glycoprotein. Although many compounds have been suggested to inhibit P-glycoprotein, however, their information on bioavailability, selectivity, ADMET properties, and molecular interactions has not been revealed. Molecular docking, ADMET analysis, molecular dynamics, Principal component analysis (PCA), and binding free energy calculations were performed. Two compounds D1 and D2 showed the best docking score against P-glycoprotein and both compounds have 4-thiazolidinone derivatives containing indolin-3 one moiety are novel anti-tumor compounds. ADMET calculation analysis predicted D1 and D2 to have acceptable pharmacokinetic properties. The MD simulation discloses that D1-P-glycoprotein and D2-P-glycoprotein complexes are in stable conformation as apo-form. Hydrophobic amino acid such as phenylalanine plays significant on the interactions of inhibitors. Principal component analysis shows that both complexes are relatively similar variables as apo-form except planarity and Columbo energy profile. In addition, Quantitative Structural Activity Relationship (QSAR) of the ligand candidates were subjected to the principal component analysis (PCA) for pattern recognition. Partial-least-square (PLS) regression analysis was further utilized to model drug candidates' QSAR for subsequent prediction of the binding energy of validated drug candidates. PCA revealed groupings of the drug candidates based on the similarity or differences in drug candidates QSAR. Moreover, the developed PLS regression accurately predicted the values of the binding energy of drug candidates, with low residual error of prediction.
Communicated by Ramaswamy H. Sarma