Structure and dynamics of the first archaeal parvulin reveal a new functionally important loop in parvulin-type prolyl isomerases

Pitfalls in measurements of R1 relaxation rates of protein backbone 15N nuclei

The dynamics of the backbone and side-chains of protein are routinely studied by interpreting experimentally determined 15N spin relaxation rates. R1(15N), the longitudinal relaxation rate, reports on fast motions and encodes, together with the transverse relaxation R2, structural information about the shape of the molecule and the orientation of the amide bond vectors in the internal diffusion frame. Determining error-free 15N longitudinal relaxation rates remains a challenge for small, disordered, and medium-sized proteins. Here, we show that mono-exponential fitting is sufficient, with no statistical preference for bi-exponential fitting up to 800 MHz. A detailed comparison of the TROSY and HSQC techniques at medium and high fields showed no statistically significant differences. The least error-prone DD/CSA interference removal technique is the selective inversion of amide signals while avoiding water resonance. The exchange of amide with solvent deuterons appears to affect the rate R1 of solvent-exposed amides in all fields tested and in each DD/CSA interference removal technique in a statistically significant manner. In summary, the most accurate R1(15N) rates in proteins are achieved by selective amide inversion, without the addition of D2O. Importantly, at high magnetic fields stronger than 800 MHz, when non-mono-exponential decay is involved, it is advisable to consider elimination of the shortest delays (typically up to 0.32 s) or bi-exponential fitting.

AlphaFold and what is next: bridging functional, systems and structural biology

INTRODUCTION

The DeepMind's AlphaFold (AF) has revolutionized biomedical and biocience research by providing both experts and non-experts with an invaluable tool for predicting protein structures. However, while AF is highly effective for predicting structures of rigid and globular proteins, it is not able to fully capture the dynamics, conformational variability, and interactions of proteins with ligands and other biomacromolecules.

AREAS COVERED

In this review, we present a comprehensive overview of the latest advancements in 3D model predictions for biomacromolecules using AF. We also provide a detailed analysis its of strengths and limitations, and explore more recent iterations, modifications, and practical applications of this strategy. Moreover, we map the path forward for expanding the landscape of AF toward predicting structures of every protein and peptide, and their interactions in the proteome in the most physiologically relevant form. This discussion is based on an extensive literature search performed using PubMed and Google Scholar.

EXPERT OPINION

While significant progress has been made to enhance AF's modeling capabilities, we argue that a combined approach integrating both various in silico and in vitro methods will be most beneficial for the future of structural biology, bridging the gaps between static and dynamic features of proteins and their functions.

Rapid Interpretation of Protein Backbone Rotation Dynamics Directly from Spin Relaxation Data

Besides their structure, dynamics is pivotal for protein functions, particularly for intrinsically disordered proteins (IDPs) that do not fold into a fixed 3D structure. However, the detection of protein dynamics is difficult for IDPs and other disordered biomolecules. NMR spin relaxation rates are sensitive to the rapid rotations of chemical bonds, but their interpretation is arduous for IDPs or molecular assemblies with a complex dynamic landscape. Here we demonstrate numerically that the dynamics of a wide range of proteins, from short peptides to partially disordered proteins and peptides in micelles, can be characterized by calculating the total effective correlation times of protein backbone N-H bond rotations, τeff, from experimentally measured transverse 15N spin relaxation rates, R2, using a linear relation. Our results enable the determination of magnetic-field-independent and intuitively understandable parameters describing protein dynamics at different regions of the sequence directly from experiments. A practical advance of the approach is demonstrated by analyzing partially disordered proteins in which rotations of disordered regions occur with timescales of 1-2 ns, independent of their size, suggesting that rotations of disordered and folded regions are uncoupled in these proteins.

Clinical Implications of MiR128, Angiotensin I Converting Enzyme and Vascular Endothelial Growth Factor Gene Abnormalities and Their Association with T2D

