Position-specific propensities of amino acids in the β-strand

RAS-dependent RAF-MAPK hyperactivation by pathogenic RIT1 is a therapeutic target in Noonan syndrome-associated cardiac hypertrophy

RIT1 is a RAS guanosine triphosphatase (GTPase) that regulates different aspects of signal transduction and is mutated in lung cancer, leukemia, and in the germline of individuals with Noonan syndrome. Pathogenic RIT1 proteins promote mitogen-activated protein kinase (MAPK) hyperactivation; however, this mechanism remains poorly understood. Here, we show that RAF kinases are direct effectors of membrane-bound mutant RIT1 necessary for MAPK activation. We identify critical residues in RIT1 that facilitate interaction with membrane lipids and show that these are necessary for association with RAF kinases and MAPK activation. Although mutant RIT1 binds to RAF kinases directly, it fails to activate MAPK signaling in the absence of classical RAS proteins. Consistent with aberrant RAF/MAPK activation as a driver of disease, we show that pathway inhibition alleviates cardiac hypertrophy in a mouse model of RIT1 mutant Noonan syndrome. These data shed light on the function of pathogenic RIT1 and identify avenues for therapeutic intervention.

Comparative Biochemical Studies of Disease-Associated Human Dicer Mutations on Processing of a Pre-microRNA and snoRNA

Dicer is an RNase III enzyme that is responsible for the maturation of small RNAs such as microRNAs. As Dicer's cleavage products play key roles in promoting cellular homeostasis through the fine-tuning of gene expression, dysregulation of Dicer activity can lead to several human diseases, including cancers. Mutations in Dicer have been found to induce tumorigenesis and lead to the development of a rare pleiotropic tumor predisposition syndrome found in children and young adults called DICER1 syndrome. These patients harbor germline and somatic mutations in Dicer that lead to defective microRNA processing and activity. While most mutations occur within Dicer's catalytic RNase III domains, alterations within the Platform-PAZ (Piwi-Argonaute-Zwille) domain also cause loss of microRNA production. Using a combination of in vitro biochemical and cellular studies, we characterized the effect of disease-relevant Platform-PAZ-associated mutations on the processing of a well-studied oncogenic microRNA, pre-microRNA-21. We then compared these results to those of a representative from another Dicer substrate class, the small nucleolar RNA, snord37. From this analysis, we provide evidence that mutations within the Platform-PAZ domain result in differential impacts on RNA binding and processing, adding new insights into the complexities of Dicer processing of small RNA substrates.

Bond clusters control rupture force limit in shear loaded histidine-Ni2+ metal-coordinated proteins

Dynamic noncovalent interactions are pivotal to the structure and function of biological proteins and have been used in bioinspired materials for similar roles. Metal-coordination bonds, in particular, are especially tunable and enable control over static and dynamic properties when incorporated into synthetic materials. Despite growing efforts to engineer metal-coordination bonds to produce strong, tough, and self-healing materials, the systematic characterization of the exact contribution of these bonds towards mechanical strength and the effect of geometric arrangements is missing, limiting the full design potential of these bonds. In this work, we engineer the cooperative rupture of metal-coordination bonds to increase the rupture strength of metal-coordinated peptide dimers. Utilizing all-atom steered molecular dynamics simulations on idealized bidentate histidine-Ni²⁺ coordinated peptides, we show that histidine-Ni²⁺ bonds can rupture cooperatively in groups of two to three bonds. We find that there is a strength limit, where adding additional coordination bonds does not contribute to the additional increase in the protein rupture strength, likely due to the highly heterogeneous rupture behavior exhibited by the coordination bonds. Further, we show that this coordination bond limit is also found natural metal-coordinated biological proteins. Using these insights, we quantitatively suggest how other proteins can be rationally designed with dynamic noncovalent interactions to exhibit cooperative bond breaking behavior. Altogether, this work provides a quantitative analysis of the cooperativity and intrinsic strength limit for metal-coordination bonds with the aim of advancing clear guiding molecular principles for the mechanical design of metal-coordinated materials.

