Experimental evaluation of single‐domain antibodies predicted by molecular dynamics simulations to have elevated thermal stability

Recently Bekker et al. [Bekker G‐J et al. Protein Sci. 2019;28:429–438.] described a computational strategy of applying molecular‐dynamics simulations to estimate the relative stabilities of single‐domain antibodies, and utilized their method to design changes with the aim of increasing the stability of a single‐domain antibody with a known crystal structure. The structure from which they generated potentially stabilizing mutations is an anti‐cholera toxin single domain antibody selected from a naïve library which has relatively low thermal stability, reflected by a melting point of 48°C. Their work was purely theoretical, so to examine their predictions, we prepared the parental and predicted stabilizing mutant single domain antibodies and examined their thermal stability, ability to refold and affinity. We found that the mutation that improved stability the most (~7°C) was one which changed an amino acid in CDR1 from an asparagine to an aspartic acid. This change unfortunately was also accompanied by a reduction in affinity. Thus, while their modeling did appear to successfully predict stabilizing mutations, introducing mutations in the binding regions is problematic. Of further interest, the mutations selected via their high temperature simulations, did improve refolding, suggesting that they were successful in stabilizing the structure at high temperatures and thereby decrease aggregation. Our result should permit them to reassess and refine their model and may one day lead to a usefulin silico approach to protein stabilization.     

Thermal stability of single‐domain antibodies estimated by molecular dynamics simulations

Single‐domain antibodies (sdAbs) function like regular antibodies, however, consist of only one domain. Because of their low molecular weight, sdAbs have advantages with respect to production and delivery to their targets and for applications such as antibody drugs and biosensors. Thus, sdAbs with high thermal stability are required. In this work, we chose seven sdAbs, which have a wide range of melting temperature (T_m) values and known structures. We applied molecular dynamics (MD) simulations to estimate their relative stability and compared them with the experimental data. High‐temperature MD simulations at 400 K and 500 K were executed with simulations at 300 K as a control. The fraction of native atomic contacts, Q, measured for the 400 K simulations showed a fairly good correlation with the T_m values. Interestingly, when the residues were classified by their hydrophobicity and size, the Q values of hydrophilic residues exhibited an even better correlation, suggesting that stabilization is correlated with favorable interactions of hydrophilic residues. Measuring the Q value on a per‐residue level enabled us to identify residues that contribute significantly to the instability and thus demonstrating how our analysis can be used in a mutant case study.     

Analysis of the binding loops configuration and surface adaptation of different crystallized single‐domain antibodies in response to various antigens

Monoclonal antibodies have revolutionized the biomedical field through their ubiquitous utilization in different diagnostics and therapeutic applications. Despite this widespread use, their large size and structural complexity have limited their versatility in specific applications. The antibody variable region that is responsible for binding antigen is embodied within domains that can be rescued individually as single‐domain antibody (sdAb) fragments. Because of the unique characteristics of sdAbs, such as low molecular weight, high physicochemical stability, and the ability to bind antigens inaccessible to conventional antibodies, they represent a viable alternative to full‐length antibodies. Consequently, 149 crystal structures of sdAbs, originating from human (VH), camelids (VHH), or sharks (VNAR), were retrieved from the Protein Data Bank, and their structures were compared. The 3 types of sdAbs displayed complementarity determining regions (CDRs) with different lengths and configurations. CDR3 of the VHH and VNAR domains were dominated by pleated and extended orientations, respectively. Although VNAR showed the smallest average molecular weight and molecular surface area compared with VHH and VH antibodies. However, the solvent accessible surface area measurements of the 3 tested sdAbs types were very similar. All the antihapten VHH antibodies showed pleated CDR3, which were sufficient to create a binding pocket to accommodate haptens (methotrexate and azo dyes) in terms of shape and electrostatic potential. The sdAbs that recognized lysozyme showed more diversity in their CDR3 orientation to enable them to recognize various topographies of lysozyme. Subsequently, the three sdAb classes were different in size and surface area and have shown distinguishable ability to optimize their CDR length and orientation to recognize different antigen classes.     

