Multi-template approach to modeling engineered disulfide bonds

Navigating the landscape of enzyme design: from molecular simulations to machine learning

Global environmental issues and sustainable development call for new technologies for fine chemical synthesis and waste valorization. Biocatalysis has attracted great attention as the alternative to the traditional organic synthesis. However, it is challenging to navigate the vast sequence space to identify those proteins with admirable biocatalytic functions. The recent development of deep-learning based structure prediction methods such as AlphaFold2 reinforced by different computational simulations or multiscale calculations has largely expanded the 3D structure databases and enabled structure-based design. While structure-based approaches shed light on site-specific enzyme engineering, they are not suitable for large-scale screening of potential biocatalysts. Effective utilization of big data using machine learning techniques opens up a new era for accelerated predictions. Here, we review the approaches and applications of structure-based and machine-learning guided enzyme design. We also provide our view on the challenges and perspectives on effectively employing enzyme design approaches integrating traditional molecular simulations and machine learning, and the importance of database construction and algorithm development in attaining predictive ML models to explore the sequence fitness landscape for the design of admirable biocatalysts.

Impact of secretin receptor homo-dimerization on natural ligand binding

Class B G protein-coupled receptors can form dimeric complexes important for high potency biological effects. Here, we apply pharmacological, biochemical, and biophysical techniques to cells and membranes expressing the prototypic secretin receptor (SecR) to gain insights into secretin binding to homo-dimeric and monomeric SecR. Spatial proximity between peptide and receptor residues, probed by disulfide bond formation, demonstrates that the secretin N-terminus moves from adjacent to extracellular loop 3 (ECL3) at wild type SecR toward ECL2 in non-dimerizing mutants. Analysis of fluorescent secretin analogs demonstrates stable engagement of the secretin C-terminal region within the receptor extracellular domain (ECD) for both dimeric and monomeric receptors, while the mid-region exhibits lower mobility while docked at the monomer. Moreover, decoupling of G protein interaction reduces mobility of the peptide mid-region at wild type receptor to levels similar to the mutant, whereas it has no further impact on the monomer. These data support a model of peptide engagement whereby the ability of SecR to dimerize promotes higher conformational dynamics of the peptide-bound receptor ECD and ECLs that likely facilitates more efficient G protein recruitment and activation, consistent with the higher observed functional potency of secretin at wild type SecR relative to the monomeric mutant receptor.

Construction and characterization of a Myceliophthora thermophila lytic polysaccharide monooxygenase mutant S174C/A93C with improved thermostability

Lytic polysaccharide monooxygenases (LPMOs) can oxidatively cleave the glycosidic bonds of crystalline polysaccharides, providing more accessible sites for polysaccharide hydrolases and promoting efficient conversion of biomass. In order to promote industrial applications of LPMOs, the stability of an LPMO of Myceliophthora thermophila C1 (MtC1LPMO) was improved by adding disulfide bonds in this study. Firstly, the structural changes of wild-type (WT) MtC1LPMO at different temperatures were explored using molecular dynamics simulations, and eight mutants were selected by combining the predicted results from Disulfide by Design (DBD), Multi agent stability prediction upon point mutations (Maestro) and Bridge disulfide (BridgeD) websites. Then, the enzymatic properties of the different mutants were determined after their expression and purification, and the mutant S174C/A93C with the highest thermal stability was obtained. The specific activities of unheated S174C/A93C and WT were 160.6 ± 1.7 U/g and 174.8 ± 7.5 U/g, respectively, while those of S174C/A93C and WT treated at 70 °C for 4 h were 77.7 ± 3.4 U/g and 46.1 ± 0.4 U/g, respectively. The transition midpoint temperature of S174C/A93C was 2.7 °C higher than that of WT. The conversion efficiency of S174C/A93C for both microcrystalline cellulose and corn straw was about 1.5 times higher than that of WT. Finally, molecular dynamics simulations revealed that the introduction of disulfide bonds increased the β-sheet content of the H1-E34 region, thus improving the rigidity of the protein. Therefore, the overall structural stability of S174C/A93C was improved, which in turn improved its thermal stability.

Alteration of Chain-Length Selectivity and Thermostability of Rhizopus oryzae Lipase via Virtual Saturation Mutagenesis Coupled with Disulfide Bond Design

