Sequence tolerance of immunoglobulin variable domain framework regions to noncanonical intradomain disulfide linkages.

Most immunoglobulin (Ig) domains bear only a single highly conserved canonical intradomain, inter-ß-sheet disulfide linkage formed between Cys23-Cys104, and incorporation of rare noncanonical disulfide linkages at other locations can enhance Ig domain stability. Here, we exhaustively surveyed the sequence tolerance of Ig variable (V) domain framework regions (FRs) to noncanonical disulfide linkages. Starting from a destabilized VH domain lacking a Cys23-Cys104 disulfide linkage, we generated and screened phage-displayed libraries of engineered VHs, bearing all possible pairwise combinations of Cys residues in neighboring ß-strands of the Ig fold FRs. This approach identified seven novel Cys pairs in VH FRs (Cys4-Cys25, Cys4-Cys118, Cys5-Cys120, Cys6-Cys119, Cys22-Cys88, Cys24-Cys86, and Cys45-Cys100; the international ImMunoGeneTics information system numbering), whose presence rescued domain folding and stability. Introduction of a subset of these noncanonical disulfide linkages (three intra-ß-sheet: Cys4-Cys25, Cys22-Cys88, and Cys24-Cys86, and one inter-ß-sheet: Cys6-Cys119) into a diverse panel of VH, VL, and VHH domains enhanced their thermostability and protease resistance without significantly impacting expression, solubility, or binding to cognate antigens. None of the noncanonical disulfide linkages identified were present in the natural human VH repertoire. These data reveal an unexpected permissiveness of Ig V domains to noncanonical disulfide linkages at diverse locations in FRs, absent in the human repertoire, whose presence is compatible with antigen recognition and improves domain stability. Our work represents the most complete assessment to date of the role of engineered noncanonical disulfide bonding within FRs in Ig V domain structure and function.

A CHO stable pool production platform for rapid clinical development of trimeric SARS-CoV-2 spike subunit vaccine antigens.

Protein expression from stably transfected Chinese hamster ovary (CHO) clones is an established but time-consuming method for manufacturing therapeutic recombinant proteins. The use of faster, alternative approaches, such as non-clonal stable pools, has been restricted due to lower productivity and longstanding regulatory guidelines. Recently, the performance of stable pools has improved dramatically, making them a viable option for quickly producing drug substance for GLP-toxicology and early-phase clinical trials in scenarios such as pandemics that demand rapid production timelines. Compared to stable CHO clones which can take several months to generate and characterize, stable pool development can be completed in only a few weeks. Here, we compared the productivity and product quality of trimeric SARS-CoV-2 spike protein ectodomains produced from stable CHO pools or clones. Using a set of biophysical and biochemical assays we show that product quality is very similar and that CHO pools demonstrate sufficient productivity to generate vaccine candidates for early clinical trials. Based on these data, we propose that regulatory guidelines should be updated to permit production of early clinical trial material from CHO pools to enable more rapid and cost-effective clinical evaluation of potentially life-saving vaccines.

Characterization of a SARS-CoV-2 spike protein reference material.

Development of diagnostic testing capability has advanced with unprecedented pace in response to the COVID-19 pandemic. An undesirable effect of such speed is a lack of standardization, often leading to unreliable test results. To assist the research community surmount this challenge, the National Research Council Canada has prepared a SARS-CoV-2 spike protein reference material, SMT1-1, as a buffered solution. Value assignment was achieved by amino acid analysis (AAA) by double isotope dilution liquid chromatography-tandem mass spectrometry (LC-ID-MS/MS) following acid hydrolysis of the protein, in combination with ultraviolet-visible spectrophotometry (UV-Vis) based on tryptophan and tyrosine absorbance at 280 nm. Homogeneity of the material was established through spectrophotometric absorbance readings at 280 nm. Transportation and long-term storage stabilities were assessed by monitoring relative changes in oligomeric state by size-exclusion liquid chromatography (LC-SEC) with UV detection. The molar concentration of the spike protein in SMT1-1 was 5.68 ± 0.22 µmol L-1 (k = 2, 95% CI), with the native trimeric form accounting for ~ 94% of the relative abundance. Reference mass concentration and mass fraction values were calculated using the protein molecular weight and density of the SMT1-1 solution. The spike protein is highly glycosylated which leads to analyte ambiguity when reporting the more commonly used mass concentration. After glycoprotein molar mass determination by LC-SEC with multi-angle light scattering detection, we thus reported mass concentration values for both the protein-only portion and intact glycoprotein as 0.813 ± 0.030 and 1.050 ± 0.068 mg mL-1 (k = 2), respectively.

