Ns Molecular Dynamics Simulations Research Articles

Reverse Mapping of Coarse Grained Polyglutamine Conformations from PRIME20 Sampling.

An inverse coarse-graining protocol is presented for generating and validating atomistic structures of large (bio-) molecules from conformations obtained via a coarse-grained sampling method. Specifically, the protocol is implemented and tested based on the (coarse-grained) PRIME20 protein model (P20/SAMC), and the resulting all-atom conformations are simulated using conventional biomolecular force fields. The phase space sampling at the coarse-grained level is performed with a stochastical approximation Monte Carlo approach. The method is applied to a series of polypeptides, specifically dimers of polyglutamine with varying chain length in aqueous solution. The majority (>70 %) of the conformations obtained from the coarse-grained peptide model can successfully be mapped back to atomistic structures that remain conformationally stable during 10 ns of molecular dynamics simulations. This work can be seen as the first step towards the overarching goal of improving our understanding of protein aggregation phenomena through simulation methods.

Computationally Predicted New Solid-State Electrolyte (Li5+x PS4+x Cl2-x : 0 ≤ x ≤ 2) and Poly Sulfide Cathodes (Li3+y PS9 or Li5+y PS9Cl2: 0 ≤ y ≤ 9) for High Performance Li Metal Anode Batteries

We used Molecular Dynamics (MD) simulations to study the structures and conductivity of sulfide-based ceramics from the superionic conductor argyrodite family for applications as solid-state electrolyte and cathode materials. These studies combined quantum mechanics (QM) and ReaxFF reactive force field based MD, where the ReaxFF parameters were developed from QM MD.[1]Currently Li10GeP2S12 (LGPS) demonstrates one of the highest Li-ion conductivity for solid state electrolyte, 12 mS/cm at 300K, but it is highly reactive with the Li-metal anode, and it is expensive due to the presence of Ge. We studied the Li ion diffusion mechanism and ionic conductivity of the promising Li6(PS4)SCl solid state electrolyte, predicted to have a conductivity of s = 6 mS/cm, close to experiment (4 mS/cm). We find that Li migration in this electrolyte occurs via conjugated substitutional type diffusion involving rearrangements of three (or more) Li-ions and ~ 3 vacant sites in a 3D matrix of anions that are essentially stationary at 298K (over 20 ns of simulation). We carried out 10 ns of MD simulation to predict a Li-ion conductivity for a single phase Li6(PS4)SCl of 5.9 mS/cm in agreement with solid state NMR measurements of 3.9 mS/cm. Our calculated activation energy of 0.24 eV is within experimentally reported range.[2]We report here the computationally predicted structure of Li5PS4Cl2, a new sulfide Li-superionic conductor with the highest Li-ion conductivity at solid state, which we predict to have Li-ion conductivity of 20 mS/cm at 298 K. We also predicted the directional ionic conductivity and Li-migration barrier in it. This Li5PS4Cl2 has not yet been synthesized and studied experimentally.We also report our results on QM and ReaxFF MD for new polysulfide cathodes, Li3+y PS4+n and Li5+y PS4+nCl2 , both based on Li3PS4. Here, n is the number of excess sulfur atoms in the fully charged polysulfide cathode per formula unit with y=0 and y is the excess Li added to the cathode during discharge processes. We evaluated the performance of these cathode electrolyte systems in terms of interfacial stability and discharge voltages.Li3PS4+5 is a potential cathode for lithium-sulfur batteries. We predicted the structures of Li3PS4+5 finding extra S3 and S7 chains attached to one S atom of each PS4 group. This leads to a density of 2.2 g/cm3. As the cathode is discharged, the lithium atoms fill up all the gaps between S atoms, leading to Li3+9 PS4+5, with 1.7 g/cm3. We studied the dynamics of the discharge as lithium ions from the metal move into the material and react with the S-S bonds to make Li2S. Our predicted discharge curve agrees well with the experimental data (Figure 1).[3]We extended these studies to design a new polysulfide cathode Li5PS4+5Cl2 to be used with our newly predicted Li-superionic conductor electrolytes. The interfacial stability of the Solid electrolyte with S-based cathode and with Li-anode were studied (Figure 2).Our work presents a new strategy for designing materials and processes for solid-state batteries using computational screening and chemical substitution. This demonstrates the potential of using lithium PS4 based superionic conductors with PS4 based cathodes for applications in solid-state batteries. We also addresses some of the key challenges in accelerating solid-state battery development, such as phase stability, transport properties, electrode compatibility etc. This work demonstrates how computational studies can guide new avenues for experimental developments to advance of solid-state battery technology for applications ranging from grid storage to electric vehicles.Acknowledgements:We acknowledge the support from Hong Kong Quantum AI Lab, AIR@InnoHK of Hong Kong Government.

