Exploring mutable conserved sites and fatal non-conserved sites by random mutation of esterase from Sulfolobus tokodaii and subtilisin from Thermococcus kodakarensis

Homologous proteins differ in their amino acid sequences at several positions. Generally, conserved sites are recognized as not suitable for amino acid substitution, and thus in evolutionary protein engineering, non-conserved sites are often selected as mutation sites. However, there have also been reports of possible mutations in conserved sites. In this study, we explored mutable conserved sites and immutable non-conserved sites by testing random mutations of two thermostable proteins, an esterase from Sulfolobus tokodaii (Sto-Est) and a subtilisin from Thermococcus kodakarensis (Tko-Sub). The subtilisin domain of Tko-Sub needs Ca²⁺ ions and the propeptide domain for stability, folding and maturation. The results from the two proteins showed that about one-third of the mutable sites were detected in conserved sites and some non-conserved sites lost enzymatic activity at high temperatures due to mutation. Of the conserved sites in Sto-Est, the sites on the loop, on the surface, and far from the active site are more resistant to mutation. In Tko-Sub, the sites flanking Ca²⁺-binding sites and propeptide were undesirable for mutation. The results presented here serve as an index for selecting mutation sites and contribute to the expansion of available sequence range by introducing mutations at conserved sites.

The development of a novel diagnostic PCR for Madurella mycetomatis using a comparative genome approach

Background

Eumycetoma is a neglected tropical disease most commonly caused by the fungus Madurella mycetomatis. Identification of eumycetoma causative agents can only be reliably performed by molecular identification, most commonly by species-specific PCR. The current M. mycetomatis specific PCR primers were recently discovered to cross-react with Madurella pseudomycetomatis. Here, we used a comparative genome approach to develop a new M. mycetomatis specific PCR for species identification.

Methodology

Predicted-protein coding sequences unique to M. mycetomatis were first identified in BLASTCLUST based on E-value, size and presence of orthologues. Primers were then developed for 16 unique sequences and evaluated against 60 M. mycetomatis isolates and other eumycetoma causing agents including the Madurella sibling species. Out of the 16, only one was found to be specific to M. mycetomatis.

Conclusion

We have discovered a predicted-protein coding sequence unique to M. mycetomatis and have developed a new species-specific PCR to be used as a novel diagnostic marker for M. mycetomatis.

Mycetoma is a neglected tropical disease characterised by tumorous swellings and grain formation. This disease can be caused by more than 70 different micro-organisms and is categorised into actinomycetoma (caused by bacteria) and eumycetoma (caused by fungi). The most common causative agent of mycetoma is the fungus Madurella mycetomatis. Diagnosis of eumycetoma is often only done clinically or by histopathological examination and culturing of the grains. Unfortunately, that often leads to misidentifications. Molecular identification is currently the most reliable method to identify the causative agents. However, we have recently discovered that the only M. mycetomatis species-specific PCR primers cross-reacts to Madurella pseudomycetomatis. Since all Madurella species cause eumycetoma and have different susceptibilities to antifungal agents, it is important to be able to accurately identify them to the species level. Here we have used a comparative genome approach to identify and design new M. mycetomatis species-specific PCR primers. These primers can be used to identify M. mycetomatis directly from grains and do not cross-react with any of the other eumycetoma causative agents tested. We, therefore, recommended the use of these primers in reference centres and local laboratories to identify M. mycetomatis to the species level.

Trends of mutation accumulation across global SARS-CoV-2 genomes: Implications for the evolution of the novel coronavirus

To understand SARS-CoV-2 microevolution, this study explored the genome-wide frequency, gene-wise distribution, and molecular nature of all point-mutations detected across its 71,703 RNA-genomes deposited in GISAID till 21 August 2020. Globally, nsp1/nsp2 and orf7a/orf3a were the most mutation-ridden non-structural and structural genes respectively. Phylogeny of 4618 spatiotemporally-representative genomes revealed that entities belonging to the early lineages are mostly spread over Asian countries, including India, whereas the recently-derived lineages are more globally distributed. Of the total 20,163 instances of polymorphism detected across global genomes, 12,594 and 7569 involved transitions and transversions, predominated by cytidine-to-uridine and guanosine-to-uridine conversions, respectively. Positive selection of nonsynonymous mutations (dN/dS >1) in most of the structural, but not the non-structural, genes indicated that SARS-CoV-2 has already harmonized its replication/transcription machineries with the host metabolism, while it is still redefining virulence/transmissibility strategies at the molecular level. Mechanistic bases and evolutionary/pathogenicity-related implications are discussed for the predominant mutation-types.

In vitro evolution of antibody affinity via insertional scanning mutagenesis of an entire antibody variable region

We report a systematic combinatorial exploration of affinity enhancement of antibodies by insertions and deletions (InDels). Transposon-based introduction of InDels via the method TRIAD (transposition-based random insertion and deletion mutagenesis) was used to generate large libraries with random in-frame InDels across the entire single-chain variable fragment gene that were further recombined and screened by ribosome display. Knowledge of potential insertion points from TRIAD libraries formed the basis of exploration of length and sequence diversity of novel insertions by insertional-scanning mutagenesis (InScaM). An overall 256-fold affinity improvement of an anti-IL-13 antibody BAK1 as a result of InDel mutagenesis and combination with known point mutations validates this approach, and suggests that the results of this InDel mutagenesis and conventional exploration of point mutations can synergize to generate antibodies with higher affinity.

Mutational landscape of K-Ras substitutions at 12th position-a systematic molecular dynamics approach

K-Ras is a small GTPase and acts as a molecular switch by recruiting GEFs and GAPs, and alternates between the inert GDP-bound and the dynamic GTP-bound forms. The amino acid at position 12 of K-Ras is a hot spot for oncogenic mutations (G12A, G12C, G12D, G12R, G12S, and G12V), disturbing the active fold of the protein, leading to cancer development. This study aimed to investigate the potential conformational changes induced by these oncogenic mutations at the 12^th position, impairing GAP-mediated GTP hydrolysis. Comprehensive computational tools (iStable, FoldX, SNPeffect, DynaMut, and CUPSAT) were used to evaluate the effect of these six mutations on the stability of wild type K-Ras protein. The docking of GTP with K-Ras was carried out using AutoDock4.2, followed by molecular dynamics simulations. Furthermore, on comparison of binding energies between the wild type K-Ras and the six mutants, we have demonstrated that the G12A and G12V mutants exhibited the strongest binding efficiency compared to the other four mutants. Trajectory analyses of these mutations revealed that G12A encountered the least deviation, fluctuation, intermolecular H-bonds, and compactness compared to the wildtype, which was supported by the lower Gibbs free energy value. Our study investigates the molecular dynamics simulations of the mutant K-Ras forms at the 12^th position, which expects to provide insights about the molecular mechanisms involved in cancer development, and may serve as a platform for targeted therapies against cancer. Communicated by Ramaswamy H. Sarma.

