What is the probability of a chance prediction of a protein structure with an rmsd of 6 A?

Inhibition of α-glucosidase activity by potential peptides derived from fermented spent coffee grounds

The control of α-glucosidase activity has been associated with managing diabetes. We previously identified three peptides with high bioactive indices derived from protein hydrolysates of fermented spent coffee grounds. In this study, the peptides YGF, GMCC, and RMYRY were synthesized and tested in vitro for their α-glucosidase inhibition activity, complemented by in silico analyses. Two of the three peptides significantly inhibited α-glucosidase activity, with the more efficient peptides being YGF and GMCC (0.42 mg/mL), resulting in decreased enzymatic activity of 95.31% and 89.79%, respectively. These peptides exhibited binding free energies with the α-glucosidase complex of -8.5 and - 6.6 kcal/mol, respectively, through hydrogen bonds and van der Waals interactions with amino acids from the active site. Pharmacokinetic analysis indicated that YGF and GMCC profiles were unrelated to toxicity. These results underscore the importance of focusing on food waste bioprocessing products to expand the range of alternatives that could aid in diabetes treatment.

In silico analysis of the interaction of de novo peptides derived from Salvia hispanica with anticancer targetsEvaluation of the anticancer potential of de novo peptides derived from Salvia hispanica through molecular docking

Cancer is one of the leading causes of death worldwide. Conventional cancer therapies are not selective to cancer cells resulting in serious side effects on patients. Thus, the need for complementary treatments that improve the patient's response to cancer therapy is highly important. To predict and evaluate the physicochemical characteristics and potential anticancer activity of the peptides identified from S. hispanica protein fraction <1 kDa through the use of in silico tools. Peptides derived from Salvia hispanica's protein fraction <1 kDa were identified and analyzed for the prediction of their physicochemical properties. The characterized peptide sequences were then submitted to a multi-criteria decision analysis to identify the peptides that possess the characteristics to potentially exert anticancer activity. Through molecular docking analysis, the potential anticancer activity of the Potentially Anticancer Peptide (PAP)-1, PAP-2, PAP-3, PAP-4, and PAP-5 was estimated by their binding interactions with cancer and apoptosis-related molecules. All five evaluated PAPs exhibited strong binding interactions (< -100 kcal/mol). However, PAP-3 showed the lowest binding free energies with several of the targets. Thus, PAP-3 shows potential to be used as a nutraceutical or ingredient for functional foods that adjuvate in cancer treatment. Conclusions: Through the molecular docking studies, the binding of the PAPs to target molecules of interest for cancer treatment was successfully simulated, from which PAP-3 exhibited the lowest binding free energies. Further in vitro and in vivo studies are required to validate the predictions obtained by the in silico analysis.Communicated by Ramaswamy H. Sarma.

In-silico study of protein-protein interactions in wheat blast using docking and molecular dynamics simulation approach

Despite advancements in agricultural research and the introduction of modern biotechnological and farming techniques, food security remains a significant issue. Although the efforts of farmers to meet the demands of a growing population, many plant diseases caused by pathogens, through their effects on cell division and tissue growth, lead to the annual loss of countless food crops. The recently emerged wheat blast fungus Magnaporthe oryzae pathotype Triticum (MoT) poses a significant danger to worldwide wheat cultivation. The fungus is a highly varied lineage of the M. oryzae, responsible for causing rice blast disease. In spite of being a significant challenge to successful wheat production in South America since 1985, the underlying biology of the wheat blast pathogen is still not fully understood. The initial outbreak of the wheat blast in South Asia had a severe impact on wheat production, resulting in a complete loss of yield in affected fields. For the purpose of enhancing disease management, it's vital to acquire a comprehensive comprehension of the infection biology of the fungus and its interaction with wheat plants on molecular levels. Host-pathogen protein interactions (HPIs) have the potential to reveal the pathogens' mechanism for overcoming the host organism. The current study delves into the interactions between the host plant wheat and MoT through protein-protein interactions, molecular docking, and 100 ns molecular dynamic simulations. This research uncovers the structural and functional basis of these proteins, leading to improved plant health and production.Communicated by Ramaswamy H. Sarma.

Diversity and selection analyses identify transmission-blocking antigens as the optimal vaccine candidates in Plasmodium falciparum

Background

A highly effective vaccine for malaria remains an elusive target, at least in part due to the under-appreciated natural parasite variation. This study aimed to investigate genetic and structural variation, and immune selection of leading malaria vaccine candidates across the Plasmodium falciparum's life cycle.

Methods

We analyzed 325 P. falciparum whole genome sequences from Zambia, in addition to 791 genomes from five other African countries available in the MalariaGEN Pf3k Rdatabase. Ten vaccine antigens spanning three life-history stages were examined for genetic and structural variations, using population genetics measures, haplotype network analysis, and 3D structure selection analysis.

Findings

Among the ten antigens analyzed, only three in the transmission-blocking vaccine category display P. falciparum 3D7 as the dominant haplotype. The antigens AMA1, CSP, MSP119 and CelTOS, are much more diverse than the other antigens, and their epitope regions are under moderate to strong balancing selection. In contrast, Rh5, a blood stage antigen, displays low diversity yet slightly stronger immune selection in the merozoite-blocking epitope region. Except for CelTOS, the transmission-blocking antigens Pfs25, Pfs48/45, Pfs230, Pfs47, and Pfs28 exhibit minimal diversity and no immune selection in epitopes that induce strain-transcending antibodies, suggesting potential effectiveness of 3D7-based vaccines in blocking transmission.

Interpretations

These findings offer valuable insights into the selection of optimal vaccine candidates against P. falciparum. Based on our results, we recommend prioritizing conserved merozoite antigens and transmission-blocking antigens. Combining these antigens in multi-stage approaches may be particularly promising for malaria vaccine development initiatives.

Funding

Purdue Department of Biological Sciences; Puskas Memorial Fellowship; National Institute of Allergy and Infectious Diseases (U19AI089680).

Computational analysis of non-synonymous single nucleotide polymorphism in the bovine PKLR geneComputational analysis of bovine PKLR gene

Pyruvate kinase (PKLR) is a potential candidate gene for milk production traits in cows. The main aim of this work is to investigate the potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the PKLR gene by using several computational tools. In silico tools including SIFT, Polyphen-2, SNAP2 and Panther indicated only 18 nsSNPs out of 170 were considered deleterious. The analysis of proteins' stability change due to amino acid substitution performed by the use of the I-mutant, MUpro, CUPSTAT, SDM and Dynamut confirmed that 9 nsSNPs decreased protein stability. ConSurf analysis predicted that all 18 nsSNPs were evolutionary moderately or highly conserved. Two different domains of PKLR protein were revealed by the InterPro tool with 12 nsSNPs positioned in the Pyruvate Kinase barrel domain and 6 nsSNP present in the Pyruvate Kinase C Terminal. The PKLR 3D model was predicted by MODELLER software and validated via Ramachandran plot and Prosa which indicated a good quality model. The analysis of energy minimizations for the native and mutated structures was performed by SWISS PDB viewer with GROMOS 96 program and showed that 3 structural and 4 functional residues had total energy higher than the native model. These findings indicate that these mutant structures (rs441424814, rs449326723, rs476805413, rs472263384, rs474320860, rs475521477, rs441633284) were less stable than the native model. Molecular Dynamics simulations were performed to confirm the impact of nsSNPs on the protein structure and function. The present study provides useful information about functional SNPs that have an impact on PKLR protein in cattle.Communicated by Ramaswamy H. Sarma.

Phylogeny and structural modeling of the transcription factor CsqR (YihW) from Escherichia coli

CsqR (YihW) is a local transcription factor that controls expression of yih genes involved in degradation of sulfoquinovose in Escherichia coli. We recently showed that expression of the respective gene cassette might be regulated by lactose. Here, we explore the phylogenetic and functional traits of CsqR. Phylogenetic analysis revealed that CsqR had a conserved Met25. Western blot demonstrated that CsqR was synthesized in the bacterial cell as two protein forms, 28.5 (CsqR-l) and 26 kDa (CsqR-s), the latter corresponding to start of translation at Met25. CsqR-s was dramatically activated during growth with sulfoquinovose as a sole carbon source, and displaced CsqR-l in the stationary phase during growth on rich medium. Molecular dynamic simulations revealed two possible states of the CsqR-s structure, with the interdomain linker being represented by either a disordered loop or an ɑ-helix. This helix allowed the hinge-like motion of the N-terminal domain resulting in a switch of CsqR-s between two conformational states, "open" and "compact". We then modeled the interaction of both CsqR forms with putative effectors sulfoquinovose, sulforhamnose, sulfoquinovosyl glycerol, and lactose, and revealed that they all preferred the same pocket in CsqR-l, while in CsqR-s there were two possible options dependent on the linker structure.