Computational Risk Assessment of Persistence, Bioaccumulation, and Toxicity of Novel Flame-Retardant Chemicals.

Novel brominated flame retardants (NBFRs) have emerged as chemicals of environmental concern, as they have been widely used as an alternative to polybrominated diphenyl ethers (PBDEs). Considering the similar structural features of NBFRs and PBDEs necessitates a comprehensive investigation to understand the physicochemical relationships of these compounds and their ability to alter biological functions. In this study, we investigated the persistent nature of NBFRs in terms of thyroid-disrupting potential by understanding the structure-stability aspects using density functional theory (DFT)-based reactivity parameters and interactions via molecular docking and molecular dynamics (MD) simulations. The results indicate that the DFT-based stability descriptor (chemical hardness) is associated with the persistent nature of NBFRs. The computed molecular interaction profile revealed prominent interactions between thyroid receptor-β (TR-β) and NBFRs. Stable trajectory and interactions with TR-β were obtained with ATE, p-TBX, PBT, PBEB, and TBBPA-DBPE during 100 ns of MD simulation. The results of these studies have suggested that the presence of a higher number of halogenated atoms increases the stability vis-à-vis the persistence and endocrine disruption potential of NBFRs.

Computational assessment and in vitro test of phytochemicals of Usnea aciculifera as potential inhibitors of Escherichia coli efflux pump AcrB

Lichens produce secondary metabolites that have many pharmaceutical activities such as antimicrobial, antioxidant, antiviral, anticancer, antigenotoxic, anti-inflammatory, analgesic and antipyretic activities. However, there is limited research on their efflux pump inhibitory activities. Twelve phytochemicals were isolated from Usnea aciculifera, and their activity of AcrAB-TolC efflux pump inhibition was evaluated. Four potential compounds, which are diffractaic acid (2), 8′ -O- methylstictic acid (5), 3-hydroxy-4-(methoxycarbonyl)-2,5-dimethylphenyl 2,4-dimethoxy-3,6-dimethylbenzoate (8) and 3-hydroxy-4-(methoxycarbonyl)-2,5-dimethylphenyl 2-hydroxy-4-methoxy-3,6-dimethylbenzoate (9), were found by virtual screening using pharmacophore and 2D-QSAR model. Compound 8 exhibited AcrB inhibition activity in vitro with an accumulation H33342 percentage compared with untreated control of 202% at a concentration of 50 µM and increased the antibacterial activity of levofloxacin by four-fold at a concentration of 200 µM. By molecular docking and molecular dynamics (MD) simulation, the binding affinity of depside and depsidone derivatives to AcrB was also clarified. Despite the poor docking score to the AcrB binding site, compound 8 was the most stable among the four complexes at 20 ns of MD simulation. The analysis of long MD at 100 ns indicated that compound 8 interacts strongly with the residues in the distal pocket, creating a stable complex with ΔGbind of −31.51 kcal.mol-1. According to the ADMETlab 2.0 web server’s predictions of pharmacokinetics and toxicities, compound 8 has the potential for drug development. Communicated by Ramaswamy H. Sarma

Identification of small-molecule glucokinase activator for type-2-diabetes treatment: a structure-based virtual screening approach