A missense mutation sheds light on a novel structure-function relationship of RANKL

The tumor necrosis factor (TNF)-like core domain of receptor activator of nuclear factor-κB ligand (RANKL) is a functional domain critical for osteoclast differentiation. One of the missense mutations identified in patients with osteoclast-poor autosomal recessive osteopetrosis (ARO) is located in residue methionine 199 that is replaced with lysine (M199K) amid the TNF-like core domain. However, the structure-function relationship of this mutation is not clear. Sequence-based alignment revealed that the fragment containing human M199 is highly conserved and equivalent to M200 in rat. Using site-directed mutagenesis, we generated three recombinant RANKL mutants M200K/A/E (M200s) by replacing the methionine 200 with lysine (M200K), alanine (M200A), and glutamic acid (M200E), representative of distinct physical properties. TRAcP staining and bone pit assay showed that M200s failed to support osteoclast formation and bone resorption, accompanied by impaired osteoclast-related signal transduction. However, no antagonistic effect was found in M200s against wild-type rat RANKL. Analysis of the crystal structure of RANKL predicted that this methionine residue is located within the hydrophobic core of the protein, thus, likely to be crucial for protein folding and stability. Consistently, differential scanning fluorimetry analysis suggested that M200s were less stable. Western blot analysis analyses further revealed impaired RANKL trimerization by M200s. Furthermore, receptor-ligand binding assay displayed interrupted interaction of M200s to its intrinsic receptors. Collectively, our studies revealed the molecular basis of human M199-induced ARO and elucidated the indispensable role of rodent residue M200 (equivalent to human M199) for the RANKL function.

Deciphering Molecular Virulence Mechanism of Mycobacterium tuberculosis Dop isopeptidase Based on Its Sequence-Structure-Function Linkage

The pupylation pathway marks proteins for prokaryotic ubiquitin-like protein (Pup)-proteasomal degradation and survival strategy of mycobacteria inside of the host macrophages. Deamidase of Pup (Dop) plays a central role in the pupylation pathway. It is still a matter of investigation to know the function of Dop in virulence of mycobacterial lineage. Hence, the present study was intended to describe the sequence-structure-function-virulence link of Dop for understanding the molecular virulence mechanism of Mycobacterium tuberculosis H37Rv (Mtb). Phylogenetic analysis of this study indicated that Dop has extensively diverged across the proteasome-harboring bacteria. The functional part of Dop was converged across the pathogenic mycobacterial lineage. The genome-wide analysis pointed out that the pupylation gene locus was identical to each other, but its genome neighborhood differed from species to species. Molecular modeling and dynamic studies proved that the predicted structure of Mtb Dop was energetically stable and low conformational freedom. Moreover, evolutionary constraints in Mtb Dop were intensively analyzed for inferring its sequence-structure-function relationships for the full virulence of Mtb. It indicated that evolutionary optimization was extensively required to stabilize its local structural environment at the side chains of mutable residues. The sequence-structure-function-virulence link of Dop might have retained in Mtb by reordering hydrophobic and hydrogen bonding patterns in the local structural environment. Thus, the results of our study provide a quest to understand the molecular virulence and pathogenesis mechanisms of Mtb during the infection process.

Analyses of displacements resulting from a point mutation in proteins

The effects of a single residue substitution on the protein backbone are frequently quite small and there are many other potential sources of structural variation for protein. We present here a methodology considering different sources of distortions in order to isolate the very effect of the mutation. To validate our methodology, we consider a well-studied family with many single mutants: the human lysozyme. Most of the perturbations are expected to be at the very localisation of the mutation, but in many cases the effects are propagated at long range. We show that the distances between the mutated residue and the 5% most disturbed residues exponentially decreases. One third of the affected residues are in direct contact with the mutated position; the remaining two thirds are potential allosteric effects. We confirm the reliability of the residues identified as significantly perturbed by comparing our results to experimental studies. We confirm with the present method all the previously identified perturbations. This study shows that mutations have long-range impact on protein backbone that can be detected, although the displacement of the affected atoms is small.

Accessing unexplored regions of sequence space in directed enzyme evolution via insertion/deletion mutagenesis

Insertions and deletions (InDels) are frequently observed in natural protein evolution, yet their potential remains untapped in laboratory evolution. Here we introduce a transposon-based mutagenesis approach (TRIAD) to generate libraries of random variants with short in-frame InDels, and screen TRIAD libraries to evolve a promiscuous arylesterase activity in a phosphotriesterase. The evolution exhibits features that differ from previous point mutagenesis campaigns: while the average activity of TRIAD variants is more compromised, a larger proportion has successfully adapted for the activity. Different functional profiles emerge: (i) both strong and weak trade-off between activities are observed; (ii) trade-off is more severe (20- to 35-fold increased kcat/KM in arylesterase with 60-400-fold decreases in phosphotriesterase activity) and (iii) improvements are present in kcat rather than just in KM, suggesting adaptive solutions. These distinct features make TRIAD an alternative to widely used point mutagenesis, accessing functional innovations and traversing unexplored fitness landscape regions.

Predicting the stability of mutant proteins by computational approaches: an overview

A very large number of computational methods to predict the change in thermodynamic stability of proteins due to mutations have been developed during the last 30 years, and many different web servers are currently available. Nevertheless, most of them suffer from severe drawbacks that decrease their general reliability and, consequently, their applicability to different goals such as protein engineering or the predictions of the effects of mutations in genetic diseases. In this review, we have summarized all the main approaches used to develop these tools, with a survey of the web servers currently available. Moreover, we have also reviewed the different assessments made during the years, in order to allow the reader to check directly the different performances of these tools, to select the one that best fits his/her needs, and to help naïve users in finding the best option for their needs.

funtrp: identifying protein positions for variation driven functional tuning

Evaluating the impact of non-synonymous genetic variants is essential for uncovering disease associations and mechanisms of evolution. An in-depth understanding of sequence changes is also fundamental for synthetic protein design and stability assessments. However, the variant effect predictor performance gain observed in recent years has not kept up with the increased complexity of new methods. One likely reason for this might be that most approaches use similar sets of gene and protein features for modeling variant effects, often emphasizing sequence conservation. While high levels of conservation highlight residues essential for protein activity, much of the variation observable in vivo is arguably weaker in its impact, thus requiring evaluation at a higher level of resolution. Here, we describe functionNeutral/Toggle/Rheostatpredictor (funtrp), a novel computational method that categorizes protein positions based on the position-specific expected range of mutational impacts: Neutral (weak/no effects), Rheostat (function-tuning positions), or Toggle (on/off switches). We show that position types do not correlate strongly with familiar protein features such as conservation or protein disorder. We also find that position type distribution varies across different protein functions. Finally, we demonstrate that position types can improve performance of existing variant effect predictors and suggest a way forward for the development of new ones.