Unraveling the binding mechanisms of SARS-CoV-2 variants through molecular simulations

The emergence of SARS-CoV-2 variants like Delta (AY.29) and Omicron (EG.5) poses continued challenges for vaccines and therapeutics. Mutations in the viral spike protein are key in altering infectivity and immune evasion. This study uses computational modeling to investigate the molecular binding mechanisms between spike protein variants and the ACE2 host receptor. Using the MARTNI force field, coarse-grained molecular dynamics (CGMD) simulations and nudged elastic band (NEB) calculations explore spike-ACE2 interactions for the wild type, Delta variant, and Omicron variant. The simulations reveal Omicron has the strongest binding affinity at -128.35 ± 10.91 kcal/mol, followed by Delta and wild type. Key mutations in Delta and Omicron, like Q493R and Q498R, optimize electrostatic contacts, enhancing ACE2 interactions. The wild-type spike has the highest transition state energy barrier at 17.87 kcal/mol, while Delta has the lowest barrier at 9.21 kcal/mol. Despite slightly higher dual barriers, Omicron's increased binding energy lowers its overall barrier to rapidly bind ACE2. These findings provide residue-level insights into mutation effects on SARS-CoV-2 infectivity. The computational modeling elucidates mechanisms underlying spike-ACE2 binding kinetics, aiding the development of vaccines and therapies targeting emerging viral strains.

Structural and biophysical properties of a [4Fe4S] ferredoxin-like protein from Synechocystis sp. PCC 6803 with a unique two domain structure

Electron carrier proteins (ECPs), binding iron-sulfur clusters, are vital components within the intricate network of metabolic and photosynthetic reactions. They play a crucial role in the distribution of reducing equivalents. In Synechocystis sp. PCC 6803, the ECP network includes at least nine ferredoxins. Previous research, including global expression analyses and protein binding studies, has offered initial insights into the functional roles of individual ferredoxins within this network. This study primarily focuses on Ferredoxin 9 (slr2059). Through sequence analysis and computational modeling, Ferredoxin 9 emerges as a unique ECP with a distinctive two-domain architecture. It consists of a C-terminal iron‑sulfur binding domain and an N-terminal domain with homology to Nil-domain proteins, connected by a structurally rigid 4-amino acid linker. Notably, in contrast to canonical [2Fe2S] ferredoxins exemplified by PetF (ssl0020), which feature highly acidic surfaces facilitating electron transfer with photosystem I reaction centers, models of Ferredoxin 9 reveal a more neutral to basic protein surface. Using a combination of electron paramagnetic resonance spectroscopy and square-wave voltammetry on heterologously produced Ferredoxin 9, this study demonstrates that the protein coordinates 2×[4Fe4S]2+/1+ redox-active and magnetically interacting clusters, with measured redox potentials of -420 ± 9 mV and - 516 ± 10 mV vs SHE. A more in-depth analysis of Fdx9's unique structure and protein sequence suggests that this type of Nil-2[4Fe4S] multi-domain ferredoxin is well conserved in cyanobacteria, bearing structural similarities to proteins involved in homocysteine synthesis in methanogens.

Computational investigation of impact of Pb(II) and Ni(II) ions on hUNG enzyme: insights from molecular dynamics simulations

Human uracil DNA glycosylase (hUNG), a crucial player in the initiation of the base excision repair pathway, is susceptible to alterations in function and conformation induced by the accumulation of toxic metals. Despite the recognized impact of toxic metals on DNA repair enzymes, there exists a notable deficiency in theoretical investigations addressing this phenomenon. This study investigates the impact of toxic heavy metal ions, Pb(II) and Ni(II), on the stability of hUNG through molecular dynamics (MD) simulations. The initial analysis involved the identification of key cavities in the hUNG enzyme. Notably, the active site cavity emerged as a promising site for ligand binding. Subsequently, AutoDockTools software was employed to dock Pb(II) and Ni(II) onto the identified cavities, followed by extensive MD simulations. The MD analysis, encompassing parameters such as root mean square deviation, radius of gyration, solvent accessible surface area, hydrogen bond variations, Ramachandran plot, principal component analysis, and root mean square fluctuations, collectively revealed distinct alterations in the behavior of the enzyme upon complexation with Pb(II) and Ni(II). Interestingly, the enzyme exhibited enhanced structural stability, reduced flexibility, and modified hydrogen bonding patterns in the presence of these toxic metal ions. The observed limitation in structural flexibility implies a more rigid and stable conformation when the enzyme complex with Pb(II) and Ni(II) compared to its free form. This structural alteration may lead to a potential reduction in enzymatic activity, suggesting that toxic metal ions influence the functional dynamics of hUNG. These computational findings offer valuable insights into the molecular interactions between metal ions and enzymes.Communicated by Ramaswamy H. Sarma.

Comprehensive characterization of coding and non-coding single nucleotide polymorphisms of the Myoneurin (MYNN) gene using molecular dynamics simulation and docking approaches

Genome-wide association studies (GWAS) identified a coding single nucleotide polymorphism, MYNN rs10936599, at chromosome 3q. MYNN gene encodes myoneurin protein, which has been associated with several cancer pathogenesis and disease development processes. However, there needed to be a more detailed characterization of this polymorphism's (and other coding and non-coding polymorphisms) structural, functional, and molecular impact. The current study addressed this gap and analyzed different properties of rs10936599 and non-coding SNPs of MYNN via a thorough computational method. The variant, rs10936599, was predicted functionally deleterious by nine functionality prediction approaches, like SIFT, PolyPhen-2, and REVEL, etc. Following that, structural modifications were estimated through the HOPE server and Mutation3D. Moreover, the mutation was found in a conserved and active residue, according to ConSurf and CPORT. Further, the secondary structures were predicted, followed by tertiary structures, and there was a significant deviation between the native and variant models. Similarly, molecular simulation also showed considerable differences in the dynamic pattern of the wildtype and mutant structures. Molecular docking revealed that the variant binds with better docking scores with ligand NOTCH2. In addition to that, non-coding SNPs located at the MYNN locus were retrieved from the ENSEMBL database. These were found to disrupt the transcription factor binding regulatory regions; nonetheless, only two affect miRNA target sites. Again, eight non-coding variants were detected in the testes with normalized expression, whereas HaploReg v4.1 unveiled annotations for non-coding variants. In summary, in silico comprehensive characterization of coding and non-coding single nucleotide polymorphisms of MYNN gene will assist researchers to work on MYNN gene and establish their association with certain types of cancers.

Prediction of interface between regions of varying degrees of order or disorderness in intrinsically disordered proteins from dihedral angles

Intrinsically disordered proteins (IDPs) are proteins that do not form uniquely defined three-dimensional (3-D) structures. Experimental research on IDPs is difficult since they go against the traditional protein structure-function paradigm. Although there are several predictors of disorder based on amino acid sequences, but very limited based on the 3-D structures of proteins. Dihedral angles have a significant role in predicting protein structure because they establish a protein's backbone, which, coupled with its side chain, establishes its overall shape. Here, we have carried out atomistic Molecular Dynamics (MD) simulations on four different proteins: one ordered protein (Monellin), two partially disordered proteins (p53-TAD and Amyloid beta (Aβ1-42) peptide), and one completely disordered protein (Histatin 5). The MD simulation trajectories for the corresponding four proteins were used to conduct dihedral angle (ϕ and ѱ) analysis. Then, the average dihedral angles for each of the residues were calculated and plotted against the residue index. We noticed steep rises or falls in the average ϕ value at certain locations in the plot. These sudden shifts in the average ϕ value reflect the interface between regions of varying degrees of order or disorderness in intrinsically disordered proteins. Using this method, the probable conformer of a protein with a higher degree of disorder can be found among the ensembles of structures sampled during the MD simulations. The results of our study offer new understandings on precisely identifying regions of various degrees of disorder in intrinsically disordered proteins.Communicated by Ramaswamy H. Sarma.

Genome-wide identification, structural homology analysis, and evolutionary diversification of the phospholipase D gene family in the venom gland of three scorpion species

BACKGROUND

Venom phospholipase D (PLDs), dermonecrotic toxins like, are the major molecules in the crude venom of scorpions, which are mainly responsible for lethality and dermonecrotic lesions during scorpion envenoming. The purpose of this study was fivefold: First, to identify transcripts coding for venom PLDs by transcriptomic analysis of the venom glands from Androctonus crassicauda, Hottentotta saulcyi, and Hemiscorpius lepturus; second, to classify them by sequence similarity to known PLDs and motif extraction method; third, to characterize scorpion PLDs; fourth to structural homology analysis with known dermonecrotic toxins; and fifth to investigate phylogenetic relationships of the PLD proteins.

RESULTS

