Pediocin and analogous bacteriocins, valued for thermal stability, serve as versatile antimicrobials in the food sector. Improving their resilience at high temperatures and deriving derivatives not only benefit food production but also offer broad-spectrum antimicrobial potential in pharmaceuticals, spanning treatments for peptic ulcers, women\'s health, and novel anticancer agents. The study aims to create mutant peptides capable of establishing a third disulfide bond or enhanced through cysteine substitutions. This involves introducing additional Cys residues into the inherent structure of pediocin PA-1 to facilitate disulfide bond formation. Five mutants (Mut 1-5) were systematically generated with double Cys substitutions and assessed for thermal stability through MD simulations across temperatures (298-394 K). The most robust mutants (Mut 1, Mut 4-5) underwent extended analysis via MD simulations, comparing their structural stability, secondary structure, and surface accessibility to the reference Pediocin PA-1 molecule. This comprehensive assessment aims to understand how Cys substitutions influence disulfide bonds and the overall thermal stability of the mutant peptides. In silico analysis indicated that Mut 1 and Mut 5, along with the reference structure, lose their helical structure and one natural disulfide bond at high temperatures, and may impacting antimicrobial activity. Conversely, Mut 4 retained its helical structure and exhibited thermal stability similar to Pediocin PA-1. Pending further experimental validation, this study implies Mut 4 may have high stability and exceptional resistance to high temperatures, potentially serving as an effective antimicrobial alternative.

The Dopa Decarboxylase (DDC) gene plays an important role in the synthesis of biogenic amines such as dopamine, serotonin, and histamine. Non-synonymous single nucleotide polymorphisms (nsSNPs) in the DDC gene have been linked with various neurodegenerative disorders. In this study, a comprehensive in silico analysis of nsSNPs in the DDC gene was conducted to assess their potential functional consequences and associations with disease outcomes. Using publicly available databases, a complete list of nsSNPs in the DDC gene was obtained. 29 computational tools and algorithms were used to characterise the effects of these nsSNPs on protein structure, function, and stability. In addition, the population-based association studies were performed to investigate possible associations between specific nsSNPs and arthritis. Our research identified four novel DDC gene nsSNPs that have a major impact on the structure and function of proteins. Through molecular dynamics simulations (MDS), we observed changes in the stability of the DDC protein induced by specific nsSNPs. Furthermore, population-based association studies have revealed potential associations between certain DDC nsSNPs and various neurological disorders, including Parkinson\'s disease and dementia. The in silico approach used in this study offers insightful information about the functional effects of nsSNPs in the DDC gene. These discoveries provide insight into the cellular processes that underlie cognitive disorders. Furthermore, the detection of disease-associated nsSNPs in the DDC gene may facilitate the development of tailored and targeted therapy approaches.Communicated by Ramaswamy H. Sarma.

Self-assembling protein nanoparticles showed promise for vaccine design due to efficient antigen presentations and safety. However, the unpredictable formations of epitopes-fused protein assemblies remain challenging in the upstream design. This study suggests employing molecular dynamic (MD) simulations to investigate the assembly properties of Hepatitis B core protein (HBc) from thermodynamic perspectives. Eight HBc derivatives were expressed in E. coli, with their self-assembly properties characterised by high-performance liquid chromatography and transmission electron microscopy. MD simulations on the dimers, based on AlphaFold-predicted 3D structures, analysed the derivative at the atomic level. Results revealed that HBc derivatives can form dissociative polymers or large multi-subunit structures due to assembly failures. The instability of the dimer in aqueous solvents or inappropriate intradimer distances could cause major assembly failures. Polar solvation energies played a vital role too in forming assemble-incompetent dimers. Importantly, our study demonstrated that MD simulations on dimers can provide preliminary predictions on the assembly properties of HBc derivatives, thus aiding vaccine design by lowering the risk of self-assembling failures in engineered proteins.Communicated by Ramaswamy H. Sarma.