ABSTRACT Glucokinase (GK, EC 2.7.1.2) is a crucial enzyme that catalyses the conversion of glucose to glucose-6-phosphate. It is used to treat type-2 diabetes (T2D), a serious metabolic disorder that is still at the forefront without proper medication. Fast Rigid Exhaustive Docking (FRED) was carried out for 400,000 compounds from the Zinc database to identify novel glucokinase activators. The hit compounds ZINC69775727, ZINC9114647, ZINC91773667, ZINC9305321, and ZINC96165848 interacted strongly with allosteric site residues and, formed hydrogen bonds with ARG 63. The hit compounds met the criterion for drug-likeness, according to the ADME prediction. The compounds were then subjected to 100 ns of molecular dynamics simulation and MM-GBSA calculation using DESMOND. The findings demonstrated that the compounds had good stability and minimal fluctuation throughout the course of the simulation, pointing to the potential of the chosen compounds for glucokinase activation. The compound ZINC69775727 in particular has the lowest binding energy of −111.1 kcal/mol, which is lower than the native ligand’s binding energy of −102.84 kcal/mol and the binding energies of the control compounds PSN-GK1 and Piragliatin, which are −102.49 kcal/mol and −107.767 kcal/mol, respectively. Therefore, the information from this work may be useful in finding novel small molecules as GKAs.

Investigation of bio-active Amaryllidaceae alkaloidal small molecules as putative SARS-CoV-2 main protease and host TMPRSS2 inhibitors: interpretation by in-silico simulation study

The novel coronavirus disease 2019 (Covid-19) outburst is still threatening global health. This highly contagious viral disease is caused by the infection of SARS-CoV-2 virus. Covid-19 and post-Covid-19 complications induce noteworthy mortality. Potential chemical hits and leads against SARS-CoV-2 for combating Covid-19 are urgently required. In the present study, a virtual-screening protocol was executed on potential Amaryllidaceae alkaloids from a pool of natural compound library against SARS-CoV-2 main protease (Mpro) and transmembrane serine protease (TMPRSS2). For the collected 1016 alkaloids from the curated library, initially, molecular docking using AutoDock Vina (ADV), and thereafter 100 ns molecular-dynamic (MD) simulation has been executed for the best top-ranked binding affinity compounds for both the viral and host proteins. Comprehensive intermolecular-binding interactions profile of Amaryllidaceae alkaloids suggested that phyto-compounds Galantamine, Lycorenine, and Neronine as potent modulators of SARS-CoV-2 Mpro and host TMPRSS2 protein. All atomistic long range 100 ns MD simulation studies of each top ranked complex in triplicates also illustrated strong binding affinity of three compounds towards Mpro and TMPRSS2. Identified compounds might be recommended as prospective anti-viral agents for future drug development selectively targeting the SARS-CoV-2 Mpro or blocking host TMPRSS2 receptor, subjected to pre-clinical and clinical assessment for a better understanding of in-vitro molecular interaction and in-vivo validation. Communicated by Ramaswamy H. Sarma

Novel pseudonucleosides and sulfamoyl-oxazolidinone β-D-glucosamine derivative as anti-COVID-19: design, synthesis, and in silico study

New pseudonucleosides containing cyclic sulfamide moiety and sulfamoyl β-D-glucosamine derivative are described. These pseudonucleosides are synthesized in good yields starting from chlorosulfonyl isocyanate and β-D-glucosamine hydrochloride in five steps; (protection, acetylation, removal of the Boc group, sulfamoylation, and cyclization). Further, novel glycosylated sulfamoyloxazolidin-2-one is prepared in three steps; carbamoylation, sulfamoylation, and intramolecular cyclization. The structures of the synthesized compounds were confirmed by usual spectroscopic and spectrometric methods NMR, IR, MS, and EA. Interesting molecular docking of the prepared pseudonucleosides and (Beclabuvir, Remdesivir) drugs with SARS-CoV-2/Mpro (PDB:5R80) was conducted using the same parameters for a fair comparison. A low binding affinity of the synthesized compounds compared to the Beclabuvir and other analysis showed that pseudonucleosides have the ability to inhibit SARS-CoV-2. After the motivating results of molecular docking study, the complex between the SARS-CoV-2 Mpro and compound 7 was subjected to 100 ns molecular dynamics (MD) simulation using Desmond module of Schrodinger suite, during which the receptor-ligand complex showed substantial stability after 10 ns of MD simulation. Also, we studied the prediction of absorption, distribution, properties of metabolism, excretion, and toxicity (ADMET) of the synthesized compounds. Communicated by Ramaswamy H. Sarma

