A three-dimensional model of the U1 small nuclear ribonucleoprotein particle

Most of the pre-mRNAs in the eukaryotic cell are comprised of protein-coding exons and non-protein-coding introns. The introns are removed and the exons are ligated together, or spliced, by a large, macromolecular complex known as the spliceosome. This RNA-protein assembly is made up of five uridine-rich small nuclear RNAs (U1-, U2-, U4-, U5- and U6-snRNA) as well over 300 proteins, which form small nuclear ribonucleoprotein particles (snRNPs). Initial recognition of the 5′ exon/intron splice site is mediated by the U1 snRNP, which is composed of the U1 snRNA as well as at least ten proteins. By combining structural informatics tools with the available biochemical and crystallographic data, we attempted to simulate a complete, three dimensional U1 snRNP from the silk moth, Bombyx mori. Comparison of our model with empirically derived crystal structures and electron micrographs pinpoints both the strengths and weaknesses in the in silico determination of macromolecular complexes. One of the most striking differences between our model and experimentally generated structures is in the positioning of the U1 snRNA stem-loops. This highlights the continuing difficulties in generating reliable, complex RNA structures; however, three-dimensional modeling of individual protein subunits by threading provided models of biological significance and the use of both automated and manual docking strategies generated a complex that closely reflects the assembly found in nature. Yet, without utilizing experimentally-derived contacts to select the most likely docking scenario, ab initio docking would fall short of providing a reliable model. Our work shows that the combination of experimental data with structural informatics tools can result in generation of near-native macromolecular complexes.     

Structural and posttranslational analysis of human calcium-binding protein,spermatid-associated 1

Recently, we detected a novel biomarker in human saliva called calcium-binding protein, spermatid-associated 1 (CABS1). CABS1 protein had previously been described only in testis, and little was known of its characteristics other than it was considered a structurally disordered protein. Levels of human CABS1 (hCABS1) in saliva correlate with stress, whereas smaller sized forms of hCABS1 in saliva are associated with resilience to stress. Interestingly, hCABS1 also has an anti-inflammatory peptide sequence near its carboxyl terminus, similar to that of a rat prohormone, submandibular rat 1. We performed phylogenetic and sequence analysis of hCABS1. We found that from 72 CABS1 sequences currently annotated in the National Center for Biotechnology Information protein database, only 14 contain the anti-inflammatory domain “TxIFELL,” all of which are primates. We performed structural unfoldability analysis using PONDER and FoldIndex and discovered three domains that are highly disordered. Predictions of three-dimensional structure of hCABS1 using RaptorX, IonCom, and I-TASSER software agreed with these findings. Predicted neutrophil elastase cleavage density also correlated with hCABS1 regions of high structural disorder. Ligand binding prediction identified Ca²⁺, Mg²⁺, Zn²⁺, leucine, and thiamine pyrophosphate, a pattern observed in enzymes associated with energy metabolism and mitochondrial localization. These new observations on hCABS1 raise intriguing questions about the interconnection between the autonomic nervous system, stress, and the immune system. However, the precise molecular mechanisms involved in the complex biology of hCABS1 remain unclear. We provide a detailed in silico analysis of relevant aspects of the structure and function of hCABS1 and postulate extracellular and intracellular roles.     

Identification of the Catalytic Site of Phospholipase D2 (PLD2) Newly Described Guanine Nucleotide Exchange Factor Activity

We have demonstrated that phospholipase D2 (PLD2) is a guanine nucleotide exchange factor (GEF) for Rac2 and determined the PLD2 domains and amino acid site(s) responsible for its GEF activity. Experiments using GST fusion proteins or GST-free counterparts, purified proteins revealed that the PX domain is sufficient to exert GEF activity similar to full-length PLD2. The PLD2-GEF catalytic site is formed by a hydrophobic pocket of residues Phe-107, Phe-129, Leu-166, and Leu-173, all of which are in the PX domain. A nearby Arg-172 is also important in the overall activity. PX mutants altering any of those five amino acids fail to have GEF activity but still bind to Rac2, while their lipase activity was mostly unaffected. In addition to the PX domain, a region in the pleckstrin homology domain (Ile-306–Ala-310) aids in the PX-mediated GEF activity by providing a docking site to hold Rac2 in place during catalysis. We conclude that PLD2 is a unique GEF, with the PX being the major catalytic domain for its GEF activity, whereas the pleckstrin homology domain assists in the PX-mediated activity. The physiological relevance of this novel GEF in cell biology is demonstrated here in chemotaxis and phagocytosis of leukocytes, as the specific PX and PH mutants abolished cell function. Thus, this study reveals for the first time the catalytic site that forms the basis for the mechanism behind the GEF activity of PLD2.     

