Sequence Conservation Analysis Research Articles

Genome-wide detection and sequence conservation analysis of long non-coding RNA during hair follicle cycle of yak

BackgroundLong non-coding RNA (lncRNA) as an important regulator has been demonstrated playing an indispensable role in the biological process of hair follicles (HFs) growth. However, their function and expression profile in the HFs cycle of yak are yet unknown. Only a few functional lncRNAs have been identified, partly due to the low sequence conservation and lack of identified conserved properties in lncRNAs. Here, lncRNA-seq was employed to detect the expression profile of lncRNAs during the HFs cycle of yak, and the sequence conservation of two datasets between yak and cashmere goat during the HFs cycle was analyzed.ResultsA total of 2884 lncRNAs were identified in 5 phases (Jan., Mar., Jun., Aug., and Oct.) during the HFs cycle of yak. Then, differential expression analysis between 3 phases (Jan., Mar., and Oct.) was performed, revealing that 198 differentially expressed lncRNAs (DELs) were obtained in the Oct.-vs-Jan. group, 280 DELs were obtained in the Jan.-vs-Mar. group, and 340 DELs were obtained in the Mar.-vs-Oct. group. Subsequently, the nearest genes of lncRNAs were searched as the potential target genes and used to explore the function of DELs by GO and KEGG enrichment analysis. Several critical pathways involved in HFs development such as Wnt signaling pathway, VEGF signaling pathway, and signaling pathways regulating pluripotency of stem cells, were enriched. To further screen key lncRNAs influencing the HFs cycle, 24 DELs with differ degree of sequence conservation were obtained via a comparative analysis of partial DELs with previously published lncRNA-seq data of cashmere goat in the HFs cycle using NCBI BLAST-2.9.0+, and 3 DELs of them were randomly selected for further detailed analysis of the sequence conservation properties.ConclusionsThis study revealed the expression pattern and potential function of lncRNAs during HFs cycle of yak, which would expand the knowledge about the role of lncRNAs in the HFs cycle. The findings related to sequence conservation properties of lncRNAs in the HFs cycle between the two species may provide valuable insights into the study of lncRNA functionality and mechanism.

In silico model of mtDNA mutations effect on secondary and 3D structure of mitochondrial rRNA and tRNA in Leber’s hereditary optic neuropathy

The Leber’s hereditary optic neuropathy (LHON) is a rare disease caused by mitochondrial DNA (mtDNA) mutations. Beside primary mutations, the effect of secondary mtDNA mutations in still unclear. We examined the effect of secondary mtDNA mutations on secondary structure of different mitochondrial RNAs. Whole mitochondrial genome sequence of LHON patients has been obtained from in six non related pedigrees by Sanger sequencing method. The effect of mutations located in mitochondrial RNA genes was examined by creating in silico models of RNA secondary and regional 3D structure, accompanied by sequence conservation analysis. All three primary LHON mutations (m.3460G>A, m.11778G>A and m.14484 T>C) were revealed in study families. Four mutations in MT-RNR1 gene (m.750A>G, m.956delC, m.1438A>G and m.1555A>G) were identified and only an m.1555A>G causes significant changes of secondary structure of mitochondrial 12S ribosomal RNA (rRNA), while it is the only mutation which does not alter its 3D structure. Five mutations (m.1811A>G, m.2706A>G, m.2831G>A, m.3010G>A and m.3197T>C) were discovered in MT-RNR2 gene and all of them induced substantial alterations of mitochondrial 16S rRNA secondary structure. Significant changes of mitochondrial 16S rRNA 3D structure are caused by m.1811A>G, m.2706A>G, m.3010G>A and m.3197T>C. A single insertion variant (m.15986insG) has been found in the MT-TP gene which encodes mitochondrial transfer RNA for Proline (tRNA Pro). This mutation does not cause substantial changes of tRNA for Proline secondary structure, while the 3D geometry remains without major changes. Most of the mutation loci exhibited high level of sequence conservation. Presence of multiple mutations in a single family appears to cause more extensive changes in mitochondrial 12S and 16S rRNA, then their individual influence. The effect of discovered mutations on in silico modelled RNA structure is in a significant correlation with the present knowledge about the potential of these mutation to participate in the pathophysiology of LHON and other human diseases. The presence of certain multiple mitochondrial RNA mutations could be a possible explanation of LHON clinical presentation in some families. All revealed mutations have been evaluated for the first time in terms of in silico structural modelling. The application of bioinformatics tools such as secondary and 3D RNA structure prediction can have a great advantage in better understanding of the molecular standpoint of the LHON pathophysiology and clinical phenotype.

