POSA: a user-driven, interactive multiple protein structure alignment server

Graphical models for identifying pore-forming proteins

Pore-forming toxins (PFTs) are proteins that form lesions in biological membranes. Better understanding of the structure and function of these proteins will be beneficial in a number of biotechnological applications, including the development of new pest control methods in agriculture. When searching for new pore formers, existing sequence homology-based methods fail to discover truly novel proteins with low sequence identity to known proteins. Search methodologies based on protein structures would help us move beyond this limitation. As the number of known structures for PFTs is very limited, it's quite challenging to identify new proteins having similar structures using computational approaches like deep learning. In this article, we therefore propose a sample-efficient graphical model, where a protein structure graph is first constructed according to consensus secondary structures. A semi-Markov conditional random fields model is then developed to perform protein sequence segmentation. We demonstrate that our method is able to distinguish structurally similar proteins even in the absence of sequence similarity (pairwise sequence identity < 0.4)-a feat not achievable by traditional approaches like HMMs. To extract proteins of interest from a genome-wide protein database for further study, we also develop an efficient framework for UniRef50 with 43 million proteins.

Force-regulated chaperone activity of BiP/ERdj3 is opposite to their homologs DnaK/DnaJ

Polypeptide chains experience mechanical tension while translocating through cellular tunnels, which are subsequently folded by molecular chaperones. However, interactions between tunnel-associated chaperones and these emerging polypeptides under force is not completely understood. Our investigation focused on mechanical chaperone activity of two tunnel-associated chaperones, BiP and ERdj3 both with and without mechanical constraints and comparing them with their cytoplasmic homologs: DnaK and DnaJ. While BiP/ERdj3 have been observed to exhibit robust foldase activity under force, DnaK/DnaJ showed holdase function. Importantly, the tunnel-associated chaperones (BiP/ERdj3) transitioned to a holdase state in the absence of force, indicating a force-dependent chaperone behavior. This chaperone-driven folding event in the tunnel generated an additional mechanical energy of up to 54 zJ, potentially aiding protein translocation. Our findings align with strain theory, where chaperones with higher intrinsic deformability act as mechanical foldases (BiP, ERdj3), while those with lower deformability serve as holdases (DnaK and DnaJ). This study thus elucidates the differential mechanically regulated chaperoning activity and introduces a novel perspective on co-translocational protein folding.

RCSB protein Data Bank: exploring protein 3D similarities via comprehensive structural alignments

MOTIVATION

Tools for pairwise alignments between 3D structures of proteins are of fundamental importance for structural biology and bioinformatics, enabling visual exploration of evolutionary and functional relationships. However, the absence of a user-friendly, browser-based tool for creating alignments and visualizing them at both 1D sequence and 3D structural levels makes this process unnecessarily cumbersome.

RESULTS

We introduce a novel pairwise structure alignment tool (rcsb.org/alignment) that seamlessly integrates into the RCSB Protein Data Bank (RCSB PDB) research-focused RCSB.org web portal. Our tool and its underlying application programming interface (alignment.rcsb.org) empowers users to align several protein chains with a reference structure by providing access to established alignment algorithms (FATCAT, CE, TM-align, or Smith-Waterman 3D). The user-friendly interface simplifies parameter setup and input selection. Within seconds, our tool enables visualization of results in both sequence (1D) and structural (3D) perspectives through the RCSB PDB RCSB.org Sequence Annotations viewer and Mol* 3D viewer, respectively. Users can effortlessly compare structures deposited in the PDB archive alongside more than a million incorporated Computed Structure Models coming from the ModelArchive and AlphaFold DB. Moreover, this tool can be used to align custom structure data by providing a link/URL or uploading atomic coordinate files directly. Importantly, alignment results can be bookmarked and shared with collaborators. By bridging the gap between 1D sequence and 3D structures of proteins, our tool facilitates deeper understanding of complex evolutionary relationships among proteins through comprehensive sequence and structural analyses.

AVAILABILITY AND IMPLEMENTATION

The alignment tool is part of the RCSB PDB research-focused RCSB.org web portal and available at rcsb.org/alignment. Programmatic access is available via alignment.rcsb.org. Frontend code has been published at github.com/rcsb/rcsb-pecos-app. Visualization is powered by the open-source Mol* viewer (github.com/molstar/molstar and github.com/molstar/rcsb-molstar) plus the Sequence Annotations in 3D Viewer (github.com/rcsb/rcsb-saguaro-3d).

Pan-cancer analyses suggest kindlin-associated global mechanochemical alterations

Kindlins serve as mechanosensitive adapters, transducing extracellular mechanical cues to intracellular biochemical signals and thus, their perturbations potentially lead to cancer progressions. Despite the kindlin involvement in tumor development, understanding their genetic and mechanochemical characteristics across different cancers remains elusive. Here, we thoroughly examined genetic alterations in kindlins across more than 10,000 patients with 33 cancer types. Our findings reveal cancer-specific alterations, particularly prevalent in advanced tumor stage and during metastatic onset. We observed a significant co-alteration between kindlins and mechanochemical proteome in various tumors through the activation of cancer-related pathways and adverse survival outcomes. Leveraging normal mode analysis, we predicted structural consequences of cancer-specific kindlin mutations, highlighting potential impacts on stability and downstream signaling pathways. Our study unraveled alterations in epithelial-mesenchymal transition markers associated with kindlin activity. This comprehensive analysis provides a resource for guiding future mechanistic investigations and therapeutic strategies targeting the roles of kindlins in cancer treatment.

Slc11 Synapomorphy: A Conserved 3D Framework Articulating Carrier Conformation Switch

Transmembrane carriers of the Slc11 family catalyze proton (H+)-dependent uptake of divalent metal ions (Me2+) such as manganese and iron-vital elements coveted during infection. The Slc11 mechanism of high-affinity Me2+ cell import is selective and conserved between prokaryotic (MntH) and eukaryotic (Nramp) homologs, though processes coupling the use of the proton motive force to Me2+ uptake evolved repeatedly. Adding bacterial piracy of Nramp genes spread in distinct environmental niches suggests selective gain of function that may benefit opportunistic pathogens. To better understand Slc11 evolution, Alphafold (AF2)/Colabfold (CF) 3D predictions for bacterial sequences from sister clades of eukaryotic descent (MCb and MCg) were compared using both native and mutant templates. AF2/CF model an array of native MCb intermediates spanning the transition from outwardly open (OO) to inwardly open (IO) carriers. In silico mutagenesis targeting (i) a set of (evolutionarily coupled) sites that may define Slc11 function (putative synapomorphy) and (ii) residues from networked communities evolving during MCb transition indicates that Slc11 synapomorphy primarily instructs a Me2+-selective conformation switch which unlocks carrier inner gate and contributes to Me2+ binding site occlusion and outer gate locking. Inner gate opening apparently proceeds from interaction between transmembrane helix (h) h5, h8 and h1a. MCg1 xenologs revealed marked differences in carrier shape and plasticity, owing partly to an altered intramolecular H+ network. Yet, targeting Slc11 synapomorphy also converted MCg1 IO models to an OO state, apparently mobilizing the same residues to control gates. But MCg1 response to mutagenesis differed, with extensive divergence within this clade correlating with MCb-like modeling properties. Notably, MCg1 divergent epistasis marks the emergence of the genus Bordetella-Achromobacter. Slc11 synapomorphy localizes to the 3D areas that deviate least among MCb and MCg1 models (either IO or OO) implying that it constitutes a 3D network of residues articulating a Me2+-selective carrier conformation switch which is maintained in fast-evolving clades at the cost of divergent epistatic interactions impacting carrier shape and dynamics.

