Prediction of the substrate for nonribosomal peptide synthetase (NRPS) adenylation domains by virtual screening

Nonribosomal peptide synthetases (NRPSs) synthesize a diverse array of bioactive small peptides, many of which are used in medicine. There is considerable interest in predicting NRPS substrate specificity in order to facilitate investigation of the many “cryptic” NRPS genes that have not been linked to any known product. However, the current sequence similarity‐based methods are unable to produce reliable predictions when there is a lack of prior specificity data, which is a particular problem for fungal NRPSs. We conducted virtual screening on the specificity‐determining domain of NRPSs, the adenylation domain, and found that virtual screening using experimentally determined structures results in good enrichment of the cognate substrate. Our results indicate that the conformation of the adenylation domain and in particular the conformation of a key conserved aromatic residue is important in determining the success of the virtual screening. When homology models of NRPS adenylation domains of known specificity, rather than experimentally determined structures, were built and used for virtual screening, good enrichment of the cognate substrate was also achieved in many cases. However, the accuracy of the models was key to the reliability of the predictions and there was a large variation in the results when different models of the same domain were used. This virtual screening approach is promising and is able to produce enrichment of the cognate substrates in many cases, but improvements in building and assessing homology models are required before the approach can be reliably applied to these models. Proteins 2015; 83:2052–2066. © 2015 Wiley Periodicals, Inc.     

The crystal structure of BlmI as a model for nonribosomal peptide synthetase peptidyl carrier proteins

Carrier proteins (CPs) play a critical role in the biosynthesis of various natural products, especially in nonribosomal peptide synthetase (NRPS) and polyketide synthase (PKS) enzymology, where the CPs are referred to as peptidyl‐carrier proteins (PCPs) or acyl‐carrier proteins (ACPs), respectively. CPs can either be a domain in large multifunctional polypeptides or standalone proteins, termed Type I and Type II, respectively. There have been many biochemical studies of the Type I PKS and NRPS CPs, and of Type II ACPs. However, recently a number of Type II PCPs have been found and biochemically characterized. In order to understand the possible interaction surfaces for combinatorial biosynthetic efforts we crystallized the first characterized and representative Type II PCP member, BlmI, from the bleomycin biosynthetic pathway from Streptomyces verticillus ATCC 15003. The structure is similar to CPs in general but most closely resembles PCPs. Comparisons with previously determined PCP structures in complex with catalytic domains reveals a common interaction surface. This surface is highly variable in charge and shape, which likely confers specificity for interactions. Previous nuclear magnetic resonance (NMR) analysis of a prototypical Type I PCP excised from the multimodular context revealed three conformational states. Comparison of the states with the structure of BlmI and other PCPs reveals that only one of the NMR states is found in other studies, suggesting the other two states may not be relevant. The state represented by the BlmI crystal structure can therefore serve as a model for both Type I and Type II PCPs. Proteins 2014; 82:1210–1218. © 2013 Wiley Periodicals, Inc.     

Recent advances in the structural biology of modular polyketide synthases and nonribosomal peptide synthetases

Polyketides and nonribosomal peptides are an important class of natural products with useful bioactivities. These compounds are similarly biosynthesized using enzymes with modular structures despite having different physicochemical properties. These enzymes are attractive targets for bioengineering to produce “unnatural” natural products owing to their modular structures. Therefore, their structures have been studied for a long time; however, the main focus was on truncated-single domains. Surprisingly, there is an increasing number of the structures of whole modules reported, most of which have been enabled through the recent advances in cryogenic electron microscopy technology. In this review, we have summarized the recent advances in the structural elucidation of whole modules.     

Plasmodium falciparum acyl carrier protein crystal structures in disulfide‐linked and reduced states and their prevalence during blood stage growth

Acyl Carrier Protein (ACP) has a single reactive sulfhydryl necessary for function in covalently binding nascent fatty acids during biosynthesis. In Plasmodium falciparum, the causative agent of the most lethal form of malaria, fatty acid biosynthesis occurs in the apicoplast organelle during the liver stage of the parasite life cycle. During the blood stage, fatty acid biosynthesis is inactive and the redox state of the apicoplast has not been determined. We solved the crystal structure of ACP from P. falciparum in reduced and disulfide‐linked forms, and observe the surprising result that the disulfide in the PfACP cross‐linked dimer is sequestered from bulk solvent in a tight molecular interface. We assessed solvent accessibility of the disulfide with small molecule reducing agents and found that the disulfide is protected from BME but less so for other common reducing agents. We examined cultured P. falciparum parasites to determine which form of PfACP is prevalent during the blood stages. We readily detected monomeric PfACP in parasite lysate, but do not observe the disulfide‐linked form, even under conditions of oxidative stress. To demonstrate that PfACP contains a free sulfhydryl and is not acylated or in the apo state, we treated blood stage parasites with the disulfide forming reagent diamide. We found that the effects of diamide are reversed with reducing agent. Together, these results suggest that the apicoplast is a reducing compartment, as suggested by models of P. falciparum metabolism, and that PfACP is maintained in a reduced state during blood stage growth. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Solution structure and proposed domain domain recognition interface of an acyl carrier protein domain from a modular polyketide synthase