Improving the thermal stability of Bacillus subtilis lipase based on multiple computational design strategies

Improving the thermal stability of enzymes is a hot and difficult point in the field of biocatalysis. Compared with the traditional directed evolution, computational assisted rational design is more efficient, and is widely used in enzyme engineering. Using Bacillus subtilis LipA as the model protein, the structure cavity of the enzyme was analyzed by Rosetta-VIP design, the mutation which was beneficial to the filling of the structure cavity (ΔΔE<0) was selected, followed by the solvent accessible surface area and evolutionary conservation analysis. The thermal stabilities of six out of sixteen designed single-point mutants were improved, with a maximum ΔTm value of 3.18 °C. These six mutations were further used for iterative combination mutation, the maximum ΔTm of the two-point and three-point combination mutants were 4.04 °C and 5.13 °C, respectively. The Tm of the four-point combination mutant M11 (F17A/L114P/I135V/M137L) was increased by 7.30 °C. The Tm of the six-point combination mutant M10 (F17A/V74I/L114P/I135V/M137A/I157L) was increased by 7.43 °C. The thermal stability of mutation with lower energy value, reduced accessible surface area, while conformed to evolutionary conservatism, was more likely to be improved. Therefore, the multiple virtual screening strategy based on the enzyme structure cavity filling, solvent accessible surface area and amino acid sequence conservation analysis can effectively improve the thermal stability of enzyme.

The Bio3D packages for structural bioinformatics.

Bio3D is a family of R packages for the analysis of biomolecular sequence, structure, and dynamics. Major functionality includes biomolecular database searching and retrieval, sequence and structure conservation analysis, ensemble normal mode analysis, protein structure and correlation network analysis, principal component, and related multivariate analysis methods. Here, we review recent package developments, including a new underlying segregation into separate packages for distinct analysis, and introduce a new method for structure analysis named ensemble difference distance matrix analysis (eDDM). The eDDM approach calculates and compares atomic distance matrices across large sets of homologous atomic structures to help identify the residue wise determinants underlying specific functional processes. An eDDM workflow is detailed along with an example application to a large protein family. As a new member of the Bio3D family, the Bio3D-eddm package supports both experimental and theoretical simulation-generated structures, is integrated with other methods for dissecting sequence-structure-function relationships, and can be used in a highly automated and reproducible manner. Bio3D is distributed as an integrated set of platform independent open source R packages available from: http://thegrantlab.org/bio3d/.

Comparative Dynamics and Functional Mechanisms of the CYP17A1 Tunnels Regulated by Ligand Binding.

As an important member of cytochrome P450 (CYP) enzymes, CYP17A1 is a dual-function monooxygenase with a critical role in the synthesis of many human steroid hormones, making it an attractive therapeutic target. The emerging structural information about CYP17A1 and the growing number of inhibitors for these enzymes call for a systematic strategy to delineate and classify mechanisms of ligand transport through tunnels that control catalytic activity. In this work, we applied an integrated computational strategy to different CYP17A1 systems with a panel of ligands to systematically study at the atomic level the mechanism of ligand-binding and tunneling dynamics. Atomistic simulations and binding free energy computations identify the dynamics of dominant tunnels and characterize energetic properties of critical residues responsible for ligand binding. The common transporting pathways including S, 3, and 2c tunnels were identified in CYP17A1 binding systems, while the 2c tunnel is a newly formed pathway upon ligand binding. We employed and integrated several computational approaches including the analysis of functional motions and sequence conservation, atomistic modeling of dynamic residue interaction networks, and perturbation response scanning analysis to dissect ligand tunneling mechanisms. The results revealed the hinge-binding and sliding motions as main functional modes of the tunnel dynamic, and a group of mediating residues as key regulators of tunnel conformational dynamics and allosteric communications. We have also examined and quantified the mutational effects on the tunnel composition, conformational dynamics, and long-range allosteric behavior. The results of this investigation are fully consistent with the experimental data, providing novel rationale to the experiments and offering valuable insights into the relationships between the structure and function of the channel networks and a robust atomistic model of activation mechanisms and allosteric interactions in CYP enzymes.