Unusual structural and functional features of TpLRR/BspA-like LRR proteins

Leucine Rich Repeat (LRR) domains, are present in hundreds of thousands of proteins across all kingdoms of life and are typically involved in protein-protein interactions and ligand recognition. LRR domains are classified into eight classes and when examined in three dimensions seven, of them form curved solenoid-like super-helices, also described as toruses, with a beta sheet on the concave (inside) and stacked alpha-helices on the convex (outside) of the torus. Here we present an overview of the least characterized 8th class of LRR proteins, the TpLRR-like LRRs, named after the Treponema pallidum protein Tp0225. Proteins from the TpLRR class differ from the proteins in all other known LRR classes by having a flipped curvature, with the beta sheet on the convex side of the torus and irregular secondary structure instead of helices on the opposite, now concave site. TpLRR proteins also present highly divergent sequence pattern of individual repeats and can associate with specific types of additional domains. Several of the characterized proteins from this class, specifically the BspA-like proteins, were found in human bacterial and protozoan pathogens, playing an important role in the interactions between the pathogens and the host immune system. In this paper we surveyed all existing experimental structures and selected AlphaFold models of the best-known proteins containing this class of LRR repeats, analyzing the relation between the pattern of conserved residues, specific structural features and functions of these proteins.

Structural modeling of cytochrome P450 51 from a deep-sea fish points to a novel structural feature in other CYP51s

Cytochromes P450 (CYP), enzymes involved in the metabolism of endogenous and xenobiotic substrates, provide an excellent model system to study how membrane proteins with unique functions have catalytically adapted through evolution. Molecular adaptation of deep-sea proteins to high hydrostatic pressure remains poorly understood. Herein, we have characterized recombinant cytochrome P450 sterol 14α-demethylase (CYP51), an essential enzyme of cholesterol biosynthesis, from an abyssal fish species, Coryphaenoides armatus. C. armatus CYP51 was heterologously expressed in Escherichia coli following N-terminal truncation and purified to homogeneity. Recombinant C. armatus CYP51 bound its sterol substrate lanosterol giving a Type I binding spectra (KD 15 μM) and catalyzed lanosterol 14α-demethylation turnover at 5.8 nmol/min/nmol P450. C. armatus CYP51 also bound the azole antifungals ketoconazole (KD 0.12 μM) and propiconazole (KD 0.54 μM) as determined by Type II absorbance spectra. Comparison of C. armatus CYP51 primary sequence and modeled structures with other CYP51s identified amino acid substitutions that may confer an ability to function under pressures of the deep sea and revealed heretofore undescribed internal cavities in human and other non-deep sea CYP51s. The functional significance of these cavities is not known. PROLOGUE: This paper is dedicated in memory of Michael Waterman and Tsuneo Omura, who as good friends and colleagues enriched our lives. They continue to inspire us.

LncRNA NNT-AS1/hsa-miR-485-5p/HSP90 axis in-silico and clinical prospect correlated-to histologic grades-based CRC stratification: A step toward ncRNA Precision

The oncogenic effects of long non-coding RNA (lncRNA) Nicotinamide Nucleotide Transhydrogenase-antisense RNA1 (NNT-AS1) role in colorectal cancer (CRC) hasn't been sufficiently inspected in relation to the Homo sapiens (hsa)-microRNA (miR)- 485-5p/ heat shock protein 90 (HSP90) axis, clinically. qRT-PCR was performed to detect lncRNA NNT-AS1 and hsa-miR-485-5p expression levels in 60 Egyptian patients' sera. HSP90 serum level was quantified using Enzyme-linked immunosorbent assay (ELISA). The relative expression level of the studied non-coding RNAs as well as the HSP90 ELISA concentration were correlated with patients clinicopathological characteristics and correlated to each other. The axis diagnostic utility in comparison with carbohydrate antigen 19-9 (CA19-9) and carcinoembryonic antigen (CEA) tumor markers (TMs) was studied by receiver operating characteristic (ROC) curve analysis. The relative lncRNA NNT-AS1 expression level fold change 56.7 (13.5-112) and HSP90 protein ELISA level 6.68 (5.14-8.77) (ng/mL) were elevated, while, for hsa-miR-485-5p 0.0474 (0.0236-0.135) expression fold change was repressed in CRC Egyptian patients' cohort sera, being compared to 28 apparently healthy control subjects. LncRNA NNT-AS1 specificity is 96.4% and a sensitivity of 91.7%, hsa-miR-485-5p showed 96.4% specificity, 90% sensitivity, and for HSP90 89.3%, 70% specificity and sensitivity, respectively. Those specificities and sensitivities were superior to the classical CRC TMs. A significant negative correlation was found between hsa-miR-485-5p with lncRNA NNT-AS1 (r = -0.933) expression fold change or with HSP90 protein blood level (r = -0.997), but, significant positive correlation was there between lncRNA NNT-AS1 and HSP90 (r = 0.927). LncRNA NNT-AS1/hsa-miR-485-5p/HSP90 axis could be a prospect for CRC development as well as diagnosis. Being correlated and related to CRC histologic grades 1-3, therefore, lncRNA NNT-AS1/hsa-miR-485-5p/HSP90 axis (not individually) expression approved clinically and in silico, could aid treatment precision.

A broadly conserved fungal alcohol oxidase (AOX) facilitates fungal invasion of plants

Alcohol oxidases (AOXs) are ecologically important enzymes that facilitate a number of plant-fungal interactions. Within Ascomycota they are primarily associated with methylotrophy, as a peroxisomal AOX catalysing the conversion of methanol to formaldehyde in methylotrophic yeast. In this study we demonstrate that AOX orthologues are phylogenetically conserved proteins that are common in the genomes of nonmethylotrophic, plant-associating fungi. Additionally, AOX orthologues are highly expressed during infection in a range of diverse pathosystems. To study the role of AOX in plant colonization, AOX knockout mutants were generated in the broad host range pathogen Sclerotinia sclerotiorum. Disease assays in soybean showed that these mutants had a significant virulence defect as evidenced by markedly reduced stem lesions and mortality rates. Chemical genomics suggested that SsAOX may function as an aromatic AOX, and growth assays demonstrated that ΔSsAOX is incapable of properly utilizing plant extract as a nutrient source. Profiling of known aromatic alcohols pointed towards the monolignol coniferyl alcohol (CA) as a possible substrate for SsAOX. As CA and other monolignols are ubiquitous among land plants, the presence of highly conserved AOX orthologues throughout Ascomycota implies that this is a broadly conserved protein used by ascomycete fungi during plant colonization.

