Amino acid neighbours and detailed conformational analysis of cysteines in proteins

Rational design of disulfide bonds to increase thermostability of Rhodococcus opacus catechol 1,2 dioxygenase

Catechol 1,2 dioxygenase is a versatile enzyme with several potential applications. However, due to its low thermostability, its industrial potential is not being met. In this study, the thermostability of a mesophilic catechol 1,2 dioxygenase from the species Rhodococcus opacus was enhanced via the introduction of disulphide bonds into its structure. Engineered designs (56) were obtained using computational prediction applications, with a set of hypothesized selection criteria narrowing the list to 9. Following recombinant production and purification, several of the designs demonstrated substantially improved protein thermostability. Notably, variant K96C-D278C yielded improvements including a 4.6°C increase in T50, a 725% increase in half-life, a 5.5°C increase in Tm, and a >10-fold increase in total turnover number compared to wild type. Stacking of best designs was not productive. Overall, current state-of-the-art prediction algorithms were effective for design of disulfide-thermostabilized catechol 1,2 dioxygenase.

Positions of cysteine residues reveal local clusters and hidden relationships to Sequons and Transmembrane domains in Human proteins

Membrane proteins often possess critical structural features, such as transmembrane domains (TMs), N-glycosylation, and disulfide bonds (SS bonds), which are essential to their structure and function. Here, we extend the study of the motifs carrying N-glycosylation, i.e. the sequons, and the Cys residues supporting the SS bonds, to the whole human proteome with a particular focus on the Cys positions in human proteins with respect to those of sequons and TMs. As the least abundant amino acid residue in protein sequences, the positions of Cys residues in proteins are not random but rather selected through evolution. We discovered that the frequency of Cys residues in proteins is length dependent, and the frequency of CC gaps formed between adjacent Cys residues can be used as a classifier to distinguish proteins with special structures and functions, such as keratin-associated proteins (KAPs), extracellular proteins with EGF-like domains, and nuclear proteins with zinc finger C2H2 domains. Most importantly, by comparing the positions of Cys residues to those of sequons and TMs, we discovered that these structural features can form dense clusters in highly repeated and mutually exclusive modalities in protein sequences. The evolutionary advantages of such complementarity among the three structural features are discussed, particularly in light of structural dynamics in proteins that are lacking from computational predictions. The discoveries made here highlight the sequence-structure-function axis in biological organisms that can be utilized in future protein engineering toward synthetic biology.

Single Deletion Unmasks Hidden Anti-Gram-Negative Bacterial Activity of an Insect Defensin-Derived Peptide

Insect defensins are a large family of antimicrobial peptides primarily active against Gram-positive bacteria. Here, we explore their hidden anti-Gram-negative bacterial potential via a nature-guided strategy inspired by natural deletion variants of Drosophila defensins. Referring to these variants, we deleted the equivalent region of an insect defensin with the first cysteine-containing N-terminus, and the last three cysteine-containing C-terminal regions remained. This 15-mer peptide exhibits low solubility and specifically targets Gram-positive bacteria. Further deletion of alanine-9 remarkably improves its solubility, unmasks its hidden anti-Gram-negative bacterial activity, and alters its states in different environments. Intriguingly, compared with the oxidized form, the 14-mer reduced peptide shows increased activity on Gram-positive and Gram-negative bacteria through a membrane-disruptive mechanism. The broad-spectrum activity and tolerance to high-salt environments and human serum, together with no toxicity to mammalian or human cells, make it a promising candidate for the design of new peptide antibiotics against Gram-negative bacterial infections.

Lessons on protein structure from interleukin-4: All disulfides are not created equal

Interleukin-4 (IL-4) is a hematopoietic cytokine composed by a four-helix bundle stabilized by an antiparallel beta-sheet and three disulfide bonds: Cys3-Cys127, Cys24-Cys65, and Cys46-Cys99. IL-4 is involved in several immune responses associated to infection, allergy, autoimmunity, and cancer. Besides its physiological relevance, IL-4 is often used as a "model" for protein design and engineering. Hence, to understand the role of each disulfide in the structure and dynamics of IL-4, we carried out several spectroscopic analyses (circular dichroism [CD], fluorescence, nuclear magnetic resonance [NMR]), and molecular dynamics (MD) simulations on wild-type IL-4 and four IL-4 disulfide mutants. All disulfide mutants showed loss of structure, altered interhelical angles, and looser core packings, showing that all disulfides are relevant for maintaining the overall fold and stability of the four-helix bundle motif, even at very low pH. In the absence of the disulfide connecting both protein termini Cys3-Cys127, C3T-IL4 showed a less packed protein core, loss of secondary structure (~9%) and fast motions on the sub-nanosecond time scale (lower S2 order parameters and larger τc correlation time), especially at the two protein termini, loops, beginning of helix A and end of helix D. In the absence of Cys24-Cys65, C24T-IL4 presented shorter alpha-helices (14% loss in helical content), altered interhelical angles, less propensity to form the small anti-parallel beta-sheet and increased dynamics. Simultaneously deprived of two disulfides (Cys3-Cys127 and Cys24-Cys65), IL-4 formed a partially folded "molten globule" with high 8-anilino-1-naphtalenesulphonic acid-binding affinity and considerable loss of secondary structure (~50%decrease), as shown by the far UV-CD, NMR, and MD data.

A computational approach to design a polyvalent vaccine against human respiratory syncytial virus

Human Respiratory Syncytial Virus (RSV) is one of the leading causes of lower respiratory tract infections (LRTI), responsible for infecting people from all age groups-a majority of which comprises infants and children. Primarily, severe RSV infections are accountable for multitudes of deaths worldwide, predominantly of children, every year. Despite several efforts to develop a vaccine against RSV as a potential countermeasure, there has been no approved or licensed vaccine available yet, to control the RSV infection effectively. Therefore, through the utilization of immunoinformatics tools, a computational approach was taken in this study, to design a multi-epitope polyvalent vaccine against two major antigenic subtypes of RSV, RSV-A and RSV-B. Potential predictions of the T-cell and B-cell epitopes were followed by extensive tests of antigenicity, allergenicity, toxicity, conservancy, homology to human proteome, transmembrane topology, and cytokine-inducing ability. The peptide vaccine was modeled, refined, and validated. Molecular docking analysis with specific Toll-like receptors (TLRs) revealed excellent interactions with suitable global binding energies. Additionally, molecular dynamics (MD) simulation ensured the stability of the docking interactions between the vaccine and TLRs. Mechanistic approaches to imitate and predict the potential immune response generated by the administration of vaccines were determined through immune simulations. Subsequent mass production of the vaccine peptide was evaluated; however, there remains a necessity for further in vitro and in vivo experiments to validate its efficacy against RSV infections.

TRAP1 S-nitrosylation as a model of population-shift mechanism to study the effects of nitric oxide on redox-sensitive oncoproteins