Molecular modelling and dynamic simulations of sequestosome 1 (SQSTM1) missense mutations linked to Paget disease of bone

The Paget disease (PDB; OMIM is 167250) is a chronic bone disease caused by pathogenic mutations in Sequestome1/p62 (SQSTM1) gene. This study has aimed to interpret the relationship of PDB linked SQSTM1 mutations with protein structure and its molecular dynamic features. The disease causative missense mutations were initially collected, and then analyzed for their, exonic and domain distribution, impact on secondary and tertiary structures, and their ability on protein-ligand interactions, using a combination of systems biology approaches. Our results show that most PDB linked SQSTM1 missense mutations affect amino acid residues clustered within or near the UBA domain (aa 389-434), which participates in the ubiquitination of substrates. We also report that the majority mutations occurred in α-helices over β-strands but their effects on the secondary structure were mostly neutral. Global tertiary structure deviations were minimal; however, at amino acid residue level minor structural changes were evident. The molecular dynamics simulation analysis showed that both PB1 and UBA domains were under constant structural fluctuations resulting in closed form conformation of SQSMT1 protein structure, when it is bound to PRKCI ligand. We also found salt bridge conformation changes in the UBA domain of SQSTM1 mutants when they bound to the PRKCI interactor protein. This finding suggests the possibility that mutations in SQSTM1 could impair its ability to ubiquitinate the substrates, eventually affecting autophagy and apoptosis, especially in mature osteoclasts. This study presents the additional insight into structure and function relationship between SQSTM1 mutations and PDB pathogenesis. Communicated by Ramaswamy H. Sarma.

Classification of VUS and unclassified variants in BRCA1 BRCT repeats by molecular dynamics simulation

Pathogenic mutation in BRCA1 gene is one of the most penetrant genetic predispositions towards cancer. Identification of the mutation provides important aspect in prevention and treatment of the mutation-caused cancer. Of the large quantity of genetic variants identified in human BRCA1, substantial portion is classified as Variant of Uncertain Significance (VUS) or unclassified variants due to the lack of functional evidence. In this study, we focused on the VUS and unclassified variants in BRCT repeat located at BRCA1 C-terminal. Utilizing the well-determined structure of BRCT repeats, we measured the influence of the variants on the structural conformations of BRCT repeats by using molecular dynamics simulation (MDS) consisting of RMSD (Root-mean-square-deviation), RMSF (Root-mean-square-fluctuations), Rg (Radius of gyration), SASA (Solvent accessible surface area), NH bond (hydrogen bond) and Covariance analysis. Using this approach, we analyzed 131 variants consisting of 89 VUS (Variant of Uncertain Significance) and 42 unclassified variants (unclassifiable by current methods) within BRCT repeats and were able to differentiate them into 78 Deleterious and 53 Tolerated variants. Comparing the results made by the saturation genome editing assay, multiple experimental assays, and BRCA1 reference databases shows that our approach provides high specificity, sensitivity and robust. Our study opens an avenue to classify VUS and unclassified variants in many cancer predisposition genes with known protein structure.

Comprehensive in silico screening and molecular dynamics studies of missense mutations in Sjogren-Larsson syndrome associated with the ALDH3A2 gene

Sjögren-Larsson syndrome (SLS) is an autoimmune disorder inherited in an autosomal recessive pattern. To date, 80 missense mutations have been identified in association with the Aldehyde Dehydrogenase 3 Family Member A2 (ALDH3A2) gene causing SLS. Disruption of the function of ALDH3A2 leads to excessive accumulation of fat in the cells, which interferes with the normal function of protective membranes or materials that are necessary for the body to function normally. We retrieved 54 missense mutations in the ALDH3A2 from the OMIM, UniProt, dbSNP, and HGMD databases that are known to cause SLS. These mutations were examined with various in silico stability tools, which predicted that the mutations p.S308N and p.R423H that are located at the protein-protein interaction domains are the most destabilizing. Furthermore, to determine the atomistic-level differences within the protein-protein interactions owing to mutations, we performed macromolecular simulation (MMS) using GROMACS to validate the motion patterns and dynamic behavior of the biological system. We found that both mutations (p.S380N and p.R423H) had significant effects on the protein-protein interaction and disrupted the dimeric interactions. The computational pipeline provided in this study helps to elucidate the potential structural and functional differences between the ALDH3A2 native and mutant homodimeric proteins, and will pave the way for drug discovery against specific targets in the SLS patients.

In silico prediction of structural changes in human papillomavirus type 16 (HPV16) E6 oncoprotein and its variants

Background

HPV16 infection is one of the main risk factors involved in the development of cervical cancer, mainly due to the high oncogenic potential of the viral proteins E6 and E7, which are involved in the different processes of malignant transformation. There is a broad spectrum of intratypical variation of E6, which is reflected in its high diversity, biological behavior, global distribution and risk of causing cervical cancer. Experimental studies have shown that the intratypical variants of the protein E6 from the European variants (E-G350, E-A176/G350, E-C188/G350) and Asian-American variants (AAa and AAc), are capable of inducing the differential expression of genes involved in the development of cervical cancer.

Results

An in silico analysis was performed to characterize the molecular effects of these variations using the structure of the HPV16 E6 oncoprotein (PDB: 4XR8; chain H) as a template. In particular, we evaluated the 3D structures of the intratypical variants by structural alignment, ERRAT, Ramachandran plots and prediction of protein disorder, which was further validated by molecular dynamics simulations. Our results, in general, showed no significant changes in the protein 3D structure. However, we observed subtle changes in protein physicochemical features and structural disorder in the N- and C-termini.

Conclusions

Our results showed that mutations in the viral oncogene E6 of six high-risk HPV16 variants are effectively neutral and do not cause significant structural changes except slight variations of structural disorder. As structural disorder is involved in rewiring protein-protein interactions, these results suggest a differential pattern of interaction of E6 with the target protein P53 and possibly different patterns of tumor aggressiveness associated with certain types of variants of the E6 oncoprotein.

Electronic supplementary material

The online version of this article (10.1186/s12860-019-0217-0) contains supplementary material, which is available to authorized users.

Structural bioinformatics analysis of variants on GPCR function

