SuperLigands – a database of ligand structures derived from the Protein Data Bank

Background  

Currently, the PDB contains approximately 29,000 protein structures comprising over 70,000 experimentally determined three-dimensional structures of over 5,000 different low molecular weight compounds. Information about these PDB ligands can be very helpful in the field of molecular modelling and prediction, particularly for the prediction of protein binding sites and function.     

Crystal Structures of an Oligopeptide-Binding Protein from the Biosynthetic Pathway of the β-Lactamase Inhibitor Clavulanic Acid

Clavulanic acid (CA) is a clinically important β-lactamase inhibitor that is produced by fermentation of Streptomyces clavuligerus. The CA biosynthesis pathway starts from arginine and glyceraldehyde-3-phosphate and proceeds via (3S,5S)-clavaminic acid, which is converted to (3R,5R)-clavaldehyde, the immediate precursor of (3R,5R)-CA. Open reading frames 7 (orf7) and 15 (orf15) of the CA biosynthesis cluster encode oligopeptide-binding proteins (OppA1 and OppA2), which are essential for CA biosynthesis. OppA1/2 are proposed to be involved in the binding and/or transport of peptides across the S. clavuligerus cell membrane. Peptide binding assays reveal that recombinant OppA1 and OppA2 bind di-/tripeptides containing arginine and certain nonapeptides including bradykinin. Crystal structures of OppA2 in its apo form and in complex with arginine or bradykinin were solved to 1.45, 1.7, and 1.7 Å resolution, respectively. The overall fold of OppA2 consists of two lobes with a deep cavity in the center, as observed for other oligopeptide-binding proteins. The large cavity creates a peptide/arginine binding cleft. The crystal structures of OppA2 in complex with arginine or bradykinin reveal that the C-terminal arginine of bradykinin binds similarly to arginine. The results are discussed in terms of the possible roles of OppA1/2 in CA biosynthesis.     

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of Cα Atoms

Background

The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the Cα atoms which are modeled with modest accuracy.

Methodology/Principal Findings

In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of Cα atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of Cα. We extend the method further to recognize potential protein-protein interface residues.

Conclusion/ Significance

Our approach to identify buried and exposed residues solely from the positions of Cα atoms resulted in an accuracy of 84%, sensitivity of 83–89% and specificity of 67–94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70–96% and specificity of 58–94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from Cα positions and all-atom models suggested that, recognition of interfacial residues using Cα atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only Cα positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.     

Identification of the Site of Inhibition of Oncogenic ras–p21-Induced Signal Transduction by a Peptide from a ras Effector Domain

We have previously found that a peptide corresponding to residues 35–47 of the ras-p21 protein, from its switch 1 effector domain region, strongly inhibits oocyte maturation induced by oncogenic p21, but not by insulin-activated cellular wild-type p21. Another ras–p21 peptide corresponding to residues 96–110 that blocks ras–jun and jun kinase (JNK) interactions exhibits a similar pattern of inhibition. We have also found that c-raf strongly induces oocyte maturation and that dominant negative c-raf strongly blocks oncogenic p21-induced oocyte maturation. We now find that the p21 35–47, but not the 96–110, peptide completely blocks c-raf-induced maturation. This finding suggests that the 35–47 peptide blocks oncogenic ras at the level of raf; that activated normal and oncogenic ras–p21 have differing requirements for raf-dependent signaling; and that the two oncogenic-ras-selective inhibitory peptides, 35–47 and 96–110, act at two different critical downstream sites, the former at raf, the latter at JNK/jun, both of which are required for oncogenic ras-p21 signaling.     

Protein disulfide isomerases – a way to tackle malaria

Protein disulfide isomerases (PDIs) ensure that specific substrate proteins are correctly folded. PDI activity plays an essential role in malaria transmission. Here we provide an overview of the role of PDIs in malaria-causing Plasmodium parasites and outline why PDI inhibition could be a novel way to treat malaria and prevent transmission.     

G Protein �¦� Subunits Regulate Cell Adhesion through Rap1a and Its Effector Radil

The activation of several G protein-coupled receptors is known to regulate the adhesive properties of cells in different contexts. Here, we reveal that Gβγ subunits of heterotrimeric G proteins regulate cell-matrix adhesiveness by activating Rap1a-dependent inside-out signals and integrin activation. We show that Gβγ subunits enter in a protein complex with activated Rap1a and its effector Radil and establish that this complex is required downstream of receptor stimulation for the activation of integrins and the positive modulation of cell-matrix adhesiveness. Moreover, we demonstrate that Gβγ and activated Rap1a promote the translocation of Radil to the plasma membrane at sites of cell-matrix contacts. These results add to the molecular understanding of how G protein-coupled receptors impinge on cell adhesion and suggest that the Gβγ·Rap1·Radil complex plays important roles in this process.     