Role of the non‐hypervariable FR3 D‐E loop in single‐domain antibody recognition of haptens and carbohydrates

Single‐domain antibodies (sdAbs), the variable domains of camelid heavy chain‐only antibodies, are generally thought to poorly recognize nonproteinaceous small molecules and carbohydrates in comparison with conventional antibodies. However, the structures of anti‐methotrexate, anti‐triclocarban and anti‐cortisol sdAbs revealed unexpected contributions of the non‐hypervariable “CDR4” loop, formed between β‐strands D and E of framework region 3, in binding. Here, we investigated the potential role of CDR4 in sdAb binding to a hapten, 15‐acetyl‐deoxynivalenol (15‐AcDON), and to carbohydrates. We constructed and panned a phage‐displayed library in which CDR4 of the 15‐AcDON‐specific sdAb, NAT‐267, was extended and randomized. From this library, we identified one sdAb, MA‐232, bearing a 14‐residue insertion in CDR4 and showing improved binding to 15‐AcDON by ELISA and surface plasmon resonance. On the basis of these results, we constructed a second set of phage‐displayed libraries in which the CDR4 and other regions of three hapten‐ or carbohydrate‐binding sdAbs were diversified. With the goal of identifying sdAbs with novel glycan‐binding specificities, we panned the library against four tumor‐associated carbohydrate antigens but were unable to enrich binding phages. Thus, we conclude that while CDR4 may play a role in binding of some rare hapten‐specific sdAbs, diversifying this region through molecular engineering is probably not a general solution to sdAb carbohydrate recognition in the absence of a paired V_L domain.     

Role of single disulfide linkages in the folding and activity of scyllatoxin‐based BH3 domain mimetics

Anti‐apoptotic Bcl‐2 proteins are implicated in pathogenic cell survival and have attracted considerable interest as therapeutic targets. We recently developed a class of synthetic peptide based on scyllatoxin (ScTx) designed to mimic the helical BH3 interaction domain of the pro‐apoptotic Bcl‐2 protein Bax. In this communication, the contribution of single disulfides in the folding and function of ScTx‐Bax peptides was investigated. We synthesized five ScTx‐Bax variants, each presenting a different combination of native disulfide linkage and evaluated their ability to directly bind Bcl‐2 in vitro. It was determined that the position of the disulfide linkage had significant implications on the structure and function of ScTx‐Bax peptides. This study underscores the importance of structural dynamics in BH3:Bcl‐2 interactions and further validates ScTx‐based ligands as potential modulators of anti‐apoptotic Bcl‐2 function. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.     

Unique single‐domain antigen binding fragments derived from naturally occurring camel heavy‐chain antibodies

The humoral immune response of camels, dromedaries and llamas includes functional antibodies formed by two heavy chains and no light chains. The amino acid sequence of the variable domain of the naturally occurring heavy‐chain antibodies reveals the necessary adaptations to compensate for the absence of the light chain. In contrast to the conventional antibodies, a large proportion of the heavy‐chain antibodies acts as competitive enzyme inhibitors. Studies on the dromedary immunoglobulin genes start to shed light on the ontogeny of these heavy‐chain antibodies. The presence of the heavy‐chain antibodies and the possibility of immunizing a dromedary allows for the production of antigen binders consisting of a single domain only. These minimal antigen‐binding fragments are well expressed in bacteria, bind the antigen with affinity in the nM range and are very stable. We expect that such camelid single domain antibodies will find their way into a number of biotechnological or medical applications. The structure of the camelid single domain is homologous to the human VH, however, the antigen‐binding loop structures deviate fundamentally from the canonical structures described for human or mouse VHs. This has two additional advantages: (1) the camel or llama derived single domain antibodies might be an ideal scaffold for anti‐idiotypic vaccinations; and (2) the development of smaller peptides or peptide mimetic drugs derived from of the antigen binding loops might be facilitated due to their less complex antigen binding site. Copyright © 1999 John Wiley & Sons, Ltd.     