G protein-coupled receptors (GPCRs) are key membrane-embedded receptor proteins, with critical roles in cellular signal transduction. In the era of precision medicine, understanding the role of natural variants on GPCR function is critical, especially from a pharmacogenomics viewpoint. Studies involved in mapping variants to GPCR structures are briefly reviewed here. The endocannabinoid system involving the central nervous system (CNS), the human cannabinoid receptor 1 (CB1), is an important drug target and its variability has implications for disease susceptibility and altered drug and pain response. We have carried out a computational study to map deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) to CB1. CB1 mutations were computationally evaluated from neutral to deleterious, and the top twelve deleterious mutations, with structural information, were found to be either close to the ligand binding region or the G-protein binding site. We have mapped these to the active and inactive CB1 X-ray crystallographic structures to correlate variants with available phenotypic information. We have also carried out molecular dynamics simulations to functionally characterize four selected mutants.

Evolution of protein kinase substrate recognition at the active site

Protein kinases catalyse the phosphorylation of target proteins, controlling most cellular processes. The specificity of serine/threonine kinases is partly determined by interactions with a few residues near the phospho-acceptor residue, forming the so-called kinase-substrate motif. Kinases have been extensively duplicated throughout evolution, but little is known about when in time new target motifs have arisen. Here, we show that sequence variation occurring early in the evolution of kinases is dominated by changes in specificity-determining residues. We then analysed kinase specificity models, based on known target sites, observing that specificity has remained mostly unchanged for recent kinase duplications. Finally, analysis of phosphorylation data from a taxonomically broad set of 48 eukaryotic species indicates that most phosphorylation motifs are broadly distributed in eukaryotes but are not present in prokaryotes. Overall, our results suggest that the set of eukaryotes kinase motifs present today was acquired around the time of the eukaryotic last common ancestor and that early expansions of the protein kinase fold rapidly explored the space of possible target motifs.

Effects of Pro1266Leu mutation on structure and function of glycoprotein Ib binding domain of von Willebrand factor

Glycoprotein Ibα (GpIbα) binding ability of A1 domain of von Willebrand factor (vWF) facilitates platelet adhesion that plays a crucial role in maintaining hemostasis and thrombosis at the site of vascular damage. There are both "loss as well as gain of function" mutations observed in this domain. Naturally occurring "gain of function" mutations leave self-activating impacts on the A1 domain which turns the normal binding to characteristic constitutive binding with GPIbα. These "gain of function" mutations are associated with the von Willebrand disease type 2B. In recent years, studies focused on understanding the mechanism and conformational patterns attached to these phenomena have been conducted, but the conformational pathways leading to such binding patterns are poorly understood as of now. To obtain a microscopic picture of such events for the better understanding of pathways, we used molecular dynamics (MD) simulations along with principal component analysis and normal mode analysis to study the effects of Pro1266Leu (Pro503Leu in structural context) mutation on the structure and function of A1 domain of vWF. MD simulations have provided atomic-level details of intermolecular motions as a function of time to understand the dynamic behavior of A1 domain of vWF. Comparative analysis of the trajectories obtained from MD simulations of both the wild type and Pro503Leu mutant suggesting appreciable conformational changes in the structure of mutant which might provide a basis for assuming the "gain of function" effects of these mutations on the A1 domain of vWF, resulting in the constitutive binding with GpIbα.

Structural and dynamical description of the enzymatic reaction of a phosphohexomutase

Enzymes are known to adopt various conformations at different points along their catalytic cycles. Here, we present a comprehensive analysis of 15 isomorphous, high resolution crystal structures of the enzyme phosphoglucomutase from the bacterium Xanthomonas citri. The protein was captured in distinct states critical to function, including enzyme-substrate, enzyme-product, and enzyme-intermediate complexes. Key residues in ligand recognition and regions undergoing conformational change are identified and correlated with the various steps of the catalytic reaction. In addition, we use principal component analysis to examine various subsets of these structures with two goals: (1) identifying sites of conformational heterogeneity through a comparison of room temperature and cryogenic structures of the apo-enzyme and (2) a priori clustering of the enzyme-ligand complexes into functionally related groups, showing sensitivity of this method to structural features difficult to detect by traditional methods. This study captures, in a single system, the structural basis of diverse substrate recognition, the subtle impact of covalent modification, and the role of ligand-induced conformational change in this representative enzyme of the α-D-phosphohexomutase superfamily.

Genetic and transcriptional variations in NRAMP-2 and OPAQUE1 genes are associated with salt stress response in wheat

SNP alleles on chromosomes 4BL and 6AL are associated with sensitivity to salt tolerance in wheat and upon validation can be exploited in the development of salt-tolerant wheat varieties. The dissection of the genetic and molecular components of salt stress response offers strong opportunities toward understanding and improving salt tolerance in crops. In this study, GWAS was employed to identify a total of 106 SNP loci (R² = 0.12-63.44%) linked to salt stress response in wheat using leaf chlorophyll fluorescence, grain quality and shoot ionic (Na⁺ and K⁺ ions) attributes. Among them, 14 SNP loci individually conferred pleiotropic effects on multiple independent salinity tolerance traits including loci at 99.04 cM (R² ≥ 14.7%) and 68.45 cM (R² ≥ 4.10%) on chromosomes 6AL and 4BL, respectively, that influenced shoot Na⁺-uptake, shoot K⁺/Na⁺ ratio, and specific energy fluxes for absorption (ABS/RC) and dissipation (DIo/RC). Analysis of the open reading frame (ORF) containing the SNP markers revealed that they are orthologous to genes involved in photosynthesis and plant stress (salt) response. Further transcript abundance and qRT-PCR analyses indicated that the genes are mostly up-regulated in salt-tolerant and down-regulated in salt-sensitive wheat genotypes including NRAMP-2 and OPAQUE1 genes on 4BL and 6AL, respectively. Both genes showed highest differential expression between contrasting genotypes when expressions of all the genes within their genetic intervals were analyzed. Possible cis-acting regulatory elements and coding sequence variation that may be involved in salt stress response were also identified in both genes. This study identified genetic and molecular components of salt stress response that are associated with Na⁺-uptake, shoot Na⁺/K⁺ ratio, ABS/RC, DIo/RC, and grain quality traits and upon functional validation would facilitate the development of gene-specific markers that could be deployed to improve salinity tolerance in wheat.

Structural Perspective on Revealing and Altering Molecular Functions of Genetic Variants Linked with Diseases