Glycoproteins Involved in Sea Urchin Temporary Adhesion

Biomedical adhesives, despite having been used increasingly in recent years, still face a major technological challenge: strong adhesion in wet environments. In this context, biological adhesives secreted by marine invertebrates have appealing characteristics to incorporate into new underwater biomimetic adhesives: water resistance, nontoxicity and biodegradability. Little is still known about temporary adhesion. Recently, a transcriptomic differential analysis of sea urchin Paracentrotus lividus tube feet pinpointed 16 adhesive/cohesive protein candidates. In addition, it has been demonstrated that the adhesive secreted by this species is composed of high molecular weight proteins associated with N-Acetylglucosamine in a specific chitobiose arrangement. As a follow-up, we aimed to investigate which of these adhesive/cohesive protein candidates were glycosylated through lectin pulldowns, protein identification by mass spectroscopy and in silico characterization. We demonstrate that at least five of the previously identified protein adhesive/cohesive candidates are glycoproteins. We also report the involvement of a third Nectin variant, the first adhesion-related protein to be identified in P. lividus. By providing a deeper characterization of these adhesive/cohesive glycoproteins, this work advances our understanding of the key features that should be replicated in future sea urchin-inspired bioadhesives.

C-terminal determinants for RNA binding motif 7 protein stability and RNA recognition

The 7SK ribonucleoprotein (RNP) is a critical regulator of eukaryotic transcription. Recently, RNA binding motif 7 (RBM7) containing an RNA recognition motif (RRM) was reported to associate with 7SK RNA and core 7SK RNP protein components in response to DNA damage. However, little is known about the mode of RBM7-7SK RNA recognition. Here, we found that RRM constructs containing extended C-termini have increased solubility compared to a minimal RRM construct, although these constructs aggregate in a temperature and concentration-dependent manner. Using solution NMR dynamics experiments, we identified additional structural features observed previously in crystal but not in solution structures. To identify potential RBM7-7SK RNA binding sites, we analyzed deposited data from in cellulo crosslinking experiments and found that RBM7 primarily crosslinks to the distal region of 7SK stem-loop 3 (SL3). Electrophoretic mobility shift assays and NMR chemical shift perturbation experiments showed weak binding to 7SK SL3 constructs in vitro. Together, these results provide new insights into RBM7 RRM folding and recognition of 7SK RNA.

Matching amino acids membrane preference profile to improve activity of antimicrobial peptides

Antimicrobial peptides (AMPs) are cationic antibiotics that can kill multidrug-resistant bacteria via membrane insertion. However, their weak activity limits their clinical use. Ironically, the cationic charge of AMPs is essential for membrane binding, but it obstructs membrane insertion. In this study, we postulate that this problem can be overcome by locating cationic amino acids at the energetically preferred membrane surface. All amino acids have an energetically preferred or less preferred membrane position profile, and this profile is strongly related to membrane insertion. However, most AMPs do not follow this profile. One exception is protegrin-1, a powerful but neglected AMP. In the present study, we found that a potent AMP, WCopW5, strongly resembles protegrin-1 and that the match between its sequence and the preferred position profile closely correlates with its antimicrobial activity. One of its derivatives, WCopW43, has antimicrobial activity comparable to that of the most effective AMPs in clinical use.

Computational Identification and Design of Complementary β-Strand Sequences

The ß-sheet is a regular secondary structure element which consists of linear segments called ß-strands. They are involved in many important biological processes, and some are known to be related to serious diseases such as neurologic disorders and amyloidosis. The self-assembly of ß-sheet peptides also has practical applications in material sciences since they can be building blocks of repeated nanostructures. Therefore, computational algorithms for identification of ß-sheet formation can offer useful insight into the mechanism of disease-prone protein segments and the construction of biocompatible nanomaterials. Despite the recent advances in structure-based methods for the assessment of atomic interactions, identifying amyloidogenic peptides has proven to be extremely difficult since they are structurally very flexible. Thus, an alternative strategy is required to describe ß-sheet formation. It has been hypothesized and observed that there are certain amino acid propensities between ß-strand pairs. Based on this hypothesis, a database search algorithm, B-SIDER, is developed for the identification and design of ß-sheet forming sequences. Given a target sequence, the algorithm identifies exact or partial matches from the structure database and constructs a position-specific score matrix. The score matrix can be utilized to design novel sequences that can form a ß-sheet specifically with the target.