We found that the venom gland of scorpions encodes two PLD isoforms: PLD1 ScoTox-beta and PLD2 ScoTox-alpha I. Two highly conserved regions shared by all PLD1s beta are GAN and HPCDC (HX2PCDC), and the most important conserved regions shared by all PLD2s alpha are two copies of the HKDG (HxKx4Dx6G) motif. We found that PLD1 beta is a 31-43 kDa acidic protein containing signal sequences, and PLD2 alpha is a 128 kDa basic protein without known signal sequences. The gene structures of PLD1 beta and PLD2 alpha contain 6 and 21 exons, respectively. Significant structural homology and similarities were found between the modeled PLD1 ScoTox-beta and the crystal structure of dermonecrotic toxins from Loxosceles intermedia.

CONCLUSIONS

This is the first report on identifying PLDs from A. crassicauda and H. saulcyi venom glands. Our work provides valuable insights into the diversity of scorpion PLD genes and could be helpful in future studies on recombinant antivenoms production.

Secretory protein Rv1987, a 'probable chitinase' from Mycobacterium tuberculosis is a novel chitin and cellulose binding protein lacking enzymatic function

Bacterial chitinases serve to hydrolyse chitin as food source or as defence mechanism. Given that chitin is not produced by mammals, it is intriguing that Mycobacterium tuberculosis, an exclusively human pathogen harbours Rv1987, a probable chitinase and secretes it. Interestingly genes annotated as chitinases are widely distributed among Mycobacterium tuberculosis complex species, clinical isolates and other human pathogens M. abscessus and M. ulcerans. However, Mycobacterial chitinases are not characterized and hence the functions remain unknown. In the present study, we show that Rv1987 is a chitin and cellulose binding protein lacking enzymatic activity in contrary to its current annotation. Further, we show Rv1987 has moon lighting functions in M. tuberculosis pathobiology signifying roles of bacterial cellulose binding clusters in infections.

Asuc_0142 of Actinobacillus succinogenes 130Z is the l-aspartate/C4-dicarboxylate exchanger DcuA

Anaerobic bacteria often use antiporters DcuB (malate/succinate antiport) or DcuA (l-aspartate/succinate antiport) for the excretion of succinate during fumarate respiration. The rumen bacterium Actinobacillus succinogenes is able to produce large amounts of succinate by fumarate respiration, using the DcuB-type transporter DcuE for l-malate/succinate antiport. Asuc_0142 was annotated as a second DcuB-type transporter. Deletion of Asuc_0142 decreased the uptake rate for l-[14C]aspartate into A. succinogenes cells. Properties of transport by heterologously expressed Asuc_0142 were investigated in an Escherichia coli mutant deficient of anaerobic C4DC transporters. Expression of Asuc_0142 resulted in high uptake activity for l-[14C]fumarate or l-[14C]aspartate, but the former showed a strong competitive inhibition by l-aspartate. In E. coli loaded with l-[14C]aspartate, [14C]succinate or [14C]fumarate, extracellular C4DCs initiated excretion of the intracellular substrates, with a preference for l-aspartateex/succinatein or l-aspartateex/fumaratein antiport. These findings indicate that Asuc_0142 represents a DcuA-type transporter for l-aspartate uptake and l-aspartateex/C4DCin antiport, differentiating it from the DcuB-type transporter DcuE for l-malateex/succinatein antiport. Sequence analysis and predicted structural characteristics confirm structural similarity of Asuc_0142 to DcuA, and Asuc_0142 was thus re-named as DcuAAs. The bovine rumen fluid contains l-aspartate (99.6 µM), whereas fumarate and l-malate are absent. Therefore, bovine rumen colonisers depend on l-aspartate as an exogenous substrate for fumarate respiration. A. succinogenes encodes HemG (protoporphyrinogen oxidase) and PyrD (dihydroorotate dehydrogenase) for haem and pyrimidine biosynthesis. The enzymes require fumarate as an electron acceptor, suggesting an essential role for l-aspartate, DcuAAs, and fumarate respiration for A. succinogenes growing in the bovine rumen.

Breaking the 'don't eat me' signal: in silico design of CD47-directed peptides for cancer immunotherapy

The leading cause of death worldwide is cancer. Although there are various therapies available to treat cancer, finding a successful one can be like searching for a needle in a haystack. Immunotherapy appears to be one of those needles in the haystack of cancer treatment. Immunotherapeutic agents enhance the immune response of the patient's body to tumor cells. One of the immunotherapeutic targets, Cluster of Differentiation 47 (CD47), releases the "don't eat me" signal when it binds to its receptor, Signal Regulatory Protein (SIRPα). Tumor cells use this signal to circumvent the immune system, rendering it ineffective. To stop tumor cells from releasing the "don't eat me" signal, the CD47-SIRPα interaction is specifically targeted in this study. To do so, in silico peptides were designed based on the structural analysis of the interaction between two proteins using point mutations on the interacting residues with the other amino acids. The peptide library was designed and docked on SIRPα using computational tools. Later on, after analyzing the docked complex, the best of them was selected for MD simulation studies of 100 ns. Further analysis after MD studies was carried out to determine the possible potential anti-SIRPα peptides.

Unveiling the Role of Sorghum RPAP3 in the Function of R2TP Complex: Insights into Protein Assembly in Plants

The chaperone R2TP has multiple subunits that assist in the proper folding, assembly, and stabilization of various protein complexes in cells and its study can offer valuable insights into the regulation and maintenance of protein assemblies in plant systems. The 'T' component of R2TP is Tah1 in yeast, consisting of 111 residues, while its counterpart in humans is RPAP3, with 665 residues. RPAP3 acts as a co-chaperone of Hsp90 and facilitates interactions between RUVBL proteins and other complex components, enhancing the recruitment of client proteins by the R2TP complex. These facts further underscore the relevance of studying this complex in different organisms. The putative gene corresponding to the RPAP3 in Sorghum bicolor, a monocotyledon plant, was cloned, and the protein (396 residues) purified for biochemical characterization. SbRPAP3 exists as a folded monomer and has a RPAP3 domain, which is present in human RPAP3 but absent in yeast Tah1. SbRPAP3 retains its functional capabilities, including binding with RUVBLs, Hsp90, and Hsp70. By elucidating the role of RPAP3 in plant R2TP complex, we can further comprehend the molecular mechanisms underlying plant-specific protein assembly and contribute to advancements in plant biology and biotechnological applications.

A computational peptide model induces cancer cells' apoptosis by docking Kringle 5 to GRP78

BACKGROUND

Cells can die through a process called apoptosis in both pathological and healthy conditions. Cancer development and progression may result from abnormal apoptosis. The 78-kDa glucose-regulated protein (GRP78) is increased on the surface of cancer cells. Kringle 5, a cell apoptosis agent, is bound to GRP78 to induce cancer cell apoptosis. Kringle 5 was docked to GRP78 using ClusPro 2.0. The interaction between Kringle 5 and GRP78 was investigated.

RESULTS

The interacting amino acids were found to be localized in three areas of Kringle 5. The proposed peptide is made up of secondary structure amino acids that contain Kringle 5 interaction residues. The 3D structure of the peptide model amino acids was created using the PEP-FOLD3 web tool.

CONCLUSIONS

The proposed peptide completely binds to the GRP78 binding site on the Kringle 5, signaling that it might be effective in the apoptosis of cancer cells.

Beta, Delta, and Omicron, Deadliest Among SARS-CoV-2 Variants: A Computational Repurposing Approach

SARS-CoV-2 has been highly susceptible to mutations since its emergence in Wuhan, China, and its subsequent propagation due to containing an RNA as its genome. The emergence of variants with improved transmissibility still poses a grave threat to global health. The spike protein mutation is mainly responsible for higher transmissibility and risk severity. This study retrieved SARS-CoV-2 variants structural and nonstructural proteins (NSPs) sequences from several geographic locations, including Africa, Asia, Europe, Oceania, and North and South America. First, multiple sequence alignments with BioEdit and protein homology modeling were performed using the SWISS Model. Then the structure visualization and structural analysis were performed by superimposing against the Wuhan sequence by Pymol to retrieve the RMSD values. Sequence alignment revealed familiar, uncommon regional among variants and, interestingly, a few unique mutations in Beta, Delta, and Omicron. Structural analysis of such unique mutations revealed that they caused structural deviations in Beta, Delta, and Omicron spike proteins. In addition, these variants were more severe in terms of hospitalization, sickness, and higher mortality, which have a substantial relationship with the structural deviations because of those unique mutations. Such evidence provides insight into the SARS-CoV-2 spike protein vulnerability toward mutation and their structural and functional deviations, particularly in Beta, Delta, and Omicron, which might be the cause of their broader coverage. This knowledge can help us with regional vaccine strain selection, virus pathogenicity testing, diagnosis, and treatment with more specific vaccines.

Molecular dynamic and bioinformatic studies of metformin-induced ACE2 phosphorylation in the presence of different SARS-CoV-2 S protein mutations

