Immune targeting of filarial glutaredoxin through a multi-epitope peptide-based vaccine: A reverse vaccinology approach.

Despite the efforts of global programme to eliminate lymphatic filariasis (GPELF), the threat of lymphatic filariasis (LF) still looms over humanity in terms of long-term disabilities, and morbidities across the globe. In light of this situation, investigators have chosen to focus on the development of immunotherapeutics targeting the physiologically important filarial-specific proteins. Glutaredoxin (16.43 kDa) plays a pivotal role in filarial redox biology, serving as a vital contributor. In the context of the intra-host survival of filarial parasites, this antioxidant helps in mitigating the oxidative stress imposed by the host immune system. Given its significant contribution, the development of a vaccine targeting glutaredoxin holds promise as a new avenue for achieving a filaria-free world. Herein, multi-epitope-based vaccine was designed using advanced immunoinformatics approach. Initially, 4B-cell epitopes and 6 T-cell epitopes (4 MHC I and 2 MHC II) were identified from the 146 amino acid long sequence of glutaredoxin of the human filarid, Wuchereria bancrofti. Subsequent clustering of these epitopes with linker peptides finalized the vaccine structure. To boost TLR-mediated innate immunity, TLR-specific adjuvants were incorporated into the designed vaccine. After that, experimental analyses confirm the designed vaccine, Vac4 as anefficient ligand of human TLR5 to elicit protective innate immunity against filarial glutaredoxin. Immune simulation further demonstrated abundant levels of IgG and IgM as crucial contributors in triggering vaccine-induced adaptive responses in the recipients. Hence, to facilitate the validation of immunogenicity of the designed vaccine, Vac4 was cloned in silico in pET28a(+) expression vector for recombinant production. Taken together, our findings suggest that vaccine-mediated targeting of filarial glutaredoxin could be a future option for intervening LF on a global scale.

Insights into the binding mechanism of ascorbic acid and violaxanthin with violaxanthin de-epoxidase (VDE) and chlorophycean violaxanthin de-epoxidase (CVDE) enzymes.

Photosynthetic organisms have evolved to work under low and high lights in photoprotection, acting as a scavenger of reactive oxygen species. The light-dependent xanthophyll cycle involved in this process is performed by a key enzyme (present in the thylakoid lumen), Violaxanthin De-Epoxidase (VDE), in the presence of violaxanthin (Vio) and ascorbic acid substrates. Phylogenetically, VDE is found to be connected with an ancestral enzyme Chlorophycean Violaxanthin De-Epoxidase (CVDE), present in the green algae on the stromal side of the thylakoid membrane. However, the structure and functions of CVDE were not known. In search of functional similarities involving this cycle, the structure, binding conformation, stability, and interaction mechanism of CVDE are explored with the two substrates compared to VDE. The structure of CVDE was determined by homology modeling and validated. In silico docking (of first-principles optimized substrates) revealed it has a larger catalytic domain than VDE. A thorough analysis of the binding affinity and stability of four enzyme-substrate complexes is performed by computing free energies and their decomposition, the root-mean-square deviation (RMSD) and fluctuation (RMSF), the radius of gyration, salt bridge, and hydrogen bonding interactions in molecular dynamics. Based on these, violaxanthin interacts with CVDE to a similar extent as that of VDE. Hence, its role is expected to be the same for both enzymes. On the contrary, ascorbic acid has a weaker interaction with CVDE than VDE. Given these interactions drive epoxidation or de-epoxidation in the xanthophyll cycle, it immediately discerns that either ascorbic acid does not participate in de-epoxidation or a different cofactor is necessary as CVDE has a weaker interaction with ascorbic acid than VDE.

