BackgroundWhite Spot Syndrome Virus (WSSV) is an enveloped double-stranded DNA virus which causes mortality of several species of shrimp, being considered one of the main pathogens that affects global shrimp farming. This virus presents a complex genome of 300kb and viral isolates that present genomes with great identity. Despite this conservation, some variable regions in the WSSV genome occur in coding regions, and these putative proteins may have some relationship with viral adaptation and virulence mechanisms. Until now, the functions of these proteins were little studied. In this work, sequences and putative proteins encoded by WSSV variable regions were characterized in silico.ResultsThe in silico approach enabled determining the variability of some sequences, as well as the identification of some domains resembling the Formin homology 2, RNA recognition motif, Xeroderma pigmentosum group D repair helicase, Hemagglutinin and Ankyrin motif. The information obtained from the sequences and the analysis of secondary and tertiary structure models allow to infer that some of these proteins possibly have functions related to protein modulation/degradation, intracellular transport, recombination and endosome fusion events.ConclusionsThe bioinformatics approaches were efficient in generating three-dimensional models and to identify domains, thereby enabling to propose possible functions for the putative polypeptides produced by the ORFs wsv129, wsv178, wsv249, wsv463a, wsv477, wsv479, wsv492, and wsv497.
Plasmodium vivax merozoite surface protein 3α (PvMSP3α) is a promising vaccine candidate which has shown strong association with immunogenicity and protectiveness. Its use is however complicated by evolutionary plasticity features which enhance immune evasion. Low complexity regions (LCRs) provide plasticity in surface proteins of Plasmodium species, but its implication in vaccine design remain unexplored. Here population genetic, comparative phylogenetic and structural biology analysis was performed on the gene encoding PvMSP3α. Three LCRs were found in PvMSP3α block II. Both the predicted tertiary structure of the protein and the phylogenetic trees based on this region were influenced by the presence of the LCRs. The LCRs were mainly B cell epitopes within or adjacent. In addition a repeat motif mimicking one of the B cell epitopes was found within the PvMSP3a block II low complexity region. This particular B cell epitope also featured rampant alanine substitutions which might impair antibody binding. The findings indicate that PvMSP3α block II possesses LCRs which might confer a strong phenotypic plasticity. The phenomenon of phenotypic plasticity and implication of LCRs in malaria immunology in general and vaccine candidate genes in particular merits further exploration.
Background Computational models of RNA 3D structure often present various inaccuracies caused by simplifications used in structure prediction methods, such as template-based modeling or coarse-grained simulations. To obtain a high-quality model, the preliminary RNA structural model needs to be refined, taking into account atomic interactions. The goal of the refinement is not only to improve the local quality of the model but to bring it globally closer to the true structure. Results We present QRNAS, a software tool for fine-grained refinement of nucleic acid structures, which is an extension of the AMBER simulation method with additional restraints. QRNAS is capable of handling RNA, DNA, chimeras, and hybrids thereof, and enables modeling of nucleic acids containing modified residues. Conclusions We demonstrate the ability of QRNAS to improve the quality of models generated with different methods. QRNAS was able to improve MolProbity scores of NMR structures, as well as of computational models generated in the course of the RNA-Puzzles experiment. The overall geometry improvement may be associated with increased model accuracy, especially on the level of correctly modeled base-pairs, but the systematic improvement of root mean square deviation to the reference structure should not be expected. The method has been integrated into a computational modeling workflow, enabling improved RNA 3D structure prediction.
BackgroundThe THAP (Thanatos Associated Proteins) protein family in humans is implicated in various important cellular processes like epigenetic regulation, maintenance of pluripotency, transposition and disorders like cancers and hemophilia. The human THAP protein family which consists of twelve members of different lengths has a well characterized amino terminal, zinc-coordinating, DNA-binding domain called the THAP domain. However, the carboxy terminus of most THAP proteins is yet to be structurally characterized. A coiled coil region is known to help in protein oligomerization in THAP1 and THAP11. It is not known if other human THAP proteins oligomerize. We have used bioinformatic tools to explore the possibility of dimerization of THAP proteins via a coiled coil region.ResultsClassification of human THAP protein into three size based groups led to the identification of an evolutionarily conserved alpha helical region, downstream of the amino terminal THAP domain. Secondary structure predictions, alpha helical wheel plots and protein models demonstrated the strong possibility of coiled coil formation in this conserved, leucine rich region of all THAP proteins except THAP10.ConclusionsThe identification of a predicted oligomerization region in the human THAP protein family opens new directions to investigate the members of this protein family.
Background It had long been thought that a protein exhibits its specific function through its own specific 3D-structure under physiological conditions. However, subsequent research has shown that there are many proteins without specific 3D-structures under physiological conditions, so-called intrinsically disordered proteins (IDPs). This study presents a new technique for predicting intrinsically disordered regions in a protein, based on our average distance map (ADM) technique. The ADM technique was developed to predict compact regions or structural domains in a protein. In a protein containing partially disordered regions, a domain region is likely to be ordered, thus it is unlikely that a disordered region would be part of any domain. Therefore, the ADM technique is expected to also predict a disordered region between domains. Results The results of our new technique are comparable to the top three performing techniques in the community-wide CASP10 experiment. We further discuss the case of p53, a tumor-suppressor protein, which is the most significant protein among cell cycle regulatory proteins. This protein exhibits a disordered character as a monomer but an ordered character when two p53s form a dimer. Conclusion Our technique can predict the location of an intrinsically disordered region in a protein with an accuracy comparable to the best techniques proposed so far. Furthermore, it can also predict a core region of IDPs forming definite 3D structures through interactions, such as dimerization. The technique in our study may also serve as a means of predicting a disordered region which would become an ordered structure when binding to another protein.
BackgroundThree transketolase genes have been identified in the human genome to date: transketolase (TKT), transketolase-like 1 (TKTL1) and transketolase-like 2 (TKTL2). Altered TKT functionality is strongly implicated in the development of diabetes and various cancers, thus offering possible therapeutic utility. It will be of great value to know whether TKTL1 and TKTL2 are, similarly, potential therapeutic targets. However, it remains unclear whether TKTL1 and TKTL2 are functional transketolases.ResultsHomology modelling of TKTL1 and TKTL2 using TKT as template, revealed that both TKTL1 and TKTL2 could assume a folded structure like TKT. TKTL1/2 presented a cleft of suitable dimensions between the homodimer surfaces that could accommodate the co-factor-substrate. An appropriate cavity and a hydrophobic nodule were also present in TKTL1/2, into which the diphosphate group fitted, and that was implicated in aminopyrimidine and thiazole ring binding in TKT, respectively. The presence of several identical residues at structurally equivalent positions in TKTL1/2 and TKT identified a network of interactions between the protein and co-factor-substrate, suggesting the functional fidelity of TKTL1/2 as transketolases.ConclusionsOur data support the hypothesis that TKTL1 and TKTL2 are functional transketolases and represent novel therapeutic targets for diabetes and cancer.
Ribose-phosphate pyrophosphokinase (EC 2.7.6.1) is an enzyme that catalyzes the ATP-dependent conversion of ribose-5-phosphate to phosphoribosyl pyrophosphate. The reaction product is a key precursor for the biosynthesis of purine and pyrimidine nucleotides. We report the 2.2 Å crystal structure of the E. coli ribose-phosphate pyrophosphobinase (EcKPRS). The protein has two type I phosphoribosyltransferase folds, related by 2-fold pseudosymmetry. The propeller-shaped homohexameric structure of KPRS is composed of a trimer of dimers, with the C-terminal domains forming the dimeric blades of the propeller and the N-terminal domains forming the hexameric core. The key, conserved active site residues are well-defined in the structure and positioned appropriately to bind substrates, adenosine monophosphate and ribose-5-phosphate. The allosteric site is also relatively well conserved but, in the EcKPRS structure, several residues from a flexible loop occupy the site where the allosteric modulator, adenosine diphosphate, is predicted to bind. The presence of the loop in the allosteric site may be an additional level of regulation, whereby low affinity molecules are precluded from binding. Overall, this study details key structural features of an enzyme that catalyzes a critical step in nucleotide metabolism. This work provides a framework for future studies of this important protein and, as nucleotides are critical for viability, may serve as a foundation for the development of novel anti-bacterial drugs.
BACKGROUND:Analysis of preferred binding regions of a ligand on a protein is important for detecting cryptic binding pockets and improving the ligand selectivity.RESULT:The enhanced sampling approach TAMD has been adapted to allow a ligand to unbind from its native binding site and explore the protein surface. This so-called re-TAMD procedure was then used to explore the interaction between the N terminal peptide of histone H3 and the YEATS domain. Depending on the length of the peptide, several regions of the protein surface were explored. The peptide conformations sampled during the re-TAMD correspond to peptide free diffusion around the protein surface.CONCLUSIONS:The re-TAMD approach permitted to get information on the relative influence of different regions of the N terminal peptide of H3 on the interaction between H3 and YEATS.
P-glycoprotein (ABCB1) is a ATP-binding cassette transporter that plays an important role in the removal of drugs and xenobiotic compounds from the cell. It is also associated with multi-drug resistance in cancer. Here we report novel features of the cryo-EM-derived structure of P-glycoprotein in the post-hydrolytic state: The cytosolic nucleotide-binding domains (NBDs) are separated despite ADP remaining bound to the NBDs. Gaps in the TMDs that connect to the inner hydrophilic cavity are back-filled by detergent head-groups from the annular detergent micelle and are close to two regions predicted to delineate two pseudo-symmetry-related drug-binding sites. In this conformation, the (newly-resolved) N-terminal extension, NBD-TMD linker region and gap-filling detergents all appear to impede NBD dimerisation. We propose a model for the mechanism of action of the exporter where ATP will be bound to the protein for most of the time, consistent with the high physiological ATP concentrations in vivo .
Background: Aquaporins (AQPs) facilitate the passage of small neutral polar molecules across membranes of the cell. In animals there are four distinct AQP subfamilies, whereof AQP8 homologues constitute one of the smallest subfamilies with just one member in man. AQP8 conducts water, ammonia, urea, glycerol and H2O2 through various membranes of animal cells. This passive channel has been connected to a number of phenomena, such as volume change of mitochondria, ammonia neurotoxicity, and mitochondrial dysfunction related to oxidative stress. Currently, there is no experimentally determined structure of an AQP8, hence the structural understanding of this subfamily is limited. The recently solved structure of the plant AQP, AtTIP2;1, which has structural and functional features in common with AQP8s, has opened up for construction of homology models that are likely to be more accurate than previous models. Results: Here we present homology models of seven vertebrate AQP8s. Modeling based on the AtTIP2;1 structure alone resulted in reasonable models except for the pore being blocked by a phenylalanine that is not present in AtTIP2;1. To achieve an open pore, these models were supplemented with models based on the bacterial water specific AQP, EcAqpZ, creating a chimeric monomeric model for each AQP8 isoform. The selectivity filter (also named the aromatic/arginine region), which defines the permeant substrate profile, comprises five amino acid residues in AtTIP2;1, including a histidine coming from loop C. Compared to AtTIP2;1, the selectivity filters of modelled AQP8s only deviates in that they are slightly more narrow and more hydrophobic due to a phenylalanine replacing the histidine from loop C. Interestingly, the models do not exclude the existence of a side pore beneath loop C similar to that described in the structure of AtTIP2;1. Conclusions: Our models concur that AQP8s are likely to have an AtTIP2;1-like selectivity filter. The detailed description of the expected configuration of residues in the selectivity filters of AQP8s provides an excellent starting point for planning of as well as rationalizing the outcome of mutational studies. Our strategy to compile hybrid models based on several templates may prove useful also for other AQPs for which structural information is limited.
Background: The Spumaretrovirinae (foamy viruses) and the Orthoretrovirinae (e.g. HIV) share many similarities both in genome structure and the sequences of the core viral encoded proteins, such as the aspartyl protease and reverse transcriptase. Similarity in the gag region of the genome is less obvious at the sequence level but has been illuminated by the recent solution of the foamy virus capsid (CA) structure. This revealed a clear structural similarity to the orthoretrovirus capsids but with marked differences that left uncertainty in the relationship between the two domains that comprise the structure.Methods: We have applied protein structure comparison methods in order to try and resolve this ambiguous relationship. These included both the DALI method and the SAP method, with rigorous statistical tests applied to the results of both methods. For this, we employed collections of artificial fold 'decoys' (generated from the pair of native structures being compared) to provide a customised background distribution for each comparison, thus allowing significance levels to be estimated.Results: We have shown that the relationship of the two domains conforms to a simple linear correspondence rather than a domain transposition. These similarities suggest that the origin of both viral capsids was a common ancestor with a double domain structure. In addition, we show that there is also a significant structural similarity between the amino and carboxy domains in both the foamy and ortho viruses.Conclusions: These results indicate that, as well as the duplication of the double domain capsid, there may have been an even more ancient gene-duplication that preceded the double domain structure. In addition, our structure comparison methodology demonstrates a general approach to problems where the components have a high intrinsic level of similarity.
After publication of the article [1], it has been brought to our attention that there is a discrepancy between the publication date on the pdf and online formats. The date on the pdf is 6th February 2018 and online is 5th February 2018. The correct publication date is the one on the pdf, 6th February 2018.
Background: Syphilis continues to be a major global health threat with 11 million new infections each year, and a global burden of 36 million cases. The causative agent of syphilis, Treponema pallidum subspecies pallidum, is a highly virulent bacterium, however the molecular mechanisms underlying T. pallidum pathogenesis remain to be definitively identified. This is due to the fact that T. pallidum is currently uncultivatable, inherently fragile and thus difficult to work with, and phylogenetically distinct with no conventional virulence factor homologs found in other pathogens. In fact, approximately 30% of its predicted protein-coding genes have no known orthologs or assigned functions. Here we employed a structural bioinformatics approach using Phyre2-based tertiary structure modeling to improve our understanding of T. pallidum protein function on a proteome-wide scale. Results: Phyre2-based tertiary structure modeling generated high-confidence predictions for 80% of the T. pallidum proteome (780/978 predicted proteins). Tertiary structure modeling also inferred the same function as primary structure-based annotations from genome sequencing pipelines for 525/605 proteins (87%), which represents 54% (525/978) of all T. pallidum proteins. Of the 175 T. pallidum proteins modeled with high confidence that were not assigned functions in the previously annotated published proteome, 167 (95%) were able to be assigned predicted functions. Twenty-one of the 175 hypothetical proteins modeled with high confidence were also predicted to exhibit significant structural similarity with proteins experimentally confirmed to be required for virulence in other pathogens. Conclusions: Phyre2-based structural modeling is a powerful bioinformatics tool that has provided insight into the potential structure and function of the majority of T. pallidum proteins and helped validate the primary structure-based annotation of more than 50% of all T. pallidum proteins with high confidence. This work represents the first T. pallidum proteome-wide structural modeling study and is one of few studies to apply this approach for the functional annotation of a whole proteome.
Background Due to dengue virus disease, half of the world population is at severe health risk. Viral encoded NS2B-NS3 protease complex causes cleavage in the nonstructural region of the viral polyprotein. The cleavage is essentially required for fully functional viral protein. It has already been reported that if function of NS2B-NS3 complex is disrupted, viral replication is inhibited. Therefore, the NS2B-NS3 is a well-characterized target for designing antiviral drug. Results In this study docking analysis was performed with active site of dengue NS2B-NS3 protein with selected plant flavonoids. More than 100 flavonoids were used for docking analysis. On the basis of docking results 10 flavonoids might be considered as the best inhibitors of NS2B-NS3 protein. The interaction studies showed resilient interactions between ligand and receptor atoms. Furthermore, QSAR and SAR studies were conducted on the basis of NS2B-NS3 protease complex docking results. The value of correlation coefficient ( r ) 0.95 shows that there was a good correlation between flavonoid structures and selected properties. Conclusion We hereby suggest that plant flavonoids could be used as potent inhibitors of dengue NS2B-NS3 protein and can be used as antiviral agents against dengue virus. Out of more than hundred plant flavonoids, ten flavonoid structures are presented in this study. On the basis of best docking results, QSAR and SAR studies were performed. These flavonoids can directly work as anti-dengue drug or with little modifications in their structures.
BACKGROUND:Dramatic progress has recently been made in cryo-electron microscopy technologies, which now make possible the reconstruction of a growing number of biomolecular structures to near-atomic resolution. However, the need persists for fitting and refinement approaches that address those cases that require modeling assistance.METHODS:In this paper, we describe algorithms to optimize the performance of such medium-resolution refinement methods. These algorithms aim to automatically optimize the parameters that define the density shape of the flexibly fitted model, as well as the time-dependent damper cutoff distance. Atomic distance constraints can be prescribed for cases where extra containment of parts of the structure is helpful, such as in regions where the density map is poorly defined. Also, we propose a simple stopping criterion that estimates the probable onset of overfitting during the simulation.RESULTS:The new set of algorithms produce more accurate fitting and refinement results, and yield a faster rate of convergence of the trajectory toward the fitted conformation. The latter is also more reliable due to the overfitting warning provided to the user.CONCLUSIONS:The algorithms described here were implemented in the new Damped-Dynamics Flexible Fitting simulation tool "DDforge" in the Situs package.
Background Eukaryotic translation initiation factor 1A (eIF1A) is universally conserved in all organisms. It has multiple functions in translation initiation, including assembly of the ribosomal pre-initiation complexes, mRNA binding, scanning, and ribosomal subunit joining. eIF1A binds directly to the small ribosomal subunit, as well as to several other translation initiation factors. The structure of an eIF1A homolog, the eIF1A domain-containing protein (eIF1AD) was recently determined but its biological functions are unknown. Since eIF1AD has a known structure, as well as a homolog, whose structure and functions have been extensively studied, it is a very attractive target for sequence and structure analysis. Results Structure/sequence analysis of eIF1AD found significant conservation in the surfaces corresponding to the ribosome-binding surfaces of its paralog eIF1A, including a nearly invariant surface-exposed tryptophan residue, which plays an important role in the interaction of eIF1A with the ribosome. These results indicate that eIF1AD may bind to the ribosome, similar to its paralog eIF1A, and could have roles in ribosome biogenenesis or regulation of translation. We identified conserved surfaces and sequence motifs in the folded domain as well as the C-terminal tail of eIF1AD, which are likely protein-protein interaction sites. The roles of these regions for eIF1AD function remain to be determined. We have also identified a set of trypanosomatid-specific surface determinants in eIF1A that could be a promising target for development of treatments against these parasites. Conclusions The results described here identify regions in eIF1A and eIF1AD that are likely to play major functional roles and are promising therapeutic targets. Our findings and hypotheses will promote new research and help elucidate the functions of eIF1AD.
Background T7 lysozyme (T7L), also known as N-acetylmuramoyl-L-alanine amidase, is a T7 bacteriophage gene product. It involves two functions: It can cut amide bonds in the bacterial cell wall and interacts with T7 RNA polymerase (T7RNAP) as a part of transcription inhibition. In this study, with the help of molecular dynamics (MD) calculations and computational interaction studies, we investigated the effect of varying pH conditions on conformational flexibilities of T7L and their influence on T7RNAP -T7L interactions. Results From the MD studies of the T7L at three different pH strengths viz. 5, neutral and 7.9 it was observed that T7L structure at pH 5 exhibited less stable nature with more residue level fluctuations, decrease of secondary structural elements and less compactness as compared to its counterparts: neutral pH and pH 7.9. The T-pad analysis of the MD trajectories identified local fluctuations in few residues that influenced the conformational differences in three pH strengths. From the docking of the minimum energy representative structures of T7L at different pH strengths (obtained from the free energy landscape analysis) with T7RNAP structures at same pH strengths, we saw strong interaction patterns at pH 7.9 and pH 5. The MD analysis of these complexes also confirmed the observations of docking study. From the combined in silico studies, it was observed that there are conformational changes in N-terminal and near helix 1 of T7L at different pH strengths, which are involved in the T7RNAP interaction, thereby varying the interaction pattern. Conclusion Since T7L has been used for developing novel therapeutics and T7RNAP one of the most biologically useful protein in both in-vitro and in vivo experiments, this in silico study of pH dependent conformational differences in T7L and the differential interaction with T7RNAP at different pH can provide a significant insight into the structural investigations on T7L and T7RNAP in varying pH environments.
BACKGROUND:Amyloidogenic proteins are most often associated with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, and Huntington's disease, but there are more than two dozen human proteins known to form amyloid fibrils associated with disease. Lysozyme is an antimicrobial protein that is used as a general model to study amyloid fibril formation. Studies aimed at elucidating the process of amyloid formation of lysozyme tend to focus on partial unfolding of the native state due to the relative instability of mutant amyloidogenic variants. While this is well supported, the data presented here suggest the native structure of the variants may also play a role in primary nucleation.RESULTS:Three-dimensional structural analysis identified lysozyme residues 21, 62, 104, and 122 as displaced in both amyloidogenic variants compared to wild type lysozyme. Residue interaction network (RIN) analysis found greater clustering of residues 112-117 in amyloidogenic variants of lysozyme compared to wild type. An analysis of the most energetically favored predicted dimers and trimers provided further evidence for a role for residues 21, 62, 104, 122, and 112-117 in amyloid formation.CONCLUSIONS:This study used lysozyme as a model to demonstrate the utility of combining 3D structural analysis with RIN analysis for studying the general process of amyloidogenesis. Results indicated that binding of two or more amyloidogenic lysozyme mutants may be involved in amyloid nucleation by placing key residues (21, 62, 104, 122, and 112-117) in proximity before partial unfolding occurs. Identifying residues in the native state that may be involved in amyloid formation could provide novel drug targets to prevent a range of amyloidoses.
Background: Myelin is a multilayered proteolipid sheath wrapped around selected axons in the nervous system. Its constituent proteins play major roles in forming of the highly regular membrane structure. P2 is a myelin-specific protein of the fatty acid binding protein (FABP) superfamily, which is able to stack lipid bilayers together, and it is a target for mutations in the human inherited neuropathy Charcot-Marie-Tooth disease. A conserved residue that has been proposed to participate in membrane and fatty acid binding and conformational changes in FABPs is Phe57. This residue is thought to be a gatekeeper for the opening of the portal region upon ligand entry and egress. Results: We performed a structural characterization of the F57A mutant of human P2. The mutant protein was crystallized in three crystal forms, all of which showed changes in the portal region and helix α2. In addition, the behaviour of the mutant protein upon lipid bilayer binding suggested more unfolding than previously observed for wild-type P2. On the other hand, membrane binding rendered F57A heat-stable, similarly to wild-type P2. Atomistic molecular dynamics simulations showed opening of the side of the discontinuous β barrel, giving important indications on the mechanism of portal region opening and ligand entry into FABPs. The results suggest a central role for Phe57 in regulating the opening of the portal region in human P2 and other FABPs, and the F57A mutation disturbs dynamic cross-correlation networks in the portal region of P2. Conclusions: Overall, the F57A variant presents similar properties to the P2 patient mutations recently linked to Charcot-Marie-Tooth disease. Our results identify Phe57 as a residue regulating conformational changes that may accompany membrane surface binding and ligand exchange in P2 and other FABPs.
Background So far, little is known about the molecular mechanisms of amyotrophic lateral sclerosis onset and progression caused by SOD1 mutations. One of the hypotheses is based on SOD1 misfolding resulting from mutations and subsequent deposition of its cytotoxic aggregates. This hypothesis is complicated by the fact that known SOD1 mutations of similar clinical effect could be distributed over the whole protein structure. Results In this work, a measure of hydrogen bond stability in conformational states was studied with elastic network analysis of 35 SOD1 mutants. Twenty-eight hydrogen bonds were detected in nine of 35 mutants with their stability being significantly different from that with the wild-type. These hydrogen bonds were formed by the amino acid residues known from the literature to be located in contact between SOD1 aggregates. Additionally, residues disposed between copper binding sites of both protein subunits were found from the models to form a stiff core, which can be involved in mechanical impulse transduction between these active centres. Conclusions The modelling highlights that both stability of the copper binding site and stability of the dimer can play an important role in ALS progression.
