MODBASE, a database of annotated comparative protein structure models, and associated resources

psnGPCRdb: The Structure-network Database of G Protein Coupled Receptors

G protein coupled receptors (GPCRs) are critical eukaryotic signal transduction gatekeepers and represent the largest protein superfamily in the human proteome, with more than 800 members. They share seven transmembrane helices organized in an up-down bundle architecture. GPCR-mediated signaling pathways have been linked to numerous human diseases, and GPCRs are the targets of approximately 35% of all drugs currently on the market. Structure network analysis, a graph theory-based approach, represents a cutting-edge tool to deeply understand GPCR function, which strongly relies on communication between the extracellular and intracellular poles of their structure. psnGPCRdb stores the structure networks (i.e., linked nodes, hubs, communities and communication pathways) computed on all updated GPCR structures in the Protein Data Bank, in their isolated states or in complex with extracellular and/or intracellular molecules. The structure communication signatures of a sub-family or family of GPCRs as well as of their small-molecule activators or inhibitors are stored as consensus networks. The database stores also all meaningful structure network-based comparisons (i.e., difference networks) of functionally different states (i.e., inactive or active) of a given receptor sub-type, or of consensus networks representative of a receptor sub-type, type, sub-family or family. Single or consensus GPCR networks hold also information on amino acid conservation. The database allows to graphically analyze 3D structure networks together with interactive data-tables. Ligand-centric networks can be analyzed as well. psnGPCRdb is unique and represents a powerful resource to unravel GPCR function with important implications in cell signaling and drug design. psnGPCRdb is freely available at: http://webpsn.hpc.unimo.it/psngpcr.php.

Spectrin: an alternate target for cytoskeletal drugs

Cytoskeletal drugs having enormous therapeutic potential act on the cytoskeletal components like actin, tubulin either by promoting polymerization or destabilizing the same. Here we present the interaction of the popular cytoskeletal drugs such as taxol, latrunculin and cytochalasin with spectrin, a huge protein with multi domains that forms the cytoskeletal network. Particularly, the actin binding domain of spectrin regulates the dynamics of the actin cytoskeleton. We followed the binding of these drugs to its actin binding domain and intact spectrin as well. These drugs bind with moderate affinity (K_b ∼ 10⁴ M^-1) and the interaction with actin binding domain is entropy driven and hydrophobic in nature as determined by Van't Hoff plot. The docking studies and molecular dynamics simulations further corroborate the experimental findings. Particularly the higher binding constants in the case of latrunculin and cytochalasin to the actin binding domain of spectrin suggest the binding sites are presumably located in its actin binding domain.Communicated by Ramaswamy H. Sarma.

Graphene BioFET sensors for SARS-CoV-2 detection: a multiscale simulation approach

Biological Field-Effect Transistors (BioFETs) have already demonstrated enormous potential for detecting minute amounts of ions and molecules. The use of two-dimensional (2D) materials has been shown to boost their performance and to enable the design of new applications. This combination deserves special interest in the current pandemic caused by the SARS-CoV-2 virus which demands fast, reliable and cheap detection methods. However, in spite of the experimental advances, there is a lack of a comprehensive and in-depth computational approach to capture the mechanisms underlying the sensor behaviour. Here, we present a multiscale platform that combines detailed atomic models of the molecules with mesoscopic device-level simulations. The fine-level description exploited in this approach accounts for the charge distribution of the receptor, its reconfiguration when the target binds to it, and the consequences in terms of sensitivity on the transduction mechanism. The results encourage the further exploration of improved sensor designs and 2D materials combined with diverse receptors selected to achieve the desired specificity.

Molecular basis of crosstalk in nuclear receptors: heterodimerization between PXR and CAR and the implication in gene regulation

The 48 human nuclear receptors (NRs) form a superfamily of transcription factors that regulate major physiological and pathological processes. Emerging evidence suggests that NR crosstalk can fundamentally change our understanding of NR biology, but detailed molecular mechanisms of crosstalk are lacking. Here, we report the molecular basis of crosstalk between the pregnane X receptor (PXR) and constitutive androstane receptor (CAR), where they form a novel heterodimer, resulting in their mutual inhibition. PXR and CAR regulate drug metabolism and energy metabolism. Although they have been broadly perceived as functionally redundant, a growing number of reports suggests a mutual inhibitory relation, but their precise mode of coordinated action remains unknown. Using methods including RNA sequencing, small-angle X-ray scattering and crosslinking mass spectrometry we demonstrate that the mutual inhibition altered gene expression globally and is attributed to the novel PXR–CAR heterodimerization via the same interface used by each receptor to heterodimerize with its functional partner, retinoid X receptor (RXR). These findings establish an unexpected functional relation between PXR, CAR and RXR, change the perceived functional relation between PXR and CAR, open new perspectives on elucidating their role and designing approaches to regulate them, and highlight the importance to comprehensively investigate nuclear receptor crosstalk.

Tryptophan Operon Diversity Reveals Evolutionary Trends among Geographically Disparate Chlamydia trachomatis Ocular and Urogenital Strains Affecting Tryptophan Repressor and Synthase Function

The obligate intracellular pathogen Chlamydia trachomatis (Ct) is the leading cause of bacterial sexually transmitted infections and blindness globally. To date, Ct urogenital strains are considered tryptophan prototrophs, utilizing indole for tryptophan synthesis within a closed-conformation tetramer comprised of two α (TrpA)- and two β (TrpB)-subunits. In contrast, ocular strains are auxotrophs due to mutations in TrpA, relying on host tryptophan pools for survival. It has been speculated that there is strong selective pressure for urogenital strains to maintain a functional operon. Here, we performed genetic, phylogenetic, and novel functional modeling analyses of 595 geographically diverse Ct ocular, urethral, vaginal, and rectal strains with complete operon sequences. We found that ocular and urogenital, but not lymphogranuloma venereum, TrpA-coding sequences were under positive selection. However, vaginal and urethral strains exhibited greater nucleotide diversity and a higher ratio of nonsynonymous to synonymous substitutions [Pi(a)/Pi(s)] than ocular strains, suggesting a more rapid evolution of beneficial mutations. We also identified nonsynonymous amino acid changes for an ocular isolate with a urogenital backbone in the intergenic region between TrpR and TrpB at the exact binding site for YtgR-the only known iron-dependent transcription factor in Chlamydia-indicating that selective pressure has disabled the response to fluctuating iron levels. In silico effects on protein stability, ligand-binding affinity, and tryptophan repressor (TrpR) affinity for single-stranded DNA (ssDNA) measured by calculating free energy changes (ΔΔG) between Ct reference and mutant tryptophan operon proteins were also analyzed. We found that tryptophan synthase function was likely suboptimal compared to other bacterial tryptophan prototrophs and that a diversity of urogenital strain mutations rendered the synthase nonfunctional or inefficient. The novel mutations identified here affected active sites in an orthosteric manner but also hindered α- and β-subunit allosteric interactions from distant sites, reducing efficiency of the tryptophan synthase. Importantly, strains with mutant proteins were inclined toward energy conservation by exhibiting an altered affinity for their respective ligands compared to reference strains, indicating greater fitness. This is not surprising as l-tryptophan is one of the most energetically costly amino acids to synthesize. Mutations in the tryptophan repressor gene (trpR) among urogenital strains were similarly detrimental to function. Our findings indicate that urogenital strains are evolving more rapidly than previously recognized with mutations that impact tryptophan operon function in a manner that is energetically beneficial, providing a novel host-pathogen evolutionary mechanism for intracellular survival.IMPORTANCE Chlamydia trachomatis (Ct) is a major global public health concern causing sexually transmitted and ocular infections affecting over 130 million and 260 million people, respectively. Sequelae include infertility, preterm birth, ectopic pregnancy, and blindness. Ct relies on available host tryptophan pools and/or substrates to synthesize tryptophan to survive. Urogenital strains synthesize tryptophan from indole using their intact tryptophan synthase (TS). Ocular strains contain a trpA frameshift mutation that encodes a truncated TrpA with loss of TS function. We found that TS function is likely suboptimal compared to other tryptophan prototrophs and that urogenital stains contain diverse mutations that render TS nonfunctional/inefficient, evolve more rapidly than previously recognized, and impact operon function in a manner that is energetically beneficial, providing an alternative host-pathogen evolutionary mechanism for intracellular survival. Our research has broad scientific appeal since our approach can be applied to other bacteria that may explain evolution/survival in host-pathogen interactions.