The SARS-CoV-2 infection activates host kinases and causes high phosphorylation in both the host and the virus. There were around 70 phosphorylation sites found in SARS-CoV-2 viral proteins. Besides, almost 15,000 host phosphorylation sites were found in SARS-CoV-2-infected cells. COVID-19 is thought to enter cells via the well-known receptor Angiotensin-Converting Enzyme 2 (ACE2) and the serine protease TMPRSS2. Substantially, the COVID-19 infection doesn't induce phosphorylation of the ACE2 receptor at Serin-680(s680). Metformin's numerous pleiotropic properties and extensive use in medicine including COVID-19, have inspired experts to call it the "aspirin of the twenty-first century". Metformin's impact on COVID-19 has been verified in clinical investigations via ACE2 receptor phosphorylation at s680. In the infection of COVID-19, sodium-dependent transporters including the major neutral amino acid (B0AT1) is regulated by ACE2. The structure of B0AT1 complexing with the COVID-19 receptor ACE2 enabled significant progress in the creation of mRNA vaccines. We aimed to study the impact of the interaction of the phosphorylation form of ACE2-s680 with wild-type (WT) and different mutations of SARS-CoV-2 infection such as delta, omicron, and gamma (γ) on their entrance of host cells as well as the regulation of B0AT1by the SARS-CoV-2 receptor ACE2. Interestingly, compared to WT SARS-CoV-2, ACE2 receptor phosphorylation at s680 produces conformational alterations in all types of SARS-CoV-2. Furthermore, our results showed for the first time that this phosphorylation significantly influences ACE2 sites K625, K676, and R678, which are key mediators for ACE2-B0AT1 complex.

Agomelatine's antiglycoxidative action-In vitro and in silico research and systematic literature review

Introduction

Agomelatine is an atypical antidepressant drug enhancing norepinephrine and dopamine liberation; nevertheless, additional mechanisms are considered for the drug's pharmacological action. Since protein glycoxidation plays a crucial role in depression pathogenesis, agomelatine's impact on carbonyl/oxidative stress was the research purpose.

Methods

Reactive oxygen species scavenging (hydroxyl radical, hydrogen peroxide, and nitrogen oxide) and antioxidant capacity (2,2-diphenyl-1-picrylhydrazyl radical and ferrous ion chelating assays) of agomelatine were marked. Agomelatine's antiglycoxidation properties were assayed in sugars (glucose, fructose, and galactose) and aldehydes- (glyoxal and methylglyoxal) glycated bovine serum albumin (BSA). Aminoguanidine and α-lipoic acid were used as standard glycation/oxidation inhibitors.

Results

Agomelatine did not show meaningful scavenging/antioxidant capacity vs. standards. Sugars/aldehydes increased glycation (↑kynurenine, ↑N-formylkynurenine, ↑dityrosine, ↑advanced glycation end products, and ↑β-amyloid) and oxidation (↑protein carbonyls and ↑advanced oxidation protein products) parameters in addition to BSA. Standards restored BSA baselines of glycation and oxidation markers, unlike agomelatine which sometimes even intensifies glycation above BSA + glycators levels. Molecular docking analysis of agomelatine in BSA demonstrated its very weak binding affinity.

Discussion

Agomelatine's very low affinity to the BSA could proclaim non-specific bonding and simplify attachment of glycation factors. Thereby, the drug may stimulate brain adaptation to carbonyl/oxidative stress as the systematic review indicates. Moreover, the drug's active metabolites could exert an antiglycoxidative effect.

Creatinase: Using Increased Entropy to Improve the Activity and Thermostability

Improving protein thermostability in mutagenesis-based enzyme engineering was often achieved by enhancing interresidue interactions via mutation to increase the enthalpy penalty of unfolding. However, this approach may trade off the functional activity due to the loss of structural flexibility of the biomolecule. Here, by performing X-ray crystallography, enzymatic kinetic experiments, neutron scattering, and thermodynamical measurements, we compared the structures, catalytic behaviors, dynamics, and thermostability between a wild-type creatinase and its four-point mutant. We found that the mutant is an entropy-driven thermostable protein with higher structural flexibility, i.e., higher conformational entropy, in the folded state compared to the wild type. The increased conformational entropy of the mutant in the folded state can reduce the entropy gain during unfolding and thus renders it greater thermostability. Moreover, the increased structural flexibility, particularly around the catalytic site, can broaden the mutant's working temperature range and considerably improve its activity at ambient conditions, which is crucial for its application in diagnosing kidney diseases. Complementary all-atom molecular dynamics simulations indicated that the four mutations replaced several of the strong interresidue interactions (electrostatic interactions and hydrogen bonds) with weak hydrophobic interactions. These substitutions not only release the structural flexibility to promote the thermostability and enzymatic activity of the protein but they also preserve the protein structure from collapsing. Our findings may pave a route for the entropy-driven strategy to design proteins with high thermostability and activity.

An Insight into the Essential Role of Carbohydrate-Binding Modules in Enzymolysis of Xanthan

To date, due to the low accessibility of enzymes to xanthan substrates, the enzymolysis of xanthan remains deficient, which hinders the industrial production of functional oligoxanthan. To enhance the enzymatic affinity against xanthan, the essential role of two carbohydrate binding modules—MiCBMx and PspCBM84, respectively, derived from Microbacterium sp. XT11 and Paenibacillus sp. 62047—in catalytic properties of endotype xanthanase MiXen were investigated for the first time. Basic characterizations and kinetic parameters of different recombinants revealed that, compared with MiCBMx, PspCBM84 dramatically increased the thermostability of endotype xanthanase, and endowed the enzyme with higher substrate affinity and catalytic efficiency. Notably, the activity of endotype xanthanase was increased by 16 times after being fused with PspCBM84. In addition, the presence of both CBMs obviously enabled endotype xanthanase to produce more oligoxanthan, and xanthan digests prepared by MiXen-CBM84 showed better antioxidant activity due to the higher content of active oligosaccharides. The results of this work lay a foundation for the rational design of endotype xanthanase and the industrial production of oligoxanthan in the future.

In silico prediction of potential inhibitors for SARS-CoV-2 Omicron variant using molecular docking and dynamics simulation-based drug repurposing

BACKGROUND

In November 2021, variant B.1.1.529 of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified by the World Health Organization (WHO) and designated Omicron. Omicron is characterized by a high number of mutations, thirty-two in total, making it more transmissible than the original virus. More than half of those mutations were found in the receptor-binding domain (RBD) that directly interacts with human angiotensin-converting enzyme 2 (ACE2). This study aimed to discover potent drugs against Omicron, which were previously repurposed for coronavirus disease 2019 (COVID-19). All repurposed anti-COVID-19 drugs were compiled from previous studies and tested against the RBD of SARS-CoV-2 Omicron.

METHODS

As a preliminary step, a molecular docking study was performed to investigate the potency of seventy-one compounds from four classes of inhibitors. The molecular characteristics of the best-performing five compounds were predicted by estimating the drug-likeness and drug score. Molecular dynamics simulations (MD) over 100 ns were performed to inspect the relative stability of the best compound within the Omicron receptor-binding site.

RESULTS

The current findings point out the crucial roles of Q493R, G496S, Q498R, N501Y, and Y505H in the RBD region of SARS-CoV-2 Omicron. Raltegravir, hesperidin, pyronaridine, and difloxacin achieved the highest drug scores compared with the other compounds in the four classes, with values of 81%, 57%, 18%, and 71%, respectively. The calculated results showed that raltegravir and hesperidin had high binding affinities and stabilities to Omicron with ΔGbinding of - 75.7304 ± 0.98324 and - 42.693536 ± 0.979056 kJ/mol, respectively. Further clinical studies should be performed for the two best compounds from this study.

Analogue and structure based approaches for modelling HIV-1 integrase inhibitors

A set of 220 inhibitors belonging to different structure classes and having HIV-1 integrase activity were collected along with their experimental pIC50 values. Geometries of all the inhibitors were fully optimized using B3LYP/6-31 + G(d) level of theory. These ligands were docked against 4 different HIV-1 integrase receptors (PDB IDs: 4LH5, 5KRS, 3ZSQ and 3ZSV). 30 docked poses were generated for all 220 inhibitors and ligand interaction of the first docked pose and the docked pose with the highest score were analysed. Residue GLU170 of 4LH5 receptor shows the highest number of interactions followed by ALA169, GLN168, HIS171 and ASP167 residues. Hydrogen bonding and stacking are mainly responsible for the interactions of these inhibitors with the receptor. We performed Molecular Dynamics (MD) simulation to observe the root-mean-square deviation (RMSD), for measure the average change of displacement between the atoms for a particular frame with respect to a reference and The Root Mean Square Fluctuation (RMSF) for characterization of local changes along the protein chain of the docked complexes. Analogue based models were generated to predict the pIC50 values for integrase inhibitors using various types of descriptors such as constitutional, geometrical, topological, quantum chemical and docking based descriptors. The best models were selected on the basis of statistical parameters and were validated by training and test set division. A few new inhibitors were designed on the basis of structure activity relationship and their pIC50 values were predicted using the generated models. All the designed new inhibitors a very high potential and may be used as potent inhibitors of HIV integrase. These models may be useful for further design and development of new and potent HIV integrase inhibitors.Communicated by Ramaswamy H. Sarma.

