Liothyronine could block the programmed death-ligand 1 (PDL1) activity: an e-Pharmacophore modeling and virtual screening study.

PURPOSE: The interaction between PD-L1 on tumor cells and the programmed death 1 (PD1) on immune cells helps them to escape the immune system elimination. Therefore, developing therapeutic agents to block this interaction has garnered a lot of attention as a therapeutic approach. In the present study, we have tried to screen for an inhibitory compound to inhibit the interaction between the PD1/PD-L1 molecules. METHODS: In this regard, the structure of PD-L1 and its inhibitor were prepared and employed to generate an e-Pharmacophore model. A library of approved compounds was prepared and toxicity analysis using Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) predictor was performed. The built e-Pharmacophore model was validated and used to screen the prepared compound library. Ligand docking and binding energy calculation were performed on the screened ligands. RESULTS: A seven-feature e-Pharmacophore model was generated using the PD-L1 complex. All of the compounds within the library passed the ADMET criteria. Performing the virtual screening, only 79 compounds have survived the criteria to fit four pharmacophoric features. The compound with the highest binding energy was the liothyronine (T3). CONCLUSION: The ability of T3 in PD1/PD-L1 checkpoint blockade along with its potential in T4 reduction could be a desirable combination in cancer treatment. These abilities of T3 could be used to restore the ability of the immune system to eliminate tumor cells.

A new multi-epitope peptide vaccine induces immune responses and protection against Leishmania infantum in BALB/c mice.

Visceral leishmaniasis (VL) is a tropical and subtropical disease which is endemic in more than eighty countries around the world. Leishmania infantum is one of the main causative agents of VL disease. Currently, there is no approved-to-market vaccine for VL therapy. In this study, we evaluated cellular and humoral immune responses induced by our newly designed multi-epitope vaccine in BALB/c mice. Four antigenic proteins, including histone H1, sterol 24-c-methyltransferase (SMT), Leishmania-specific hypothetical protein (LiHy), and Leishmania-specific antigenic protein (LSAP) were chosen for the prediction of potential immunodominant epitopes. Moreover, to enhance vaccine immunogenicity, two toll-like receptors 4 (TLR4) agonists, resuscitation-promoting factors of Mycobacterium tuberculosis (RpfE and RpfB), were employed as the built-in adjuvants. Immunization with the designed multi-epitope vaccine elicited a robust Th1-type immune response, compared to other groups, as shown by increased levels of IL-2, IFN-γ, TNF-α, and IgG2a. Furthermore, a significant decrease was observed in Th-2-type-related cytokines such as IL-4 in immunized mice. The designed construct also induced a significant reduction in parasite load (p < 0.0001), conferring protection against L. infantum challenge. This study could be promising in gaining insight towards the potential of peptide epitope-based vaccines as effective protective approaches against Leishmania species.

Design of an engineered ACE2 as a novel therapeutics against COVID-19.

The interaction between the angiotensin-converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of the spike protein from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plays a pivotal role in virus entry into the host cells. Since recombinant ACE2 protein has been suggested as an anti-SARS-CoV-2 therapeutic agent, this study was conducted to design an ACE2 protein with more desirable properties. In this regard, the amino acids with central roles in enzymatic activity of the ACE2 were substituted. Moreover, saturation mutagenesis at the interaction interface between the ACE2 and RBD was performed to increase their interaction affinity. The best mutations to increase the structural and thermal stability of the ACE2 were also selected based on B factors and mutation effects. The obtained resulted revealed that the Arg273Gln and Thr445Gly mutation have drastically reduced the binding affinity of the angiotensin-II into the active site of ACE2. The Thr27Arg mutation was determined to be the most potent mutation to increase the binding affinity. The Asp427Arg mutation was done to decrease the flexibility of the region with high B factor. The Pro451Met mutation along with the Gly448Trp mutation was predicted to increase the thermodynamic stability and thermostability of the ACE2. The designed therapeutic ACE2 would have no enzymatic activity while it could bear stronger interaction with Spike glycoprotein of the SARS-CoV-2. Moreover, decreased in vivo enzymatic degradation would be anticipated due to increased thermostability. This engineered ACE2 could be exploited as a novel therapeutic agent against COVID-19 after necessary evaluations.