Dynamical Rearrangement of Human Epidermal Growth Factor Receptor 2 upon Antibody Binding: Effects on the Dimerization

Human epidermal growth factor 2 (HER2) is a ligand-free tyrosine kinase receptor of the HER family that is overexpressed in some of the most aggressive tumours. Although it is known that HER2 dimerization involves a specific region of its extracellular domain, the so-called “dimerization arm”, the mechanism of dimerization inhibition remains uncertain. However, uncovering how antibody interactions lead to inhibition of HER2 dimerization is of key importance in understanding its role in tumour progression and therapy. Herein, we employed several computational modelling techniques for a molecular-level understanding of the interactions between HER and specific anti-HER2 antibodies, namely an antigen-binding (Fab) fragment (F0178) and a single-chain variable fragment from Trastuzumab (scFv). Specifically, we investigated the effects of antibody-HER2 interactions on the key residues of “dimerization arm” from molecular dynamics (MD) simulations of unbound HER (in a total of 1 µs), as well as ScFv:HER2 and F0178:HER2 complexes (for a total of 2.5 µs). A deep surface analysis of HER receptor revealed that the binding of specific anti-HER2 antibodies induced conformational changes both in the interfacial residues, which was expected, and in the ECDII (extracellular domain), in particular at the “dimerization arm”, which is critical in establishing protein–protein interface (PPI) interactions. Our results support and advance the knowledge on the already described trastuzumab effect on blocking HER2 dimerization through synergistic inhibition and/or steric hindrance. Furthermore, our approach offers a new strategy for fine-tuning target activity through allosteric ligands.

3D spatial organization and network-guided comparison of mutation profiles in Glioblastoma reveals similarities across patients

Glioblastoma multiforme (GBM) is the most aggressive type of brain tumor. Molecular heterogeneity is a hallmark of GBM tumors that is a barrier in developing treatment strategies. In this study, we used the nonsynonymous mutations of GBM tumors deposited in The Cancer Genome Atlas (TCGA) and applied a systems level approach based on biophysical characteristics of mutations and their organization in patient-specific subnetworks to reduce inter-patient heterogeneity and to gain potential clinically relevant insights. Approximately 10% of the mutations are located in "patches" which are defined as the set of residues spatially in close proximity that are mutated across multiple patients. Grouping mutations as 3D patches reduces the heterogeneity across patients. There are multiple patches that are relatively small in oncogenes, whereas there are a small number of very large patches in tumor suppressors. Additionally, different patches in the same protein are often located at different domains that can mediate different functions. We stratified the patients into five groups based on their potentially affected pathways that are revealed from the patient-specific subnetworks. These subnetworks were constructed by integrating mutation profiles of the patients with the interactome data. Network-guided clustering showed significant association between the groups and patient survival (P-value = 0.0408). Also, each group carries a set of signature 3D mutation patches that affect predominant pathways. We integrated drug sensitivity data of GBM cell lines with the mutation patches and the patient groups to analyze the possible therapeutic outcome of these patches. We found that Pazopanib might be effective in Group 3 by targeting CSF1R. Additionally, inhibiting ATM that is a mediator of PTEN phosphorylation may be ineffective in Group 2. We believe that from mutations to networks and eventually to clinical and therapeutic data, this study provides a novel perspective in the network-guided precision medicine.

Clinical Persistence of Chlamydia trachomatis Sexually Transmitted Strains Involves Novel Mutations in the Functional αββα Tetramer of the Tryptophan Synthase Operon

Clinical persistence of Chlamydia trachomatis (Ct) sexually transmitted infections (STIs) is a major public health concern. In vitro persistence is known to develop through interferon gamma (IFN-γ) induction of indoleamine 2,3-dioxygenase (IDO), which catabolizes tryptophan, an essential amino acid for Ct replication. The organism can recover from persistence by synthesizing tryptophan from indole, a substrate for the enzyme tryptophan synthase. The majority of Ct strains, except for reference strain B/TW-5/OT, contain an operon comprised of α and β subunits that encode TrpA and TrpB, respectively, and form a functional αββα tetramer. However, trpA mutations in ocular Ct strains, which are responsible for the blinding eye disease known as trachoma, abrogate tryptophan synthesis from indole. We examined serial urogenital samples from a woman who had recurrent Ct infections over 4 years despite antibiotic treatment. The Ct isolates from each infection episode were genome sequenced and analyzed for phenotypic, structural, and functional characteristics. All isolates contained identical mutations in trpA and developed aberrant bodies within intracellular inclusions, visualized by transmission electron microscopy, even when supplemented with indole following IFN-γ treatment. Each isolate displayed an altered αββα structure, could not synthesize tryptophan from indole, and had significantly lower trpBA expression but higher intracellular tryptophan levels compared with those of reference Ct strain F/IC-Cal3. Our data indicate that emergent mutations in the tryptophan operon, which were previously thought to be restricted only to ocular Ct strains, likely resulted in in vivo persistence in the described patient and represents a novel host-pathogen adaptive strategy for survival.IMPORTANCEChlamydia trachomatis (Ct) is the most common sexually transmitted bacterium with more than 131 million cases occurring annually worldwide. Ct infections are often asymptomatic, persisting for many years despite treatment. In vitro recovery from persistence occurs when indole is utilized by the organism's tryptophan synthase to synthesize tryptophan, an essential amino acid for replication. Ocular but not urogenital Ct strains contain mutations in the synthase that abrogate tryptophan synthesis. Here, we discovered that the genomes of serial isolates from a woman with recurrent, treated Ct STIs over many years were identical with a novel synthase mutation. This likely allowed long-term in vivo persistence where active infection resumed only when tryptophan became available. Our findings indicate an emerging adaptive host-pathogen evolutionary strategy for survival in the urogenital tract that will prompt the field to further explore chlamydial persistence, evaluate the genetics of mutant Ct strains and fitness within the host, and their implications for disease pathogenesis.

Supersecondary Structures and Fragment Libraries

The use of smotifs and fragment libraries has proven useful to both simplify and increase the quality of protein models. Here, we present Profrager, a tool that automatically generates putative structural fragments to reproduce local motifs of proteins given a target sequence. Profrager is highly customizable, allowing the user to select the number of fragments per library, the ranking method is able to generate fragments of all sizes, and it was recently modified to include the possibility of output exclusively smotifs.

A novel approach of virulome based reverse vaccinology for exploring and validating peptide-based vaccine candidates against the most troublesome nosocomial pathogen: Acinetobacter baumannii