Genome-wide analysis of the Populus trichocarpa laccase gene family and functional identification of PtrLAC23

INTRODUCTION

Biofuel is a kind of sustainable, renewable and environment friendly energy. Lignocellulose from the stems of woody plants is the main raw material for "second generation biofuels". Lignin content limits fermentation yield and is therefore a major obstacle in biofuel production. Plant laccase plays an important role in the final step of lignin formation, which provides a new strategy for us to obtain ideal biofuels by regulating the expression of laccase genes to directly gain the desired lignin content or change the composition of lignin.

METHODS

Multiple sequence alignment and phylogenetic analysis were used to classify PtrLAC genes; sequence features of PtrLACs were revealed by gene structure and motif composition analysis; gene duplication, interspecific collinearity and Ka/Ks analysis were conducted to identify ancient PtrLACs; expression levels of PtrLAC genes were measured by RNA-Seq data and qRT-PCR; domain analysis combine with cis-acting elements prediction together showed the potential function of PtrLACs. Furthermore, Alphafold2 was used to simulate laccase 3D structures, proLAC23::LAC23-eGFP transgenic Populus stem transects were applied to fluorescence observation.

RESULTS

A comprehensive analysis of the P. trichocarpa laccase gene (PtLAC) family was performed. Some ancient PtrLAC genes such as PtrLAC25, PtrLAC19 and PtrLAC41 were identified. Gene structure and distribution of conserved motifs clearly showed sequence characteristics of each PtrLAC. Combining published RNA-Seq data and qRT-PCR analysis, we revealed the expression pattern of PtrLAC gene family. Prediction results of cis-acting elements show that PtrLAC gene regulation was closely related to light. Through above analyses, we selected 5 laccases and used Alphafold2 to simulate protein 3D structures, results showed that PtrLAC23 may be closely related to the lignification. Fluorescence observation of proLAC23::LAC23-eGFP transgenic Populus stem transects and qRT-PCR results confirmed our hypothesis again.

DISCUSSION

In this study, we fully analyzed the Populus trichocarpa laccase gene family and identified key laccase genes related to lignification. These findings not only provide new insights into the characteristics and functions of Populus laccase, but also give a new understanding of the broad prospects of plant laccase in lignocellulosic biofuel production.

In silico screening and epitope mapping of leptospiral outer membrane protein-Lsa46

Leptospirosis is one of the neglected diseases caused by the spirochete, Leptospira interrogans. Leptospiral surface adhesion (Lsa) proteins are surface exposed outer membrane proteins present in the pathogen. It acts as laminin and plasminogen binding proteins which enable them to infect host cells. The major target for the development of vaccine in the current era focuses on surface exposed outer membrane proteins, as they can induce strong and fast immune response in hosts. Therefore, the present study mapped the potential epitopes of the Leptospiral outer membrane proteins, mainly the surface adhesion proteins. Protein sequence analysis of Lsa proteins was done by in silico methods. The primary protein sequence analysis revealed Lsa46 as a suitable target which can be a potent Leptospiral vaccine candidate. Its structure was modelled by threading based method in I-TASSER server and validated by Ramachandran plot. The predicted epitope's interactions with human IgG, IgM(Fab) and T-cell receptor TCR(αβ) were performed by molecular docking studies using Biovia Discovery studio 2018. One of the predicted B-cell epitopes and the IgG showed desirable binding interactions, while four of the predicted B-cell epitopes and T-cell epitopes showed desirable binding interactions with IgM and TCR respectively. The molecular dynamic simulation studies carried out with the molecular docked complexes gave minimized energies indicating stable interactions. The structural analysis of the entire simulated complex showed a stable nature except for one of the Epitope-IgM complex. Further the binding free energy calculation of eight receptor-ligand complex predicted them energetically stable. The results of the study help in elucidating the structural and functional characterization of Lsa46 for epitope-based vaccine design.Communicated by Ramaswamy H. Sarma.

The Effects of One-Point Mutation on the New Delhi Metallo Beta-Lactamase-1 Resistance toward Carbapenem Antibiotics and β-Lactamase Inhibitors: An In Silico Systematic Approach

Antibiotic resistance has been becoming more and more critical due to bacteria's evolving hydrolysis enzymes. The NDM-1 enzyme could hydrolyze not only carbapenems but also most of β-lactam's antibiotics and inhibitors. In fact, variant strains could impose a high impact on the resistance of bacteria producing NDM-1. Although previous studies showed the effect of some variants toward antibiotics and inhibitors binding, there has been no research systematically evaluating the effects of alternative one-point mutations on the hydrolysis capacity of NDM-1. This study aims to identify which mutants could increase or decrease the effectiveness of antibiotics and β-lactamase inhibitors toward bacteria. Firstly, 35 different variants with a high probability of emergence based on the PAM-1 matrix were constructed and then docked with 5 ligands, namely d-captopril, l-captopril, thiorphan, imipenem, and meropenem. The selected complexes underwent molecular dynamics simulation and free energy binding estimation, with the results showing that the substitutions at residues 122 and 124 most influenced the binding ability of NDM-1 toward inhibitors and antibiotics. The H122R mutant decreases the binding ability between d-captopril and NDM-1 and diminishes the effectiveness of this antibiotic toward Enterobacteriaceae. However, the H122R mutant has a contrary impact on thiorphan, which should be tested in vitro and in vivo in further experiments.

Collaborative Assessment of Molecular Geometries and Energies from the Open Force Field

Force fields form the basis for classical molecular simulations, and their accuracy is crucial for the quality of, for instance, protein-ligand binding simulations in drug discovery. The huge diversity of small-molecule chemistry makes it a challenge to build and parameterize a suitable force field. The Open Force Field Initiative is a combined industry and academic consortium developing a state-of-the-art small-molecule force field. In this report, industry members of the consortium worked together to objectively evaluate the performance of the force fields (referred to here as OpenFF) produced by the initiative on a combined public and proprietary dataset of 19,653 relevant molecules selected from their internal research and compound collections. This evaluation was important because it was completely blind; at most partners, none of the molecules or data were used in force field development or testing prior to this work. We compare the Open Force Field "Sage" version 2.0.0 and "Parsley" version 1.3.0 with GAFF-2.11-AM1BCC, OPLS4, and SMIRNOFF99Frosst. We analyzed force-field-optimized geometries and conformer energies compared to reference quantum mechanical data. We show that OPLS4 performs best, and the latest Open Force Field release shows a clear improvement compared to its predecessors. The performance of established force fields such as GAFF-2.11 was generally worse. While OpenFF researchers were involved in building the benchmarking infrastructure used in this work, benchmarking was done entirely in-house within industrial organizations and the resulting assessment is reported here. This work assesses the force field performance using separate benchmarking steps, external datasets, and involving external research groups. This effort may also be unique in terms of the number of different industrial partners involved, with 10 different companies participating in the benchmark efforts.

Potential antiviral peptides targeting the SARS-CoV-2 spike protein

BACKGROUND

The coronavirus disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection became an international pandemic and created a public health crisis. The binding of the viral Spike glycoprotein to the human cell receptor angiotensin-converting enzyme 2 (ACE2) initiates viral infection. The development of efficient treatments to combat coronavirus disease is considered essential.

METHODS

An in silico approach was employed to design amino acid peptide inhibitor against the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein. The designed inhibitor (SARS-CoV-2 PEP 49) consists of amino acids with the α1 helix and the β4 - β5 sheets of ACE2. The PEP-FOLD3 web tool was used to create the 3D structures of the peptide amino acids. Analyzing the interaction between ACE2 and the RBD of the Spike protein for three protein data bank entries (6M0J, 7C8D, and 7A95) indicated that the interacting amino acids were contained inside two regions of ACE2: the α1 helical protease domain (PD) and the β4 - β5 sheets.

RESULTS

Molecular docking analysis of the designed inhibitor demonstrated that SARS-CoV-2 PEP 49 attaches directly to the ACE2 binding site of the Spike protein with a binding affinity greater than the ACE2, implying that the SARS-CoV-2 PEP 49 model may be useful as a potential RBD binding blocker.

In silico analysis of SARS-CoV-2 spike protein N501Y and N501T mutation effects on human ACE2 binding