Prediction of Protein–Protein Interaction Sites in Sequences and 3D Structures by Random Forests

Identifying interaction sites in proteins provides important clues to the function of a protein and is becoming increasingly relevant in topics such as systems biology and drug discovery. Although there are numerous papers on the prediction of interaction sites using information derived from structure, there are only a few case reports on the prediction of interaction residues based solely on protein sequence. Here, a sliding window approach is combined with the Random Forests method to predict protein interaction sites using (i) a combination of sequence- and structure-derived parameters and (ii) sequence information alone. For sequence-based prediction we achieved a precision of 84% with a 26% recall and an F-measure of 40%. When combined with structural information, the prediction performance increases to a precision of 76% and a recall of 38% with an F-measure of 51%. We also present an attempt to rationalize the sliding window size and demonstrate that a nine-residue window is the most suitable for predictor construction. Finally, we demonstrate the applicability of our prediction methods by modeling the Ras–Raf complex using predicted interaction sites as target binding interfaces. Our results suggest that it is possible to predict protein interaction sites with quite a high accuracy using only sequence information.     

Selection of Protein–Protein Interactions of Desired Affinities with a Bandpass Circuit

We have developed a genetic circuit in Escherichia coli that can be used to select for protein–protein interactions of different strengths by changing antibiotic concentrations in the media. The genetic circuit links protein–protein interaction strength to β-lactamase activity while simultaneously imposing tuneable positive and negative selection pressure for β-lactamase activity. Cells only survive if they express interacting proteins with affinities that fall within set high- and low-pass thresholds; i.e. the circuit therefore acts as a bandpass filter for protein–protein interactions. We show that the circuit can be used to recover protein–protein interactions of desired affinity from a mixed population with a range of affinities. The circuit can also be used to select for inhibitors of protein–protein interactions of defined strength.     

Comparison of the Predicted Structures of Loops in the ras–SOS Protein Bound to a Single ras-p21 Protein with the Crystallographically Determined Structures in SOS Bound to Two ras-p21 Proteins

We have previously computed the structures of three loops, residues 591–596, 654–675 and 742–751, in the ras-p21 protein-binding domain (residues 568–1044) of the guanine nucleotide-exchange-promoting SOS protein that were crystallographically undefined when one molecule of ras-p21 (unbound to nucleotide) binds to SOS. Based on our computational results, we synthesized three peptides corresponding to sequences of each of these three loops and found that all three peptides strongly inhibit ras-p21 signaling. More recently, a new crystal structure of SOS has been determined in which this protein binds to two molecules of ras-p21, one unbound to GTP and one bound to GTP. In this structure, the 654–675 loop and residues 742–743 and 750–751 are now crystallographically defined. We have superimposed our energy-minimized structure of the ras-binding domain of SOS bound to one molecule of ras-p21 on the X-ray structure for SOS bound to two molecules of ras-p21. We find that, while the two structures are superimposable, there are large deviations of the residues 673 and 676 and 741 and 752, flanking the two loop segments. This suggests that the binding of the extra ras-p21 molecule, which is far from each of the three loops, induces conformational changes in these domains and further supports their role in signal transduction. In spite of these differences, we have superimposed our computed structures for the loop residues on those from the more recent X-ray structure. Our structure for the 654–675 segment is an anti-parallel beta-sheet with a reverse turn at residues 663–665; in the X-ray structure residues 655–662 adopt an alpha-helical conformation; on the other hand, our computed structure for residues 663–675 superimpose on the X-ray structure for these residues. We further find that our computed structures for residues 742–743 and 750–751 are superimposable on the X-ray structure for these residues.     

Protein dynamics of a β-sheet protein