S-nitrosylation is a post-translational modification in which nitric oxide (NO) binds to the thiol group of cysteine, generating an S-nitrosothiol (SNO) adduct. S-nitrosylation has different physiological roles, and its alteration has also been linked to a growing list of pathologies, including cancer. SNO can affect the function and stability of different proteins, such as the mitochondrial chaperone TRAP1. Interestingly, the SNO site (C501) of TRAP1 is in the proximity of another cysteine (C527). This feature suggests that the S-nitrosylated C501 could engage in a disulfide bridge with C527 in TRAP1, resembling the well-known ability of S-nitrosylated cysteines to resolve in disulfide bridge with vicinal cysteines. We used enhanced sampling simulations and in-vitro biochemical assays to address the structural mechanisms induced by TRAP1 S-nitrosylation. We showed that the SNO site induces conformational changes in the proximal cysteine and favors conformations suitable for disulfide bridge formation. We explored 4172 known S-nitrosylated proteins using high-throughput structural analyses. Furthermore, we used a coarse-grained model for 44 protein targets to account for protein flexibility. This resulted in the identification of up to 1248 proximal cysteines, which could sense the redox state of the SNO site, opening new perspectives on the biological effects of redox switches. In addition, we devised two bioinformatic workflows ( https://github.com/ELELAB/SNO_investigation_pipelines ) to identify proximal or vicinal cysteines for a SNO site with accompanying structural annotations. Finally, we analyzed mutations in tumor suppressors or oncogenes in connection with the conformational switch induced by S-nitrosylation. We classified the variants as neutral, stabilizing, or destabilizing for the propensity to be S-nitrosylated and undergo the population-shift mechanism. The methods applied here provide a comprehensive toolkit for future high-throughput studies of new protein candidates, variant classification, and a rich data source for the research community in the NO field.

Designing Formulation Strategies for Enhanced Stability of Therapeutic Peptides in Aqueous Solutions: A Review

Over the past few decades, there has been a tremendous increase in the utilization of therapeutic peptides. Therapeutic peptides are usually administered via the parenteral route, requiring an aqueous formulation. Unfortunately, peptides are often unstable in aqueous solutions, affecting stability and bioactivity. Although a stable and dry formulation for reconstitution might be designed, from a pharmaco-economic and practical convenience point of view, a peptide formulation in an aqueous liquid form is preferred. Designing formulation strategies that optimize peptide stability may improve bioavailability and increase therapeutic efficacy. This literature review provides an overview of various degradation pathways and formulation strategies to stabilize therapeutic peptides in aqueous solutions. First, we introduce the major peptide stability issues in liquid formulations and the degradation mechanisms. Then, we present a variety of known strategies to inhibit or slow down peptide degradation. Overall, the most practical approaches to peptide stabilization are pH optimization and selecting the appropriate type of buffer. Other practical strategies to reduce peptide degradation rates in solution are the application of co-solvency, air exclusion, viscosity enhancement, PEGylation, and using polyol excipients.

Computational designing of a novel subunit vaccine for human cytomegalovirus by employing the immunoinformatics framework

Human cytomegalovirus (HCMV) is a widespread virus that can cause serious and irreversible neurological damage in newborns and even death in children who do not have the access to much-needed medications. While some vaccines and drugs are found to be effective against HCMV, their extended use has given rise to dose-limiting toxicities and the development of drug-resistant mutants among patients. Despite half a century's worth of research, the lack of a licensed HCMV vaccine heightens the need to develop newer antiviral therapies and vaccine candidates with improved effectiveness and reduced side effects. In this study, the immunoinformatics approach was utilized to design a potential polyvalent epitope-based vaccine effective against the four virulent strains of HCMV. The vaccine was constructed using seven CD8+ cytotoxic T lymphocytes epitopes, nine CD4+ helper T lymphocyte epitopes, and twelve linear B-cell lymphocyte epitopes that were predicted to be antigenic, non-allergenic, non-toxic, fully conserved, and non-human homologous. Subsequently, molecular docking study, protein-protein interaction analysis, molecular dynamics simulation (including the root mean square fluctuation (RMSF) and root mean square deviation (RMSD)), and immune simulation study rendered promising results assuring the vaccine to be stable, safe, and effective. Finally, in silico cloning was conducted to develop an efficient mass production strategy of the vaccine. However, further in vitro and in vivo research studies on the proposed vaccine are required to confirm its safety and efficacy.Communicated by Ramaswamy H. Sarma.

Functionally annotating cysteine disulfides and metal binding sites in the plant kingdom using AlphaFold2 predicted structures

Deep learning algorithms such as AlphaFold2 predict three-dimensional protein structure with high confidence. The recent release of more than 200 million structural models provides an unprecedented resource for functional protein annotation. Here, we used AlphaFold2 predicted structures of fifteen plant proteomes to functionally and evolutionary analyze cysteine residues in the plant kingdom. In addition to identification of metal ligands coordinated by cysteine residues, we systematically analyzed cysteine disulfides present in these structural predictions. Our analysis demonstrates most of these predicted disulfides are trustworthy due their high agreement (∼96%) with those present in X-ray and NMR protein structures, their characteristic disulfide stereochemistry, the biased subcellular distribution of their proteins and a higher degree of oxidation of their respective cysteines as measured by proteomics. Adopting an evolutionary perspective, zinc binding sites are increasingly present at the expense of iron-sulfur clusters in plants. Interestingly, disulfide formation is increased in secreted proteins of land plants, likely promoting sequence evolution to adapt to changing environments encountered by plants. In summary, Alphafold2 predicted structural models are a rich source of information for studying the role of cysteines residues in proteins of interest and for protein redox biology in general.

Immunoinformatics-guided designing and in silico analysis of epitope-based polyvalent vaccines against multiple strains of human coronavirus (HCoV)

OBJECTIVES

The group of human coronaviruses (HCoVs) consists of some highly pathogenic viruses that have caused several outbreaks in the past. The newly emerged strain of HCoV, the SARS-CoV-2 is responsible for the recent global pandemic that has already caused the death of hundreds of thousands of people due to the lack of effective therapeutic options.

METHODS

In this study, immunoinformatics methods were used to design epitope-based polyvalent vaccines which are expected to be effective against four different pathogenic strains of HCoV i.e., HCoV-OC43, HCoV-SARS, HCoV-MERS, and SARS-CoV-2.

RESULTS

The constructed vaccines consist of highly antigenic, non-allergenic, nontoxic, conserved, and non-homologous T-cell and B-cell epitopes from all the four viral strains. Therefore, they should be able to provide strong protection against all these strains. Protein-protein docking was performed to predict the best vaccine construct. Later, the MD simulation and immune simulation of the best vaccine construct also predicted satisfactory results. Finally, in silico cloning was performed to develop a mass production strategy of the vaccine.

CONCLUSION

If satisfactory results are achieved in further in vivo and in vitro studies, then the vaccines designed in this study might be effective as preventative measures against the selected HCoV strains.

Immunoinformatic Design of a Multivalent Peptide Vaccine Against Mucormycosis: Targeting FTR1 Protein of Major Causative Fungi

Mucormycosis is a potentially fatal illness that arises in immunocompromised people due to diabetic ketoacidosis, neutropenia, organ transplantation, and elevated serum levels of accessible iron. The sudden spread of mucormycosis in COVID-19 patients engendered massive concern worldwide. Comorbidities including diabetes, cancer, steroid-based medications, long-term ventilation, and increased ferritin serum concentration in COVID-19 patients trigger favorable fungi growth that in turn effectuate mucormycosis. The necessity of FTR1 gene-encoded ferrous permease for host iron acquisition by fungi has been found in different studies recently. Thus, targeting the transit component could be a potential solution. Unfortunately, no appropriate antifungal vaccine has been constructed as of yet. To date, mucormycosis has been treated with antiviral therapy and surgical treatment only. Thus, in this study, the FTR1 protein has been targeted to design a convenient and novel epitope-based vaccine with the help of immunoinformatics against four different virulent fungal species. Furthermore, the vaccine was constructed using 8 CTL, 2 HTL, and 1 LBL epitopes that were found to be highly antigenic, non-allergenic, non-toxic, and fully conserved among the fungi under consideration. The vaccine has very reassuring stability due to its high pI value of 9.97, conclusive of a basic range. The vaccine was then subjected to molecular docking, molecular dynamics, and immune simulation studies to confirm the biological environment’s safety, efficacy, and stability. The vaccine constructs were found to be safe in addition to being effective. Finally, we used in-silico cloning to develop an effective strategy for vaccine mass production. The designed vaccine will be a potential therapeutic not only to control mucormycosis in COVID-19 patients but also be effective in general mucormycosis events. However, further in vitro, and in vivo testing is needed to confirm the vaccine’s safety and efficacy in controlling fungal infections. If successful, this vaccine could provide a low-cost and effective method of preventing the spread of mucormycosis worldwide.

The use of bromopyridazinedione derivatives in chemical biology

Tools that facilitate the chemical modification of peptides and proteins are gaining an increasing amount of interest across many avenues of chemical biology as they enable a plethora of therapeutic, imaging and diagnostic applications. Cysteine residues and disulfide bonds have been highlighted as appealing targets for modification due to the highly homogenous nature of the products that can be formed through their site-selective modification. Amongst the reagents available for the site-selective modification of cysteine(s)/disulfide(s), pyridazinediones (PDs) have played a particularly important and enabling role. In this review, we outline the unique chemical features that make PDs especially well-suited to cysteine/disulfide modification on a wide variety of proteins and peptides, as well as provide context as to the problems solved (and applications enabled) by this technology.

The Effect of Length and Structure of Attached Polyethylene Glycol Chain on Hydrodynamic Radius, and Separation of PEGylated Human Serum Albumin by Chromatography

Purpose: This study focuses on the effect of length and structure of attached polyethylene glycol (PEG) chain on hydrodynamic radius (R_h ) and chromatographic retention of PEGylated protein. To this aim human serum albumin (HSA) as a standard protein was PEGylated site specifically with mPEG-maleimide. Methods: Separated PEG_HSA fractions were analyzed by size exclusion and anion exchange chromatography (AExC). The purity of fractions and the relative mobility of PEGylated and native proteins were analyzed by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). Hydrodynamic radius was determined based on the retention time of fractions on size exclusion chromatography (SEC), and also according to the previously reported equations. Results: A linear relation was shown between the molecular weight of attached PEG and R_h of PEGylated HSA. No significant difference between R_h of proteins modified with linear and branched PEG was shown. In SDS-PAGE, the delaying effect of branched PEG on movement of PEGylated protein was higher than that of linear PEG. Conclusion: As PEGylated HSA and dimer HSA have almost the same size and in SEC they elute at very close retention times, so in this case ion exchange chromatography (IExC) is more effective than SEC in separation of PEGylated HSA. Branched PEG- HSA showed earlier elution on anion exchange chromatography compared to linear PEG-HSA, that this can explain the different shielding effect of various structures of attached PEGs. The smaller size of PEGylated HSA in compare to the sum of the hydrodynamic radiuses of native HSA and attached PEG could be as a result of shielded attachment of polymer around protein.