Acinetobacter baumannii is one of the major cause of nosocomial infections around the globe. The emergence of hyper-virulent strains of the pathogen greatly narrows down therapeutic options for patients infected with this red alert superbug. Development of a peptide-based vaccine can offers an alternative, attractive, and cost-effective remedy for multidrug-resistant A. baumannii associated complications. Herein, we introduced a novel virulome based Reverse Vaccinology for screening peptide based vaccine candidates against A. baumannii and its validation using a negative control. The pipeline screened "FYLNDQPVS" of polysaccharide export outer membrane protein (EpsA) and "LQNNTRRMK" of chaperone-usher pathway protein B (CsuB) as broad-spectrum peptides for induction of targeted immune responses. The 9-mer epitope of both proteins was rendered virulent, antigenic, non-allergen, and highly conserved among thirty-four completely annotated strains. Interactome examination unravels peptides protein direct and indirect interactions with biological significant pathways, essential for A. baumannii pathogenesis and survival. Protein-peptide docking aids in addition by unveiling deep binding of the epitopes in the active site of the most prevalent binding allele in the human population-the DRB1*0101. Both the proteins till to date are not characterized for immunoprotective efficacy and desirable to be deciphered experimentally. The designed series of in silico filters rejected few recently reported peptide and non-peptide vaccine targets and has delivered outcomes, which we believe will enrich the existing knowledge of vaccinology against this life-threatening human pathogen.

An effector from the Huanglongbing-associated pathogen targets citrus proteases

The citrus industry is facing an unprecedented challenge from Huanglongbing (HLB). All cultivars can be affected by the HLB-associated bacterium 'Candidatus Liberibacter asiaticus' (CLas) and there is no known resistance. Insight into HLB pathogenesis is urgently needed in order to develop effective management strategies. Here, we use Sec-delivered effector 1 (SDE1), which is conserved in all CLas isolates, as a molecular probe to understand CLas virulence. We show that SDE1 directly interacts with citrus papain-like cysteine proteases (PLCPs) and inhibits protease activity. PLCPs are defense-inducible and exhibit increased protein accumulation in CLas-infected trees, suggesting a role in citrus defense responses. We analyzed PLCP activity in field samples, revealing specific members that increase in abundance but remain unchanged in activity during infection. SDE1-expressing transgenic citrus also exhibit reduced PLCP activity. These data demonstrate that SDE1 inhibits citrus PLCPs, which are immune-related proteases that enhance defense responses in plants.

Molecular mechanisms and structural features of cardiomyopathy-causing troponin T mutants in the tropomyosin overlap region

Point mutations in genes encoding sarcomeric proteins are the leading cause of inherited primary cardiomyopathies. Among them are mutations in the TNNT2 gene that encodes cardiac troponin T (TnT). These mutations are clustered in the tropomyosin (Tm) binding region of TnT, TNT1 (residues 80-180). To understand the mechanistic changes caused by pathogenic mutations in the TNT1 region, six hypertrophic cardiomyopathy (HCM) and two dilated cardiomyopathy (DCM) mutants were studied by biochemical approaches. Binding assays in the absence and presence of actin revealed changes in the affinity of some, but not all, TnT mutants for Tm relative to WT TnT. HCM mutants were hypersensitive and DCM mutants were hyposensitive to Ca²⁺ in regulated actomyosin ATPase activities. To gain better insight into the disease mechanism, we modeled the structure of TNT1 and its interactions with Tm. The stability predictions made by the model correlated well with the affinity changes observed in vitro of TnT mutants for Tm. The changes in Ca²⁺ sensitivity showed a strong correlation with the changes in binding affinity. We suggest the primary reason by which these TNNT2 mutations between residues 92 and 144 cause cardiomyopathy is by changing the affinity of TnT for Tm within the TNT1 region.

An in silico strategy for identification of novel drug targets against Plasmodium falciparum

The apicomplexan parasite Plasmodium falciparum is responsible for global malaria burden. With the reported resistance to artemisinin chemotherapy, there is an urgent need to maintain early phase drug discovery and identify novel drug targets for successful eradication of the pathogen from the host. In our previous work on comparative genomics study for identification of putative essential genes and therapeutic candidates in P. falciparum, we predicted 11 proteins as anti-malarial drug targets from PlasmoDB database. In this paper, we made an attempt for identification of novel drug targets in P. falciparum genome using a sequence of computational methods from Malaria Parasite Metabolic Pathway database. The study reported the identification of 71 proteins as potential drug targets for anti-malarial interventions. Furthermore, homology modeling and molecular dynamic simulation study of one of the potential drug targets, aminodeoxychorismate lyase, was carried to predict the 3D structure of the protein. Structure and ligand-based drug designing reported MMV019742 from Pathogen Box and TCAMS-141515 from GSK-TCAMS library as potential hits. The reliability of the binding mode of the inhibitors is confirmed by GROMACS for a simulation time of 20 ns in water environment. This will be helpful for experimental validation of the small-molecule inhibitor.

Genome-wide associations of CD46 and IFI44L genetic variants with neutralizing antibody response to measles vaccine

Population-based studies have revealed 2-10% measles vaccine failure rate even after two vaccine doses. While the mechanisms behind this remain unknown, we hypothesized that host genetic factors are likely to be involved. We performed a genome-wide association study of measles specific neutralizing antibody and IFNγ ELISPOT response in a combined sample of 2872 subjects. We identified two distinct chromosome 1 regions (previously associated with MMR-related febrile seizures), associated with vaccine-induced measles neutralizing antibody titers. The 1q32 region contained 20 significant SNPs in/around the measles virus receptor-encoding CD46 gene, including the intronic rs2724384 (p value = 2.64 × 10^-09) and rs2724374 (p value = 3.16 × 10^-09) SNPs. The 1q31.1 region contained nine significant SNPs in/around IFI44L, including the intronic rs1333973 (p value = 1.41 × 10^-10) and the missense rs273259 (His73Arg, p value = 2.87 × 10^-10) SNPs. Analysis of differential exon usage with mRNA-Seq data and RT-PCR suggests the involvement of rs2724374 minor G allele in the CD46 STP region exon B skipping, resulting in shorter CD46 isoforms. Our study reveals common CD46 and IFI44L SNPs associated with measles-specific humoral immunity, and highlights the importance of alternative splicing/virus cellular receptor isoform usage as a mechanism explaining inter-individual variation in immune response after live measles vaccine.

Whole-genome sequencing identifies common-to-rare variants associated with human blood metabolites

Genetic factors modifying the blood metabolome have been investigated through genome-wide association studies (GWAS) of common genetic variants and through exome sequencing. We conducted a whole-genome sequencing study of common, low-frequency and rare variants to associate genetic variations with blood metabolite levels using comprehensive metabolite profiling in 1,960 adults. We focused the analysis on 644 metabolites with consistent levels across three longitudinal data collections. Genetic sequence variations at 101 loci were associated with the levels of 246 (38%) metabolites (P ≤ 1.9 × 10^-11). We identified 113 (10.7%) among 1,054 unrelated individuals in the cohort who carried heterozygous rare variants likely influencing the function of 17 genes. Thirteen of the 17 genes are associated with inborn errors of metabolism or other pediatric genetic conditions. This study extends the map of loci influencing the metabolome and highlights the importance of heterozygous rare variants in determining abnormal blood metabolic phenotypes in adults.

PB2 substitutions V598T/I increase the virulence of H7N9 influenza A virus in mammals

PB2 is one of the subunits of the influenza A virus (IAV) polymerase complex. By bioinformatics analysis we identified PB2 substitutions at positions 389 and 598 among IAV isolates from humans, which might associate with viral pathogenicity. To evaluate the biological significance of these substitutions, PB2-K389R and -V598T/I mutant viruses of avian H7N9 IAVs were generated by reverse genetics. Compared to the wild type, the mutant viruses displayed an enhanced growth capacity in human and mammalian cells. Meanwhile, they presented increased transcription and replication by producing higher levels of viral mRNA, cRNA and vRNA. Minireplicon assays indicated that the polymerase activity was elevated by these substitutions. Notably, the PB2-V598T/I substitutions substantially increased virus replication and virulence in mice. Together, we demonstrated that the substitutions PB2-V598T/I contributed to higher IAV replication and virulence in mammals, which added to the knowledge of IAV virulence determinants and benefited the surveillance of IAVs.