Rhodnius prolixus Nitrophorin 4 (abbreviated NP4) is an almost pure β-sheet heme protein. Its dynamics is investigated by X-ray structure determination at eight different temperatures from 122 to 304 K and by means of Mössbauer spectroscopy. A comparison of this β-sheet protein with the pure α-helical protein myoglobin (abbreviated Mbmet) is performed. The mean square displacement derived from the Mössbauer spectra increases linearly with temperature below a characteristic temperature T _c. It is about 10 K larger than that of myoglobin. Above T _c the mean square displacements increase dramatically. The Mössbauer spectra are analyzed by a two state model. The increased mean square displacements are caused by very slow motions occurring on a time scale faster than 140 ns. With respect to these motions NP4 shows the same protein specific modes as Mbmet. There is, however, a difference in the fast vibration regime. The B values found in the X-ray structures vary linearly over the entire temperature range. The mean square displacements in NP4 increase with slopes which are 60% larger than those observed for Mbmet. This indicates that nitrophorin has a larger structural distribution which makes it more flexible than myoglobin.     

Characterization of a Single-Stranded DNA-Binding-Like Protein from Nanoarchaeum equitans—A Nucleic Acid Binding Protein with Broad Substrate Specificity

BackgroundSSB (single-stranded DNA-binding) proteins play an essential role in all living cells and viruses, as they are involved in processes connected with ssDNA metabolism. There has recently been an increasing interest in SSBs, since they can be applied in molecular biology techniques and analytical methods. Nanoarchaeum equitans, the only known representative of Archaea phylum Nanoarchaeota, is a hyperthermophilic, nanosized, obligatory parasite/symbiont of Ignicoccus hospitalis.ResultsThis paper reports on the ssb-like gene cloning, gene expression and characterization of a novel nucleic acid binding protein from Nanoarchaeum equitans archaeon (NeqSSB-like protein). This protein consists of 243 amino acid residues and one OB fold per monomer. It is biologically active as a monomer like as SSBs from some viruses. The NeqSSB-like protein displays a low sequence similarity to the Escherichia coli SSB, namely 10% identity and 29% similarity, and is the most similar to the Sulfolobus solfataricus SSB (14% identity and 32% similarity). The NeqSSB-like protein binds to ssDNA, although it can also bind mRNA and, surprisingly, various dsDNA forms, with no structure-dependent preferences as evidenced by gel mobility shift assays. The size of the ssDNA binding site, which was estimated using fluorescence spectroscopy, is 7±1 nt. No salt-dependent binding mode transition was observed. NeqSSB-like protein probably utilizes a different model for ssDNA binding than the SSB proteins studied so far. This protein is highly thermostable; the half-life of the ssDNA binding activity is 5 min at 100°C and melting temperature (T_m) is 100.2°C as shown by differential scanning calorimetry (DSC) analysis.ConclusionNeqSSB-like protein is a novel highly thermostable protein which possesses a unique broad substrate specificity and is able to bind all types of nucleic acids.     

The AlphaFold Database of Protein Structures: A Biologist’s Guide

AlphaFold, the deep learning algorithm developed by DeepMind, recently released the three-dimensional models of the whole human proteome to the scientific community. Here we discuss the advantages, limitations and the still unsolved challenges of the AlphaFold models from the perspective of a biologist, who may not be an expert in structural biology.     

Microalgae as a source of bioactive molecules – experience from cyanophyte research

Detrimental and beneficial properties of prokaryoticand eukaryotic microalgae may be qualifications givingthem posibillities for their biotechnologicalutilization. Phycotoxins and related products mayserve as material for useful drugs. Cyanotoxins arean example. Bioassay screening methods are used todetect specific biological activity (e.g. antibioticeffects). A study demonstrating antibacterialproperties among cyanophytes showed that these effectswere caused by substances produced by both toxigenicand non-toxigenic strains. Differences inbiosynthetic activity are influenced by the type ofclones and by the culture conditions (genetic andbiochemical variation). Experience of recentscreening (1997–1998) of compounds from cyanophytes –including 25 species from 14 genera – display thediscovery of more than 30 different antibioticsubstances. The further isolation of newpharmaceuticals from natural resources is becomingdifficult, if not uncommon species of producerorganisms are taken into consideration.     

An Experimental Tool to Estimate the Probability of a Nucleotide Presence in the Crystal Structures of the Nucleotide–Protein Complexes

A correlation between the ligand–protein affinity and the identification of the ligand in the experimental electron density maps obtained by X-ray crystallography has been tested for a number of RNA-binding proteins. Bacterial translation regulators ProQ, TRAP, Rop, and Hfq together with their archaeal homologues SmAP have been used. The equilibrium dissociation constants for the N-methyl-anthraniloyl-labelled adenosine and guanosine monophosphates titrated by the proteins have been determined by the fluorescent anisotropy measurements. The estimated stability of the nucleotide–protein complexes has been matched with a presence of the nucleotides in the structures of the proposed nucleotide–protein complexes. It has been shown that the ribonucleotides can be definitely identified in the experimental electron density maps at equilibrium dissociation constant <10 μM. At K_D of 20–40 μM, long incubation of the protein crystals in the nucleotide solution is required to obtain the structures of the complexes. The complexes with K_D value higher than 50 μM are not stable enough to survive in crystallization conditions.     