Potential targets of severe acute respiratory syndrome coronavirus 2 of clinical drug fluvoxamine: Docking and molecular dynamics studies to elucidate viral action.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued evolving for survival and adaptation by mutating itself into different variants of concern, including omicron. Several studies and clinical trials found fluvoxamine, an Food and Drug Administration-approved antidepressant drug, to be effective at preventing mild coronavirus disease 2019 (COVID-19) from progressing to severe diseases. However, the mechanism of fluvoxamine's direct antiviral action against COVID-19 is still unknown. Fluvoxamine was docked with 11 SARS-CoV-2 targets and subjected to stability, conformational changes, and binding free energy analyses to explore its mode of action. Of the targets, nonstructural protein 14 (NSP14), main protease (Mpro), and papain-like protease (PLpro) had the best docking scores with fluvoxamine. Consistent with the docking results, it was confirmed by molecular dynamics simulations that the NSP14 N7-MTase ((N7-guanine)-methyltransferase)-fluvoxamine, Mpro-fluvoxamine, and PLpro-fluvoxamine complexes are stable, with the lowest binding free energies of -105.1, -82.7, and - 38.5 kJ/mol, respectively. A number of hotspot residues involved in the interaction were also identified. These include Glu166, Asp187, His41, and Cys145 in Mpro, Gly163 and Arg166 in PLpro, and Glu302, Gly333, and Phe426 in NSP14, which could aid in the development of better antivirals against SARS-CoV-2.

Reverse vaccinology assisted design of a novel multi-epitope vaccine to target Wuchereria bancrofti cystatin: An immunoinformatics approach.

Proteases are the critical mediators of immunomodulation exerted by the filarial parasites to bypass and divert host immunity. Cystatin is a small (∼15 kDa) immunomodulatory filarial protein and known to contribute in the immunomodulation strategy by inducing anti-inflammatory response through alternative activation of macrophages. Recently, Wuchereria bancrofti cystatin has been discovered as a ligand of human toll-like receptor 4 which is key behind the cystatin-induced anti-inflammatory response in major human antigen-presenting cells. Considering the pivotal role of cystatin in the immunobiology of filariasis, cystatin could be an efficacious target for developing vaccine. Herein, we present the design and in-silico analyses of a multi-epitope-based peptide vaccine to target W. bancrofti cystatin through immune-informatics approaches. The 262 amino acid long antigen construct comprises 9 MHC-I epitopes and MHC-II epitopes linked together by GPGPG peptide alongside an adjuvant (50S ribosomal protein L7/L12) at N terminus and 6 His tags at C terminus. Molecular docking study reveals that the peptide could trigger TLR4-MD2 to induce protective innate immune responses while the induced adaptive responses were found to be mediated by IgG, IgM and Th1 mediated responses. Notably, the designed vaccine exhibits high stability and no allergenicity in-silico. Furthermore, the muti epitope-vaccine was also predicted for its RNA structure and cloned in pET30ax for further experimental validation. Taken together, this study presents a novel multi-epitope peptide vaccine for triggering efficient innate and adaptive immune responses against W. bancrofti to intervene LF through immunotherapy.

Designing efficient multi-epitope peptide-based vaccine by targeting the antioxidant thioredoxin of bancroftian filarial parasite.

Thioredoxin is a low molecular weight redox-active protein of filarial parasite that plays a crucial role in downregulating the host immune response to prolong the survival of the parasite within the host body. It has the ability to cope up with the oxidative challenges posed by the host. Hence, the antioxidant protein of the filarial parasite has been suggested to be a useful target for immunotherapeutic intervention of human filariasis. In this study, we have designed a multi-epitope peptide-based vaccine using thioredoxin of Wuchereria bancrofti. Different MHC-I and MHC-II epitopes were predicted using various web servers to construct the vaccine model as MHC-I and MHC-II epitopes are crucial for the development of both humoral and cellular immune responses. Moreover, TLRs specific adjuvants were also incorporated into the vaccine candidates as TLRs are the key immunomodulator to execute innate immunity. Protein-protein molecular docking and simulation analysis between the vaccine and human TLR was performed. TLR5 is the most potent receptor to convey the vaccine-mediated inductive signal for eliciting an innate immune response. A satisfactory immunogenic report from an in-silico immune simulation experiment directed us to propose our vaccine model for experimental and clinical validation. The reverse translated vaccine sequence was also cloned in pET28a(+) to apply the concept in a wet lab experiment in near future. Taken together, this in-silico study on the design of a vaccine construct to target W. bancrofti thioredoxin is predicted to be a future hope in saving human-being from the threat of filariasis.