Synthesis, in silico study, and biological evaluation of cyclosulfamide derivatives as new anticholinesterase inhibitors

After the great spread of Alzheimer's disease (AD) in the past few years, many researchers directed their work toward developing an effective treatment for this disease and to discovering very advanced drugs. Cyclosulfamides are considered as versatile pharmacophores in the construction of new molecules with excellent activities. Therefore, a series of cyclosulfamide have been synthesized and evaluated as anti-Alzheimer agents through in vitro acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) inhibition. Most of tested compounds showed an average inhibitory activity against two targeted enzymes compared with the reference ligand. Specifically, 10a showed the best inhibition of AChE enzyme with an IC50 = 45.30±0.88 µM; while 4a exhibited the most potent BuChE enzyme inhibition with an IC50 = 52.87±3.73 µM.Further, we were able to determine a feasible binding mode for cyclosulfamide derivatives owing to molecular docking studies, which offered prospective evidence to identify significant interactions between the active site of AChE and the synthesized ligands.Following the encouraging findings of the molecular docking investigation, the complex AChE-10a was put through a 100 ns Desmond of Schrodinger simulation of molecular dynamics (MD), during which the receptor-ligand combination showed significant stability after 10 ns of MD simulation.

Potential SARS-CoV-2 RdRp inhibitors of cytidine derivatives: Molecular docking, molecular dynamic simulations, ADMET, and POM analyses for the identification of pharmacophore sites.

The RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 is one of the optimum targets for antiviral drug design and development. The hydroxyl groups of cytidine structures were modified with different aliphatic and aromatic groups to obtain 5´-O-acyl and 2´,3´-di-O-acyl derivatives, and then, these derivatives were employed in molecular modeling, antiviral prediction, molecular docking, molecular dynamics, pharmacological and POM studies. Density functional theory (DFT) at the B3LYP/6-31G++ level analyzed biochemical behavior and molecular electrostatic potential (MESP) of the modified cytidine derivatives. The antiviral parameters of the mutated derivatives revealed promising drug properties compared with those of standard antiviral drugs. Molecular docking has determined binding affinities and interactions between the cytidine derivatives and SARS-CoV-2 RdRp. The modified derivatives strongly interacted with prime Pro620 and Lys621 residues. The binding conformation and interactions stability were investigated by 200 ns of molecular dynamics simulations and predicted the compounds to firmly dock inside the RdRp binding pocket. Interestingly, the binding residues of the derivatives were revealed in high equilibrium showing an enhanced binding affinity for the molecules. Intermolecular interactions are dominated by both Van der Waals and electrostatic energies. Finally, the pharmacokinetic characterization of the optimized inhibitors confirmed the safety of derivatives due to their improved kinetic properties. The selected cytidine derivatives can be suggested as potential inhibitors against SARS-CoV-2. The POM Theory supports the hypothesis above by confirming the existence of an antiviral (Oδ--O'δ-) pharmacophore site of Hits.

Modification and application of highly active alkaline pectin lyase

Alkaline pectate lyase has developmental prospects in the textile, pulp, paper, and food industries. In this study, we selected BacPelA, the pectin lyase with the highest expression activity from Bacillus clausii, modified and expressed in Escherichia coli BL21(DE3). Through fragment replacement, the catalytic activity of the enzyme was significantly improved. The optimum pH and temperature of the modified pectin lyase (PGLA-rep4) were 11.0 and 70 °C, respectively. It also exhibited a superior ability to cleave methylated pectin. The enzyme activity of PGLA-rep4, measured at 235 nm with 0.2% apple pectin as the substrate, was 554.0 U/mL, and the specific enzyme activity after purification using a nickel column was 822.9 U/mg. After approximately 20 ns of molecular dynamics simulation, the structure of the pectin lyase PGLA-rep4 tended to be stable. The root mean square fluctuation (RMSF) values at the key catalytically active site, LYS168, were higher than those of the wildtype PGLA. In addition, PGLA-rep4 was relatively stable in the presence of metal ions. PGLA-rep4 has good enzymatic properties and activities and maintains a high pH and temperature. This study provides a successful strategy for enhancing the catalytic activity of PGLA-rep4, making it the ultimate candidate for degumming and various uses in the pulp, paper, and textile industries.

Backbone N-Amination Promotes the Folding of β-Hairpin Peptides via a Network of Hydrogen Bonds.

Molecular dynamics (MD) simulations have been used to characterize the effects of backbone N-amination of residues in a model β-hairpin peptide. This modification is of considerable interest as N-aminated peptides have been shown to inhibit amyloid-type aggregation. Six derivatives of the β-hairpin peptide, which contain one, two, or four N-aminated residues, have been studied. For each peptide 100 ns MD simulations starting from the folded β-hairpin structure were performed. The effects of the N-amination prove to be very sequence dependent. N-Amination of a residue involved in interstrand hydrogen bonding (Val3) leads to unfolding of the β-hairpin, whereas N-amination of a residue toward the C-terminus (Leu11) gives fraying at the termini of the peptide. In the other derivatives the peptide remains folded, with increasing levels of N-amination reducing the right-handed twist of the β-hairpin and favoring population of a type II' rather than a type I' β-turn. MD simulations (100 ns) have also been run for each peptide starting from an unfolded extended chain. Here, the peptide with four N-aminated residues shows the most folding into the β-hairpin (34%). Analysis of the simulations shows that N-amination favors the population of β (φ, ψ) conformations by the preceding residue due to, at least in part, a network of weak NH2(i)-CO(i) and NH2(i)-CO(i-2) hydrogen bonds. It also leads to a reduction of misfolding because of changes in the hydrogen-bonding potential. Both of these features help funnel the peptide to the folded β-hairpin structure. The conformational insights provided through this work give a firm foundation for the design of N-aminated peptide inhibitors for modulating protein-protein interactions and aggregation.

A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds.

Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has fomented a climate of fear worldwide due to its rapidly spreading nature, and high mortality rate. The World Health Organization declared it a global pandemic on 11 March 2020. Many endeavors have been made to find appropriate medications to restrain the SARS-CoV-2 infection from spreading but there is no specific antiviral therapy to date. However, a computer-aided drug design approach can be an alternative to identify probable drug candidates within a short time. SARS-CoV-2 main protease is a proven drug target, and it plays a pivotal role in viral replication and transcription. Methods: In this study, we identified a total of 114 essential oil compounds as a feasible anti-SARS-CoV-2 agent from several online reservoirs. These compounds were screened by incorporating absorption, distribution, metabolism, excretion, and toxicity (ADMET) profiling, molecular docking, and 50 ns of molecular dynamics simulation to identify potential drug candidates. The crystallized SARS-CoV-2 main protease structure was collected from the Research Collaboratory for Structural Bioinformatics Protein Data Bank database (Protein Data Bank ID 6LU7). Results: According to the results of the ADMET study, none of the compounds have any side effects that could reduce their druglikeness or pharmacokinetic properties. Among 114 compounds, we selected bisabololoxide B, eremanthin, and leptospermone as top drug candidates based on their higher binding affinity scores, and strong interaction with the Cys 145-His 41 catalytic dyad. Finally, the molecular dynamics simulation was implemented to evaluate the structural stability of the ligand-receptor complex. Molecular dynamics simulation disclosed that all the hits showed conformational stability compared to the positive control α-ketoamide. Conclusions: Our study showed that the top three hits might work as potential anti-SARS-CoV-2 agents, which can pave the way for discovering new drugs, but further in vivo trials will require for experimental validation.

Development of novel monoamine oxidase B (MAO-B) inhibitors by combined application of docking-based alignment, 3D-QSAR, ADMET prediction, molecular dynamics simulation, and MM_GBSA binding free energy

Unsaturated ketone derivatives are known as inhibitors of monoamine oxidase B (MAO-B), a potential drug target of Parkinson’s disease. Here, docking-based alignment, 3 D-QSAR (three-dimensional quantitative structure-activity relationship) studies, ADMET (absorption, distribution, metabolism, excretion, and toxicity) prediction, molecular dynamics (MD) simulation, and MM_GBSA binding free energy were performed on a novel series of MAO-B inhibitors. The objective is to predict new MAO-B inhibitors with high potency activity. The 3 D-QSAR models were created using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). Molecular docking findings indicated that compounds with strong inhibitory efficacy also had a high binding affinity. 3 D-QSAR studies showed the importance of steric, electrostatic, and H-bond acceptor fields on the inhibitory activity of MAO-B. Based on the appropriate 3 D-QSAR model, a new series of MAO-B inhibitors were predicted and their pharmacokinetic characteristics were evaluated using in silico ADMET prediction. All screened compounds show good oral bioavailability without any side effects. Moreover, the dynamic behavior and stability of the most active compounds were evaluated using MD simulations. The results showed that unsaturated ketone derivatives are stable and compact during the 100 ns of MD simulation. Finally, the binding free energy of complexes was determined using the MM_GBSA method; the findings indicated that the T1 compound is more stable (ΔGbinding = −409.506 KJ/mol) than the data set's highest active compound (ΔGbinding = −31.883 KJ/mol). Communicated by Ramaswamy H. Sarma