In silico analyses of Wilms׳ tumor protein to designing a novel multi-epitope DNA vaccine against cancer.

Predefined and pre-weighted objective criteria and essential role of Wilms׳ tumor wild type gene (WT1) for maintaining transformed features of cancer cells confirm the high potency of WT1 as a valuable cancer antigen. The antigen was at the top of the ranking among 75 representative cancer antigens. In the present study, an in silico approach was launched to characterized novel CTL epitopes and design a novel multi-epitope DNA vaccine to elicit a desirable immune response against cancers over expressing WT1. Forty-four novel epitopes were described. A multi-epitope construct was designed based on predicted epitopes which is 310 residues in length. The vaccine candidate designed here displays acceptable population coverage (>65%) in different ethnicities as well as high probability of eliciting WT1 antibodies which both are pertinent goals in the context of appropriate multi-epitope vaccines. Various in silico analyses indicate that final vaccine is a qualified immunotherapy candidate capable of eliciting both CD4+ and CD8+ T cell responses.

Production of a novel multi-epitope peptide vaccine for cancer immunotherapy in TC-1 tumor-bearing mice.

In our previous research, several bioinformatic strategies were utilized to design an efficient multi-epitope peptide vaccine (MEV) against cancer. The designed vaccine consists of Wilms tumor-1 (WT-1) and human papillomavirus (HPV) E7 cytotoxic T lymphocyte (CTL) epitopes, tetanus toxin fragment C (TTFrC) and HLA-DR epitope (PADRE) helper T lymphocyte (HTL) epitopes and heparin-binding hemagglutinin (HBHA) as an immunostimulatory adjuvant. All segments were fused together by suitable linkers. In the current study, we cloned and expressed the designed MEV in E. coli. We subsequently performed in vivo preventative and therapeutic assays to evaluate antitumor efficacy of the vaccine against the HPV-16 E7-expressing murine tumor cell line TC-1 as a model for cancer immunotherapy. The results showed that in preventive experiments, vaccination with MEV significantly augmented the IgG antibody titer and the percentage of tumor-free mice compared to control groups (PBS and E7). Moreover, in therapeutic experiments, vaccination with MEV led to a reduction in the number of metastatic nodules, lung weights and the ratio of lung weights to body weights compared to other groups. In sum, our epitope vaccine could efficiently induce preventive and therapeutic antitumor immunity in TC-1 tumor bearing mice.

In Silico Design and Evaluation of a Novel Therapeutic Agent Against the Spike Protein as a Novel Treatment Strategy for COVID-19 Treatment.

BACKGROUND: Coronavirus disease 2019 (COVID-19) is a viral respiratory disease that is associated with severe damage to other human organs. It causes by a novel coronavirus, and it is spreading all over the world. To date, there is some approved vaccine or therapeutic agent which could be effective against this disease. But their effectiveness against mutated strains is not studied completely. The spike glycoprotein on the surface of the coronaviruses gives the virus the ability to bind to host cell receptors and enter cells. Inhibition of attachment of these spikes can lead to virus neutralization by inhibiting viral entrance. AIMS: In this study, we tried to use the virus entrance strategy against itself by utilizing virus receptor (ACE-2) in order to design an engineered protein consisting of a human Fc antibody fragment and a part of ACE-2, which reacts with virus RBD, and we also evaluated this interaction by computational methods and in silico methods. Subsequently, we have designed a new protein structure to bind with this site and inhibit the virus from attaching to its cell receptor, mechanically or chemically. METHODS: Various in silico software, bioinformatics, and patent databases were used to retrieve the requested gene and protein sequences. The physicochemical properties and possibility of allergenicity were also examined. Three-dimensional structure prediction and molecular docking were also performed to develop the most suitable therapeutic protein. RESULTS: The designed protein consisted of a total of 256 amino acids with a molecular weight of 28984.62 and 5.92 as a theoretical isoelectric point. Instability and aliphatic index and grand average of hydropathicity are 49.99, 69.57 and -0.594, respectively. CONCLUSIONS: In silico studies can provide a good opportunity to study viral proteins and new drugs or compounds since they do not need direct exposure to infectious agents or equipped laboratories. The suggested therapeutic agent should be further characterized in vitro and in vivo.