Structural information of biological macromolecules is crucial and necessary to deliver predictions about the effects of mutations-whether polymorphic or deleterious (i.e., disease causing), wherein, thermodynamic parameters, namely, folding and binding free energies potentially serve as effective biomarkers. It may be emphasized that the effect of a mutation depends on various factors, including the type of protein (globular, membrane or intrinsically disordered protein) and the structural context in which it occurs. Such information may positively aid drug-design. Furthermore, due to the intrinsic plasticity of proteins, even mutations involving radical change of the structural and physico⁻chemical properties of the amino acids (native vs. mutant) can still have minimal effects on protein thermodynamics. However, if a mutation causes significant perturbation by either folding or binding free energies, it is quite likely to be deleterious. Mitigating such effects is a promising alternative to the traditional approaches of designing inhibitors. This can be done by structure-based in silico screening of small molecules for which binding to the dysfunctional protein restores its wild type thermodynamics. In this review we emphasize the effects of mutations on two important biophysical properties, stability and binding affinity, and how structures can be used for structure-based drug design to mitigate the effects of disease-causing variants on the above biophysical properties.

Concordant Changes in Gene Expression and Nucleotides Underlie Independent Adaptation to Hydrogen-Sulfide-Rich Environments

The colonization of novel environments often involves changes in gene expression, protein coding sequence, or both. Studies of how populations adapt to novel conditions, however, often focus on only one of these two processes, potentially missing out on the relative importance of different parts of the evolutionary process. In this study, our objectives were 1) to better understand the qualitative concordance between conclusions drawn from analyses of differential expression and changes in genic sequence and 2) to quantitatively test whether differentially expressed genes were enriched for sites putatively under positive selection within gene regions. To achieve this, we compared populations of fish (Poecilia mexicana) that have independently adapted to hydrogen-sulfide-rich environments in southern Mexico to adjacent populations residing in nonsulfidic waters. Specifically, we used RNA-sequencing data to compare both gene expression and DNA sequence differences between populations. Analyzing these two different data types led to similar conclusions about which biochemical pathways (sulfide detoxification and cellular respiration) were involved in adaptation to sulfidic environments. Additionally, we found a greater overlap between genes putatively under selection and differentially expressed genes than expected by chance. We conclude that considering both differential expression and changes in DNA sequence led to a more comprehensive understanding of how these populations adapted to extreme environmental conditions. Our results imply that changes in both gene expression and DNA sequence-sometimes at the same loci-may be involved in adaptation.

Contributions of substitutions and indels to the structural variations in ancient protein superfamilies

Background

Quantitative evaluation of protein structural evolution is important for our understanding of protein biological functions and their evolutionary adaptation, and is useful in guiding protein engineering. However, compared to the models for sequence evolution, the quantitative models for protein structural evolution received less attention. Ancient protein superfamilies are often considered versatile, allowing genetic and functional diversifications during long-term evolution. In this study, we investigated the quantitative impacts of sequence variations on the structural evolution of homologues in 68 ancient protein superfamilies that exist widely in sequenced eukaryotic, bacterial and archaeal genomes.

Results

We found that the accumulated structural variations within ancient superfamilies could be explained largely by a bilinear model that simultaneously considers amino acid substitution and insertion/deletion (indel). Both substitutions and indels are essential for explaining the structural variations within ancient superfamilies. For those ancient superfamilies with high bilinear multiple correlation coefficients, the influence of each unit of substitution or indel on structural variations is almost constant within each superfamily, but varies greatly among different superfamilies. The influence of each unit indel on structural variations is always larger than that of each unit substitution within each superfamily, but the accumulated contributions of indels to structural variations are lower than those of substitutions in most superfamilies. The total contributions of sequence indels and substitutions (46% and 54%, respectively) to the structural variations that result from sequence variations are slightly different in ancient superfamilies.

Conclusions

Structural variations within ancient protein superfamilies accumulated under the significantly bilinear influence of amino acid substitutions and indels in sequences. Both substitutions and indels are essential for explaining the structural variations within ancient superfamilies. For those structural variations resulting from sequence variations, the total contribution of indels is slightly lower than that of amino acid substitutions. The regular clock exists not only in protein sequences, but also probably in protein structures.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5178-8) contains supplementary material, which is available to authorized users.

Topology of active, membrane-embedded Bax in the context of a toroidal pore

Bax is a Bcl-2 protein critical for apoptosis induction. In healthy cells, Bax is mostly a monomeric, cytosolic protein, while upon apoptosis initiation it inserts into the outer mitochondrial membrane, oligomerizes, and forms pores that release proapoptotic factors like Cytochrome c into the cytosol. The structures of active Bax and its homolog Bak are only partially understood and the topology of the proteins with respect to the membrane bilayer is controversially described in the literature. Here, we systematically review and examine the protein-membrane, protein-water, and protein-protein contacts of the nine helices of active Bax and Bak, and add a new set of topology data obtained by fluorescence and EPR methods. We conclude based on the consistent part of the datasets that the core/dimerization domain of Bax (Bak) is water exposed with only helices 4 and 5 in membrane contact, whereas the piercing/latch domain is in peripheral membrane contact, with helix 9 being transmembrane. Among the available structural models, those considering the dimerization/core domain at the rim of a toroidal pore are the most plausible to describe the active state of the proteins, although the structural flexibility of the piercing/latch domain does not allow unambiguous discrimination between the existing models.

Trade-offs with stability modulate innate and mutationally acquired drug resistance in bacterial dihydrofolate reductase enzymes

Structural stability is a major constraint on the evolution of protein sequences. However, under strong directional selection, mutations that confer novel phenotypes but compromise structural stability of proteins may be permissible. During the evolution of antibiotic resistance, mutations that confer drug resistance often have pleiotropic effects on the structure and function of antibiotic-target proteins, usually essential metabolic enzymes. In the present study, we show that trimethoprim (TMP)-resistant alleles of dihydrofolate reductase from Escherichia coli (EcDHFR) harboring the Trp30Gly, Trp30Arg or Trp30Cys mutations are significantly less stable than the wild-type, making them prone to aggregation and proteolysis. This destabilization is associated with a lower expression level, resulting in a fitness cost and negative epistasis with other TMP-resistant mutations in EcDHFR. Using structure-based mutational analysis, we show that perturbation of critical stabilizing hydrophobic interactions in wild-type EcDHFR enzyme explains the phenotypes of Trp30 mutants. Surprisingly, though crucial for the stability of EcDHFR, significant sequence variation is found at this site among bacterial dihydrofolate reductases (DHFRs). Mutational and computational analyses in EcDHFR and in DHFR enzymes from Staphylococcus aureus and Mycobacterium tuberculosis demonstrate that natural variation at this site and its interacting hydrophobic residues modulates TMP resistance in other bacterial DHFRs as well, and may explain the different susceptibilities of bacterial pathogens to TMP. Our study demonstrates that trade-offs between structural stability and function can influence innate drug resistance as well as the potential for mutationally acquired drug resistance of an enzyme.