UNcleProt (Universal Nuclear Protein database of barley): The first nuclear protein database that distinguishes proteins from different phases of the cell cycle

Proteins are the most abundant component of the cell nucleus, where they perform a plethora of functions, including the assembly of long DNA molecules into condensed chromatin, DNA replication and repair, regulation of gene expression, synthesis of RNA molecules and their modification. Proteins are important components of nuclear bodies and are involved in the maintenance of the nuclear architecture, transport across the nuclear envelope and cell division. Given their importance, the current poor knowledge of plant nuclear proteins and their dynamics during the cell's life and division is striking. Several factors hamper the analysis of the plant nuclear proteome, but the most critical seems to be the contamination of nuclei by cytosolic material during their isolation. With the availability of an efficient protocol for the purification of plant nuclei, based on flow cytometric sorting, contamination by cytoplasmic remnants can be minimized. Moreover, flow cytometry allows the separation of nuclei in different stages of the cell cycle (G1, S, and G2). This strategy has led to the identification of large number of nuclear proteins from barley (Hordeum vulgare), thus triggering the creation of a dedicated database called UNcleProt, http://barley.gambrinus.ueb.cas.cz/.

Osteocalcin and Sex Hormone Binding Globulin Compete on a Specific Binding Site of GPRC6A

The undercarboxylated form of osteocalcin (ucOC) regulates male fertility and energy metabolism, acting through the G protein-coupled receptor (GPRC)6A, thus forming a new pancreas-bone-testis axis. Recently, GPRC6A has also been suggested to mediate the nongenomic responses of free testosterone (T). However, these data did not consider the physiological scenario, where circulating T is mainly bound to sex hormone-binding globulin (SHBG) and only a small percentage circulates freely in the blood. Here, by the use of computational modelling, we document the existence of similar structural moieties between ucOC and SHBG that are predicted to bind to GPRC6A at docking analysis. This hypothesis of competition was assessed by binding experiments on human embryonic kidney-293 cells transfected with human GPRC6A gene. Unliganded SHBG specifically bound the membrane of human embryonic kidney-293 cells transfected with GPRC6A and was displaced by ucOC when coincubated at 100-fold molar excess. Furthermore, specific downstream Erk1/2 phosphorylation after stimulation of GPRC6A with ucOC was significantly blunted by 100-fold molar excess of unliganded SHBG. Intriguingly previous incubation with unliganded SHBG, followed by incubation with T, induced Erk1/2 phosphorylation in a dose-dependent manner. Neither binding nor stimulating activities were shown for SHBG saturated with T. Experiments on mutation constructs of GPRC6A strengthened the hypothesis of a common binding site of ucOC and SHBG. Given the role of GPRC6A on energy metabolism, these data agree with epidemiological association between SHBG levels and insulin sensitivity, suggest GPRC6A as a likely SHBG receptor, and add bases for the possible regulation of androgen activity in a nonsteroidal manner.

In-silico structural analysis of E509K mutation in LARGE and T192M mutation in Alpha Dystroglycan in the inhibition of glycosylation of Alpha Dystroglycan by LARGE

Impaired glycosylation of cellular receptor Alpha Dystroglycan (α-DG) leads to dystroglycanopathy. Glycoprotein α-DG is the receptor protein in the Dystrophin Associated Protein Complex (DAPC), a macromolecular gathering on muscle cell membrane to form a bridge between extracellular matrix (ECM) and cellular actin cytoskeleton. Proper glycosylation of α-DG is mediated by the glycosylating enzyme LARGE. Mutations either in α-DG or in LARGE lead to improper glycosylations of α-DG thereby hampering the formation of final Laminin binding form α-DG resulting in dystroglycanopathy. In our current work, we explored the structural changes associated with the presence of mutations in α-DG as well as in the enzyme LARGE. We further extended our research to understand the effect of the mutations onto protein-enzyme interactions. Moreover, since LARGE transfers the sugar moiety (glucuronic acid; GlcA) onto α-DG, we tried to analyze what effect the mutation in LARGE confers on this enzyme ligand interaction. This work for the first time addressed the molecular changes occurring in the structures α-DG, LARGE and their interactions and shed lights on the as yet poorly understood mechanism behind the dystroglycanopathy onset.

BALL-SNPgp-from genetic variants toward computational diagnostics

In medical research, it is crucial to understand the functional consequences of genetic alterations, for example, non-synonymous single nucleotide variants (nsSNVs). NsSNVs are known to be causative for several human diseases. However, the genetic basis of complex disorders such as diabetes or cancer comprises multiple factors. Methods to analyze putative synergetic effects of multiple such factors, however, are limited. Here, we concentrate on nsSNVs and present BALL-SNPgp, a tool for structural and functional characterization of nsSNVs, which is aimed to improve pathogenicity assessment in computational diagnostics. Based on annotated SNV data, BALL-SNPgp creates a three-dimensional visualization of the encoded protein, collects available information from different resources concerning disease relevance and other functional annotations, performs cluster analysis, predicts putative binding pockets and provides data on known interaction sites.

AVAILABILITY AND IMPLEMENTATION

BALL-SNPgp is based on the comprehensive C ++ framework Biochemical Algorithms Library (BALL) and its visualization front-end BALLView. Our tool is available at www.ccb.uni-saarland.de/BALL-SNPgp

CONTACT

ballsnp@milaman.cs.uni-saarland.de.

Negative regulation of bacterial killing and inflammation by two novel CD16 ligands

Sepsis, a leading cause of death worldwide, involves exacerbated proinflammatory responses and inefficient bacterial clearance. Phagocytic cells play a crucial part in the prevention of sepsis by clearing bacteria through host innate receptors. Here, we used a phage display library to identify two peptides in Escherichia coli that interact with host innate receptors. One of these peptides, encoded by the wzxE gene of E. coli K-12, was involved in the transbilayer movement of a trisaccharide-lipid intermediate in the assembly of enterobacterial common antigen. Peptide-receptor interactions induced CD16-mediated inhibitory immunoreceptor tyrosine-based activating motif signaling, blocking the production of ROS and bacterial killing. This CD16-mediated inhibitory signaling was abrogated in a WzxE(-/-) mutant of E. coli K-12, restoring the production of ROS and bacterial killing. Taken together, the two novel CD16 ligands identified negatively regulate bacterial killing and inflammation. Our findings may contribute toward the development of new immunotherapies for E. coli-mediated infectious diseases and inflammation.

Computational Framework for Prediction of Peptide Sequences That May Mediate Multiple Protein Interactions in Cancer-Associated Hub Proteins

A considerable proportion of protein-protein interactions (PPIs) in the cell are estimated to be mediated by very short peptide segments that approximately conform to specific sequence patterns known as linear motifs (LMs), often present in the disordered regions in the eukaryotic proteins. These peptides have been found to interact with low affinity and are able bind to multiple interactors, thus playing an important role in the PPI networks involving date hubs. In this work, PPI data and de novo motif identification based method (MEME) were used to identify such peptides in three cancer-associated hub proteins—MYC, APC and MDM2. The peptides corresponding to the significant LMs identified for each hub protein were aligned, the overlapping regions across these peptides being termed as overlapping linear peptides (OLPs). These OLPs were thus predicted to be responsible for multiple PPIs of the corresponding hub proteins and a scoring system was developed to rank them. We predicted six OLPs in MYC and five OLPs in MDM2 that scored higher than OLP predictions from randomly generated protein sets. Two OLP sequences from the C-terminal of MYC were predicted to bind with FBXW7, component of an E3 ubiquitin-protein ligase complex involved in proteasomal degradation of MYC. Similarly, we identified peptides in the C-terminal of MDM2 interacting with FKBP3, which has a specific role in auto-ubiquitinylation of MDM2. The peptide sequences predicted in MYC and MDM2 look promising for designing orthosteric inhibitors against possible disease-associated PPIs. Since these OLPs can interact with other proteins as well, these inhibitors should be specific to the targeted interactor to prevent undesired side-effects. This computational framework has been designed to predict and rank the peptide regions that may mediate multiple PPIs and can be applied to other disease-associated date hub proteins for prediction of novel therapeutic targets of small molecule PPI modulators.