Functional and Structural Characterization of OXA-935, a Novel OXA-10-Family β-Lactamase from Pseudomonas aeruginosa

Resistance to antipseudomonal penicillins and cephalosporins is often driven by the overproduction of the intrinsic β-lactamase AmpC. However, OXA-10-family β-lactamases are a rich source of resistance in Pseudomonas aeruginosa. OXA β-lactamases have a propensity for mutation that leads to extended spectrum cephalosporinase and carbapenemase activity. In this study, we identified isolates from a subclade of the multidrug-resistant (MDR) high risk P. aeruginosa clonal complex CC446 with a resistance to ceftazidime. A genomic analysis revealed that these isolates harbored a plasmid containing a novel allele of blaOXA-10, named blaOXA-935, which was predicted to produce an OXA-10 variant with two amino acid substitutions: an aspartic acid instead of a glycine at position 157 and a serine instead of a phenylalanine at position 153. The G157D mutation, present in OXA-14, is associated with the resistance of P. aeruginosa to ceftazidime. Compared to OXA-14, OXA-935 showed increased catalytic efficiency for ceftazidime. The deletion of blaOXA-935 restored the sensitivity to ceftazidime, and susceptibility profiling of P. aeruginosa laboratory strains expressing blaOXA-935 revealed that OXA-935 conferred ceftazidime resistance. To better understand the impacts of the variant amino acids, we determined the crystal structures of OXA-14 and OXA-935. Compared to OXA-14, the F153S mutation in OXA-935 conferred increased flexibility in the omega (Ω) loop. Amino acid changes that confer extended spectrum cephalosporinase activity to OXA-10-family β-lactamases are concerning, given the rising reliance on novel β-lactam/β-lactamase inhibitor combinations, such as ceftolozane-tazobactam and ceftazidime-avibactam, to treat MDR P. aeruginosa infections.

Nramp: Deprive and conquer?

Solute carriers 11 (Slc11) evolved from bacterial permease (MntH) to eukaryotic antibacterial defense (Nramp) while continuously mediating proton (H⁺)-dependent manganese (Mn²⁺) import. Also, Nramp horizontal gene transfer (HGT) toward bacteria led to mntH polyphyly. Prior demonstration that evolutionary rate-shifts distinguishing Slc11 from outgroup carriers dictate catalytic specificity suggested that resolving Slc11 family tree may provide a function-aware phylogenetic framework. Hence, MntH C (MC) subgroups resulted from HGTs of prototype Nramp (pNs) parologs while archetype Nramp (aNs) correlated with phagocytosis. PHI-Blast based taxonomic profiling confirmed MntH B phylogroup is confined to anaerobic bacteria vs. MntH A (MA)’s broad distribution; suggested niche-related spread of MC subgroups; established that MA-variant MH, which carries ‘eukaryotic signature’ marks, predominates in archaea. Slc11 phylogeny shows MH is sister to Nramp. Site-specific analysis of Slc11 charge network known to interact with the protonmotive force demonstrates sequential rate-shifts that recapitulate Slc11 evolution. 3D mapping of similarly coevolved sites across Slc11 hydrophobic core revealed successive targeting of discrete areas. The data imply that pN HGT could advantage recipient bacteria for H⁺-dependent Mn²⁺ acquisition and Alphafold 3D models suggest conformational divergence among MC subgroups. It is proposed that Slc11 originated as a bacterial stress resistance function allowing Mn²⁺-dependent persistence in conditions adverse for growth, and that archaeal MH could contribute to eukaryogenesis as a Mn²⁺ sequestering defense perhaps favoring intracellular growth-competent bacteria.

Large-Size Subunit Catalases Are Chimeric Proteins: A H2O2 Selecting Domain with Catalase Activity Fused to a Hsp31-Derived Domain Conferring Protein Stability and Chaperone Activity

Bacterial and fungal large-size subunit catalases (LSCs) are like small-size subunit catalases (SSCs) but have an additional C-terminal domain (CT). The catalytic domain is conserved at both primary sequence and structural levels and its amino acid composition is optimized to select H₂O₂ over water. The CT is structurally conserved, has an amino acid composition similar to very stable proteins, confers high stability to LSCs, and has independent molecular chaperone activity. While heat and denaturing agents increased Neurospora crassa catalase-1 (CAT-1) activity, a CAT-1 version lacking the CT (C63) was no longer activated by these agents. The addition of catalase-3 (CAT-3) CT to the CAT-1 or CAT-3 catalase domains prevented their heat denaturation in vitro. Protein structural alignments indicated CT similarity with members of the DJ-1/PfpI superfamily and the CT dimers present in LSCs constitute a new type of symmetric dimer within this superfamily. However, only the bacterial Hsp31 proteins show sequence similarity to the bacterial and fungal catalase mobile coil (MC) and are phylogenetically related to MC_CT sequences. LSCs might have originated by fusion of SSC and Hsp31 encoding genes during early bacterial diversification, conferring at the same time great stability and molecular chaperone activity to the novel catalases.

Structure-based pharmacophore mapping and virtual screening of natural products to identify polypharmacological inhibitor against c-MET/EGFR/VEGFR-2

Three receptor tyrosine kinases (RTKs), c-MET, EGFR, and VEGFR-2 have been identified as potential oncogenic targets involved in tumor development, metastasis, and invasion. Designing inhibitors that can simultaneously interact with multiple targets is a promising approach, therefore, inhibiting these three RTKs with a single chemical component might give an effective chemotherapeutic strategy for addressing the disease while limiting adverse effects. The in-silico methods have been developed to identify the polypharmacological inhibitors particularly for drug repurposing and multitarget drug design. Here, to find a viable inhibitor from natural source against these three RTKs, structure-based pharmacophore mapping and virtual screening of SN-II database were carried out. The filtered compound SN00020821, identified as Cedeodarin, from different computational approaches, demonstrated good interactions with all the three targets, c-MET/EGFR/VEGFR-2, with interaction energies of -42.35 kcal/mol, -49.32 kcal/mol and -44.83 kcal/mol, respectively. SN00020821displayed stable key interactions with critical amino acids of all the three receptors' kinase catalytic domains including "DFG motif" explored through the MD simulations. Furthermore, it also met the ADMET requirements and was determined to be drug-like as predicted from the Lipinski's rule of five and Veber's rule. Finally, SN00020821 provides a novel molecular scaffold that could be investigated further as a polypharmacological anticancer therapeutic candidate that targets the three RTKs.Communicated by Ramaswamy H. Sarma.

Everything Is Everywhere: Physiological Responses of the Mediterranean Sea and Eastern Pacific Ocean Epiphyte Cobetia Sp. to Varying Nutrient Concentration

