Native and modeled disulfide bonds in proteins: knowledge-based approaches toward structure prediction of disulfide-rich polypeptides

DSDBASE 2.0: updated version of DiSulphide dataBASE, a database on disulphide bonds in proteins

Disulphide bonds are stabilizing crosslinks in proteins and serve to enhance their thermal stability. In proteins that are small and rich in disulphide bonds, they could be the major determining factor for the choice of conformational state since their constraints on appropriate backbone conformation can be substantial. Such crosslinks and their positional conservation could itself enable protein family and functional association. Despite the importance of the field, there is no comprehensive database on disulphide crosslinks that is available to the public. Herein we provide information on disulphides in DSDBASE2.0, an updated and significantly expanded database that is freely available, fully annotated and manually curated database on native and modelled disulphides. The web interface also provides several useful computational tools that have been specifically developed for proteins containing disulphide crosslinks. The modelling of disulphide crosslinks is performed using stereochemical criteria, coded within our Modelling of Disulphides in Proteins (MODIP) algorithm. The inclusion of modelled disulphides potentially enhances the loop database substantially, thereby permitting the recognition of compatible polypeptide segments that could serve as templates for immediate modelling. The DSDBASE2.0 database has been updated to include 153,944 PDB entries, 216,096 native and 20,153,850 modelled disulphide bond segments from PDB January 2021 release. The current database also provides a resource to user-friendly search for multiple disulphide bond containing loops, along with annotation of their function using GO and subcellular localization of the query. Furthermore, it is possible to obtain the three-dimensional models of disulphide-rich small proteins using an independent algorithm, RANMOD, that generates and examines random, but allowed backbone conformations of the polypeptide. DSDBASE2.0 still remains the largest open-access repository that organizes all disulphide bonds of proteins on a single platform. The database can be accessed from http://caps.ncbs.res.in/dsdbase2.

Enhancement of Pharmaceutical Urate Oxidase Thermostability by Rational Design of De Novo Disulfide Bridge

Background and Purpose

As a therapeutic enzyme, urate oxidase is utilized in the reduction of uric acid in various conditions such as gout or tumor syndrome lysis. However, even bearing kinetical advantage over other counterparts, it suffers from structural instability most likely due to its subcellular and fungal origin.

Objectives

In this research, by using rational design and introduction of de novo disulfide bridge in urate oxidase structure, we designed and created a thermostable urate oxidase for the first time.

Materials and Methods

Utilizing site-directed mutagenesis and only with one point mutation we constructed two separate mutants: Ala6Cys and Ser282Cys which covalently linked subunits of enzyme each other. Single mutation to cysteine created three inter-chain disulfide bridges and one hydrogen bond in Ala6Cys and two disulfide bridges in Ser282Cys.

Results

Both mutants showed 10 °C increase in optimum activity compared to wild-type enzyme while the K_m values for both increased by 50% and their specific activity compromised. The thermal stability of Ser282Cys increased remarkably by comparing Ala6Cys and wild-type enzymes. Estimation of half life for wild-type enzyme demonstrated 38.5 min, while for Ala6Cys and Ser282Cys were 138 and 115 min, respectively. Interestingly, the optimal pH of both mutants was broaden from 7 to 10, which could make them candidates for industrial applications.

Conclusion

It seemed that introducing disulfide bridges resulted in local and overall rigidity by bringing two adjacent sites of enzyme together and decreasing the conformational entropy of unfolding state is responsible for the enhancement of thermostability.

Recombinant Production of a Novel Fusion Protein: Listeriolysin O Fragment Fused to S1 Subunit of Pertussis Toxin

Background:

Some resources have suggested that genetically inactivated PTs bear a more protective effect than chemically inactivated products. This study aimed to produce new version of PT, by cloning an inactive PTS1 in a fusion form with N-terminal half of the LLO pore-forming toxin.

Methods:

Deposited pdb structure file of the PT was used to model an extra disulfide bond. Codon-optimized ORF of the PTS1 was used to make recombinant constructs of PTS1 and LLO-PTS1 in the pPSG-IBA35 vector. The recombinant PTS1 and LLO-PTS1 proteins were expressed in BL21 DE3 and SHuffle T7 strains of E. coli and purified by affinity chromatography. Cytotoxic effects of the recombinant proteins were examined in the MCF-7 cell line.

Results:

The purity of the products proved to be more than 85%, and the efficiency of the disulfide bond formation in SHuffle T7 strain was higher than BL21 DE3 strain. No cytotoxicity of the recombinant proteins was observed in MCF-7 cells. Soluble recombinant PTS1 and LLO-PTS1 proteins were produced in SHuffle T7 strain of E. coli with high efficiency of disulfide bonds formation.