Molecular, Evolutionary, and Structural Analysis of the Terminal Protein Domain of Hepatitis B Virus Polymerase, a Potential Drug Target.

Approximately 250 million people are living with chronic hepatitis B virus (HBV) infections, which claim nearly a million lives annually. The target of all current HBV drug therapies (except interferon) is the viral polymerase; specifically, the reverse transcriptase domain. Although no high-resolution structure exists for the HBV polymerase, several recent advances have helped to map its functions to specific domains. The terminal protein (TP) domain, unique to hepadnaviruses such as HBV, has been implicated in the binding and packaging of the viral RNA, as well as the initial priming of and downstream synthesis of viral DNA—all of which make the TP domain an attractive novel drug target. This review encompasses three types of analysis: sequence conservation analysis, secondary structure prediction, and the results from mutational studies. It is concluded that the TP domain of HBV polymerase is comprised of seven subdomains (three unstructured loops and four helical regions) and that all three loop subdomains and Helix 5 are the major determinants of HBV function within the TP domain. Further studies, such as modeling inhibitors of these critical TP subdomains, will advance the TP domain of HBV polymerase as a therapeutic drug target in the progression towards a cure.

Molecular insight into human P-glycoprotein allosteric transition from outward- to inward-facing state

Human P-glycoprotein (P-gp) belongs to ATP-Binding Cassette (ABC) transporter, which can export varieties of anti-cancer drugs out of tumor cells, resulting in multidrug resistance (MDR) of tumor cells. We utilized targeted molecular dynamics simulation to explore the allosteric transition of human P-gp from outward- to inward-facing state. The results show that NBDs experience translational and rotational movement, caused by the successive hydrolysis of ATP molecules. The reorientation of the TMDs begins with the closing of the extracellular gate, followed by the opening of the cytoplasmic gate. The opening delay of the cytoplasmic gate is due to the extensive electrostatic interactions here. The sequence conservation analysis reveals that charged residues in TM3, 4 are highly conserved for achieving the unidirectional transport property. Our results give detailed movement and energy analyses to the reverse allosteric process of human P-gp and are helpful to the understanding of the working mechanism of human P-gp.

Whole exome sequencing of a Saudi family and systems biology analysis identifies CPED1 as a putative causative gene to Celiac Disease.

Celiac disease (CD) is a gastrointestinal disorder whose genetic basis is not fully understood. Therefore, we studied a Saudi family with two CD affected siblings to discover the causal genetic defect. Through whole exome sequencing (WES), we identified that both siblings have inherited an extremely rare and deleterious CPED1 genetic variant (c.241 A > G; p.Thr81Ala) segregating as autosomal recessive mutation, suggesting its putative causal role in the CD. Saudi population specific minor allele frequency (MAF) analysis has confirmed its extremely rare prevalence in homozygous condition (MAF is 0.0004). The Sanger sequencing analysis confirmed the absence of this homozygous variant in 100 sporadic Saudi CD cases. Genotype-Tissue Expression (GTEx) data has revealed that CPED1 is abundantly expressed in gastrointestinal mucosa. By using a combination of systems biology approaches like protein 3D modeling, stability analysis and nucleotide sequence conservation analysis, we have further established that this variant is deleterious to the structural and functional aspects of CPED1 protein. To the best of our knowledge, this variant has not been previously reported in CD or any other gastrointestinal disease. The cell culture and animal model studies could provide further insight into the exact role of CPED1 p.Thr81Ala variant in the pathophysiology of CD. In conclusion, by using WES and systems biology analysis, present study for the first-time reports CPED1 as a potential causative gene for CD in a Saudi family with potential implications to both disease diagnosis and genetic counseling.

