In Silico Analysis of Putative Paralytic Shellfish Poisoning Toxins Export Proteins in Cyanobacteria

Paralytic shellfish poisoning toxins (PSTs) are a family of more than 30 natural alkaloids synthesized by dinoflagellates and cyanobacteria whose toxicity in animals is mediated by voltage-gated Na⁺ channel blocking. The export of PST analogues may be through SxtF and SxtM, two putative MATE (multidrug and toxic compound extrusion) family transporters encoded in PSTs biosynthetic gene cluster (sxt). sxtM is present in every sxt cluster analyzed; however, sxtF is only present in the Cylindrospermopsis-Raphidiopsis clade. These transporters are energetically coupled with an electrochemical gradient of proton (H⁺) or sodium (Na⁺) ions across membranes. Because the functional role of PSTs remains unknown and methods for genetic manipulation in PST-producing organisms have not yet been developed, protein structure analyses will allow us to understand their function. By analyzing the sxt cluster of eight PST-producing cyanobacteria, we found no correlation between the presence of sxtF or sxtM and a specific PSTs profile. Phylogenetic analyses of SxtF/M showed a high conservation of SxtF in the Cylindrospermopsis-Raphidiopsis clade, suggesting conserved substrate affinity. Two domains involved in Na⁺ and drug recognition from NorM proteins (MATE family) of Vibrio parahaemolyticus and V. cholerae are present in SxtF/M. The Na⁺ recognition domain was conserved in both SxtF/M, indicating that Na⁺ can maintain the role as a cation anti-transporter. Consensus motifs for toxin binding differed between SxtF and SxtM implying differential substrate binding. Through protein modeling and docking analysis, we found that there is no marked affinity between the recognition domain and a specific PST analogue. This agrees with our previous results of PST export in R. brookii D9, where we observed that the response to Na⁺ incubation was similar to different analogues. These results reassert the hypothesis regarding the involvement of Na⁺ in toxin export, as well as the motifs L³⁹⁸XGLQD⁴⁰³ (SxtM) and L³⁹⁰VGLRD³⁹⁵ (SxtF) in toxin recognition.     

In Silico Prediction and In Vivo Validation of Daphnia pulex Micrornas

Daphnia pulex, the crustacean with the first sequenced genome, is an important organism that has been widely used in ecological and toxicological research. MicroRNAs (miRNAs) are 21–25 nucleotide small non-coding RNAs that are involved in a myriad of physiological processes. In this research, we predicted 75 D. pulex miRNAs by sequence homology and secondary structure identification from the full genome sequence. Fourteen predicted miRNAs were selected for quantitative real time polymerase chain reaction (RT-PCR) validation. Out of these, eight (mir-8, mir-9, mir-12, mir-92, mir-100, mir-133, mir-153 and mir-283) were successfully amplified and validated. Next, expression levels were quantified at three different life stages (days 4, 8 and 12 of age) using U6 spliceosomal RNA as a reference gene. The expression of mir-8, mir-9, mir-12, mir-92 and mir-100 significantly differed across time suggesting these microRNAs might play a critical role during D. pulex development. This is the first study to identify and validate miRNAs in D. pulex, which is an important first step in further studies that evaluate their roles in development and response to environmental and ecological stimuli.     

In Vitro,In Silico and In Vivo Studies of Ursolic Acid as an Anti-Filarial Agent