Structure-guided design of a potent Clostridiodes difficile toxin A inhibitor.

Crystal structures of camelid heavy-chain antibody variable domains (VHHs) bound to fragments of the combined repetitive oligopeptides domain of Clostridiodes difficile toxin A (TcdA) reveal that the C-terminus of VHH A20 was located 30 Å away from the N-terminus of VHH A26. Based on this observation, we generated a biparatopic fusion protein with A20 at the N-terminus, followed by a (GS)6 linker and A26 at the C-terminus. This A20-A26 fusion protein shows an improvement in binding affinity and a dramatic increase in TcdA neutralization potency (>330-fold [IC 50]; ≥2,700-fold [IC 99]) when compared to the unfused A20 and A26 VHHs. A20-A26 also shows much higher binding affinity and neutralization potency when compared to a series of control antibody constructs that include fusions of two A20 VHHs, fusions of two A26 VHHs, a biparatopic fusion with A26 at the N-terminus and A20 at the C-terminus (A26-A20), and actoxumab. In particular, A20-A26 displays a 310-fold (IC 50) to 29,000-fold (IC 99) higher neutralization potency than A26-A20. Size-exclusion chromatography-multiangle light scattering (SEC-MALS) analyses further reveal that A20-A26 binds to TcdA with 1:1 stoichiometry and simultaneous engagement of both A20 and A26 epitopes as expected based on the biparatopic design inspired by the crystal structures of TcdA bound to A20 and A26. In contrast, the control constructs show varied and heterogeneous binding modes. These results highlight the importance of molecular geometric constraints in generating highly potent antibody-based reagents capable of exploiting the simultaneous binding of more than one paratope to an antigen.

Corrigendum: Structure-guided design of a potent Clostridioides difficile toxin A inhibitor.

[This corrects the article DOI: 10.3389/fmicb.2023.1110541.].

Influence of variant-specific mutations, temperature and pH on conformations of a large set of SARS-CoV-2 spike trimer vaccine antigen candidates.

SARS-CoV-2 subunit vaccines continue to be the focus of intense clinical development worldwide. Protein antigens in these vaccines most commonly consist of the spike ectodomain fused to a heterologous trimerization sequence, designed to mimic the compact, prefusion conformation of the spike on the virus surface. Since 2020, we have produced dozens of such constructs in CHO cells, consisting of spike variants with different mutations fused to different trimerization sequences. This set of constructs displayed notable conformational heterogeneity, with two distinct trimer species consistently detected by analytical size exclusion chromatography. A recent report showed that spike ectodomain fusion constructs can adopt an alternative trimer conformation consisting of loosely associated ectodomain protomers. Here, we applied multiple biophysical and immunological techniques to demonstrate that this alternative conformation is formed to a significant extent by several SARS-CoV-2 variant spike proteins. We have also examined the influence of temperature and pH, which can induce inter-conversion of the two forms. The substantial structural differences between these trimer types may impact their performance as vaccine antigens.

An anti-HER2 biparatopic antibody that induces unique HER2 clustering and complement-dependent cytotoxicity.

Human epidermal growth factor receptor 2 (HER2) is a receptor tyrosine kinase that plays an oncogenic role in breast, gastric and other solid tumors. However, anti-HER2 therapies are only currently approved for the treatment of breast and gastric/gastric esophageal junction cancers and treatment resistance remains a problem. Here, we engineer an anti-HER2 IgG1 bispecific, biparatopic antibody (Ab), zanidatamab, with unique and enhanced functionalities compared to both trastuzumab and the combination of trastuzumab plus pertuzumab (tras + pert). Zanidatamab binds adjacent HER2 molecules in trans and initiates distinct HER2 reorganization, as shown by polarized cell surface HER2 caps and large HER2 clusters, not observed with trastuzumab or tras + pert. Moreover, zanidatamab, but not trastuzumab nor tras + pert, elicit potent complement-dependent cytotoxicity (CDC) against high HER2-expressing tumor cells in vitro. Zanidatamab also mediates HER2 internalization and downregulation, inhibition of both cell signaling and tumor growth, antibody-dependent cellular cytotoxicity (ADCC) and phagocytosis (ADCP), and also shows superior in vivo antitumor activity compared to tras + pert in a HER2-expressing xenograft model. Collectively, we show that zanidatamab has multiple and distinct mechanisms of action derived from the structural effects of biparatopic HER2 engagement.