Conclusion:

The LLO-PTS1 with corrected disulfide bonds was successfully expressed in E. coli SHuffleT7 strain. Due to the safety for human cells, this chimeric molecule can be an option to prevent pertussis disease if its immunostimulatory effects would be confirmed in the future.

Structure-based discovery of antiviral inhibitors targeting the E dimer interface of Japanese encephalitis virus

Flaviviruses are emerging arthropod-borne viruses posing a great threat to human beings worldwide. The E dimer configuration of the flavivirus was prominent during viral assembly, maturation and entry. Neutralization antibodies targeting E dimer played the important role in controlling the flavivirus infection. Previously, the ideal drug target of small molecular inhibitors of JEV was viral proteases and polymerases. The crystal structure of JEV E protein showed a conserved pocket in it is important at membrane fusion step. Recently, a set of anti-virus drugs has been found by virtual screening. Here, we show that the fusion-loop pocket of JEV E protein was a conservative region and an ideal drug target. ChemDiv-3 from virtual screening as the lead compound was found to show a relatively modest inhibition effect for JEV in vitro and in vivo test and could interfere with the formation of JEV sE dimer. ChemDiv-3 interacts with the amino acid residues ASN 313, PRO 314, ALA 315, and VAL 323 in E protein via hydrogen bonds for occupation of the fusion-loop pocket. The key binding sites LYS 312, ALA 513 and THR 317 forming the fusion-loop pocket are the same and other auxiliary sites are similar among the flavivirus. Taken together, the fusion-loop pocket of the flavivirus could be one promising target for drug discovery.

Covalently linked dengue virus envelope glycoprotein dimers reduce exposure of the immunodominant fusion loop epitope

A problem in the search for an efficient vaccine against dengue virus is the immunodominance of the fusion loop epitope (FLE), a segment of the envelope protein E that is buried at the interface of the E dimers coating mature viral particles. Anti-FLE antibodies are broadly cross-reactive but poorly neutralizing, displaying a strong infection enhancing potential. FLE exposure takes place via dynamic ‘breathing' of E dimers at the virion surface. In contrast, antibodies targeting the E dimer epitope (EDE), readily exposed at the E dimer interface over the region of the conserved fusion loop, are very potent and broadly neutralizing. We here engineer E dimers locked by inter-subunit disulfide bonds, and show by X-ray crystallography and by binding to a panel of human antibodies that these engineered dimers do not expose the FLE, while retaining the EDE exposure. These locked dimers are strong immunogen candidates for a next-generation vaccine.

The immunodominant epitope of dengue virus envelope protein (E) induces poorly neutralizing antibodies, which poses a problem for vaccine development. Here, the authors engineer covalently locked E dimers exposing an epitope that has been shown to induce potent and broadly neutralizing antibodies.

Engineering and kinetic stabilization of the therapeutic enzyme Anabeana variabilis phenylalanine ammonia lyase

Anabeana variabilis phenylalanine ammonia lyase has just recently been discovered and introduced in clinical trials of phenylketonuria enzyme replacement therapy for its outstanding kinetic properties. In the present study, kinetic stabilization of this therapeutically important enzyme has been explored by introduction of a disulfide bond into the structure. Site-directed mutagenesis was performed with quick-change PCR method. Recombinant wild-type and mutated enzymes were expressed in Escherichia coli, and his-tagged proteins were affinity purified. Formation of disulfide bond was confirmed by Ellman's method, and then chemical unfolding, kinetic behavior, and thermal inactivation of mutated enzyme were compared with the wild type. Based on our results, the Q292C mutation resulted in a significant improvement in kinetic stability and resistance against chemical unfolding of the enzyme while kinetic parameters and pH profile of enzyme activity were remained unaffected. The results of the present study provided an insight towards designing phenylalanine ammonia lyases with higher stability.

Comparative genomics of odorant binding proteins in Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus

About 1 million people in the world die each year from diseases spread by mosquitoes, and understanding the mechanism of host identification by the mosquitoes through olfaction is at stake. The role of odorant binding proteins (OBPs) in the primary molecular events of olfaction in mosquitoes is becoming an important focus of biological research in this area. Here, we present a comprehensive comparative genomics study of OBPs in the three disease-transmitting mosquito species Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus starting with the identification of 110 new OBPs in these three genomes. We have characterized their genomic distribution and orthologous and phylogenetic relationships. The diversity and expansion observed with respect to the Aedes and Culex genomes suggests that the OBP gene family acquired functional diversity concurrently with functional constraints posed on these two species. Sequences with unique features have been characterized such as the "two-domain OBPs" (previously known as Atypical OBPs) and "MinusC OBPs" in mosquito genomes. The extensive comparative genomics featured in this work hence provides useful primary insights into the role of OBPs in the molecular adaptations of mosquito olfactory system and could provide more clues for the identification of potential targets for insect repellants and attractants.

Analysis of factors that induce cysteine bonding state

Regarding the fact that the protein structure is principally encoded in its sequence, investigating the bonding state of cysteine has gained a great deal of attention due to its significance in the formation of protein structure. Due to lack of evident influence of free cysteines on the protein structure, it may be expected that only half-cystines convey encoded information. The results obtained from the analysis of amino acid distribution in proximity of both states of cysteines explicitly indicated that perquisite information for inducing cysteine bonding state is present even in the flanking amino acid sequences of free cysteines.

Structure and folding of disulfide-rich miniproteins: insights from molecular dynamics simulations and MM-PBSA free energy calculations

The fold of small disulfide-rich proteins largely relies on two or more disulfide bridges that are main components of the hydrophobic core. Because of the small size of these proteins and their high cystine content, the cysteine connectivity has been difficult to ascertain in some cases, leading to uncertainties and debates in the literature. Here, we use molecular dynamics simulations and MM-PBSA free energy calculations to compare similar folds with different disulfide pairings in two disulfide-rich miniprotein families, namely the knottins and the short-chain scorpion toxins, for which the connectivity has been discussed. We first show that the MM-PBSA approach is able to discriminate the correct knotted topology of knottins from the laddered one. Interestingly, a comparison of the free energy components for kalata B1 and MCoTI-II suggests that cyclotides and squash inhibitors, although sharing the same scaffold, are stabilized through different interactions. Application to short-chain scorpion toxins suggests that the conventional cysteine pairing found in many homologous toxins is significantly more stable than the unconventional pairing reported for maurotoxin and for spinoxin. This would mean that native maurotoxin and spinoxin are not at the lowest free energy minimum and might result from kinetically rather than thermodynamically driven oxidative folding processes. For both knottins and toxins, the correct or conventional disulfide connectivities provide lower flexibilities and smaller deviations from the initial conformations. Overall, our work suggests that molecular dynamics simulations and the MM-PBSA approach to estimate free energies are useful tools to analyze and compare disulfide bridge connectivities in miniproteins.

Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background

Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue.

Results

Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues.

Conclusion

We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Predicting disulfide connectivity from protein sequence using multiple sequence feature vectors and secondary structure

MOTIVATION

Disulfide bonds are primary covalent crosslinks between two cysteine residues in proteins that play critical roles in stabilizing the protein structures and are commonly found in extracy-toplasmatic or secreted proteins. In protein folding prediction, the localization of disulfide bonds can greatly reduce the search in conformational space. Therefore, there is a great need to develop computational methods capable of accurately predicting disulfide connectivity patterns in proteins that could have potentially important applications.

RESULTS

We have developed a novel method to predict disulfide connectivity patterns from protein primary sequence, using a support vector regression (SVR) approach based on multiple sequence feature vectors and predicted secondary structure by the PSIPRED program. The results indicate that our method could achieve a prediction accuracy of 74.4% and 77.9%, respectively, when averaged on proteins with two to five disulfide bridges using 4-fold cross-validation, measured on the protein and cysteine pair on a well-defined non-homologous dataset. We assessed the effects of different sequence encoding schemes on the prediction performance of disulfide connectivity. It has been shown that the sequence encoding scheme based on multiple sequence feature vectors coupled with predicted secondary structure can significantly improve the prediction accuracy, thus enabling our method to outperform most of other currently available predictors. Our work provides a complementary approach to the current algorithms that should be useful in computationally assigning disulfide connectivity patterns and helps in the annotation of protein sequences generated by large-scale whole-genome projects.

AVAILABILITY

The prediction web server and Supplementary Material are accessible at http://foo.maths.uq.edu.au/~huber/disulfide

HARMONY: a server for the assessment of protein structures