Prediction of VH–VL domain orientation for antibody variable domain modeling

The antigen‐binding site of antibodies forms at the interface of their two variable domains, VH and VL, making VH–VL domain orientation a factor that codetermines antibody specificity and affinity. Preserving VH–VL domain orientation in the process of antibody engineering is important in order to retain the original antibody properties, and predicting the correct VH–VL orientation has also been recognized as an important factor in antibody homology modeling. In this article, we present a fast sequence‐based predictor that predicts VH–VL domain orientation with Q² values ranging from 0.54 to 0.73 on the evaluation set. We describe VH–VL orientation in terms of the six absolute ABangle parameters that have recently been proposed as a means to separate the different degrees of freedom of VH–VL domain orientation. In order to assess the impact of adjusting VH–VL orientation according to our predictions, we use the set of antibody structures of the recently published Antibody Modeling Assessment (AMA) II study. In comparison to the original AMAII homology models, we find an improvement in the accuracy of VH–VL orientation modeling, which also translates into an improvement in the average root‐mean‐square deviation with regard to the crystal structures. Proteins 2015; 83:681–695. © 2015 Wiley Periodicals, Inc.     

Trapping a 96 degrees domain rotation in two distinct conformations by engineered disulfide bridges

Engineering disulfide bridges is a common technique to lock a protein movement in a defined conformational state. We have designed two double mutants of Escherichia coli 5'-nucleotidase to trap the enzyme in both an open (S228C, P513C) and a closed (P90C, L424C) conformation by the formation of disulfide bridges. The mutant proteins have been expressed, purified, and crystallized, to structurally characterize the designed variants. The S228C, P513C is a double mutant crystallized in two different crystal forms with three independent conformers, which differ from each other by a rotation of up to 12 degrees of the C-terminal domain with respect to the N-terminal domain. This finding, as well as an analysis of the domain motion in the crystal, indicates that the enzyme still exhibits considerable residual domain flexibility. In the double mutant that was designed to trap the enzyme in the closed conformation, the structure analysis reveals an unexpected intermediate conformation along the 96 degrees rotation trajectory between the open and closed enzyme forms. A comparison of the five independent conformers analyzed in this study shows that the domain movement of the variant enzymes is characterized by a sliding movement of the residues of the domain interface along the interface, which is in contrast to a classical closure motion where the residues of the domain interface move perpendicular to the interface.     

Chemical synthesis and evaluation of a backbone‐cyclized minimized 2‐helix Z‐domain

The Z‐molecule is a small, engineered IgG‐binding affinity protein derived from the immunoglobulin‐binding domain B of Staphylococcus aureus protein A. The Z‐domain consists of 58 amino acids forming a well‐defined antiparallel three‐helix structure. Two of the three helices are involved in ligand binding, whereas the third helix provides structural support to the three‐helix bundle. The small size and the stable three‐helix structure are two attractive properties comprised in the Z‐domain, but a further reduction in size of the protein is valuable for several reasons. Reduction in size facilitates synthetic production of any protein‐based molecule, which is beneficial from an economical viewpoint. In addition, a smaller protein is easier to manipulate through chemical modifications. By omitting the third stabilizing helix from the Z‐domain and joining the N‐ and C‐termini by a native peptide bond, the affinity protein obtains the advantageous properties of a smaller scaffold and in addition becomes resistant to exoproteases. We here demonstrate the synthesis and evaluation of a novel cyclic two‐helix Z‐domain. The molecule has retained affinity for its target protein, is resistant to heat treatment, and lacks both N‐ and C‐termini. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Methods to identify the substrates of thiol‐disulfide oxidoreductases