Missense Mutations Modify the Conformational Ensemble of the α-Synuclein Monomer Which Exhibits a Two-Phase Characteristic

α-Synuclein is an intrinsically disordered protein occurring in different conformations and prone to aggregate in β-sheet structures, which are the hallmark of the Parkinson disease. Missense mutations are associated with familial forms of this neuropathy. How these single amino-acid substitutions modify the conformations of wild-type α-synuclein is unclear. Here, using coarse-grained molecular dynamics simulations, we sampled the conformational space of the wild type and mutants (A30P, A53P, and E46K) of α-synuclein monomers for an effective time scale of 29.7 ms. To characterize the structures, we developed an algorithm, CUTABI (CUrvature and Torsion based of Alpha-helix and Beta-sheet Identification), to identify residues in the α-helix and β-sheet from C^α-coordinates. CUTABI was built from the results of the analysis of 14,652 selected protein structures using the Dictionary of Secondary Structure of Proteins (DSSP) algorithm. DSSP results are reproduced with 93% of success for 10 times lower computational cost. A two-dimensional probability density map of α-synuclein as a function of the number of residues in the α-helix and β-sheet is computed for wild-type and mutated proteins from molecular dynamics trajectories. The density of conformational states reveals a two-phase characteristic with a homogeneous phase (state B, β-sheets) and a heterogeneous phase (state HB, mixture of α-helices and β-sheets). The B state represents 40% of the conformations for the wild-type, A30P, and E46K and only 25% for A53T. The density of conformational states of the B state for A53T and A30P mutants differs from the wild-type one. In addition, the mutant A53T has a larger propensity to form helices than the others. These findings indicate that the equilibrium between the different conformations of the α-synuclein monomer is modified by the missense mutations in a subtle way. The α-helix and β-sheet contents are promising order parameters for intrinsically disordered proteins, whereas other structural properties such as average gyration radius, R_g, or probability distribution of R_g cannot discriminate significantly the conformational ensembles of the wild type and mutants. When separated in states B and HB, the distributions of R_g are more significantly different, indicating that global structural parameters alone are insufficient to characterize the conformational ensembles of the α-synuclein monomer.

Secondary Structure of the Novel Myosin Binding Domain WYR and Implications within Myosin Structure

Structural changes in the myosin II light meromyosin (LMM) that influence thick filament mechanical properties and muscle function are modulated by LMM-binding proteins. Flightin is an LMM-binding protein indispensable for the function of Drosophila indirect flight muscle (IFM). Flightin has a three-domain structure that includes WYR, a novel 52 aa domain conserved throughout Pancrustacea. In this study, we (i) test the hypothesis that WYR binds the LMM, (ii) characterize the secondary structure of WYR, and (iii) examine the structural impact WYR has on the LMM. Circular dichroism at 260-190 nm reveals a structural profile for WYR and supports an interaction between WYR and LMM. A WYR-LMM interaction is supported by co-sedimentation with a stoichiometry of ~2.4:1. The WYR-LMM interaction results in an overall increased coiled-coil content, while curtailing ɑ helical content. WYR is found to be composed of 15% turns, 31% antiparallel β, and 48% 'other' content. We propose a structural model of WYR consisting of an antiparallel β hairpin between Q92-K114 centered on an ASX or β turn around N102, with a G1 bulge at G117. The Drosophila LMM segment used, V1346-I1941, encompassing conserved skip residues 2-4, is found to possess a traditional helical profile but is interpreted as having <30% helical content by multiple methods of deconvolution. This low helicity may be affiliated with the dynamic behavior of the structure in solution or the inclusion of a known non-helical region in the C-terminus. Our results support the hypothesis that WYR binds the LMM and that this interaction brings about structural changes in the coiled-coil. These studies implicate flightin, via the WYR domain, for distinct shifts in LMM secondary structure that could influence the structural properties and stabilization of the thick filament, scaling to modulation of whole muscle function.