Computational Design To Reduce Conformational Flexibility and Aggregation Rates of an Antibody Fab Fragment

Computationally guided semirational design has significant potential for improving the aggregation kinetics of protein biopharmaceuticals. While improvement in the global conformational stability can stabilize proteins to aggregation under some conditions, previous studies suggest that such an approach is limited, because thermal transition temperatures ( T_m) and the fraction of protein unfolded ( f_T) tend to only correlate with aggregation kinetics where the protein is incubated at temperatures approaching the T_m. This is because under these conditions, aggregation from globally unfolded protein becomes dominant. However, under native conditions, the aggregation kinetics are presumed to be dependent on local structural fluctuations or partial unfolding of the native state, which reveal regions of high propensity to form protein-protein interactions that lead to aggregation. In this work, we have targeted the design of stabilizing mutations to regions of the A33 Fab surface structure, which were predicted to be more flexible. This Fab already has high global stability, and global unfolding is not the main cause of aggregation under most conditions. Therefore, the aim was to reduce the conformational flexibility and entropy of the native protein at various locations and thus identify which of those regions has the greatest influence on the aggregation kinetics. Highly dynamic regions of structure were identified through both molecular dynamics simulation and B-factor analysis of related X-ray crystal structures. The most flexible residues were mutated into more stable variants, as predicted by Rosetta, which evaluates the ΔΔ G_ND for each potential point mutation. Additional destabilizing variants were prepared as controls to evaluate the prediction accuracy and also to assess the general influence of conformational stability on aggregation kinetics. The thermal conformational stability, and aggregation rates of 18 variants at 65 °C, were each examined at pH 4, 200 mM ionic strength, under which conditions the initial wild-type protein was <5% unfolded. Variants with decreased T_m values led to more rapid aggregation due to an increase in the fraction of protein unfolded under the conditions studied. As expected, no significant improvements were observed in the global conformational stability as measured by T_m. However, 6 of the 12 stable variants led to an increase in the cooperativity of unfolding, consistent with lower conformational flexibility and entropy in the native ensemble. Three of these had 5-11% lower aggregation rates, and their structural clustering indicated that the local dynamics of the C-terminus of the heavy chain had a role in influencing the aggregation rate.

The Staphylococcus aureus Two-Component System AgrAC Displays Four Distinct Genomic Arrangements That Delineate Genomic Virulence Factor Signatures

Two-component systems (TCSs) consist of a histidine kinase and a response regulator. Here, we evaluated the conservation of the AgrAC TCS among 149 completely sequenced Staphylococcus aureus strains. It is composed of four genes: agrBDCA. We found that: (i) AgrAC system (agr) was found in all but one of the 149 strains, (ii) the agr positive strains were further classified into four agr types based on AgrD protein sequences, (iii) the four agr types not only specified the chromosomal arrangement of the agr genes but also the sequence divergence of AgrC histidine kinase protein, which confers signal specificity, (iv) the sequence divergence was reflected in distinct structural properties especially in the transmembrane region and second extracellular binding domain, and (v) there was a strong correlation between the agr type and the virulence genomic profile of the organism. Taken together, these results demonstrate that bioinformatic analysis of the agr locus leads to a classification system that correlates with the presence of virulence factors and protein structural properties.

Labeling for Quantitative Comparison of Imaging Measurements in Vitro and in Cells

Qualitative imaging of biomolecular localization and distribution inside cells has revolutionized cell biology. Most of these powerful techniques require modifications to the target biomolecule. Over the past 10 years, these techniques have been extended to quantitative measurements, from in-cell protein folding rates to complex dissociation constants to RNA lifetimes. Such measurements can be affected even when a target molecule is just mildly perturbed by its labels. Here, the impact of labeling on protein (and RNA) structure, stability, and function in cells is discussed via practical examples from the recent literature. General guidelines for selecting and validating modification sites are provided to bring the best from cell biology and imaging to quantitative biophysical experiments inside cells.

Site-Mutation of Hydrophobic Core Residues Synchronically Poise Super Interleukin 2 for Signaling: Identifying Distant Structural Effects through Affordable Computations

A superkine variant of interleukin-2 with six site mutations away from the binding interface developed from the yeast display technique has been previously characterized as undergoing a distal structure alteration which is responsible for its super-potency and provides an elegant case study with which to get insight about how to utilize allosteric effect to achieve desirable protein functions. By examining the dynamic network and the allosteric pathways related to those mutated residues using various computational approaches, we found that nanosecond time scale all-atom molecular dynamics simulations can identify the dynamic network as efficient as an ensemble algorithm. The differentiated pathways for the six core residues form a dynamic network that outlines the area of structure alteration. The results offer potentials of using affordable computing power to predict allosteric structure of mutants in knowledge-based mutagenesis.

Atomic mutagenesis of stop codon nucleotides reveals the chemical prerequisites for release factor-mediated peptide release

Termination of protein synthesis is triggered by the recognition of a stop codon at the ribosomal A site and is mediated by class I release factors (RFs). Whereas in bacteria, RF1 and RF2 promote termination at UAA/UAG and UAA/UGA stop codons, respectively, eukaryotes only depend on one RF (eRF1) to initiate peptide release at all three stop codons. Based on several structural as well as biochemical studies, interactions between mRNA, tRNA, and rRNA have been proposed to be required for stop codon recognition. In this study, the influence of these interactions was investigated by using chemically modified stop codons. Single functional groups within stop codon nucleotides were substituted to weaken or completely eliminate specific interactions between the respective mRNA and RFs. Our findings provide detailed insight into the recognition mode of bacterial and eukaryotic RFs, thereby revealing the chemical groups of nucleotides that define the identity of stop codons and provide the means to discriminate against noncognate stop codons or UGG sense codons.

In silico identification of rescue sites by double force scanning

Motivation

A deleterious amino acid change in a protein can be compensated by a second-site rescue mutation. These compensatory mechanisms can be mimicked by drugs. In particular, the location of rescue mutations can be used to identify protein regions that can be targeted by small molecules to reactivate a damaged mutant.

Results

We present the first general computational method to detect rescue sites. By mimicking the effect of mutations through the application of forces, the double force scanning (DFS) method identifies the second-site residues that make the protein structure most resilient to the effect of pathogenic mutations. We tested DFS predictions against two datasets containing experimentally validated and putative evolutionary-related rescue sites. A remarkably good agreement was found between predictions and experimental data. Indeed, almost half of the rescue sites in p53 was correctly predicted by DFS, with 65% of remaining sites in contact with DFS predictions. Similar results were found for other proteins in the evolutionary dataset.