The formation of a disulfide bond is a critical step in the folding of numerous secretory and membrane proteins and catalyzed in vivo. A variety of mechanisms and protein structures have evolved to catalyze oxidative protein folding. Those enzymes that directly interact with a folding protein to accelerate its oxidative folding are mostly thiol‐disulfide oxidoreductases that belong to the thioredoxin superfamily. The enzymes of this class often use a CXXC active‐site motif embedded in their thioredoxin‐like fold to promote formation, isomerization, and reduction of a disulfide bond in their target proteins. Over the past decade or so, an increasing number of substrates of the thiol‐disulfide oxidoreductases that are present in the ER of mammalian cells have been discovered, revealing that the enzymes play unexpectedly diverse physiological functions. However, functions of some of these enzymes still remain unclear due to the lack of information on their substrates. Here, we review the methods used by researchers to identify the substrates of these enzymes and provide data that show the importance of using trichloroacetic acid in sample preparation for the substrate identification, hoping to aid future studies. We particularly focus on successful studies that have uncovered physiological substrates and functions of the enzymes in the periplasm of Gram‐negative bacteria and the endoplasmic reticulum of mammalian cells. Similar approaches should be applicable to enzymes in other cellular compartments or in other organisms.     

Kinetically trapped metastable intermediate of a disulfide‐deficient mutant of the starch‐binding domain of glucoamylase

Refolding of a thermally unfolded disulfide‐deficient mutant of the starch‐binding domain of glucoamylase was investigated using differential scanning calorimetry, isothermal titration calorimetry, CD, and ¹H NMR. When the protein solution was rapidly cooled from a higher temperature, a kinetic intermediate was formed during refolding. The intermediate was unexpectedly stable compared with typical folding intermediates that have short half‐lives. It was shown that this intermediate contained substantial secondary structure and tertiary packing and had the same binding ability with β‐cyclodextrin as the native state, suggesting that the intermediate is highly‐ordered and native‐like on the whole. These characteristics differ from those of partially folded intermediates such as molten globule states. Far‐UV CD spectra showed that the secondary structure was once disrupted during the transition from the intermediate to the native state. These results suggest that the intermediate could be an off‐pathway type, possibly a misfolded state, that has to undergo unfolding on its way to the native state.     

Enhanced transport of plant‐produced rabies single‐chain antibody‐RVG peptide fusion protein across an in cellulo blood–brain barrier device

The biomedical applications of antibody engineering are developing rapidly and have been expanded to plant expression platforms. In this study, we have generated a novel antibody molecule in planta for targeted delivery across the blood–brain barrier (BBB). Rabies virus (RABV) is a neurotropic virus for which there is no effective treatment after entry into the central nervous system. This study investigated the use of a RABV glycoprotein peptide sequence to assist delivery of a rabies neutralizing single‐chain antibody (ScFv) across an in cellulo model of human BBB. The 29 amino acid rabies virus peptide (RVG) recognizes the nicotinic acetylcholine receptor (nAchR) at neuromuscular junctions and the BBB. ScFv and ScFv‐RVG fusion proteins were produced in Nicotiana benthamiana by transient expression. Both molecules were successfully expressed and purified, but the ScFv expression level was significantly higher than that of ScFv‐RVG fusion. Both ScFv and ScFv‐RVG fusion molecules had potent neutralization activity against RABVin cellulo. The ScFv‐RVG fusion demonstrated increased binding to nAchR and entry into neuronal cells, compared to ScFv alone. Additionally, a human brain endothelial cell line BBB model was used to demonstrate that plant‐produced ScFv‐RVG^P fusion could translocate across the cells. This study indicates that the plant‐produced ScFv‐RVG^P fusion protein was able to cross the in celluloBBB and neutralize RABV.     

Discovery of an intermolecular disulfide bond required for the thermostability of a heterodimeric protein from the thermophile Hydrogenobacter thermophilus

