PDBsum new things

Insights into docking in megasynthases from the investigation of the toblerol trans-AT polyketide synthase: many α-helical means to an end

The fidelity of biosynthesis by modular polyketide synthases (PKSs) depends on specific moderate affinity interactions between successive polypeptide subunits mediated by docking domains (DDs). These sequence elements are notably portable, allowing their transplantation into alternative biosynthetic and metabolic contexts. Herein, we use integrative structural biology to characterize a pair of DDs from the toblerol trans-AT PKS. Both are intrinsically disordered regions (IDRs) that fold into a 3 α-helix docking complex of unprecedented topology. The C-terminal docking domain (CDD) resembles the 4 α-helix type (4HB) CDDs, which shows that the same type of DD can be redeployed to form complexes of distinct geometry. By carefully re-examining known DD structures, we further extend this observation to type 2 docking domains, establishing previously unsuspected structural relations between DD types. Taken together, these data illustrate the plasticity of α-helical DDs, which allow the formation of a diverse topological spectrum of docked complexes. The newly identified DDs should also find utility in modular PKS genetic engineering.

Identification of novel T-cell epitopes on monkeypox virus and development of multi-epitopes vaccine using immunoinformatics approaches

Monkeypox virus (MPV) is closely related to the smallpox virus, and previous data from Africa suggest that the smallpox vaccine (VARV) is at least 85% effective in preventing MPV. No multi-epitope vaccine has yet been developed to prevent MPV infection. In this work, we used in silico structural biology and advanced immunoinformatic strategies to design a multi-epitope subunit vaccine against MPV infection. The designed vaccine sequence is adjuvanted with CpG-ODN and includes HTL/CTL epitopes for similar proteins between vaccinia virus (VACV) that induced T-cell production in vaccinated volunteers and the first draft sequence of the MPV genome associated with the suspected outbreak in several countries, May 2022. In addition, the specific binding of the modified vaccine and the immune Toll-like receptor 9 (TLR9) was estimated by molecular interaction studies. Strong interaction in the binding groove as well as good docking scores confirmed the stringency of the modified vaccine. The stability of the interaction was confirmed by a classical molecular dynamics simulation and normal mode analysis. Then, the immune simulation also indicated the ability of this vaccine to induce an effective immune response against MPV. Codon optimization and in silico cloning of the vaccine into the pET-28a (+) vector also showed its expression potential in the E. coli K12 system. The promising data obtained from the various in silico studies indicate that this vaccine is effective against MPV. However, additional in vitro and in vivo studies are still needed to confirm its efficacy.Communicated by Ramaswamy H. Sarma.

Subtractive proteomics-based vaccine targets annotation and reverse vaccinology approaches to identify multiepitope vaccine against Plesiomonas shigelloides

Plesiomonas shigelloides, an aquatic bacterium belonging to the Enterobacteriaceae family, is a frequent cause of gastroenteritis with diarrhea and gastrointestinal severe disease. Despite decades of research, discovering a licensed and globally accessible vaccine is still years away. Developing a putative vaccine that can combat the Plesiomonas shigelloides infection by boosting population immunity against P. shigelloides is direly needed. In the framework of the current study, the entire proteome of P. shigelloides was explored using subtractive genomics integrated with the immunoinformatics approach for designing an effective vaccine construct against P. shigelloides. The overall stability of the vaccine construct was evaluated using molecular docking, which demonstrated that MEV showed higher binding affinities with toll-like receptors (TLR4: 51.5 ± 10.3, TLR2: 60.5 ± 9.2) and MHC receptors(MHCI: 79.7 ± 11.2 kcal/mol, MHCII: 70.4 ± 23.7). Further, the therapeutic efficacy of the vaccine construct for generating an efficient immune response was evaluated by computational immunological simulation. Finally, computer-based cloning and improvement in codon composition without altering amino acid sequence led to the development of a proposed vaccine. In a nutshell, the findings of this study add to the existing knowledge about the pathogenesis of this infection. The schemed MEV can be a possible prophylactic agent for individuals infected with P. shigelloides. Nevertheless, further authentication is required to guarantee its safeness and immunogenic potential.

Design and Synthesis of 2-Phenylindolizine Acetamides: Molecular Docking, in Vitro Antimicrobial and Anticancer Activity Evaluation

In the present work, we synthesized a small library of 2-phenylindolizine acetamide derivatives 7a-i and studied their biological activity. The synthesis was accomplished starting with easily available starting material phenacyl bromide 1 proceeding through the key intermediate 6-methyl-7-nitro-2-phenylindolizine 4. All the compounds 7a-i were characterized using spectroscopy viz., 1H-NMR, 13C NMR, FTIR, and mass spectrometry. Interestingly, 2-phenylindolizine scaffolds 7c, 7f and 7g revealed a remarkable antibacterial activity against relevant organisms S. aureus, E. coli, S. pneumoniae, P. aeruginosa. The target compounds 7e and 7h showed excellent anticancer activity against Colo-205 and MDA-MB-231 cell lines with IC50 values of 68.62, 62.91, 54.23 and 46.34 μM respectively. Additionally, all the 2-phenylindolizine acetamide derivatives 7a-i were subjected to molecular docking prediction by Autodock 4.2. Compounds 7a, 7f and 7c exhibited very good hydrogen bonding amino acid interactions Asp83 (2.23 Å), Asp83 (2.08 Å), His74 (2.05 Å), His76 (1.71 Å), Ser80 (1.05 Å) with active site of Topoisomerase-IV from S. pneumoniae (4KPE). Further, the compounds 7a-i have revealed acceptable ranges for drug-likeliness properties upon evaluation using SwissADME for ADMET and physiochemical properties.

Cryo-EM analysis of Pseudomonas phage Pa193 structural components

The World Health Organization has designated Pseudomonas aeruginosa as a critical pathogen for the development of new antimicrobials. Bacterial viruses, or bacteriophages, have been used in various clinical settings, commonly called phage therapy, to address this growing public health crisis. Here, we describe a high-resolution structural atlas of a therapeutic, contractile-tailed Pseudomonas phage, Pa193. We used bioinformatics, proteomics, and cryogenic electron microscopy single particle analysis to identify, annotate, and build atomic models for 21 distinct structural polypeptide chains forming the icosahedral capsid, neck, contractile tail, and baseplate. We identified a putative scaffolding protein stabilizing the interior of the capsid 5-fold vertex. We also visualized a large portion of Pa193 ~ 500 Å long tail fibers and resolved the interface between the baseplate and tail fibers. The work presented here provides a framework to support a better understanding of phages as biomedicines for phage therapy and inform engineering opportunities.

Structure of glyceraldehyde-3-phosphate dehydrogenase from Paracoccidioides lutzii in complex with an aldonic sugar acid

The pathogen Paracoccidioides lutzii (Pb01) is found in South America countries Colombia, Ecuador, Venezuela and Brazil, especially in the central, west, and north regions of the latter. It belongs to the Ajellomycetaceae family, Onygenales order, and is typically thermodimorphic, presenting yeast cells when it grows in animal tissues, but mycelia when in the environment, where it produces the infectious propagule. This fungus is one of the etiologic agents of Paracoccidioidomycosis (PCM), the most important endemic fungal infection in Latin America. Investigations on its genome have contributed to a better understanding about its metabolism and revealed the complexity of several metabolic glycolytic pathways. Glyceraldehyde-3-Phosphate Dehydrogenase from Paracoccidioides lutzii (PlGAPDH) is considered a moonlighting protein and participates in several biological processes of this pathogen. The enzyme was expressed and purified, as seen in SDS-PAGE gel, crystallized and had its three dimensional structure (3D) determined in complex with NAD+, a sulphate ion and d-galactonic acid, therefore, a type of 'GAA site'. It is the first GAPDH structure to show this chemical type in this site and how this protein can bind an acid derived from oxidation of a linear hexose.

Proteomics-based vaccine targets annotation and design of multi-epitope vaccine against antibiotic-resistant Streptococcus gallolyticus