Bacteria are essential in the maintenance and sustainment of marine environments (e.g., benthic systems), playing a key role in marine food webs and nutrient cycling. These microorganisms can live associated as epiphytic or endophytic populations with superior organisms with valuable ecological functions, e.g., seagrasses. Here, we isolated, identified, sequenced, and exposed two strains of the same species (i.e., identified as Cobetia sp.) from two different marine environments to different nutrient regimes using batch cultures: (1) Cobetia sp. UIB 001 from the endemic Mediterranean seagrass Posidonia oceanica and (2) Cobetia sp. 4B UA from the endemic Humboldt Current System (HCS) seagrass Heterozostera chilensis. From our physiological studies, both strains behaved as bacteria capable to cope with different nutrient and pH regimes, i.e., N, P, and Fe combined with different pH levels, both in long-term (12 days (d)) and short-term studies (4 d/96 h (h)). We showed that the isolated strains were sensitive to the N source (inorganic and organic) at low and high concentrations and low pH levels. Low availability of phosphorus (P) and Fe had a negative independent effect on growth, especially in the long-term studies. The strain UIB 001 showed a better adaptation to low nutrient concentrations, being a potential N₂-fixer, reaching higher growth rates (μ) than the HCS strain. P-acquisition mechanisms were deeply investigated at the enzymatic (i.e., alkaline phosphatase activity, APA) and structural level (e.g., alkaline phosphatase D, PhoD). Finally, these results were complemented with the study of biochemical markers, i.e., reactive oxygen species (ROS). In short, we present how ecological niches (i.e., MS and HCS) might determine, select, and modify the genomic and phenotypic features of the same bacterial species (i.e., Cobetia spp.) found in different marine environments, pointing to a direct correlation between adaptability and oligotrophy of seawater.

Neofunctionalization of Glycolytic Enzymes: An Evolutionary Route to Plant Parasitism in the Oomycete Phytophthora nicotianae

Oomycetes, of the genus Phytophthora, comprise of some of the most devastating plant pathogens. Parasitism of Phytophthora results from evolution from an autotrophic ancestor and adaptation to a wide range of environments, involving metabolic adaptation. Sequence mining showed that Phytophthora spp. display an unusual repertoire of glycolytic enzymes, made of multigene families and enzyme replacements. To investigate the impact of these gene duplications on the biology of Phytophthora and, eventually, identify novel functions associated to gene expansion, we focused our study on the first glycolytic step on P. nicotianae, a broad host range pathogen. We reveal that this step is committed by a set of three glucokinase types that differ by their structure, enzymatic properties, and evolutionary histories. In addition, they are expressed differentially during the P. nicotianae life cycle, including plant infection. Last, we show that there is a strong association between the expression of a glucokinase member in planta and extent of plant infection. Together, these results suggest that metabolic adaptation is a component of the processes underlying evolution of parasitism in Phytophthora, which may possibly involve the neofunctionalization of metabolic enzymes.

Mn2+ coordinates Cap-0-RNA to align substrates for efficient 2'-O-methyl transfer by SARS-CoV-2 nsp16

Virally encoded 2′-O-methyltransferases catalyze the last step in the capping of viral RNAs, which protects the RNAs from degradation and prevents them from triggering host defenses. Minasov et al. report structures of the SARS-CoV-2 methyltransferase, a heterodimeric complex of the enzyme nsp16 and its coactivator nsp10, in complex with a short, capped RNA (instead of the RNA cap analogs used to generate previous structures), the methyl donor SAM, and divalent metal cations. The metal ions and a four-residue insert of nsp16 were important for precisely aligning the RNA substrate in the active site for efficient catalysis. This insert is present in coronavirus but not in mammalian methyltransferases, suggesting this site as a potential target for the design of coronavirus-specific methyltransferase inhibitors.

Capping of viral messenger RNAs is essential for efficient translation, for virus replication, and for preventing detection by the host cell innate response system. The SARS-CoV-2 genome encodes the 2′-O-methyltransferase nsp16, which, when bound to the coactivator nsp10, uses S-adenosylmethionine (SAM) as a donor to transfer a methyl group to the first ribonucleotide of the mRNA in the final step of viral mRNA capping. Here, we provide biochemical and structural evidence that this reaction requires divalent cations, preferably Mn²⁺, and a coronavirus-specific four-residue insert. We determined the x-ray structures of the SARS-CoV-2 2′-O-methyltransferase (the nsp16-nsp10 heterodimer) in complex with its reaction substrates, products, and divalent metal cations. These structural snapshots revealed that metal ions and the insert stabilize interactions between the capped RNA and nsp16, resulting in the precise alignment of the ribonucleotides in the active site. Comparison of available structures of 2′-O-methyltransferases with capped RNAs from different organisms revealed that the four-residue insert unique to coronavirus nsp16 alters the backbone conformation of the capped RNA in the binding groove, thereby promoting catalysis. This insert is highly conserved across coronaviruses, and its absence in mammalian methyltransferases makes this region a promising site for structure-guided drug design of selective coronavirus inhibitors.

Semi-Automated Cell Panning for Efficient Isolation of FGFR3-Targeting Antibody

Phage display technology is a widely used practical tool for isolating binding molecules against the desired targets in phage libraries. In the case of targeting the membrane protein with its natural conformation, conventional bio-panning has limitations on the efficient screening of the functionally relevant antibodies. To enrich the single-chain variable fragment (scFv) pools for recognizing the natural conformation of the membrane targets, the conventional bio-panning and screening process was modified to include the semi-automated cell panning protocol. Using FGFR3-overexpressing patient-derived cancer cells, biotin-X-DHPE was introduced and coupled to Streptavidin-coated magnetic beads for use in the solution-phage bio-panning procedure. The resulting clones of scFv were compared to the diversity of the binding region, especially on CDR-H3. The clones enriched further by cell-based panning procedure possessed a similar binding site and the CDR-H3 loop structure. The resulting antibodies inhibited cell growth and induced target degradation. This process may be a useful tool for screening biologically related antibodies that recognize natural conformational structure on cell membrane protein. Furthermore, cell-based panning has the potential to further expand to a high-throughput screening (HTS) system and automation process.

Cancer-Associated circRNA-miRNA-mRNA Regulatory Networks: A Meta-Analysis

Recent advances in sequencing technologies and the discovery of non-coding RNAs (ncRNAs) have provided new insights in the molecular pathogenesis of cancers. Several studies have implicated the role of ncRNAs, including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and recently discovered circular RNAs (circRNAs) in tumorigenesis and metastasis. Unlike linear RNAs, circRNAs are highly stable and closed-loop RNA molecules. It has been established that circRNAs regulate gene expression by controlling the functions of miRNAs and RNA-binding protein (RBP) or by translating into proteins. The circRNA-miRNA-mRNA regulatory axis is associated with human diseases, such as cancers, Alzheimer's disease, and diabetes. In this study, we explored the interaction among circRNAs, miRNAs, and their target genes in various cancers using state-of-the-art bioinformatics tools. We identified differentially expressed circRNAs, miRNAs, and mRNAs on multiple cancers from publicly available data. Furthermore, we identified many crucial drivers and tumor suppressor genes in the circRNA-miRNA-mRNA regulatory axis in various cancers. Together, this study data provide a deeper understanding of the circRNA-miRNA-mRNA regulatory mechanisms in cancers.