Intrinsic basis of thermostability of prolyl oligopeptidase from Pyrococcus furiosus.

Salt-bridges play a key role in the thermostability of proteins adapted in stress environments whose intrinsic basis remains to be understood. We find that the higher hydrophilicity of PfP than that of HuP is due to the charged but not the polar residues. The primary role of these residues is to enhance the salt-bridges and their ME. Unlike HuP, PfP has made many changes in its intrinsic property to strengthen the salt-bridge. First, the desolvation energy is reduced by directing the salt-bridge towards the surface. Second, it has made bridge-energy more favorable by recruiting energetically advantageous partners with high helix-propensity among the six possible salt-bridge pairs. Third, ME-residues that perform intricate interactions have increased their energy contribution by making major changes in their binary properties. The use of salt-bridge partners as ME-residues, and ME-residues' overlapping usage, predominant in helices, and energetically favorable substitution are some of the favorable features of PfP compared to HuP. These changes in PfP reduce the unfavorable, increase the favorable ME-energy. Thus, the per salt-bridge stability of PfP is greater than that of HuP. Further, unfavorable target ME-residues can be identified whose mutation can increase the stability of salt-bridge. The study applies to other similar systems.

Designing of a novel multi-epitope peptide based vaccine against Brugia malayi: An in silico approach.

In spite of the tremendous efforts of the World Health Organization, scientific and medical community to eradicate lymphatic filariasis (LF) within 2020, the disease is still taking a huge toll on mankind throughout the globe. The current therapeutic strategies and solution measures against this alarming condition are suffering from a number of limitations such as inadequate effectiveness of the drugs against the adult-stage parasites, low bioavailability, and emergence of resistance. Considering this situation, development of the new therapeutics are urgently needed to combat human LF, especially targeting the adult filarial nematodes. Brugia malayi, the causative parasite for the human brugian filariasis majorly found in the countries of the South-Asia. In this study, we have designed a vaccine candidate using B-cell and T-cell epitopes derived from the aspartic protease of B. malayi (BmASP-1) and found to display significant humoral and cell mediated immune responses using in-silico approaches. Protein-protein docking between the human Toll-like receptor 4 (TLR4) and the vaccine candidate helped us to predict the way of inductive signaling that leads to immune-response. Molecular dynamics (MD) simulation studies further confirmed the proper docking between the TLR4 and vaccine candidate. Moreover, in-silico cloning of the vaccine element within the expression vector was found useful to optimize the restriction sites as well as to determine the primer location. Taken together, the in-silico vaccine candidate depicted in this study promises could be a useful therapeutic option for treating LF and experimental validation of this study is expected to strengthen the candidature of the said vaccine in the future.

Anti-HIV potential of diarylpyrimidine derivatives as non-nucleoside reverse transcriptase inhibitors: design, synthesis, docking, TOPKAT analysis and molecular dynamics simulations.