Streptococcus gallolyticus is a non-motile, gram-positive bacterium that causes infective endocarditis. S. gallolyticus has developed resistance to existing antibiotics, and no vaccine is currently available. Therefore, it is essential to develop an effective S. gallolyticus vaccine. Core proteomics was used in this study together with subtractive proteomics and reverse vaccinology approach to find antigenic proteins that could be utilized for the design of the S. gallolyticus multi-epitope vaccine. The pipeline identified two antigenic proteins as potential vaccine targets: penicillin-binding protein and the ATP synthase subunit. T and B cell epitopes from the specific proteins were forecasted employing several immunoinformatics and bioinformatics resources. A vaccine (360 amino acids) was created using a combination of seven cytotoxic T cell lymphocyte (CTL), three helper T cell lymphocyte (HTL), and five linear B cell lymphocyte (LBL) epitopes. To increase immune responses, the vaccine was paired with a cholera enterotoxin subunit B (CTB) adjuvant. The developed vaccine was highly antigenic, non-allergenic, and stable for human use. The vaccine's binding affinity and molecular interactions with the human immunological receptor TLR4 were studied using molecular mechanics/generalized Born surface area (MMGBSA), molecular docking, and molecular dynamic (MD) simulation analyses. Escherichia coli (strain K12) plasmid vector pET-28a ( +) was used to examine the ability of the vaccine to be expressed. According to the outcomes of these computer experiments, the vaccine is quite promising in terms of developing a protective immunity against diseases. However, in vitro and animal research are required to validate our findings.

Revolutionizing and identifying novel drug targets in Citrobacter koseri via subtractive proteomics and development of a multi-epitope vaccine using reverse vaccinology and immuno-informatics

Citrobacter koseri is a gram-negative rod that has been linked to infections in people with significant comorbidities and immunocompromised immune systems. It is most commonly known to cause urinary tract infections. Thus, the development of an efficacious C. koseri vaccine is imperative, as the pathogen has acquired resistance to current antibiotics. Subtractive proteomics was employed during this research to identify potential antigenic proteins to design an effective vaccine against C. koseri. The pipeline identified two antigenic proteins as potential vaccine targets: DP-3-O-acyl-N-acetylglucosamine deacetylase and Arabinose 5-phosphate isomerase. B and T cell epitopes from the specific proteins were forecasted employing several immunoinformatic and bioinformatics resources. A vaccine was created using a combination of seven cytotoxic T cell lymphocytes (CTL), five helper T cell lymphocyte (HTL), and seven linear B cell lymphocyte (LBL) epitopes. An adjuvant (β-defensin) was added to the vaccine to enhance immunological responses. The created vaccine was stable for use in humans, highly antigenic, and non-allergenic. The vaccine's molecular and interactions binding affinity with the human immunological receptor TLR3 were studied using MMGBSA, molecular dynamics (MD) simulations, and molecular docking analyses. E. coli (strain-K12) plasmid vector pET-28a (+) was used to examine the ability of the vaccine to be expressed. The vaccine shows great promise in terms of developing protective immunity against diseases, based on the results of these computer experiments. However, in vitro and animal research are required to validate our findings.Communicated by Ramaswamy H. Sarma.

Integrative structural analysis of Pseudomonas phage DEV reveals a genome ejection motor

DEV is an obligatory lytic Pseudomonas phage of the N4-like genus, recently reclassified as Schitoviridae. The DEV genome encodes 91 ORFs, including a 3,398 amino acid virion-associated RNA polymerase. Here, we describe the complete architecture of DEV, determined using a combination of cryo-electron microscopy localized reconstruction, biochemical methods, and genetic knockouts. We built de novo structures of all capsid factors and tail components involved in host attachment. We demonstrate that DEV long tail fibers are essential for infection of Pseudomonas aeruginosa and dispensable for infecting mutants with a truncated lipopolysaccharide devoid of the O-antigen. We identified DEV ejection proteins and, unexpectedly, found that the giant DEV RNA polymerase, the hallmark of the Schitoviridae family, is an ejection protein. We propose that DEV ejection proteins form a genome ejection motor across the host cell envelope and that these structural principles are conserved in all Schitoviridae.

Evaluation of mutations on O6-methylguanine methyl transferase structure and its interactions: molecular dynamics simulation study

O6-methylguanine DNA methyl transferase (MGMT) is a significant vehicle for the cellular clearance of alkyl lesions, particularly the methyl group of the O-6 and O-4 positions of guanine and thymine, respectively. Many publications have studied the correlation between polymorphisms in MGMT and susceptibility to various cancers. In the present study, we investigated the consequence of L84F, common single-nucleotide polymorphism, K125E, site-specific mutagenesis, and L84F/K125E on conformation, stability, and behavior of MGMT in the free form and interaction with proliferating cell nuclear antigen (PCNA) and DNA as partners in the biochemical network by using molecular dynamics simulation method. Our results showed that all free variants of MGMT differed from the native form. However, among all free variants of MGMT, the L84F/K125E variant exhibited similar properties compared with the wild-type. In contrast, in complex modes, only amino acid residues of the L84F variant are involved in the interactions with PCNA and DNA somewhat differently relative to the wild-type. Furthermore, L84F SNP showed the highest binding free energy compared to other variants and native forms. These alterations in the amino acids and binding free energy of L84F relative to the native are the reasons for changing its region connection compared to the native form. Therefore, we propose conducting further investigations into the impact of inhibitors or chemotherapeutic agents to assess their effectiveness on MGMT variants compared to the wild-type, aiming to reduce the cost of cancer treatment that will depend on inhibiting native MGMT protein.Communicated by Ramaswamy H. Sarma.

Uncovering the structural stability of Magnaporthe oryzae effectors: a secretome-wide in silico analysis

Rice blast, caused by the ascomycete fungus Magnaporthe oryzae, is a deadly disease and a major threat to global food security. The pathogen secretes small proteinaceous effectors, virulence factors, inside the host to manipulate and perturb the host immune system, allowing the pathogen to colonize and establish a successful infection. While the molecular functions of several effectors are characterized, very little is known about the structural stability of these effectors. We analyzed a total of 554 small secretory proteins (SSPs) from the M. oryzae secretome to decipher key features of intrinsic disorder (ID) and the structural dynamics of the selected putative effectors through thorough and systematic in silico studies. Our results suggest that out of the total SSPs, 66% were predicted as effector proteins, released either into the apoplast or cytoplasm of the host cell. Of these, 68% were found to be intrinsically disordered effector proteins (IDEPs). Among the six distinct classes of disordered effectors, we observed peculiar relationships between the localization of several effectors in the apoplast or cytoplasm and the degree of disorder. We determined the degree of structural disorder and its impact on protein foldability across all the putative small secretory effector proteins from the blast pathogen, further validated by molecular dynamics simulation studies. This study provides definite clues toward unraveling the mystery behind the importance of structural distortions in effectors and their impact on plant-pathogen interactions. The study of these dynamical segments may help identify new effectors as well.Communicated by Ramaswamy H. Sarma.

Molecular Architecture of Salmonella Typhimurium Virus P22 Genome Ejection Machinery

Bacteriophage P22 is a prototypical member of the Podoviridae superfamily. Since its discovery in 1952, P22 has become a paradigm for phage transduction and a model for icosahedral viral capsid assembly. Here, we describe the complete architecture of the P22 tail apparatus (gp1, gp4, gp10, gp9, and gp26) and the potential location and organization of P22 ejection proteins (gp7, gp20, and gp16), determined using cryo-EM localized reconstruction, genetic knockouts, and biochemical analysis. We found that the tail apparatus exists in two equivalent conformations, rotated by ∼6° relative to the capsid. Portal protomers make unique contacts with coat subunits in both conformations, explaining the 12:5 symmetry mismatch. The tail assembles around the hexameric tail hub (gp10), which folds into an interrupted β-propeller characterized by an apical insertion domain. The tail hub connects proximally to the dodecameric portal protein and head-to-tail adapter (gp4), distally to the trimeric tail needle (gp26), and laterally to six trimeric tailspikes (gp9) that attach asymmetrically to gp10 insertion domain. Cryo-EM analysis of P22 mutants lacking the ejection proteins gp7 or gp20 and biochemical analysis of purified recombinant proteins suggest that gp7 and gp20 form a molecular complex associated with the tail apparatus via the portal protein barrel. We identified a putative signal transduction pathway from the tailspike to the tail needle, mediated by three flexible loops in the tail hub, that explains how lipopolysaccharide (LPS) is sufficient to trigger the ejection of the P22 DNA in vitro.

The DNA-binding induced (de)AMPylation activity of a Coxiella burnetii Fic enzyme targets Histone H3

The intracellular bacterial pathogen Coxiella burnetii evades the host response by secreting effector proteins that aid in establishing a replication-friendly niche. Bacterial filamentation induced by cyclic AMP (Fic) enzymes can act as effectors by covalently modifying target proteins with the posttranslational AMPylation by transferring adenosine monophosphate (AMP) from adenosine triphosphate (ATP) to a hydroxyl-containing side chain. Here we identify the gene product of C. burnetii CBU_0822, termed C. burnetii Fic 2 (CbFic2), to AMPylate host cell histone H3 at serine 10 and serine 28. We show that CbFic2 acts as a bifunctional enzyme, both capable of AMPylation as well as deAMPylation, and is regulated by the binding of DNA via a C-terminal helix-turn-helix domain. We propose that CbFic2 performs AMPylation in its monomeric state, switching to a deAMPylating dimer upon DNA binding. This study unveils reversible histone modification by a specific enzyme of a pathogenic bacterium.