As part of our drug discovery program for anti-filarial agents from Indian medicinal plants, leaves of Eucalyptus tereticornis were chemically investigated, which resulted in the isolation and characterization of an anti-filarial agent, ursolic acid (UA) as a major constituent. Antifilarial activity of UA against the human lymphatic filarial parasite Brugia malayi using in vitro and in vivo assays, and in silico docking search on glutathione-s-transferase (GST) parasitic enzyme were carried out. The UA was lethal to microfilariae (mf; LC₁₀₀: 50; IC₅₀: 8.84 µM) and female adult worms (LC₁₀₀: 100; IC50: 35.36 µM) as observed by motility assay; it exerted 86% inhibition in MTT reduction potential of the adult parasites. The selectivity index (SI) of UA for the parasites was found safe. This was supported by the molecular docking studies, which showed adequate docking (LibDock) scores for UA (−8.6) with respect to the standard antifilarial drugs, ivermectin (IVM −8.4) and diethylcarbamazine (DEC-C −4.6) on glutathione-s-transferase enzyme. Further, in silico pharmacokinetic and drug-likeness studies showed that UA possesses drug-like properties. Furthermore, UA was evaluated in vivo in B. malayi-M. coucha model (natural infection), which showed 54% macrofilaricidal activity, 56% female worm sterility and almost unchanged microfilaraemia maintained throughout observation period with no adverse effect on the host. Thus, in conclusion in vitro, in silico and in vivo results indicate that UA is a promising, inexpensive, widely available natural lead, which can be designed and developed into a macrofilaricidal drug. To the best of our knowledge this is the first ever report on the anti-filarial potential of UA from E. tereticornis, which is in full agreement with the Thomson Reuter''s ‘Metadrug’ tool screening predictions.     

In Silico Adoption of an Orphan Nuclear Receptor NR4A1

A 4.1μs molecular dynamics simulation of the NR4A1 (hNur77) apo-protein has been undertaken and a previously undetected druggable pocket has become apparent that is located remotely from the ‘traditional’ nuclear receptor ligand-binding site. A NR4A1/bis-indole ligand complex at this novel site has been found to be stable over 1 μs of simulation and to result in an interesting conformational transmission to a remote loop that has the capacity to communicate with a NBRE within a RXR-α/NR4A1 heterodimer. Several features of the simulations undertaken indicate how NR4A1 can be affected by alternate-site modulators.     

In Silico Target-Specific siRNA Design Based on Domain Transfer in Heterogeneous Data

RNA interference via exogenous small interference RNAs (siRNA) is a powerful tool in gene function study and disease treatment. Designing efficient and specific siRNA on target gene remains the key issue in RNAi. Although various in silico models have been proposed for rational siRNA design, most of them focus on the efficiencies of selected siRNAs, while limited effort has been made to improve their specificities targeted on specific mRNAs, which is related to reducing off-target effects (OTEs) in RNAi. In our study, we propose for the first time that the enhancement of target specificity of siRNA design can be achieved computationally by domain transfer in heterogeneous data sources from different siRNA targets. A transfer learning based method i.e., heterogeneous regression (HEGS) is presented for target-specific siRNA efficacy modeling and feature selection. Based on the model, (1) the target regression model can be built by extracting information from related data in other targets/experiments, thus increasing the target specificity in siRNA design with the help of information from siRNAs binding to other homologous genes, and (2) the potential features correlated to the current siRNA design can be identified even when there is lack of experimental validated siRNA affinity data on this target. In summary, our findings present useful instructions for a better target-specific siRNA design, with potential applications in genome-wide high-throughput screening of effective siRNA, and will provide further insights on the mechanism of RNAi.     

In Silico Screening of the Key Cellular Remodeling Targets in Chronic Atrial Fibrillation

Chronic atrial fibrillation (AF) is a complex disease with underlying changes in electrophysiology, calcium signaling and the structure of atrial myocytes. How these individual remodeling targets and their emergent interactions contribute to cell physiology in chronic AF is not well understood. To approach this problem, we performed in silico experiments in a computational model of the human atrial myocyte. The remodeled function of cellular components was based on a broad literature review of in vitro findings in chronic AF, and these were integrated into the model to define a cohort of virtual cells. Simulation results indicate that while the altered function of calcium and potassium ion channels alone causes a pronounced decrease in action potential duration, remodeling of intracellular calcium handling also has a substantial impact on the chronic AF phenotype. We additionally found that the reduction in amplitude of the calcium transient in chronic AF as compared to normal sinus rhythm is primarily due to the remodeling of calcium channel function, calcium handling and cellular geometry. Finally, we found that decreased electrical resistance of the membrane together with remodeled calcium handling synergistically decreased cellular excitability and the subsequent inducibility of repolarization abnormalities in the human atrial myocyte in chronic AF. We conclude that the presented results highlight the complexity of both intrinsic cellular interactions and emergent properties of human atrial myocytes in chronic AF. Therefore, reversing remodeling for a single remodeled component does little to restore the normal sinus rhythm phenotype. These findings may have important implications for developing novel therapeutic approaches for chronic AF.     