In view of the low toxicity of NNRTIs in comparison to NRTIs, a new series of diarylpyrimidine derivatives has been designed as NNRTIs against HIV-1. In silico studies using DS 3.0 software have shown that these compounds behaved as NNRTIs while interacting at the allosteric site of HIV-RT. The designed compounds have shown promising docking results, which revealed that all compounds formed hydrogen bonds with Lys101, Lys103, Tyr181, Tyr318 and π- interactions with Tyr181, Tyr188, Phe227 and Trp229 amino acid residues located in the non-nucleoside inhibitor binding pocket (NNIBP) of HIV-RT protein. The intended molecules have shown high binding affinity with HIV-1 RT, analogous to standard drug molecule-etravirine. TOPKAT results confirmed that the designed compounds were found to be less toxic than the reference drug. Further, employing molecular dynamics simulations, the complexes of the best screened compound 6 and etravirine with the HIV-1 RT protein were analyzed by calculating the RMSD, RMSF, Rg, number of hydrogen bonds, principal components of the coordinates, molecular mechanics-Poisson-Boltzmann surface area-based binding free energy and their decomposition for different interactions. The analysis demonstrated the higher stability of compound 6 than the standard drug etravirine with HIV-1 RT. The interactions like hydrogen-bonding, van-der-Waals, electrostatic and the solvent accessible surface energy have favorable contributions to the complex stability. Thus, the shortlisted designed compound has great promise as a potential inhibitor against HIV-1 RT.

Anion Sensing by Novel Triarylboranes Containing Boraanthracene: DFT Functional Assessment, Selective Interactions, and Mechanism Demonstration.

Analytical methods often involve expensive instrumentation and tedious sample pretreatment for an analyte detection. Being toxic and detrimental to human health, sensing of cyanide (CN-), fluoride (F-), chloride (Cl-), bromide (Br-), nitrate (NO3 -), acetate (CH3COO-), and bisulfate (HSO4 -) is performed by a boron-based molecular receptor, N,N,N,3,5-pentamethyl-4-{2-thia-9-boratricyclo[8.4.0.03,8]tetradeca-1(10),3(8),4,6,11,13-hexaen-9-yl}anili-nium (1), and the three newly designed receptors from it. Thermodynamics, electronic structure, and photophysical properties are computed by employing density functional theory (DFT) and time-dependent density functional theory (TD-DFT) to explore selective sensing of these anions and its mechanism. Free-energy changes (ΔG) and binding energies (ΔE) suggest that among these anions, only binding of CN- and F- is thermodynamically feasible with a very strong binding affinity with the receptors. Boron atoms containing positive natural charges act as the electrophilic centers to bind the anions involving a 2p-2p orbital overlap resulting in charge transfer. In the receptor-analyte complexes with CN- and F-, fluorescence is quenched due to the intramolecular charge transfer transitions (π-π* transitions in the case of the receptors lead to fluorescence), internal conversion, and associated configurational changes. Among the six tested functionals, CAM-B3LYP/6-31G(d) is found to be the most accurate one. The designed receptors are better fluorescent probes for F- and CN-, demonstrating their importance for the practical utility.

PROPAB: Computation of Propensities and Other Properties from Segments of 3D structure of Proteins.

Residues in allelic positions, in the local segment of aligned sequences of proteins show wide variations. Here, we describe PROPAB that computes the propensity tables for helix, strand and coil types from multiple 3D structure files following ab initio statistical procedure. It also classifies them in range specific and chain specific manners. It further computes percentage composition and physicochemical properties along with residues propensities. It also prepares FASTA files for different segments (helix, strand and coil) in the exact order that they follow in the sequence. Representative analyses on orthologous (homologous across species) proteins demonstrate wide segmental variations of physicochemical properties. Such variations provide insights to relate the adaptation of these proteins in a given functional constraint under diverse environmental conditions. Thus, the program finds applications in the structural and evolutionary analysis of proteins. AVAILABILITY: PROPAB is freely available at http://sourceforge.net/projects/propab/for worldwide user.

POWAINDv1.0: A Program for Protein-Water Interactions Determination.

Protein is the most exposed biomolecule in the aqueous environment of the cell. Its structure maintains a delicate balance between the rigidity and the flexibility that imparts binding specificity to its substrate/ligand, etc. Intramolecular interactions of polar and non-polar groups of amino acid residues and intermolecular weak interactions between these groups and shell-waters may contribute to the overall stability of the tertiary structure. However, the question as to what are the dynamics of interactions of shell-water with respect to weak forces and atom-groups of protein (AGP), requires systematic investigations. In this end, we have developed a procedure POWAINDv1.0 that analyzes interactions of crystallographic shell-waters (CSH) in residues and AGP specific manner. The shell-water and AGP specific bridge-interactions are also extracted. Further, the program analyzes favorable and unfavorable nature of each interaction based on the actual and 75% of the sum of van der Waals (vdW) radii of interacting atoms. The EXCEL-outputs are useful in understanding the profile for AGP-CSH interactions and contribution of each component in AGP. Taken together, the program provides intricate details on CSHprotein interactions and finds application in the structural Bioinformatics.

