Single-cell multi-omics analysis identifies context-specific gene regulatory gates and mechanisms.

There is a growing interest in inferring context specific gene regulatory networks from single-cell RNA sequencing (scRNA-seq) data. This involves identifying the regulatory relationships between transcription factors (TFs) and genes in individual cells, and then characterizing these relationships at the level of specific cell types or cell states. In this study, we introduce scGATE (single-cell gene regulatory gate) as a novel computational tool for inferring TF-gene interaction networks and reconstructing Boolean logic gates involving regulatory TFs using scRNA-seq data. In contrast to current Boolean models, scGATE eliminates the need for individual formulations and likelihood calculations for each Boolean rule (e.g. AND, OR, XOR). By employing a Bayesian framework, scGATE infers the Boolean rule after fitting the model to the data, resulting in significant reductions in time-complexities for logic-based studies. We have applied assay for transposase-accessible chromatin with sequencing (scATAC-seq) data and TF DNA binding motifs to filter out non-relevant TFs in gene regulations. By integrating single-cell clustering with these external cues, scGATE is able to infer context specific networks. The performance of scGATE is evaluated using synthetic and real single-cell multi-omics data from mouse tissues and human blood, demonstrating its superiority over existing tools for reconstructing TF-gene networks. Additionally, scGATE provides a flexible framework for understanding the complex combinatorial and cooperative relationships among TFs regulating target genes by inferring Boolean logic gates among them.

Prospects and challenges of cancer systems medicine: from genes to disease networks.

It is becoming evident that holistic perspectives toward cancer are crucial in deciphering the overwhelming complexity of tumors. Single-layer analysis of genome-wide data has greatly contributed to our understanding of cellular systems and their perturbations. However, fundamental gaps in our knowledge persist and hamper the design of effective interventions. It is becoming more apparent than ever, that cancer should not only be viewed as a disease of the genome but as a disease of the cellular system. Integrative multilayer approaches are emerging as vigorous assets in our endeavors to achieve systemic views on cancer biology. Herein, we provide a comprehensive review of the approaches, methods and technologies that can serve to achieve systemic perspectives of cancer. We start with genome-wide single-layer approaches of omics analyses of cellular systems and move on to multilayer integrative approaches in which in-depth descriptions of proteogenomics and network-based data analysis are provided. Proteogenomics is a remarkable example of how the integration of multiple levels of information can reduce our blind spots and increase the accuracy and reliability of our interpretations and network-based data analysis is a major approach for data interpretation and a robust scaffold for data integration and modeling. Overall, this review aims to increase cross-field awareness of the approaches and challenges regarding the omics-based study of cancer and to facilitate the necessary shift toward holistic approaches.

wpLogicNet: logic gate and structure inference in gene regulatory networks.

MOTIVATION: The gene regulatory process resembles a logic system in which a target gene is regulated by a logic gate among its regulators. While various computational techniques are developed for a gene regulatory network (GRN) reconstruction, the study of logical relationships has received little attention. Here, we propose a novel tool called wpLogicNet that simultaneously infers both the directed GRN structures and logic gates among genes or transcription factors (TFs) that regulate their target genes, based on continuous steady-state gene expressions. RESULTS: wpLogicNet proposes a framework to infer the logic gates among any number of regulators, with a low time-complexity. This distinguishes wpLogicNet from the existing logic-based models that are limited to inferring the gate between two genes or TFs. Our method applies a Bayesian mixture model to estimate the likelihood of the target gene profile and to infer the logic gate a posteriori. Furthermore, in structure-aware mode, wpLogicNet reconstructs the logic gates in TF-gene or gene-gene interaction networks with known structures. The predicted logic gates are validated on simulated datasets of TF-gene interaction networks from Escherichia coli. For the directed-edge inference, the method is validated on datasets from E.coli and DREAM project. The results show that compared to other well-known methods, wpLogicNet is more precise in reconstructing the network and logical relationships among genes. AVAILABILITY AND IMPLEMENTATION: The datasets and R package of wpLogicNet are available in the github repository, https://github.com/CompBioIPM/wpLogicNet. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Contribution of the intra-hippocampal orexin system in the regulation of restraint stress response to pain-related behaviors in the formalin test.