Semi-Automatic In Silico Gap Closure Enabled De Novo Assembly of Two Dehalobacter Genomes from Metagenomic Data

Typically, the assembly and closure of a complete bacterial genome requires substantial additional effort spent in a wet lab for gap resolution and genome polishing. Assembly is further confounded by subspecies polymorphism when starting from metagenome sequence data. In this paper, we describe an in silico gap-resolution strategy that can substantially improve assembly. This strategy resolves assembly gaps in scaffolds using pre-assembled contigs, followed by verification with read mapping. It is capable of resolving assembly gaps caused by repetitive elements and subspecies polymorphisms. Using this strategy, we realized the de novo assembly of the first two Dehalobacter genomes from the metagenomes of two anaerobic mixed microbial cultures capable of reductive dechlorination of chlorinated ethanes and chloroform. Only four additional PCR reactions were required even though the initial assembly with Newbler v. 2.5 produced 101 contigs within 9 scaffolds belonging to two Dehalobacter strains. By applying this strategy to the re-assembly of a recently published genome of Bacteroides, we demonstrate its potential utility for other sequencing projects, both metagenomic and genomic.     

Ih Tunes Theta/Gamma Oscillations and Cross-Frequency Coupling In an In Silico CA3 Model

channels are uniquely positioned to act as neuromodulatory control points for tuning hippocampal theta (4–12 Hz) and gamma (25 Hz) oscillations, oscillations which are thought to have importance for organization of information flow. contributes to neuronal membrane resonance and resting membrane potential, and is modulated by second messengers. We investigated oscillatory control using a multiscale computer model of hippocampal CA3, where each cell class (pyramidal, basket, and oriens-lacunosum moleculare cells), contained type-appropriate isoforms of . Our model demonstrated that modulation of pyramidal and basket allows tuning theta and gamma oscillation frequency and amplitude. Pyramidal also controlled cross-frequency coupling (CFC) and allowed shifting gamma generation towards particular phases of the theta cycle, effected via ''s ability to set pyramidal excitability. Our model predicts that in vivo neuromodulatory control of allows flexibly controlling CFC and the timing of gamma discharges at particular theta phases.     

In Silico Analysis of Antibiotic Resistance Genes in the Gut Microflora of Individuals from Diverse Geographies and Age-Groups

The spread of antibiotic resistance, originating from the rampant and unrestrictive use of antibiotics in humans and livestock over the past few decades has emerged as a global health problem. This problem has been further compounded by recent reports implicating the gut microbial communities to act as reservoirs of antibiotic resistance. We have profiled the presence of probable antibiotic resistance genes in the gut flora of 275 individuals from eight different nationalities. For this purpose, available metagenomic data sets corresponding to 275 gut microbiomes were analyzed. Sequence similarity searches of the genomic fragments constituting each of these metagenomes were performed against genes conferring resistance to around 240 antibiotics. Potential antibiotic resistance genes conferring resistance against 53 different antibiotics were detected in the human gut microflora analysed in this study. In addition to several geography/country-specific patterns, four distinct clusters of gut microbiomes, referred to as ‘Resistotypes’, exhibiting similarities in their antibiotic resistance profiles, were identified. Groups of antibiotics having similarities in their resistance patterns within each of these clusters were also detected. Apart from this, mobile multi-drug resistance gene operons were detected in certain gut microbiomes. The study highlighted an alarmingly high abundance of antibiotic resistance genes in two infant gut microbiomes. The results obtained in the present study presents a holistic ‘big picture’ on the spectra of antibiotic resistance within our gut microbiota across different geographies. Such insights may help in implementation of new regulations and stringency on the existing ones.     