Conformational dynamics in peptide toxins: Implications for receptor interactions and molecular design

Peptide toxins are potent and often exquisitely selective probes of the structure and function of ion channels and receptors, and are therefore of significant interest to the pharmaceutical and biotech industries as both pharmacological tools and therapeutic leads. The three-dimensional structures of peptide toxins are essential as a basis for understanding their structure-activity relationships and their binding to target receptors, as well as in guiding the design of analogues with modified potency and/or selectivity for key targets. NMR spectroscopy has played a key role in elucidating the structures of peptide toxins and probing their structure-function relationships. In this article, we highlight the additional important contribution of NMR to characterising the dynamics of peptide toxins. We also compare the information available from NMR measurements with that afforded by molecular dynamics simulations. We describe several examples of the importance of dynamics measurements over a range of timescales for understanding the structure-function relationships of peptide toxins and their receptor engagement. Peptide toxins that inhibit the voltage-gated potassium channel K_V1.3 with pM affinities display different degrees of conformational flexibility, even though they contain multiple disulfide bonds, and this flexibility can affect the relative orientation of residues that have been shown to be critical for channel binding. Information on the dynamic properties of peptide toxins is important in the design of analogues or mimetics where receptor-bound structures are not available.

3D interaction homology: Hydropathic interaction environments of serine and cysteine are strikingly different and their roles adapt in membrane proteins

Atomic-resolution protein structural models are prerequisites for many downstream activities like structure-function studies or structure-based drug discovery. Unfortunately, this data is often unavailable for some of the most interesting and therapeutically important proteins. Thus, computational tools for building native-like structural models from less-than-ideal experimental data are needed. To this end, interaction homology exploits the character, strength and loci of the sets of interactions that define a structure. Each residue type has its own limited set of backbone angle-dependent interaction motifs, as defined by their environments. In this work, we characterize the interactions of serine, cysteine and S-bridged cysteine in terms of 3D hydropathic environment maps. As a result, we explore several intriguing questions. Are the environments different between the isosteric serine and cysteine residues? Do some environments promote the formation of cystine S-S bonds? With the increasing availability of structural data for water-insoluble membrane proteins, are there environmental differences for these residues between soluble and membrane proteins? The environments surrounding serine and cysteine residues are dramatically different: serine residues are about 50% solvent exposed, while cysteines are only 10% exposed; the latter are more involved in hydrophobic interactions although there are backbone angle-dependent differences. Our analysis suggests that one driving force for -S-S- bond formation is a rather substantial increase in burial and hydrophobic interactions in cystines. Serine and cysteine become less and more, respectively, solvent-exposed in membrane proteins. 3D hydropathic environment maps are an evolving structure analysis tool showing promise as elements in a new protein structure prediction paradigm.

Stress-dependent conformational changes of artemin: Effects of heat and oxidant

Artemin is an abundant thermostable protein in Artemia embryos and it is considered as a highly efficient molecular chaperone against extreme environmental stress conditions. The conformational dynamics of artemin have been suggested to play a critical role in its biological functions. In this study, we have investigated the conformational and functional changes of artemin under heat and oxidative stresses to identify the relationship between its structure and function. The tertiary and quaternary structures of artemin were evaluated by fluorescence measurements, protein cross-linking analysis, and dynamic light scattering. Based on the structural analysis, artemin showed irreversible substantial conformational lability in responses to heat and oxidant, which was mainly mediated through the hydrophobic interactions and dimerization of the chaperone. In addition, the chaperone-like activity of heated and oxidized artemin was examined using lysozyme refolding assay and the results showed that although both factors, i.e. heat and oxidant, at specific levels improved artemin potency, simultaneous incubation with both stressors significantly triggered the chaperone activation. Moreover, the heat-induced dimerization of artemin was found to be the most critical factor for its activation. It was suggested that oxidation presumably acts through stabilizing the dimer structures of artemin through formation of disulfide bridges between the subunits and strengthens its chaperoning efficacy. Accordingly, it is proposed that artemin probably exists in a monomer–oligomer equilibrium in Artemia cysts and environmental stresses and intracellular portion of protein substrates may shift the equilibrium towards the active dimer forms of the chaperone.