Stress-induced antinociception (SIA) is due to the activation of several neural pathways and neurotransmitters that often suppress pain perception. Studies have shown that the orexin neuropeptide system is essential in pain modulation. Therefore, this study aimed to investigate the role of orexinergic receptors in the hippocampal CA1 region in modulating SIA response during the formalin test as an animal model of inflammatory pain. The orexin-1 receptor (OX1r) antagonist, SB334867, at 1, 3, 10, and 30 nmol or TCS OX2 29 as an orexin-2 receptor (OX2r) antagonist at the same doses were microinjected into the CA1 region in rats. Five minutes later, rats were exposed to restraint stress (RS) for 3 h, and pain-related behaviors were monitored in 5-min blocks for the 60-min test period in the formalin test. Results showed that applying RS for 3 h reduced pain responses in the early and late phases of the formalin test. The main findings showed that intra-CA1 injection of orexin receptor antagonists reduced the antinociception caused by stress in both phases of the formalin test. In addition, the contribution of OX2r in mediating the antinociceptive effect of stress was more prominent than that of OX1r in the early phase of the formalin test. However, in the late phase, both receptors worked similarly. Accordingly, the orexin system and its two receptors in the CA1 region of the hippocampus regulate SIA response to this animal model of pain in formalin test.

Design, construction and in vivo functional assessment of a hinge truncated sFLT01.

Gene therapy for the treatment of ocular neovascularization has reached clinical trial phases. The AAV2-sFLT01 construct was already evaluated in a phase 1 open-label trial administered intravitreally to patients with advanced neovascular age-related macular degeneration. SFLT01 protein functions by binding to VEGF and PlGF molecules and inhibiting their activities simultaneously. It consists of human VEGFR1/Flt-1 (hVEGFR1), a polyglycine linker, and the Fc region of human IgG1. The IgG1 upper hinge region of the sFLT01 molecule makes it vulnerable to radical attacks and prone to causing immune reactions. This study pursued two goals: (i) minimizing the immunogenicity and vulnerability of the molecule by designing a truncated molecule called htsFLT01 (hinge truncated sFLT01) that lacked the IgG1 upper hinge and lacked 2 amino acids from the core hinge region; and (ii) investigating the structural and functional properties of the aforesaid chimeric molecule at different levels (in silico, in vitro, and in vivo). Molecular dynamics simulations and molecular mechanics energies combined with Poisson-Boltzmann and surface area continuum solvation calculations revealed comparable free energy of binding and binding affinity for sFLT01 and htsFLT01 to their cognate ligands. Conditioned media from human retinal pigment epithelial (hRPE) cells that expressed htsFLT01 significantly reduced tube formation in HUVECs. The AAV2-htsFLT01 virus suppressed vascular development in the eyes of newborn mice. The htsFLT01 gene construct is a novel anti-angiogenic tool with promising improvements compared to existing treatments.

Identification of UBA1 as the causative gene of an X-linked non-Kennedy spinal-bulbar muscular atrophy.