In Silico Models for Dynamic Connected Cell Cultures Mimicking Hepatocyte-Endothelial Cell-Adipocyte Interaction Circle

The biochemistry of a system made up of three kinds of cell is virtually impossible to work out without the use of in silico models. Here, we deal with homeostatic balance phenomena from a metabolic point of view and we present a new computational model merging three single-cell models, already available from our research group: the first model reproduced the metabolic behaviour of a hepatocyte, the second one represented an endothelial cell, and the third one described an adipocyte. Multiple interconnections were created among these three models in order to mimic the main physiological interactions that are known for the examined cell phenotypes. The ultimate aim was to recreate the accomplishment of the homeostatic balance as it was observed for an in vitro connected three-culture system concerning glucose and lipid metabolism in the presence of the medium flow. The whole model was based on a modular approach and on a set of nonlinear differential equations implemented in Simulink, applying Michaelis-Menten kinetic laws and some energy balance considerations to the studied metabolic pathways. Our in silico model was then validated against experimental datasets coming from literature about the cited in vitro model. The agreement between simulated and experimental results was good and the behaviour of the connected culture system was reproduced through an adequate parameter evaluation. The developed model may help other researchers to investigate further about integrated metabolism and the regulation mechanisms underlying the physiological homeostasis.     

In Vitro HIV-1 Evolution in Response to Triple Reverse Transcriptase Inhibitors & In Silico Phenotypic Analysis

Background

Effectiveness of ART regimens strongly depends upon complex interactions between the selective pressure of drugs and the evolution of mutations that allow or restrict drug resistance.

Methods

Four clinical isolates from NRTI-exposed, NNRTI-naive subjects were passaged in increasing concentrations of NVP in combination with 1 µM 3 TC and 2 µM ADV to assess selective pressures of multi-drug treatment. A novel parameter inference procedure, based on a stochastic viral growth model, was used to estimate phenotypic resistance and fitness from in vitro combination passage experiments.

Results

Newly developed mathematical methods estimated key phenotypic parameters of mutations arising through selective pressure exerted by 3 TC and NVP. Concentrations of 1 µM 3 TC maintained the M184V mutation, which was associated with intrinsic fitness deficits. Increasing NVP concentrations selected major NNRTI resistance mutations. The evolutionary pathway of NVP resistance was highly dependent on the viral genetic background, epistasis as well as stochasticity. Parameter estimation indicated that the previously unrecognized mutation L228Q was associated with NVP resistance in some isolates.

Conclusion

Serial passage of viruses in the presence of multiple drugs may resemble the selection of mutations observed among treated individuals and populations in vivo and indicate evolutionary preferences and restrictions. Phenotypic resistance estimated here “in silico” from in vitro passage experiments agreed well with previous knowledge, suggesting that the unique combination of “wet-” and “dry-lab” experimentation may improve our understanding of HIV-1 resistance evolution in the future.     

In Silico Analysis of Functional Single Nucleotide Polymorphisms in the Human TRIM22 Gene

Tripartite motif protein 22 (TRIM22) is an evolutionarily ancient protein that plays an integral role in the host innate immune response to viruses. The antiviral TRIM22 protein has been shown to inhibit the replication of a number of viruses, including HIV-1, hepatitis B, and influenza A. TRIM22 expression has also been associated with multiple sclerosis, cancer, and autoimmune disease. In this study, multiple in silico computational methods were used to identify non-synonymous or amino acid-changing SNPs (nsSNP) that are deleterious to TRIM22 structure and/or function. A sequence homology-based approach was adopted for screening nsSNPs in TRIM22, including six different in silico prediction algorithms and evolutionary conservation data from the ConSurf web server. In total, 14 high-risk nsSNPs were identified in TRIM22, most of which are located in a protein interaction module called the B30.2 domain. Additionally, 9 of the top high-risk nsSNPs altered the putative structure of TRIM22''s B30.2 domain, particularly in the surface-exposed v2 and v3 regions. These same regions are critical for retroviral restriction by the closely-related TRIM5α protein. A number of putative structural and functional residues, including several sites that undergo post-translational modification, were also identified in TRIM22. This study is the first extensive in silico analysis of the highly polymorphic TRIM22 gene and will be a valuable resource for future targeted mechanistic and population-based studies.     