Engineering a novel subunit vaccine against SARS-CoV-2 by exploring immunoinformatics approach

As the number of infections and deaths caused by the recent COVID-19 pandemic is increasing dramatically day-by-day, scientists are rushing towards developing possible countermeasures to fight the deadly virus, SARS-CoV-2. Although many efforts have already been put forward for developing potential vaccines; however, most of them are proved to possess negative consequences. Therefore, in this study, immunoinformatics methods were exploited to design a novel epitope-based subunit vaccine against the SARS-CoV-2, targeting four essential proteins of the virus i.e., spike glycoprotein, nucleocapsid phosphoprotein, membrane glycoprotein, and envelope protein. The highly antigenic, non-allergenic, non-toxic, non-human homolog, and 100% conserved (across other isolates from different regions of the world) epitopes were used for constructing the vaccine. In total, fourteen CTL epitopes and eighteen HTL epitopes were used to construct the vaccine. Thereafter, several in silico validations i.e., the molecular docking, molecular dynamics simulation (including the RMSF and RMSD studies), and immune simulation studies were also performed which predicted that the designed vaccine should be quite safe, effective, and stable within the biological environment. Finally, in silico cloning and codon adaptation studies were also conducted to design an effective mass production strategy of the vaccine. However, more in vitro and in vivo studies are required on the predicted vaccine to finally validate its safety and efficacy.

Blueprint of epitope-based multivalent and multipathogenic vaccines: targeted against the dengue and zika viruses

Both dengue virus (DENV) and zika virus (ZIKV) belong to the highly infectious Flaviviridae family that has already caused several outbreaks and epidemics in many countries. DENV and ZIKV cause two of the most wide spread mosquito-borne viral diseases in the world, dengue fever (DENF) and zika fever (ZIKF), respectively. In many regions around the world, both of these diseases can outbreak together and can be lethal as well as life-threatening. Unfortunately, there is no functional and satisfactory vaccine available to combat these viruses. Therefore, in this study, we have attempted to design a blue print of potential multivalent and multipathogenic vaccines using immunoinformatics approach, which can combat both the DENV and ZIKV infections, simultaneously. Initially, three vaccines were designed; containing highly antigenic, non-allergenic, and non-toxic epitopes of T-cell (100% conserved) and B-cell from all the four DENV serotypes and ZIKV. In total, nine cytotoxic T-lymphocytic (CTL), nine helper T-lymphocytic (HTL), and seven B-cell lymphocytic (BCL) epitopes were used to construct three vaccines using three different adjuvants, designated as 'V1', 'V2', and 'V3'. Later, V3 was found to be the best vaccine construct, determined by molecular docking analysis. Thereafter, several in silico validation studies including molecular dynamics simulation and immune simulation were performed which indicated that V3 might be quite stable and should generate substantial immune response in the biological environment. However, further in vivo and in vitro validation might be required to finally confirm the safety and efficacy of our suggested vaccine constructs.Communicated by Ramaswamy H. Sarma.

Exploiting the reverse vaccinology approach to design novel subunit vaccines against Ebola virus

Ebola virus is a highly pathogenic RNA virus that causes the Ebola haemorrhagic fever in human. This virus is considered as one of the dangerous viruses in the world with very high mortality rate. To date, no epitope-based subunit vaccine has yet been discovered to fight against Ebola although the outbreaks of this deadly virus took many lives in the past. In this study, approaches of reverse vaccinology were utilized in combination with different tools of immunoinformatics to design subunit vaccines against Ebola virus strain Mayinga-76. Three potential antigenic proteins of this virus i.e., matrix protein VP40, envelope glycoprotein and nucleoprotein were selected to construct the subunit vaccine. The MHC class-I, MHC class-II and B-cell epitopes were determined initially and after some robust analysis i.e., antigenicity, allergenicity, toxicity, conservancy and molecular docking study, EV-1, EV-2 and EV-3 were constructed as three potential vaccine constructs. These vaccine constructs are also expected to be effective on few other strains of Ebola virus since the highly conserved epitopes were used for vaccine construction. Thereafter, molecular docking study was conducted on these vaccines and EV-1 emerged as the best vaccine construct. Afterward, molecular dynamics simulation study revealed the good performances and stability of the intended vaccine protein. Finally, codon adaptation and in silico cloning were carried out to design a possible plasmid (pET-19b plasmid vector was used) for large scale production of the EV-1 vaccine. However, further in vitro and in vivo studies might be required on the predicted vaccines for final validation.

Involvement of the α-helical and Src homology 3 domains in the molecular assembly and enzymatic activity of human α1,6-fucosyltransferase, FUT8

BACKGROUND

Previous structural analyses showed that human α1,6-fucosyltransferase, FUT8 contains a catalytic domain along with two additional domains, N-terminal α-helical domain and C-terminal Src homology 3 domain, but these domains are unique to FUT8 among glycosyltransferases. The role that these domains play in formation of the active form of FUT8 has not been investigated. This study reports on attempts to determine the involvement of these domains in the functions of FUT8.

METHODS

Based on molecular modeling, the domain mutants were constructed by truncation and site-directed mutagenesis, and were heterologously expressed in Sf21 or COS-1 cells. The mutants were analyzed by SDS-PAGE and assayed for enzymatic activity. In vivo cross-linking experiments by introducing disulfide bonds were also carried out to examine the orientation of the domains in the molecular assembly.

RESULTS

Mutagenesis and molecular modeling findings suggest that human FUT8 potentially forms homodimer in vivo via intermolecular hydrophobic interactions involving α-helical domains. Truncation or site-directed mutagenesis findings indicated that α-helical and SH3 domains are all required for enzymatic activity. In addition, in vivo cross-linking experiments clearly indicated that the SH3 domain located in close proximity to the α-helical domain in an intermolecular manner.

CONCLUSIONS

α-Helical and SH3 domains are required for a fully active enzyme, and are also involved in homophilic dimerization, which probably results in the formation of the active form of human FUT8.

GENERAL SIGNIFICANCE

α-Helical and SH3 domains, which are not commonly found in glycosyltransferases, play roles in the formation of the functional quaternary structure of human FUT8.

The biophysics of superoxide dismutase-1 and amyotrophic lateral sclerosis