Computational Modelling, Functional Characterization and Molecular Docking to Lead Compounds of Bordetella pertussis Diaminopimelate Epimerase

Bordetella pertussis, the causative agent of whooping cough, is an opportunistic virulent bacterial pathogen that is resistant to a wide range of antibiotics due to a variety of resistance mechanisms. Looking at the increasing number of infections caused by B. pertussis and its resistance to diverse antibiotics, it is essential to develop alternative strategies to fight against B. pertussis. Diaminopimelate epimerase (DapF) is an important enzyme of the lysine biosynthesis pathway in B. pertussis that catalyzes the formation of meso-2, 6-diaminoheptanedioate (meso-DAP), which is an important step in lysine metabolism. Therefore, Bordetella pertussis diaminopimelate epimerase (DapF) becomes an ideal target for antimicrobial drug development. In the present study, computational modelling, functional characterization, binding studies, and docking studies of BpDapF with lead compounds were carried out using different in silico tools. In silico prediction results in the secondary structure, 3-D structure analysis, and protein-protein interaction analysis of BpDapF. Docking studies further showed the respective amino acid residues for ligands in the phosphate‑binding loop of BpDapF play a vital role in the formation of H‑bonds with these ligands. The site where the ligand was bound is a deep groove, which is regarded as the binding cavity of the protein. Biochemical studies indicated that Limonin (binding energy - 8.8 kcal/mol), Ajmalicine (binding energy - 8.7 kcal/mol), Clinafloxacin (binding energy - 8.3 kcal/mol), Dexamethasone (binding energy - 8.2 kcal/mol), and Tetracycline (binding energy - 8.1 kcal/mol) exhibited promising binding towards the drug target DapF of B. pertussis in comparison with the binding between other drugs and act as the potential inhibitors of BpDapF that eventually can reduce the catalytic activity of BpDapF.

Design and characterization of a multi-epitope vaccine targeting Chlamydia abortus using immunoinformatics approach

Chlamydiosis is a widespread ailment affecting humans, livestock, and wildlife, caused by C. abortus, a member of the Chlamydia genus. This disease leads to reproductive disorders in bovines and poses a zoonotic risk, resulting in adverse outcomes such as abortion, stillbirths, weak offspring, endometritis, repeat breeding, and perinatal mortality. However, current chlamydiosis vaccines have limitations in terms of safety, efficacy, and stability, necessitating the development of effective and safe alternatives. In this study, our objective was to design a multi-epitope vaccine (MEV) targeting all strains of C. abortus using bioinformatics and immunoinformatics approaches. We identified highly antigenic and non-allergic proteins (yidC, yajC, secY, CAB503, and CAB746) using VaxiJen and AlgPred tools. Physicochemical analyses and secondary structure predictions confirmed protein stability through ProtParam and SOPMA methods. Furthermore, we employed IEDB-AR, NETMHCpan, and ToxinPred2 tools to predict cytotoxic T lymphocyte (CTL), helper T lymphocyte (HTL), and B-cell epitopes, resulting in the identification of conserved epitopes for further analysis. The MEV construct, consisting of 545 amino acids, incorporated the adjuvant Beta defensin-3, along with 9 CTL epitopes and 21 HTL epitopes linked by EAAAK, KK, and AAY linkers. We assessed the safety and immunogenicity of the vaccine through comprehensive evaluations of antigenicity, toxicity, allergenicity, and physicochemical properties. Structural stability and quality were examined using 3D modeling via the ab initio approach with the Robetta platform. Molecular docking analysis explored the compatibility of the MEV with Toll-like receptor 9 (TLR9) using ClusPro, while molecular dynamics simulation with the DESMOND Maestro software predicted the stability and flexibility of the docked complex. Despite promising in silico findings, further wet lab investigations are crucial to validate the safety and efficacy of the MEV. Successful development and validation of this MEV hold significant potential in combatting chlamydiosis in both animal and human populations.Communicated by Ramaswamy H. Sarma.

Structural and functional analyses of Burkholderia pseudomallei BPSL1038 reveal a Cas-2/VapD nuclease sub-family

Burkholderia pseudomallei is a highly versatile pathogen with ~25% of its genome annotated to encode hypothetical proteins. One such hypothetical protein, BPSL1038, is conserved across seven bacterial genera and 654 Burkholderia spp. Here, we present a 1.55 Å resolution crystal structure of BPSL1038. The overall structure folded into a modified βαββαβα ferredoxin fold similar to known Cas2 nucleases. The Cas2 equivalent catalytic aspartate (D11) pairs are conserved in BPSL1038 although B. pseudomallei has no known CRISPR associated system. Functional analysis revealed that BPSL1038 is a nuclease with endonuclease activity towards double-stranded DNA. The DNase activity is divalent ion independent and optimum at pH 6. The concentration of monovalent ions (Na+ and K+) is crucial for nuclease activity. An active site with a unique D11(X20)SST motif was identified and proposed for BPSL1038 and its orthologs. Structure modelling indicates the catalytic role of the D11(X20)SST motif and that the arginine residues R10 and R30 may interact with the nucleic acid backbone. The structural similarity of BPSL1038 to Cas2 proteins suggests that BPSL1038 may represent a sub-family of nucleases that share a common ancestor with Cas2.

Ultraviolet Resonance Raman Spectra of Serum Albumins

The ultraviolet resonance Raman (UVRR) spectra of the two proteins bovine serum albumin (BSA) and human serum albumin (HSA) in an aqueous solution are compared with the aim to distinguish between them based on their very similar amino acid composition and structure and to obtain signals from tryptophan that has only very few residues. Comparison of the protein spectra with solutions of tryptophan, tyrosine, and phenylalanine in comparative ratios as in the two proteins shows that at an excitation wavelength of 220 nm, the spectra are dominated by the strong resonant contribution from these three amino acids. While the strong enhancement of two and one single tryptophan residue in BSA and HSA, respectively, results in pronounced bands assigned to fundamental vibrations of tryptophan, its weaker overtones and combination bands do not play a major role in the spectral range above 1800 cm-1. There, the protein spectra clearly reveal the signals of overtones and combination bands of phenylalanine and tyrosine. Assignments of spectral features in the range of Raman shifts from 3800 to 5100 cm-1 to combinations comprising fundamentals and overtones of tyrosine were supported by spectra of amino acid mixtures that contain deuterated tyrosine. The information in the high-frequency region of the UVRR spectra could provide information that is complementary to near-infrared absorption spectroscopy of the proteins.

Cupressus sempervirens L. flavonoids as potent inhibitors to xanthine oxidase: in vitro, molecular docking, ADMET and PASS studies

Excessive intake of purine-rich foods such as seafood and red meat leads to excess xanthine oxidase activity and provokes gout attacks. The aim of this paper is to evaluate in vitro and in silico, the inhibition effect of Cupressus sempervirens plant extracts (flavonoids (Cae) and alkaloids (CaK)) and its six derivative compounds on bovine xanthine oxidase (BXO). The in silico study consists of molecular docking with GOLD v4.0 based on the best PLPchem score (PLP) and prediction of biological activity with the PASS server tool. The inhibitors used were lignan (cp1), Amentoflavone (cp2), Cupressuflavone (cp3), Isocryptomerin (cp4), Hinokiflavone (cp5), and Neolignan (cp6). The in vitro results showed that CaK gives an IC50 of 3.52 ± 0.04 μg/ml. Similarly, Cae saved an IC50 of 8.46 ± 1.98 μg/ml compared with the control (2.82 ± 0.10 μg/ml). The in silico results show that cp1 was the best inhibitor model (PLP of 88.09) with approved pharmacokinetics. These findings suggest that cp1 and cp2 may offer good alternatives for the treatment of hyperuricemia; cp3 was moderate, while the others (cp4 to cp6) were considered weak inhibitors according to their PLP.Communicated by Ramaswamy H. Sarma.

Machine learning-based modulation of Ca2+-binding affinity in EF-hand proteins and comparative structural insights into site-specific cooperative binding