The SARS-CoV-2 is an RNA-based virus and the most vital step of its survival is the attachment to hACE2 through its spike protein. Although SARS-CoV-2 has the ability to maintain high accurate replication and it can be accepted as a low mutation risked virus, it already showed more than nine thousand mutations in spike protein, of which 44 mutations are located within a 3.2 Å interacting distance from the hACE2 receptor. Mutations on spike protein, N501Y and N501T raised serious concerns for higher transmissibility and resistance towards current vaccines. In the current study, the mutational outcomes of N501Y and N501T on the hACE2-SARS CoV-2 spike protein complexation were analyzed by employing all-atom classic molecular dynamics (MD) simulations. These simulations revealed that both N501Y and N501T mutations increased the binding strength of spike protein to the host hACE2, predicted by binding free energy analysis via MM/GBSA rescoring scheme. This study highlights the importance of energy-based analysis for identifying mutational outcomes and will shed light on handling long-term and effective treatment strategies including repurposing anti-viral drugs, anti-SARS-CoV-2 antibodies, vaccines, and antisense based-therapies.

Targeting Natural Plant Metabolites for Hunting SARS-CoV-2 Omicron BA.1 Variant Inhibitors: Extraction, Molecular Docking, Molecular Dynamics, and Physicochemical Properties Study

(1) Background: SARS-CoV-2 Omicron BA.1 is the most common variation found in most countries and is responsible for 99% of cases in the United States. To overcome this challenge, there is an urgent need to discover effective inhibitors to prevent the emerging BA.1 variant. Natural products, particularly flavonoids, have had widespread success in reducing COVID-19 prevalence. (2) Methods: In the ongoing study, fifteen compounds were annotated from Echium angustifolium and peach (Prunus persica), which were computationally analyzed using various in silico techniques. Molecular docking calculations were performed for the identified phytochemicals to investigate their efficacy. Molecular dynamics (MD) simulations over 200 ns followed by molecular mechanics Poisson–Boltzmann surface area calculations (MM/PBSA) were performed to estimate the binding energy. Bioactivity was also calculated for the best components in terms of drug likeness and drug score. (3) Results: The data obtained from the molecular docking study demonstrated that five compounds exhibited remarkable potency, with docking scores greater than −9.0 kcal/mol. Among them, compounds 1, 2 and 4 showed higher stability within the active site of Omicron BA.1, with ΔG_binding values of −49.02, −48.07, and −67.47 KJ/mol, respectively. These findings imply that the discovered phytoconstituents are promising in the search for anti-Omicron BA.1 drugs and should be investigated in future in vitro and in vivo research.

Structural analysis of PpSP15 and PsSP9 sand fly salivary proteins designed with a self-cleavable linker as a live vaccine candidate against cutaneous leishmaniasis

Background

Leishmania parasites are deposited in the host through sand fly bites along with sand fly saliva. Therefore, salivary proteins are promising vaccine candidates for controlling leishmaniasis. Herein, two immunogenic salivary proteins, PpSP15 from Phlebotomus papatasi and PsSP9 from Phlebotomus sergenti, were selected as vaccine candidates to be delivered by live Leishmania tarentolae as vector. The stepwise in silico protocol advantaged in this study for multi-protein design in L. tarentolae is then described in detail.

Methods

All possible combinations of two salivary proteins, PpSP15 and PsSP9, with or without T2A peptide were designed at the mRNA and protein levels. Then, the best combination for the vaccine candidate was selected based on mRNA and protein stability along with peptide analysis.

Results

At the mRNA level, the most favored secondary structure was PpSP15-T2A-PsSP9. At the protein level, the refined three-dimensional models of all combinations were structurally valid; however, local quality estimation showed that the PpSp15-T2A-PsSP9 fusion had higher stability for each amino acid position, with low root-mean-square deviation (RMSD), compared with the original proteins. In silico evaluation confirmed the PpSP15-T2A-PsSP9 combination as a good Th1-polarizing candidate in terms of high IFN-γ production and low IL-10/TGF-β ratio in response to three consecutive immunizations. Potential protein expression was then confirmed by Western blotting.

Conclusions

The approach presented herein is among the first studies to have privileged protein homology modeling along with mRNA analysis for logical live vaccine design-coding multi-proteins.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05437-x.

Exploring CRD mobility during RAS/RAF engagement at the membrane

During the activation of mitogen-activated protein kinase (MAPK) signaling, the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF bind to active RAS at the plasma membrane. The orientation of RAS at the membrane may be critical for formation of the RAS-RBDCRD complex and subsequent signaling. To explore how RAS membrane orientation relates to the protein dynamics within the RAS-RBDCRD complex, we perform multiscale coarse-grained and all-atom molecular dynamics (MD) simulations of KRAS4b bound to the RBD and CRD domains of RAF-1, both in solution and anchored to a model plasma membrane. Solution MD simulations describe dynamic KRAS4b-CRD conformations, suggesting that the CRD has sufficient flexibility in this environment to substantially change its binding interface with KRAS4b. In contrast, when the ternary complex is anchored to the membrane, the mobility of the CRD relative to KRAS4b is restricted, resulting in fewer distinct KRAS4b-CRD conformations. These simulations implicate membrane orientations of the ternary complex that are consistent with NMR measurements. While a crystal structure-like conformation is observed in both solution and membrane simulations, a particular intermolecular rearrangement of the ternary complex is observed only when it is anchored to the membrane. This configuration emerges when the CRD hydrophobic loops are inserted into the membrane and helices α3-5 of KRAS4b are solvent exposed. This membrane-specific configuration is stabilized by KRAS4b-CRD contacts that are not observed in the crystal structure. These results suggest modulatory interplay between the CRD and plasma membrane that correlate with RAS/RAF complex structure and dynamics, and potentially influence subsequent steps in the activation of MAPK signaling.

Molecular Characterization and In Silico Analyses of Maurolipin Structure as a Secretory Phospholipase A 2 ( sPLA 2 ) from Venom Glands of Iranian Scorpio maurus (Arachnida: Scorpionida)

The venom is a mixture of various compounds with specific biological activities, such as the phospholipase  A₂ (PLA₂) enzyme present in scorpion venom. PLA₂ plays a key role in inhibiting ryanodine receptor channels and has neurotoxic activity. This study is the first investigation of molecular characterization, cloning, and in silico analyses of PLA₂ from Iranian Scorpio maurus, named Maurolipin. After RNA extraction from S. maurus venom glands, cDNA was synthesized and amplified through RT-PCR using specific primers. Amplified Maurolipin was cloned in TA cloning vector, pTG19. For in silico analyses, the characterized gene was analyzed utilizing different software. Maurolipin coding gene with 432 base pair nucleotide length encoded a protein of 144 amino acid residues and 16.34 kilodaltons. Comparing the coding sequence of Maurolipin with other characterized PLA₂ from different species of scorpions showed that this protein was a member of the PLA₂ superfamily. According to SWISS-MODEL prediction, Maurolipin had 38.83% identity with bee venom PLA₂ with 100% confidence and 39% identity with insect phospholipase A₂ family, which Phyre2 predicted. According to the three-dimensional structure prediction, Maurolipin with five disulfide bonds has a very high similarity to the structure of PLA₂ that belonged to the group III subfamily. The in silico analyses showed that phospholipase A₂ coding gene and protein structure is different based on scorpion species and geographical condition in which they live.

Glycosylated Flavonoid Compounds as Potent CYP121 Inhibitors of Mycobacterium tuberculosis

Due to the concerning rise in the number of multiple- and prolonged-drug-resistant (MDR and XDR) Mycobacterium tuberculosis (Mtb) strains, unprecedented demand has been created to design and develop novel therapeutic drugs with higher efficacy and safety. In this study, with a focused view on implementing an in silico drug design pipeline, a diverse set of glycosylated flavonoids were screened against the Mtb cytochrome-P450 enzyme 121 (CYP121), which is established as an approved drug target for the treatment of Mtb infection. A total of 148 glycosylated flavonoids were screened using structure-based virtual screening against the crystallized ligand, i.e., the L44 inhibitor, binding pocket in the Mtb CYP121 protein. Following this, only the top six compounds with the highest binding scores (kcal/mol) were considered for further intermolecular interaction and dynamic stability using 100 ns classical molecular dynamics simulation. These results suggested a considerable number of hydrogen and hydrophobic interactions and thermodynamic stability in comparison to the reference complex, i.e., the CYP121-L44 inhibitor. Furthermore, binding free energy via the MMGBSA method conducted on the last 10 ns interval of MD simulation trajectories revealed the substantial affinity of glycosylated compounds with Mtb CYP121 protein against reference complex. Notably, both the docked poses and residual energy decomposition via the MMGBSA method demonstrated the essential role of active residues in the interactions with glycosylated compounds by comparison with the reference complex. Collectively, this study demonstrates the viability of these screened glycosylated flavonoids as potential inhibitors of Mtb CYP121 for further experimental validation to develop a therapy for the treatment of drug-resistant Mtb strains.

Structural Characteristics of the 5′-Terminal Region of Mouse p53 mRNA and Identification of Proteins That Bind to This mRNA Region