Stable intronic sequence RNAs (sisRNAs) are selected regions in introns with distinct properties

BackgroundStable introns and intronic fragments make up the largest population of RNA in the oocyte nucleus of the frog Xenopus tropicalis. These stable intronic sequence RNAs (sisRNAs) persist through the onset of zygotic transcription when synchronous cell division has ended, and the developing embryo consists of approximately 8000 cells. Despite their abundance, the sequence properties and biological function of sisRNAs are just beginning to be understood.ResultsTo characterize this population of non-coding RNA, we identified all of the sisRNAs in the X. tropicalis oocyte nucleus using published high-throughput RNA sequencing data. Our analysis revealed that sisRNAs, have an average length of ~ 360 nt, are widely expressed from genes with multiple introns, and are derived from specific regions of introns that are GC and TG rich, while CpG poor. They are enriched in introns at both ends of transcripts but preferentially at the 3′ end. The consensus binding sites of specific transcription factors such as Stat3 are enriched in sisRNAs, suggesting an association between sisRNAs and transcription factors involved in early development. Evolutionary conservation analysis of sisRNA sequences in seven vertebrate genomes indicates that sisRNAs are as conserved as other parts of introns, but much less conserved than exons.ConclusionIn total, our results indicate sisRNAs are selected intron regions with distinct properties and may play a role in gene expression regulation.

Rhapsody: predicting the pathogenicity of human missense variants.

MotivationThe biological effects of human missense variants have been studied experimentally for decades but predicting their effects in clinical molecular diagnostics remains challenging. Available computational tools are usually based on the analysis of sequence conservation and structural properties of the mutant protein. We recently introduced a new machine learning method that demonstrated for the first time the significance of protein dynamics in determining the pathogenicity of missense variants.ResultsHere, we present a new interface (Rhapsody) that enables fully automated assessment of pathogenicity, incorporating both sequence coevolution data and structure- and dynamics-based features. Benchmarked against a dataset of about 20 000 annotated variants, the methodology is shown to outperform well-established and/or advanced prediction tools. We illustrate the utility of Rhapsody by in silico saturation mutagenesis studies of human H-Ras, phosphatase and tensin homolog and thiopurine S-methyltransferase.Availability and implementationThe new tool is available both as an online webserver at http://rhapsody.csb.pitt.edu and as an open-source Python package (GitHub repository: https://github.com/prody/rhapsody; PyPI package installation: pip install prody-rhapsody). Links to additional resources, tutorials and package documentation are provided in the 'Python package' section of the website.Supplementary information Supplementary data are available at Bioinformatics online.

Identification and characterization of protective CD8+ T-epitopes in a malaria vaccine candidate SLTRiP.

IntroductionEfforts are required at developing an effective vaccine that can inhibit malaria prevalence and transmission. Identifying the critical immunogenic antigens and understanding their interactions with host proteins forms a major focus of subunit vaccine development. Previously, our laboratory showed that SLTRiP conferred protection to the liver stage of Plasmodium growth in rodents. In the follow‐up of earlier research, we demonstrate that SLTRiP‐mediated protection is majorly concentrated in specific regions of protein.MethodTo identify particular protective regions of protein, we synthesized multiple nonoverlapping fragments from SLTRiP protein. From this, we designed a panel of 8‐20mer synthetic peptides, which were predicted using T‐epitope‐based prediction algorithm. We utilized the IFN‐γ enzyme‐linked immunosorbent spot assay to identify immunodominant peptides. The latter were used to immunize mice, and these mice were challenged to assess protection.ResultsThe protective polypeptide fragment SLTRiP C3 and SLTRiP C4 were identified, by expressing and testing multiple fragments of PbSLTRiP protein. The immune responses generated by these fragments were compared to identify the immunodominant fragment. The T‐epitopes were predicted from SLTRiP protein using computer‐based algorithms. The in vitro immune responses generated by these peptides were compared with each other to identify the immunodominant T‐epitope. Immunization using these peptides showed significant reduction in parasite numbers during liver stage.ConclusionOur findings show that the protective efficacy shown by SLTRiP is localized in particular protein fragments. The peptides designed from such regions showed protective efficacy equivalent to whole protein. The sequence conservation analysis with human Plasmodium species also showed that these peptides were conserved. In conclusion, these peptides or their equivalent from other Plasmodium species could impart protection against malaria in their respective hosts too. Our studies provide a basis for the inclusion of these peptides in clinical vaccine constructs against malaria.