Different active-site loop orientation in serine hydrolases versus acyltransferases.

Acyl transfer is a key reaction in biosynthesis, including synthesis of antibiotics and polyesters. Although researchers have long recognized the similar protein fold and catalytic machinery in acyltransferases and hydrolases, the molecular basis for the different reactivity has been a long-standing mystery. By comparison of X-ray structures, we identified a different oxyanion-loop orientation in the active site. In esterases/lipases a carbonyl oxygen points toward the active site, whereas in acyltransferases a NH of the main-chain amide points toward the active site. Amino acid sequence comparisons alone cannot identify such a difference in the main-chain orientation. To identify how this difference might change the reaction mechanism, we solved the X-ray crystal structure of Pseudomonas fluorescens esterase containing a sulfonate transition-state analogue bound to the active-site serine. This structure mimics the transition state for the attack of water on the acyl-enzyme and shows a bridging water molecule between the carbonyl oxygen mentioned above and the sulfonyl oxygen that mimics the attacking water. A possible mechanistic role for this bridging water molecule is to position and activate the attacking water molecule in hydrolases, but to deactivate the attacking water molecule in acyl transferases.

The crystal structure and dimerization interface of GADD45gamma.

Gadd45 proteins are recognized as tumor and autoimmune suppressors whose expression can be induced by genotoxic stresses. These proteins are involved in cell cycle control, growth arrest, and apoptosis through interactions with a wide variety of binding partners. We report here the crystal structure of Gadd45gamma, which reveals a fold comprising an alphabetaalpha sandwich with a central five-stranded mixed beta-sheet with alpha-helices packed on either side. Based on crystallographic symmetry we identified the dimer interface of Gadd45gamma dimers by generating point mutants that compromised dimerization while leaving the tertiary structure of the monomer intact. The dimer interface comprises a four-helix bundle involving residues that are the most highly conserved among Gadd45 isoforms. Cell-based assays using these point mutants demonstrate that dimerization is essential for growth inhibition. This structural information provides a new context for evaluation of the plethora of protein-protein interactions that govern the many functions of the Gadd45 family of proteins.

Structure of the noncatalytic domains and global fold of the protein disulfide isomerase ERp72.

Protein disulfide isomerases are a family of proteins that catalyze the oxidation and isomerization of disulfide bonds in newly synthesized proteins in the endoplasmic reticulum. The family includes general enzymes such as PDI that recognize unfolded proteins, and others that are selective for specific classes of proteins. Here, we report the X-ray crystal structure of central non-catalytic domains of a specific isomerase, ERp72 (also called CaBP2 and protein disulfide-isomerase A4) from Rattus norvegicus. The structure reveals strong similarity to ERp57, a PDI-family member that interacts with the lectin-like chaperones calnexin and calreticulin but, unexpectedly, ERp72 does not interact with calnexin as shown by isothermal titration calorimetry and nuclear magnetic resonance (NMR) spectroscopy. Small-angle X-ray scattering (SAXS) of ERp72 was used to develop models of the full-length protein using both rigid body refinement and ab initio simulated annealing of dummy atoms. The two methods show excellent agreement and define the relative positions of the five thioredoxin-like domains of ERp72 and potential substrate or chaperone binding sites.

Switching catalysis from hydrolysis to perhydrolysis in Pseudomonas fluorescens esterase.

Many serine hydrolases catalyze perhydrolysis, the reversible formation of peracids from carboxylic acids and hydrogen peroxide. Recently, we showed that a single amino acid substitution in the alcohol binding pocket, L29P, in Pseudomonas fluorescens (SIK WI) aryl esterase (PFE) increased the specificity constant of PFE for peracetic acid formation >100-fold [Bernhardt et al. (2005) Angew. Chem., Int. Ed. 44, 2742]. In this paper, we extend this work to address the three following questions. First, what is the molecular basis of the increase in perhydrolysis activity? We previously proposed that the L29P substitution creates a hydrogen bond between the enzyme and hydrogen peroxide in the transition state. Here we report two X-ray structures of L29P PFE that support this proposal. Both structures show a main chain carbonyl oxygen closer to the active site serine as expected. One structure further shows acetate in the active site in an orientation consistent with reaction by an acyl-enzyme mechanism. We also detected an acyl-enzyme intermediate in the hydrolysis of epsilon-caprolactone by mass spectrometry. Second, can we further increase perhydrolysis activity? We discovered that the reverse reaction, hydrolysis of peracetic acid to acetic acid and hydrogen peroxide, occurs at nearly the diffusion limited rate. Since the reverse reaction cannot increase further, neither can the forward reaction. Consistent with this prediction, two variants with additional amino acid substitutions showed 2-fold higher k(cat), but K(m) also increased so the specificity constant, k(cat)/K(m), remained similar. Third, how does the L29P substitution change the esterase activity? Ester hydrolysis decreased for most esters (75-fold for ethyl acetate) but not for methyl esters. In contrast, L29P PFE catalyzed hydrolysis of epsilon-caprolactone five times more efficiently than wild-type PFE. Molecular modeling suggests that moving the carbonyl group closer to the active site blocks access for larger alcohol moieties but binds epsilon-caprolactone more tightly. These results are consistent with the natural function of perhydrolases being either hydrolysis of peroxycarboxylic acids or hydrolysis of lactones.