Targeting Efficient Features of Urate Oxidase to Increase Its Solubility.

With the demand for mass production of protein drugs, solubility has become a serious issue. Extrinsic and intrinsic factors both affect this property. A homotetrameric cofactor-free urate oxidase (UOX) is not sufficiently soluble. To engineer UOX for optimum solubility, it is important to identify the most effective factor that influences solubility. The most effective feature to target for protein engineering was determined by measuring various solubility-related factors of UOX. A large library of homologous sequences was obtained from the databases. The data was reduced to six enzymes from different organisms. On the basis of various sequence- and structure-derived elements, the most and the least soluble enzymes were defined. To determine the best protein engineering target for modification, features of the most and least soluble enzymes were compared. Metabacillus fastidiosus UOX was the most soluble enzyme, while Agrobacterium globiformis UOX was the least soluble. According to the comparison-constant method, positive surface patches caused by arginine residue distribution are appropriate targets for modification. Two Arg to Ala mutations were introduced to the least soluble enzyme to test this hypothesis. These mutations significantly enhanced the mutant's solubility. While different algorithms produced conflicting results, it was difficult to determine which proteins were most and least soluble. Solubility prediction requires multiple algorithms based on these controversies. Protein surfaces should be investigated regionally rather than globally, and both sequence and structural data should be considered. Several other biotechnological products could be engineered using the data reduction and comparison-constant methods used in this study.

A potential in silico antibody-antigen based diagnostic test for precise identification of Acinetobacter baumannii.

Acinetobacter baumannii is a problematic nosocomial pathogen. The resistance to a wide range of antimicrobial agents, attributable to its biofilm phenotype, makes the treatment very difficult. Biofilm is a common feature of most pathogens. Biofilm associated proteins (Bap) are cellular surface components directly involved in biofilm formation process. The dearth of a fast precise diagnostic test and versatility of Bap sequences in A. baumannii were intuitions to design this study. In silico analysis is a reliable alternative to laborious experimental work in this connection. Databases were searched for an antigenic conserved region of Bap specific to A. baumannii. The region was selected based on alignments and propensity scales. Tertiary structure for this region was built and predicted B-cell epitopes were mapped on the surface of the built model. Our protein subunit was found to be a potential antigen, possessing several antigenic determinants, eliciting antibody. Hence this subunit could be used as a suitable agent for antibody-antigen based diagnostic test. This specific antigen can minimize laboratory errors in identification of A. baumannii and thus help clinicians to quick and precise diagnosis of the bacteria and initiatives to the treatment of the infection. Antigenicity of the region could also be explored for elicitation of antibody to protect the individuals exposed to A. baumannii.

Interaction of indole-3-acetic acid with horseradish peroxidase as a potential anticancer agent: from docking to molecular dynamics simulation.

The oxidation process, catalyzed by the peroxidase enzymes, occurs in all domains of life to detoxify the hydrogen peroxide toxicity. The most well-known, applicable and vastly studied member of the peroxidases family is horseradish peroxidase (HRP), especially the isoenzyme C (HRP C). HRP (primarily HRP C) is commercially available and applicable in biotechnology and diagnosis. Recently, a novel application of HRP has been introduced in cancer therapy as the combination of HRP with indole-3-acetic acid (IAA). The anticancer activity of HRP/IAA complex is through oxidation of IAA by HRP in hypoxic tumor condition, which leads to apoptosis and cancerous cell death. However, the molecular interaction of HRP/IAA has not been elucidated. Identifying the interaction of IAA with HRP would provide a better insight into its function and applications. In this study, molecular docking and molecular dynamics (MD) simulation were applied to determine the molecular interaction of the IAA/HRP complex. The docking study represented that IAA bound at the 'exposed' heme edge of the HRP enzyme, and the IAA entrance to the enzyme was situated at the carboxymethyl side-chain of the selected structure. Our computational results showed the HRP/IAA complex structure stability. While hydrogen bond formation with ARG38 and HIS42 stabilized the substrate, hydrophobic interactions with Phe68, Gly69, Leu138, Pro139, Pro141 and Phe179 contributed to IAA/HRP complex stability. The results can help to better understand peroxidase enzyme activity and would pave the way for future development of new therapeutics with improved anticancer efficacy.Communicated by Ramaswamy H. Sarma.