Conservation analysis of SARS-CoV-2 spike suggests complicated viral adaptation history from bat to human.

BackgroundThe current coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome (SARS)-CoV-2, has become the most devastating public health emergency in the 21st century and one of the most influential plagues in history. Studies on the origin of SARS-CoV-2 have generally agreed that the virus probably comes from bat, closely related to a bat CoV named BCoV-RaTG13 taken from horseshoe bat (Rhinolophus affinis), with Malayan pangolin (Manis javanica) being a plausible intermediate host. However, due to the relatively low number of SARS-CoV-2-related strains available in public domain, the evolutionary history remains unclear.MethodologyNine hundred ninety-five coronavirus sequences from NCBI Genbank and GISAID were obtained and multiple sequence alignment was carried out to categorize SARS-CoV-2 related groups. Spike sequences were analyzed using similarity analysis and conservation analyses. Mutation analysis was used to identify variations within receptor-binding domain (RBD) in spike for SARS-CoV-2-related strains.ResultsWe identified a family of SARS-CoV-2-related strains, including the closest relatives, bat CoV RaTG13 and pangolin CoV strains. Sequence similarity analysis and conservation analysis on spike sequence identified that N-terminal domain, RBD and S2 subunit display different degrees of conservation with several coronavirus strains. Mutation analysis on contact sites in SARS-CoV-2 RBD reveals that human-susceptibility probably emerges in pangolin.Conclusion and implicationWe conclude that the spike sequence of SARS-CoV-2 is the result of multiple recombination events during its transmission from bat to human, and we propose a framework of evolutionary history that resolve the relationship of BCoV-RaTG13 and pangolin coronaviruses with SARS-CoV-2.Lay SummaryThis study analyses whole-genome and spike sequences of coronavirus from NCBI using phylogenetic and conservation analyses to reconstruct the evolutionary history of severe acute respiratory syndrome (SARS)-CoV-2 and proposes an evolutionary history of spike in the progenitors of SARS-CoV-2 from bat to human through mammal hosts before they recombine into the current form.

Techniques for Small Non-Coding RNA Analysis in Seeds of Forest Tree Species.

In recent years, the scientific community has become aware that epigenetic mechanisms play a more important role in gene regulatory networks (GRNs) than was hitherto thought, as accumulating evidence has shown that changes in epigenetics without genetic variation can affect complex traits over multiple generations. Within the epigenetic machinery, small non-coding RNAs (sRNAs, 18-24 nucleotides in length) are evolutionarily conserved RNA molecules that target mRNAs for deregulation or translational repression. They commonly have high-level regulatory functions in GRNs by mediating DNA and/or histone methylation and gene silencing essential for plant developmental programs and adaptability. Local adaptation enables plants to acquire a high fitness by, for example, properly timing developmental transitions to match plant growth stages with organism's favorable seasons. In particular, the seed represents a key evolutionary adaptation of seed plants that facilitates dispersal and reinitiates the development coupled in time with suitable environmental conditions. With the advent of high-throughput sequencing for sRNAs and computational approaches for sRNA detection and categorization, it is now feasible to unravel how sRNAs contribute to the fitness of tree species that can survive hundreds of years (e.g., conifers). Of particular interest is to disentangle the roles of sRNAs from complex genomic information in tree species with intimidating genomic sizes (commonly 20-30Gb in conifers) and abundant nongenic components (e.g., >60% transposable elements). In this chapter, we use seeds of the conifer Picea glauca as a study system to describe the methods and protocols we used or have recently updated, from high-quality RNA isolation to sRNA identification, sequence conservation, abundance comparison, and functional analysis.