Ca2+-binding proteins are present in almost all living organisms and different types display different levels of binding affinities for the cation. Here, we report two new scoring schemes enabling the user to estimate and manipulate the calcium binding affinities in EF hand containing proteins. To validate this, we designed a unique EF-hand loop capable of binding calcium with high affinity by altering five residues. The N-terminal domain of Entamoeba histolytica calcium-binding protein1 (NtEhCaBP1) is used for site-directed mutagenesis to incorporate the designed loop sequence into the second EF-hand motif of this protein, referred as Nt-EhCaBP1-EF2 mutant. The binding isotherms calculated using ITC calorimetry showed that Nt-EhCaBP1-EF2 mutant site binds Ca2+ with higher affinity than Wt-Nt-EhCaBP1, by ∼600 times. The crystal structure of the mutant displayed more compact Ca2+-coordination spheres in both of its EF loops than the structure of the wildtype protein. The compact coordination sphere of EF-2 causes the bend in the helix-3, which leads to the formation of unexpected hexamer of NtEhCaBP1-EF2 mutant structure. Further dynamic correlation analysis revealed that the mutation in the second EF loop changed the entire residue network of the monomer, resulting in stronger coordination of Ca2+ even in another EF-hand loop.

Unraveling the impact of ORF3a Q57H mutation on SARS-CoV-2: insights from molecular dynamics

ORF3a is a conserved accessory protein of SARS-CoV-2, linked to viral infection and pathogenesis, with acquired mutations at various locations. Previous studies have shown that the occurrence of the Q57H mutation is higher in comparison to other positions in ORF3a. This mutation is known to induce conformational changes, yet the extent of structural alteration and its role in the viral adaptation process remain unknown. Here we performed molecular dynamics (MD) simulations of wt-ORF3a, Q57H, and Q57A mutants to analyze structural changes caused by mutations compared to the native protein. The MD analysis revealed that Q57H and Q57A mutants show significant structural changes in the dimer conformation than the wt-ORF3a. This dimer conformer narrows down the ion channel cavity, which reduces Na + or K + permeability leading to decrease the antigenic response that can help the virus to escape the host immune system. Non-bonding interaction analysis shows the Q57H mutant has more interacting residues, resulting in more stability within dimer conformation than the wt-ORF3a and Q57A. Moreover, both mutant dimers (Q57H and Q57A) form a novel salt-bridge interaction at the same position between A:Asp142 and B:Lys61, whereas such an interaction is absent in the wt-ORF3a dimer. We have also noticed that the TM3 domain's flexibility in Q57H is increased because of strong inter-domain interactions of TM1 and TM2 within the dimer conformation. These unusual interactions and flexibility of Q57H mutant can have significant impacts on the SARS-CoV-2 adaptations, virulence, transmission, and immune system evasion. Our findings are consistent with the previous experimental data and provided details information on the structural perturbation in ORF3a caused by mutations, which can help better understand the structural change at the molecular level as well as the reason for the high virulence properties of this variant.Communicated by Ramaswamy H. Sarma.

Whole Exome-Wide Association Identifies Rare Variants in GALNT9 Associated with Middle Eastern Papillary Thyroid Carcinoma Risk

Papillary thyroid carcinoma (PTC) is the commonest thyroid cancer. The majority of inherited causes of PTC remain elusive. However, understanding the genetic underpinnings and origins remains a challenging endeavor. An exome-wide association study was performed to identify rare germline variants in coding regions associated with PTC risk in the Middle Eastern population. By analyzing exome-sequencing data from 249 PTC patients (cases) and 1395 individuals without any known cancer (controls), GALNT9 emerged as being strongly associated with rare inactivating variants (RIVs) (4/249 cases vs. 1/1395 controls, OR = 22.75, p = 5.09 × 10-5). Furthermore, three genes, TRIM40, ARHGAP23, and SOX4, were enriched for rare damaging variants (RDVs) at the exome-wide threshold (p < 2.5 × 10-6). An additional seven genes (VARS1, ZBED9, PRRC2A, VWA7, TRIM31, TRIM40, and COL8A2) were associated with a Middle Eastern PTC risk based on the sequence kernel association test (SKAT). This study underscores the potential of GALNT9 and other implicated genes in PTC predisposition, illuminating the need for large collaborations and innovative approaches to understand the genetic heterogeneity of PTC predisposition.

Biochemical and structural characterization of the first-discovered metazoan DNA cytosine-N4 methyltransferase from the bdelloid rotifer Adineta vaga

Much is known about the generation, removal, and roles of 5-methylcytosine (5mC) in eukaryote DNA, and there is a growing body of evidence regarding N6-methyladenine, but very little is known about N4-methylcytosine (4mC) in the DNA of eukaryotes. The gene for the first metazoan DNA methyltransferase generating 4mC (N4CMT) was reported and characterized recently by others, in tiny freshwater invertebrates called bdelloid rotifers. Bdelloid rotifers are ancient, apparently asexual animals, and lack canonical 5mC DNA methyltransferases. Here, we characterize the kinetic properties and structural features of the catalytic domain of the N4CMT protein from the bdelloid rotifer Adineta vaga. We find that N4CMT generates high-level methylation at preferred sites, (a/c)CG(t/c/a), and low-level methylation at disfavored sites, exemplified by ACGG. Like the mammalian de novo 5mC DNA methyltransferase 3A/3B (DNMT3A/3B), N4CMT methylates CpG dinucleotides on both DNA strands, generating hemimethylated intermediates and eventually fully methylated CpG sites, particularly in the context of favored symmetric sites. In addition, like DNMT3A/3B, N4CMT methylates non-CpG sites, mainly CpA/TpG, though at a lower rate. Both N4CMT and DNMT3A/3B even prefer similar CpG-flanking sequences. Structurally, the catalytic domain of N4CMT closely resembles the Caulobacter crescentus cell cycle-regulated DNA methyltransferase. The symmetric methylation of CpG, and similarity to a cell cycle-regulated DNA methyltransferase, together suggest that N4CMT might also carry out DNA synthesis-dependent methylation following DNA replication.

In Silico Analysis of Honey Bee Peptides as Potential Inhibitors of Capripoxvirus DNA-Directed RNA Polymerase

The genus Capripoxvirus belongs to the Poxviridae family. The sheeppox, goatpox, and lumpy skin disease viruses are three species of this genus with 96% identity in their genomes. These are financially devastating viral infections among cattle, which cause a reduction in animal products and lead to a loss in livestock industries. In the current study, the phylogenetic analysis was carried out to reveal the evolutionary relationships of Capripoxvirus species (i.e., sheeppox virus (SPPV), goatpox virus (GTPV), and lumpy skin disease virus (LSDV)) with other viruses from the Poxviridae family with >96% query coverage to find the similarity index among all members. The three viruses (i.e., SPPV, GTPV, and LSDV) joined the clade of Capripoxvirus of the Poxviridae family in the phylogenetic tree and exhibited close evolutionary relationships. The multiple sequence alignment using ClustalOmega revealed significant variations in the protein sequences of the DNA-dependent RNA polymerase of SPPV, GTPV, and LSDV. The three-dimensional structures of five selected bee peptides and DNA-directed RNA polymerase of SPPV, GTPV, and LSDV were predicted using trRosetta and I-TASSER and used for molecular docking and simulation studies. The protein-protein docking was carried out using HADDOCK server to explore the antiviral activity of peptides as honey bee proteins against SPPV, GTPV, and LSDV. In total, five peptides were docked to DNA-directed RNA polymerase of these viruses. The peptides mellitin and secapin-1 displayed the lowest binding scores (-106.9 +/- 7.2 kcal/mol and -101.4 +/- 11.3 kcal/mol, respectively) and the best patterns with stable complexes. The molecular dynamics simulation indicated that the complex of protein DNA-dependent RNA polymerase and the peptide melittin stayed firmly connected and the peptide binding to the receptor protein was stable. The findings of this study provide the evidence of bee peptides as potent antimicrobial agents against sheeppox, goatpox, and lumpy skin disease viruses with no complexity.

Exploring the Binding Mechanism of NRG1-ERBB3 Complex and Discovery of Potent Natural Products to Reduce Diabetes-Assisted Breast Cancer Progression