Rhizopus oryzae lipase (ROL) is one of the most important enzymes used in the food, biofuel, and pharmaceutical industries. However, the highly demanding conditions of industrial processes can reduce its stability and activity. To seek a feasible method to improve both the catalytic activity and the thermostability of this lipase, first, the structure of ROL was divided into catalytic and noncatalytic regions by identifying critical amino acids in the crevice-like binding pocket. Second, a mutant screening library aimed at improvement of ROL catalytic performance by virtual saturation mutagenesis of residues in the catalytic region was constructed based on Rosetta's Cartesian_ddg protocol. A double mutant, E265V/S267W (with an E-to-V change at residue 265 and an S-to-W change at residue 267), with markedly improved catalytic activity toward diverse chain-length fatty acid esters was identified. Then, computational design of disulfide bonds was conducted for the noncatalytic amino acids of E265V/S267W, and two potential disulfide bonds, S61C-S115C and E190C-E238C, were identified as candidates. Experimental data validated that the variant E265V/S267W/S61C-S115C/E190C-E238C had superior stability, with an increase of 8.5°C in the melting temperature and a half-life of 31.7 min at 60°C, 4.2-fold longer than that of the wild-type enzyme. Moreover, the variant improved the lipase activity toward five 4-nitrophenyl esters by 1.5 to 3.8 times, exhibiting a potential to modify the catalytic efficiency. IMPORTANCE Rhizopus oryzae lipase (ROL) is very attractive in biotechnology and industry as a safe and environmentally friendly biocatalyst. Functional expression of ROL in Escherichia coli facilitates effective high-throughput screening for positive variants. This work highlights a method to improve both selectivity and thermostability based on a combination of virtual saturation mutagenesis in the substrate pocket and disulfide bond prediction in the noncatalytic region. Using the method, ROL thermostability and activity to diverse 4-nitrophenyl esters could be substantially improved. The strategy of rational introduction of multiple mutations in different functional domains of the enzyme is a great prospect in the modification of biocatalysts.

Engineering the Thermostability of the Mono- and Diacylglycerol Lipase SMG1 for the Synthesis of Diacylglycerols

Diacylglycerols (DAGs) display huge application prospectives in food industries. Therefore, new strategies to produce diacylglycerides are needed. Malassezia globose lipase (SMG1) could be used to synthesize DAGs. However, the poor thermostability of SMG1 seriously hampers its application. Herein, a rational design was used to generate a more thermostable SMG1. Compared with the wild type (WT), the M5D mutant (Q34P/A37P/M176V/G177A/M294R/ G28C-P206C), which contains five single-point mutations and one additional disulfide bond, displayed a 14.0 °C increase in the melting temperature (T_m), 5 °C in the optimal temperature, and 1154.3-fold in the half-life (t_1/2) at 55 °C. Meanwhile, the specific activity towards DAGs of the M5D variant was improved by 3.0-fold compared to the WT. Molecular dynamics (MD) simulations revealed that the M5D mutant showed an improved rigid structure. Additionally, the WT and the M5D variants were immobilized and used for the production of DAGs. Compared with the WT, the immobilized M5D-catalyzed esterification showed a 9.1% higher DAG content and a 22.9% increase in residual activity after nine consecutive cycles. This study will pave the way for the industrial application of SMG1.

Improving Both the Thermostability and Catalytic Efficiency of Phospholipase D from Moritella sp. JT01 through Disulfide Bond Engineering Strategy

Mining of Phospholipase D (PLD) with high activity and stability has attracted strong interest for investigation. A novel PLD from marine Moritella sp. JT01 (MsPLD) was biochemically and structurally characterized in our previous study; however, the short half-life time (t_1/2) under its optimum reaction temperature seriously hampered its further applications. Herein, the disulfide bond engineering strategy was applied to improve its thermostability. Compared with wild-type MsPLD, mutant S148C-T206C/D225C-A328C with the addition of two disulfide bonds exhibited a 3.1-fold t_1/2 at 35 °C and a 5.7 °C increase in melting temperature (T_m). Unexpectedly, its specific activity and catalytic efficiency (k_cat/K_m) also increased by 22.7% and 36.5%, respectively. The enhanced activity might be attributed to an increase in the activation entropy by displacing more water molecules by the transition state. The results of molecular dynamics simulations (MD) revealed that the introduction of double disulfide bonds rigidified the global structure of the mutant, which might cause the enhanced thermostability. Finally, the synthesis capacity of the mutant to synthesize phosphatidic acid (PA) was evaluated. The conversion rate of PA reached about 80% after 6 h reaction with wild-type MsPLD but reached 78% after 2 h with mutant S148C-T206C/D225C-A328C, which significantly reduced the time needed for the reaction to reach equilibrium. The present results pave the way for further application of MsPLD in the food and pharmaceutical industries.

Harnessing the potential of bacterial oxidative folding to aid protein production