Availability and implementation

The DFS code is available under GPL at https://fornililab.github.io/dfs/.

Supplementary information

Supplementary data are available at Bioinformatics online.

The Hypertension Pandemic: An Evolutionary Perspective

Hypertension affects over 1.2 billion individuals worldwide and has become the most critical and expensive public health problem. Hypertension is a multifactorial disease involving environmental and genetic factors together with risk-conferring behaviors. The cause of the disease is identified in ∼10% of the cases (secondary hypertension), but in 90% of the cases no etiology is found (primary or essential hypertension). For this reason, a better understanding of the mechanisms controlling blood pressure in normal and hypertensive patients is the aim of very active experimental and clinical research. In this article, we review the importance of the renin-angiotensin-aldosterone system (RAAS) for the control of blood pressure, focusing on the evolution of the system and its critical importance for adaptation of vertebrates to a terrestrial and dry environment. The evolution of blood pressure control during the evolution of primates, hominins, and humans is discussed, together with the role of common genetic factors and the possible causes of the current hypertension pandemic in the light of evolutionary medicine.

Global substitution of hemeproteins with noncanonical amino acids in Escherichia coli with intact cofactor maturation machinery

Global substitution of canonical amino acids (cAAs) with noncanonical (ncAAs) counterparts in proteins whose function is dependent on post-translational events such as cofactor binding is still a methodically challenging and difficult task as ncAA insertion generally interferes with the cofactor biosynthesis machinery. Here, we report a technology for the expression of fully substituted and functionally active cofactor-containing hemeproteins. The maturation process which yields an intact cofactor is timely separated from cAA→ncAA substitutions. This is achieved by an optimised expression and fermentation procedure which includes pre-induction of the heme cofactor biosynthesis followed by an incorporation experiment at multiple positions in the protein sequence. This simple strategy can be potentially applied for engineering of other cofactor-containing enzymes.

Associations of ERAP1 coding variants and domain specific interaction with HLA-C∗06 in the early onset psoriasis patients of India

Interferon-γ-induced aminopeptidase ERAP1 trims peptides within the endoplasmic reticulum so that they can be loaded onto MHC class I and presented to the CD8⁺ T-cells. ERAP1 association and its interaction with HLA-C∗06 is controversial across different populations. We have investigated the association and possible functional role of non-synonymous SNPs at different exons of ERAP1 (rs26653: Arg127Pro, rs30187: Lys528Arg and rs27044: Gln730Glu) and their interactions with HLA-C∗06 in psoriasis. Significant associations of HLA-C∗06 (OR=5.47, P<2.2×10^-16), rs30187 (OR 1.35, P=7.4×10^-4) and rs27044 (OR=1.24, P=5.8×10^-3) were observed. All three ERAP1 SNPs showed significant association only for HLA-C∗06 positive patients, while rs30187 and rs27044 showed significant association only for early onset patients (rs30187: OR=1.47, P=9.6×10^-5; rs27044: OR=1.36, P=3.3×10^-4). No differential expression of ERAP1 was observed either between paired uninvolved and involved skin tissues of psoriasis patients or between non-risk and risk variants in the involved skin. Significant epistatic interaction was observed between HLA-C∗06 and the SNP (rs27044) located at the peptide-binding cavity of ERAP1. Evolutionary conservation analysis among mammals showed confinement of Lys528 and Gln730 within highly conserved regions of ERAP1 and suggested the possible detrimental effect of this allele in ERAP1 regulation.

Effects of short indels on protein structure and function in human genomes

Insertions and deletions (indels) represent the second most common type of genetic variations in human genomes. Indels can be deleterious and contribute to disease susceptibility as recent genome sequencing projects revealed a large number of indels in various cancer types. In this study, we investigated the possible effects of small coding indels on protein structure and function, and the baseline characteristics of indels in 2504 individuals of 26 populations from the 1000 Genomes Project. We found that each population has a distinct pattern in genes with small indels. Frameshift (FS) indels are enriched in olfactory receptor activity while non-frameshift (NFS) indels are enriched in transcription-related proteins. Structural analysis of NFS indels revealed that they predominantly adopt coil or disordered conformations, especially in proteins with transcription-related NFS indels. These results suggest that the annotated coding indels from the 1000 Genomes Project, while contributing to genetic variations and phenotypic diversity, generally do not affect the core protein structures and have no deleterious effect on essential biological processes. In addition, we found that a number of reference genome annotations might need to be updated due to the high prevalence of annotated homozygous indels in the general population.

Computational analysis for GNAQ mutations: New insights on the molecular etiology of Sturge-Weber syndrome

Somatic activating mutations in the GNAQ have been recently associated with several congenital genetic disorders and tumors; however, the molecular mechanism/etiology that leads to GNAQ somatic mosaic mutation are unknown. Here, we reported a case of Sturge-Weber Syndrome (SWS) manifesting cutaneous vascular malformations (hemifacial Port-wine stain), cerebral and ocular vascular abnormalities (including epilepsy and glaucoma) and harboring a c.548G>A (p.R183Q) somatic mosaic mutation in GNAQ. Computational modeling studies were performed to assistant with the comprehension of the functional impact of p.R183Q and p.Q209L mutations in GNAQ, which encodes a G protein subunit alpha q (Gαq). The p.R183Q mutation was predicted to abolish hydrogen bonds between R183 residue and GDP molecule, destabilizing the inactive GDP-bound conformation of the Gαq mutants. Furthermore, replacement of R183 by Q183 residue was predicted to promote conformation changes in protein surface features affecting the switch I region, a key region that undergoes conformational changes triggered by receptor binding during signal transduction. In addition, replacement of Q209 by L209 residue was predicted to affect the molecular interaction between Gαq and Gβ subunit, impairing formation of the inactive heterotrimeric complex. These findings, in association with PPI network analysis, indicate that p.R183Q and p.Q209L mutations result in the over-activation of different downstream effectors, which in turn will determine the distinct cell responses and phenotype. These findings bring new insights on molecular etiology of vascular malformations associated to SWS and on different mechanisms underlying hyperactivation of downstream pathways to Gαq.

Biomolecular engineering for nanobio/bionanotechnology