Theta-Defensins to Counter COVID-19 as Furin Inhibitors: In Silico Efficiency Prediction and Novel Compound Design.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was characterized as a pandemic by the World Health Organization (WHO) in Dec. 2019. SARS-CoV-2 binds to the cell membrane through spike proteins on its surface and infects the cell. Furin, a host-cell enzyme, possesses a binding site for the spike protein. Thus, molecules that block furin could potentially be a therapeutic solution. Defensins are antimicrobial peptides that can hypothetically inhibit furin because of their arginine-rich structure. Theta-defensins, a subclass of defensins, have attracted attention as drug candidates due to their small size, unique structure, and involvement in several defense mechanisms. Theta-defensins could be a potential treatment for COVID-19 through furin inhibition and an anti-inflammatory mechanism. Note that inflammatory events are a significant and deadly condition that could happen at the later stages of COVID-19 infection. Here, the potential of theta-defensins against SARS-CoV-2 infection was investigated through in silico approaches. Based on docking analysis results, theta-defensins can function as furin inhibitors. Additionally, a novel candidate peptide against COVID-19 with optimal properties regarding antigenicity, stability, electrostatic potential, and binding strength was proposed. Further in vitro/in vivo investigations could verify the efficiency of the designed novel peptide.

A unique antigen against SARS-CoV-2, Acinetobacter baumannii, and Pseudomonas aeruginosa.

The recent outbreak of COVID-19 has increased hospital admissions, which could elevate the risk of nosocomial infections, such as A. baumannii and P. aeruginosa infections. Although effective vaccines have been developed against SARS-CoV-2, no approved treatment option is still available against antimicrobial-resistant strains of A. baumannii and P. aeruginosa. In the current study, an all-in-one antigen was designed based on an innovative, state-of-the-art strategy. In this regard, experimentally validated linear epitopes of spike protein (SARS-CoV-2), OmpA (A. baumannii), and OprF (P. aeruginosa) were selected to be harbored by mature OmpA as a scaffold. The selected epitopes were used to replace the loops and turns of the barrel domain in OmpA; OprF311-341 replaced the most similar sequence within the OmpA, and three validated epitopes of OmpA were retained intact. The obtained antigen encompasses five antigenic peptides of spike protein, which are involved in SARS-CoV-2 pathogenicity. One of these epitopes, viz. QTQTNSPRRARSV could trigger antibodies preventing super-antigenic characteristics of spike and alleviating probable autoimmune responses. The designed antigen could raise antibodies neutralizing emerging variants of SARS-CoV-2 since at least two epitopes are consensus. In conclusion, the designed antigen is expected to raise protective antibodies against SARS-CoV-2, A. baumannii, and P. aeruginosa.

Impact of SNPs, off-targets, and passive permeability on efficacy of BCL6 degrading drugs assigned by virtual screening and 3D-QSAR approach.