Protein folding is central to both biological function and recombinant protein production. In bacterial expression systems, which are easy to use and offer high protein yields, production of the protein of interest in its native fold can be hampered by the limitations of endogenous posttranslational modification systems. Disulfide bond formation, entailing the covalent linkage of proximal cysteine amino acids, is a fundamental posttranslational modification reaction that often underpins protein stability, especially in extracytoplasmic environments. When these bonds are not formed correctly, the yield and activity of the resultant protein are dramatically decreased. Although the mechanism of oxidative protein folding is well understood, unwanted or incorrect disulfide bond formation often presents a stumbling block for the expression of cysteine-containing proteins in bacteria. It is therefore important to consider the biochemistry of prokaryotic disulfide bond formation systems in the context of protein production, in order to take advantage of the full potential of such pathways in biotechnology applications. Here, we provide a critical overview of the use of bacterial oxidative folding in protein production so far, and propose a practical decision-making workflow for exploiting disulfide bond formation for the expression of any given protein of interest.

Capturing Peptide-GPCR Interactions and Their Dynamics

Many biological functions of peptides are mediated through G protein-coupled receptors (GPCRs). Upon ligand binding, GPCRs undergo conformational changes that facilitate the binding and activation of multiple effectors. GPCRs regulate nearly all physiological processes and are a favorite pharmacological target. In particular, drugs are sought after that elicit the recruitment of selected effectors only (biased ligands). Understanding how ligands bind to GPCRs and which conformational changes they induce is a fundamental step toward the development of more efficient and specific drugs. Moreover, it is emerging that the dynamic of the ligand–receptor interaction contributes to the specificity of both ligand recognition and effector recruitment, an aspect that is missing in structural snapshots from crystallography. We describe here biochemical and biophysical techniques to address ligand–receptor interactions in their structural and dynamic aspects, which include mutagenesis, crosslinking, spectroscopic techniques, and mass-spectrometry profiling. With a main focus on peptide receptors, we present methods to unveil the ligand–receptor contact interface and methods that address conformational changes both in the ligand and the GPCR. The presented studies highlight a wide structural heterogeneity among peptide receptors, reveal distinct structural changes occurring during ligand binding and a surprisingly high dynamics of the ligand–GPCR complexes.

Structure and dynamics of the active Gs-coupled human secretin receptor

The class B secretin GPCR (SecR) has broad physiological effects, with target potential for treatment of metabolic and cardiovascular disease. Molecular understanding of SecR binding and activation is important for its therapeutic exploitation. We combined cryo-electron microscopy, molecular dynamics, and biochemical cross-linking to determine a 2.3 Å structure, and interrogate dynamics, of secretin bound to the SecR:Gs complex. SecR exhibited a unique organization of its extracellular domain (ECD) relative to its 7-transmembrane (TM) core, forming more extended interactions than other family members. Numerous polar interactions formed between secretin and the receptor extracellular loops (ECLs) and TM helices. Cysteine-cross-linking, cryo-electron microscopy multivariate analysis and molecular dynamics simulations revealed that interactions between peptide and receptor were dynamic, and suggested a model for initial peptide engagement where early interactions between the far N-terminus of the peptide and SecR ECL2 likely occur following initial binding of the peptide C-terminus to the ECD.

Yosshi: a web-server for disulfide engineering by bioinformatic analysis of diverse protein families

Disulfide bonds play a significant role in protein stability, function or regulation but are poorly conserved among evolutionarily related proteins. The Yosshi can help to understand the role of S–S bonds by comparing sequences and structures of homologs with diverse properties and different disulfide connectivity patterns within a common structural fold of a superfamily, and assist to select the most promising hot-spots to improve stability of proteins/enzymes or modulate their functions by introducing naturally occurring crosslinks. The bioinformatic analysis is supported by the integrated Mustguseal web-server to construct large structure-guided sequence alignments of functionally diverse protein families that can include thousands of proteins based on all available information in public databases. The Yosshi+Mustguseal is a new integrated web-tool for a systematic homology-driven analysis and engineering of S–S bonds that facilitates a broader interpretation of disulfides not just as a factor of structural stability, but rather as a mechanism to implement functional diversity within a superfamily. The results can be downloaded as a content-rich PyMol session file or further studied online using the HTML5-based interactive analysis tools. Both web-servers are free and open to all users at https://biokinet.belozersky.msu.ru/yosshi and there is no login requirement.

Enhancing the Thermostability of Rhizomucor miehei Lipase with a Limited Screening Library by Rational-Design Point Mutations and Disulfide Bonds