Diabetes mellitus significantly contributes to breast cancer progression, where hyperglycemia upregulates specific genes, leading to more aggressive tumor growth. In patients with BC that develop diabetes, neuregulin 1 (NRG1) and epidermal growth factor receptor 3 (ERBB3) overexpression exacerbate tumor growth and progression. Since the interaction between NRG1 and ERBB3 is critical for tumor growth, understanding the molecular mechanisms underlying NRG1-ERBB3 complex formation is essential for elucidating diabetes-assisted breast cancer progression. However, the key residues forming the NRG1-ERBB3 complex remain unknown. Here, we substituted specific residues in NRG1 with alanine and studied its interactions with ERBB3 using computational structural biology tools. We further screened the South African natural compounds database to target the complex's interface residues to discover potential inhibitors. The conformational stability and dynamic features of NRG1-WT, -H2A, -L3A, and -K35A complexed with ERBB3 were subjected to 400 ns molecular dynamics simulations. The free binding energies of all NRG1-ERBB3 complexes were calculated using the molecular mechanics-generalized Born surface area (MM/GBSA). The H2 and L3 alanine substitutions caused a loss of interaction with ERBB3 residue D73, weakening the interaction with ERBB3. Screening 1300 natural compounds identified four (SANC00643, SANC00824, SANC00975, and SANC00335) with the best potential to inhibit ERRB3-NRG1 coupling. The binding free energies for each complex were - 48.55 kcal/mol for SANC00643, - 47.68 kcal/mol for SANC00824, - 46.04 kcal/mol for SANC00975, and - 45.29 kcal/mol for SANC00335, showing their overall stronger binding with ERBB3 than NRG1 and their potential to act as ERBB3-NRG1 complex inhibitors. In conclusion, this complex may represent a residue-specific drug target to inhibit BC progression.

Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA

Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.

In Silico Investigation of a Chimeric IL24-LK6 Fusion Protein as a Potent Candidate Against Breast Cancer

Targeted delivery of therapeutic anticancer chimeric molecules enhances the efficacy of drug by improving cellular uptake and circulation time. Engineering the molecules to facilitate the specific interaction between chimeric protein and its receptor is critical to elucidate biological mechanism as well as accuracy in modeling of complexes. A theoretically designed novel protein-protein interfaces can serve as a bottom-up method for comprehensive understanding of interacting protein residues. This study was aimed for in silico analyses of a chimeric fusion protein against breast cancer. The amino acid sequences of the interleukin 24 (IL-24) and LK-6 peptide were used to design the chimeric fusion protein via a rigid linker. The secondary and tertiary structures along with physicochemical properties by ProtParam and solubility were predicted using online software. The validation and quality of the fusion protein was confirmed by Rampage and ERRAT2. The newly designed fusion construct has a total length of 179 amino acids. The top-ranked structure from alpha fold2 showed 18.1 KD molecular weight by ProtParam, quality factor of 94.152 by ERRAT, and a valid structure by a Ramachandran plot with 88.5% residues in the favored region. Finally, the docking and simulation studies were performed using HADDOCK and Desmond module of Schrodinger. The quality, validity, interaction analysis, and stability of the fusion protein depict a functional molecule. The fusion gene IL24-LK6 after cloning and expression in a suitable prokaryotic cell might be a useful candidate for developing a novel anticancer therapy.

Evaluation of therapeutic potentials of selected phytochemicals against Nipah virus, a multi-dimensional in silico study

The current study attempted to evaluate the potential of fifty-three (53) natural compounds as Nipah virus attachment glycoprotein (NiV G) inhibitors through in silico molecular docking study. Pharmacophore alignment of the four (4) selected compounds (Naringin, Mulberrofuran B, Rutin and Quercetin 3-galactoside) through Principal Component Analysis (PCA) revealed that common pharmacophores, namely four H bond acceptors, one H bond donor and two aromatic groups were responsible for the residual interaction with the target protein. Out of these four compounds, Naringin was found to have the highest inhibitory potential ( - 9.19 kcal mol^-1) against the target protein NiV G, when compared to the control drug, Ribavirin ( - 6.95 kcal mol^-1). The molecular dynamic simulation revealed that Naringin could make a stable complex with the target protein in the near-native physiological condition. Finally, MM-PBSA (Molecular Mechanics-Poisson-Boltzmann Solvent-Accessible Surface Area) analysis in agreement with our molecular docking result, showed that Naringin ( - 218.664 kJ mol^-1) could strongly bind with the target protein NiV G than the control drug Ribavirin ( - 83.812 kJ mol^-1).

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03595-y.

Design of a multi-epitope Zika virus vaccine candidate - an in-silico study

Zika virus (ZIKV), an RNA virus, rapidly spreads Aedes mosquito-borne sickness. Currently, there are neither effective vaccines nor therapeutics available to prevent or treat ZIKV infection. In this study, to address these unmet medical needs, we aimed to design B- and T-cell candidate multi-epitope-based subunit against ZIKV using an in silico approach. In this study we applied immunoinformatics, molecular docking, and dynamic simulation assessments targeting the most immunogenic proteins; the capsid (C), envelope (E) proteins and the non-stuctural protein (NS1), described in our previous study, and which predicted immunodominant B and T cell epitopes. The final non-allergenic and highly antigenic multi-epitope was constituted of immunogenic screened-epitopes (3 CTL and 3 HTL) and the β-defensin as an adjuvant that have been linked using EAAAK, AAY, and GPGPG linkers, respectively. The final construct containing 143 amino acids was characterized for its allergenicity, antigenicity, and physiochemical properties; and found to be safe and immunogenic with a good prediction of solubility. The existence of IFN-γ epitopes asserts the capacity to trigger strong immune responses. Subsequently, the molecular docking among vaccine and immune receptors (TLR2/TLR4) was revealed with a good binding affinity with and stable molecular interactions. Molecular dynamics simulation confirmed the stability of the complexes. Finally, the construct was subjected to in silico cloning demonstrating the efficiently of its expression in E.coli. However, this study needs the experimental validation to demonstrate vaccine safety and efficacy.Communicated by Ramaswamy H. Sarma.

Cryo-EM structure of the folded-back state of human β-cardiac myosin

To save energy and precisely regulate cardiac contractility, cardiac muscle myosin heads are sequestered in an 'off' state that can be converted to an 'on' state when exertion is increased. The 'off' state is equated with a folded-back structure known as the interacting-heads motif (IHM), which is a regulatory feature of all class-2 muscle and non-muscle myosins. We report here the human β-cardiac myosin IHM structure determined by cryo-electron microscopy to 3.6 Å resolution, providing details of all the interfaces stabilizing the 'off' state. The structure shows that these interfaces are hot spots of hypertrophic cardiomyopathy mutations that are thought to cause hypercontractility by destabilizing the 'off' state. Importantly, the cardiac and smooth muscle myosin IHM structures dramatically differ, providing structural evidence for the divergent physiological regulation of these muscle types. The cardiac IHM structure will facilitate development of clinically useful new molecules that modulate IHM stability.

Cryo-EM structure of the folded-back state of human β-cardiac myosin

During normal levels of exertion, many cardiac muscle myosin heads are sequestered in an off-state even during systolic contraction to save energy and for precise regulation. They can be converted to an on-state when exertion is increased. Hypercontractility caused by hypertrophic cardiomyopathy (HCM) myosin mutations is often the result of shifting the equilibrium toward more heads in the on-state. The off-state is equated with a folded-back structure known as the interacting head motif (IHM), which is a regulatory feature of all muscle myosins and class-2 non-muscle myosins. We report here the human β-cardiac myosin IHM structure to 3.6 Å resolution. The structure shows that the interfaces are hot spots of HCM mutations and reveals details of the significant interactions. Importantly, the structures of cardiac and smooth muscle myosin IHMs are dramatically different. This challenges the concept that the IHM structure is conserved in all muscle types and opens new perspectives in the understanding of muscle physiology. The cardiac IHM structure has been the missing puzzle piece to fully understand the development of inherited cardiomyopathies. This work will pave the way for the development of new molecules able to stabilize or destabilize the IHM in a personalized medicine approach. *This manuscript was submitted to Nature Communications in August 2022 and dealt efficiently by the editors. All reviewers received this version of the manuscript before 9 208 August 2022. They also received coordinates and maps of our high resolution structure on the 18 208 August 2022. Due to slowness of at least one reviewer, this contribution was delayed for acceptance by Nature Communications and we are now depositing in bioRxiv the originally submitted version written in July 2022 for everyone to see. Indeed, two bioRxiv contributions at lower resolution but adding similar concepts on thick filament regulation were deposited this week in bioRxiv, one of the contributions having had access to our coordinates. We hope that our data at high resolution will be helpful for all readers that appreciate that high resolution information is required to build accurate atomic models and discuss implications for sarcomere regulation and the effects of cardiomyopathy mutations on heart muscle function.

Evolutionary Aspects of Selenium Binding Protein (SBP)