B-cell lymphoma 6 (BCL6) regulates various genes and is reported to be overexpressed in lymphomas and other malignancies. Thus, BCL6 inhibition or its tagging for degradation would be an amenable therapeutic approach. A library of 2500 approved drugs was employed to find BCL6 inhibitory molecules via virtual screening. Moreover, the 3D core structure of 170 BCL6 inhibitors was used to build a 3D QSAR model and predict the biological activity. The SNP database was analyzed to study the impact on the destabilization of BCL6/drug interactions. Structural similarity search and molecular docking analyses were used to assess the interaction between possible off-targets and BCL6 inhibitors. The tendency of drugs for passive membrane permeability was also analyzed. Lifitegrast (DB11611) had favorable binding properties and biological activity compared to the BI-3802. Missense SNPs were located at the essential interaction sites of the BCL6. Structural similarity search resulted in five BTB-domain containing off-target proteins. BI-3802 and Lifitegrast had similar chemical behavior and binding properties against off-target candidates. More interestingly, the binding affinity of BI-3802 (against off-targets) was higher than Lifitegrast. Energetically, Lifitegrast was less favorable for passive membrane permeability. The interaction between BCL6 and BI-3802 is more prone to SNP-derived variations. On the other hand, higher nonspecific binding of BI-3802 to off-target proteins could bring about higher undesirable properties. It should also be noted that energetically less desirable passive membrane translocation of Lifitegrast would demand drug delivery vehicles. However, further empirical evaluation of Lifitegrast would unveil its true potential.

Pierce into Structural Changes of Interactions Between Mutated Spike Glycoproteins and ACE2 to Evaluate Its Potential Biological and Therapeutic Consequences.

The structural consequences of ongoing mutations on the SARS-CoV-2 spike-protein remains to be fully elucidated. These mutations could change the binding affinity between the virus and its target cell. Moreover, obtaining new mutations would also change the therapeutic efficacy of the designed drug candidates. To evaluate these consequences, 3D structure of a mutant spike protein was predicted and checked for stability, cavity sites, and residue depth. The docking analyses were performed between the 3D model of the mutated spike protein and the ACE2 protein and an engineered therapeutic ACE2 against COVID-19. The obtained results revealed that the N501Y substitution has altered the interaction orientation, augmented the number of interface bonds, and increased the affinity against the ACE2. On the other hand, the P681H mutation contributed to the increased cavity size and relatively higher residue depth. The binding affinity between the engineered therapeutic ACE2 and the mutant spike was significantly higher with a distinguished binding orientation. It could be concluded that the mutant spike protein increased the affinity, preserved the location, changed the orientation, and altered the interface amino acids of its interaction with both the ACE2 and its therapeutic engineered version. The obtained results corroborate the more aggressive nature of mutated SARS-CoV-2 due to their higher binding affinity. Moreover, designed ACe2-baased therapeutics would be still highly effective against covid-19, which could be the result of conserved nature of cellular ACE2. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10989-021-10346-1.

Protection against Acinetobacter baumannii infection via its functional deprivation of biofilm associated protein (Bap).

Acinetobacter baumannii, a major nosocomial pathogen, has remarkable capacity to acquire antimicrobial resistance attributable to its biofilm formation ability. Biofilm associated protein (Bap), a specific cell surface protein, is directly involved in biofilm formation by A. baumannii and plays a major role in bacterial infectious processes. In the present study we cloned, expressed and purified a 371 amino acid subunit of Bap. Mice were immunized using recombinant Bap subunit. They were then challenged with A. baumannii to evaluate the immunogenicity and protectivity of Bap subunit. Humoral immune response to Bap was determined by ELISA. Injection of Bap subunit resulted in high antibody titers. Decrease in bacterial cell counts of the immunized mice was evident 18 h after challenge. Reaction of antibodies against Bap with several strains suggests that not only immunodominant regions of Bap in A. baumannii strains are conserved but also have the same epitope presenting pattern in different strains. Immunodominant region of Bap possesses target sites for a protective humoral immune response to A. baumannii. This seems to be a conserved region erecting efficacy of Bap as an appropriate vaccine candidate.

SARS-CoV-2 Proteome Harbors Peptides Which Are Able to Trigger Autoimmunity Responses: Implications for Infection, Vaccination, and Population Coverage.