Type 2 DM (T2D) results from the interaction of the genetic and environmental risk factors. Vascular endothelial growth factor (VEGF), angiotensin I-converting enzyme (ACE), and MicroRNAs (MiRNAs) are involved in important physiological processes. Gene variations in VEGF, ACE and MiRNA genes are associated with diseases. In this study we investigated the associations of the VEGF-2578 C/A (rs699947), VEGF-2549 insertion/deletion (I/D), and ACE I/D rs4646994 and Mir128a (rs11888095) gene variations with T2D using the amplification refractory mutation system PCR (ARMS-PCR) and mutation specific PCR (MSP). We screened 122 T2D cases and 126 healthy controls (HCs) for the rs699947, and 133 T2D cases and 133 HCs for the VEGF I/D polymorphism. For the ACE I/D we screened 152 cases and 150 HCs, and we screened 129 cases and 112 HCs for the Mir128a (rs11888095). The results showed that the CA genotype of the VEGF rs699947 and D allele of the VEGF I/D polymorphisms were associated with T2D with OR =2.01, p-value = 0.011, and OR = 2.42, p-value = 0.010, respectively. The result indicated the D allele of the ACE ID was protective against T2D with OR = 0.10, p-value = 0.0001, whereas the TC genotype and the T allele of the Mir128a (rs11888095) were associated with increased risk to T2D with OR = 3.16, p-value = 0.0001, and OR = 1.68, p-value = 0.01, respectively. We conclude that the VEGF (rs699947), VEGF I/D and Mir128a (rs11888095) are potential risk loci for T2D, and that the D allele of the ACE ID polymorphism may be protective against T2D. These results help in identification and stratification for the individuals that at risk for T2D. However, future well-designed studies in different populations and with larger sample sizes are required. Moreover, studies to examine the effects of these polymorphisms on VEGF and ACE proteins are recommended.

Distribution of Peptidyl-Prolyl Isomerase (PPIase) in the Archaea

Cis-trans isomerization of the peptide bond prior to proline is an intrinsically slow process but plays an essential role in protein folding. In vivo cis-trans isomerization reaction is catalyzed by Peptidyl-prolyl isomerase (PPIases), a category of proteins widely distributed among all the three domains of life. The present study is majorly focused on the distribution of different types of PPIases in the archaeal domain. All the three hitherto known families of PPIases (namely FKBP, Cyclophilin and parvulin) were studied to identify the evolutionary conservation across the phylum archaea. The basic function of cyclophilin, FKBP and parvulin has been conserved whereas the sequence alignment suggested variations in each clade. The conserved residues within the predicted motif of each family are unique. The available protein structures of different PPIase across various domains were aligned to ascertain the structural variation in the catalytic site. The structural alignment of native PPIase proteins among various groups suggested that the apo-protein may have variable conformations but when bound to their specific inhibitors, they attain similar active site configuration. This is the first study of its kind which explores the distribution of archaeal PPIases, along with detailed structural and functional analysis of each type of PPIase found in archaea.

Tracing Eukaryotic Ribosome Biogenesis Factors Into the Archaeal Domain Sheds Light on the Evolution of Functional Complexity

Ribosome assembly is an essential and carefully choreographed cellular process. In eukaryotes, several 100 proteins, distributed across the nucleolus, nucleus, and cytoplasm, co-ordinate the step-wise assembly of four ribosomal RNAs (rRNAs) and approximately 80 ribosomal proteins (RPs) into the mature ribosomal subunits. Due to the inherent complexity of the assembly process, functional studies identifying ribosome biogenesis factors and, more importantly, their precise functions and interplay are confined to a few and very well-established model organisms. Although best characterized in yeast (Saccharomyces cerevisiae), emerging links to disease and the discovery of additional layers of regulation have recently encouraged deeper analysis of the pathway in human cells. In archaea, ribosome biogenesis is less well-understood. However, their simpler sub-cellular structure should allow a less elaborated assembly procedure, potentially providing insights into the functional essentials of ribosome biogenesis that evolved long before the diversification of archaea and eukaryotes. Here, we use a comprehensive phylogenetic profiling setup, integrating targeted ortholog searches with automated scoring of protein domain architecture similarities and an assessment of when search sensitivity becomes limiting, to trace 301 curated eukaryotic ribosome biogenesis factors across 982 taxa spanning the tree of life and including 727 archaea. We show that both factor loss and lineage-specific modifications of factor function modulate ribosome biogenesis, and we highlight that limited sensitivity of the ortholog search can confound evolutionary conclusions. Projecting into the archaeal domain, we find that only few factors are consistently present across the analyzed taxa, and lineage-specific loss is common. While members of the Asgard group are not special with respect to their inventory of ribosome biogenesis factors (RBFs), they unite the highest number of orthologs to eukaryotic RBFs in one taxon. Using large ribosomal subunit maturation as an example, we demonstrate that archaea pursue a simplified version of the corresponding steps in eukaryotes. Much of the complexity of this process evolved on the eukaryotic lineage by the duplication of ribosomal proteins and their subsequent functional diversification into ribosome biogenesis factors. This highlights that studying ribosome biogenesis in archaea provides fundamental information also for understanding the process in eukaryotes.