Evidence of colorectal cancer risk associated variant Lys25Ser in the proximity of human bone morphogenetic protein 2

Colorectal cancer (CRC) is the third most prevalent cancer and fourth leading cause of cancer-related deaths globally. It has been shown that the nsSNP variants play an important role in diseases, however it remained unclear how these variants are associated with the disease. Recently, several CRC risk associated SNPs have been discovered, however rs961253 (Lys25Arg at 20p12.3) located in the proximity of bone morphogenetic protein 2 (Bmp2) and fermitin family homolog 1 Fermt1 genes have been reported to be highly associated with the CRC risk. Here we provide evidence for the first time in silico biological functional and structural implications of non-synonymous (nsSNPs) CRC disease-associated variant Lys25Arg via molecular dynamic (MD) simulation. Protein structural analysis was performed with a particular variant allele (A/C, Lys25Arg) and compared with the predicted native protein structure. Our results showed that this nsSNP will cause changes in the protein structure and as a result is associated with the disease. In addition to the native and mutant 3D structures of CRC associated risk allele protein domain (CRAPD), they were also analyzed using solvent accessibility models for further protein stability confirmation. Taken together, this study confirmed that this variant has functional effect and structural impact on the CRAPD and may play an important role in CRC disease progression; hence it could be a reasonable approach for studying the effect of other deleterious variants in future studies.     

Antiviral potential of 4-hydroxypanduratin A,secondary metabolite of Fingerroot,Boesenbergia pandurata (Schult.), towards Japanese Encephalitis virus NS2B/NS3 protease

4-hydroxypanduratin A is a secondary metabolite of Boesenbergia pandurata Schult. (Fingerroot) plant with various pharmacological activities such as neuroprotective, potent antioxidant, antibacterial and antifungal. Flaviviral NS2B/NS3 protease activity is essential for polyprotein processing and viral replication for Japanese Encephalitis Virus (JEV), a major cause of Acute Encephaltis in Asia. Inhibition of formation of this complex by arresting the binding of NS2B with NS3 would reduce the enzyme''s activity to meager proportions and hence would prevent further viral proliferation. The automated 3D structure of NS2B protein of the JEV GP78 was predicted based on the sequence-to-structure-to-function paradigm using I-TASSER and the function of NS2B protein was inferred by matching to other known proteins. The stereochemical quality of predicted structure was checked by PROCHECK. The antiviral activity of 4-hydroxypanduratin A against NS2B protein as a potential drug has been elucidated in this paper. Docking simulation analysis showed 4-hydroxypanduratin A as potential inhibitor of NS2B protein/cofactor which is necessary for NS3 protease activity. 220 derivatives of 4-hydroxypanduratin A were virtually screened with rigid criteria of Lipinski''s rule of 5 using Autodock4.2. 4-hydroxypanduratin A was found interacting with target hydrophilic domain in NS2B protein by two Hbonds (Gly80 and Asp81) with active residues, several hydrophobic interactions, Log P value of 5.6, inhibition constant (Ki) of 51.07nM and lowest binding energy of -9.95Kcal/Mol. Hence, 4-hydroxypanduratin A targeted to Site 2 will have sufficient profound effect to inhibit protease activity to abrogate viral replication. It could be a promising potential drug candidate for JEV infections using NS2B Site 2 as a Drug target.     

An in silico approach to understand the structure–function properties of a serine protease (Bacifrinase) from Bacillus cereus and experimental evidence to support the interaction of Bacifrinase with fibrinogen and thrombin

Microbial fibrinogenolytic serine proteases find therapeutic applications in the treatment of thrombosis- and hyperfibrinogenemia-associated disorders. However, analysis of structure–function properties of an enzyme is utmost important before its commercial application. In this study, an attempt has been made to understand the structure of a fibrinogenolytic protease enzyme, “Bacifrinase” from Bacillus cereus strain AB01. From the molecular dynamics trajectory analysis, the modelled three-dimensional structure of the protease was found to be stable and the presence of a catalytic triad made up of Asp102, His83 and Ser195 suggests that it is a serine protease. To understand the mechanism of enzyme–substrate and enzyme–inhibitor interactions, the equilibrated protein was docked with human fibrinogen (the physiological substrate of this enzyme), human thrombin and with ten selective protease inhibitors. The Bacifrinase–chymostatin interaction was the strongest among the selected protease inhibitors. The serine protease inhibitor phenyl methane sulphonyl fluoride was found to interact with the Ser134 residue of Bacifrinase. Furthermore, protein–protein docking study revealed the fibrinogenolytic property of Bacifrinase and its interaction with Aα-, Bβ- and Cγ-chains human fibrinogen to a different extent. However, biochemical analysis showed that Bacifrinase did not hydrolyse the γ-chain of fibrinogen. The in silico and spectrofluorometric studies also showed interaction of Bacifrinase with thrombin as well as fibrinogen with a Kd value of 16.5 and .81 nM, respectively. Our findings have shed light on the salient structural features of Bacifrinase and confirm that it is a fibrinogenolytic serine protease.     