Rhizomucor miehei lipase (RML), as a kind of eukaryotic protein catalyst, plays an important role in the food, organic chemical, and biofuel industries. However, RML retains its catalytic activity below 50°C, which limits its industrial applications at higher temperatures. Soluble expression of this eukaryotic protein in Escherichia coli not only helps to screen for thermostable mutants quickly but also provides the opportunity to develop rapid and effective ways to enhance the thermal stability of eukaryotic proteins. Therefore, in this study, RML was engineered using multiple computational design methods, followed by filtration via conservation analysis and functional region assessment. We successfully obtained a limited screening library (only 36 candidates) to validate thermostable single point mutants, among which 24 of the candidates showed higher thermostability and 13 point mutations resulted in an apparent melting temperature ([Formula: see text]) of at least 1°C higher. Furthermore, both of the two disulfide bonds predicted from four rational-design algorithms were further introduced and found to stabilize RML. The most stable mutant, with T18K/T22I/E230I/S56C-N63C/V189C-D238C mutations, exhibited a 14.3°C-higher [Formula: see text] and a 12.5-fold increase in half-life at 70°C. The catalytic efficiency of the engineered lipase was 39% higher than that of the wild type. The results demonstrate that rationally designed point mutations and disulfide bonds can effectively reduce the number of screened clones to enhance the thermostability of RML.IMPORTANCER. miehei lipase, whose structure is well established, can be widely applied in diverse chemical processes. Soluble expression of R. miehei lipase in E. coli provides an opportunity to explore efficient methods for enhancing eukaryotic protein thermostability. This study highlights a strategy that combines computational algorithms to predict single point mutations and disulfide bonds in RML without losing catalytic activity. Through this strategy, an RML variant with greatly enhanced thermostability was obtained. This study provides a competitive alternative for wild-type RML in practical applications and further a rapid and effective strategy for thermostability engineering.

Engineering of a membrane-triggered activity switch in coagulation factor VIIa

Coagulation factor VIIa (FVIIa) is an intrinsically poor serine protease that requires assistance from its cofactor tissue factor (TF) to trigger the extrinsic pathway of blood coagulation. TF stimulates FVIIa through allosteric maturation of its active site and by facilitating substrate recognition. The surface dependence of the latter property allowed us to design a potent membrane-triggered activity switch in FVIIa by engineering a disulfide cross-link between an allosterically silent FVIIa variant and soluble TF. These results show that optimization of substrate recognition remote from the active site represents a promising new route to simultaneously enhance and localize the procoagulant activity of FVIIa for therapeutic purposes.

Recombinant factor VIIa (FVIIa) variants with increased activity offer the promise to improve the treatment of bleeding episodes in patients with inhibitor-complicated hemophilia. Here, an approach was adopted to enhance the activity of FVIIa by selectively optimizing substrate turnover at the membrane surface. Under physiological conditions, endogenous FVIIa engages its cell-localized cofactor tissue factor (TF), which stimulates activity through membrane-dependent substrate recognition and allosteric effects. To exploit these properties of TF, a covalent complex between FVIIa and the soluble ectodomain of TF (sTF) was engineered by introduction of a nonperturbing cystine bridge (FVIIa Q64C-sTF G109C) in the interface. Upon coexpression, FVIIa Q64C and sTF G109C spontaneously assembled into a covalent complex with functional properties similar to the noncovalent wild-type complex. Additional introduction of a FVIIa-M306D mutation to uncouple the sTF-mediated allosteric stimulation of FVIIa provided a final complex with FVIIa-like activity in solution, while exhibiting a two to three orders-of-magnitude increase in activity relative to FVIIa upon exposure to a procoagulant membrane. In a mouse model of hemophilia A, the complex normalized hemostasis upon vascular injury at a dose of 0.3 nmol/kg compared with 300 nmol/kg for FVIIa.

Disulfide Trapping for Modeling and Structure Determination of Receptor: Chemokine Complexes

Despite the recent breakthrough advances in GPCR crystallography, structure determination of protein-protein complexes involving chemokine receptors and their endogenous chemokine ligands remains challenging. Here, we describe disulfide trapping, a methodology for generating irreversible covalent binary protein complexes from unbound protein partners by introducing two cysteine residues, one per interaction partner, at selected positions within their interaction interface. Disulfide trapping can serve at least two distinct purposes: (i) stabilization of the complex to assist structural studies and/or (ii) determination of pairwise residue proximities to guide molecular modeling. Methods for characterization of disulfide-trapped complexes are described and evaluated in terms of throughput, sensitivity, and specificity toward the most energetically favorable crosslinks. Due to abundance of native disulfide bonds at receptor:chemokine interfaces, disulfide trapping of their complexes can be associated with intramolecular disulfide shuffling and result in misfolding of the component proteins; because of this, evidence from several experiments is typically needed to firmly establish a positive disulfide crosslink. An optimal pipeline that maximizes throughput and minimizes time and costs by early triage of unsuccessful candidate constructs is proposed.