Positive selection as a key player for SARS-CoV-2 pathogenicity: Insights into ORF1ab, S and E genes

The human β-coronavirus SARS-CoV-2 epidemic started in late December 2019 in Wuhan, China. It causes Covid-19 disease which has become pandemic. Each of the five-known human β-coronaviruses has four major structural proteins (E, M, N and S) and 16 non-structural proteins encoded by ORF1a and ORF1b together (ORF1ab) that are involved in virus pathogenicity and infectivity. Here, we performed detailed positive selection analyses for those six genes among the four previously known human β-coronaviruses and within 38 SARS-CoV-2 genomes to assess signatures of adaptive evolution using maximum likelihood approaches. Our results suggest that three genes (E, S and ORF1ab genes) are under strong signatures of positive selection among human β-coronavirus, influencing codons that are located in functional important protein domains. The E protein-coding gene showed signatures of positive selection in two sites, Asp 66 and Ser 68, located inside a putative transmembrane α-helical domain C-terminal part, which is preferentially composed by hydrophilic residues. Such Asp and Ser sites substitutions (hydrophilic residues) increase the stability of the transmembrane domain in SARS-CoV-2. Moreover, substitutions in the spike (S) protein S1 N-terminal domain have been found, all of them were located on the S protein surface, suggesting their importance in viral transmissibility and survival. Furthermore, evidence of strong positive selection was detected in three of the SARS-CoV-2 nonstructural proteins (NSP1, NSP3, NSP16), which are encoded by ORF1ab and play vital roles in suppressing host translation machinery, viral replication and transcription and inhibiting the host immune response. These results are insightful to assess the role of positive selection in the SARS-CoV-2 encoded proteins, which will allow to better understand the virulent pathogenicity of the virus and potentially identifying targets for drug or vaccine strategy design

Sequence and evolutionary analysis of bacterial ribosomal S1 proteins

The multi-domain bacterial S1 protein is the largest and most functionally important ribosomal protein of the 30S subunit, which interacts with both mRNA and proteins. The family of ribosomal S1 proteins differs in the classical sense from a protein with tandem repeats and has a "bead-on-string" organization, where each repeat is folded into a globular domain. Based on our recent data, the study of evolutionary relationships for the bacterial phyla will provide evidence for one of the proposed theories of the evolutionary development of proteins with structural repeats: from multiple repeats of assembles to single repeats, or vice versa. In this comparative analysis of 1333 S1 sequences that were identified in 24 different phyla, we demonstrate how such phyla can form independently/dependently during evolution. To the best of our knowledge, this work is the first study of the evolutionary history of bacterial ribosomal S1 proteins. The collected and structured data can be useful to computer biologists as a resource for determining percent identity, amino acid composition and logo motifs, as well as dN/dS ratio in bacterial S1 protein. The obtained research data indicate that the evolutionary development of bacterial ribosomal S1 proteins evolved from multiple assemblies to single repeat. The presented data are integrated into the server, which can be accessed at http://oka.protres.ru:4200.

Local sequence-structure relationships in proteins

We seek to understand the interplay between amino acid sequence and local structure in proteins. Are some amino acids unique in their ability to fit harmoniously into certain local structures? What is the role of sequence in sculpting the putative native state folds from myriad possible conformations? In order to address these questions, we represent the local structure of each Cα atom of a protein by just two angles, θ and μ, and we analyze a set of more than 4,000 protein structures from the PDB. We use a hierarchical clustering scheme to divide the 20 amino acids into six distinct groups based on their similarity to each other in fitting local structural space. We present the results of a detailed analysis of patterns of amino acid specificity in adopting local structural conformations and show that the sequence-structure correlation is not very strong compared with a random assignment of sequence to structure. Yet, our analysis may be useful to determine an effective scoring rubric for quantifying the match of an amino acid to its putative local structure.