Insights from the docking and molecular dynamics simulation of the Phosphopantetheinyl transferase (PptT) structural model from Mycobacterium tuberculosis

A great challenge is posed to the treatment of tuberculosis due to the evolution of multidrug-resistant (MDR) and extensively drugresistant (XDR) strains of Mycobacterium tuberculosis in recent times. The complex cell envelope of the bacterium contains unusual structures of lipids which protects the bacterium from host enzymes and escape immune response. To overcome the drug resistance, targeting “drug targets” which have a critical role in growth and virulence factor is a novel approach for better tuberculosis treatment. The enzyme Phosphopantetheinyl transferase (PptT) is an attractive drug target as it is primarily involved in post translational modification of various types-I polyketide synthases and assembly of mycobactin, which is required for lipid virulence factors. Our in silico studies reported that the structural model of M.tuberculosis PptT characterizes the structure-function activity. The refinement of the model was carried out with molecular dynamics simulations and was analyzed with root mean square deviation (RMSD), and radius of gyration (Rg). This confirmed the structural behavior of PptT in dynamic system. Molecular docking with substrate coenzyme A (CoA) identified the binding pocket and key residues His93, Asp114 and Arg169 involved in PptT-CoA binding. In conclusion, our results show that the M.tuberculosis PptT model and critical CoA binding pocket initiate the inhibitor design of PptT towards tuberculosis treatment.     

A protocol for computer-based protein structure and function prediction

Genome sequencing projects have ciphered millions of protein sequence, which require knowledge of their structure and function to improve the understanding of their biological role. Although experimental methods can provide detailed information for a small fraction of these proteins, computational modeling is needed for the majority of protein molecules which are experimentally uncharacterized. The I-TASSER server is an on-line workbench for high-resolution modeling of protein structure and function. Given a protein sequence, a typical output from the I-TASSER server includes secondary structure prediction, predicted solvent accessibility of each residue, homologous template proteins detected by threading and structure alignments, up to five full-length tertiary structural models, and structure-based functional annotations for enzyme classification, Gene Ontology terms and protein-ligand binding sites. All the predictions are tagged with a confidence score which tells how accurate the predictions are without knowing the experimental data. To facilitate the special requests of end users, the server provides channels to accept user-specified inter-residue distance and contact maps to interactively change the I-TASSER modeling; it also allows users to specify any proteins as template, or to exclude any template proteins during the structure assembly simulations. The structural information could be collected by the users based on experimental evidences or biological insights with the purpose of improving the quality of I-TASSER predictions. The server was evaluated as the best programs for protein structure and function predictions in the recent community-wide CASP experiments. There are currently >20,000 registered scientists from over 100 countries who are using the on-line I-TASSER server.     

In silico structural, functional and pathogenicity evaluation of a novel mutation: an overview of HSD3B2 gene mutations

Mutations of 3 beta hydroxysteroid dehydrogenase type II (HSD3B2) gene result in different clinical consequences. We explain a patient who demonstrated a salt wasting form of 3βHSD deficiency in infancy. Signs of hyponatremia and hyperkalemia were recognized in the infant with ambiguous genitalia and perineal hypospadias. The 46,XY male was genotyped by direct sequencing of HSD3B2 gene. Steroid profiles showed elevated concentration of 17 hydroxyprogesterone, and decrease in concentration of cortisol, and testosterone. Dehydroepiandrotone (DHEA) to androstenedione ratio had 6 fold increases. Direct sequencing of the patient revealed homozygous missense A82P mutation in exon 3. This mutation was confirmed by segregation analysis of the parents. Bioinformatic tools were used for in silico structural and functional analyses. Also, the pathological effect of the mutation was validated by different software. Alanine is a conserved amino acid in the membrane binding domain of the enzyme and proline substitution was predicted to destabilize the protein. This report may highlight the importance of the screening programs of the disorder in Iran.