Substitutional Analysis of Orthologous Protein Families Using BLOCKS.

Orthologous proteins, form due to divergence of parental sequence, perform similar function under different environmental and biological conditions. Amino acid changes at locus specific positions form hetero-pairs whose role in BLOCK evolution is yet to be understood. We involve eight protein BLOCKs of known divergence rate to gain insight into the role of hetero-pairs in evolution. Our procedure APBEST uses BLOCK-FASTA file to extract BLOCK specific evolutionary parameters such as dominantly used hetero-pair (D), usage of hetero-pairs (E), non-conservative to conservative substitution ratio (R), maximally-diverse residue (MDR), residue (RD) and class (CD) specific diversity. All these parameters show BLOCK specific variation. Conservative nature of D points towards restoration of function of BLOCK. While E sets the upper-limit of usage of hereto-pairs, strong correlation of R with divergence-rate indicates that the later is directly dependent on non-conservative substitutions. The observation that MDR, measure of positional diversity, occupy very limited positions in BLOCK indicates accommodation of diversity is positionally restricted. Overall, the study extract observed hetero-pair related quantitative and multi-parametric details of BLOCK, which finds application in evolutionary biology.

ADSBET2: Automated Determination of Salt-Bridge Energy-Terms version 2.

UNLABELLED: Component (bridge: ΔΔGbrd , background: ΔΔGprot , desolvation: ΔΔGdsolv ) and net (ΔΔGnet ) energy-terms of salt-bridge-structure (SBS) are auto-generated by the program ADSBET that makes use of general purpose Adaptive Poison Boltzmann Solver (APBS) method. While the procedure reports gross energy terms (Kcal Mol(-1) ), report on bond-multiplicity corrected normalized energyterms (Kcal Mol(-1) Bond(-1) ) along with their accessibility (ASA) in monomer, isolated-SBS (ISBS) and networked-SBS (NSBS) format would be very useful for statistical comparison among SBSs and understanding their location in protein structure. In this end, ADSBET2 potentially incorporates these features along with additional model for side-chain. Gross and normalized energy-terms are redirected in monomer, ISBS and NSBS format along with their ASA informations. It works on any number of SBSs for any number of structure files present in a database. Taken together, ADSBET2 has been suitable for statistical analyses of SBSs energetics and finds applications in protein engineering and structural bioinformatics. AVAILABILITY: ADSBET2 is freely available at http://sourceforge.net/projects/ADSBET2/ for all users.

PHYSICO2: an UNIX based standalone procedure for computation of physicochemical, window-dependent and substitution based evolutionary properties of protein sequences along with automated block preparation tool, version 2.

UNLABELLED: Automated genome sequencing procedure is enriching the sequence database very fast. To achieve a balance between the entry of sequences in the database and their analyses, efficient software is required. In this end PHYSICO2, compare to earlier PHYSICO and other public domain tools, is most efficient in that it i] extracts physicochemical, window-dependent and homologousposition-based-substitution (PWS) properties including positional and BLOCK-specific diversity and conservation, ii] provides users with optional-flexibility in setting relevant input-parameters, iii] helps users to prepare BLOCK-FASTA-file by the use of Automated Block Preparation Tool of the program, iv] performs fast, accurate and user-friendly analyses and v] redirects itemized outputs in excel format along with detailed methodology. The program package contains documentation describing application of methods. Overall the program acts as efficient PWS-analyzer and finds application in sequence-bioinformatics. AVAILABILITY: PHYSICO2: is freely available at http://sourceforge.net/projects/physico2/ along with its documentation at https://sourceforge.net/projects/physico2/files/Documentation.pdf/download for all users.