Structure-based engineering of pH-dependent antibody binding for selective targeting of solid-tumor microenvironment.

Recent development of monoclonal antibodies as mainstream anticancer agents demands further optimization of their safety for use in humans. Potent targeting and/or effector activities on normal tissues is an obvious toxicity concern. Optimization of specific tumor targeting could be achieved by taking advantage of the extracellular acidity of solid tumors relative to normal tissues. Here, we applied a structure-based computational approach to engineer anti-human epidermal growth factor receptor 2 (Her2) antibodies with selective binding in the acidic tumor microenvironment. We used an affinity maturation platform in which dual-pH histidine-scanning mutagenesis was implemented for pH selectivity optimization. Testing of a small set of designs for binding to the recombinant Her2 ectodomain led to the identification of antigen-binding fragment (Fab) variants with the desired pH-dependent binding behavior. Binding selectivity toward acidic pH was improved by as much as 25-fold relative to the parental bH1-Fab. In vitro experiments on cells expressing intact Her2 confirmed that designed variants formatted as IgG1/k full-size antibodies have high affinity and inhibit the growth of tumor spheroids at a level comparable to that of the benchmark anti-Her2 antibody trastuzumab (Herceptin®) at acidic pH, whereas these effects were significantly reduced at physiological pH. In contrast, both Herceptin and the parental bH1 antibody exhibited strong cell binding and growth inhibition irrespective of pH. This work demonstrates the feasibility of computational optimization of antibodies for selective targeting of the acidic environment such as that found in many solid tumors.

Lectin control of protein folding and sorting in the secretory pathway.

Glycan moieties are essential for folding, sorting and targeting of glycoproteins through the secretory pathway to various cellular compartments. The molecular mechanisms that underlie these processes, however, are only now coming to light. Recent crystallographic and NMR studies of proteins located in the endoplasmic reticulum (ER), Golgi complex and ER-Golgi intermediate compartment have illuminated their roles in glycoprotein folding and secretion. Calnexin and calreticulin, both ER-resident proteins, have lectin domains that are crucial for their function as chaperones. The crystal structure of the carbohydrate-recognition domain of ER-Golgi intermediate compartment (ERGIC)-53 complements the biochemical and functional characterization of the protein, confirming that a lectin domain is essential for the role of this protein in sorting and transfer of glycoproteins from the ER to the Golgi complex. The lectin domains of calnexin and ERGIC-53 are structurally similar, although there is little primary sequence similarity. By contrast, sequence similarity between ERGIC-53 and vesicular integral membrane protein (VIP36), a Golgi-resident protein, leaves little doubt that a similar lectin domain is central to the transport and/or sorting functions of VIP36. The theme emerging from these studies is that carbohydrate recognition and modification are central to mediation of glycoprotein folding and secretion.

Binding symmetry and surface flexibility mediate antibody self-association.