Molecular Determination of mirRNA-126 rs4636297, Phosphoinositide-3-Kinase Regulatory Subunit 1-Gene Variability rs7713645, rs706713 (Tyr73Tyr), rs3730089 (Met326Ile) and Their Association with Susceptibility to T2D

Type 2 diabetes is a metabolic disease characterized by elevated blood sugar. It has serious complications and socioeconomic impact. The MicroRNAs are short single-stranded and non-coding RNA molecules. They regulate gene expression at the post-transcriptional levels. They are important for many physiological processes including metabolism, growth, and others. The phosphoinositide 3-kinase (PI3K) is important for insulin signaling and glucose uptake. The genome wide association studies have identified the association of certain loci with diseases including T2D. In this study we have examined the association of miR126 rs4636297 and Phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1) gene Variations rs7713645, rs706713 (Tyr73Tyr), and rs3730089 (Met326Ile) with T2D using the amplification refractory mutation system PCR. Results indicated that there was a significant different (p-value < 0.05) in the Mir126 rs4636297 genotypes distribution between cases and controls, and the minor allele of the rs4636297 was also associated with T2D with OR = 0.58, p-value < 0.05. In addition results showed that there were significant differences (p-value < 0.05) of rs4636297 genotype distribution of patients with normal and patient with abnormal lipid profile. Results also showed that the PIK3R1 rs7713645 and rs3730089 genotype distribution was significantly different between cases and controls with a p-values < 0.05. In addition, the minor allele of the rs7713645 and rs3730089 were associated with T2D with OR = 0.58, p-value < 0.05. We conclude that the Mir126 rs4636297 and PIK3R1 SNPs (rs7713645 and rs3730089) were associated with T2D. These results need verification in future studies with larger sample sizes and in different populations. Protein-protein interaction and enzyme assay studies are also required to uncover the effect of the SNPs on the PI3K regulatory subunit (PI3KR1) and PI3K catalytic activity.

Structural Analysis of the SANT/Myb Domain of FLASH and YARP Proteins and Their Complex with the C-Terminal Fragment of NPAT by NMR Spectroscopy and Computer Simulations

FLICE-associated huge protein (FLASH), Yin Yang 1-Associated Protein-Related Protein (YARP) and Nuclear Protein, Ataxia-Telangiectasia Locus (NPAT) localize to discrete nuclear structures called histone locus bodies (HLBs) where they control various steps in histone gene expression. Near the C-terminus, FLASH and YARP contain a highly homologous domain that interacts with the C-terminal region of NPAT. Structural aspects of the FLASH-NPAT and YARP-NPAT complexes and their role in histone gene expression remain largely unknown. In this study, we used multidimensional NMR spectroscopy and in silico modeling to analyze the C-terminal domain in FLASH and YARP in an unbound form and in a complex with the last 31 amino acids of NPAT. Our results demonstrate that FLASH and YARP domains share the same fold of a triple α -helical bundle that resembles the DNA binding domain of Myb transcriptional factors and the SANT domain found in chromatin-modifying and remodeling complexes. The NPAT peptide contains a single α -helix that makes multiple contacts with α -helices I and III of the FLASH and YARP domains. Surprisingly, in spite of sharing a significant amino acid similarity, each domain likely binds NPAT using a unique network of interactions, yielding two distinct complexes. In silico modeling suggests that both complexes are structurally compatible with DNA binding, raising the possibility that they may function in identifying specific sequences within histone gene clusters, hence initiating the assembly of HLBs and regulating histone gene expression during cell cycle progression.