Selenium-binding proteins represent a ubiquitous protein family and recently SBP1 was described as a new stress response regulator in plants. SBP1 has been characterized as a methanethiol oxidase, however its exact role remains unclear. Moreover, in mammals, it is involved in the regulation of anti-carcinogenic growth and progression as well as reduction/oxidation modulation and detoxification. In this work, we delineate the functional potential of certain motifs of SBP in the context of evolutionary relationships. The phylogenetic profiling approach revealed the absence of SBP in the fungi phylum as well as in most non eukaryotic organisms. The phylogenetic tree also indicates the differentiation and evolution of characteristic SBP motifs. Main evolutionary events concern the CSSC motif for which Acidobacteria, Fungi and Archaea carry modifications. Moreover, the CC motif is harbored by some bacteria and remains conserved in Plants, while modified to CxxC in Animals. Thus, the characteristic sequence motifs of SBPs mainly appeared in Archaea and Bacteria and retained in Animals and Plants. Our results demonstrate the emergence of SBP from bacteria and most likely as a methanethiol oxidase.

Interactions of heparin derivatives with recombinant human keratinocyte growth factor: Structural stability and bioactivity effect study

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.

Integrated ATAC-Seq and RNA-Seq Data Analysis to Reveal OsbZIP14 Function in Rice in Response to Heat Stress

Transcription factors (TFs) play critical roles in mediating the plant response to various abiotic stresses, particularly heat stress. Plants respond to elevated temperatures by modulating the expression of genes involved in diverse metabolic pathways, a regulatory process primarily governed by multiple TFs in a networked configuration. Many TFs, such as WRKY, MYB, NAC, bZIP, zinc finger protein, AP2/ERF, DREB, ERF, bHLH, and brassinosteroids, are associated with heat shock factor (Hsf) families, and are involved in heat stress tolerance. These TFs hold the potential to control multiple genes, which makes them ideal targets for enhancing the heat stress tolerance of crop plants. Despite their immense importance, only a small number of heat-stress-responsive TFs have been identified in rice. The molecular mechanisms underpinning the role of TFs in rice adaptation to heat stress still need to be researched. This study identified three TF genes, including OsbZIP14, OsMYB2, and OsHSF7, by integrating transcriptomic and epigenetic sequencing data analysis of rice in response to heat stress. Through comprehensive bioinformatics analysis, we demonstrated that OsbZIP14, one of the key heat-responsive TF genes, contained a basic-leucine zipper domain and primarily functioned as a nuclear TF with transcriptional activation capability. By knocking out the OsbZIP14 gene in the rice cultivar Zhonghua 11, we observed that the knockout mutant OsbZIP14 exhibited dwarfism with reduced tiller during the grain-filling stage. Under high-temperature treatment, it was also demonstrated that in the OsbZIP14 mutant, the expression of the OsbZIP58 gene, a key regulator of rice seed storage protein (SSP) accumulation, was upregulated. Furthermore, bimolecular fluorescence complementation (BiFC) experiments uncovered a direct interaction between OsbZIP14 and OsbZIP58. Our results suggested that OsbZIP14 acts as a key TF gene through the concerted action of OsbZIP58 and OsbZIP14 during rice filling under heat stress. These findings provide good candidate genes for genetic improvement of rice but also offer valuable scientific insights into the mechanism of heat tolerance stress in rice.

B-box1 Domain of MID1 Interacts with the Ube2D1 E2 Enzyme Differently Than RING E3 Ligases

The MID1 TRIM protein is important for ventral midline development in vertebrates, and mutations of its B-box1 domain result in several birth defects. The B-box1 domain of the human MID1 protein binds two zinc atoms and adopt a similar ββα-RING structure. This domain is required for the efficient ubiquitination of protein phosphatase 2A, alpha4, and fused kinase. Considering the structural similarity, the MID1 B-box1 domain exhibits mono-autoubiquitination activity, in contrast to poly-autoubiquitination observed for RING E3 ligases. To understand its mechanism of action, the interaction of the B-box1 domain with Ube2D1 (UbcH5a, E2), a preferred E2 ligase, is investigated. Using isothermal titration calorimetry, the MID1 RING and B-box1 domains were observed to have similar binding affinities with the Ube2D1 protein. However, NMR ¹⁵N-¹H Heteronuclear Single Quantum Coherence titration, ¹⁵N relaxation data, and High Ambiguity Driven protein-protein DOCKing (HADDOCK) calculations show the B-box1 domain binding on a surface distinct from where RING domains bind. The novel binding interaction shows the B-box1 domain partially overlapping the noncovalent Ube2D1 and a ubiquitin binding site that is necessary for poly-autoubiquitination activity. The B-box1 domain also displaces the ubiquitin from the Ube2D1 protein. These studies reveal a novel binding interaction between the zinc-binding ββα-fold B-box1 domain and the Ube2D enzyme family and that this difference in binding, compared to RING E3 ligases, provides a rationale for its auto-monoubiquitination E3 ligase activity.

Prediction and Evaluation of Bioactive Properties of Cowpea Protein Hydrolysates

Cowpea protein hydrolysates were prepared using thermolysin, alcalase, and trypsin and analysed for bioactive properties, and then, the release of bioactive peptides was investigated in silico. It was found that the degree of hydrolysis reached 48% after 24 h hydrolysis with alcalase. The hydrolysate prepared using alcalase showed higher ACE inhibitory (62%) and DPPH scavenging activity (19%). SDS-PAGE analysis revealed that vignin was the major protein in cowpea protein isolate. In silico analysis indicated the presence of potential bioactive peptides with potent bioactivity in the primary structure of proteins. The 3D structure of proteins was built, upon which bioactive peptides were mapped using their location in the primary structure. The secondary structure and solvent accessible surface around each bioactive peptide were then calculated. On this basis, the higher degree of hydrolysis and bioactive properties of cowpea protein hydrolysate prepared by alcalase were explained, and structural factors influencing the release of bioactive peptides were investigated.

Immunoinformatics Study: Multi-Epitope Based Vaccine Design from SARS-CoV-2 Spike Glycoprotein

The coronavirus disease 2019 outbreak has become a huge challenge in the human sector for the past two years. The coronavirus is capable of mutating at a higher rate than other viruses. Thus, an approach for creating an effective vaccine is still needed to induce antibodies against multiple variants with lower side effects. Currently, there is a lack of research on designing a multiepitope of the COVID-19 spike protein for the Indonesian population with comprehensive immunoinformatic analysis. Therefore, this study aimed to design a multiepitope-based vaccine for the Indonesian population using an immunoinformatic approach. This study was conducted using the SARS-CoV-2 spike glycoprotein sequences from Indonesia that were retrieved from the GISAID database. Three SARS-CoV-2 sequences, with IDs of EIJK-61453, UGM0002, and B.1.1.7 were selected. The CD8+ cytotoxic T-cell lymphocyte (CTL) epitope, CD4+ helper T lymphocyte (HTL) epitope, B-cell epitope, and IFN-γ production were predicted. After modeling the vaccines, molecular docking, molecular dynamics, in silico immune simulations, and plasmid vector design were performed. The designed vaccine is antigenic, non-allergenic, non-toxic, capable of inducing IFN-γ with a population reach of 86.29% in Indonesia, and has good stability during molecular dynamics and immune simulation. Hence, this vaccine model is recommended to be investigated for further study.

Structural Consequences of IRS-2 nsSNPs and Implication for Insulin Receptor Substrate-2 Protein Stability

Single-Nucleotide Polymorphisms (SNPs) are common genetic variations implicated in human diseases. The non-synonymous SNPs (nsSNPs) affect the proteins' structures and their molecular interactions with other interacting proteins during the accomplishment of biochemical processes. This ultimately causes proteins functional perturbation and disease phenotypes. The Insulin receptor substrate-2 (IRS-2) protein promotes glucose absorption and participates in the biological regulation of glucose metabolism and energy production. Several IRS-2 SNPs are reported in association with type 2 diabetes and obesity in human populations. However, there are no comprehensive reports about the protein structural consequences of these nsSNPs. Keeping in view the pathophysiological consequences of the IRS-2 nsSNPs, we designed the current study to understand their possible structural impact on coding protein. The prioritized list of the deleterious IRS-2 nsSNPs was acquired from multiple bioinformatics resources, including VEP (SIFT, PolyPhen, and Condel), PROVEAN, SNPs&GO, PMut, and SNAP2. The protein structure stability assessment of these nsSNPs was performed by MuPro and I-Mutant-3.0 servers via structural modeling approaches. The atomic-level structural and molecular dynamics (MD) impact of these nsSNPs were examined using GROMACS 2019.2 software package. The analyses initially predicted 8 high-risk nsSNPs located in the highly conserved regions of IRS-2. The MD simulation analysis eventually prioritized the N232Y, R218C, and R104H nsSNPs that predicted to significantly compromise the structure stability and may affect the biological function of IRS-2. These nsSNPs are predicted as high-risk candidates for diabetes and obesity. The validation of protein structural impact of these shortlisted nsSNPs may provide biochemical insight into the IRS-2-mediated type-2 diabetes.