Solution stability is an important factor in the optimization of engineered biotherapeutic candidates such as monoclonal antibodies because of its possible effects on manufacturability, pharmacology, efficacy and safety. A detailed atomic understanding of the mechanisms governing self-association of natively folded protein monomers is required to devise predictive tools to guide screening and re-engineering along the drug development pipeline. We investigated pairs of affinity-matured full-size antibodies and observed drastically different propensities to aggregate from variants differing by a single amino-acid. Biophysical testing showed that antigen-binding fragments (Fabs) from the aggregating antibodies also reversibly associated with equilibrium dissociation constants in the low-micromolar range. Crystal structures (PDB accession codes 6MXR, 6MXS, 6MY4, 6MY5) and bottom-up hydrogen-exchange mass spectrometry revealed that Fab self-association occurs in a symmetric mode that involves the antigen complementarity-determining regions. Subtle local conformational changes incurred upon point mutation of monomeric variants foster formation of complementary polar interactions and hydrophobic contacts to generate a dimeric Fab interface. Testing of popular in silico tools generally indicated low reliabilities for predicting the aggregation propensities observed. A structure-aggregation data set is provided here in order to stimulate further improvements of in silico tools for prediction of native aggregation. Incorporation of intermolecular docking, conformational flexibility, and short-range packing interactions may all be necessary features of the ideal algorithm.

Crystal structure of the bb' domains of the protein disulfide isomerase ERp57.

The synthesis of proteins in the endoplasmic reticulum (ER) is limited by the rate of correct disulfide bond formation. This process is carried out by protein disulfide isomerases, a family of ER proteins which includes general enzymes such as PDI that recognize unfolded proteins and others that are selective for specific proteins or classes. Using small-angle X-ray scattering and X-ray crystallography, we report the structure of a selective isomerase, ERp57, and its interactions with the lectin chaperone calnexin. Using isothermal titration calorimetry and NMR spectroscopy, we show that the b' domain of ERp57 binds calnexin with micromolar affinity through a conserved patch of basic residues. Disruption of this binding site by mutagenesis abrogates folding of RNase B in an in vitro assay. The relative positions of the ERp57 catalytic sites and calnexin binding site suggest that activation by calnexin is due to substrate recruitment rather than a direct stimulation of ERp57 oxidoreductase activity.

The structure of the exopolyphosphatase (PPX) from Escherichia coli O157:H7 suggests a binding mode for long polyphosphate chains.

Polyphosphate (polyP) is a linear polymer consisting of tens to hundreds of phosphate molecules joined together by high-energy anhydride bonds. These polymers are found in virtually all prokaryotic and eukaryotic cells and perform many functions; prominent among them are the responses to many stresses. Polyphosphate is synthesized by polyP kinase (PPK), using the terminal phosphate of ATP as the substrate, and degraded to inorganic phosphate by both endo- and exopolyphosphatases. Here we report the crystal structure and analysis of the polyphosphate phosphatase PPX from Escherichia coli O157:H7 refined at 2.2 Angstroms resolution. PPX is made of four domains. Domains I and II display structural similarity with one another and share the ribonuclease-H-like fold. Domain III bears structural similarity to the N-terminal, HD domain of SpoT. Domain IV, the smallest domain, has structural counterparts in cold-shock associated RNA-binding proteins but is of unknown function in PPX. The putative PPX active site is located at the interface between domains I and II. In the crystal structure of PPX these two domains are close together and represent the "closed" state. Comparison with the crystal structure of PPX/GPPA from Aquifex aeolicus reveals close structural similarity between domains I and II of the two enzymes, with the PPX/GPPA representing an "open" state. A striking feature of the dimer is a deep S-shaped canyon extending along the dimer interface and lined with positively charged residues. The active site region opens to this canyon. We postulate that this is a likely site of polyP binding.

Stability-Diversity Tradeoffs Impose Fundamental Constraints on Selection of Synthetic Human VH/VL Single-Domain Antibodies from In Vitro Display Libraries.

Human autonomous VH/VL single-domain antibodies (sdAbs) are attractive therapeutic molecules, but often suffer from suboptimal stability, solubility and affinity for cognate antigens. Most commonly, human sdAbs have been isolated from in vitro display libraries constructed via synthetic randomization of rearranged VH/VL domains. Here, we describe the design and characterization of three novel human VH/VL sdAb libraries through a process of: (i) exhaustive biophysical characterization of 20 potential VH/VL sdAb library scaffolds, including assessment of expression yield, aggregation resistance, thermostability and tolerance to complementarity-determining region (CDR) substitutions; (ii) in vitro randomization of the CDRs of three VH/VL sdAb scaffolds, with tailored amino acid representation designed to promote solubility and expressibility; and (iii) systematic benchmarking of the three VH/VL libraries by panning against five model antigens. We isolated ≥1 antigen-specific human sdAb against four of five targets (13 VHs and 7 VLs in total); these were predominantly monomeric, had antigen-binding affinities ranging from 5 nM to 12 µM (average: 2-3 µM), but had highly variable expression yields (range: 0.1-19 mg/L). Despite our efforts to identify the most stable VH/VL scaffolds, selection of antigen-specific binders from these libraries was unpredictable (overall success rate for all library-target screens: ~53%) with a high attrition rate of sdAbs exhibiting false positive binding by ELISA. By analyzing VH/VL sdAb library sequence composition following selection for monomeric antibody expression (binding to protein A/L followed by amplification in bacterial cells), we found that some VH/VL sdAbs had marked growth advantages over others, and that the amino acid composition of the CDRs of this set of sdAbs was dramatically restricted (bias toward Asp and His and away from aromatic and hydrophobic residues). Thus, CDR sequence clearly dramatically impacts the stability of human autonomous VH/VL immunoglobulin domain folds, and sequence-stability tradeoffs must be taken into account during the design of such libraries.