ConBind: motif-aware cross-species alignment for the identification of functional transcription factor binding sites

Eukaryotic gene expression is regulated by transcription factors (TFs) binding to promoter as well as distal enhancers. TFs recognize short, but specific binding sites (TFBSs) that are located within the promoter and enhancer regions. Functionally relevant TFBSs are often highly conserved during evolution leaving a strong phylogenetic signal. While multiple sequence alignment (MSA) is a potent tool to detect the phylogenetic signal, the current MSA implementations are optimized to align the maximum number of identical nucleotides. This approach might result in the omission of conserved motifs that contain interchangeable nucleotides such as the ETS motif (IUPAC code: GGAW). Here, we introduce ConBind, a novel method to enhance alignment of short motifs, even if their mutual sequence similarity is only partial. ConBind improves the identification of conserved TFBSs by improving the alignment accuracy of TFBS families within orthologous DNA sequences. Functional validation of the Gfi1b + 13 enhancer reveals that ConBind identifies additional functionally important ETS binding sites that were missed by all other tested alignment tools. In addition to the analysis of known regulatory regions, our web tool is useful for the analysis of TFBSs on so far unknown DNA regions identified through ChIP-sequencing.

Antiprotozoal Nitazoxanide Derivatives: Synthesis, Bioassays and QSAR Study Combined with Docking for Mechanistic Insight

In view of the serious health problems concerning infectious diseases in heavily populated areas, we followed the strategy of lead compound diversification to evaluate the near-by chemical space for new organic compounds. To this end, twenty derivatives of nitazoxanide (NTZ) were synthesized and tested for activity against Entamoeba histolytica parasites. To ensure drug-likeliness and activity relatedness of the new compounds, the synthetic work was assisted by a quantitative structure-activity relationships study (QSAR). Many of the inherent downsides – well-known to QSAR practitioners – we circumvented thanks to workarounds which we proposed in prior QSAR publication. To gain further mechanistic insight on a molecular level, ligand-enzyme docking simulations were carried out since NTZ is known to inhibit the protozoal pyruvate ferredoxin oxidoreductase (PFOR) enzyme as its biomolecular target.

Building a virtual ligand screening pipeline using free software: a survey

Virtual screening, the search for bioactive compounds via computational methods, provides a wide range of opportunities to speed up drug development and reduce the associated risks and costs. While virtual screening is already a standard practice in pharmaceutical companies, its applications in preclinical academic research still remain under-exploited, in spite of an increasing availability of dedicated free databases and software tools. In this survey, an overview of recent developments in this field is presented, focusing on free software and data repositories for screening as alternatives to their commercial counterparts, and outlining how available resources can be interlinked into a comprehensive virtual screening pipeline using typical academic computing facilities. Finally, to facilitate the set-up of corresponding pipelines, a downloadable software system is provided, using platform virtualization to integrate pre-installed screening tools and scripts for reproducible application across different operating systems.

Robust classification of protein variation using structural modelling and large-scale data integration

Existing methods for interpreting protein variation focus on annotating mutation pathogenicity rather than detailed interpretation of variant deleteriousness and frequently use only sequence-based or structure-based information. We present VIPUR, a computational framework that seamlessly integrates sequence analysis and structural modelling (using the Rosetta protein modelling suite) to identify and interpret deleterious protein variants. To train VIPUR, we collected 9477 protein variants with known effects on protein function from multiple organisms and curated structural models for each variant from crystal structures and homology models. VIPUR can be applied to mutations in any organism's proteome with improved generalized accuracy (AUROC .83) and interpretability (AUPR .87) compared to other methods. We demonstrate that VIPUR's predictions of deleteriousness match the biological phenotypes in ClinVar and provide a clear ranking of prediction confidence. We use VIPUR to interpret known mutations associated with inflammation and diabetes, demonstrating the structural diversity of disrupted functional sites and improved interpretation of mutations associated with human diseases. Lastly, we demonstrate VIPUR's ability to highlight candidate variants associated with human diseases by applying VIPUR to de novo variants associated with autism spectrum disorders.

Comparative analyses of the proteins from Mycobacterium tuberculosis and human genomes: Identification of potential tuberculosis drug targets

Tuberculosis, one of the major infectious diseases affecting human beings is caused by the bacillus Mycobacterium tuberculosis. Increased resistance to known drugs commonly used for the treatment of tuberculosis has created an urgent need to identify new targets for validation and to develop drugs. In this study, we have used various bioinformatics tools, to compare the protein sequences from twenty-three M. tuberculosis genome strains along with the known human protein sequences, in order to identify the 'conserved' M. tuberculosis proteins absent in human. Further, based on the analysis of protein interaction networks, we selected one-hundred and forty proteins that were predicted as potential M. tuberculosis drug targets and prioritized according to the ranking of 'clusters' of interacting proteins. Comparison of the predicted 140 TB targets with literature indicated that 46 of them were previously reported, thereby increasing the confidence in our predictions of the remaining 94 targets too. The analyses of the structures and functions corresponding to the predicted potential TB drug targets indicated a diverse range of proteins that included ten 'druggable' targets with some of the known drugs.

The neuronal porosome complex in health and disease

Cup-shaped secretory portals at the cell plasma membrane called porosomes mediate the precision release of intravesicular material from cells. Membrane-bound secretory vesicles transiently dock and fuse at the base of porosomes facing the cytosol to expel pressurized intravesicular contents from the cell during secretion. The structure, isolation, composition, and functional reconstitution of the neuronal porosome complex have greatly progressed, providing a molecular understanding of its function in health and disease. Neuronal porosomes are 15 nm cup-shaped lipoprotein structures composed of nearly 40 proteins, compared to the 120 nm nuclear pore complex composed of >500 protein molecules. Membrane proteins compose the porosome complex, making it practically impossible to solve its atomic structure. However, atomic force microscopy and small-angle X-ray solution scattering studies have provided three-dimensional structural details of the native neuronal porosome at sub-nanometer resolution, providing insights into the molecular mechanism of its function. The participation of several porosome proteins previously implicated in neurotransmission and neurological disorders, further attest to the crosstalk between porosome proteins and their coordinated involvement in release of neurotransmitter at the synapse.

Comparisons of Allergenic and Metazoan Parasite Proteins: Allergy the Price of Immunity

Allergic reactions can be considered as maladaptive IgE immune responses towards environmental antigens. Intriguingly, these mechanisms are observed to be very similar to those implicated in the acquisition of an important degree of immunity against metazoan parasites (helminths and arthropods) in mammalian hosts. Based on the hypothesis that IgE-mediated immune responses evolved in mammals to provide extra protection against metazoan parasites rather than to cause allergy, we predict that the environmental allergens will share key properties with the metazoan parasite antigens that are specifically targeted by IgE in infected human populations. We seek to test this prediction by examining if significant similarity exists between molecular features of allergens and helminth proteins that induce an IgE response in the human host. By employing various computational approaches, 2712 unique protein molecules that are known IgE antigens were searched against a dataset of proteins from helminths and parasitic arthropods, resulting in a comprehensive list of 2445 parasite proteins that show significant similarity through sequence and structure with allergenic proteins. Nearly half of these parasite proteins from 31 species fall within the 10 most abundant allergenic protein domain families (EF-hand, Tropomyosin, CAP, Profilin, Lipocalin, Trypsin-like serine protease, Cupin, BetV1, Expansin and Prolamin). We identified epitopic-like regions in 206 parasite proteins and present the first example of a plant protein (BetV1) that is the commonest allergen in pollen in a worm, and confirming it as the target of IgE in schistosomiasis infected humans. The identification of significant similarity, inclusive of the epitopic regions, between allergens and helminth proteins against which IgE is an observed marker of protective immunity explains the ‘off-target’ effects of the IgE-mediated immune system in allergy. All these findings can impact the discovery and design of molecules used in immunotherapy of allergic conditions.