BACKGROUND AND PURPOSE: Spinal-bulbar muscular atrophy (SBMA) (Kennedy's disease) is a motor neuron disease. Kennedy's disease is nearly exclusively caused by mutations in the androgen receptor encoding gene (AR). The results of studies aimed at identification of the genetic cause of a disease that best approximates SBMA in a pedigree (four patients) without mutations in AR are reported. METHODS: Clinical investigations included thorough neurological and non-neurological examinations and testing. Genetic analysis was performed by exome sequencing using standard protocols. UBA1 mutations were modeled on the crystal structure of UBA1. RESULTS: The clinical features of the patients are described in detail. A missense mutation in UBA1 (c.T1499C; p.Ile500Thr) was identified as the probable cause of the non-Kennedy SBMA in the pedigree. Like AR, UBA1 is positioned on chromosome X. UBA1 is a highly conserved gene. It encodes ubiquitin-like modifier activating enzyme 1 (UBA1) which is the major E1 enzyme of the ubiquitin-proteasome system. Interestingly, UBA1 mutations can also cause infantile-onset X-linked spinal muscular atrophy (XL-SMA). The mutation identified here and the XL-SMA causative mutations were shown to affect amino acids positioned in the vicinity of UBA1's ATP binding site and to cause structural changes. CONCLUSION: UBA1 was identified as a novel SBMA causative gene. The gene affects protein homeostasis which is one of most important components of the pathology of neurodegeneration. The contribution of this same gene to the etiology of XL-SMA is discussed.

Network analysis and the impact of Aflibercept on specific mediators of angiogenesis in HUVEC cells.

Angiogenesis, inflammation and endothelial cells' migration and proliferation exert fundamental roles in different diseases. However, more studies are needed to identify key proteins and pathways involved in these processes. Aflibercept has received the approval of the US Food and Drug Administration (FDA) for the treatment of wet AMD and colorectal cancer. Moreover, the effect of Aflibercept on VEGFR2 downstream signalling pathways has not been investigated yet. Here, we integrated text mining data, protein-protein interaction networks and multi-experiment microarray data to specify candidate genes that are involved in VEGFA/VEGFR2 signalling pathways. Network analysis of candidate genes determined the importance of the nominated genes via different centrality parameters. Thereupon, several genes-with the highest centrality indexes-were recruited to investigate the impact of Aflibercept on their expression pattern in HUVEC cells. Real-time PCR was performed, and relative expression of the specific genes revealed that Aflibercept modulated angiogenic process by VEGF/PI3KA/AKT/mTOR axis, invasion by MMP14/MMP9 axis and inflammation-related angiogenesis by IL-6-STAT3 axis. Data showed Aflibercept simultaneously affected these processes and determined the nominated axes that had been affected by the drug. Furthermore, integrating the results of Aflibercept on expression of candidate genes with the current network analysis suggested that resistance against the Aflibercept effect is a plausible process in HUVEC cells.

Optogenetic control of neural differentiation in Opto-mGluR6 engineered retinal pigment epithelial cell line and mesenchymal stem cells.

In retinal degenerative disorders, when neural retinal cells are damaged, cell transplantation is one of the most promising therapeutic approaches. Optogenetic technology plays an essential role in the neural differentiation of stem cells via membrane depolarization. This study explored the efficacy of blue light stimulation in neuroretinal differentiation of Opto-mGluR6-engineered mouse retinal pigment epithelium (mRPE) and bone marrow mesenchymal stem cells (BMSCs). mRPE and BMSCs were selected for optogenetic study due to their capability to differentiate into retinal-specific neurons. BMSCs were isolated and phenotypically characterized by the expression of mesenchymal stem cell-specific markers, CD44 (99%) and CD105 (98.8%). mRPE culture identity was confirmed by expression of RPE-specific marker, RPE65, and epithelial cell marker, ZO-1. mRPE cells and BMSCs were transduced with AAV-MCS-IRES-EGFP-Opto-mGluR6 viral vector and stimulated for 5 days with blue light (470 nm). RNA and protein expression of Opto-mGluR6 were verified. Optogenetic stimulation-induced elevated intracellular Ca2+ levels in mRPE- and BMS-treated cells. Significant increase in cell growth rate and G1/S phase transition were detected in mRPE- and BMSCs-treated cultures. Pou4f1, Dlx2, Eomes, Barlh2, Neurod2, Neurod6, Rorb, Rxrg, Nr2f2, Ascl1, Hes5, and Sox8 were overexpressed in treated BMSCs and Barlh2, Rorb, and Sox8 were overexpressed in treated mRPE cells. Expression of Rho, Thy1, OPN1MW, Recoverin, and CRABP, as retinal-specific neuron markers, in mRPE and BMS cell cultures were demonstrated. Differentiation of ganglion, amacrine, photoreceptor cells, and bipolar and Muller precursors were determined in BMSCs-treated culture and were compared with mRPE. mRPE cells represented more abundant terminal Muller glial differentiation compared with BMSCs. Our results also demonstrated that optical stimulation increased the intracellular Ca2+ level and proliferation and differentiation of Opto-mGluR6-engineered BMSCs. It seems that optogenetic stimulation of mRPE- and BMSCs-engineered cells would be a potential therapeutic approach for retinal degenerative disorders.

Population-level heterogeneity as a reflection of mixed strategy: A computational perspective on the Crabtree effect.

Is population-level heterogeneity a reflection of distinct subpopulations, exhibiting different metabolic functions, or dynamic metabolism of individuals within the population? This fascinating question has remained a subject of great interest in studying metabolic specialization in microorganisms. The Crabtree effect - i.e., the ability of some microorganisms to switch from respiration to fermentation in the presence of oxygen - is an appropriate case study to address the aforementioned question. Game-theoretical approaches have been routinely used to examine and explain the way a microorganism, such as yeast, would switch between the two ATP-producing pathways, i.e., respiration and fermentation. Here we attempt to explain the switch between respiration and fermentation in yeast by constructing a simple metabolic switch. We then utilize an individual-based model, in which each individual is equipped with all the relevant chemical reactions, to see how cells equipped with such metabolic switch would behave in different conditions. We further investigate our proposed metabolic switch using the game-theoretical approach. Based on this model, we postulate that the population-level metabolic heterogeneity in microorganisms can simply arise from individuals utilizing a mixed strategy.

Orexin-2 receptor antagonism in the cornu ammonis 1 region of hippocampus prevented the antinociceptive responses induced by chemical stimulation of the lateral hypothalamus in the animal model of persistent pain.

Orexins are excitatory neuropeptides, mainly produced by neurons located in the lateral hypothalamus, which project to many brain areas. The orexinergic system plays a fundamental role in arousal, sleep/wakefulness, feeding, energy homeostasis, motivation, reward, stress and pain modulation. As a prominent part of the limbic system, the hippocampus has been involved in formalin-induced nociception modulation. Moreover, hippocampus regions express both orexin-1 (OX1) and orexin-2 (OX2) receptors. The present study investigated the role of OX2 receptors (OX2R) within the cornu ammonis 1 (CA1) region of the hippocampus in the mediation of lateral hypothalamus-induced antinociception. Fifty-three male Wistar rats were unilaterally implanted with two separate cannulae into the lateral hypothalamus and CA1. Animals were pretreated with intra-CA1 TCS OX2 29 as an OX2R antagonist before intra-lateral hypothalamus administration of carbachol (250 nM) as a muscarinic agonist for chemical stimulation of orexinergic neurons. Formalin test was used as an animal model of persistent pain, following intra-lateral hypothalamus carbachol microinjection. Results showed that the chemical stimulation of the lateral hypothalamus significantly attenuated formalin-evoked nociceptive behaviors during both phases of the formalin test, and administration of TCS OX2 29 into the CA1 blocked these antinociceptive responses in both phases, especially in the late phase. These findings suggest that OX2 receptors in the CA1 partially mediate the lateral hypothalamus-induced antinociceptive responses in persistent inflammatory pain.

Inter-protein residue covariation information unravels physically interacting protein dimers.

BACKGROUND: Predicting physical interaction between proteins is one of the greatest challenges in computational biology. There are considerable various protein interactions and a huge number of protein sequences and synthetic peptides with unknown interacting counterparts. Most of co-evolutionary methods discover a combination of physical interplays and functional associations. However, there are only a handful of approaches which specifically infer physical interactions. Hybrid co-evolutionary methods exploit inter-protein residue coevolution to unravel specific physical interacting proteins. In this study, we introduce a hybrid co-evolutionary-based approach to predict physical interplays between pairs of protein families, starting from protein sequences only. RESULTS: In the present analysis, pairs of multiple sequence alignments are constructed for each dimer and the covariation between residues in those pairs are calculated by CCMpred (Contacts from Correlated Mutations predicted) and three mutual information based approaches for ten accessible surface area threshold groups. Then, whole residue couplings between proteins of each dimer are unified into a single Frobenius norm value. Norms of residue contact matrices of all dimers in different accessible surface area thresholds are fed into support vector machine as single or multiple feature models. The results of training the classifiers by single features show no apparent different accuracies in distinct methods for different accessible surface area thresholds. Nevertheless, mutual information product and context likelihood of relatedness procedures may roughly have an overall higher and lower performances than other two methods for different accessible surface area cut-offs, respectively. The results also demonstrate that training support vector machine with multiple norm features for several accessible surface area thresholds leads to a considerable improvement of prediction performance. In this context, CCMpred roughly achieves an overall better performance than mutual information based approaches. The best accuracy, sensitivity, specificity, precision and negative predictive value for that method are 0.98, 1, 0.962, 0.96, and 0.962, respectively. CONCLUSIONS: In this paper, by feeding norm values of protein dimers into support vector machines in different accessible surface area thresholds, we demonstrate that even small number of proteins in pairs of multiple alignments could allow one to accurately discriminate between positive and negative dimers.

LogicNet: probabilistic continuous logics in reconstructing gene regulatory networks.

BACKGROUND: Gene Regulatory Networks (GRNs) have been previously studied by using Boolean/multi-state logics. While the gene expression values are usually scaled into the range [0, 1], these GRN inference methods apply a threshold to discretize the data, resulting in missing information. Most of studies apply fuzzy logics to infer the logical gene-gene interactions from continuous data. However, all these approaches require an a priori known network structure. RESULTS: Here, by introducing a new probabilistic logic for continuous data, we propose a novel logic-based approach (called the LogicNet) for the simultaneous reconstruction of the GRN structure and identification of the logics among the regulatory genes, from the continuous gene expression data. In contrast to the previous approaches, the LogicNet does not require an a priori known network structure to infer the logics. The proposed probabilistic logic is superior to the existing fuzzy logics and is more relevant to the biological contexts than the fuzzy logics. The performance of the LogicNet is superior to that of several Mutual Information-based and regression-based tools for reconstructing GRNs. CONCLUSIONS: The LogicNet reconstructs GRNs and logic functions without requiring prior knowledge of the network structure. Moreover, in another application, the LogicNet can be applied for logic function detection from the known regulatory genes-target interactions. We also conclude that computational modeling of the logical interactions among the regulatory genes significantly improves the GRN reconstruction accuracy.

Mid-infrared hybrid Si/VO2 modulator electrically driven by graphene electrodes.

Silicon photonic platforms are of significant interest for a variety of applications that operate in the mid-infrared regime. However, the realization of efficient mid-IR modulators, key components in any integrated optics platform, is still a challenging topic. Here, an ultra-compact high-speed hybrid Si/VO2 modulator operating at a mid-IR wavelength of 3.8 µm is presented. Electrical properties of graphene are employed to achieve a reversible insulating-metal phase transition in VO2 by electrical actuation. The thermal characteristics of graphene are employed to improve the response time of the VO2 phase transition through speed up heating and dissipation processes, thus enhancing the modulation speed. Optical and thermal simulations show an extinction ratio of 4.4 dB/µm, an insertion loss of 0.1 dB/µm, and high modulation speed of 23 ns. A larger modulation depth as high as 10 dB/µm can be achieved at the cost of lower modulation speed.

From indication to decision: A hierarchical approach to model the chemotactic behavior of Escherichia coli.

Reducing the complex behavior of living entities to its underlying physical and chemical processes is a formidable task in biology. Complex behaviors can be characterized as decision making: the ability to process the incoming information via an intracellular network and act upon this information to choose appropriate strategies. Motility is one such behavior that has been the focus many modeling efforts in the past. Our aim is to reduce the chemotactic behavior in Escherichia coli to its molecular constituents in order to paint a comprehensive and end-to-end picture of this intricate behavior. We utilize a hierarchical approach, consisting of three layers, to achieve this goal: at the first level, chemical reactions involved in chemotaxis are simulated. In the second level, the chemical reactions give rise to the mechanical movement of six independent flagella. At the last layer, the two lower layers are combined to allow a digital bacterium to receive information from its environment and swim through it with verve. Our results are in concert with the experimental studies concerning the motility of E.coli cells. In addition, we show that our detailed model of chemotaxis is reducible to a non-homogeneous Markov process.

A Stochastic Model of DNA Double-Strand Breaks Repair Throughout the Cell Cycle.

Cell cycle phase is a decisive factor in determining the repair pathway of DNA double-strand breaks (DSBs) by non-homologous end joining (NHEJ) or homologous recombination (HR). Recent experimental studies revealed that 53BP1 and BRCA1 are the key mediators of the DNA damage response (DDR) with antagonizing roles in choosing the appropriate DSB repair pathway in G1, S, and G2 phases. Here, we present a stochastic model of biochemical kinetics involved in detecting and repairing DNA DSBs induced by ionizing radiation during the cell cycle progression. A three-dimensional stochastic process is defined to monitor the cell cycle phase and DSBs repair at times after irradiation. To estimate the model parameters, a Metropolis Monte Carlo method is applied to perform maximum likelihood estimation utilizing the kinetics of γ-H2AX and RAD51 foci formation in G1, S, and G2 phases. The recruitment of DSB repair proteins is verified by comparing our model predictions with the corresponding experimental data on human cells after exposure to X and γ-radiation. Furthermore, the interaction between 53BP1 and BRCA1 is simulated for G1 and S/G2 phases determining the competition between NHEJ and HR pathways in repairing induced DSBs throughout the cell cycle. In accordance with recent biological data, the numerical results demonstrate that the maximum proportion of HR occurs in S phase cells and the high level of NHEJ takes place in G1 and G2 phases. Moreover, the stochastic realizations of the total yield of simple and complex DSBs ligation are compared for G1 and S/G2 damaged cells. Finally, the proposed stochastic model is validated when DSBs induced by different particle radiation such as iron, silicon, oxygen, proton, and carbon.

Enhanced quasi-phase-matched second-harmonic generation in a silicon-organic hybrid microring resonator.

An enhanced second-harmonic generation process in a dually resonant silicon-organic hybrid microring resonator is presented. The necessary phase matching is fulfilled through quasi-phase matching by periodically altering the spatial distribution of the second-order nonlinearity ${\chi ^{(2)}}$χ(2) of polymer infiltrated into a slot layer along the microring resonator propagation length. Absolute conversion efficiency up to ${P_{2\omega }}/{P_\omega } = {38}\% $P2ω/Pω=38% with only 10 mW pump power at $\lambda = 3.1\;{\unicode{x00B5}{\rm m}}$λ=3.1µm is obtained in an optimized structure, which shows a ${\sim}{76}$∼76-fold improvement compared to that of the straight waveguide. Also, our efficiency is comparable to that of previous works, even with the much lower pump power. We expect that our device will provide a potential application for bridging the mid-infrared and near-infrared bands on silicon photonics platforms.

Molecular Basis of Bicalutamide Response Alteration of Androgen Receptor Caused by Single Nucleotide Polymorphisms: An In Silico Investigation.

The vast majority of drugs act through binding to their protein targets. Prediction of the interaction between small molecules and these receptors is a key element in the process of drug discovery. Advances in structural biology have enabled us to resolve the three-dimensional structure of proteins, which are the targets of the drugs. Pharmacogenetics also helped researchers to study the structural variations arise from the single nucleotide polymorphisms (SNPs) and to survey the effects these variations in drug design and development. These improvements led to the identification of structural changes caused by SNPs, which affect the drug interaction with their receptors, called drug response. In this study, the interaction between androgen receptor and bicalutamide was investigated using a computational analysis. The results of these analyses were then used for identification of nonsynonymous SNPs that are potentially involved in drug response alterations. The data show that amino acids Met895, Trp741, Arg752, Ile899, Leu707, Gly708, Gln711, Met745, Met749, Thr877, Phe764, Met742, Asn705 and Leu704 are the main residues involved in the interaction between androgen receptor and bicalutamide. The occurrence of nonsynonymous polymorphisms I843T, L708R, H690P, I870M, N757S, L713F, G744E, L678P, M788V, M781I, A722T, H875Y, I842V, and F827L in this receptor greatly affected its interaction with bicalutamide, and they were able to cause drug resistance. The results of this study could be useful in predicting the response to treatment in patients receiving bicalutamide.

Identification of a missense variant in CLDN2 in obstructive azoospermia.

Obstructive azoospermia (OA), defined as an obstruction in any region of the male genital tract, accounts for 40% of all azoospermia cases. Of all OA cases, ~30% are thought to have a genetic origin, however, hitherto, the underlying genetic etiology of the majority of these cases remain unknown. To address this, we took a family-based whole-exome sequencing approach to identify causal variants of OA in a multiplex family with epidydimal obstruction. A novel gain-of-function missense variant in CLDN2 (c.481G>C; p.Gly161Arg) was found to co-segregate with the phenotype, consistent with the X-linked inheritance pattern observed in the pedigree. To assess the pathogenicity of this variant, the wild and mutant protein structures were modeled and their potential for strand formation in multimeric form was assessed and compared. The results showed that dimeric and tetrameric arrangements of Claudin-2 were not only reduced, but were also significantly altered by this single residue change. We, therefore, envisage that this amino acid change likely forms a polymeric discontinuous strand, which may lead to the disruption of tight junctions among epithelial cells. This missense variant is thus likely to be responsible for the disruption of the blood-epididymis barrier, causing dislodged epithelial cells to clog the genital tract, hence causing OA. This study not only sheds light on the underlying pathobiology of OA, but also provides a basis for more efficient diagnosis in the clinical setting.

New insights into resistance of Helicobacter pylori against third- and fourth-generation fluoroquinolones: A molecular docking study of prevalent GyrA mutations.

BACKGROUND: Fluoroquinolones hinder bacterial DNA replication by inhibiting DNA gyrase. However, mutations, in the QRDR segment of its A subunit (GyrA), cause antibiotic resistance. Here, the interactions of levofloxacin (LVX), gemifloxacin (GXN), and moxifloxacin (MXN) with Helicobacter pylori GyrA, in LVX-resistant vs -sensitive strains, were studied. METHODS: Levoflixacin-sensitive (n = 4) and -resistant (n = 9) H pylori strains, randomly selected from another antibiotic susceptibility study, underwent PCR amplification of gyrA gene, spanning the QRDR segment. The amplified gene fragments were sequenced and aligned. The homology model of H pylori GyrA was built based on that of Escherichia coli, and energy minimization was done. The interaction patterns of LVX, GXN, and MXN with GyrA were analyzed via molecular docking studies. RESULTS: Sequence alignment of the 13 studied strains, created 5 categories of strains: (A) wild type-like (H pylori ATCC26695), (B) N87K-only, (C) D91N-only, (D) N87K + V94L, and (E) D91N + A97V mutations. The minimum inhibitory concentrations (MIC) for LVX-sensitive (category A) and -resistant (categories B-E) strains were <1 mg/L and ≥32 mg/L, respectively. The binding mode of GyrA in category A with LVX identified G35/N87/Y90/D91/V94/G114/S115/I116/D117/G118/D119, as key residues, some residing outside the QRDR segment. Category B strains lost only one interaction (G35), which led to elevated binding free energy (∆G) and full LVX resistance. Categories C-E lost more contacts, with higher ∆G and again full LVX resistance. GXN bound to GyrA of categories A and B via a different set of key residues, while MXN retained the lost contact (G35) in LVX-resistant, category B strains. CONCLUSION: Using molecular docking tools, we identified the key residues responsible for interaction of GyrA with LVX, GXN, and MXN. In the presence of N87K-only mutation, the loss of one of these contacts (ie, G35) led to full LVX resistance. Yet, GXN and MXN overcame this mutation, by retaining all key contacts with GyrA.

Assessing the impact of exact reads on reducing the error rate of read mapping.

BACKGROUND: Nowadays, according to valuable resources of high-quality genome sequences, reference-based assembly methods with high accuracy and efficiency are strongly required. Many different algorithms have been designed for mapping reads onto a genome sequence which try to enhance the accuracy of reconstructed genomes. In this problem, one of the challenges occurs when some reads are aligned to multiple locations due to repetitive regions in the genomes. RESULTS: In this paper, our goal is to decrease the error rate of rebuilt genomes by resolving multi-mapping reads. To achieve this purpose, we reduce the search space for the reads which can be aligned against the genome with mismatches, insertions or deletions to decrease the probability of incorrect read mapping. We propose a pipeline divided to three steps: ExactMapping, InExactMapping, and MergingContigs, where exact and inexact reads are aligned in two separate phases. We test our pipeline on some simulated and real data sets by applying some read mappers. The results show that the two-step mapping of reads onto the contigs generated by a mapper such as Bowtie2, BWA and Yara is effective in improving the contigs in terms of error rate. CONCLUSIONS: Assessment results of our pipeline suggest that reducing the error rate of read mapping, not only can improve the genomes reconstructed by reference-based assembly in a reasonable running time, but can also have an impact on improving the genomes generated by de novo assembly. In fact, our pipeline produces genomes comparable to those of a multi-mapping reads resolution tool, namely MMR by decreasing the number of multi-mapping reads. Consequently, we introduce EIM as a post-processing step to genomes reconstructed by mappers.