Molecular modelling and in vitro studies of Daruharidra as a potent alpha-amylase inhibitor

The present study aims at exploring the potential of the Daruharidra plant (stem and bark) for inhibition of alpha-amylase. Aqueous and ethanolic extraction yielded the highest total phenolic content (TPC) of 101.4 and 111.8 mcg of gallic acid equivalent. Methanol and ethanol extract had Total flavonoid content (TFC) of 319.6 and 288.3 mcg of quercetin equivalent, respectively. In contrast, petroleum ether extraction resulted in the lowest TPC of 23.6 and TFC of 8.33 mcg, respectively. Methanol (5.554 mg/ml), acetone (6.576 mg/ml), and ethanol (7.321 mg/ml) extract had the lowest IC50 values in alpha-amylase inhibition with the mode of inhibition being non-competitive inhibition. HR-LCMS was used for comprehension of phytoconstituents present in the extract. Amongst hundreds of hits observed 10 ligands of alkaloid nature were used for docking studies. Berbamine, alloxanthine, protopine and benazepril along with reference molecule (Acarbose) were subjected to Molecular dynamics (MD) simulation to analyze the stability of the docked protein-ligand complex. The values of RMSD, RMSF, RG, H-Bond and SASA, the interaction energy of all protein-ligand complexes were calculated after 150 ns of MD simulation. The results of screened complexes revealed good stability as compared to reference Acarbose. These screened ligands used for simulation have the most negative binding energies that interacted with alpha amylase enzyme having −9.28 kcal/mol, −7.51 kcal/mol, −7.73 kcal/mol and −8.00 kcal/mol, energies respectively. The results show significant alpha-amylase inhibitory activity and interaction of ligands targeting this enzyme, which can be used for cross-validation, in vitro Communicated by Ramaswamy H. Sarma

Different compounds against Angiotensin-Converting Enzyme 2 (ACE2) receptor potentially containing the infectivity of SARS-CoV-2: an in silico study.