Allergy is an increasingly widespread clinical problem that leads to various conditions such as allergic asthma and susceptibility to anaphylactic shock. These conditions arise from exposure to a range of environmental and food proteins (‘allergens’) that are recognised by a form of immune system antibody called IgE. This part of the immune system is thought to have evolved to provide mammals with additional rapid response mechanisms to combat metazoan parasites. Here, we address the pertinent question, ‘what makes an Allergen an Allergen’ as, although they constitute a very small percentage of known proteins, they appear to be diverse and unrelated. Using computational studies, we have established molecular similarity between parasite proteins and allergens that affect the nature of immune response and are able to predict the regions of parasite proteins that potentially share similarity with the IgE-binding region(s) of the allergens. Our experimental studies support the computational predictions, and we can present the first confirmed example of a plant pollen-like protein in a worm that is targeted by IgE. The results of this study will enable us to predict likely allergens in food and environmental organisms and to help design protein molecules to treat allergy in the future.

Lhx6 and Lhx8 promote palate development through negative regulation of a cell cycle inhibitor gene, p57Kip2

Cleft palate is a common birth defect in humans. Therefore, understanding the molecular genetics of palate development is important from both scientific and medical perspectives. Lhx6 and Lhx8 encode LIM homeodomain transcription factors, and inactivation of both genes in mice resulted in profound craniofacial defects including cleft secondary palate. The initial outgrowth of the palate was severely impaired in the mutant embryos, due to decreased cell proliferation. Through genome-wide transcriptional profiling, we discovered that p57(Kip2) (Cdkn1c), encoding a cell cycle inhibitor, was up-regulated in the prospective palate of Lhx6(-/-);Lhx8(-/-) mutants. p57(Kip2) has been linked to Beckwith-Wiedemann syndrome and IMAGe syndrome in humans, which are developmental disorders with increased incidents of palate defects among the patients. To determine the molecular mechanism underlying the regulation of p57(Kip2) by the Lhx genes, we combined chromatin immunoprecipitation, in silico search for transcription factor-binding motifs, and in vitro reporter assays with putative cis-regulatory elements. The results of these experiments indicated that LHX6 and LHX8 regulated p57(Kip2) via both direct and indirect mechanisms, with the latter mediated by Forkhead box (FOX) family transcription factors. Together, our findings uncovered a novel connection between the initiation of palate development and a cell cycle inhibitor via LHX. We propose a model in which Lhx6 and Lhx8 negatively regulate p57(Kip2) expression in the prospective palate area to allow adequate levels of cell proliferation and thereby promote normal palate development. This is the first report elucidating a molecular genetic pathway downstream of Lhx in palate development.

Computational Studies of Benzoxazinone Derivatives as Antiviral Agents against Herpes Virus Type 1 Protease

Herpes simplex virus infections have been described in the medical literature for centuries, yet the the drugs available nowadays for therapy are largely ineffective and low oral bioavailability plays an important role on the inefficacy of the treatments. Additionally, the details of the inhibition of Herpes Virus type 1 are still not fully understood. Studies have shown that several viruses encode one or more proteases required for the production new infectious virions. This study presents an analysis of the interactions between HSV-1 protease and benzoxazinone derivatives through a combination of structure-activity relationships, comparative modeling and molecular docking studies. The structure activity relationship results showed an important contribution of hydrophobic and polarizable groups and limitations for bulky groups in specific positions. Two Herpes Virus type 1 protease models were constructed and compared to achieve the best model which was obtained by MODELLER. Molecular docking results pointed to an important interaction between the most potent benzoxazinone derivative and Ser129, consistent with previous mechanistic data. Moreover, we also observed hydrophobic interactions that may play an important role in the stabilization of inhibitors in the active site. Finally, we performed druglikeness and drugscore studies of the most potent derivatives and the drugs currently used against Herpes virus.

Allosteric regulation of focal adhesion kinase by PIP₂ and ATP

Focal adhesion kinase (FAK) is a nonreceptor tyrosine kinase that regulates cell signaling, proliferation, migration, and development. A major mechanism of regulation of FAK activity is an intramolecular autoinhibitory interaction between two of its domains--the catalytic and FERM domains. Upon cell adhesion to the extracellular matrix, FAK is being translocated toward focal adhesion sites and activated. Interactions of FAK with phosphoinositide phosphatidylinsositol-4,5-bis-phosphate (PIP₂) are required to activate FAK. However, the molecular mechanism of the activation remains poorly understood. Recent fluorescence resonance energy transfer experiments revealed a closure of the FERM-kinase interface upon ATP binding, which is reversed upon additional binding of PIP₂. Here, we addressed the allosteric regulation of FAK by performing all-atom molecular-dynamics simulations of a FAK fragment containing the catalytic and FERM domains, and comparing the dynamics in the absence or presence of ATP and PIP₂. As a major conformational change, we observe a closing and opening motion upon ATP and additional PIP₂ binding, respectively, in good agreement with the fluorescence resonance energy transfer experiments. To reveal how the binding of the regulatory PIP₂ to the FERM F2 lobe is transduced to the very distant F1/N-lobe interface, we employed force distribution analysis. We identified a network of mainly charged residue-residue interactions spanning from the PIP₂ binding site to the distant interface between the kinase and FERM domains, comprising candidate residues for mutagenesis to validate the predicted mechanism of FAK activation.

Sequence and structural difference favors a distinct preference of Wnt3a binding with co-receptor LRP6

Wnt signaling pathway plays a key role in a wide array of development and physiological processes. Wnt proteins interact with two different co-receptors LRP5/6 and ROR 2, leading to different signal transductions in the cell. Though the Wnt family of proteins has high sequence similarity the specificity for particular co-receptor is not well understood. The choice of pathway is attributed to the binding of Wnt complex to the co-receptor. Our current study is a novel approach using homology modeling, docking, and structural alignment to unravel the structural differences between Wnt3a and Wnt5b binding to LRP6. The conservation of a protruding loop has been identified in Wnt3a protein indicating an enhanced ability of Wnt3a to bind to LRP5/6 against its counter parts. The docking studies have further substantiated the findings. This could potentially help us design and develop novel inhibitors targeting Wnt3a-LRP6 complex in specific tissues or disease states.

Methods for generating and applying pharmacophore models as virtual screening filters and for bioactivity profiling

Biological effects of small molecules in an organism result from favorable interactions between the molecules and their target proteins. These interactions depend on chemical functionalities, bonds, and their 3D-orientations towards each other. These 3D-arrangements of chemical functionalities that make a small molecule active towards its target can be described by pharmacophore models. In these models, chemical functionalities are represented as so-called features. Commonly, they are obtained either from a set of active compounds or directly from the observed protein-ligand interactions as present in X-ray crystal structures, NMR structures, or docking poses. In this review, we explain the basics of pharmacophore modeling including dataset generation, 3D-representations and conformational analysis of small molecules, pharmacophore model construction, model validation, and its benefits to virtual screening and other applications.

Neuronal Porosome-The Secretory Portal at the Nerve Terminal: It's Structure-Function, Composition, and Reconstitution