A systematic analysis of the beta hairpin motif in the Protein Data Bank

The beta hairpin motif is a ubiquitous protein structural motif that can be found in molecules across the tree of life. This motif, which is also popular in synthetically designed proteins and peptides, is known for its stability and adaptability to broad functions. Here, we systematically probe all 49,000 unique beta hairpin substructures contained within the Protein Data Bank (PDB) to uncover key characteristics correlated with stable beta hairpin structure, including amino acid biases and enriched interstrand contacts. We find that position specific amino acid preferences, while seen throughout the beta hairpin structure, are most evident within the turn region, where they depend on subtle turn dynamics associated with turn length and secondary structure. We also establish a set of broad design principles, such as the inclusion of aspartic acid residues at a specific position and the careful consideration of desired secondary structure when selecting residues for the turn region, that can be applied to the generation of libraries encoding proteins or peptides containing beta hairpin structures.

Control over the fibrillization yield by varying the oligomeric nucleation propensities of self-assembling peptides

Self-assembling peptides are an exemplary class of supramolecular biomaterials of broad biomedical utility. Mechanistic studies on the peptide self-assembly demonstrated the importance of the oligomeric intermediates towards the properties of the supramolecular biomaterials being formed. In this study, we demonstrate how the overall yield of the supramolecular assemblies are moderated through subtle molecular changes in the peptide monomers. This strategy is exemplified with a set of surfactant-like peptides (SLPs) with different β-sheet propensities and charged residues flanking the aggregation domains. By integrating different techniques, we show that these molecular changes can alter both the nucleation propensity of the oligomeric intermediates and the thermodynamic stability of the fibril structures. We demonstrate that the amount of assembled nanofibers are critically defined by the oligomeric nucleation propensities. Our findings offer guidance on designing self-assembling peptides for different biomedical applications, as well as insights into the role of protein gatekeeper sequences in preventing amyloidosis.

Molecular determinants of chaperone interactions on MHC-I for folding and antigen repertoire selection

The interplay between a highly polymorphic set of MHC-I alleles and molecular chaperones shapes the repertoire of peptide antigens displayed on the cell surface for T cell surveillance. Here, we demonstrate that the molecular chaperone TAP-binding protein related (TAPBPR) associates with a broad range of partially folded MHC-I species inside the cell. Bimolecular fluorescence complementation and deep mutational scanning reveal that TAPBPR recognition is polarized toward the α₂ domain of the peptide-binding groove, and depends on the formation of a conserved MHC-I disulfide epitope in the α₂ domain. Conversely, thermodynamic measurements of TAPBPR binding for a representative set of properly conformed, peptide-loaded molecules suggest a narrower MHC-I specificity range. Using solution NMR, we find that the extent of dynamics at "hotspot" surfaces confers TAPBPR recognition of a sparsely populated MHC-I state attained through a global conformational change. Consistently, restriction of MHC-I groove plasticity through the introduction of a disulfide bond between the α₁/α₂ helices abrogates TAPBPR binding, both in solution and on a cellular membrane, while intracellular binding is tolerant of many destabilizing MHC-I substitutions. Our data support parallel TAPBPR functions of 1) chaperoning unstable MHC-I molecules with broad allele-specificity at early stages of their folding process, and 2) editing the peptide cargo of properly conformed MHC-I molecules en route to the surface, which demonstrates a narrower specificity. Our results suggest that TAPBPR exploits localized structural adaptations, both near and distant to the peptide-binding groove, to selectively recognize discrete conformational states sampled by MHC-I alleles, toward editing the repertoire of displayed antigens.

B-SIDER: Computational Algorithm for the Design of Complementary β-Sheet Sequences