A mouse model has often been used in studies of p53 gene expression. Detailed interpretation of functional studies is, however, hampered by insufficient knowledge of the impact of mouse p53 mRNA’s structure and its interactions with proteins in the translation process. In particular, the 5′-terminal region of mouse p53 mRNA is an important region which takes part in the regulation of the synthesis of p53 protein and its N-truncated isoform Δ41p53. In this work, the spatial folding of the 5′-terminal region of mouse p53 mRNA and its selected sub-fragments was proposed based on the results of the SAXS method and the RNAComposer program. Subsequently, RNA-assisted affinity chromatography was used to identify proteins present in mouse fibroblast cell lysates that are able to bind the RNA oligomer, which corresponds to the 5′-terminal region of mouse p53 mRNA. Possible sites to which the selected, identified proteins can bind were proposed. Interestingly, most of these binding sites coincide with the sites determined as accessible to hybridization of complementary oligonucleotides. Finally, the high binding affinity of hnRNP K and PCBP2 to the 5′-terminal region of mouse p53 mRNA was confirmed and their possible binding sites were proposed.

Impact of Antibiotics as Waste, Physical, Chemical, and Enzymatical Degradation: Use of Laccases

The first traces of Tetracycline (TE) were detected in human skeletons from Sudan and Egypt, finding that it may be related to the diet of the time, the use of some dyes, and the use of soils loaded with microorganisms, such as Streptomyces spp., among other microorganisms capable of producing antibiotics. However, most people only recognise authors dating between 1904 and 1940, such as Ehrlich, Domagk, and Fleming. Antibiotics are the therapeutic option for countless infections treatment; unfortunately, they are the second most common group of drugs in wastewaters worldwide due to failures in industrial waste treatments (pharmaceutics, hospitals, senior residences) and their irrational use in humans and animals. The main antibiotics problem lies in delivered and non-prescribed human use, use in livestock as growth promoters, and crop cultivation as biocides (regulated activities that have not complied in some places). This practice has led to the toxicity of the environment as antibiotics generate eutrophication, water pollution, nutrient imbalance, and press antibiotic resistance. In addition, the removal of antibiotics is not a required process in global wastewater treatment standards. This review aims to raise awareness of the negative impact of antibiotics as residues and physical, chemical, and biological treatments for their degradation. We discuss the high cost of physical and chemical treatments, the risk of using chemicals that worsen the situation, and the fact that each antibiotic class can be transformed differently with each of these treatments and generate new compounds that could be more toxic than the original ones; also, we discuss the use of enzymes for antibiotic degradation, with emphasis on laccases.

Exploration of chalcones as 3-chymotrypsin-like protease (3CLpro) inhibitors of SARS-CoV-2 using computational approaches

The main protease 3CL^pro is one of the potential targets against coronavirus. Inhibiting this enzyme leads to the interruption of viral replication. Chalcone and its derivatives were reported to possess the ability to bind to 3CL^pro protease in the binding pocket. This study explored an in-house database of 269 chalcones as 3CL^pro inhibitors using in silico screening models, including molecular docking, molecular dynamics simulation, binding free energy calculation, and ADME prediction. C264 and C235 stand out as the two most potential structures. The top hit compound C264 was with the Jamda score of −2.8329 and the MM/GBSA binding energy mean value of −28.23 ± 3.53 kcal/mol, which was lower than the reference ligand. Despite the lower mean binding energy (−22.07 ± 3.39 kcal/mol), in-depth analysis of binding interaction suggested C235 could be another potential candidate. Further, in vitro and in vivo experiments are required to confirm the inhibitory ability.

Computational Reconstruction and Analysis of Structural Models of Avian Cryptochrome 4

A recent study by Xu et al. (Nature, 2021, 594, 535-540) provided strong evidence that cryptochrome 4 (Cry4) is a key protein to endow migratory birds with the magnetic compass sense. The investigation compared the magnetic field response of Cry4 from migratory and nonmigratory bird species and suggested that a difference in magnetic sensitivity could exist. This finding prompted an in-depth investigation into Cry4 protein differences on the structural and dynamic levels. In the present study, the pigeon Cry4 (ClCry4) crystal structure was used to reconstruct the missing avian Cry4 protein structures via homology modeling for carefully selected bird species. The reconstructed Cry4 structure from European robin, Eurasian blackcap, zebra finch, chicken, and pigeon were subsequently simulated dynamically and analyzed. The studied avian Cry4 structures show flexibility in analogous regions pointing to similar activation mechanisms and/or signaling interaction partners. It can be concluded that the experimentally recorded difference in the magnetic field sensitivity of Cry4 from different birds is unlikely to be due to solely intrinsic dynamics of the proteins but requires additional factors that have not yet been identified.

Seeking antiviral drugs to inhibit SARS-CoV-2 RNA dependent RNA polymerase: A molecular docking analysis

COVID-19 outbreak associated with the severe acute respiratory syndrome coronavirus (SARS-CoV-2) raised health concerns across the globe and has been considered highly transmissible between people. In attempts for finding therapeutic treatment for the new disease, this work has focused on examining the polymerase inhibitors against the SARS-CoV-2 nsp12 and co-factors nsp8 and nsp7. Several polymerase inhibitors were examined against PDB ID: 6M71 using computational analysis evaluating the ligand's binding affinity to replicating groove to the active site. The findings of this analysis showed Cytarabine of -5.65 Kcal/mol with the highest binding probability (70%) to replicating groove of 6M71. The complex stability was then examined over 19 ns molecular dynamics simulation suggesting that Cytarabine might be possible potent inhibitor for the SARS-CoV-2 RNA Dependent RNA Polymerase.

Evidence that nuclear receptors are related to terpene synthases

Ligand-activated nuclear receptors (NRs) orchestrate development, growth, and reproduction across all animal lifeforms - the Metazoa - but how NRs evolved remains mysterious. Given the NR ligands including steroids and retinoids are predominantly terpenoids, we asked whether NRs might have evolved from enzymes that catalyze terpene synthesis and metabolism. We provide evidence suggesting that NRs may be related to the terpene synthase (TS) enzyme superfamily. Based on over 10,000 3D structural comparisons, we report that the NR ligand-binding domain and TS enzymes share a conserved core of seven α-helical segments. In addition, the 3D locations of the major ligand-contacting residues are also conserved between the two protein classes. Primary sequence comparisons reveal suggestive similarities specifically between NRs and the subfamily of cis-isoprene transferases, notably with dehydrodolichyl pyrophosphate synthase and its obligate partner, NUS1/NOGOB receptor. Pharmacological overlaps between NRs and TS enzymes add weight to the contention that they share a distant evolutionary origin, and the combined data raise the possibility that a ligand-gated receptor may have arisen from an enzyme antecedent. However, our findings do not formally exclude other interpretations such as convergent evolution, and further analysis will be necessary to confirm the inferred relationship between the two protein classes.

Identification of potential interleukin-8 inhibitors acting on the interactive site between chemokine and CXCR2 receptor: A computational approach

Interactions between interleukin (IL)-8 and its receptors, CXCR1, and CXCR2, serve crucial roles in inflammatory conditions and various types of cancers. Inhibition of this signaling pathway has been exploited as a promising strategy in treating these diseases. However, most studies only focused on the design of allosteric antagonists-bound receptors on the intracellular side of IL-8 receptors. Recently, the first cryo-EM structures of IL-8-CXCR2-Gi complexes have been solved, revealing the unique binding and activation modes of the endogenous chemokine IL-8. Hence, we set to identify small molecule inhibitors for IL-8 using critical protein-protein interaction between IL-8 and CXCR2 at the orthosteric binding site. The pharmacophore models and molecular docking screened compounds from DrugBank and NCI databases. The oral bioavailability of the top 23 ligands from the screening was then predicted by the SwissAMDE tool. Molecular dynamics simulation and free binding energy calculation were performed for the best compounds. The result indicated that DB14770, DB12121, and DB03916 could form strong interactions and stable protein-ligand complexes with IL-8. These three candidates are potential IL-8 inhibitors that can be further evaluated by in vitro experiments in the next stage.

Atomistic mechanism of leptin and leptin-receptor association

The leptin-leptin receptor complex is at the very core of energy homeostasis and immune system regulation, among many other functions. In this work, we built homology models of leptin and the leptin binding domain (LBD) of the receptor from humans and mice. Docking analyses were used to obtain the coordinates of the native leptin-LBD complexes and a mixed heterodimer formed by human leptin and mouse LBD. Molecular dynamics (MD) simulations were performed using all models (monomers and heterodimers) as initial coordinates and the GROMACS program. The overall structural and dynamical behaviors are similar for the three complexes. Upon MD simulations, several new interactions appear. In particular, hydrophobic interactions, with more than 90% persistence, seem to be the most relevant for the stability of the dimers, as well as the pair formed by Asp85_Lep and Arg468_LBD. This in silico analysis provides structural and dynamical information, at the atomistic level, about the mechanism of leptin-LBD complex formation and leptin receptor activation. This knowledge might be used in the rational drug design of therapeutics to modulate leptin signaling.Communicated by Ramaswamy H. Sarma.