Factors that increase protein thermostability are of considerable interest in both scientific and industrial fields. Disulfide bonds are one of such factors that increase thermostability, but are rarely found in intracellular proteins because of the reducing environment of the cytosol. Here, we report the first example of an intermolecular disulfide bond between heteromeric subunits of a novel-type phosphoserine phosphatase from a thermophilic bacterium Hydrogenobacter thermophilus, which contributes to the protein thermostability at the physiological temperature. Comparison of remaining soluble proteins between wild-type and cysteine-deleted mutant using SDS-PAGE revealed that the disulfide bond increases the thermostability of the whole protein by tightly connecting a subunit with low solubility to the partner with higher solubility. Furthermore, it was strongly suggested that the disulfide bond is formed and contributes to the stability in vivo. This finding will open new avenues for the design of proteins with increased thermostability.     

A phage‐displayed single‐chain Fab library optimized for rapid production of single‐chain IgGs

Phage‐displayed synthetic antibody (Ab) repertoires have become a major source of affinity reagents for basic and clinical research. Specific Abs identified from such libraries are often screened as fragments antigen binding (Fabs) produced in bacteria, and those with desired biochemical characteristics are reformatted for production as full‐length immunoglobulin G (IgG) in mammalian cells. The conversion of Fabs to IgGs is a cumbersome and often rate‐limiting step in the development of Abs. Moreover, biochemical properties required for lead IgG development are not always shared by the Fabs, and these issues are not uncovered until a significant effort has been spent on Abs that ultimately will not be useful. Thus, there is a need for simple and rapid techniques to convert phage‐displayed Fabs to IgGs at an early stage of the Ab screening process. We report the generation of a highly diverse phage‐displayed synthetic single‐chain Fab (scFab) library, in which the light and heavy chains were tethered with an optimized linker. Following selection, pools of scFabs were converted to single‐chain IgGs (scIgGs) en masse, enabling facile screening of hundreds of phage‐derived scIgGs. We show that this approach can be used to rapidly screen for and select scIgGs that target cell‐surface receptors, and scIgGs behave the same as conventional IgGs.     

X‐ray crystallographic validation of structure predictions used in computational design for protein stabilization

Protein engineering aimed at enhancing enzyme stability is increasingly supported by computational methods for calculation of mutant folding energies and for the design of disulfide bonds. To examine the accuracy of mutant structure predictions underlying these computational methods, crystal structures of thermostable limonene epoxide hydrolase variants obtained by computational library design were determined. Four different predicted effects indeed contributed to the obtained stabilization: (i) enhanced interactions between a flexible loop close to the N‐terminus and the rest of the protein; (ii) improved interactions at the dimer interface; (iii) removal of unsatisfied hydrogen bonding groups; and (iv) introduction of additional positively charged groups at the surface. The structures of an eightfold and an elevenfold mutant showed that most mutations introduced the intended stabilizing interactions, and side‐chain conformations were correctly predicted for 72–88% of the point mutations. However, mutations that introduced a disulfide bond in a flexible region had a larger influence on the backbone conformation than predicted. The enzyme active sites were unaltered, in agreement with the observed preservation of catalytic activities. The structures also revealed how a c‐Myc tag, which was introduced for facile detection and purification, can reduce access to the active site and thereby lower the catalytic activity. Finally, sequence analysis showed that comprehensive mutant energy calculations discovered stabilizing mutations that are not proposed by the consensus or B‐FIT methods. Proteins 2015; 83:940–951. © 2015 Wiley Periodicals, Inc.     

Microfabricated arrays of cylindrical wells facilitate single‐molecule enzymology of α‐chymotrypsin

Single‐molecule enzymology allows scientists to examine the distributions of kinetic rates among members of a population. We describe a simple method for the analysis of single‐molecule enzymatic kinetics and provide comparisons to ensemble‐averaged kinetics. To isolate our model enzyme, α‐chymotrypsin, into single molecules, we use an array of cylindrical poly(dimethylsiloxane) wells 2 μm in diameter and 1.35 μm in height. Inside the wells, a protease assay with a profluorescent substrate detects α‐chymotrypsin activity. We hold the concentration of α‐chymotrypsin at 0.39 nM in a given well with an enzyme‐to‐substrate ratio of 1:6,666 molecules. Fluorescence emitted by the substrate is proportional to enzyme activity and detectable by a charge‐coupled device. This method allows for the simultaneous real‐time characterization of hundreds of individual enzymes. We analyze single‐molecule kinetics by recording and observing their intensity trajectories over time. By testing our method with our current instruments, we confirm that our methodology is useful for the analysis of single enzymes for extracting static inhomogeneity. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2009     