The β-sheet is an element of protein secondary structure, and intra-/intermolecular β-sheet interactions play pivotal roles in biological regulatory processes including scaffolding, transporting, and oligomerization. In nature, a β-sheet formation is tightly regulated because dysregulated β-stacking often leads to severe diseases such as Alzheimer's, Parkinson's, systemic amyloidosis, or diabetes. Thus, the identification of intrinsic β-sheet-forming propensities can provide valuable insight into protein designs for the development of novel therapeutics. However, structure-based design methods may not be generally applicable to such amyloidogenic peptides mainly owing to high structural plasticity and complexity. Therefore, an alternative design strategy based on complementary sequence information is of significant importance. Herein, we developed a database search method called β-Stacking Interaction DEsign for Reciprocity (B-SIDER) for the design of complementary β-strands. This method makes use of the structural database information and generates target-specific score matrices. The discriminatory power of the B-SIDER score function was tested on representative amyloidogenic peptide substructures against a sequence-based score matrix (PASTA 2.0) and two popular ab initio protein design score functions (Rosetta and FoldX). B-SIDER is able to distinguish wild-type amyloidogenic β-strands as favored interactions in a more consistent manner than other methods. B-SIDER was prospectively applied to the design of complementary β-strands for a splitGFP scaffold. Three variants were identified to have stronger interactions than the original sequence selected through a directed evolution, emitting higher fluorescence intensities. Our results indicate that B-SIDER can be applicable to the design of other β-strands, assisting in the development of therapeutics against disease-related amyloidogenic peptides.

Physicochemical Position-Dependent Properties in the Protein Secondary Structures

BACKGROUND

Establishing theories for designing arbitrary protein structures is complicated and depends on understanding the principles for protein folding, which is affected by applied features. Computer algorithms can reach high precision and stability in computationally designed enzymes and binders by applying informative features obtained from natural structures.

METHODS

In this study, a position-specific analysis of secondary structures (α-helix, β-strand, and tight turn) was performed to reveal novel features for protein structure prediction and protein design.

RESULTS

Our results showed that the secondary structures in the N-termini region tend to be more compact than C-termini. Decoying periodicity in length and distribution of amino acids in α-helices is deciphered using the curve-fitting methods. Compared with α-helix, β-strands do not show distinct periodicities in length. Also, significant differences in position-dependent distribution of physicochemical properties are shown in secondary structures.

CONCLUSION

Position-specific propensities in our study underline valuable parameters that could be used by researchers in the field of structural biology, particularly protein design through site-directed mutagenesis.

Effect of site-directed point mutations on protein misfolding: A simulation study

A Monte Carlo simulation based sequence design method is proposed to investigate the role of site-directed point mutations in protein misfolding. Site-directed point mutations are incorporated in the designed sequences of selected proteins. While most mutated sequences correctly fold to their native conformation, some of them stabilize in other nonnative conformations and thus misfold/unfold. The results suggest that a critical number of hydrophobic amino acid residues must be present in the core of the correctly folded proteins, whereas proteins misfold/unfold if this number of hydrophobic residues falls below the critical limit. A protein can accommodate only a particular number of hydrophobic residues at the surface, provided a large number of hydrophilic residues are present at the surface and critical hydrophobicity of the core is preserved. Some surface sites are observed to be equally sensitive toward site-directed point mutations as the core sites. Point mutations with highly polar and charged amino acids increases the misfold/unfold propensity of proteins. Substitution of natural amino acids at sites with different number of nonbonded contacts suggests that both amino acid identity and its respective site-specificity determine the stability of a protein. A clash-match method is developed to calculate the number of matching and clashing interactions in the mutated protein sequences. While misfolded/unfolded sequences have a higher number of clashing and a lower number of matching interactions, the correctly folded sequences have a lower number of clashing and a higher number of matching interactions. These results are valid for different SCOP classes of proteins.

Insight into the Self-Assembling Properties of Peptergents: A Molecular Dynamics Simulation Study