An allosteric model for control of pore opening by substrate binding in the EutL microcompartment shell protein

The ethanolamine utilization (Eut) microcompartment is a protein-based metabolic organelle that is strongly associated with pathogenesis in bacteria that inhabit the human gut. The exterior shell of this elaborate protein complex is composed from a few thousand copies of BMC-domain shell proteins, which form a semi-permeable diffusion barrier that provides the interior enzymes with substrates and cofactors while simultaneously retaining metabolic intermediates. The ability of this protein shell to regulate passage of substrate and cofactor molecules is critical for microcompartment function, but the details of how this diffusion barrier can allow the passage of large cofactors while still retaining small intermediates remain unclear. Previous work has revealed two conformations of the EutL shell protein, providing substantial evidence for a gated pore that might allow the passage of large cofactors. Here we report structural and biophysical evidence to show that ethanolamine, the substrate of the Eut microcompartment, acts as a negative allosteric regulator of EutL pore opening. Specifically, a series of X-ray crystal structures of EutL from Clostridium perfringens, along with equilibrium binding studies, reveal that ethanolamine binds to EutL at a site that exists in the closed-pore conformation and which is incompatible with opening of the large pore for cofactor transport. The allosteric mechanism we propose is consistent with the cofactor requirements of the Eut microcompartment, leading to a new model for EutL function. Furthermore, our results suggest the possibility of redox modulation of the allosteric mechanism, opening potentially new lines of investigation.

Experiment-Guided Molecular Modeling of Protein-Protein Complexes Involving GPCRs

Experimental structure determination for G protein-coupled receptors (GPCRs) and especially their complexes with protein and peptide ligands is at its infancy. In the absence of complex structures, molecular modeling and docking play a large role not only by providing a proper 3D context for interpretation of biochemical and biophysical data, but also by prospectively guiding experiments. Experimentally confirmed restraints may help improve the accuracy and information content of the computational models. Here we present a hybrid molecular modeling protocol that integrates heterogeneous experimental data with force field-based calculations in the stochastic global optimization of the conformations and relative orientations of binding partners. Some experimental data, such as pharmacophore-like chemical fields or disulfide-trapping restraints, can be seamlessly incorporated in the protocol, while other types of data are more useful at the stage of solution filtering. The protocol was successfully applied to modeling and design of a stable construct that resulted in crystallization of the first complex between a chemokine and its receptor. Examples from this work are used to illustrate the steps of the protocol. The utility of different types of experimental data for modeling and docking is discussed and caveats associated with data misinterpretation are highlighted.

Conformational dynamics of individual antibodies using computational docking and AFM

Molecular recognition between a receptor and a ligand requires a certain level of flexibility in macromolecules. In this study, we aimed at analyzing the conformational variability of receptors portrayed by monoclonal antibodies that have been individually imaged using atomic force microscopy (AFM). Individual antibodies were chemically coupled to activated mica surface, and they have been imaged using AFM in ambient conditions. The resulting topographical surface of antibodies was used to assemble the three subunits constituting antibodies: two antigen-binding fragments and one crystallizable fragment using a surface-constrained computational docking approach. Reconstructed structures based on 10 individual topographical surfaces of antibodies are presented for which separation and relative orientation of the subunits were measured. When compared with three X-ray structures of antibodies present in the protein data bank database, results indicate that several arrangements of the reconstructed subunits are comparable with those of known structures. Nevertheless, no reconstructed structure superimposes adequately to any particular X-ray structure consequence of the antibody flexibility. We conclude that high-resolution AFM imaging with appropriate computational reconstruction tools is adapted to study the conformational dynamics of large individual macromolecules deposited on mica.

Computationally designed libraries for rapid enzyme stabilization

The ability to engineer enzymes and other proteins to any desired stability would have wide-ranging applications. Here, we demonstrate that computational design of a library with chemically diverse stabilizing mutations allows the engineering of drastically stabilized and fully functional variants of the mesostable enzyme limonene epoxide hydrolase. First, point mutations were selected if they significantly improved the predicted free energy of protein folding. Disulfide bonds were designed using sampling of backbone conformational space, which tripled the number of experimentally stabilizing disulfide bridges. Next, orthogonal in silico screening steps were used to remove chemically unreasonable mutations and mutations that are predicted to increase protein flexibility. The resulting library of 64 variants was experimentally screened, which revealed 21 (pairs of) stabilizing mutations located both in relatively rigid and in flexible areas of the enzyme. Finally, combining 10-12 of these confirmed mutations resulted in multi-site mutants with an increase in apparent melting temperature from 50 to 85°C, enhanced catalytic activity, preserved regioselectivity and a >250-fold longer half-life. The developed Framework for Rapid Enzyme Stabilization by Computational libraries (FRESCO) requires far less screening than conventional directed evolution.