Amphiphilic peptides, such as Aß amyloids, can adsorb at an interface between two immiscible electrolyte solutions (ITIES). Based on previous work (vide infra), a hydrophilic/hydrophobic interface is used as a simple biomimetic system for studying drug interactions. The ITIES provides a 2D interface to study ion-transfer processes associated with aggregation, as a function of Galvani potential difference. Here, the aggregation/complexation behaviour of Aβ(1-42) is studied in the presence of Cu (II) ions, together with the effect of a multifunctional peptidomimetic inhibitor (P6). Cyclic and differential pulse voltammetry proved to be particularly sensitive to the detection of the complexation and aggregation of Aβ(1-42), enabling estimations of changes in lipophilicity upon binding to Cu (II) and P6. At a 1:1 ratio of Cu (II):Aβ(1-42), fresh samples showed a single DPV (Differential Pulse Voltammetry) peak half wave transfer potential (E1/2) at 0.40 V. Upon increasing the ratio of Cu (II) two-fold, fluctuations were observed in the DPVs, indicating aggregation. The approximate stoichiometry and binding properties of Aβ(1-42) during complexation with Cu (II) were determined by performing a differential pulse voltammetry (DPV) standard addition method, which showed two binding regimes. A pKa of 8.1 was estimated, with a Cu:Aβ1-42 ratio~1:1.7. Studies using molecular dynamics simulations of peptides at the ITIES show that Aβ(1-42) strands interact through the formation of β-sheet stabilised structures. In the absence of copper, binding/unbinding is dynamic, and interactions are relatively weak, leading to the observation of parallel and anti-parallel arrangements of β-sheet stabilised aggregates. In the presence of copper ions, strong binding occurs between a copper ion and histidine residues on two peptides. This provides a convenient geometry for inducing favourable interactions between folded β-sheet structures. Circular Dichroism spectroscopy (CD spectroscopy) was used to support the aggregation behaviour of the Aβ(1-42) peptides following the addition of Cu (II) and P6 to the aqueous phase.

DNA polymerases create complementary DNA strands in living cells and are crucial to genome transmission and maintenance. These enzymes possess similar human right-handed folds which contain thumb, fingers, and palm subdomains and contribute to polymerization activities. These enzymes are classified into seven evolutionary families, A, B, C, D, X, Y, and RT, based on amino acid sequence analysis and biochemical characteristics. Family A DNA polymerases exist in an extended range of organisms including mesophilic, thermophilic, and hyper-thermophilic bacteria, participate in DNA replication and repair, and have a broad application in molecular biology and biotechnology. In this study, we attempted to detect factors that play a role in the thermostability properties of this family member despite their remarkable similarities in structure and function. For this purpose, similarities and differences in amino acid sequences, structure, and dynamics of these enzymes have been inspected. Our results demonstrated that thermophilic and hyper-thermophilic enzymes have more charged, aromatic, and polar residues than mesophilic ones and consequently show further electrostatic and cation-pi interactions. In addition, in thermophilic enzymes, aliphatic residues tend to position in buried states more than mesophilic enzymes. These residues within their aliphatic parts increase hydrophobic core packing and therefore enhance the thermostability of these enzymes. Furthermore, a decrease in thermophilic cavities volumes assists in the protein compactness enhancement. Moreover, molecular dynamic simulation results revealed that increasing temperature impacts mesophilic enzymes further than thermophilic ones that reflect on polar and aliphatic residues surface area and hydrogen bonds changes.

Food preservation is a schematic and scientific procedure employed for the maintenance and improvement of food\'s quality, shelf life, and nutritional value. Although, on one hand, ancient conventional methods such as freezing, pasteurization, canning, and chemical methods have the potential to lengthen the shelf life of edible substances, but on the other hand, they can deteriorate its nutritional value as well. Present research focuses on the identification of promising bacteriocins against Pseudomonas fragi via subtractive proteomics pipeline as an alternative approach for food preservation. Bacteriocins are small peptides produced by certain microbes to naturally defend themselves by destroying other closely related bacteria residing in their neighborhood. P. fragi lies among the most notable microbes responsible for the elicitation of food spoilage. Due to increasing emergence and prevalence of multidrug resistance bacteria, there is a need to unravel novel drug targets, crucially involved in food decay process. Based on subtractive scrutinization, UDP-N-acetylglucosamine O-acyltransferase (LpxA) was chosen as promising therapeutic protein target that could play a significant role in progression of food spoilage. Subtilosin A, thuricin-CD, and mutacin B-NY266 were found as the most robust inhibitors of LpxA according to the molecular docking assay results. Molecular dynamic simulations and binding energy calculations via MM/PBSA method of LpxA and three top hit docked complexes, i.e., LpxA-subtilosin A, LpxA-thuricin-CD, and LpxA-mutacin B-NY266, revealed stability throughout simulations and ensured that shortlisted bacteriocins had strong affinity for LpxA.