The glomerulus – a view from the outside – the podocyte

In the past decade, podocyte research has been greatly aided by the development of powerful new molecular, cellular and animal tools, leading to elucidation of an increasing number of proteins involved in podocyte function and identification of mutated genes in hereditary glomerulopathies. Accumulating evidence indicates that podocyte disorders may not only underlie these hereditary glomerulopathies but also play crucial role in a broad spectrum of acquired glomerular diseases. Genetic susceptibility, environmental influence and systemic responses are all involved in the mediation of the pathogenesis of podocytopathies. Injured podocytes may predisopose to further injury of other podocytes and other adjacent/distant renal cells in a vicious cycle, leading to inexorable progression of glomerular injury. The classic view is that podocytes have a limited ability to proliferate in the normal mature kidney. However, recent research in rodents has provided suggestive evidence for podocyte regeneration resulting from differentiation of progenitor cells within Bowman's capsule.     

Environmental Impact of the Production of Mealworms as a Protein Source for Humans – A Life Cycle Assessment

The demand for animal protein is expected to rise by 70–80% between 2012 and 2050, while the current animal production sector already causes major environmental degradation. Edible insects are suggested as a more sustainable source of animal protein. However, few experimental data regarding environmental impact of insect production are available. Therefore, a lifecycle assessment for mealworm production was conducted, in which greenhouse gas production, energy use and land use were quantified and compared to conventional sources of animal protein. Production of one kg of edible protein from milk, chicken, pork or beef result in higher greenhouse gas emissions, require similar amounts of energy and require much more land. This study demonstrates that mealworms should be considered a more sustainable source of edible protein.     

NMR Structures of a Nonapeptide from DNA Binding Domain of Human Polymerase-α Determined by Iterative Complete-Relaxation-Matrix Approach

Abstract

Nuclear magnetic resonance structures of a nonapeptide, ERFKCPCPT, selected from the DNA binding domain of human polymerase-a, were determined by complete relaxation matrix analysis of transverse NOE data. The structures exhibit a type III turn with residues KCPC, and the remaining residues exhibit non-ordered structures. The turn was confirmed by α, N (i,i+3) connectivity, a low temperature coefficient of NH chemical shift (?3.1 × 1O^?3) of the fourth residue, ³J_NHα coupling constants, and characteristic CD peaks at 228 and 200 nm. Furthermore, ø and ψ dihedral angles for the i + 1, and i + 2 residues of the tum are found to be?80 and?41 and?60 and?40 degrees. The first proline residue is trans- while the second exists in both eis- and trans- configurations, with trans- being more than 80% populated. The trans-configuration was established from C5α-P6α correlation and ø and ψ angles of the proline. The five-membered proline ring is in DOWN puckered (C-β-exo/C-γ-endo) conformation. The structure of the peptide reveals that the two cysteine thiols are?5 A° apart and appropriately positioned to covalently bind cis-diamminedichloroplatinum(II), a widely used anti-cancer drug.     

Design of a Zinc-Finger Hydrolase with a Synthetic αββ Protein

Recent advances in protein design have opened avenues for the creation of artificial enzymes needed for biotechnological and pharmaceutical applications. However, designing efficient enzymes remains an unrealized ambition, as the design must incorporate a catalytic apparatus specific for the desired reaction. Here we present a de novo design approach to evolve a minimal carbonic anhydrase mimic. We followed a step-by-step design of first folding the main chain followed by sequence variation for substrate binding and catalysis. To optimize the fold, we designed an αββ protein based on a Zn-finger. We then inverse-designed the sequences to provide stability to the fold along with flexibility of linker regions to optimize Zn binding and substrate hydrolysis. The resultant peptides were synthesized and assessed for Zn and substrate binding affinity by fluorescence and ITC followed by evaluation of catalytic efficiency with UV-based enzyme kinetic assays. We were successful in mimicking carbonic anhydrase activity in a peptide of twenty two residues, using p-nitrophenyl acetate as a CO₂ surrogate. Although our design had modest activity, being a simple structure is an advantage for further improvement in efficiency. Our approach opens a way forward to evolving an efficient biocatalyst for any industrial reaction of interest.