Few proteins have come under such intense scrutiny as superoxide dismutase-1 (SOD1). For almost a century, scientists have dissected its form, function and then later its malfunction in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). We now know SOD1 is a zinc and copper metalloenzyme that clears superoxide as part of our antioxidant defence and respiratory regulation systems. The possibility of reduced structural integrity was suggested by the first crystal structures of human SOD1 even before deleterious mutations in the sod1 gene were linked to the ALS. This concept evolved in the intervening years as an impressive array of biophysical studies examined the characteristics of mutant SOD1 in great detail. We now recognise how ALS-related mutations perturb the SOD1 maturation processes, reduce its ability to fold and reduce its thermal stability and half-life. Mutant SOD1 is therefore predisposed to monomerisation, non-canonical self-interactions, the formation of small misfolded oligomers and ultimately accumulation in the tell-tale insoluble inclusions found within the neurons of ALS patients. We have also seen that several post-translational modifications could push wild-type SOD1 down this toxic pathway. Recently we have come to view ALS as a prion-like disease where both the symptoms, and indeed SOD1 misfolding itself, are transmitted to neighbouring cells. This raises the possibility of intervention after the initial disease presentation. Several small-molecule and biologic-based strategies have been devised which directly target the SOD1 molecule to change the behaviour thought to be responsible for ALS. Here we provide a comprehensive review of the many biophysical advances that sculpted our view of SOD1 biology and the recent work that aims to apply this knowledge for therapeutic outcomes in ALS.

The unusual conformation of cross-strand disulfide bonds is critical to the stability of β-hairpin peptides

The cross-strand disulfides (CSDs) found in β-hairpin antimicrobial peptides (β-AMPs) show a unique disulfide geometry that is characterized by unusual torsion angles and a short Cα-Cα distance. While the sequence and disulfide bond connectivity of disulfide-rich peptides is well studied, much less is known about the disulfide geometry found in CSDs and their role in the stability of β-AMPs. To address this, we solved the nuclear magnetic resonance (NMR) structure of the β-AMP gomesin (Gm) at 278, 298, and 310 K, examined the disulfide bond geometry of over 800 disulfide-rich peptides, and carried out extensive molecular dynamics (MD) simulation of the peptides Gm and protegrin. The NMR data suggests Cα-Cα distances characteristic for CSDs are independent of temperature. Analysis of disulfide-rich peptides from the Protein Data Bank revealed that right-handed and left-handed rotamers are equally likely in CSDs. The previously reported preference for right-handed rotamers was likely biased by restricting the analysis to peptides and proteins solved using X-ray crystallography. Furthermore, data from MD simulations showed that the short Cα-Cα distance is critical for the stability of these peptides. The unique disulfide geometry of CSDs poses a challenge to biomolecular force fields and to retain the stability of β-hairpin fold over long simulation times, restraints on the torsion angles might be required.

Rational Design of Alginate Lyase from Microbulbifer sp. Q7 to Improve Thermal Stability

Alginate lyase degrades alginate by the β-elimination mechanism to produce oligosaccharides with special bioactivities. The low thermal stability of alginate lyase limits its industrial application. In this study, introducing the disulfide bonds while using the rational design methodology enhanced the thermal stability of alginate lyase cAlyM from Microbulbifer sp. Q7. Enzyme catalytic sites, secondary structure, spatial configuration, and molecular dynamic simulation were comprehensively analyzed. When compared with cAlyM, the mutants D102C-A300C and G103C-T113C showed an increase by 2.25 and 1.16 h, respectively, in half-life time at 45 °C, in addition to increases by 1.7 °C and 0.4 °C in the melting temperature, respectively. The enzyme-specific activity and k_cat/K_m values of D102C-A300C were 1.8- and 1.5-times higher than those of cAlyM, respectively. The rational design strategy that was used in this study provides a valuable method for improving the thermal stability of the alginate lyase.

Association of BUD13-ZNF259-APOA5-APOA1-SIK3 cluster polymorphism in 11q23.3 and structure of APOA5 with increased plasma triglyceride levels in a Korean population

In this association study on chromosome 11, the data from 12,537 Korean individuals within the Health Examinee (HEXA) and the Korea Association Resource (KARE) projects were analysed to identify genetic loci correlating with increased and decreased plasma triglyceride (TG) levels. We identified a locus in chromosomal region 11q23.3 that harbours genes BUD13, ZNF259, APOA5, APOA1, and SIK3, which may be associated with plasma TG levels. In this locus, 13 relevant single-nucleotide polymorphisms (SNPs) were found: rs184616707, rs118175510, rs60954647, rs79408961, and rs180373 (near BUD13); rs11604424 (in ZNF259); rs2075291, rs651821, and rs7123666 (in or near APOA5); rs525028 (near APOA1), and rs645258, rs10160754, and rs142395187 (in or near SIK3). All 13 SNPs satisfied the genome-wide significance level (P < 5.0 × 10^-8) in both meta-analysis and conditional analysis. Haplotype analysis of six SNPs (rs79408961, rs180373, rs2075291, rs651821, rs525028, and rs10160754) that were selected based on the β coefficient and conditional P values, revealed nine common haplotypes (with frequency 0.02-0.34) associated with both increased and reduced TG levels. Furthermore, to shed light on possible structural implications, we modelled and simulated the G185C variant of APOA5 (corresponding to rs2075291), which showed the strongest association. Molecular dynamics simulation results showed that this polymorphic variant of APOA5 has a different hydrogen bond network, increased average distance between chains, and an ability to form distinct clusters. Owing to the orientation of cysteine, the possibility of disulphide bond formation with other proteins is evident. In summary, our association and modelling analyses provided evidence that genetic variations in chromosomal region 11q23.3 are associated with elevated TG levels.

Allostery in Coagulation Factor VIIa Revealed by Ensemble Refinement of Crystallographic Structures

A critical step in injury-induced initiation of blood coagulation is the formation of the complex between the trypsin-like protease coagulation factor VIIa (FVIIa) and its cofactor tissue factor (TF), which converts FVIIa from an intrinsically poor enzyme to an active protease capable of activating zymogens of downstream coagulation proteases. Unlike its constitutively active ancestor trypsin, FVIIa is allosterically activated (by TF). Here, ensemble refinement of crystallographic structures, which uses multiple copies of the entire structure as a means of representing structural flexibility, is applied to explore the impacts of inhibitor binding to trypsin and FVIIa, as well as cofactor binding to FVIIa. To assess the conformational flexibility and its role in allosteric pathways in these proteases, main-chain hydrogen bond networks are analyzed by calculating the hydrogen-bond propensity. Mapping pairwise propensity differences between relevant structures shows that binding of the inhibitor benzamidine to trypsin has a minor influence on the protease flexibility. For FVIIa, in contrast, the protease domain is "locked" into the catalytically competent trypsin-like configuration upon benzamidine binding as indicated by the stabilization of key structural features: the nonprime binding cleft and the oxyanion hole are stabilized, and the effect propagates from the active site region to the calcium-binding site and to the vicinity of the disulphide bridge connecting with the light chain. TF binding to FVIIa furthermore results in stabilization of the 170 loop, which in turn propagates an allosteric signal from the TF-binding region to the active site. Analyses of disulphide bridge energy and flexibility reflect the striking stability difference between the unregulated enzyme and the allosterically activated form after inhibitor or cofactor binding. The ensemble refinement analyses show directly, for the first time to our knowledge, whole-domain structural footprints of TF-induced allosteric networks present in x-ray crystallographic structures of FVIIa, which previously only have been hypothesized or indirectly inferred.

Detection of miRNA cancer biomarkers using light activated Molecular Beacons