In recent years, humanity has had to face a critical pandemic due to SARS-CoV-2. In the rapid search for effective drugs against this RNA-positive virus, the repurposing of already existing nucleotide/nucleoside analogs able to stop RNA replication by inhibiting the RNA-dependent RNA polymerase enzyme has been evaluated. In this process, a valid contribution has been the use of in silico experiments, which allow for a rapid evaluation of the possible effectiveness of the proposed drugs. Here we propose a molecular dynamic study to provide insight into the inhibition mechanism of Penciclovir, a nucleotide analog on the RNA-dependent RNA polymerase enzyme. Besides the presented results, in this article, for the first time, molecular dynamic simulations have been performed considering not only the RNA-dependent RNA polymerase protein, but also its cofactors (fundamental for RNA replication) and double-strand RNA.

Heparin and heparan sulfate are important glycosaminoglycans that can regulate the activities of many vital proteins, especially the fibroblast growth factor (FGF) family. Because FGF7 (KGF) has an important role in tissue repair and maintaining the integrity of the mucosal barrier, recombinant human keratinocyte growth factor (rhKGF, palifermin) has been approved for the treatment of wound healing and oral cavity. Due to heparin plays an important role in the KGF signaling pathway, a more detailed study of the drug-drug interactions (DDIs) between rhKGF and heparin at the atomic level and investigating their synergistic effect on each other in terms of biology, especially in silico, is necessary for a better understanding of DDIs. In this study, DDIs between rhKGF and low-molecular weight heparin types (LMWH) were investigated. In this regard, scrutiny of the influence of the synergistic heparin types on the structure and biostability of rhKGF is accomplished using computational methods such as molecular docking and molecular dynamic simulations (MDs). Subsequently, the motion behavior of rhKGF in interaction with LMWHs was evaluated based on eigenvectors by using principal component analysis (PCA). Also, the binding free energies of rhKGF-LMWH complexes were calculated by the molecular mechanics/Poisson-Boltzmann surface area (MM-BPSA) method. The result showed that rhKGF-idraparinux (-6.9 kcal/mol) and rhKGF-heparin (-6.0 kcal/mol) complexes had significant binding affinity as well as they had a more stable binding to rhKGF than to other LMWH during 100 ns simulation. However, in order to confirm the curative effect of these drugs, clinical trials must be done.

Evodiamine (EVO) and rutaecarpine (RUT) are the main active compounds of the traditional Chinese medicinal herb Evodia rutaecarpa. Here, we fully optimized the molecular geometries of EVO and RUT at the B3LYP/6-311++G (d, p) level of density functional theory. The natural population analysis (NPA) charges, frontier molecular orbitals, molecular electrostatic potentials, and the chemical reactivity descriptors for EVO and RUT were also investigated. Furthermore, molecular docking, molecular dynamics simulations, and the analysis of the binding free energies of EVO and RUT were carried out against the anticancer target topoisomerase 1 (TOP1) to clarify their anticancer mechanisms. The docking results indicated that they could inhibit TOP1 by intercalating into the cleaved DNA-binding site to form a TOP1−DNA−ligand ternary complex, suggesting that they may be potential TOP1 inhibitors. Molecular dynamics (MD) simulations evaluated the binding stability of the TOP1−DNA−ligand ternary complex. The calculation of binding free energy showed that the binding ability of EVO with TOP1 was stronger than that of RUT. These results elucidated the structure−activity relationship and the antitumor mechanism of EVO and RUT at the molecular level. It is suggested that EVO and RUT may be potential compounds for the development of new anticancer drugs.

The SARS-CoV-2 targets were evaluated for a set of FDA-approved drugs using a combination of drug repositioning and rigorous computational modeling methodologies such as molecular docking and molecular dynamics (MD) simulations followed by binding free energy calculations. Six FDA-approved drugs including, Ouabain, Digitoxin, Digoxin, Proscillaridin, Salinomycin and Niclosamide with promising anti-SARS-CoV-2 activity were screened in silico against four SARS-CoV-2 proteins-papain-like protease (PLpro), RNA-dependent RNA polymerase (RdRp), SARS-CoV-2 main protease (Mpro), and adaptor-associated kinase 1 (AAK1)-in an attempt to define their promising targets. The applied computational techniques suggest that all the tested drugs exhibited excellent binding patterns with higher scores and stable complexes compared to the native protein cocrystallized inhibitors. Ouabain was suggested to act as a dual inhibitor for both PLpro and Mpro enzymes, while Digitoxin bonded perfectly to RdRp. In addition, Salinomycin targeted PLpro. Particularly, Niclosamide was found to target AAK1 with greater affinity compared to the reference drug. Our study provides comprehensive molecular-level insights for identifying or designing novel anti-COVID-19 drugs.