Identification and phylogenetic analysis of RNA binding domain abundant in apicomplexans or RAP proteins

The RNA binding domain abundant in apicomplexans (RAP) is a protein domain identified in a diverse group of proteins, called RAP proteins, many of which have been shown to be involved in RNA binding. To understand the expansion and potential function of the RAP proteins, we conducted a hidden Markov model based screen among the proteomes of 54 eukaryotes, 17 bacteria and 12 archaea. We demonstrated that the domain is present in closely and distantly related organisms with particular expansions in Alveolata and Chlorophyta, and are not unique to Apicomplexa as previously believed. All RAP proteins identified can be decomposed into two parts. In the N-terminal region, the presence of variable helical repeats seems to participate in the specific targeting of diverse RNAs, while the RAP domain is mostly identified in the C-terminal region and is highly conserved across the different phylogenetic groups studied. Several conserved residues defining the signature motif could be crucial to ensure the function(s) of the RAP proteins. Modelling of RAP domains in apicomplexan parasites confirmed an ⍺/β structure of a restriction endonuclease-like fold. The phylogenetic trees generated from multiple alignment of RAP domains and full-length proteins from various distantly related eukaryotes indicated a complex evolutionary history of this family. We further discuss these results to assess the potential function of this protein family in apicomplexan parasites.

Thermostable D-amino acid decarboxylases derived from Thermotoga maritima diaminopimelate decarboxylase

Diaminopimelate decarboxylases (DAPDCs) are highly selective enzymes that catalyze the common final step in different lysine biosynthetic pathways, i.e. the conversion of meso-diaminopimelate (DAP) to L-lysine. We examined the modification of the substrate specificity of the thermostable decarboxylase from Thermotoga maritima with the aim to introduce activity with 2-aminopimelic acid (2-APA) since its decarboxylation leads to 6-aminocaproic acid (6-ACA), a building block for the synthesis of nylon-6. Structure-based mutagenesis of the distal carboxylate binding site resulted in a set of enzyme variants with new activities toward different D-amino acids. One of the mutants (E315T) had lost most of its activity toward DAP and primarily acted as a 2-APA decarboxylase. We next used computational modeling to explain the observed shift in catalytic activities of the mutants. The results suggest that predictive computational protocols can support the redesign of the catalytic properties of this class of decarboxylating PLP-dependent enzymes.

Two bi-functional cytochrome P450 CYP72 enzymes from olive (Olea europaea) catalyze the oxidative C-C bond cleavage in the biosynthesis of secoxy-iridoids - flavor and quality determinants in olive oil

Olive (Olea europaea) is an important crop in Europe, with high cultural, economic and nutritional significance. Olive oil flavor and quality depend on phenolic secoiridoids, but the biosynthetic pathway of these iridoids remains largely uncharacterized. We discovered two bifunctional cytochrome P450 enzymes, catalyzing the rare oxidative C-C bond cleavage of 7-epi-loganin to produce oleoside methyl ester (OeOMES) and secoxyloganin (OeSXS), both through a ketologanin intermediary. Although these enzymes are homologous to the previously reported Catharanthus roseus secologanin synthase (CrSLS), the substrate and product profiles differ. Biochemical assays provided mechanistic insights into the two-step OeOMES and CrSLS reactions. Model-guided mutations of OeOMES changed the product profile in a predictable manner, revealing insights into the molecular basis for this change in product specificity. Our results suggest that, in contrast to published hypotheses, in planta production of secoxy-iridoids is secologanin-independent. Notably, sequence data of cultivated and wild olives point to a relation between domestication and OeOMES expression. Thus, the discovery of this key biosynthetic gene suggests a link between domestication and secondary metabolism, and could potentially be used as a genetic marker to guide next-generation breeding programs.

ModFlex: Towards Function Focused Protein Modeling