Modeling SARS-CoV-2 spike/ACE2 protein-protein interactions for predicting the binding affinity of new spike variants for ACE2, and novel ACE2 structurally related human protein targets, for COVID-19 handling in the 3PM context

Aims

The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM.

Methods

In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein–protein interface.

Results

Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike “boat-shaped” receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the “stern” or “bow” regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the “hull” region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities.

Conclusion

In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-021-00267-w.

SARS-CoV-2 variant surge and vaccine breakthrough infection: A computational analysis

Coronavirus Delta variant was first detected in India in October of 2020, and it led to a massive second wave of COVID-19 cases in the country. Since then, the highly infectious Delta strain has been spreading globally. The Delta variant and its sub-lineages showed an increased infection rate with a reduced effect of the potential antibody neutralization. The current work is a modeled computational analysis of the mutated receptor-binding domain (RBD) of the SARS-CoV-2 B.1.617 lineage binding with ACE2 and GRP78 to understand the increased strain transmissibility. The cell-surface Glucose Regulated Protein 78 (GRP78) attached to the mutated ACE2-SARS-CoV-2 Spike RBD complex is modeled. The results showed that GRP78 β-substrate-binding domain weakly binds to the wild-type RBD combined with angiotensin-converting enzyme 2 (ACE2) within the SARS-CoV-2 Spike RBD-ACE2 complex. Both GRP78 and ACE2 bind approximately in the same region on the wild-type SARS-CoV-2 Spike RBD surface. On the other hand, GRP78 strongly binds to the mutated SARS-CoV-2 Spike RBD in the RBD-ACE2 complex through the α-substrate-binding domain instead of β-substrate-binding domain in a different region from that of ACE2. The current findings suggest that blocking the main ACE2 pathway may not prevent the interactions between GRP78 and the mutated SARS-CoV-2 Spike RBD, which might introduce an additional avenue into the virus invasion for the host cell if the ACE2 pathway is blocked by the neutralized antibodies. Hence, the peptide satpdb10668 has been proposed as a potential inhibitor of SARS-CoV-2 attachment and virus invasion into the host cell.

In silico assessment of chronic toxicity of a combination drug namely ‘Olmesartan medoxomil and Hydrochlorothiazide’, marketed in Bangladesh

Background

Nowadays combination therapy has become more popular due to their additional effect, synergistic effect and antagonistic effect. Any of these can influence the treatment profile. Combination therapy is used to treat some chronic diseases like diabetes, hypertension, cancer etc. But recently India has banned some fixed dose drug combinations due to their increased chances of adverse drug effects and drug interactions. So it is the time to take a look on the present drug combinations available in Bangladesh. An in silico study may provide important information about their probable toxicities. Drugs available in the combination may deposit slowly in the body and may lead to toxicities. Here an antihypertensive drug combination ‘Olmesartan medoxomil and Hydrochlorothiazide’ had been studied.

Results

Olmesartan medoxomil and Hydrochlorothiazide have not been found to comply any similar protein to interact with each other, thus no possible chance of additional toxicity of the combination in case of long term use.

Conclusions

At first, using PubChem the ligand was searched for a canonical SMILE. By inputting the canonical SMILE in Protox, a basic information about toxicities was predicted. From Swiss Target Prediction, target proteins responsible for both efficacy and toxicity were identified. These protein structures were downloaded from Protein Data Bank and edited with Flare. Undesired amino acid, ligand–ligand complex, fatty acid, and water molecules were removed by PyMOL. Structurally modified proteins and ligands were inputted in Swiss PDB viewer for energy minimization. Energy minimization is a very important step because unfavorable bond length, bond strength and torsion angle between protein and ligand may interfere with docking procedure. Then docking between Olmesartan medoxomil (ligand) and the proteins responsible for efficacy and toxicity was performed by PyRx. Vina binding affinity provided the value of binding strength between the ligand and the proteins, which determines how strong the bond is. The more negative the vina binding affinity, the stronger the bond. Discovery studio software was used to visualize the docking complexes. Same steps were followed for Hydrochlorothiazide to identify proteins responsible for desired and undesired effects, but no toxic effect was found from protox.

Molecular dynamic study of SARS-CoV-2 with various S protein mutations and their effect on thermodynamic properties

Studying the structure of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) protein is important to understand the infection process. The S protein is necessary in completing the virus life cycle and is responsible for the appearance of new variants and drug and vaccine resistance. Understanding the structure and dynamics of biological macromolecules is essential for understanding how they function. In this work, we investigated the effects of mutations on S protein stability and solubility through molecular dynamic (MD) simulation in a 100 ns (ns) period. We screened four variants in addition to the wild type (WT). Results show that changes on MD simulation parameters of S protein indicate fluctuations and changes in the conformation, especially in the area between 300 and 600 amino acids (aa). This provides us an image of how the virus protein can reshape itself to adapt to any changes that occur in human angiotensin-converting enzyme 2 or drugs that can target the protein region. Our results also show that the Brazil variant has high fluctuations and unstable folding at some stages compared with other variants.

In silico molecular docking analysis for repurposing approved antiviral drugs against SARS-CoV-2 main protease

Developing a safe and effective antiviral treatment takes a decade, however, when it comes to the coronavirus disease (COVID-19), time is a sensitive matter to slow the spread of the pandemic. Screening approved antiviral drugs against COVID-19 would speed the process of finding therapeutic treatment. The current study examines commercially approved drugs to repurpose them against COVID-19 virus main protease using structure-based in-silico screening. The main protease of the coronavirus is essential in the viral replication and is involved in polyprotein cleavage and immune regulation, making it an effective target when developing the treatment. A Number of approved antiviral drugs were tested against COVID-19 virus using molecular docking analysis by calculating the free natural affinity of the binding ligand to the active site pocket and the catalytic residues without forcing the docking of the ligand to active site. COVID-19 virus protease solved structure (PDB ID: 6LU7) is targeted by repurposed drugs. The molecular docking analysis results have shown that the binding of Remdesivir and Mycophenolic acid acyl glucuronide with the protein drug target has optimal binding features supporting that Remdesivir and Mycophenolic acid acyl glucuronide can be used as potential anti-viral treatment against COVID-19 disease.

Fluorescence-based screening for engineered aldo-keto reductase KmAKR with improved catalytic performance and extended substrate scope

BACKGROUND

Aldo-keto reductases-catalyzed transformations of ketones to chiral alcohols have become an established biocatalytic process step in the pharmaceutical industry. Previously, we have discovered an aldo-keto reductase (AKR) from Kluyveromyces marxianus that is active to the aliphatic tert-butyl 6-substituted (5R/S)-hydroxy-3-oxohexanoates, but it is inactive to aromatic ketones. In order to acquire an excellent KmAKRmutant for ensuring the simultaneous improvement of activity-thermostability toward tert-butyl 6-cyano-(5R)-hydroxy-3-oxohexanoate ((5R)-1) and broadening the universal application prospects toward more substrates covering both aliphatic and aromatic ketones, a fluorescence-based high-throughput (HT) screening technique was established.

MAIN METHODS AND MAJOR RESULTS

The directed evolution of KmAKR variant M5 (KmAKR-W297H/Y296W/K29H/Y28A/T63M) produced the "best" variant M5-Q213A/T23V. It exhibited enhanced activity-thermostability toward (5R)-1, improved activity toward all 18 test substrates and strict R-stereoselectivity toward 10 substrates in comparison to M5. The enhancement of enzymatic activity and the extension of substrate scope covering aromatic ketones are proposed to be largely attributed to pushing the binding pocket of M5-Q213A/T23V to the enzyme surface, decreasing the steric hindrance at the entrance and enhancing the flexibility of loops surrounding the active center. In addition, combined with 0.94 g dry cell weight (DCW) L^-1 glucose dehydrogenase from Exiguobacterium sibiricum (EsGDH) for NADPH regeneration, 2.81 g DCW L^-1 M5-Q213A/T23V completely converted (5R)-1 of up to 450 g L^-1 at 120 g g^-1 substrates/catalysts (S/C), yielding the corresponding optically pure tert-butyl 6-cyano-(3R,5R)-dihydroxyhexanoate ((3R,5R)-2, > 99.5% d.e._p ) with a space-time yield (STY) of 1.08 kg L^-1 day^-1 .

CONCLUSIONS

A fluorescence-based HT screening system was developed to tailor KmAKR's activity, thermostability and substrate scope. The "best" variant M5-Q213A/T23V holds great potential application for the synthesis of aliphatic/aromatic R-configuration alcohols.