By manipulating the various physicochemical properties of amino acids, the design of peptides with specific self-assembling properties has been emerging for more than a decade. In this context, short peptides possessing detergent properties (so-called "peptergents") have been developed to self-assemble into well-ordered nanostructures that can stabilize membrane proteins for crystallization. In this study, the peptide with "peptergency" properties, called ADA8 and extensively described by Tao et al., is studied by molecular dynamic simulations for its self-assembling properties in different conditions. In water, it spontaneously forms beta sheets with a β barrel-like structure. We next simulated the interaction of this peptide with a membrane protein, the bacteriorhodopsin, in the presence or absence of a micelle of dodecylphosphocholine. According to the literature, the peptergent ADA8 is thought to generate a belt of β structures around the hydrophobic helical domain that could help stabilize purified membrane proteins. Molecular dynamic simulations are here used to image this mechanism and provide further molecular details for the replacement of detergent molecules around the protein. In addition, we generalized this behavior by designing an amphipathic peptide with beta propensity, which was called ABZ12. Both peptides are able to surround the membrane protein and displace surfactant molecules. To our best knowledge, this is the first molecular mechanism proposed for "peptergency".

Synthetic single domain antibodies for the conformational trapping of membrane proteins

Mechanistic and structural studies of membrane proteins require their stabilization in specific conformations. Single domain antibodies are potent reagents for this purpose, but their generation relies on immunizations, which impedes selections in the presence of ligands typically needed to populate defined conformational states. To overcome this key limitation, we developed an in vitro selection platform based on synthetic single domain antibodies named sybodies. To target the limited hydrophilic surfaces of membrane proteins, we designed three sybody libraries that exhibit different shapes and moderate hydrophobicity of the randomized surface. A robust binder selection cascade combining ribosome and phage display enabled the generation of conformation-selective, high affinity sybodies against an ABC transporter and two previously intractable human SLC transporters, GlyT1 and ENT1. The platform does not require access to animal facilities and builds exclusively on commercially available reagents, thus enabling every lab to rapidly generate binders against challenging membrane proteins.

Stabilization of the Reductase Domain in the Catalytically Self-Sufficient Cytochrome P450BM3 by Consensus-Guided Mutagenesis

The multidomain, catalytically self-sufficient cytochrome P450 BM-3 from Bacillus megaterium (P450BM3 ) constitutes a versatile enzyme for the oxyfunctionalization of organic molecules and natural products. However, the limited stability of the diflavin reductase domain limits the utility of this enzyme for synthetic applications. In this work, a consensus-guided mutagenesis approach was applied to enhance the thermal stability of the reductase domain of P450BM3 . Upon phylogenetic analysis of a set of distantly related P450s (>38 % identity), a total of 14 amino acid substitutions were identified and evaluated in terms of their stabilizing effects relative to the wild-type reductase domain. Recombination of the six most stabilizing mutations generated two thermostable variants featuring up to tenfold longer half-lives at 50 °C and increased catalytic performance at elevated temperatures. Further characterization of the engineered P450BM3 variants indicated that the introduced mutations increased the thermal stability of the FAD-binding domain and that the optimal temperature (Topt ) of the enzyme had shifted from 25 to 40 °C. This work demonstrates the effectiveness of consensus mutagenesis for enhancing the stability of the reductase component of a multidomain P450. The stabilized P450BM3 variants developed here could potentially provide more robust scaffolds for the engineering of oxidation biocatalysts.

Structure–mechanical property correlations of hydrogel forming β-sheet peptides

This review discusses about β-sheet peptide structure at the molecular level and the bulk mechanical properties of the corresponding hydrogels.

Analyses of the general rule on residue pair frequencies in local amino acid sequences of soluble, ordered proteins