Characterization of recombinant glutathione reductase from Antarctic yeast Rhodotorula mucilaginosa

Glutathione reductase (GR) catalyzes the reduction of glutathione disulfide, which helps to maintain a cellular reducing environment during stress in organisms. However, GR of polar yeast has not been well-characterized so far. To fully understand the molecular and enzymatic properties of GR from extreme area and broaden its knowledge, a cDNA-encoding GR from Antarctic sea-ice yeast Rhodotorula mucilaginosa (designated as RmGR) was cloned and expressed in Escherichia coli. The open reading frame of RmGR was 1500 bp, encoding an enzyme of 499 amino acids with a predicted pI of 6.07 and molecular weight of 54.8 kDa. SDS-PAGE and gel filtration analysis results showed that the RmGR was a homodimer. Conserved sequence analysis revealed that RmGR behaved typical characteristics of GR, containing a pyridine nucleotide-disulfide oxidoreductase active site, a flavin adenine dinucleotide (FAD), and reduced nicotinamide adenine dinucleotide phosphate (NADPH) binding motifs and two glutathione oxidized binding motifs. The recombinant enzyme displayed relatively high stability at pH 3.0–6.0 and 20–40 °C with optimal enzymatic activity at 30 °C and pH 7.5. Real-time quantitative PCR data showed that the expression of RmGR gene was upregulated after copper treatment of yeast cells. Moreover, RmGR-transformed E. coli exhibited stronger growth profiles than cells with an empty vector after copper and cadmium treatment. Our results demonstrated the possible function of RmGR in adaptation to heavy metals and its potential application in heterologous expression systems.

Identification of an Intronic Regulatory Element Necessary for Tissue-Specific Expression of Foxn1 in Thymic Epithelial Cells.

The thymus is critical for the establishment of the adaptive immune system and the development of a diverse T cell repertoire. T cell development depends upon cell-cell interactions with epithelial cells in the thymus. The thymus is composed of two different types of epithelial cells: cortical and medullary epithelial cells. Both of these express and critically depend on the transcription factor Foxn1 Foxn1 is also expressed in the hair follicle, and disruption of Foxn1 function in mice results in severe thymic developmental defects and the hairless (nude) phenotype. Despite its importance, little is known about the direct regulation of Foxn1 expression. In this study, we identify a cis-regulatory element (RE) critical for expression of Foxn1 in mouse thymic epithelial cells but dispensable for expression in hair follicles. Analysis of chromatin accessibility, histone modifications, and sequence conservation identified regions within the first intron of Foxn1 that possessed the characteristics of REs. Systematic knockout of candidate regions lead us to identify a 1.6 kb region that, when deleted, results in a near total disruption of thymus development. Interestingly, Foxn1 expression and function in the hair follicle were unaffected. RNA fluorescent in situ hybridization showed a near complete loss of Foxn1 mRNA expression in the embryonic thymic bud. Our studies have identified a genomic RE with thymic-specific control of Foxn1 gene expression.

Use of Transcriptomic Analysis for the Design of Antibodies to Serve as Neuronal and Synaptic Markers in Macrobrachium crustaceans