Early detection of cancer biomarkers can reduce cancer mortality rate. miRNAs are small non-coding RNAs whose expression changes upon the onset of various types of cancer. Biosensors that specifically detect such biomarkers can be engineered and integrated into point-of-care devices (POC) using label-free detection, high sensibility and compactness. In this paper, a new engineered Molecular Beacon (MB) construct used to detect miRNAs is presented. Such a construct is immobilized onto biosensor surfaces in a covalent and spatially oriented way using the photonic technology Light Assisted Molecular Immobilization (LAMI). The construct consists of a Cy3 labelled MB covalently attached to a light-switchable peptide. One MB construct contains a poly-A sequence in its loop region while the other contains a sequence complementary to the cancer biomarker miRNA-21. The constructs have been characterized by UV-Vis spectroscopy, mass spectrometry and HPLC. LAMI led to the successful immobilization of the engineered constructs onto thiol functionalized optically flat quartz slides and Silicon on Insulator (SOI) sensor surfaces. The immobilized Cy3 labelled MB construct has been imaged using confocal fluorescence microscopy (CFM). The bioavailability of the immobilized engineered MB biosensors was confirmed through specific hybridization with the Cy5 labelled complementary sequence and imaged by CFM and FRET. Hybridization kinetics have been monitored using steady state fluorescence spectroscopy. The label-free detection of miRNA-21 was also achieved by using integrated photonic sensing structures. The engineered light sensitive constructs can be immobilized onto thiol reactive surfaces and are currently being integrated in a POC device for the detection of cancer biomarkers.

Photonic based detection strategies of cancer miRNA biomarkers after Light Assisted Molecular Immobilization (LAMI) of peptide-MB biosensor constructs.

Functional protein nanostructures: a chemical toolbox

Nature has evolved an optimal synthetic factory in the form of translational and posttranslational processes by which millions of proteins with defined primary sequences and 3D structures can be built. Nature's toolkit gives rise to protein building blocks, which dictates their spatial arrangement to form functional protein nanostructures that serve a myriad of functions in cells, ranging from biocatalysis, formation of structural networks, and regulation of biochemical processes, to sensing. With the advent of chemical tools for site-selective protein modifications and recombinant engineering, there is a rapid development to develop and apply synthetic methods for creating structurally defined, functional protein nanostructures for a broad range of applications in the fields of catalysis, materials and biomedical sciences. In this review, design principles and structural features for achieving and characterizing functional protein nanostructures by synthetic approaches are summarized. The synthetic customization of protein building blocks, the design and introduction of recognition units and linkers and subsequent assembly into structurally defined protein architectures are discussed herein. Key examples of these supramolecular protein nanostructures, their unique functions and resultant impact for biomedical applications are highlighted.

Bioconjugation with Maleimides: A Useful Tool for Chemical Biology

Maleimide chemistry stands out in the bioconjugation toolbox by virtue of its synthetic accessibility, excellent reactivity, and practicability. The second-generation of clinically approved antibody-drug conjugates (ADC) and much of the current ADC pipeline in clinical trials contain the maleimide linkage. However, thiosuccinimide linkages are now known to be less robust than once thought, and ergo, are correlated with suboptimal pharmacodynamics, pharmacokinetics, and safety profiles in some ADC constructs. Rational design of novel generations of maleimides and maleimide-type reagents have been reported to address the shortcomings of classical maleimides, allowing for the formation of robust bioconjugate linkages. This review highlights the main strategies for rational reagent design that have allowed irreversible bioconjugations in cysteines, reversible labelling strategies and disulfide re-bridging.

Aromatic interactions in β-hairpin scaffold stability: A historical perspective

Non-covalent interactions between naturally occurring aromatic residues have been widely exploited as scaffold stabilizing agents in de novo designed peptides and in Nature - inspired structures. Our understanding of the factors driving aromatic interactions and their observed interaction geometries have advanced remarkably with improvements in conventional structural studies, availability of novel molecular methods and in silico studies, which have together provided atomistic information on aromatic interactions and interaction strengths. This review attempts to recapitulate the early advances in our understanding of aromatic interactions as stabilizing agents of peptide β-hairpins.

Structure of an H1-Bound 6-Nucleosome Array Reveals an Untwisted Two-Start Chromatin Fiber Conformation

Chromatin adopts a diversity of regular and irregular fiber structures in vitro and in vivo. However, how an array of nucleosomes folds into and switches between different fiber conformations is poorly understood. We report the 9.7 Å resolution crystal structure of a 6-nucleosome array bound to linker histone H1 determined under ionic conditions that favor incomplete chromatin condensation. The structure reveals a flat two-start helix with uniform nucleosomal stacking interfaces and a nucleosome packing density that is only half that of a twisted 30-nm fiber. Hydroxyl radical footprinting indicates that H1 binds the array in an on-dyad configuration resembling that observed for mononucleosomes. Biophysical, cryo-EM, and crosslinking data validate the crystal structure and reveal that a minor change in ionic environment shifts the conformational landscape to a more compact, twisted form. These findings provide insights into the structural plasticity of chromatin and suggest a possible assembly pathway for a 30-nm fiber.

Structure-function analyses reveal a triple β-turn receptor-bound conformation of adrenomedullin 2/intermedin and enable peptide antagonist design

The cardioprotective vasodilator peptide adrenomedullin 2/intermedin (AM2/IMD) and the related adrenomedullin (AM) and calcitonin gene-related peptide (CGRP) signal through three heterodimeric receptors comprising the calcitonin receptor-like class B G protein-coupled receptor (CLR) and a variable receptor activity-modifying protein (RAMP1, -2, or -3) that determines ligand selectivity. The CGRP receptor (RAMP1:CLR) favors CGRP binding, whereas the AM₁ (RAMP2:CLR) and AM₂ (RAMP3:CLR) receptors favor AM binding. How AM2/IMD binds the receptors and how RAMPs modulate its binding is unknown. Here, we show that AM2/IMD binds the three purified RAMP-CLR extracellular domain (ECD) complexes with a selectivity profile that is distinct from those of CGRP and AM. AM2/IMD bound all three ECD complexes but preferred the CGRP and AM₂ receptor complexes. A 2.05 Å resolution crystal structure of an AM2/IMD antagonist fragment-bound RAMP1-CLR ECD complex revealed that AM2/IMD binds the complex through a unique triple β-turn conformation that was confirmed by peptide and receptor mutagenesis. Comparisons of the receptor-bound conformations of AM2/IMD, AM, and a high-affinity CGRP analog revealed differences that may have implications for biased signaling. Guided by the structure, enhanced-affinity AM2/IMD antagonist variants were developed, including one that discriminates the AM₁ and AM₂ receptors with ∼40-fold difference in affinities and one stabilized by an intramolecular disulfide bond. These results reveal differences in how the three peptides engage the receptors, inform development of AM2/IMD-based pharmacological tools and therapeutics, and provide insights into RAMP modulation of receptor pharmacology.

Self-homodimerization of an actinoporin by disulfide bridging reveals implications for their structure and pore formation

The Trp111 to Cys mutant of sticholysin I, an actinoporin from Stichodactyla helianthus sea anemone, forms a homodimer via a disulfide bridge. The purified dimer is 193 times less hemolytic than the monomer. Ultracentrifugation, dynamic light scattering and size-exclusion chromatography demonstrate that monomers and dimers are the only independent oligomeric states encountered. Indeed, circular dichroism and fluorescence spectroscopies showed that Trp/Tyr residues participate in homodimerization and that the dimer is less thermostable than the monomer. A homodimer three-dimensional model was constructed and indicates that Trp147/Tyr137 are at the homodimer interface. Spectroscopy results validated the 3D-model and assigned 85° to the disulfide bridge dihedral angle responsible for dimerization. The homodimer model suggests that alterations in the membrane/carbohydrate-binding sites in one of the monomers, as result of dimerization, could explain the decrease in the homodimer ability to form pores.