Cup-shaped secretory portals at the cell plasma membrane called porosomes mediate secretion from cells. Membrane bound secretory vesicles transiently dock and fuse at the cytosolic compartment of the porosome base to expel intravesicular contents to the outside during cell secretion. In the past decade, the structure, isolation, composition, and functional reconstitution of the neuronal porosome complex has been accomplished providing a molecular understanding of its structure-function. Neuronal porosomes are 15 nm cup-shaped lipoprotein structures composed of nearly 40 proteins. Being a membrane-associated supramolecular complex has precluded determination of the atomic structure of the porosome. However recent studies using small-angle X-ray solution scattering (SAXS), provide at sub-nanometer resolution, the native 3D structure of the neuronal porosome complex associated with docked synaptic vesicle at the nerve terminal. Additionally, results from the SAXS study and earlier studies using atomic force microscopy, provide the possible molecular mechanism involved in porosome-mediated neurotransmitter release at the nerve terminal.

In-silico characterization of Formin Binding Protein 4 Family of proteins

Members of the Formin Binding Protein 4 Family or the FNBP4 were indirectly reported to be associated with many of the biological processes. These proteins possess two WW domains. So far there are practically no reports regarding the characterization and classification of the protein by any means. Keeping in mind the importance of the proteins from this FNBP4 family, we have tried an in silico approach to come up with a comprehensive analysis of the proteins. We have analyzed the proteins by considering their sequence conservation, their phylogenetic distributions among the different organisms. We have also investigated the functional properties of the WW domains in the proteins. Finally, we have made an attempt to elucidate the structural details of the domains and predicted the possible modes of their interactions. Our findings show that FNBP4 is eukaryotic in its distribution and follows a trend of evolution where animal and plant homologues have evolved in an independent manner. While the WW domain is the only common motif present across the FNBP4 family of proteins, there are different classes (mainly two) of WW domains that are found among different FNBP4 proteins. Structure function predictions indicate a possible role of FNBP4 in either protein stabilization control or transcript processing. Our study on FNBP4 may therefore open up new avenues to generate new interest in this highly important but largely unexplored class of proteins. Future studies with proteins from this family may answer many important questions of protein-protein interactions in different biologically important processes.

FixingTIM: interactive exploration of sequence and structural data to identify functional mutations in protein families

BACKGROUND

Knowledge of the 3D structure and functionality of proteins can lead to insight into the associated cellular processes, speed up the creation of pharmaceutical products, and develop drugs that are more effective in combating disease.

METHODS

We present the design and implementation of a visual mining and analysis tool to help identify protein mutations across a family of structural models and to help discover the effect of these mutations on protein function. We integrate 3D structure and sequence information in a common visual interface; multiple linked views and a computational backbone allow comparison at the molecular and atomic levels, while a novel trend-image visual abstraction allows for the sorting and mining of large collections of sequences and of their residues.

RESULTS

We evaluate our approach on the triosephosphate isomerase (TIM) family structural models and sequence data and show that our tool provides an effective, scalable way to navigate a family of proteins, as well as a means to inspect the structure and sequence of individual proteins.

CONCLUSIONS

The TIM application shows that our tool can assist in the navigation of families of proteins, as well as in the exploration of individual protein structures. In conjunction with domain expert knowledge, this interactive tool can help provide biophysical insight into why specific mutations affect function and potentially suggest additional modifications to the protein that could be used to rescue functionality.

Comparative modeling and virtual screening for the identification of novel inhibitors for myo-inositol-1-phosphate synthase

Myo-inositol-1-phosphate (MIP) synthase is a key enzyme in the myo-inositol biosynthesis pathway. Disruption of the inositol signaling pathway is associated with bipolar disorders. Previous work suggested that MIP synthase could be an attractive target for the development of anti-bipolar drugs. Inhibition of this enzyme could possibly help in reducing the risk of a disease in patients. With this objective, three dimensional structure of the protein was modeled followed by the active site prediction. For the first time, computational studies were carried out to obtain structural insights into the interactive behavior of this enzyme with ligands. Virtual screening was carried out using FILTER, ROCS and EON modules of the OpenEye scientific software. Natural products from the ZINC database were used for the screening process. Resulting compounds were docked into active site of the target protein using FRED (Fast Rigid Exhaustive Docking) and GOLD (Genetic Optimization for Ligand Docking) docking programs. The analysis indicated extensive hydrogen bonding network and hydrophobic interactions which play a significant role in ligand binding. Four compounds are shortlisted and their binding assay analysis is underway.

Smotifs as structural local descriptors of supersecondary elements: classification, completeness and applications

Protein structures are made up of periodic and aperiodic structural elements (i.e., α-helices, β-strands and loops). Despite the apparent lack of regular structure, loops have specific conformations and play a central role in the folding, dynamics, and function of proteins. In this article, we reviewed our previous works in the study of protein loops as local supersecondary structural motifs or Smotifs. We reexamined our works about the structural classification of loops (ArchDB) and its application to loop structure prediction (ArchPRED), including the assessment of the limits of knowledge-based loop structure prediction methods. We finalized this article by focusing on the modular nature of proteins and how the concept of Smotifs provides a convenient and practical approach to decompose proteins into strings of concatenated Smotifs and how can this be used in computational protein design and protein structure prediction.

Protein structure modeling in the proteomics era

Structural proteomics aims to understand the structural basis of protein interactions and functions. A prerequisite for this is the availability of 3D protein structures that mediate the biochemical interactions. The explosion in the number of available gene sequences set the stage for the next step in genome-scale projects -- to obtain 3D structures for each protein. To achieve this ambitious goal, the slow and costly structure determination experiments are supplemented with theoretical approaches. The current state and recent advances in structure modeling approaches are reviewed here, with special emphasis on comparative protein structure modeling techniques.

Hypoglycosylation of dystroglycan due to T192M mutation: a molecular insight behind the fact

Abnormal glycosylation of dystroglycan (DG), a transmembrane glycoprotein, results in a group of diseases known as dystroglycanopathy. A severe dystroglycanopathy known as the limb girdle disease MDDGC9 [OMIM: 613818] occurs as a result of hypoglycosylation of alpha subunit of DG. Reasons behind this has been traced back to a point mutation (T192M) in DG that leads to weakening of interactions of DG protein with laminin and subsequent loss of signal flow through the DG protein. In this work we have tried to analyze the molecular details of the interactions between DG and laminin1 in order to propose a mechanism about the onset of the disease MDDGC9. We have observed noticeable changes between the modeled structures of wild type and mutant DG proteins. We also have employed molecular docking techniques to study and compare the binding interactions between laminin1 and both the wild type and mutant DG proteins. The docking simulations have revealed that the mutant DG has weaker interactions with laminin1 as compared to the wild type DG. Till date there are no previous reports that deal with the elucidation of the interactions of DG with laminin1 from the molecular level. Our study is therefore the first of its kind which analyzes the differences in binding patterns of laminin1 with both the wild type and mutant DG proteins. Our work would therefore facilitate analysis of the molecular mechanism of the disease MDDGC9. Future work based on our results may be useful for the development of suitable drugs against this disease.

Solving nucleic acid structures by molecular replacement: examples from group II intron studies

Structured RNA molecules are key players in ensuring cellular viability. It is now emerging that, like proteins, the functions of many nucleic acids are dictated by their tertiary folds. At the same time, the number of known crystal structures of nucleic acids is also increasing rapidly. In this context, molecular replacement will become an increasingly useful technique for phasing nucleic acid crystallographic data in the near future. Here, strategies to select, create and refine molecular-replacement search models for nucleic acids are discussed. Using examples taken primarily from research on group II introns, it is shown that nucleic acids are amenable to different and potentially more flexible and sophisticated molecular-replacement searches than proteins. These observations specifically aim to encourage future crystallographic studies on the newly discovered repertoire of noncoding transcripts.