Hyperthin nanochains composed of self-polymerizing protein shackles

Protein fibrils are expected to have applications as functional nanomaterials because of their sophisticated structures; however, nanoscale ordering of the functional units of protein fibrils remains challenging. Here we design a series of self-polymerizing protein monomers, referred to as protein shackles, derived from modified recombinant subunits of pili from Streptococcus pyogenes. The monomers polymerize into nanochains through spontaneous irreversible covalent bond formation. We design the protein shackles so that their reactions can be controlled by altering redox conditions, which affect disulphide bond formation between engineered cysteine residues. The interaction between the monomers improves their polymerization reactivity and determines morphologies of the polymers. In addition, green fluorescent protein-tagged protein shackles can polymerize, indicating proteins can be stably attached to the nanochains with its functionality preserved. Furthermore we demonstrate that a molecular-recognizable nanochain binds to its partner with an enhanced binding ability in solution. These characteristics are expected to be applied for novel protein nanomaterials.

Understanding the evolutionary structural variability and target specificity of tick salivary Kunitz peptides using next generation transcriptome data

BACKGROUND

Ticks are blood-sucking arthropods and a primary function of tick salivary proteins is to counteract the host's immune response. Tick salivary Kunitz-domain proteins perform multiple functions within the feeding lesion and have been classified as venoms; thereby, constituting them as one of the important elements in the arms race with the host. The two main mechanisms advocated to explain the functional heterogeneity of tick salivary Kunitz-domain proteins are gene sharing and gene duplication. Both do not, however, elucidate the evolution of the Kunitz family in ticks from a structural dynamic point of view. The Red Queen hypothesis offers a fruitful theoretical framework to give a dynamic explanation for host-parasite interactions. Using the recent salivary gland Ixodes ricinus transcriptome we analyze, for the first time, single Kunitz-domain encoding transcripts by means of computational, structural bioinformatics and phylogenetic approaches to improve our understanding of the structural evolution of this important multigenic protein family.

RESULTS

Organizing the I. ricinus single Kunitz-domain peptides based on their cysteine motif allowed us to specify a putative target and to relate this target specificity to Illumina transcript reads during tick feeding. We observe that several of these Kunitz peptide groups vary in their translated amino acid sequence, secondary structure, antigenicity, and intrinsic disorder, and that the majority of these groups are subject to a purifying (negative) selection. We finalize by describing the evolution and emergence of these Kunitz peptides. The overall interpretation of our analyses discloses a rapidly emerging Kunitz group with a distinct disulfide bond pattern from the I. ricinus salivary gland transcriptome.

CONCLUSIONS

We propose a model to explain the structural and functional evolution of tick salivary Kunitz peptides that we call target-oriented evolution. Our study reveals that combining analytical approaches (transcriptomes, computational, bioinformatics and phylogenetics) improves our understanding of the biological functions of important salivary gland mediators during tick feeding.

Protein disulfide engineering

Improving the stability of proteins is an important goal in many biomedical and industrial applications. A logical approach is to emulate stabilizing molecular interactions found in nature. Disulfide bonds are covalent interactions that provide substantial stability to many proteins and conform to well-defined geometric conformations, thus making them appealing candidates in protein engineering efforts. Disulfide engineering is the directed design of novel disulfide bonds into target proteins. This important biotechnological tool has achieved considerable success in a wide range of applications, yet the rules that govern the stabilizing effects of disulfide bonds are not fully characterized. Contrary to expectations, many designed disulfide bonds have resulted in decreased stability of the modified protein. We review progress in disulfide engineering, with an emphasis on the issue of stability and computational methods that facilitate engineering efforts.

Nanoscale structural features determined by AFM for single virus particles