Autoimmune diseases (ADs) could occur due to infectious diseases and vaccination programs. Since millions of people are expected to be infected with SARS-CoV-2 and vaccinated against it, autoimmune consequences seem inevitable. Therefore, we have investigated the whole proteome of the SARS-CoV-2 for its ability to trigger ADs. In this regard, the entire proteome of the SARS-CoV-2 was chopped into more than 48000 peptides. The produced peptides were searched against the entire human proteome to find shared peptides with similar experimentally confirmed T-cell and B-cell epitopes. The obtained peptides were checked for their ability to bind to HLA molecules. The possible population coverage was calculated for the most potent peptides. The obtained results indicated that the SARS-CoV-2 and human proteomes share 23 peptides originated from ORF1ab polyprotein, nonstructural protein NS7a, Surface glycoprotein, and Envelope protein of SARS-CoV-2. Among these peptides, 21 peptides had experimentally confirmed equivalent epitopes. Amongst, only nine peptides were predicted to bind to HLAs with known global allele frequency data, and three peptides were able to bind to experimentally confirmed HLAs of equivalent epitopes. Given the HLAs which have already been reported to be associated with ADs, the ESGLKTIL, RYPANSIV, NVAITRAK, and RRARSVAS were determined to be the most harmful peptides of the SARS-CoV-2 proteome. It would be expected that the COVID-19 pandemic and the vaccination against this pathogen could significantly increase the ADs incidences, especially in populations harboring HLA-B*08:01, HLA-A*024:02, HLA-A*11:01 and HLA-B*27:05. The Southeast Asia, East Asia, and Oceania are at higher risk of AD development.

Hotspots for mutations in the SARS-CoV-2 spike glycoprotein: a correspondence analysis.

Spike glycoprotein (Sgp) is liable for binding of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to the host receptors. Since Sgp is the main target for vaccine and drug designing, elucidating its mutation pattern could help in this regard. This study is aimed at investigating the correspondence of specific residues to the SgpSARS-CoV-2 functionality by explorative interpretation of sequence alignments. Centrality analysis of the Sgp dissects the importance of these residues in the interaction network of the RBD-ACE2 (receptor-binding domain) complex and furin cleavage site. Correspondence of RBD to threonine500 and asparagine501 and furin cleavage site to glutamine675, glutamine677, threonine678, and alanine684 was observed; all residues are exactly located at the interaction interfaces. The harmonious location of residues dictates the RBD binding property and the flexibility, hydrophobicity, and accessibility of the furin cleavage site. These species-specific residues can be assumed as real targets of evolution, while other substitutions tend to support them. Moreover, all these residues are parts of experimentally identified epitopes. Therefore, their substitution may affect vaccine efficacy. Higher rate of RBD maintenance than furin cleavage site was predicted. The accumulation of substitutions reinforces the probability of the multi-host circulation of the virus and emphasizes the enduring evolutionary events.

In silico Approaches for the Design and Optimization of Interfering Peptides Against Protein-Protein Interactions.

Large contact surfaces of protein-protein interactions (PPIs) remain to be an ongoing issue in the discovery and design of small molecule modulators. Peptides are intrinsically capable of exploring larger surfaces, stable, and bioavailable, and therefore bear a high therapeutic value in the treatment of various diseases, including cancer, infectious diseases, and neurodegenerative diseases. Given these promising properties, a long way has been covered in the field of targeting PPIs via peptide design strategies. In silico tools have recently become an inevitable approach for the design and optimization of these interfering peptides. Various algorithms have been developed to scrutinize the PPI interfaces. Moreover, different databases and software tools have been created to predict the peptide structures and their interactions with target protein complexes. High-throughput screening of large peptide libraries against PPIs; "hotspot" identification; structure-based and off-structure approaches of peptide design; 3D peptide modeling; peptide optimization strategies like cyclization; and peptide binding energy evaluation are among the capabilities of in silico tools. In the present study, the most recent advances in the field of in silico approaches for the design of interfering peptides against PPIs will be reviewed. The future perspective of the field and its advantages and limitations will also be pinpointed.

In silico analysis of antibody triggering biofilm associated protein in Acinetobacter baumannii.