SBION2: Analyses of Salt Bridges from Multiple Structure Files, Version 2.

UNLABELLED: Specific electrostatics (i.e. salt-bridge) includes both local and non-local interactions that contribute to the overall stability of proteins. It has been shown that a salt-bridge could either be buried or exposed, networked or isolated, hydrogen-bonded or nonhydrogen bonded, in secondary-structure or in coil, formed by single or multiple bonds. Further it could also participates either in intra- or inter-dipole interactions with preference in orientation either for basic residue at N-terminal (orientation-I) or acidic residue at N-terminal (orientation-II). In this context SBION2 is unique in that it reports above mentioned binary items in excel format along with details on intra and inter-dipole interactions and orientations. These results are suitable for post run statistical analyses involving large datasets. Reports are also made on protein-protein interactions, intervening residue distances and general residue specific salt-bridge details. A ready to use compact supplementary table is also produced. The program runs in three alternative modes. Each mode works on any number of structure files with any number of chains at any given atomic distance of ion-pair. Thus SBION2 provides intricate details on salt-bridges and finds application in structural bioinformatics. AVAILABILITY: SBION2 is freely available at http://sourceforge.net/projects/sbion2/ for academic users.

SBION: A Program for Analyses of Salt-Bridges from Multiple Structure Files.

UNLABELLED: Salt-bridge and network salt-bridge are specific electrostatic interactions that contribute to the overall stability of proteins. In hierarchical protein folding model, these interactions play crucial role in nucleation process. The advent and growth of protein structure database and its availability in public domain made an urgent need for context dependent rapid analysis of salt-bridges. While these analyses on single protein is cumbersome and time-consuming, batch analyses need efficient software for rapid topological scan of a large number of protein for extracting details on (i) fraction of salt-bridge residues (acidic and basic). (ii) Chain specific intra-molecular salt-bridges, (iii) inter-molecular salt-bridges (protein-protein interactions) in all possible binary combinations (iv) network salt-bridges and (v) secondary structure distribution of salt-bridge residues. To the best of our knowledge, such efficient software is not available in public domain. At this juncture, we have developed a program i.e. SBION which can perform all the above mentioned computations for any number of protein with any number of chain at any given distance of ion-pair. It is highly efficient, fast, error-free and user friendly. Finally we would say that our SBION indeed possesses potential for applications in the field of structural and comparative bioinformatics studies. AVAILABILITY: SBION is freely available for non-commercial/academic institutions on formal request to the corresponding author (akbanerjee@biotech.buruniv.ac.in).

PHYSICO: An UNIX based Standalone Procedure for Computation of Individual and Group Properties of Protein Sequences.

UNLABELLED: In the genomic and proteomic era, efficient and automated analyses of sequence properties of protein have become an important task in bioinformatics. There are general public licensed (GPL) software tools to perform a part of the job. However, computations of mean properties of large number of orthologous sequences are not possible from the above mentioned GPL sets. Further, there is no GPL software or server which can calculate window dependent sequence properties for a large number of sequences in a single run. With a view to overcome above limitations, we have developed a standalone procedure i.e. PHYSICO, which performs various stages of computation in a single run based on the type of input provided either in RAW-FASTA or BLOCK-FASTA format and makes excel output for: a) Composition, Class composition, Mean molecular weight, Isoelectic point, Aliphatic index and GRAVY, b) column based compositions, variability and difference matrix, c) 25 kinds of window dependent sequence properties. The program is fast, efficient, error free and user friendly. Calculation of mean and standard deviation of homologous sequences sets, for comparison purpose when relevant, is another attribute of the program; a property seldom seen in existing GPL softwares. AVAILABILITY: PHYSICO is freely available for non-commercial/academic user in formal request to the corresponding author akbanerjee@biotech.buruniv.ac.in.