Conformational and mechanical stability of the isolated large subunit of membrane-bound [NiFe]-hydrogenase from Cupriavidus necator

Comprising at least a bipartite architecture, the large subunit of [NiFe]-hydrogenase harbors the catalytic nickel-iron site while the small subunit houses an array of electron-transferring Fe-S clusters. Recently, some [NiFe]-hydrogenase large subunits have been isolated showing an intact and redox active catalytic cofactor. In this computational study we have investigated one of these metalloproteins, namely the large subunit HoxG of the membrane-bound hydrogenase from Cupriavidus necator (CnMBH), targeting its conformational and mechanical stability using molecular modelling and long all-atom Gaussian accelerated molecular dynamics (GaMD). Our simulations predict that isolated HoxG is stable in aqueous solution and preserves a large portion of its mechanical properties, but loses rigidity in regions around the active site, in contrast to the MBH heterodimer. Inspired by biochemical data showing dimerization of the HoxG protein and IR measurements revealing an increased stability of the [NiFe] cofactor in protein preparations with higher dimer content, corresponding simulations of homodimeric forms were also undertaken. While the monomeric subunit contains several flexible regions, our data predicts a regained rigidity in homodimer models. Furthermore, we computed the electrostatic properties of models obtained by enhanced sampling with GaMD, which displays a significant amount of positive charge at the protein surface, especially in solvent-exposed former dimer interfaces. These data offer novel insights on the way the [NiFe] core is protected from de-assembly and provide hints for enzyme anchoring to surfaces, which is essential information for further investigations on these minimal enzymes.

Characterisation of Bioactive Peptides from Red Alga Gracilariopsis chorda

In this study, we studied the bioactive peptides produced by thermolysin hydrolysis of a water-soluble protein (WSP) from the red alga Gracilariopsis chorda, whose major components are phycobiliproteins and Ribulose-1,5-bisphosphate carboxylase-oxygenase (RuBisCo). The results showed that WSP hydrolysate exhibited significantly higher ACE inhibitory activity (92% inhibition) compared to DPP-IV inhibitory activity and DPPH scavenging activity. The phycobiliproteins and RuBisCo of G. chorda contain a high proportion of hydrophobic (31.0-46.5%) and aromatic (5.1-46.5%) amino acid residues, which was considered suitable for the formation of peptides with strong ACE inhibitory activity. Therefore, we searched for peptides with strong ACE inhibitory activity and identified two novel peptides (IDHY and LVVER). Then, their interaction with human ACE was evaluated by molecular docking, and IDHY was found to be a promising inhibitor. In silico analysis was then performed on the structural factors affecting ACE inhibitory peptide release, using the predicted 3D structures of phycobiliproteins and RuBisCo. The results showed that most of the ACE inhibitory peptides are located in the highly solvent accessible α-helix. Therefore, it was suggested that G. chorda is a good source of bioactive peptides, especially ACE-inhibitory peptides.

Assessment of potassium ion channel during electric signalling in biofilm formation of Acinetobacter baumannii for finding antibiofilm molecule

Acinetobacter baumannii is an opportunistic ESKAPE pathogen which causes nosocomial infections and can produce biofilms that act as resistant determinants. The role of quorum sensing (chemical signaling) in biofilm establishment has already been studied extensively, but the existence of electrochemical signaling during biofilm formation by A. baumannii has not yet been investigated. The current study evaluated the presence of electrical signaling, types of ion channels involved, and their role in biofilm formation using spectroscopic and microbiological methods. The findings suggest that the potassium ion channel has a significant role in the electrical signaling during the biofilm formation by A. baumannii. Further, in-silico screening, molecular mechanics, and molecular dynamic simulation studies identify a potential lead, ZINC12496555(a specific inhibitor), which targets the potassium ion channel protein of A. baumannii. Mutational analysis of the interacting residues showed alterations in the unfolding rate of this protein after the selected mutation, which shows its role in the stability of this protein. It was also observed that identified lead has high antibiofilm activity, no human off-targets, and non-cytotoxicity to cell lines. Thus, identified lead against the potassium channel of A baumannii may be used as an effective therapeutic for the treatment of A. baumannii infections after further experimental validation.

Exploring the structural basis to develop efficient multi-epitope vaccines displaying interaction with HLA and TAP and TLR3 molecules to prevent NIPAH infection, a global threat to human health

Nipah virus (NiV) is an emerging zoonotic virus that caused several serious outbreaks in the south asian region with high mortality rates ranging from 40 to 90% since 2001. NiV infection causes lethal encephalitis and respiratory disease with the symptom of endothelial cell-cell fusion. No specific and effective vaccine has yet been reported against NiV. To address the urgent need for a specific and effective vaccine against NiV infection, in the present study, we have designed two Multi-Epitope Vaccines (MEVs) composed of 33 Cytotoxic T lymphocyte (CTL) epitopes and 38 Helper T lymphocyte (HTL) epitopes. Out of those CTL and HTL combined 71 epitopes, 61 novel epitopes targeting nine different NiV proteins were not used before for vaccine design. Codon optimization for the cDNA of both the designed MEVs might ensure high expression potential in the human cell line as stable proteins. Both MEVs carry potential B cell linear epitope overlapping regions, B cell discontinuous epitopes as well as IFN-γ inducing epitopes. Additional criteria such as sequence consensus amongst CTL, HTL and B Cell epitopes was implemented for the design of final constructs constituting MEVs. Hence, the designed MEVs carry the potential to elicit cell-mediated as well as humoral immune response. Selected overlapping CTL and HTL epitopes were validated for their stable molecular interactions with HLA class I and II alleles and in case of CTL epitopes with human Transporter Associated with antigen Processing (TAP) cavity. The structure based epitope cross validation for interaction with TAP cavity was used as another criteria choosing final epitopes for NiV MEVs. Finally, human Beta-defensin 2 and Beta-defensin 3 were used as adjuvants to enhance the immune response of both the MEVs. Molecular dynamics simulation studies of MEVs-TLR3 ectodomain (Human Toll-Like Receptor 3) complex indicated the stable molecular interaction. We conclude that the MEVs designed and in silico validated here could be highly potential vaccine candidates to combat NiV infections, with great effectiveness, high specificity and large human population coverage worldwide.

Integrative Bioinformatics Approaches Indicate a Particular Pattern of Some SARS-CoV-2 and Non-SARS-CoV-2 Proteins

Pattern recognition plays a critical role in integrative bioinformatics to determine the structural patterns of proteins of viruses such as SARS-CoV-2. This study identifies the pattern of SARS-CoV-2 proteins to depict the structure-function relationships of the protein alphabets of SARS-CoV-2 and COVID-19. The assembly enumeration algorithm, Anisotropic Network Model, Gaussian Network Model, Markovian Stochastic Model, and image comparison protein-like alphabets were used. The distance score was the lowest with 22 for "I" and highest with 40 for "9". For post-processing and decision, two protein alphabets "C" (PDB ID: 6XC3) and "S" (PDB ID: 7OYG) were evaluated to understand the structural, functional, and evolutionary relationships, and we found uniqueness in the functionality of proteins. Here, models were constructed using "SARS-CoV-2 proteins" (12 numbers) and "non-SARS-CoV-2 proteins" (14 numbers) to create two words, "SARS-CoV-2" and "COVID-19". Similarly, we developed two slogans: "Vaccinate the world against COVID-19" and "Say no to SARS-CoV-2", which were made with the proteins structure. It might generate vaccine-related interest to broad reader categories. Finally, the evolutionary process appears to enhance the protein structure smoothly to provide suitable functionality shaped by natural selection.

High-resolution cryo-EM structure of the Shigella virus Sf6 genome delivery tail machine

Sf6 is a bacterial virus that infects the human pathogen Shigella flexneri. Here, we describe the cryo-electron microscopy structure of the Sf6 tail machine before DNA ejection, which we determined at a 2.7-angstrom resolution. We built de novo structures of all tail components and resolved four symmetry-mismatched interfaces. Unexpectedly, we found that the tail exists in two conformations, rotated by ~6° with respect to the capsid. The two tail conformers are identical in structure but differ solely in how the portal and head-to-tail adaptor carboxyl termini bond with the capsid at the fivefold vertex, similar to a diamond held over a five-pronged ring in two nonidentical states. Thus, in the mature Sf6 tail, the portal structure does not morph locally to accommodate the symmetry mismatch but exists in two energetic minima rotated by a discrete angle. We propose that the design principles of the Sf6 tail are conserved across P22-like Podoviridae.