In Silico Repositioning-Chemogenomics Strategy Identifies New Drugs with Potential Activity against Multiple Life Stages of Schistosoma mansoni

Morbidity and mortality caused by schistosomiasis are serious public health problems in developing countries. Because praziquantel is the only drug in therapeutic use, the risk of drug resistance is a concern. In the search for new schistosomicidal drugs, we performed a target-based chemogenomics screen of a dataset of 2,114 proteins to identify drugs that are approved for clinical use in humans that may be active against multiple life stages of Schistosoma mansoni. Each of these proteins was treated as a potential drug target, and its amino acid sequence was used to interrogate three databases: Therapeutic Target Database (TTD), DrugBank and STITCH. Predicted drug-target interactions were refined using a combination of approaches, including pairwise alignment, conservation state of functional regions and chemical space analysis. To validate our strategy, several drugs previously shown to be active against Schistosoma species were correctly predicted, such as clonazepam, auranofin, nifedipine, and artesunate. We were also able to identify 115 drugs that have not yet been experimentally tested against schistosomes and that require further assessment. Some examples are aprindine, gentamicin, clotrimazole, tetrabenazine, griseofulvin, and cinnarizine. In conclusion, we have developed a systematic and focused computer-aided approach to propose approved drugs that may warrant testing and/or serve as lead compounds for the design of new drugs against schistosomes.     

In Silico,In Vitro and In Vivo Analysis of Binding Affinity between N and C-Domains of Clostridium perfringens Alpha Toxin

Clostridium perfringens alpha toxin/phospholipase C (CP-PLC) is one of the most potent bacterial toxins known to cause soft tissue infections like gas gangrene in humans and animals. It is the first bacterial toxin demonstrated to be an enzyme with phospholipase, sphingomyelinase and lecithinase activities. The toxin is comprised of an enzymatic N-domain and a binding C-domain interconnected by a flexible linker. The N-domain alone is non-toxic to mammalian cells, but incubation with C-domain restores the toxicity, the mechanism of which is still not elucidated. The objectives of the current study were to investigate the formation of a stable N and C-domain complex, to determine possible interactions between the two domains in silico and to characterize the in vitro and in vivo correlates of the interaction. To establish the existence of a stable N and C-domain hybrid, in vitro pull down assay and dot-Far Western blotting assays were employed, where it was clearly revealed that the two domains bound to each other to form an intermediate. Using bioinformatics tools like MetaPPISP, PatchDock and FireDock, we predicted that the two domains may interact with each other through electrostatic interactions between at least six pairs of amino acids. This N and C-domains interacted with each other in 1:1 ratio and the hybrid lysed mouse erythrocytes in a slower kinetics when compared with wild type native Cp-PLC. BALB/c mice when challenged with N and C-domain hybrid demonstrated severe myonecrosis at the site of injection while no death was observed. Our results provide further insight into better understanding the mechanism for the toxicity of Cp-PLC N and C-domain mixture.     

In Silico Investigations of the Anti-Catabolic Effects of Pamidronate and Denosumab on Multiple Myeloma-Induced Bone Disease