A robust cosolvent-compatible halohydrin dehalogenase by computational library design

To improve the applicability of halohydrin dehalogenase as a catalyst for reactions in the presence of organic cosolvents, we explored a computational library design strategy (Framework for Rapid Enzyme Stabilization by Computational libraries) that involves discovery and in silico evaluation of stabilizing mutations. Energy calculations, disulfide bond predictions and molecular dynamics simulations identified 218 point mutations and 35 disulfide bonds with predicted stabilizing effects. Experiments confirmed 29 stabilizing point mutations, most of which were located in two distinct regions, whereas introduction of disulfide bonds was not effective. Combining the best mutations resulted in a 12-fold mutant (HheC-H12) with a 28°C higher apparent melting temperature and a remarkable increase in resistance to cosolvents. This variant also showed a higher optimum temperature for catalysis while activity at low temperature was preserved. Mutant H12 was used as a template for the introduction of mutations that enhance enantioselectivity or activity. Crystal structures showed that the structural changes in the H12 mutant mostly agreed with the computational predictions and that the enhanced stability was mainly due to mutations that redistributed surface charges and improved interactions between subunits, the latter including better interactions of water molecules at the subunit interfaces.

Enzymatic modification by point mutation and functional analysis of an omega-6 fatty acid desaturase from Arctic Chlamydomonas sp

Arctic Chlamydomonas sp. is a dominant microalgal strain in cold or frozen freshwater in the Arctic region. The full-length open reading frame of the omega-6 fatty acid desaturase gene (AChFAD6) was obtained from the transcriptomic database of Arctic Chlamydomonas sp. from the KOPRI culture collection of polar micro-organisms. Amino acid sequence analysis indicated the presence of three conserved histidine-rich segments as unique characteristics of omega-6 fatty acid desaturases, and three transmembrane regions transported to plastidic membranes by chloroplast transit peptides in the N-terminal region. The AChFAD6 desaturase activity was examined by expressing wild-type and V254A mutant (Mut-AChFAD6) heterologous recombinant proteins. Quantitative gas chromatography indicated that the concentration of linoleic acids in AChFAD6-transformed cells increased more than 3-fold [6.73 ± 0.13 mg g^-1 dry cell weight (DCW)] compared with cells transformed with vector alone. In contrast, transformation with Mut-AChFAD6 increased the concentration of oleic acid to 9.23 ± 0.18 mg g^-1 DCW, indicating a change in enzymatic activity to mimic that of stearoyl-CoA desaturase. These results demonstrate that AChFAD6 of Arctic Chlamydomonas sp. increases membrane fluidity by enhancing denaturation of C18 fatty acids and facilitates production of large quantities of linoleic fatty acids in prokaryotic expression systems.

Elimination of a Free Cysteine by Creation of a Disulfide Bond Increases the Activity and Stability of Candida boidinii Formate Dehydrogenase

NAD+-dependent formate dehydrogenase (FDH; EC 1.2.1.2) is an industrial enzyme widely used for NADH regeneration. However, enzyme inactivation caused by the oxidation of cysteine residues is a flaw of native FDH. In this study, we relieved the oxidation of the free cysteine of FDH from Candida boidinii (CboFDH) through the construction of disulfide bonds between A10 and C23 as well as I239 and C262. Variants A10C, I239C, and A10C/I239C were obtained by the site-directed mutagenesis and their properties were studied. Results showed that there were no significant changes in the optimum temperature and pH between variants and wild-type CboFDH. However, the stabilities of all variant enzymes were improved. Specifically, the CboFDH variant A10C (A10Cfdh) showed a significant increase in copper ion resistance and acid resistance, a 6.7-fold increase in half-life at 60°C, and a 1.4-fold increase in catalytic efficiency compared with the wild type. Asymmetric synthesis of l-tert-leucine indicated that the process time was reduced by 40% with variant A10Cfdh, which benefited from the increase in catalytic efficiency. Circular dichroism analysis and molecular dynamics simulation indicated that variants that contained disulfide bonds lowered the overall root mean square deviation (RMSD) and consequently increased the protein rigidity without affecting the secondary structure of enzyme. This work is expected to provide a viable strategy to avoid the microbial enzyme inactivation caused by the oxidation of the free cysteine residues and improving their performances.

IMPORTANCE

FDH is widely used for NADH regeneration in dehydrogenase-based synthesis of optically active compounds to decrease the cost of production. This study highlighted a viable strategy that was used to eliminate the oxidation of free cysteine residues of FDH from Candida boidinii by the introduction of disulfide bonds. Using this strategy, we obtained a variant FDH with improved activity and stability. The improvement of activity and stability of FDH is expected to reduce its price and then further to decrease the cost of its application.

Site-Selective Disulfide Modification of Proteins: Expanding Diversity beyond the Proteome

The synthetic transformation of polypeptides with molecular accuracy holds great promise for providing functional and structural diversity beyond the proteome. Consequently, the last decade has seen an exponential growth of site-directed chemistry to install additional features into peptides and proteins even inside living cells. The disulfide rebridging strategy has emerged as a powerful tool for site-selective modifications since most proteins contain disulfide bonds. In this Review, we present the chemical design, advantages and limitations of the disulfide rebridging reagents, while summarizing their relevance for synthetic customization of functional protein bioconjugates, as well as the resultant impact and advancement for biomedical applications.

Conservation and diversity in the ultralong third heavy-chain complementarity-determining region of bovine antibodies

A subset of bovine antibodies have an exceptionally long third heavy-chain complementarity determining region (CDR H3) that is highly variable in sequence and includes multiple cysteines. These long CDR H3s (up to 69 residues) fold into a long stalk atop which sits a knob domain that is located far from the antibody surface. Three new bovine Fab crystal structures have been determined to decipher the conserved and variable features of ultralong CDR H3s that lead to diversity in antigen recognition. Despite high sequence variability, the stalks adopt a conserved β-ribbon structure, while the knob regions share a conserved β-sheet that serves as a scaffold for two connecting loops of variable length and conformation, as well as one conserved disulfide. Variation in patterns and connectivity of the remaining disulfides contribute to the knob structural diversity. The unusual architecture of these ultralong bovine CDR H3s for generating diversity is unique in adaptive immune systems.

Rational Design of Disulfide Bonds Increases Thermostability of a Mesophilic 1,3-1,4-β-Glucanase from Bacillus terquilensis