Comparative computational study to augment UbiA prenyltransferases inherent in purple photosynthetic bacteria cultured from mangrove microbial mats in Qatar for coenzyme Q10 biosynthesis

Coenzyme Q10 (CoQ10) is a powerful antioxidant with a myriad of applications in healthcare and cosmetic industries. The most effective route of CoQ10 production is microbial biosynthesis. In this study, four CoQ10 biosynthesizing purple photosynthetic bacteria: Rhodobacter blasticus, Rhodovulum adriaticum, Afifella pfennigii and Rhodovulum marinum, were identified using 16S rRNA sequencing of enriched microbial mat samples obtained from Purple Island mangroves (Qatar). The membrane bound enzyme 4-hydroxybenzoate octaprenyltransferase (UbiA) is pivotal for bacterial biosynthesis of CoQ10. The identified bacteria could be inducted as efficient industrial bio-synthesizers of CoQ10 by engineering their UbiA enzymes. Therefore, the mutation sites and substitution residues for potential functional enhancement were determined by comparative computational study. Two mutation sites were identified within the two conserved Asp-rich motifs, and the effect of proposed mutations in substrate binding affinity of the UbiA enzymes was assessed using multiple ligand simultaneous docking (MLSD) studies, as a groundwork for experimental studies.

Immunoinformatics design of multi-epitope peptide-based vaccine against Haemophilus influenzae strain using cell division protein

Haemophilus influenzae is a pathogen that causes invasive bacterial infections in humans. The highest prevalence lies in both young children and adults. Generally, there are no vaccines available that target all the strains of Haemophilus influenzae. Hence, the purpose of this research is to employ bioinformatics and immunoinformatics approaches to design a Multi-Epitope Vaccine candidate employing the pathogenic cell division protein FtsN that specifically combat all the Haemophilus influenzae strains. The current research focuses on developing subunit vaccine in contrast to vaccines generated from the entire pathogen. This will be accomplished by combining multiple bioinformatics and immunoinformatics approaches. As a result, prospective T cells (helper T lymphocyte and cytotoxic T lymphocytes) and B cells epitopes were investigated. The human leukocyte antigen allele having strong associations with the antigenic and overlapping epitopes were chosen, with 70% of the total coverage of the world population. To construct a linked vaccine design, multiple linkers were used. To increase the immunogenic profile, an adjuvant was linked using EAAAK linker. The final vaccine construct with 149 amino acids was obtained after adjuvants and linkers were added. The developed Multi-Epitope Vaccine has a high antigenicity as well as viable physiochemical features. The 3D conformation was modeled and undergoes refinement and validation using bioinformatics methods. Furthermore, protein-protein molecular docking analysis was performed to predict the effective binding poses of Multi-Epitope Vaccine with the Toll-like receptor 4 protein. Besides, vaccine underwent the codon translational optimization and computational cloning to verify the reliability and proper Multi-Epitope Vaccine expression. In addition, it is necessary to conduct experiments and research in the laboratory to demonstrate that the vaccine that has been developed is immunogenic and protective.

Novel Copper Oxide Bio-Nanocrystals to Target Outer Membrane Lectin of Vancomycin-Resistant Enterococcus faecium (VREfm): In Silico, Bioavailability, Antimicrobial, and Anticancer Potential

In present study, we used Olea europaea leaf extract to biosynthesize in situ Copper Oxide nanocrystals (CuO @OVLe NCs) with powerful antibacterial and anti-cancer capabilities. Physio-chemical analyses, such as UV/Vis, FTIR, XRD, EDX, SEM, and TEM, were applied to characterize CuO @OVLe NCs. The UV/Vis spectrum demonstrated a strong peak at 345 nm. Furthermore, FTIR, XRD, and EDX validated the coating operation's contact with colloidal CuO @OVLe NCs. According to TEM and SEM analyses, CuO @OVLe NCs exhibited a spherical shape and uniform distribution of size with aggregation, for an average size of ~75 nm. The nanoparticles demonstrated a considerable antibacterial effect against E. faecium bacterial growth, as well as an increased inhibition rate in a dose-dependent manner on the MCF-7, PC3, and HpeG2 cancer cell lines and a decreased inhibition rate on WRL-68. Molecular docking and MD simulation were used to demonstrate the high binding affinity of a ligand (Oleuropein) toward the lectin receptor complex of the outer membrane to vancomycin-resistant E. faecium (VREfm) via amino acids (Leu 195, Thr 288, His 165, and Ser 196). Hence, our results expand the accessibility of OVLe's bioactive components as a promising natural source for the manufacture of physiologically active components and the creation of green biosynthesis of metal nanocrystals.

Parallel functional annotation of cancer-associated missense mutations in histone methyltransferases

Using exome sequencing for biomarker discovery and precision medicine requires connecting nucleotide-level variation with functional changes in encoded proteins. However, for functionally annotating the thousands of cancer-associated missense mutations, or variants of uncertain significance (VUS), purifying variant proteins for biochemical and functional analysis is cost-prohibitive and inefficient. We describe parallel functional annotation (PFA) of large numbers of VUS using small cultures and crude extracts in 96-well plates. Using members of a histone methyltransferase family, we demonstrate high-throughput structural and functional annotation of cancer-associated mutations. By combining functional annotation of paralogs, we discovered two phylogenetic and clustering parameters that improve the accuracy of sequence-based functional predictions to over 90%. Our results demonstrate the value of PFA for defining oncogenic/tumor suppressor functions of histone methyltransferases as well as enhancing the accuracy of sequence-based algorithms in predicting the effects of cancer-associated mutations.

Development of CD163 receptor-based enzyme-linked immunosorbent assay for diagnosis of porcine reproductive and respiratory syndrome virus

Porcine reproductive and respiratory syndrome (PRRS) is an important economical disease in the global swine industry. The accurate detection of the PRRS virus (PRRSV) antigen is essential for the disease control and prevention programme. In this study, an indirect enzyme-linked immunosorbent test (PRRSVCD163-iELISA) was developed for the detection of the PRRSV antigen in samples of post-mortem swine tissue using the recombinant pig CD163 receptor protein as the capture ligand. The test was found to be specific for PRRSV, with no cross-reactions with other prevalent pig viral pathogens. The assay was validated by testing 217 post-mortem porcine tissue samples and the results were found to be satisfactory with a relative accuracy of 88.88%. Our assay is also quite precise, with intra- and inter-assay CVs of 6% and 10%, respectively. These findings imply that the PRRSVCD163-iELISA developed is capable of detecting the PRRSV antigen in swine post-mortem tissue samples. This research showed that porcine CD163, the PRRSV cellular receptor, can be exploited to build a diagnostic technique for the detection of PRRSV antigen.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03376-z.

Terminase Subunits from the Pseudomonas-Phage E217

Pseudomonas phages are increasingly important biomedicines for phage therapy, but little is known about how these viruses package DNA. This paper explores the terminase subunits from the Myoviridae E217, a Pseudomonas-phage used in an experimental cocktail to eradicate P. aeruginosa in vitro and in animal models. We identified the large (TerL) and small (TerS) terminase subunits in two genes ∼58 kbs away from each other in the E217 genome. TerL presents a classical two-domain architecture, consisting of an N-terminal ATPase and C-terminal nuclease domain arranged into a bean-shaped tertiary structure. A 2.05 Å crystal structure of the C-terminal domain revealed an RNase H-like fold with two magnesium ions in the nuclease active site. Mutations in TerL residues involved in magnesium coordination had a dominant-negative effect on phage growth. However, the two ions identified in the active site were too far from each other to promote two-metal-ion catalysis, suggesting a conformational change is required for nuclease activity. We also determined a 3.38 Å cryo-EM reconstruction of E217 TerS that revealed a ring-like decamer, departing from the most common nonameric quaternary structure observed thus far. E217 TerS contains both N-terminal helix-turn-helix motifs enriched in basic residues and a central channel lined with basic residues large enough to accommodate double-stranded DNA. Overexpression of TerS caused a more than a 4-fold reduction of E217 burst size, suggesting a catalytic amount of the protein is required for packaging. Together, these data expand the molecular repertoire of viral terminase subunits to Pseudomonas-phages used for phage therapy.