In this work, we propose "single-image analysis", as opposed to multi-image averaging, for extracting valuable information from AFM images of single bio-particles. This approach allows us to study molecular systems imaged by AFM under general circumstances without restrictions on their structural forms. As feature exhibition is a resolution correlation, we have performed AFM imaging on surfaces of tobacco mosaic virus (TMV) to demonstrate variations of structural patterns with probing resolution. Two AFM images were acquired with the same tip at different probing resolutions in terms of pixel width, i.e., 1.95 and 0.49 nm per pixel. For assessment, we have constructed an in silico topograph based on the three-dimensional crystal structure of TMV as a reference. The prominent artifacts observed in the AFM-determined shape of TMV were attributed to tip convolutions. The width of TMV rod was systematically overestimated by ~10 nm at both probing resolutions of AFM. Nevertheless, the effects of tip convolution were less severe in vertical orientation so that the estimated height of TMV by AFM imaging was in close agreement with the in silico X-ray topograph. Using dedicated image processing algorithms, we found that at low resolution (i.e., 1.95 nm per pixel), the extracted surface features of TMV can be interpreted as a partial or full helical repeat (three complete turns with ~7.0 nm in length), while individual protein subunits (~2.5 nm) were perceivable only at high resolution. The present study shows that the scales of revealed structural features in AFM images are subject to both probing resolution and processing algorithms for image analysis.

Proteolysin, a novel highly thermostable and cosolvent-compatible protease from the thermophilic bacterium Coprothermobacter proteolyticus

Through genome mining, we identified a gene encoding a putative serine protease of the thermitase subgroup of subtilases (EC 3.4.21.66) in the thermophilic bacterium Coprothermobacter proteolyticus. The gene was functionally expressed in Escherichia coli, and the enzyme, which we called proteolysin, was purified to near homogeneity from crude cell lysate by a single heat treatment step. Proteolysin has a broad pH tolerance and is active at temperatures of up to 80°C. In addition, the enzyme shows good activity and stability in the presence of organic solvents, detergents, and dithiothreitol, and it remains active in 6 M guanidinium hydrochloride. Based on its stability and activity profile, proteolysin can be an excellent candidate for applications where resistance to harsh process conditions is required.

Eukaryotic GPN-loop GTPases paralogs use a dimeric assembly reminiscent of archeal GPN

GTPases are molecular switches that regulate a wide-range of cellular processes. The GPN-loop GTPase (GPN) is a sub-family of P-loop NTPase that evolved from a single gene copy in archaea to triplicate paralog genes in eukaryotes, each having a non-redundant essential function in cell. In Saccharomyces cerevisiae, yGPN1 and yGPN2 are involved in sister chromatid cohesion mechanism, whereas nothing is known regarding yGPN3 function. Previous high-throughput experiments suggested that GPN paralogs interaction may occur. In this work, GPN|GPN contact was analyzed in details using TAP-Tag approach, yeast two-hybrid assay, in silico energy computation and site-directed mutagenesis of a conserved Glu residue located at the center of the interaction interface. It is demonstrated that this residue is essential for cell viability. A chromatid cohesion assay revealed that, like yGPN1 and yGPN2, yGPN3 also plays a role in sister chromatid cohesion. These results suggest that all three GPN proteins act at the molecular level in sister chromatid cohesion mechanism as a GPN|GPN complex reminiscent of the homodimeric structure of PAB0955, an archaeal member of GPN-loop GTPase.

CING: an integrated residue-based structure validation program suite