Mirror-image packing in enantiomer discrimination molecular basis for the enantioselectivity of B.cepacia lipase toward 2-methyl-3-phenyl-1-propanol.

Synthetic chemists often exploit the high enantioselectivity of lipases to prepare pure enantiomers of primary alcohols, but the molecular basis for this enantioselectivity is unknown. The crystal structures of two phosphonate transition-state analogs bound to Burkholderia cepacia lipase reveal this molecular basis for a typical primary alcohol: 2-methyl-3-phenyl-1-propanol. The enantiomeric alcohol moieties adopt surprisingly similar orientations, with only subtle differences that make it difficult to predict how to alter enantioselectivity. These structures, along with a survey of previous structures of enzyme bound enantiomers, reveal that binding of enantiomers does not involve an exchange of two substituent positions as most researchers assumed. Instead, the enantiomers adopt mirror-image packing, where three of the four substituents at the stereocenter lie in similar positions. The fourth substituent, hydrogen, points in opposite directions.

Mutational analysis of the gephyrin-related molybdenum cofactor biosynthetic gene cnxE from the lower eukaryote Aspergillus nidulans.

We report the identification of a number of mutations that result in amino acid replacements (and their phenotypic characterization) in either the MogA-like domain or domains 2 and 3 of the MoeA-like region of the Aspergillus nidulans cnxE gene. These domains are functionally required since mutations that result in amino acid substitutions in any one domain lead to the loss or to a substantial reduction in all three identified molybdoenzyme activities (i.e., nitrate reductase, xanthine dehydrogenase, and nicotinate hydroxylase). Certain cnxE mutants that show partial growth with nitrate as the nitrogen source in contrast do not grow on hypoxanthine or nicotinate. Complementation between mutants carrying lesions in the MogA-like domain or the MoeA-like region, respectively, most likely occurs at the protein level. A homology model of CnxE based on the dimeric structure of E. coli MoeA is presented and the position of inactivating mutations (due to amino acid replacements) in the MoeA-like functional region of the CnxE protein is mapped to this model. Finally, the activity of nicotinate hydroxylase, unlike that of nitrate reductase and xanthine dehydrogenase, is not restored in cnxE mutants grown in the presence of excess molybdate.

Antibody light chain variable domains and their biophysically improved versions for human immunotherapy.

We set out to gain deeper insight into the potential of antibody light chain variable domains (VLs) as immunotherapeutics. To this end, we generated a naïve human VL phage display library and, by using a method previously shown to select for non-aggregating antibody heavy chain variable domains (VHs), we isolated a diversity of VL domains by panning the library against B cell super-antigen protein L. Eight domains representing different germline origins were shown to be non-aggregating at concentrations as high as 450 µM, indicating VL repertoires are a rich source of non-aggregating domains. In addition, the VLs demonstrated high expression yields in E. coli, protein L binding and high reversibility of thermal unfolding. A side-by-side comparison with a set of non-aggregating human VHs revealed that the VLs had similar overall profiles with respect to melting temperature (T(m)), reversibility of thermal unfolding and resistance to gastrointestinal proteases. Successful engineering of a non-canonical disulfide linkage in the core of VLs did not compromise the non-aggregation state or protein L binding properties. Furthermore, the introduced disulfide bond significantly increased their T(m)s, by 5.5-17.5 ° C, and pepsin resistance, although it somewhat reduced expression yields and subtly changed the structure of VLs. Human VLs and engineered versions may make suitable therapeutics due to their desirable biophysical features. The disulfide linkage-engineered VLs may be the preferred therapeutic format because of their higher stability, especially for oral therapy applications that necessitate high resistance to the stomach's acidic pH and pepsin.