Polyketides are a medicinally important class of natural products. The architecture of modular polyketide synthases (PKSs), composed of multiple covalently linked domains grouped into modules, provides an attractive framework for engineering novel polyketide-producing assemblies. However, impaired domain-domain interactions can compromise the efficiency of engineered polyketide biosynthesis. To facilitate the study of these domain-domain interactions, we have used nuclear magnetic resonance (NMR) spectroscopy to determine the first solution structure of an acyl carrier protein (ACP) domain from a modular PKS, 6-deoxyerythronolide B synthase (DEBS). The tertiary fold of this 10-kD domain is a three-helical bundle; an additional short helix in the second loop also contributes to the core helical packing. Superposition of residues 14-94 of the ensemble on the mean structure yields an average atomic RMSD of 0.64 +/- 0.09 Angstrom for the backbone atoms (1.21 +/- 0.13 Angstrom for all non-hydrogen atoms). The three major helices superimpose with a backbone RMSD of 0.48 +/- 0.10 Angstrom (0.99 +/- 0.11 Angstrom for non-hydrogen atoms). Based on this solution structure, homology models were constructed for five other DEBS ACP domains. Comparison of their steric and electrostatic surfaces at the putative interaction interface (centered on helix II) suggests a model for protein-protein recognition of ACP domains, consistent with the previously observed specificity. Site-directed mutagenesis experiments indicate that two of the identified residues influence the specificity of ACP recognition.     

Structure of the terminal PCP domain of the non‐ribosomal peptide synthetase in teicoplanin biosynthesis

The biosynthesis of the glycopeptide antibiotics, of which teicoplanin and vancomycin are representative members, relies on the combination of non‐ribosomal peptide synthesis and modification of the peptide by cytochrome P450 (Oxy) enzymes while the peptide remains bound to the peptide synthesis machinery. We have structurally characterized the final peptidyl carrier protein domain of the teicoplanin non‐ribosomal peptide synthetase machinery: this domain is believed to mediate the interactions with tailoring Oxy enzymes in addition to its function as a shuttle for intermediates between multiple non‐ribosomal peptide synthetase domains. Using solution state NMR, we have determined structures of this PCP domain in two states, the apo and the post‐translationally modified holo state, both of which conform to a four‐helix bundle assembly. The structures exhibit the same general fold as the majority of known carrier protein structures, in spite of the complex biosynthetic role that PCP domains from the final non‐ribosomal peptide synthetase module must play in glycopeptide antibiotic biosynthesis. These structures thus support the hypothesis that it is subtle rearrangements, rather than dramatic conformational changes, which govern carrier protein interactions and selectivity during non‐ribosomal peptide synthesis. Proteins 2015; 83:711–721. © 2015 Wiley Periodicals, Inc.     

Accurate prediction of protein structural classes using functional domains and predicted secondary structure sequences

Protein structural class prediction is one of the challenging problems in bioinformatics. Previous methods directly based on the similarity of amino acid (AA) sequences have been shown to be insufficient for low-similarity protein data-sets. To improve the prediction accuracy for such low-similarity proteins, different methods have been recently proposed that explore the novel feature sets based on predicted secondary structure propensities. In this paper, we focus on protein structural class prediction using combinations of the novel features including secondary structure propensities as well as functional domain (FD) features extracted from the InterPro signature database. Our comprehensive experimental results based on several benchmark data-sets have shown that the integration of new FD features substantially improves the accuracy of structural class prediction for low-similarity proteins as they capture meaningful relationships among AA residues that are far away in protein sequence. The proposed prediction method has also been tested to predict structural classes for partially disordered proteins with the reasonable prediction accuracy, which is a more difficult problem comparing to structural class prediction for commonly used benchmark data-sets and has never been done before to the best of our knowledge. In addition, to avoid overfitting with a large number of features, feature selection is applied to select discriminating features that contribute to achieve high prediction accuracy. The selected features have been shown to achieve stable prediction performance across different benchmark data-sets.     

A structural explanation for the mechanism and specificity of plant branching enzymes I and IIb