1,3-1,4-β-glucanase is an important biocatalyst in brewing industry and animal feed industry, while its low thermostability often reduces its application performance. In this study, the thermostability of a mesophilic β-glucanase from Bacillus terquilensis was enhanced by rational design and engineering of disulfide bonds in the protein structure. Protein spatial configuration was analyzed to pre-exclude the residues pairs which negatively conflicted with the protein structure and ensure the contact of catalytic center. The changes in protein overall and local flexibility among the wild-type enzyme and the designated mutants were predicted to select the potential disulfide bonds for enhancement of thermostability. Two residue pairs (N31C-T187C and P102C-N125C) were chosen as engineering targets and both of them were proved to significantly enhance the protein thermostability. After combinational mutagenesis, the double mutant N31C-T187C/P102C-N125C showed a 48.3% increase in half-life value at 60°C and a 4.1°C rise in melting temperature (Tm) compared to wild-type enzyme. The catalytic property of N31C-T187C/P102C-N125C mutant was similar to that of wild-type enzyme. Interestingly, the optimal pH of double mutant was shifted from pH6.5 to pH6.0, which could also increase its industrial application. By comparison with mutants with single-Cys substitutions, the introduction of disulfide bonds and the induced new hydrogen bonds were proved to result in both local and overall rigidification and should be responsible for the improved thermostability. Therefore, the introduction of disulfide bonds for thermostability improvement could be rationally and highly-effectively designed by combination with spatial configuration analysis and molecular dynamics simulation.

Hyperstabilization of Tetrameric Bacillus sp. TB-90 Urate Oxidase by Introducing Disulfide Bonds through Structural Plasticity

Bacillus sp. TB-90 urate oxidase (BTUO) is one of the most thermostable homotetrameric enzymes. We previously reported [Hibi, T., et al. (2014) Biochemistry 53, 3879-3888] that specific binding of a sulfate anion induced thermostabilization of the enzyme, because the bound sulfate formed a salt bridge with two Arg298 residues, which stabilized the packing between two β-barrel dimers. To extensively characterize the sulfate-binding site, Arg298 was substituted with cysteine by site-directed mutagenesis. This substitution markedly increased the protein melting temperature by ∼ 20 °C compared with that of the wild-type enzyme, which was canceled by reduction with dithiothreitol. Calorimetric analysis of the thermal denaturation suggested that the hyperstabilization resulted from suppression of the dissociation of the tetramer into the two homodimers. The crystal structure of R298C at 2.05 Å resolution revealed distinct disulfide bond formation between the symmetrically related subunits via Cys298, although the Cβ distance between Arg298 residues of the wild-type enzyme (5.4 Å apart) was too large to predict stable formation of an engineered disulfide cross-link. Disulfide bonding was associated with local disordering of interface loop II (residues 277-300), which suggested that the structural plasticity of the loop allowed hyperstabilization by disulfide formation. Another conformational change in the C-terminal region led to intersubunit hydrogen bonding between Arg7 and Asp312, which probably promoted mutant thermostability. Knowledge of the disulfide linkage of flexible loops at the subunit interface will help in the development of new strategies for enhancing the thermostabilization of multimeric proteins.

Metal-Directed Design of Supramolecular Protein Assemblies

Owing to their central roles in cellular signaling, construction, and biochemistry, protein-protein interactions (PPIs) and protein self-assembly have become a major focus of molecular design and synthetic biology. In order to circumvent the complexity of constructing extensive noncovalent interfaces, which are typically involved in natural PPIs and protein self-assembly, we have developed two design strategies, metal-directed protein self-assembly (MDPSA) and metal-templated interface redesign (MeTIR). These strategies, inspired by both the proposed evolutionary roles of metals and their prevalence in natural PPIs, take advantage of the favorable properties of metal coordination (bonding strength, directionality, and reversibility) to guide protein self-assembly with minimal design and engineering. Using a small, monomeric protein (cytochrome cb562) as a model building block, we employed MDPSA and MeTIR to create a diverse array of functional supramolecular architectures which range from structurally tunable oligomers to metalloprotein complexes that can properly self-assemble in living cells into novel metalloenzymes. The design principles and strategies outlined herein should be readily applicable to other protein systems with the goal of creating new PPIs and protein assemblies with structures and functions not yet produced by natural evolution.

Cysteine-Rich Atrial Secretory Protein from the Snail Achatina achatina: Purification and Structural Characterization

Despite extensive studies of cardiac bioactive peptides and their functions in molluscs, soluble proteins expressed in the heart and secreted into the circulation have not yet been reported. In this study, we describe an 18.1-kDa, cysteine-rich atrial secretory protein (CRASP) isolated from the terrestrial snail Achatina achatina that has no detectable sequence similarity to any known protein or nucleotide sequence. CRASP is an acidic, 158-residue, N-glycosylated protein composed of eight alpha-helical segments stabilized with five disulphide bonds. A combination of fold recognition algorithms and ab initio folding predicted that CRASP adopts an all-alpha, right-handed superhelical fold. CRASP is most strongly expressed in the atrium in secretory atrial granular cells, and substantial amounts of CRASP are released from the heart upon nerve stimulation. CRASP is detected in the haemolymph of intact animals at nanomolar concentrations. CRASP is the first secretory protein expressed in molluscan atrium to be reported. We propose that CRASP is an example of a taxonomically restricted gene that might be responsible for adaptations specific for terrestrial pulmonates.

Engineering Streptavidin and a Streptavidin-Binding Peptide with Infinite Binding Affinity and Reversible Binding Capability: Purification of a Tagged Recombinant Protein to High Purity via Affinity-Driven Thiol Coupling

To extend and improve the utility of the streptavidin-binding peptide tag (SBP-tag) in applications ranging from affinity purification to the reversible immobilization of recombinant proteins, a cysteine residue was introduced to the streptavidin mutein SAVSBPM18 and the SBP-tag to generate SAVSBPM32 and SBP(A18C), respectively. This pair of derivatives is capable of forming a disulfide bond through the newly introduced cysteine residues. SAVSBPM32 binds SBP-tag and biotin with binding affinities (K_d ~ 10^-8M) that are similar to SAVSBPM18. Although SBP(A18C) binds to SAVSBPM32 more weakly than SBP-tag, the binding affinity is sufficient to bring the two binding partners together efficiently before they are locked together via disulfide bond formation–a phenomenon we have named affinity-driven thiol coupling. Under the condition with SBP(A18C) tags in excess, two SBP(A18C) tags can be captured by a tetrameric SAVSBPM32. The stoichiometry of the disulfide-bonded SAVSBPM32-SBP(A18C) complex was determined using a novel two-dimensional electrophoresis method which has general applications for analyzing the composition of disulfide-bonded protein complexes. To illustrate the application of this reversible immobilization technology, optimized conditions were established to use the SAVSBPM32-affinity matrix for the purification of a SBP(A18C)-tagged reporter protein to high purity. Furthermore, we show that the SAVSBPM32-affinity matrix can also be applied to purify a biotinylated protein and a reporter protein tagged with the unmodified SBP-tag. The dual (covalent and non-covalent) binding modes possible in this system offer great flexibility to many different applications which need reversible immobilization capability.

Developments and recent advancements in the field of endogenous amino acid selective bond forming reactions for bioconjugation

Bioconjugation methodologies have proven to play a central enabling role in the recent development of biotherapeutics and chemical biology approaches. Recent endeavours in these fields shed light on unprecedented chemical challenges to attain bioselectivity, biocompatibility, and biostability required by modern applications. In this review the current developments in various techniques of selective bond forming reactions of proteins and peptides were highlighted. The utility of each endogenous amino acid-selective conjugation methodology in the fields of biology and protein science has been surveyed with emphasis on the most relevant among reported transformations; selectivity and practical use have been discussed.