There is a wide, and continuously widening, gap between the number of proteins known only by their amino acid sequence versus those structurally characterized by direct experiment. To close this gap, we mostly rely on homology-based inference and modeling to reason about the structures of the uncharacterized proteins by using structures of homologous proteins as templates. With the rapidly growing size of the Protein Data Bank, there are often multiple choices of templates, including multiple sets of coordinates from the same protein. The substantial conformational differences observed between different experimental structures of the same protein often reflect function related structural flexibility. Thus, depending on the questions being asked, using distant homologs, or coordinate sets with lower resolution but solved in the appropriate functional form, as templates may be more informative. The ModFlex server (https://modflex.org/) addresses this seldom mentioned gap in the standard homology modeling approach by providing the user with an interface with multiple options and tools to select the most relevant template and explore the range of structural diversity in the available templates. ModFlex is closely integrated with a range of other programs and servers developed in our group for the analysis and visualization of protein structural flexibility and divergence.

Dissecting protein domain variability in the core RNA interference machinery of five insect orders

RNA interference (RNAi)-mediated gene silencing can be used to control specific insect pest populations. Unfortunately, the variable efficiency in the knockdown levels of target genes has narrowed the applicability of this technology to a few species. Here, we examine the current state of knowledge regarding the miRNA (micro RNA) and siRNA (small interfering RNA) pathways in insects and investigate the structural variability at key protein domains of the RNAi machinery. Our goal was to correlate domain variability with mechanisms affecting the gene silencing efficiency. To this end, the protein domains of 168 insect species, encompassing the orders Coleoptera, Diptera, Hemiptera, Hymenoptera, and Lepidoptera, were analysed using our pipeline, which takes advantage of meticulous structure-based sequence alignments. We used phylogenetic inference and the evolutionary rate coefficient (K) to outline the variability across domain regions and surfaces. Our results show that four domains, namely dsrm, Helicase, PAZ and Ribonuclease III, are the main contributors of protein variability in the RNAi machinery across different insect orders. We discuss the potential roles of these domains in regulating RNAi-mediated gene silencing and the role of loop regions in fine-tuning RNAi efficiency. Additionally, we identified several order-specific singularities which indicate that lepidopterans have evolved differently from other insect orders, possibly due to constant coevolution with plants and viruses. In conclusion, our results highlight several variability hotspots that deserve further investigation in order to improve the application of RNAi technology in the control of insect pests.

The Cytotoxicity of RNase-Derived Peptides

Bacterial ribonuclease binase exhibits a cytotoxic effect on tumor cells possessing certain oncogenes. The aim of this study was to identify the structural parts of the binase molecule that exert cytotoxicity. Out of five designed peptides, the peptides representing the binase regions 21-50 and 74-94 have the highest cytotoxic potential toward human cervical HeLa and breast BT-20 and MCF-7 cancer cells. The peptides B21-50 and B74-94 were not able to enter human lung adenocarcinoma A549 cells, unlike BT-20 cells, explaining their failure to inhibit A549 cell proliferation. The peptide B74-94 shares similarities with epidermal growth factor (EGF), suggesting the peptide's specificity for EGF receptor overexpressed in BT-20 cells. Thus, the binase-derived peptides have the potential of being further developed as tumor-targeting peptides.

FATCAT 2.0: towards a better understanding of the structural diversity of proteins

FATCAT 2.0 server (http://fatcat.godziklab.org/), provides access to a flexible protein structure alignment algorithm developed in our group. In such an alignment, rotations and translations between elements in the structure are allowed to minimize the overall root mean square deviation (RMSD) between the compared structures. This allows to effectively compare protein structures even if they underwent structural rearrangements in different functional forms, different crystallization conditions or as a result of mutations. The major update for the server introduces a new graphical interface, much faster database searches and several new options for visualization of the structural differences between proteins

Comparison of metal-bound and unbound structures of aminopeptidase B proteins from Escherichia coli and Yersinia pestis

Protein degradation by aminopeptidases is involved in bacterial responses to stress. Escherichia coli produces two metal-dependent M17 family leucine aminopeptidases (LAPs), aminopeptidase A (PepA) and aminopeptidase B (PepB). Several structures have been solved for PepA as well as other bacterial M17 peptidases. Herein, we report the first structures of a PepB M17 peptidase. The E. coli PepB protein structure was determined at a resolution of 2.05 and 2.6 Å. One structure has both Zn²⁺ and Mn²⁺ , while the second structure has two Zn²⁺ ions bound to the active site. A 2.75 Å apo structure is also reported for PepB from Yersinia pestis. Both proteins form homohexamers, similar to the overall arrangement of PepA and other M17 peptidases. However, the divergent N-terminal domain in PepB is much larger resulting in a tertiary structure that is more expanded. Modeling of a dipeptide substrate into the C-terminal LAP domain reveals contacts that account for PepB to uniquely cleave after aspartate.

The crystal structure of nsp10-nsp16 heterodimer from SARS-CoV-2 in complex with S-adenosylmethionine

SARS-CoV-2 is a member of the coronaviridae family and is the etiological agent of the respiratory Coronavirus Disease 2019. The virus has spread rapidly around the world resulting in over two million cases and nearly 150,000 deaths as of April 17, 2020. Since no treatments or vaccines are available to treat COVID-19 and SARS-CoV-2, respiratory complications derived from the infections have overwhelmed healthcare systems around the world. This virus is related to SARS-CoV-1, the virus that caused the 2002-2004 outbreak of Severe Acute Respiratory Syndrome. In January 2020, the Center for Structural Genomics of Infectious Diseases implemented a structural genomics pipeline to solve the structures of proteins essential for coronavirus replication-transcription. Here we show the first structure of the SARS-CoV-2 nsp10-nsp16 2'-O-methyltransferase complex with S-adenosylmethionine at a resolution of 1.80 Å. This heterodimer complex is essential for capping viral mRNA transcripts for efficient translation and to evade immune surveillance.

Bioinformatic Analysis of the Flavin-Dependent Amine Oxidase Superfamily: Adaptations for Substrate Specificity and Catalytic Diversity

The flavin-dependent amine oxidase (FAO) superfamily consists of over 9000 nonredundant sequences represented in all domains of life. Of the thousands of members identified, only 214 have been functionally annotated to date, and 40 unique structures are represented in the Protein Data Bank. The few functionally characterized members share a catalytic mechanism involving the oxidation of an amine substrate through transfer of a hydride to the FAD cofactor, with differences observed in substrate specificities. Previous studies have focused on comparing a subset of superfamily members. Here, we present a comprehensive analysis of the FAO superfamily based on reaction mechanism and substrate recognition. Using a dataset of 9192 sequences, a sequence similarity network, and subsequently, a genome neighborhood network were constructed, organizing the superfamily into eight subgroups that accord with substrate type. Likewise, through phylogenetic analysis, the evolutionary relationship of subgroups was determined, delineating the divergence between enzymes based on organism, substrate, and mechanism. In addition, using sequences and atomic coordinates of 22 structures from the Protein Data Bank to perform sequence and structural alignments, active-site elements were identified, showing divergence from the canonical aromatic-cage residues to accommodate large substrates. These specificity determinants are held in a structural framework comprising a core domain catalyzing the oxidation of amines with an auxiliary domain for substrate recognition. Overall, analysis of the FAO superfamily reveals a modular fold with cofactor and substrate-binding domains allowing for diversity of recognition via insertion/deletions. This flexibility allows facile evolution of new activities, as shown by reinvention of function between subfamilies.

In silico identification and evaluation of Bacillus subtilis cold shock protein B (cspB)-like plant RNA chaperones

Cold shock domain (CSD) proteins with nucleic acid binding properties are well conserved from bacteria to higher organisms. In bacteria, the cold shock proteins (CSPs) are single domain RNA chaperones, whereas in animals and plants, CSDs are accompanied by additional domains with roles in transcription regulation. Bacterial CSPs (Escherischia coli-cspA and Bacilus subtilis-cspB) have successfully imparted drought tolerance in transgenic plants; however, these cannot be deployed in food crops due to their low public acceptance of transgenics with bacterial genes. Therefore, this study aimed to identify CSPB-like proteins from plants that can be used for developing drought tolerant transgenic crops. Twelve single domain plant CSPs presenting >40% sequence identity with CSPB were identified. All 12 plant CSPs were modeled by homology modeling and refined by molecular dynamics simulation for 10 ns. Selected plant CSPs and CSPB exhibited high structural similarity (Tm-score: 0.63-0.86). Structure based phylogenetic analysis revealed that Triticum aestivum-csp1 and Aegilops tauschii-cspE are structurally closer to CSPB compared to their orthologs and paralogs. Molecular docking with three RNA molecules (5U, UC3U, and C2UC) indicates that Ricinus communis-csd1 and T. aestivum-csp1 have a binding pattern and docking scores similar to those of CSPB. Furthermore, MD simulations for 20 ns and analysis of RMSD, RMSF, Rg as well as the number of hydrogen bonds in all the three complexes revealed that plant CSP-RNA complexes behave in a similar manner to that of the CSPB-RNA complex, making them highly potential candidate genes for developing drought tolerance in transgenic plants. Communicated by Ramaswamy H. Sarma.

Investigation of the Relationship between the S1 Domain and Its Molecular Functions Derived from Studies of the Tertiary Structure

S1 domain, a structural variant of one of the "oldest" OB-folds (oligonucleotide/oligosaccharide-binding fold), is widespread in various proteins in three domains of life: Bacteria, Eukaryotes, and Archaea. In this study, it was shown that S1 domains of bacterial, eukaryotic, and archaeal proteins have a low percentage of identity, which indicates the uniqueness of the scaffold and is associated with protein functions. Assessment of the predisposition of tertiary flexibility of S1 domains using computational and statistical tools showed similar structural features and revealed functional flexible regions that are potentially involved in the interaction of natural binding partners. In addition, we analyzed the relative number and distribution of S1 domains in all domains of life and established specific features based on sequences and structures associated with molecular functions. The results correlate with the presence of repeats of the S1 domain in proteins containing the S1 domain in the range from one (bacterial and archaeal) to 15 (eukaryotic) and, apparently, are associated with the need for individual proteins to increase the affinity and specificity of protein binding to ligands.

The Role of Iron in the P-Acquisition Mechanisms of the Unicellular N2-Fixing Cyanobacteria Halothece sp., Found in Association With the Mediterranean Seagrass Posidonia oceanica

Posidonia oceanica, an endemic seagrass of the Mediterranean Sea harbors a high diversity of N₂-fixing prokaryotes. One of these is Halothece sp., a unicellular N₂-fixing cyanobacteria detected through nifH analysis from the epiphytes of P. oceanica. The most related strain in culture is Halothece sp. PCC 7418 and this was used as the test organism in this study. In the Mediterranean Sea, phosphorus (P) and iron (Fe) can be the major limiting nutrients for N₂ fixation. However, information about the mechanisms of P-acquisition and the role of metals (i.e., Fe) in these processes for N₂-fixing bacteria is scarce. From our genomic analyses of the test organism and other phylogenetically related N₂-fixing strains, Halothece sp. PCC 7418 is one of the strains with the greatest number of gene copies (eight copies) of alkaline phosphatases (APases). Our structural analysis of PhoD (alkaline phosphatase type D) and PhoU (phosphate acquisition regulator) of Halothece sp. PCC 7418 showed the connection among metals (Ca²⁺ and Fe³⁺), and the P-acquisition mechanisms. Here, we measured the rates of alkaline phosphatase activity (APA) through MUF-P hydrolysis under different combinations of concentrations of inorganic P (PO43−) and Fe in experiments under N₂-fixing (low NO3− availability) and non-N₂ fixing (high NO3− availability) conditions. Our results showed that APA rates were enhanced by the increase in Fe availability under low levels of PO43−, especially under N₂-fixing conditions. Moreover, the increased PO43−-uptake was reflected in the increased of the P-cellular content of the cells under N₂ fixation conditions. We also found a positive significant relationship between cellular P and cellular Fe content of the cells (r² = 0.71, p < 0.05). Our results also indicated that Fe-uptake in Halothece sp. PCC 7418 was P and Fe-dependent. This study gives first insights of P-acquisition mechanisms in the N₂-fixing cyanobacteria (Halothece sp.) found in P. oceanica and highlights the role of Fe in these processes.

mTM-align: a server for fast protein structure database search and multiple protein structure alignment

With the rapid increase of the number of protein structures in the Protein Data Bank, it becomes urgent to develop algorithms for efficient protein structure comparisons. In this article, we present the mTM-align server, which consists of two closely related modules: one for structure database search and the other for multiple structure alignment. The database search is speeded up based on a heuristic algorithm and a hierarchical organization of the structures in the database. The multiple structure alignment is performed using the recently developed algorithm mTM-align. Benchmark tests demonstrate that our algorithms outperform other peering methods for both modules, in terms of speed and accuracy. One of the unique features for the server is the interplay between database search and multiple structure alignment. The server provides service not only for performing fast database search, but also for making accurate multiple structure alignment with the structures found by the search. For the database search, it takes about 2-5 min for a structure of a medium size (∼300 residues). For the multiple structure alignment, it takes a few seconds for ∼10 structures of medium sizes. The server is freely available at: http://yanglab.nankai.edu.cn/mTM-align/.

Automated shape-based clustering of 3D immunoglobulin protein structures in chronic lymphocytic leukemia

Background

Although the etiology of chronic lymphocytic leukemia (CLL), the most common type of adult leukemia, is still unclear, strong evidence implicates antigen involvement in disease ontogeny and evolution. Primary and 3D structure analysis has been utilised in order to discover indications of antigenic pressure. The latter has been mostly based on the 3D models of the clonotypic B cell receptor immunoglobulin (BcR IG) amino acid sequences. Therefore, their accuracy is directly dependent on the quality of the model construction algorithms and the specific methods used to compare the ensuing models. Thus far, reliable and robust methods that can group the IG 3D models based on their structural characteristics are missing.

Results

Here we propose a novel method for clustering a set of proteins based on their 3D structure focusing on 3D structures of BcR IG from a large series of patients with CLL. The method combines techniques from the areas of bioinformatics, 3D object recognition and machine learning. The clustering procedure is based on the extraction of 3D descriptors, encoding various properties of the local and global geometrical structure of the proteins. The descriptors are extracted from aligned pairs of proteins. A combination of individual 3D descriptors is also used as an additional method. The comparison of the automatically generated clusters to manual annotation by experts shows an increased accuracy when using the 3D descriptors compared to plain bioinformatics-based comparison. The accuracy is increased even more when using the combination of 3D descriptors.

Conclusions

The experimental results verify that the use of 3D descriptors commonly used for 3D object recognition can be effectively applied to distinguishing structural differences of proteins. The proposed approach can be applied to provide hints for the existence of structural groups in a large set of unannotated BcR IG protein files in both CLL and, by logical extension, other contexts where it is relevant to characterize BcR IG structural similarity. The method does not present any limitations in application and can be extended to other types of proteins.

First insights into a type II toxin-antitoxin system from the clinical isolate Mycobacterium sp. MHSD3, similar to epsilon/zeta systems

A putative type II toxin-antitoxin (TA) system was found in the clinical isolate Mycobacterium sp. MHSD3, a strain closely related to Mycobacterium chelonae. Further analyses of the protein sequences of the two genes revealed the presence of domains related to a TA system. BLAST analyses indicated the presence of closely related proteins in the genomes of other recently published M. chelonae strains. The functionality of both elements of the TA system was demonstrated when expressed in Escherichia coli cells, and the predicted structure of the toxin is very similar to those of well-known zeta-toxins, leading to the definition of a type II TA system similar to epsilon/zeta TA systems in strains that are closely related to M. chelonae.

Computational Analysis of Gynura bicolor Bioactive Compounds as Dipeptidyl Peptidase-IV Inhibitor

The inhibition of dipeptidyl peptidase-IV (DPPIV) is a popular route for the treatment of type-2 diabetes. Commercially available gliptin-based drugs such as sitagliptin, anagliptin, linagliptin, saxagliptin, and alogliptin were specifically developed as DPPIV inhibitors for diabetic patients. The use of Gynura bicolor in treating diabetes had been reported in various in vitro experiments. However, an understanding of the inhibitory actions of G. bicolor bioactive compounds on DPPIV is still lacking and this may provide crucial information for the development of more potent and natural sources of DPPIV inhibitors. Evaluation of G. bicolor bioactive compounds for potent DPPIV inhibitors was computationally conducted using Lead IT and iGEMDOCK software, and the best free-binding energy scores for G. bicolor bioactive compounds were evaluated in comparison with the commercial DPPIV inhibitors, sitagliptin, anagliptin, linagliptin, saxagliptin, and alogliptin. Drug-likeness and absorption, distribution, metabolism, and excretion (ADME) analysis were also performed. Based on molecular docking analysis, four of the identified bioactive compounds in G. bicolor, 3-caffeoylquinic acid, 5-O-caffeoylquinic acid, 3,4-dicaffeoylquinic acid, and trans-5-p-coumaroylquinic acid, resulted in lower free-binding energy scores when compared with two of the commercially available gliptin inhibitors. The results revealed that bioactive compounds in G. bicolor are potential natural inhibitors of DPPIV.

XSuLT: a web server for structural annotation and representation of sequence-structure alignments

The web server XSuLT, an enhanced version of the protein alignment annotation program JoY, formats a submitted multiple-sequence alignment using three-dimensional (3D) structural information in order to assist in the comparative analysis of protein evolution and in the optimization of alignments for comparative modelling and construct design. In addition to the features analysed by JoY, which include secondary structure, solvent accessibility and sidechain hydrogen bonds, XSuLT annotates each amino acid residue with residue depth, chain and ligand interactions, inter-residue contacts, sequence entropy, root mean square deviation and secondary structure and disorder prediction. It is also now integrated with built-in 3D visualization which interacts with the formatted alignment to facilitate inspection and understanding. Results can be downloaded as stand-alone HTML for the formatted alignment and as XML with the underlying annotation data. XSuLT is freely available at http://structure.bioc.cam.ac.uk/xsult/.

Similarity Search of Flexible 3D Molecules Combining Local and Global Shape Descriptors

In this paper, a framework for shape-based similarity search of 3D molecular structures is presented. The proposed framework exploits simultaneously the discriminative capabilities of a global, a local, and a hybrid local-global shape feature to produce a geometric descriptor that achieves higher retrieval accuracy than each feature does separately. Global and hybrid features are extracted using pairwise computations of diffusion distances between the points of the molecular surface, while the local feature is based on accumulating pairwise relations among oriented surface points into local histograms. The local features are integrated into a global descriptor vector using the bag-of-features approach. Due to the intrinsic property of its constituting shape features to be invariant to articulations of the 3D objects, the framework is appropriate for similarity search of flexible 3D molecules, while at the same time it is also accurate in retrieving rigid 3D molecules. The proposed framework is evaluated in flexible and rigid shape matching of 3D protein structures as well as in shape-based virtual screening of large ligand databases with quite promising results.

Structural basis for DNA recognition by the transcription regulator MetR

MetR, a LysR-type transcriptional regulator (LTTR), has been extensively studied owing to its role in the control of methionine biosynthesis in proteobacteria. A MetR homodimer binds to a 24-base-pair operator region of the met genes and specifically recognizes the interrupted palindromic sequence 5'-TGAA-N5-TTCA-3'. Mechanistic details underlying the interaction of MetR with its target DNA at the molecular level remain unknown. In this work, the crystal structure of the DNA-binding domain (DBD) of MetR was determined at 2.16 Å resolution. MetR-DBD adopts a winged-helix-turn-helix (wHTH) motif and shares significant fold similarity with the DBD of the LTTR protein BenM. Furthermore, a data-driven macromolecular-docking strategy was used to model the structure of MetR-DBD bound to DNA, which revealed that a bent conformation of DNA is required for the recognition helix α3 and the wing loop of the wHTH motif to interact with the major and minor grooves, respectively. Comparison of the MetR-DBD-DNA complex with the crystal structures of other LTTR-DBD-DNA complexes revealed residues that may confer operator-sequence binding specificity for MetR. Taken together, the results show that MetR-DBD uses a combination of direct base-specific interactions and indirect shape recognition of the promoter to regulate the transcription of met genes.

Calculating and scoring high quality multiple flexible protein structure alignments

MOTIVATION

Calculating multiple protein structure alignments (MSAs) is important for understanding functional and evolutionary relationships between protein families, and for modeling protein structures by homology. While incorporating backbone flexibility promises to circumvent many of the limitations of rigid MSA algorithms, very few flexible MSA algorithms exist today. This article describes several novel improvements to the Kpax algorithm which allow high quality flexible MSAs to be calculated. This article also introduces a new Gaussian-based MSA quality measure called 'M-score', which circumvents the pitfalls of RMSD-based quality measures.

RESULTS

As well as calculating flexible MSAs, the new version of Kpax can also score MSAs from other aligners and from previously aligned reference datasets. Results are presented for a large-scale evaluation of the Homstrad, SABmark and SISY benchmark sets using Kpax and Matt as examples of state-of-the-art flexible aligners and 3DCOMB as an example of a state-of-the-art rigid aligner. These results demonstrate the utility of the M-score as a measure of MSA quality and show that high quality MSAs may be achieved when structural flexibility is properly taken into account.

AVAILABILITY AND IMPLEMENTATION

Kpax 5.0 may be downloaded for academic use at http://kpax.loria.fr/

CONTACT

dave.ritchie@inria.fr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Crystal structure of HINT from Helicobacter pylori

Proteins belonging to the histidine triad (HIT) superfamily bind nucleotides and use the histidine triad motif to carry out dinucleotidyl hydrolase, nucleotidyltransferase and phosphoramidite hydrolase activities. Five different branches of this superfamily are known to exist. Defects in these proteins in humans are linked to many diseases such as ataxia, diseases of RNA metabolism and cell-cycle regulation, and various types of cancer. The histidine triad nucleotide protein (HINT) is nearly identical to proteins that have been classified as protein kinase C-interacting proteins (PKCIs), which also have the ability to bind and inhibit protein kinase C. The structure of HINT, which exists as a homodimer, is highly conserved from humans to bacteria and shares homology with the product of fragile histidine triad protein (FHit), a tumour suppressor gene of this superfamily. Here, the structure of HINT from Helicobacter pylori (HpHINT) in complex with AMP is reported at a resolution of 3 Å. The final model has R and Rfree values of 26 and 28%, respectively, with good electron density. Structural comparison with previously reported homologues and phylogenetic analysis shows H. pylori HINT to be the smallest among them, and suggests that it branched out separately during the course of evolution. Overall, this structure has contributed to a better understanding of this protein across the animal kingdom.

PDBFlex: exploring flexibility in protein structures

The PDBFlex database, available freely and with no login requirements at http://pdbflex.org, provides information on flexibility of protein structures as revealed by the analysis of variations between depositions of different structural models of the same protein in the Protein Data Bank (PDB). PDBFlex collects information on all instances of such depositions, identifying them by a 95% sequence identity threshold, performs analysis of their structural differences and clusters them according to their structural similarities for easy analysis. The PDBFlex contains tools and viewers enabling in-depth examination of structural variability including: 2D-scaling visualization of RMSD distances between structures of the same protein, graphs of average local RMSD in the aligned structures of protein chains, graphical presentation of differences in secondary structure and observed structural disorder (unresolved residues), difference distance maps between all sets of coordinates and 3D views of individual structures and simulated transitions between different conformations, the latter displayed using JSMol visualization software.