It is unclear whether the new anti-catabolic agent denosumab represents a viable alternative to the widely used anti-catabolic agent pamidronate in the treatment of Multiple Myeloma (MM)-induced bone disease. This lack of clarity primarily stems from the lack of sufficient clinical investigations, which are costly and time consuming. However, in silico investigations require less time and expense, suggesting that they may be a useful complement to traditional clinical investigations. In this paper, we aim to (i) develop integrated computational models that are suitable for investigating the effects of pamidronate and denosumab on MM-induced bone disease and (ii) evaluate the responses to pamidronate and denosumab treatments using these integrated models. To achieve these goals, pharmacokinetic models of pamidronate and denosumab are first developed and then calibrated and validated using different clinical datasets. Next, the integrated computational models are developed by incorporating the simulated transient concentrations of pamidronate and denosumab and simulations of their actions on the MM-bone compartment into the previously proposed MM-bone model. These integrated models are further calibrated and validated by different clinical datasets so that they are suitable to be applied to investigate the responses to the pamidronate and denosumab treatments. Finally, these responses are evaluated by quantifying the bone volume, bone turnover, and MM-cell density. This evaluation identifies four denosumab regimes that potentially produce an overall improved bone-related response compared with the recommended pamidronate regime. This in silico investigation supports the idea that denosumab represents an appropriate alternative to pamidronate in the treatment of MM-induced bone disease.     

In Silico Identification of Multi-target Anti-SARS-CoV-2 Peptides from Quinoa Seed Proteins

International Journal of Peptide Research and Therapeutics - Peptides are promising antagonists against severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2). To expedite drug...     

In Silico Molecular Comparisons of C. elegans and Mammalian Pharmacology Identify Distinct Targets That Regulate Feeding

Phenotypic screens can identify molecules that are at once penetrant and active on the integrated circuitry of a whole cell or organism. These advantages are offset by the need to identify the targets underlying the phenotypes. Additionally, logistical considerations limit screening for certain physiological and behavioral phenotypes to organisms such as zebrafish and C. elegans. This further raises the challenge of elucidating whether compound-target relationships found in model organisms are preserved in humans. To address these challenges we searched for compounds that affect feeding behavior in C. elegans and sought to identify their molecular mechanisms of action. Here, we applied predictive chemoinformatics to small molecules previously identified in a C. elegans phenotypic screen likely to be enriched for feeding regulatory compounds. Based on the predictions, 16 of these compounds were tested in vitro against 20 mammalian targets. Of these, nine were active, with affinities ranging from 9 nM to 10 µM. Four of these nine compounds were found to alter feeding. We then verified the in vitro findings in vivo through genetic knockdowns, the use of previously characterized compounds with high affinity for the four targets, and chemical genetic epistasis, which is the effect of combined chemical and genetic perturbations on a phenotype relative to that of each perturbation in isolation. Our findings reveal four previously unrecognized pathways that regulate feeding in C. elegans with strong parallels in mammals. Together, our study addresses three inherent challenges in phenotypic screening: the identification of the molecular targets from a phenotypic screen, the confirmation of the in vivo relevance of these targets, and the evolutionary conservation and relevance of these targets to their human orthologs.     

In Silico Insights into Protein-Protein Interactions and Folding Dynamics of the Saposin-Like Domain of Solanum tuberosum Aspartic Protease

The plant-specific insert is an approximately 100-residue domain found exclusively within the C-terminal lobe of some plant aspartic proteases. Structurally, this domain is a member of the saposin-like protein family, and is involved in plant pathogen defense as well as vacuolar targeting of the parent protease molecule. Similar to other members of the saposin-like protein family, most notably saposins A and C, the recently resolved crystal structure of potato (Solanum tuberosum) plant-specific insert has been shown to exist in a substrate-bound open conformation in which the plant-specific insert oligomerizes to form homodimers. In addition to the open structure, a closed conformation also exists having the classic saposin fold of the saposin-like protein family as observed in the crystal structure of barley (Hordeum vulgare L.) plant-specific insert. In the present study, the mechanisms of tertiary and quaternary conformation changes of potato plant-specific insert were investigated in silico as a function of pH. Umbrella sampling and determination of the free energy change of dissociation of the plant-specific insert homodimer revealed that increasing the pH of the system to near physiological levels reduced the free energy barrier to dissociation. Furthermore, principal component analysis was used to characterize conformational changes at both acidic and neutral pH. The results indicated that the plant-specific insert may adopt a tertiary structure similar to the characteristic saposin fold and suggest a potential new structural motif among saposin-like proteins. To our knowledge, this acidified PSI structure presents the first example of an alternative saposin-fold motif for any member of the large and diverse SAPLIP family.