Biophysical insight into the interaction of levocabastine with human serum albumin: spectroscopy and molecular docking approach

Interaction of levocabastine with human serum albumin (HSA) is investigated by applying fluorescence spectroscopy, circular dichroism spectroscopy and molecular docking methods. Levocabastine is an important drug in treatment of allergy and currently a target drug for drug repurposing to treat other diseases like vernal keratoconjuctivitis. Fluorescence quenching data revealed that levocabastine bind weakly to protein with binding constant in the order of 10³ M^-1. Förster resonance energy transfer results indicated the binding distance of 2.28 nm for levocabastine. Synchronous fluorescence result suggest slight blue shift for tryptophan upon levocabastine binding, binding of levocabastine impelled rise in α-helical structure in protein, while there are minimal changes in tertiary structure in protein. Moreover, docking results indicate levocabastine binds to pocket near to the drug site-I in HSA via hydrogen bonding and hydrophobic interactions. Understanding the interaction of levocabastine with HSA is significant for the advancement of therapeutic and diagnostic strategies for optimal treatment results.Communicated by Ramaswamy H. Sarma.

Coarse-grain simulations on NMR conformational ensembles highlight functional residues in proteins

Dynamics are a key feature of protein function, and this is especially true of gating residues, which occupy cavity or tunnel lining positions in the protein structure, and will reversibly switch between open and closed conformations in order to control the diffusion of small molecules within a protein's internal matrix. Earlier work on globins and hydrogenases have shown that these gating residues can be detected using a multiscale scheme combining all-atom classic molecular dynamics simulations and coarse-grain calculations of the resulting conformational ensemble mechanical properties. Here, we show that the structural variations observed in the conformational ensembles produced by NMR spectroscopy experiments are sufficient to induce noticeable mechanical changes in a protein, which in turn can be used to identify residues important for function and forming a mechanical nucleus in the protein core. This new approach, which combines experimental data and rapid coarse-grain calculations and no longer needs to resort to time-consuming all-atom simulations, was successfully applied to five different protein families.

H-NS uses an autoinhibitory conformational switch for environment-controlled gene silencing

As an environment-dependent pleiotropic gene regulator in Gram-negative bacteria, the H-NS protein is crucial for adaptation and toxicity control of human pathogens such as Salmonella, Vibrio cholerae or enterohaemorrhagic Escherichia coli. Changes in temperature affect the capacity of H-NS to form multimers that condense DNA and restrict gene expression. However, the molecular mechanism through which H-NS senses temperature and other physiochemical parameters remains unclear and controversial. Combining structural, biophysical and computational analyses, we show that human body temperature promotes unfolding of the central dimerization domain, breaking up H-NS multimers. This unfolding event enables an autoinhibitory compact H-NS conformation that blocks DNA binding. Our integrative approach provides the molecular basis for H-NS-mediated environment-sensing and may open new avenues for the control of pathogenic multi-drug resistant bacteria.

Fine-tuning the extent and dynamics of binding cleft opening as a potential general regulatory mechanism in parvulin-type peptidyl prolyl isomerases

Parvulins or rotamases form a distinct group within peptidyl prolyl cis-trans isomerases. Their exact mode of action as well as the role of conserved residues in the family are still not unambiguously resolved. Using backbone S² order parameters and NOEs as restraints, we have generated dynamic structural ensembles of three distinct parvulins, SaPrsA, TbPin1 and CsPinA. The resulting ensembles are in good agreement with the experimental data but reveal important differences between the three enzymes. The largest difference can be attributed to the extent of the opening of the substrate binding cleft, along which motional mode the three molecules occupy distinct regions. Comparison with a wide range of other available parvulin structures highlights structural divergence along the bottom of the binding cleft acting as a hinge during the opening-closing motion. In the prototype WW-domain containing parvulin, Pin1, this region is also important in forming contacts with the WW domain known to modulate enzymatic activity of the catalytic domain. We hypothesize that modulation of the extent and dynamics of the identified ‘breathing motion’ might be one of the factors responsible for functional differences in the distinct parvulin subfamilies.

Discovery, Molecular Mechanisms, and Industrial Applications of Cold-Active Enzymes