Structures of single‐layer β‐sheet proteins evolved from β‐hairpin repeats

Free‐standing single‐layer β‐sheets are extremely rare in naturally occurring proteins, even though β‐sheet motifs are ubiquitous. Here we report the crystal structures of three homologous, single‐layer, anti‐parallel β‐sheet proteins, comprised of three or four twisted β‐hairpin repeats. The structures reveal that, in addition to the hydrogen bond network characteristic of β‐sheets, additional hydrophobic interactions mediated by small clusters of residues adjacent to the turns likely play a significant role in the structural stability and compensate for the lack of a compact hydrophobic core. These structures enabled identification of a family of secreted proteins that are broadly distributed in bacteria from the human gut microbiome and are putatively involved in the metabolism of complex carbohydrates. A conserved surface patch, rich in solvent‐exposed tyrosine residues, was identified on the concave surface of the β‐sheet. These new modular single‐layer β‐sheet proteins may serve as a new model system for studying folding and design of β‐rich proteins.     

Clustering of multi‐domain protein sequences

The overall function of a multi‐domain protein is determined by the functional and structural interplay of its constituent domains. Traditional sequence alignment‐based methods commonly utilize domain‐level information and provide classification only at the level of domains. Such methods are not capable of taking into account the contributions of other domains in the proteins, and domain‐linker regions and classify multi‐domain proteins. An alignment‐free protein sequence comparison tool, CLAP (CLAssification of Proteins) was previously developed in our laboratory to especially handle multi‐domain protein sequences without a requirement of defining domain boundaries and sequential order of domains. Through this method we aim to achieve a biologically meaningful classification scheme for multi‐domain protein sequences. In this article, CLAP‐based classification has been explored on 5 datasets of multi‐domain proteins and we present detailed analysis for proteins containing (1) Tyrosine phosphatase and (2) SH3 domain. At the domain‐level CLAP‐based classification scheme resulted in a clustering similar to that obtained from an alignment‐based method. CLAP‐based clusters obtained for full‐length datasets were shown to comprise of proteins with similar functions and domain architectures. Our study demonstrates that multi‐domain proteins could be classified effectively by considering full‐length sequences without a requirement of identification of domains in the sequence.     

Expression of recombinant anti‐colorectal cancer large single‐chain monoclonal antibody in insect cells

Colorectal cancer is the third most commonly diagnosed cancer in the world. Monoclonal antibody (mAb) CO17‐1A recognizes the tumor‐associated antigen GA733‐2, a cell surface glycoprotein highly expressed in colorectal carcinoma cells which is applicable for preventing and curing colorectal cancer. In this study, we tried to produce a new recombinant anti‐colorectal cancer large single chain (lsc) mAb based on mAb CO17‐1A in the baculovirus‐insect cell protein expression system. Two kinds of recombinant lsc mAbs were generated where variable light chain (VL) and heavy chain (HC) of mAb CO17‐1A were fused together by an interchain linker. The only difference between the two mAbs is based on fusion of an ER retention signal (KDEL) at its C‐terminus of HC. Polymerase chain reaction analysis verified the presence of both recombinant genes in the bacmid for generating viral expression vectors in insect cells. Western blot confirmed the expression of lsc mAbs in baculovirus‐infected insect cells. Cell enzyme linked immunosorbent assay (ELISA) showed that the mAbs from cell lysates bound to SW480 and SW620 human colorectal cancer cells. These results indicate that the baculovirus insect expression system can produce anti‐colorectal lsc mAb recognizing human colorectal cancer cells.