Acinetobacter baumannii surface protein, commonly known as biofilm associated protein (Bap), is involved in biofilm formation. A high propensity among the clinical isolates to form biofilm and a significant association of biofilms with multiple drug resistance has been demonstrated. Production of antibodies can be used for inhibition of biofilm and control of the diseases caused by A. baumannii. Large molecular mass of Bap justifies an approach to identifying A. baumannii effective antigens. It has a core domain of seven repeat modules A-G. With the large number of available biofilm gene sequences, bioinformatic tools are needed to identify the genes encoding the antigens. Proteins containing these tandem repeats of Bap domains have high propensities to attach to each other to form biofilm. We hypothesized that conserved and functional domains of tandem repeat could be identified with a search and alignment of the repeats for evaluation of antigenic determinants. Here we demonstrate the results of bioinformatics screening and gene scan of the gene sequence database of homolog sequences to identify conserved domains. Higher scoring hits were found in repeat modules mostly D, B, C and A, respectively. Upon the analysis four regions of highly structural and functional conserved regions from Bap sequence of A. baumannii were selected. 3D structure, antigenicity and solubility predictions revealed that these regions were appropriate candidates for antibody production.

Non-adaptive Evolution of Trimeric Autotransporters in Brucellaceae.

Brucella species are Gram-negative, facultative intracellular pathogens. They are the main cause of brucellosis, which has led to a global health burden. Adherence of the pathogen to the host cells is the first step in the infection process. The bacteria can adhere to various biotic and abiotic surfaces using their outer membrane proteins. Trimeric autotransporter adhesins (TAAs) are modular homotrimers of various length and domain complexity. They are a diverse, and widespread gene family constituting the type Vc secretion pathway. These adhesins have been established as virulence factors in Brucellaceae. To date, no comprehensive and exhaustive study has been performed on the trimeric autotransporter family in the genus. In the present study, various bioinformatics tools were used to provide a novel evolutionary insight into the sequence and structure of this protein family in Brucellaceae. To this end, a dataset of all trimeric autotransporters from the Brucella genomes was built. Analyses included but were not limited to sequence alignment, phylogenetic tree constructions, codon-based test for selection, clustering of the sequences, and structure (primary to quaternary) predictions. Batch analyzes of the dataset suggested the existence of a few structural domains within the whole population. BatA from the B. abortus 2308 genome was selected as a reference to describe the features of these structural domains. Furthermore, we examined the structural basis for the observed rigidity and resiliency of the protein structure through a molecular dynamics evaluation, which led us to deduce that the random drift results in the non-adaptive evolution of the trimeric autotransporter genes in the Brucella genus. Notably, the modifications have occurred across the genus without interference of gene transmission.

Trimeric autotransporter adhesins in Acinetobacter baumannii, coincidental evolution at work.

Trimeric autotransporter (TAA), also known as type Vc secretion system, is expressed by many strains of Acinetobacter baumannii, an opportunistic pathogen, which is responsible for nosocomial infections worldwide. TAAs, are modular homotrimeric virulence factors, containing a signal peptide, complex stalk, and conserved membrane anchoring domain. The evolutionary mechanisms underlying the evolvement of these adhesins are not clear. Here, we showed that TAA genes were laterally acquired and underwent gene duplication and recombination. The heterogeneity of TAA nucleotide sequences, GC content, codon usage, and the probability of recombination and duplication events were assessed by MEGA7. Given the heterogeneity of sequences, we used all-against-all BLAST for clustering the TAAs. The pattern of distribution of TAAs are highly scattered; GC content and codon usage for these genes are variable. Multiple events of lateral gene transfer from the early history of Acinetobacter and the occurrence of gene duplication, gene loss, and recombination after acquiring the alien genes may explain the scattered pattern of distribution of TAAs. Additionally, this gene is not present in many clinical isolates of A. baumannii, thus is not a single virulence factor attributing to the infection. The advantage of harboring such genes might be adopting to different environments by developing the biofilm communities. We suggested that TAA genes were laterally acquired in the environmental context and incidentally provided some benefits at the infection site. Thus, coincidental evolution theory may be better suited for describing the evolution of TAA genes in A. baumannii genomes.