Crustaceans are major constituents of aquatic ecosystems and, as such, changes in their behavior and the structure and function of their bodies can serve as indicators of alterations in their immediate environment, such as those associated with climate change and anthropogenic contamination. We want to develop markers of crustacean neural circuits and terminal synaptic structure that can be used for quantifiable characterization with confocal microscopy. When compared with other invertebrate species such as Drosophila and C. elegans, little is known regarding crustacean genomic information, particularly of the animal's nervous system. We have used bioinformatics and a de novo transcriptome assembly approach to identify potential targets for developing specific antibodies to serve as nervous system markers for freshwater prawns of the Macrobrachium spp. These species of prawn are found in rivers of Puerto Rico and are being used to study the impact of exposure to emergent contaminants, such as phthalates and metabolites of pharmaceuticals, in natural sources of freshwater near urban areas.Total RNA was extracted from brain ganglia of Macrobrachium carcinus freshwater prawns and Illumina Next Generation Sequencing was performed using an Eel Pond mRNA Seq Protocol to construct a de novo transcriptome. Transcript sequences related to glutamate receptors were ran through multiple sequence conservation analysis using the MUSCLE webserver. Conserved sequences were analyzed for antigenicity to design specific antibodies. Raw data from the Illumina Sequencing yielded 97,202,662 sequences: 47,630,546 paired and 1,941,570 singletons. Assembly with Trinity resulted in 197,898 assembled sequences, of which 30,576 were annotated: 9,600 annotated by orthology, 17,197 annotated by homology, and 3,779 annotated by transcript families. In the search for glutamate receptor sequences, we found 138 sequences related to ionotropic receptors, 32 sequences related to metabotropic receptors, and 18 sequences to unidentified receptors. Using this information, we were able to develop antibodies for prawn AMPA ionotropic glutamate receptor 1 (GluRIA1), metabotropic glutamate receptor subunit 1 and 4 (mGluRG1S1 and mGluRG3S4), and NMDA ionotropic receptor 2B (GluRNMDA2B). Immunohistochemistry of the ventral nerve cord has thus far shown both punctuate and cell body staining of mGluRG1S1 located throughout all the central nervous system, including the eyestalks, brain, and the circumesophageal, subesophageal, thoracic and abdominal ganglia. The other antibodies are presently being tested in the prawn central nervous system as well as in peripheral synapses such as the neuromuscular junction.Support or Funding InformationNSF DBI‐1337284NSF CREST NSF HRD‐1736019NIH NIMHD 8G12‐MD007600 (RCMI)NIGMS‐RISE R25 GM061838 NSF (USA): CREST HRD‐1137725This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Structural characterization and molecular dynamics simulations of the caprine and bovine solute carrier family 11 A1 (SLC11A1).

Natural Resistance-Associated Macrophage Proteins are a family of transmembrane divalent metal ion transporters, with important implications in life of both bacteria and mammals. Among them, the Solute Carrier family 11 member A1 (SLC11A1) has been implicated with susceptibility to infection by Mycobacterium avium subspecies paratuberculosis (MAP), potentially causing Crohn's disease in humans and paratuberculosis (PTB) in ruminants. Our previous research had focused on sequencing the mRNA of the caprine slc11a1 gene and pinpointed polymorphisms that contribute to caprine SLC11A1's susceptibility to infection by MAP in PTB. Despite its importance, little is known on the structural/dynamic features of mammalian SLC11A1 that may influence its function under normal or pathological conditions at the protein level. In this work we studied the structural architecture of SLC11A1 in Capra hircus and Bos taurus through molecular modeling, molecular dynamics simulations in different, functionally relevant configurations, free energy calculations of protein-metal interactions and sequence conservation analysis. The results of this study propose a three dimensional structure for SLC11A1 with conserved sequence and structural features and provide hints for a potential mechanism through which divalent metal ion transport is conducted. Given the importance of SLC11A1 in susceptibility to PTB, this study provides a framework for further studies on the structure and dynamics of SLC11A1 in other organisms, to gain 3D structural insight into the macromolecular arrangements of SLC11A1 but also suggesting a potential mechanism which divalent metal ion transport is conducted.

The highly conserved domain of RND multidrug efflux pumps in pathogenic Gram-negative bacteria