Cold-active enzymes constitute an attractive resource for biotechnological applications. Their high catalytic activity at temperatures below 25°C makes them excellent biocatalysts that eliminate the need of heating processes hampering the quality, sustainability, and cost-effectiveness of industrial production. Here we provide a review of the isolation and characterization of novel cold-active enzymes from microorganisms inhabiting different environments, including a revision of the latest techniques that have been used for accomplishing these paramount tasks. We address the progress made in the overexpression and purification of cold-adapted enzymes, the evolutionary and molecular basis of their high activity at low temperatures and the experimental and computational techniques used for their identification, along with protein engineering endeavors based on these observations to improve some of the properties of cold-adapted enzymes to better suit specific applications. We finally focus on examples of the evaluation of their potential use as biocatalysts under conditions that reproduce the challenges imposed by the use of solvents and additives in industrial processes and of the successful use of cold-adapted enzymes in biotechnological and industrial applications.

NmPin from the marine thaumarchaeote Nitrosopumilus maritimus is an active membrane associated prolyl isomerase

Background

Peptidyl-prolyl isomerases (PPIases) are present in all forms of life and play a crucial role in protein folding and regulation. They catalyze the cis-trans isomerization of the peptide bond that precedes proline residues in numerous proteins. The parvulins, which is one family of PPIases, have been extensively investigated in several eukaryotes. However, nothing is known about their expression, function and localization in archaea.

Results

Here, we describe the endogenous expression, molecular structure, function and cellular localization of NmPin, a single-domain parvulin-type PPIase from Nitrosopumilus maritimus. This marine chemolithoautotrophic archaeon belongs to the globally abundant phylum Thaumarchaeota. Using high resolution NMR spectroscopy we demonstrate that the 3D structure of NmPin adopts a parvulin fold and confirmed its peptidyl-prolyl isomerase activity by protease-coupled assays and mutagenesis studies. A detailed topological analysis revealed a positively charged lysine-rich patch on the protein surface, which is conserved in all known parvulin sequences of thaumarchaeotes and targets NmPin to lipids in vitro. Immunofluorescence microscopy confirms that the protein is attached to the outer archaeal cell membrane in vivo. Transmission electron microscopy uncovered that NmPin has a uniform distribution at the membrane surface, which is correlated with a native cell shape of the prokaryote.

Conclusion

We present a novel solution structure of a catalytically active thaumarchaeal parvulin. Our results reveal that a lysine-rich patch in NmPin mediates membrane localization. These findings provide a model whereby NmPin is located between the archaeal membrane and the surface layer and hence suggest proteins of the S-layer as the key target substrates of this parvulin.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0274-1) contains supplementary material, which is available to authorized users.

Parvulin 17-catalyzed Tubulin Polymerization Is Regulated by Calmodulin in a Calcium-dependent Manner

Recently we have shown that the peptidyl-prolyl cis/trans isomerase parvulin 17 (Par17) interacts with tubulin in a GTP-dependent manner, thereby promoting the formation of microtubules. Microtubule assembly is regulated by Ca(2+)-loaded calmodulin (Ca(2+)/CaM) both in the intact cell and under in vitro conditions via direct interaction with microtubule-associated proteins. Here we provide the first evidence that Ca(2+)/CaM interacts also with Par17 in a physiologically relevant way, thus preventing Par17-promoted microtubule assembly. In contrast, parvulin 14 (Par14), which lacks only the first 25 N-terminal residues of the Par17 sequence, does not interact with Ca(2+)/CaM, indicating that this interaction is exclusive for Par17. Pulldown experiments and chemical shift perturbation analysis with (15)N-labeled Par17 furthermore confirmed that calmodulin (CaM) interacts in a Ca(2+)-dependent manner with the Par17 N terminus. The reverse experiment with (15)N-labeled Ca(2+)/CaM demonstrated that the N-terminal Par17 segment binds to both CaM lobes simultaneously, indicating that Ca(2+)/CaM undergoes a conformational change to form a binding channel between its two lobes, apparently similar to the structure of the CaM-smMLCK(796-815) complex. In vitro tubulin polymerization assays furthermore showed that Ca(2+)/CaM completely suppresses Par17-promoted microtubule assembly. The results imply that Ca(2+)/CaM binding to the N-terminal segment of Par17 causes steric hindrance of the Par17 active site, thus interfering with the Par17/tubulin interaction. This Ca(2+)/CaM-mediated control of Par17-assisted microtubule assembly may provide a mechanism that couples Ca(2+) signaling with microtubule function.