Protein models: the Grand Challenge of protein docking

Characterization of life processes at the molecular level requires structural details of protein-protein interactions (PPIs). The number of experimentally determined protein structures accounts only for a fraction of known proteins. This gap has to be bridged by modeling, typically using experimentally determined structures as templates to model related proteins. The fraction of experimentally determined PPI structures is even smaller than that for the individual proteins, due to a larger number of interactions than the number of individual proteins, and a greater difficulty of crystallizing protein-protein complexes. The approaches to structural modeling of PPI (docking) often have to rely on modeled structures of the interactors, especially in the case of large PPI networks. Structures of modeled proteins are typically less accurate than the ones determined by X-ray crystallography or nuclear magnetic resonance. Thus the utility of approaches to dock these structures should be assessed by thorough benchmarking, specifically designed for protein models. To be credible, such benchmarking has to be based on carefully curated sets of structures with levels of distortion typical for modeled proteins. This article presents such a suite of models built for the benchmark set of the X-ray structures from the Dockground resource (http://dockground.bioinformatics.ku.edu) by a combination of homology modeling and Nudged Elastic Band method. For each monomer, six models were generated with predefined C(α) root mean square deviation from the native structure (1, 2, …, 6 Å). The sets and the accompanying data provide a comprehensive resource for the development of docking methodology for modeled proteins.

An estimated 5% of new protein structures solved today represent a new Pfam family

High-resolution structural knowledge is key to understanding how proteins function at the molecular level. The number of entries in the Protein Data Bank (PDB), the repository of all publicly available protein structures, continues to increase, with more than 8000 structures released in 2012 alone. The authors of this article have studied how structural coverage of the protein-sequence space has changed over time by monitoring the number of Pfam families that acquired their first representative structure each year from 1976 to 2012. Twenty years ago, for every 100 new PDB entries released, an estimated 20 Pfam families acquired their first structure. By 2012, this decreased to only about five families per 100 structures. The reasons behind the slower pace at which previously uncharacterized families are being structurally covered were investigated. It was found that although more than 50% of current Pfam families are still without a structural representative, this set is enriched in families that are small, functionally uncharacterized or rich in problem features such as intrinsically disordered and transmembrane regions. While these are important constraints, the reasons why it may not yet be time to give up the pursuit of a targeted but more comprehensive structural coverage of the protein-sequence space are discussed.

Stability Mechanisms of Laccase Isoforms using a Modified FoldX Protocol Applicable to Widely Different Proteins

A recent computational protocol that accurately predicts and rationalizes protein multisite mutant stabilities has been extended to handle widely different isoforms of laccases. We apply the protocol to four isoenzymes of Trametes versicolor laccase (TvL) with variable lengths (498-503 residues) and thermostability (Topt ∼ 45-80 °C) and with 67-77% sequence identity. The extended protocol uses (i) statistical averaging, (ii) a molecular-dynamics-validated "compromise" homology model to minimize bias that causes proteins close in sequence to a structural template to be too stable due to having the benefits of the better sampled template (typically from a crystal structure), (iii) correction for hysteresis that favors the input template to overdestabilize, and (iv) a preparative protocol to provide robust input sequences of equal length. The computed ΔΔG values are in good agreement with the major trends in experimental stabilities; that is, the approach may be applicable for fast estimates of the relative stabilities of proteins with as little as 70% identity, something that is currently extremely challenging. The computed stability changes associated with variations are Gaussian-distributed, in good agreement with experimental distributions of stability effects from mutation. The residues causing the differential stability of the four isoforms are consistent with a range of compiled laccase wild type data, suggesting that we may have identified general drivers of laccase stability. Several sites near Cu, notably 79, 241, and 245, or near substrate, mainly 265, are identified that contribute to stability-function trade-offs, of relevance to the search for new proficient and stable variants of these important industrial enzymes.

Global analysis of phosphorylation and ubiquitylation cross-talk in protein degradation

Cross-talk between different types of post-translational modifications on the same protein molecule adds specificity and combinatorial logic to signal processing, but it has not been characterized on a large-scale basis. We developed two methods to identify protein isoforms that are both phosphorylated and ubiquitylated in the yeast Saccharomyces cerevisiae, identifying 466 proteins with 2,100 phosphorylation sites co-occurring with 2,189 ubiquitylation sites. We applied these methods quantitatively to identify phosphorylation sites that regulate protein degradation via the ubiquitin-proteasome system. Our results demonstrate that distinct phosphorylation sites are often used in conjunction with ubiquitylation and that these sites are more highly conserved than the entire set of phosphorylation sites. Finally, we investigated how the phosphorylation machinery can be regulated by ubiquitylation. We found evidence for novel regulatory mechanisms of kinases and 14-3-3 scaffold proteins via proteasome-independent ubiquitylation.

Modulation of TET2 expression and 5-methylcytosine oxidation by the CXXC domain protein IDAX

TET (ten-eleven-translocation) proteins are Fe(ii)- and α-ketoglutarate-dependent dioxygenases that modify the methylation status of DNA by successively oxidizing 5-methylcytosine to 5-hydroxymethylcytosine, 5-formylcytosine and 5-carboxycytosine, potential intermediates in the active erasure of DNA-methylation marks. Here we show that IDAX (also known as CXXC4), a reported inhibitor of Wnt signalling that has been implicated in malignant renal cell carcinoma and colonic villous adenoma, regulates TET2 protein expression. IDAX was originally encoded within an ancestral TET2 gene that underwent a chromosomal gene inversion during evolution, thus separating the TET2 CXXC domain from the catalytic domain. The IDAX CXXC domain binds DNA sequences containing unmethylated CpG dinucleotides, localizes to promoters and CpG islands in genomic DNA and interacts directly with the catalytic domain of TET2. Unexpectedly, IDAX expression results in caspase activation and TET2 protein downregulation, in a manner that depends on DNA binding through the IDAX CXXC domain, suggesting that IDAX recruits TET2 to DNA before degradation. IDAX depletion prevents TET2 downregulation in differentiating mouse embryonic stem cells, and short hairpin RNA against IDAX increases TET2 protein expression in the human monocytic cell line U937. Notably, we find that the expression and activity of TET3 is also regulated through its CXXC domain. Taken together, these results establish the separate and linked CXXC domains of TET2 and TET3, respectively, as previously unknown regulators of caspase activation and TET enzymatic activity.

On the immortality of television sets: "function" in the human genome according to the evolution-free gospel of ENCODE

A recent slew of ENCyclopedia Of DNA Elements (ENCODE) Consortium publications, specifically the article signed by all Consortium members, put forward the idea that more than 80% of the human genome is functional. This claim flies in the face of current estimates according to which the fraction of the genome that is evolutionarily conserved through purifying selection is less than 10%. Thus, according to the ENCODE Consortium, a biological function can be maintained indefinitely without selection, which implies that at least 80 - 10 = 70% of the genome is perfectly invulnerable to deleterious mutations, either because no mutation can ever occur in these "functional" regions or because no mutation in these regions can ever be deleterious. This absurd conclusion was reached through various means, chiefly by employing the seldom used "causal role" definition of biological function and then applying it inconsistently to different biochemical properties, by committing a logical fallacy known as "affirming the consequent," by failing to appreciate the crucial difference between "junk DNA" and "garbage DNA," by using analytical methods that yield biased errors and inflate estimates of functionality, by favoring statistical sensitivity over specificity, and by emphasizing statistical significance rather than the magnitude of the effect. Here, we detail the many logical and methodological transgressions involved in assigning functionality to almost every nucleotide in the human genome. The ENCODE results were predicted by one of its authors to necessitate the rewriting of textbooks. We agree, many textbooks dealing with marketing, mass-media hype, and public relations may well have to be rewritten.