Branching enzymes (BEs) are essential in the biosynthesis of starch and glycogen and play critical roles in determining the fine structure of these polymers. The substrates of these BEs are long carbohydrate chains that interact with these enzymes via multiple binding sites on the enzyme’s surface. By controlling the branched-chain length distribution, BEs can mediate the physiological properties of starch and glycogen moieties; however, the mechanism and structural determinants of this specificity remain mysterious. In this study, we identify a large dodecaose binding surface on rice BE I (BEI) that reaches from the outside of the active site to the active site of the enzyme. Mutagenesis activity assays confirm the importance of this binding site in enzyme catalysis, from which we conclude that it is likely the acceptor chain binding site. Comparison of the structures of BE from Cyanothece and BE1 from rice allowed us to model the location of the donor-binding site. We also identified two loops that likely interact with the donor chain and whose sequences diverge between plant BE1, which tends to transfer longer chains, and BEIIb, which transfers exclusively much shorter chains. When the sequences of these loops were swapped with the BEIIb sequence, rice BE1 also became a short-chain transferring enzyme, demonstrating the key role these loops play in specificity. Taken together, these results provide a more complete picture of the structure, selectivity, and activity of BEs.     

Analysis of two linked genes coding for the acyl carrier protein (ACP) from Arabidopsis thaliana (columbia)

Two linked genes, A1 and A2, coding for nearly identical isoforms of the acyl carrier protein (ACP) were isolated from an Arabidopsis thaliana (columbia) genomic library and sequenced. The amino acids deduced from the nucleotide sequence of the two genes indicate they encode distinct transit peptides, but the mature proteins are the same except for residue 79. Both genes are predicted to contain three introns in similar positions, although they differ in sequence and length. The introns interrupt regions coding for a) the transit peptide, b) the junction of the transit peptide and mature protein, and c) the highly conserved domain surrounding serine 38 to which the phosphopantetheine is attached. Primer extension analysis indicates that at least A1 is active in young plants.     

Single mutations toggle the substrate selectivity of multifunctional Camptotheca secologanic acid synthases

Terpene indole alkaloids (TIAs) are plant-derived specialized metabolites with widespread use in medicine. Species-specific pathways derive various TIAs from common intermediates, strictosidine or strictosidinic acid, produced by coupling tryptamine with secologanin or secologanic acid. The penultimate reaction in this pathway is catalyzed by either secologanin synthase (SLS) or secologanic acid synthase (SLAS) according to whether plants produce secologanin from loganin or secologanic acid from loganic acid. Previous work has identified SLSs and SLASs from different species, but the determinants of selectivity remain unclear. Here, combining molecular modeling, ancestral sequence reconstruction, and biochemical methodologies, we identified key residues that toggle SLS and SLAS selectivity in two CYP72A (cytochrome P450) subfamily enzymes from Camptotheca acuminata. We found that the positions of foremost importance are in substrate recognition sequence 1 (SRS1), where mutations to either of two adjacent histidine residues switched selectivity; His131Phe selects for and increases secologanin production whereas His132Asp selects for secologanic acid production. Furthermore, a change in SRS3 in the predicted substrate entry channel (Arg/Lys270Thr) and another in SRS4 at the start of the I-helix (Ser324Glu) decreased enzyme activity toward either substrate. We propose that the Camptotheca SLASs have maintained the broadened activities found in a common asterid ancestor, even as the Camptotheca lineage lost its ability to produce loganin while the campanulid and lamiid lineages specialized to produce secologanin by acquiring mutations in SRS1. The identification here of the residues essential for the broad substrate scope of SLASs presents opportunities for more tailored heterologous production of TIAs.     

Prediction of protein secondary structure content using amino acid composition and evolutionary information

Knowing protein structure and inferring its function from the structure are one of the main issues of computational structural biology, and often the first step is studying protein secondary structure. There have been many attempts to predict protein secondary structure contents. Previous attempts assumed that the content of protein secondary structure can be predicted successfully using the information on the amino acid composition of a protein. Recent methods achieved remarkable prediction accuracy by using the expanded composition information. The overall average error of the most successful method is 3.4%. Here, we demonstrate that even if we only use the simple amino acid composition information alone, it is possible to improve the prediction accuracy significantly if the evolutionary information is included. The idea is motivated by the observation that evolutionarily related proteins share the similar structure. After calculating the homolog-averaged amino acid composition of a protein, which can be easily obtained from the multiple sequence alignment by running PSI-BLAST, those 20 numbers are learned by a multiple linear regression, an artificial neural network and a support vector regression. The overall average error of method by a support vector regression is 3.3%. It is remarkable that we obtain the comparable accuracy without utilizing the expanded composition information such as pair-coupled amino acid composition. This work again demonstrates that the amino acid composition is a fundamental characteristic of a protein. It is anticipated that our novel idea can be applied to many areas of protein bioinformatics where the amino acid composition information is utilized, such as subcellular localization prediction, enzyme subclass prediction, domain boundary prediction, signal sequence prediction, and prediction of unfolded segment in a protein sequence, to name a few.     

Structures of the first and second double-stranded RNA-binding domains of human TAR RNA-binding protein