The Ess1 prolyl isomerase: traffic cop of the RNA polymerase II transcription cycle

Ess1 is a prolyl isomerase that regulates the structure and function of eukaryotic RNA polymerase II. Ess1 works by catalyzing the cis/trans conversion of pSer5-Pro6 bonds, and to a lesser extent pSer2-Pro3 bonds, within the carboxy-terminal domain (CTD) of Rpb1, the largest subunit of RNA pol II. Ess1 is conserved in organisms ranging from yeast to humans. In budding yeast, Ess1 is essential for growth and is required for efficient transcription initiation and termination, RNA processing, and suppression of cryptic transcription. In mammals, Ess1 (called Pin1) functions in a variety of pathways, including transcription, but it is not essential. Recent work has shown that Ess1 coordinates the binding and release of CTD-binding proteins that function as co-factors in the RNA pol II complex. In this way, Ess1 plays an integral role in writing (and reading) the so-called CTD code to promote production of mature RNA pol II transcripts including non-coding RNAs and mRNAs.

Transient domain interactions enhance the affinity of the mitotic regulator Pin1 toward phosphorylated peptide ligands

The mitotic regulator Pin1 plays an important role in protein quality control and age-related medical conditions such as Alzheimer disease and Parkinson disease. Although its cellular role has been thoroughly investigated during the past decade, the molecular mechanisms underlying its function remain elusive. We provide evidence for interactions between the two domains of Pin1. Several residues displayed unequivocal peak splits in nuclear magnetic resonance spectra, indicative of two different conformational states in equilibrium. Pareto analysis of paramagnetic relaxation enhancement data demonstrates that the two domains approach each other upon addition of a nonpeptidic ligand. Titration experiments with phosphorylated peptides monitored by fluorescence anisotropy and chemical shift perturbation indicate that domain interactions increase Pin1's affinity toward peptide ligands. We propose this interplay of the domains and ligands to be a general mechanism for a large class of two-domain proteins.

TSAR: a program for automatic resonance assignment using 2D cross-sections of high dimensionality, high-resolution spectra

While NMR studies of proteins typically aim at structure, dynamics or interactions, resonance assignments represent in almost all cases the initial step of the analysis. With increasing complexity of the NMR spectra, for example due to decreasing extent of ordered structure, this task often becomes both difficult and time-consuming, and the recording of high-dimensional data with high-resolution may be essential. Random sampling of the evolution time space, combined with sparse multidimensional Fourier transform (SMFT), allows for efficient recording of very high dimensional spectra (≥4 dimensions) while maintaining high resolution. However, the nature of this data demands for automation of the assignment process. Here we present the program TSAR (Tool for SMFT-based Assignment of Resonances), which exploits all advantages of SMFT input. Moreover, its flexibility allows to process data from any type of experiments that provide sequential connectivities. The algorithm was tested on several protein samples, including a disordered 81-residue fragment of the δ subunit of RNA polymerase from Bacillus subtilis containing various repetitive sequences. For our test examples, TSAR achieves a high percentage of assigned residues without any erroneous assignments.

Solution structural analysis of the single-domain parvulin TbPin1

Background

Pin1-type parvulins are phosphorylation-dependent peptidyl-prolyl cis-trans isomerases. Their functions have been widely reported to be involved in a variety of cellular responses or processes, such as cell division, transcription, and apoptosis, as well as in human diseases including Alzheimer's disease and cancers. TbPin1 was identified as a novel class of Pin1-type parvulins from Trypanosoma brucei, containing a unique PPIase domain, which can catalyze the isomerization of phosphorylated Ser/Thr-Pro peptide bond.

Methodology/Principal Findings