Biomolecular engineering can be used to purposefully manipulate biomolecules, such as peptides, proteins, nucleic acids and lipids, within the framework of the relations among their structures, functions and properties, as well as their applicability to such areas as developing novel biomaterials, biosensing, bioimaging, and clinical diagnostics and therapeutics. Nanotechnology can also be used to design and tune the sizes, shapes, properties and functionality of nanomaterials. As such, there are considerable overlaps between nanotechnology and biomolecular engineering, in that both are concerned with the structure and behavior of materials on the nanometer scale or smaller. Therefore, in combination with nanotechnology, biomolecular engineering is expected to open up new fields of nanobio/bionanotechnology and to contribute to the development of novel nanobiomaterials, nanobiodevices and nanobiosystems. This review highlights recent studies using engineered biological molecules (e.g., oligonucleotides, peptides, proteins, enzymes, polysaccharides, lipids, biological cofactors and ligands) combined with functional nanomaterials in nanobio/bionanotechnology applications, including therapeutics, diagnostics, biosensing, bioanalysis and biocatalysts. Furthermore, this review focuses on five areas of recent advances in biomolecular engineering: (a) nucleic acid engineering, (b) gene engineering, (c) protein engineering, (d) chemical and enzymatic conjugation technologies, and (e) linker engineering. Precisely engineered nanobiomaterials, nanobiodevices and nanobiosystems are anticipated to emerge as next-generation platforms for bioelectronics, biosensors, biocatalysts, molecular imaging modalities, biological actuators, and biomedical applications.

TK, C. GPD, R. S, Zayed H. Determining the role of missense mutations in the POU domain of HNF1A that reduce the DNA-binding affinity: A computational approach

Maturity-onset diabetes of the young type 3 (MODY3) is a non-ketotic form of diabetes associated with poor insulin secretion. Over the past years, several studies have reported the association of missense mutations in the Hepatocyte Nuclear Factor 1 Alpha (HNF1A) with MODY3. Missense mutations in the POU homeodomain (POUH) of HNF1A hinder binding to the DNA, thereby leading to a dysfunctional protein. Missense mutations of the HNF1A were retrieved from public databases and subjected to a three-step computational mutational analysis to identify the underlying mechanism. First, the pathogenicity and stability of the mutations were analyzed to determine whether they alter protein structure and function. Second, the sequence conservation and DNA-binding sites of the mutant positions were assessed; as HNF1A protein is a transcription factor. Finally, the biochemical properties of the biological system were validated using molecular dynamic simulations in Gromacs 4.6.3 package. Two arginine residues (131 and 203) in the HNF1A protein are highly conserved residues and contribute to the function of the protein. Furthermore, the R131W, R131Q, and R203C mutations were predicted to be highly deleterious by in silico tools and showed lower binding affinity with DNA when compared to the native protein using the molecular docking analysis. Triplicate runs of molecular dynamic (MD) simulations (50ns) revealed smaller changes in patterns of deviation, fluctuation, and compactness, in complexes containing the R131Q and R131W mutations, compared to complexes containing the R203C mutant complex. We observed reduction in the number of intermolecular hydrogen bonds, compactness, and electrostatic potential, as well as the loss of salt bridges, in the R203C mutant complex. Substitution of arginine with cysteine at position 203 decreases the affinity of the protein for DNA, thereby destabilizing the protein. Based on our current findings, the MD approach is an important tool for elucidating the impact and affinity of mutations in DNA-protein interactions and understanding their function.

Infectious bronchitis vaccine virus detection and part-S1 genetic variation following single or dual inoculation in broiler chicks

An investigation was undertaken of the extent of genetic variation occurring within infectious bronchitis virus (IBV) vaccine strains following vaccination of day-old broiler chicks. Chicks were divided into seven groups, with two groups receiving single Massachusetts (Mass) vaccinations while the other four were inoculated with combinations of different IBV serotypes; Mass, 793B, D274 and Arkansas (Ark). The remaining group was maintained as an unvaccinated control. Following vaccination, swabs and tissues collected at intervals were pooled and RNA was extracted for detection of IBV by reverse transcription polymerase chain reaction. Positive amplicons were sequenced for the part-S1 gene and compared to the original vaccine strain sequences. Single nucleotide polymorphisms, amino acid variations and hydrophobicity changes were identified and recorded for each sampling point. A total of 106 single nucleotide polymorphisms were detected within 28 isolates. The average single nucleotide polymorphism counts of swab isolates were greater than those found in tissue samples. This translated into 64 amino acid changes; however only six resulted in a change to the hydrophobicity properties. All hydrophobic alterations occurred within swab isolates and the majority were recovered at 3 days post vaccination suggesting such changes to be detrimental to early virus survival. Nucleotide deletions were seen only in the group given the combination of Mass and Ark. Of the 16 sequenced samples in this group, 13 contained the same AAT deletion at position 1033 1035 in the Ark strains. Findings presented in this study demonstrate alteration in the S1 nucleotide sequence following co-administration of live IBV vaccines.

Site-selective tagging of proteins by pnictogen-mediated self-assembly

Trivalent pnictogens (Pn) enable the selective self-assembly between an engineered di-cysteine motif in a protein and a thiol-containing lanthanide (Ln) probe.

Chromosomal rearrangements and protein globularity changes in Mycobacterium tuberculosis isolates from cerebrospinal fluid

Background

Meningitis is a major cause of mortality in tuberculosis (TB). It is not clear what factors promote central nervous system invasion and pathology but it has been reported that certain strains of Mycobacterium tuberculosis (Mtb) might have genetic traits associated with neurotropism.

Methods

In this study, we generated whole genome sequences of eight clinical strains of Mtb that were isolated from the cerebrospinal fluid (CSF) of patients presenting with tuberculous meningitis (TBM) in Malaysia, and compared them to the genomes of H37Rv and other respiratory Mtb genomes either downloaded from public databases or extracted from local sputum isolates. We aimed to find genomic features that might be distinctly different between CSF-derived and respiratory Mtb.

Results

Genome-wide comparisons revealed rearrangements (translocations, inversions, insertions and deletions) and non-synonymous SNPs in our CSF-derived strains that were not observed in the respiratory Mtb genomes used for comparison. These rearranged segments were rich in genes for PE (proline-glutamate)/PPE (proline-proline-glutamate), transcriptional and membrane proteins. Similarly, most of the ns SNPs common in CSF strains were noted in genes encoding PE/PPE proteins. Protein globularity differences were observed among mycobacteria from CSF and respiratory sources and in proteins previously reported to be associated with TB meningitis. Transcription factors and other transcription regulators featured prominently in these proteins. Homologs of proteins associated with Streptococcus pneumoniae meningitis and Neisseria meningitidis virulence were identified in neuropathogenic as well as respiratory mycobacterial spp. examined in this study.

Discussion

The occurrence of in silico genetic differences in CSF-derived but not respiratory Mtb suggests their possible involvement in the pathogenesis of TBM. However, overall findings in this comparative analysis support the postulation that TB meningeal infection is more likely to be related to the expression of multiple virulence factors on interaction with host defences than to CNS tropism associated with specific genetic traits.