RND (Resistance-Nodulation-Division) family transporters have a vital role in both intrinsic and acquired multi-drug resistance in Gram-negative bacteria. It is important to find a conserved domain in the RND family between different pathogenic bacteria for diagnostic and therapeutic purpose. Total sequences of three-component system RND efflux pumps were retrieved from NCBI nucleotide and protein database and were subjected to conservation and variation analysis using the multiple sequence alignment feature of the CLC workbench. The phylogenetic tree for main transporters was drawn and the three-dimensional structure was also evaluated. From the sequence conservation analysis, highly conserved residues with 282 base pair (94 amino acid) long were identified. The location of the highly conserved domain is positioned in the domain 1 crystallographic structure of AcrB Escherichia coli and MexB Pseudomonas aeruginosa. The main transporter component phylogenetic tree shows the clusters of different genotypes and their evolutionary association. Each of three components of RND proteins is crucial for drug efflux, and the absence of even one component makes the entire complex totally nonfunctional. Therefore, this highly conserved region can be used to disable the RND multidrug efflux pumps. In addition, this highly conserved can also be used for diagnostic aspects.

Engineering methionine γ-lyase from Citrobacter freundii for anticancer activity

Methionine deprivation of cancer cells, which are deficient in methionine biosynthesis, has been envisioned as a therapeutic strategy to reduce cancer cell viability. Methionine γ-lyase (MGL), an enzyme that degrades methionine, has been exploited to selectively remove the amino acid from cancer cell environment. In order to increase MGL catalytic activity, we performed sequence and structure conservation analysis of MGLs from various microorganisms. Whereas most of the residues in the active site and at the dimer interface were found to be conserved, residues located in the C-terminal flexible loop, forming a wall of the active site entry channel, were found to be variable. Therefore, we carried out site-saturation mutagenesis at four independent positions of the C-terminal flexible loop, P357, V358, P360 and A366 of MGL from Citrobacter freundii, generating libraries that were screened for activity. Among the active variants, V358Y exhibits a 1.9-fold increase in the catalytic rate and a 3-fold increase in KM, resulting in a catalytic efficiency similar to wild type MGL. V358Y cytotoxic activity was assessed towards a panel of cancer and nonmalignant cell lines and found to exhibit IC50 lower than the wild type. The comparison of the 3D-structure of V358Y MGL with other MGL available structures indicates that the C-terminal loop is either in an open or closed conformation that does not depend on the amino acid at position 358. Nevertheless, mutations at this position allosterically affects catalysis.

Characterization of mutation spectrum and identification of novel mutations in ATP7B gene from a cohort of Wilson disease patients: Functional and therapeutic implications

Wilson disease (WD), a copper metabolism disorder, occurs due to the presence of mutations in the gene encoding ATP7B, a protein that primarily facilitates hepatic copper excretion. A better understanding of spectrum and functional significance of ATP7B variants is critical to formulating targeted and personalized therapies. Henceforth, we screened and sequenced 21 exons of ATP7B gene from 50 WD patients and 60 healthy subjects. We identified 28 variants comprising, seven novels in 20% alleles, while eight variations affecting 23% alleles were first time reported in Indian cohort. The c.813C>A, p.(Cys271*) (10%) was the most frequent mutation. Bioinformatics analysis revealed five of seven novel variants viz. c.1600C>A, p.(Pro534Thr); c.1616C>A, p.(Pro539His); c.1924G>T, p.(Asp642Tyr); c.2168G>C, p.(Arg723Thr); c.2174G>C, p.(Arg725Thr) resulted in protein misfolding. Sequence conservation analysis of ATP7B regions containing novel variants documented an evolutionarily conserved nature. Functional analysis of these novel variants in five different cell lines lacking inherent ATP7B expression demonstrated sensitivity to CuCl2 -treatment, experiencing augmented cellular copper retention and decreased copper excretion as well as ceruloplasmin secretion to that of wildtype-ATP7B expressing cells. Interestingly, pharmacological chaperone 4-phenylbutyrate, a clinically approved compound, partially restored protein function of ATP7B mutants. These findings might enable novel treatment strategies in WD by clinically enhancing the protein expression of mutant ATP7B with residual copper export activity.