The amino acid sequences of soluble, ordered proteins with stable structures have evolved due to biological and physical requirements, thus distinguishing them from random sequences. Previous analyses have focused on extracting the features that frequently appear in protein substructures, such as α-helix and β-sheet, but the universal features of protein sequences have not been addressed. To clarify the differences between native protein sequences and random sequences, we analyzed 7368 soluble, ordered protein sequences, by inspecting the observed and expected occurrences of 400 amino acid pairs in local proximity, up to 10 residues along the sequence in comparison with their expected occurrence in random sequence. We found the trend that the hydrophobic residue pairs and the polar residue pairs are significantly decreased, whereas the pairs between a hydrophobic residue and a polar residue are increased. This trend was universally observed regardless of the secondary structure content but was not observed in protein sequences that include intrinsically disordered regions, indicating that it can be a general rule of protein foldability. The possible benefits of this rule are discussed from the viewpoints of protein aggregation and disorder, which are both caused by low-complexity regions of hydrophobic or polar residues.

Thermodynamics of ideal proteinogenic homopolymer chains as a function of the energy spectrum E, helical propensity ω and enthalpic energy barrier

A reformulation and generalization of the Zwanzig model (ZW model) for ideal homopolymer chains poly-X, where X represents any of the twenty naturally occurring proteinogenic amino acid residues is presented. This reformulation and generalization provides a direct connection between coarse-grained parameters originally proposed in the ZW model with variables from the Lifson-Roig (LR) theory, such as the helical propensity per residue ω, and new variables introduced here, such as the energy gap Δ between unfolded and folded structures, as well as the ratio f of the energy scales involved. This enables us to discover the relevance of the energy spectrum E to the onset of configurational phase transitions. From the configurational partition function Q, thermodynamic properties such as the configurational entropy S, specific heat v and average energy <E> are calculated in terms of the number of residues K, temperature T, helical propensity ω and energy barrier ΔH for different poly-X chains in vacuo. Results obtained here provide substantial evidence that configurational phase transitions for ideal poly-X chains correspond to first-order phase transitions. An anomalous behavior of the thermodynamic functions <E>, Cv, S with respect to the number K of residues is also highlighted. On-going methods of solution are outlined.

Dependence of α-helical and β-sheet amino acid propensities on the overall protein fold type

Background

A large number of studies have been carried out to obtain amino acid propensities for α-helices and β-sheets. The obtained propensities for α-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the β-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects β-sheet propensity.

Results

We calculated amino acid propensities for α-helices and β-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that α-helix propensities do not differ significantly by fold, but β-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for α-helix, but such is not the case for the β-sheet propensities. We also found some fold dependence on amino acid frequency in β-strands. Folds with a high Ser, Thr and Asn content at exposed sites in β-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = −0.90) and to have flat β-sheets. At buried sites in β-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = −0.93). "All-β" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas "α/β" proteins tend to have a higher content of Val, Ile and Leu.

Conclusions

The α-helix propensities are similar for all folds and for exposed and buried residues. However, β-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in β-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a β-strand or association between two β sheets.

Improved sequence-based prediction of strand residues

Accurate identification of strand residues aids prediction and analysis of numerous structural and functional aspects of proteins. We propose a sequence-based predictor, BETArPRED, which improves prediction of strand residues and β-strand segments. BETArPRED uses a novel design that accepts strand residues predicted by SSpro and predicts the remaining positions utilizing a logistic regression classifier with nine custom-designed features. These are derived from the primary sequence, the secondary structure (SS) predicted by SSpro, PSIPRED and SPINE, and residue depth as predicted by RDpred. Our features utilize certain local (window-based) patterns in the predicted SS and combine information about the predicted SS and residue depth. BETArPRED is evaluated on 432 sequences that share low identity with the training chains, and on the CASP8 dataset. We compare BETArPRED with seven modern SS predictors, and the top-performing automated structure predictor in CASP8, the ZHANG-server. BETArPRED provides statistically significant improvements over each of the SS predictors; it improves prediction of strand residues and β-strands, and it finds β-strands that were missed by the other methods. When compared with the ZHANG-server, we improve predictions of strand segments and predict more actual strand residues, while the other predictor achieves higher rate of correct strand residue predictions when under-predicting them.