Protein structure validation is an important step in computational modeling and structure determination. Stereochemical assessment of protein structures examine internal parameters such as bond lengths and Ramachandran (varphi,psi) angles. Gross structure prediction methods such as inverse folding procedure and structure determination especially at low resolution can sometimes give rise to models that are incorrect due to assignment of misfolds or mistracing of electron density maps. Such errors are not reflected as strain in internal parameters. HARMONY is a procedure that examines the compatibility between the sequence and the structure of a protein by assigning scores to individual residues and their amino acid exchange patterns after considering their local environments. Local environments are described by the backbone conformation, solvent accessibility and hydrogen bonding patterns. We are now providing HARMONY through a web server such that users can submit their protein structure files and, if required, the alignment of homologous sequences. Scores are mapped on the structure for subsequent examination that is useful to also recognize regions of possible local errors in protein structures. HARMONY server is located at http://caps.ncbs.res.in/harmony/

Analysis of disulphide bond connectivity patterns in protein tertiary structure

The analysis of disulphide bond containing proteins in the Protein Data Bank (PDB) revealed that out of 27,209 protein structures analyzed, 12,832 proteins contain at least one intra-chain disulphide bond and 811 proteins contain at least one inter-chain disulphide bond. The intra-chain disulphide bond containing proteins can be grouped into 256 categories based on the number of disulphide bonds and the disulphide bond connectivity patterns (DBCPs) that were generated according to the position of half-cystine residues along the protein chain. The PDB entries corresponding to these 256 categories represent 509 unique SCOP superfamilies. A simple web-based computational tool is made freely available at the website that allows flexible queries to be made on the database in order to retrieve useful information on the disulphide bond containing proteins in the PDB. The database is useful to identify the different SCOP superfamilies associated with a particular disulphide bond connectivity pattern or vice versa. It is possible to define a query based either on a single field or a combination of the following fields, i.e., PDB code, protein name, SCOP superfamily name, number of disulphide bonds, disulphide bond connectivity pattern and the number of amino acid residues in a protein chain and retrieve information that match the criterion. Thereby, the database may be useful to select suitable protein structural templates in order to model the more distantly related protein homologs/analogs using the comparative modeling methods.

Construction of a 3D model of CP12, a protein linker

The chloroplast protein CP12 is known to play a leading role in a complex formation with the enzymes GAPDH and PRK. As a preliminary step towards the understanding of the complex formation mechanism and the exact role of this protein linker, a comparative modelling of the CP12 protein of the green alga Chlamydomonas reinhardtii was performed. Because of the very few structural information and poor template similarities, the derivation of the model consisted in an iterative trial-and-error procedure using the comparative modelling program MODELLER, the following three structure validation programs PROCHECK, PROSA, and WHATIF, and molecular mechanics energy refinement of the model using the program CHARMM. The analysis of the final model reveals a scaffold of key residues that is believed to be essential in the folding mechanism and that coincides with the residues conserved throughout the CP12 family. Our results suggest that this protein is a typical disordered protein. Finally, the various mechanisms by which the CP12 protein can self-interact or binds to other enzymes are discussed in light of its modelled structure and characteristics.

Limited tendency of α-helical residues to form disulfide bridges: a structural explanation

Disulfide bridges have an enormous impact on the structure of a large number of proteins and polypeptides. Understanding the structural basis that regulates their formation may be important for the design of novel peptide-based molecules with a specific fold and stability. Here we report a statistical analysis of the relationships between secondary structure and disulfide bond formation, carried out using a large database of protein structures. Our analyses confirm the observation sporadically reported in previous investigations that cysteine residues located in alpha-helices display a limited tendency to form disulfide bridges. The very low occurrence of the disulfide bond in all alpha-chains compared to all beta-chains indicates that this property is also evident when proteins with different topologies are investigated. Taking advantage of the large database that endorsed the analysis on relatively rare motifs, we demonstrate that cysteine residues embedded in 3(10) helices present a good tendency to form disulfide bonds. This result is somewhat surprising since 3(10) helices are commonly assimilated into alpha-helices. A plausible structural explanation for the observed data has been derived combining analyses of disulfide bond sequence separation and of the length of the different secondary structure elements.

Prediction of disulfide connectivity from protein sequences

The difficulties in predicting disulfide connectivity from protein sequences lie in the nonlocal properties of the disulfide bridges that involve cysteine pairs at large sequence separation. Though some progress has been recently made in the prediction of disulfide connectivity, the current methods predict less than half of the disulfide patterns for the data set sharing less than 30% sequence identity. In this report, we use the support vector machines based on sequence features such as the coupling between the local sequence environments of cysteine pair, the cysteines sequence separations, and the global sequence descriptor, such as amino acid content. Our approach is able to predict 55% of the disulfide patterns of proteins with two to five disulfide bridges, which is 11-26% higher than other methods in the literature.