We determined the solution structure of TbPin1 and performed ¹⁵N relaxation measurements to analyze its backbone dynamics using multi-dimensional heteronuclear NMR spectroscopy. The average RMSD values of the 20 lowest energy structures are 0.50±0.05 Å for backbone heavy atoms and 0.85±0.08 Å for all heavy atoms. TbPin1 adopts the typical catalytic tertiary structure of Pin1-type parvulins, which comprises a globular fold with a four-stranded anti-parallel β-sheet core surrounded by three α-helices and one 3₁₀-helix. The global structure of TbPin1 is relatively rigid except the active site. The 2D EXSY spectra illustrate that TbPin1 possesses a phosphorylation-dependent PPIase activity. The binding sites of TbPin1 for a phosphorylated peptide substrate {SSYFSG[p]TPLEDDSD} were determined by the chemical shift perturbation approach. Residues Ser15, Arg18, Asn19, Val21, Ser22, Val32, Gly66, Ser67, Met83, Asp105 and Gly107 are involved in substantial contact with the substrate.

Conclusions/Significance

The solution structure of TbPin1 and the binding sites of the phosphorylated peptide substrate on TbPin1 were determined. The work is helpful for further understanding the molecular basis of the substrate specificity for Pin1-type parvulin family and enzyme catalysis.

Single-domain parvulins constitute a specific marker for recently proposed deep-branching archaeal subgroups

Peptidyl-prolyl cis/trans isomerases (PPIases) are enzymes assisting protein folding and protein quality control in organisms of all kingdoms of life. In contrast to the other sub-classes of PPIases, the cyclophilins and the FK-506 binding proteins, little was formerly known about the parvulin type of PPIase in Archaea. Recently, the first solution structure of an archaeal parvulin, the PinA protein from Cenarchaeum symbiosum, was reported. Investigation of occurrence and frequency of PPIase sequences in numerous archaeal genomes now revealed a strong tendency for thermophilic microorganisms to reduce the number of PPIases. Single-domain parvulins were mostly found in the genomes of recently proposed deep-branching archaeal subgroups, the Thaumarchaeota and the ARMANs (archaeal Richmond Mine acidophilic nanoorganisms). Hence, we used the parvulin sequence to reclassify available archaeal metagenomic contigs, thereby, adding new members to these subgroups. A combination of genomic background analysis and phylogenetic approaches of parvulin sequences suggested that the assigned sequences belong to at least two distinct groups of Thaumarchaeota. Finally, machine learning approaches were applied to identify amino acid residues that separate archaeal and bacterial parvulin proteins from each other. When mapped onto the recent PinA solution structure, most of these positions form a cluster at one site of the protein possibly indicating a different functionality of the two groups of parvulin proteins.

FK506-Binding protein 22 from a psychrophilic bacterium, a cold shock-inducible peptidyl prolyl isomerase with the ability to assist in protein folding

Adaptation of microorganisms to low temperatures remains to be fully elucidated. It has been previously reported that peptidyl prolyl cis-trans isomerases (PPIases) are involved in cold adaptation of various microorganisms whether they are hyperthermophiles, mesophiles or phsycrophiles. The rate of cis-trans isomerization at low temperatures is much slower than that at higher temperatures and may cause problems in protein folding. However, the mechanisms by which PPIases are involved in cold adaptation remain unclear. Here we used FK506-binding protein 22, a cold shock protein from the psychrophilic bacterium Shewanella sp. SIB1 (SIB1 FKBP22) as a model protein to decipher the involvement of PPIases in cold adaptation. SIB1 FKBP22 is homodimer that assumes a V-shaped structure based on a tertiary model. Each monomer consists of an N-domain responsible for dimerization and a C-catalytic domain. SIB1 FKBP22 is a typical cold-adapted enzyme as indicated by the increase of catalytic efficiency at low temperatures, the downward shift in optimal temperature of activity and the reduction in the conformational stability. SIB1 FKBP22 is considered as foldase and chaperone based on its ability to catalyze refolding of a cis-proline containing protein and bind to a folding intermediate protein, respectively. The foldase and chaperone activites of SIB1 FKBP22 are thought to be important for cold adaptation of Shewanella sp. SIB1. These activities are also employed by other PPIases for being involved in cold adaptation of various microorganisms. Despite other biological roles of PPIases, we proposed that foldase and chaperone activities of PPIases are the main requirement for overcoming the cold-stress problem in microorganisms due to folding of proteins.