Novel SARS coronavirus or SARS-CoV-2 is a novel coronavirus that was identified and spread from Wuhan in 2019. On January 30th, the World Health Organization declared the coronavirus outbreak as a Global Public Health Emergency. Although Remdesivir and Molnupiravir are FDA-approved drugs for COVID-19, finding new efficient and low-cost antiviral drugs against COVID-19 for applying in more countries can still be helpful. One of the potential sources for finding new and low-cost drugs is the herbal compounds in addition to repurposing FDA-approved drugs. So, in this study, we focused on finding effective drug candidates against COVID-19 based on the computational approaches. As ACE2 serves as a critical receptor for cell entry of this virus. Inhibiting the binding site of SARS-CoV-2 on human ACE2 provides a promising therapeutic approach for developing drugs against SARS-CoV-2. Herein, we applied a bioinformatics approach to identify possible potential inhibitors of SARS-CoV-2. A library of FDA-approved compounds and five natural compounds was screened using Smina docking. Top-docking compounds are then applied in Molecular Dynamics (MD) simulation to assess the stability of ACE2-inhibitor complexes. Results indicate that Luteolin and Chrysin represent high conformation stability with ACE2 during 120 ns of Molecular Dynamics simulation. The binding free energies of Luteolin and Chrysin were calculated by the Molecular Mechanics/Poisson–Boltzmann Surface Area method (MM/PBSA) which confirmed the relative binding free energy of these drugs to ACE2 in favor of the effective binding. So, Luteolin and Chrysin could sufficiently interact with ACE2 and block the Spike binding pocket of ACE2 and can be a potential inhibitor against the binding of SARS-CoV-2 to ACE2 receptor which is an early stage of infection. Luteolin and Chrysin could be suggestive as beneficial compounds for preventing or reducing SARS-CoV-2 transmission and infection which need experimental work to prove.Graphical abstract Supplementary InformationThe online version contains supplementary material available at 10.1007/s00894-022-05059-1.

Development a novel multiepitope DNA vaccine against human SARS coronavirus-2: an immunoinformatic designing study.

Human SARS coronavirus 2 (SARS-CoV-2) causes the current global COVID-19 pandemic. The production of an efficient vaccine against COVID-19 is under heavy investigation. In this study, we have designed a novel multiepitope DNA vaccine against SARS-CoV-2 using reverse vaccinology and DNA vaccine approaches. Applying these strategies led to reduce the time and costs of vaccine development and also improve the immune protective characteristics of the vaccine. For this purpose, epitopes of nucleocapsid, membrane glycoprotein, and ORF8 proteins of SARS-CoV-2 chose as targets for B and T-cell receptors. Accordingly, DNA sequences of selected epitopes have optimized for protein expression in the eukaryotic system. To this end, the Kozak and tissue plasminogen activator sequences were added into the epitope sequences for proper protein expression and secretion, respectively. Furthermore, interleukin-2 and beta-defensin 1 preproprotein sequences were incorporated to the designed DNA vaccine as an adjuvant. Modeling and refinement of fused protein composed of SARS-CoV-2 multiepitope antigens (fuspMA) have performed based on homology modeling of orthologous peptides, then constructed 3D model of fuspMA was more investigated during 50 ns of molecular dynamics simulation. Further bioinformatics predictions demonstrated that fuspMA is a stable protein with acceptable antigenic features and no allergenicity or toxicity characteristics. Finally, the affinity of fuspMA to the MHC I and II and TLRs molecules validated by the molecular docking procedure. In conclusion, it seems the designed multiepitope DNA vaccine could have a chance to be introduced as an efficient vaccine against COVID-19 after more in vivo evaluations.

Diversity of bacteriocins in the microbiome of the Tucuruí Hydroelectric Power Plant water reservoir and three-dimensional structure prediction of a zoocin.

Bacteriocins are antimicrobial peptides expressed by bacteria through ribosomal activity. In this study, we analyzed the diversity of bacteriocin-like genes in the Tucuruí-HPP using a whole-metagenome shotgun sequencing approach. Three layers of the water column were analyzed (photic, aphotic and sediment). Detection of bacteriocin-like genes was performed with blastx using the BAGEL4 database as subject sequences. In order to calculate the abundance of bacteriocin-like genes we also determined the number of 16S rRNA genes using blastn. Taxonomic analysis was performed using RAST server and the metagenome was assembled using IDBA-UD in order to recover the full sequence of a zoocin which had its three-dimensional structure determined. The photic zone presented the highest number of reads affiliated to bacteriocins. The most abundant bacteriocins were sonorensin, Klebicin D , pyocin and colicin. The zoocin model was composed of eight anti-parallel β-sheets and two α-helices with a Zn2+ ion in the active site. This model was considerably stable during 10 ns of molecular dynamics simulation. We observed a high diversity of bacteriocins in the Tucuruí-HPP, demonstrating that the environment is an inexhaustible source for prospecting these molecules. Finally, the zoocin model can be used for further studies of substrate binding and molecular mechanisms involving peptidoglycan degradation.