The TAR RNA-binding Protein (TRBP) is a double-stranded RNA (dsRNA)-binding protein, which binds to Dicer and is required for the RNA interference pathway. TRBP consists of three dsRNA-binding domains (dsRBDs). The first and second dsRBDs (dsRBD1 and dsRBD2, respectively) have affinities for dsRNA, whereas the third dsRBD (dsRBD3) binds to Dicer. In this study, we prepared the single domain fragments of human TRBP corresponding to dsRBD1 and dsRBD2 and solved the crystal structure of dsRBD1 and the solution structure of dsRBD2. The two structures contain an α-β-β-β-α fold, which is common to the dsRBDs. The overall structures of dsRBD1 and dsRBD2 are similar to each other, except for a slight shift of the first α helix. The residues involved in dsRNA binding are conserved. We examined the small interfering RNA (siRNA)-binding properties of these dsRBDs by isothermal titration colorimetry measurements. The dsRBD1 and dsRBD2 fragments both bound to siRNA, with dissociation constants of 220 and 113 nM, respectively. In contrast, the full-length TRBP and its fragment with dsRBD1 and dsRBD2 exhibited much smaller dissociation constants (0.24 and 0.25 nM, respectively), indicating that the tandem dsRBDs bind simultaneously to one siRNA molecule. On the other hand, the loop between the first α helix and the first β strand of dsRBD2, but not dsRBD1, has a Trp residue, which forms hydrophobic and cation-π interactions with the surrounding residues. A circular dichroism analysis revealed that the thermal stability of dsRBD2 is higher than that of dsRBD1 and depends on the Trp residue.     

Model peptide studies of sequence regions in the elastomeric biomineralization protein, Lustrin A. I. The C-domain consensus-PG-, -NVNCT-motif

The lustrin superfamily represents a unique group of biomineralization proteins localized between layered aragonite mineral plates (i.e., nacre layer) in mollusk shell. Recent atomic force microscopy (AFM) pulling studies have demonstrated that the lustrin‐containing organic nacre layer in the abalone, Haliotis rufescens, exhibits a typical sawtooth force‐extension curve with hysteretic recovery. This force extension behavior is reminiscent of reversible unfolding and refolding in elastomeric proteins such as titin and tenascin. Since secondary structure plays an important role in force‐induced protein unfolding and refolding, the question is, What secondary structure(s) exist within the major domains of Lustrin A? Using a model peptide (FPGKNVNCTSGE) representing the 12‐residue consensus sequence found near the N‐termini of the first eight cysteine‐rich domains (C‐domains) within the Lustrin A protein, we employed CD, NMR spectroscopy, and simulated annealing/minimization to determine the secondary structure preferences for this sequence. At pH 7.4, we find that the 12‐mer sequence adopts a loop conformation, consisting of a “bend” or “turn” involving residues G3–K4 and N7–C8–T9, with extended conformations arising at F1–G3; K4–V6; T9–S10–G11 in the sequence. Minor pH‐dependent conformational effects were noted for this peptide; however, there is no evidence for a salt‐bridge interaction between the K4 and E12 side chains. The presence of a loop conformation within the highly conserved —PG—, —NVNCT— sequence of C1–C8 domains may have important structural and mechanistic implications for the Lustrin A protein with regard to elastic behavior. © 2002 Wiley Periodicals, Inc. Biopolymers 63: 358–369, 2002     

Structural mapping of the coiled‐coil domain of a bacterial condensin and comparative analyses across all domains of life suggest conserved features of SMC proteins

The structural maintenance of chromosomes (SMC) proteins form the cores of multisubunit complexes that are required for the segregation and global organization of chromosomes in all domains of life. These proteins share a common domain structure in which N‐ and C‐ terminal regions pack against one another to form a globular ATPase domain. This “head” domain is connected to a central, globular, “hinge” or dimerization domain by a long, antiparallel coiled coil. To date, most efforts for structural characterization of SMC proteins have focused on the globular domains. Recently, however, we developed a method to map interstrand interactions in the 50‐nm coiled‐coil domain of MukB, the divergent SMC protein found in γ‐proteobacteria. Here, we apply that technique to map the structure of the Bacillus subtilis SMC (BsSMC) coiled‐coil domain. We find that, in contrast to the relatively complicated coiled‐coil domain of MukB, the BsSMC domain is nearly continuous, with only two detectable coiled‐coil interruptions. Near the middle of the domain is a break in coiled‐coil structure in which there are three more residues on the C‐terminal strand than on the N‐terminal strand. Close to the head domain, there is a second break with a significantly longer insertion on the same strand. These results provide an experience base that allows an informed interpretation of the output of coiled‐coil prediction algorithms for this family of proteins. A comparison of such predictions suggests that these coiled‐coil deviations are highly conserved across SMC types in a wide variety of organisms, including humans. Proteins 2015; 83:1027–1045. © 2015 Wiley Periodicals, Inc.