We present a suite of programs, named CING for Common Interface for NMR Structure Generation that provides for a residue-based, integrated validation of the structural NMR ensemble in conjunction with the experimental restraints and other input data. External validation programs and new internal validation routines compare the NMR-derived models with empirical data, measured chemical shifts, distance- and dihedral restraints and the results are visualized in a dynamic Web 2.0 report. A red-orange-green score is used for residues and restraints to direct the user to those critiques that warrant further investigation. Overall green scores below ~20 % accompanied by red scores over ~50 % are strongly indicative of poorly modelled structures. The publically accessible, secure iCing webserver ( https://nmr.le.ac.uk ) allows individual users to upload the NMR data and run a CING validation analysis.

Mapping spatial approximations between the amino terminus of secretin and each of the extracellular loops of its receptor using cysteine trapping

While it is evident that the carboxyl-terminal region of natural peptide ligands bind to the amino-terminal domain of class B GPCRs, how their biologically critical amino-terminal regions dock to the receptor is unclear. We utilize cysteine trapping to systematically explore spatial approximations among residues in the first five positions of secretin and in every position within the receptor extracellular loops (ECLs). Only Cys(2) and Cys(5) secretin analogues exhibited full activity and retained moderate binding affinity (IC(50): 92±4 and 83±1 nM, respectively). When these peptides probed 61 human secretin receptor cysteine-replacement mutants, a broad network of receptor residues could form disulfide bonds consistent with a dynamic ligand-receptor interface. Two distinct patterns of disulfide bond formation were observed: Cys(2) predominantly labeled residues in the amino terminus of ECL2 and ECL3 (relative labeling intensity: Ser(340), 94±7%; Pro(341), 84±9%; Phe(258), 73±5%; Trp(274) 62±8%), and Cys(5) labeled those in the carboxyl terminus of ECL2 and ECL3 (Gln(348), 100%; Ile(347), 73±12%; Glu(342), 59±10%; Phe(351), 58±11%). These constraints were utilized in molecular modeling, providing improved understanding of the structure of the transmembrane bundle and interconnecting loops, the orientation between receptor domains, and the molecular basis of ligand docking. Key spatial approximations between peptide and receptor predicted by this model (H(1)-W(274), D(3)-N(268), G(4)-F(258)) were supported by mutagenesis and residue-residue complementation studies.

Functional mapping of factor VIII C2 domain

The factor VIII (FVIII) is a cofactor of the coagulation cascade. The FVIII C2 domain is a critical domain that participates in the interactions with the von Willebrand factor and the phospholipidic surfaces. To assess the importance of each residue of this domain in the maintenance of the structure and the function of FVIII, a number (n=139) of mutants were generated by substituting the original residues, from Ser2173 to Gly2325, by an alanine. Mutants were built within a complete B domain-deleted FVIII and expressed in COS-1 cells. Mutant antigen levels and procoagulant activities were measured. Two in silico analyses, a sliding average procedure and an analysis of the mutation energy cost were conducted in parallel on the FVIII structure. Both results were in agreement with the functional data, and illustrated the benefit of using such strategies prior to targeting specific residues in the aim of generating active recombinant molecules. The functional assays identify the residues that are important to maintaining the structure of the C2 domain, mainly those forming β-sheet, and those that can afford substitution, establishing a detailed functional relation with the available crystallographic data. This study provided a comprehensive functional mapping of the FVIII C2 domain and discussed the implication of specific residues in respect to the maintenance in the activity and structure stability, the efficiency in secretion, the binding to phospholipids and the formation of epitope.

Engineered human antibody constant domains with increased stability

The immunoglobulin (Ig) constant CH2 domain is critical for antibody effector functions. Isolated CH2 domains are promising as scaffolds for construction of libraries containing diverse binders that could also confer some effector functions. However, previous work has shown that an isolated murine CH2 domain is relatively unstable to thermally induced unfolding. To explore unfolding mechanisms of isolated human CH2 and increase its stability gamma1 CH2 was cloned and a panel of cysteine mutants was constructed. Human gamma1 CH2 unfolded at a higher temperature (T(m) = 54.1 degrees C, as measured by circular dichroism) than that previously reported for a mouse CH2 (41 degrees C). One mutant (m01) was remarkably stable (T(m) = 73.8 degrees C). Similar results were obtained by differential scanning calorimetry. This mutant was also significantly more stable than the wild-type CH2 against urea induced unfolding (50% unfolding at urea concentration of 6.8 m versus 4.2 m). The m01 was highly soluble and monomeric. The existence of the second disulfide bond in m01 and its correct position were demonstrated by mass spectrometry and nuclear magnetic resonance spectroscopy, respectively. The loops were on average more flexible than the framework in both CH2 and m01, and the overall secondary structure was not affected by the additional disulfide bond. These data suggest that a human CH2 domain is relatively stable to unfolding at physiological temperature, and that both CH2 and the highly stable mutant m01 are promising new scaffolds for the development of therapeutics against human diseases.

Structures and functional implications of an AMP-binding cystathionine beta-synthase domain protein from a hyperthermophilic archaeon

Cystathionine beta-synthase domains are found in a myriad of proteins from organisms across the tree of life and have been hypothesized to function as regulatory modules that sense the energy charge of cells. Here we characterize the structure and stability of PAE2072, a dimeric tandem cystathionine beta-synthase domain protein from the hyperthermophilic crenarchaeon Pyrobaculum aerophilum. Crystal structures of the protein in unliganded and AMP-bound forms, determined at resolutions of 2.10 and 2.35 A, respectively, reveal remarkable conservation of key functional features seen in the gamma subunit of the eukaryotic AMP-activated protein kinase. The structures also confirm the presence of a suspected intermolecular disulfide bond between the two subunits that is shown to stabilize the protein. Our AMP-bound structure represents a first step in investigating the function of a large class of uncharacterized prokaryotic proteins. In addition, this work extends previous studies that have suggested that, in certain thermophilic microbes, disulfide bonds play a key role in stabilizing intracellular proteins and protein-protein complexes.