Self-Assembly of Black Cumin Oil-Based Nanoemulsion on Various Surfactants: A Molecular Dynamics Study

Black cumin is commonly used as traditional medicine due to its wide range of pharmacological potential. Black cumin oil (BCO) was often prepared as nanoemulsion to improve its solubility, stability, and bioavailability. This study was conducted to investigate the molecular behavior as well as structural evolution of BCO-surfactant systems during self-assembly micellization using molecular dynamics (MD) simulations. Several BCO constituents and variations of surfactants were employed to model BCO-surfactant systems. 50 ns of MD simulations were performed to elucidate their evolution of structures and physicochemical properties during formation. Results showed that BCO-tween20 and BCO-lecithin were able to form spherical-shaped micelles with the effective radii of 10.20 and 8.67 nm at the end of the simulation. Also, from the root mean square deviation and radius of gyration profile, it is showed that BCO-tween20 system was able to maintain the stability of its structure throughout the simulation. Results also revealed that self-assembly of BCO-surfactant systems were exothermic processes, confirming spontaneous nature upon formation

A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds.

Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has fomented a climate of fear worldwide due to its rapidly spreading nature, and high mortality rate. The World Health Organization (WHO) declared it as a global pandemic on 11th March, 2020. Many endeavors have been made to find appropriate medications to restrain the SARS CoV-2 infection from spreading but there is no specific antiviral therapy to date. However, a computer-aided drug design approach can be an alternative to identify probable drug candidates within a short time. SARS-CoV-2 main protease is a proven drug target, and it plays a pivotal role in viral replication and transcription. Methods: In this study, we identified a total of 114 essential oil compounds as a feasible anti-SARS-CoV-2 agent from several online reservoirs. These compounds were screened by incorporating ADMET profiling, molecular docking, and 50 ns of molecular dynamics simulation to identify potential drug candidates against the SARS-CoV-2 main protease. The crystallized SARS-CoV-2 main protease structure was collected from the RCSB PDB database (PDB ID 6LU7). Results: According to the results of the ADMET study, none of the compounds have any side effects that could reduce their druglikeness or pharmacokinetic properties. Out of 114 compounds, we selected bisabololoxide B, eremanthin, and leptospermone as our top drug candidates based on their higher binding affinity scores, and strong interaction with the Cys 145-His 41 catalytic dyad. Finally, the molecular dynamics simulation was implemented to evaluate the structural stability of the ligand-receptor complex. MD simulations disclosed that all the hits showed conformational stability compared to the positive control α-ketoamide. Conclusions: Our study showed that the top three hits might work as potential anti-SARS-CoV-2 agents, which can pave the way for discovering new drugs, but for experimental validation, they will require more in vivo trials.

In-silico study of seaweed secondary metabolites as AXL kinase inhibitors

AXL kinase is an attractive cancer target for drug design and it is involved in different cancers. A set of molecule databases with 1072 secondary metabolites from seaweeds were screened against the AXL kinase active site and eight molecules were shortlisted for further studies. From the docking analysis of the complexes, four molecules GA011, BE005, BC010, and BC005 are showing prominent binging. From the 100 ns of molecular dynamics simulations and ligand-bound complex MM-PBSA free energy analysis, two molecules BC010 (ΔG = −135.38 kJ/mol) and BE005 (ΔG = −141.72 kJ/mol) are showing molecule stability in the active site also showing very strong binding free energies. It suggests these molecules could be the potent molecules for AXL kinase.