Structure of a low-melting-temperature anti-cholera toxin: llama V(H)H domain

TEMPRO: nanobody melting temperature estimation model using protein embeddings

Single-domain antibodies (sdAbs) or nanobodies have received widespread attention due to their small size (~ 15 kDa) and diverse applications in bio-derived therapeutics. As many modern biotechnology breakthroughs are applied to antibody engineering and design, nanobody thermostability or melting temperature (Tm) is crucial for their successful utilization. In this study, we present TEMPRO which is a predictive modeling approach for estimating the Tm of nanobodies using computational methods. Our methodology integrates various nanobody biophysical features to include Evolutionary Scale Modeling (ESM) embeddings, NetSurfP3 structural predictions, pLDDT scores per sdAb region from AlphaFold2, and each sequence's physicochemical characteristics. This approach is validated with our combined dataset containing 567 unique sequences with corresponding experimental Tm values from a manually curated internal data and a recently published nanobody database, NbThermo. Our results indicate the efficacy of protein embeddings in reliably predicting the Tm of sdAbs with mean absolute error (MAE) of 4.03 °C and root mean squared error (RMSE) of 5.66 °C, thus offering a valuable tool for the optimization of nanobodies for various biomedical and therapeutic applications. Moreover, we have validated the models' performance using experimentally determined Tms from nanobodies not found in NbThermo. This predictive model not only enhances nanobody thermostability prediction, but also provides a useful perspective of using embeddings as a tool for facilitating a broader applicability of downstream protein analyses.

Thermal Stability Estimation of Single Domain Antibodies Using Molecular Dynamics Simulations

In this chapter, we describe a protocol to estimate the thermal stability of single domain antibodies (sdAbs) using molecular dynamics (MD) simulations. This method measures the Q-value, the fraction of the native contacts, along the trajectory of high-temperature MD simulations starting from the experimental X-ray structure. We show a good correlation between the Q-value and the experimental melting temperature (T_m) in seven sdAbs. Assessing the Q-value on a per-residue level enabled us to identify residues that contribute to the instability and thus demonstrate which residues could be mutated to improve the stability and have later been validated by experiments. Our protocol extends beyond the application on sdAbs, as it is also suitable for other proteins and to determine the interfacial stability between protein and ligand.

Enhanced Conformational Sampling of Nanobody CDR H3 Loop by Generalized Replica-Exchange with Solute Tempering

The variable domains of heavy-chain antibodies, known as nanobodies, are potential substitutes for IgG antibodies. They have similar affinities to antigens as antibodies, but are more heat resistant. Their small size allows us to exploit computational approaches for structural modeling or design. Here, we investigate the applicability of an enhanced sampling method, a generalized replica-exchange with solute tempering (gREST) for sampling CDR-H3 loop structures of nanobodies. In the conventional replica-exchange methods, temperatures of only a whole system or scaling parameters of a solute molecule are selected for temperature or parameter exchange. In gREST, we can flexibly select a part of a solute molecule and a part of the potential energy terms as a parameter exchange region. We selected the CDR-H3 loop and investigated which potential energy term should be selected for the efficient sampling of the loop structures. We found that the gREST with dihedral terms can explore a global conformational space, but the relaxation to the global equilibrium is slow. On the other hand, gREST with all the potential energy terms can sample the equilibrium distribution, but the structural exploration is slower than with dihedral terms. The lessons learned from this study can be applied to future studies of loop modeling.

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.

Chromatibody, a novel non-invasive molecular tool to explore and manipulate chromatin in living cells

Chromatin function is involved in many cellular processes, its visualization or modification being essential in many developmental or cellular studies. Here, we present the characterization of chromatibody, a chromatin-binding single-domain, and explore its use in living cells. This non-intercalating tool specifically binds the heterodimer of H2A–H2B histones and displays a versatile reactivity, specifically labeling chromatin from yeast to mammals. We show that this genetically encoded probe, when fused to fluorescent proteins, allows non-invasive real-time chromatin imaging. Chromatibody is a dynamic chromatin probe that can be modulated. Finally, chromatibody is an efficient tool to target an enzymatic activity to the nucleosome, such as the DNA damage-dependent H2A ubiquitylation, which can modify this epigenetic mark at the scale of the genome and result in DNA damage signaling and repair defects. Taken together, these results identify chromatibody as a universal non-invasive tool for either in vivo chromatin imaging or to manipulate the chromatin landscape.

Summary: Chromatibody is a chromatin-binding single-domain antibody, derived from llama nanobodies, that can be used as a novel non-invasive molecular tool to explore and manipulate chromatin in living cells.

Structural and mutational analysis of a monomeric and dimeric form of a single domain antibody with implications for protein misfolding

Camelid single domain antibodies (sdAb) are known for their thermal stability and reversible refolding. We have characterized an unusually stable sdAb recognizing Staphylococcal enterotoxin B with one of the highest reported melting temperatures (T(m) = 85°C). Unexpectedly, ∼10-20% of the protein formed a dimer in solution. Three other cases where <20% of the sdAb dimerized have been reported; however, this is the first report of both the monomeric and dimeric X-ray crystal structures. Concentration of the monomer did not lead to the formation of new dimer suggesting a stable conformationally distinct species in a fraction of the cytoplasmically expressed protein. Comparison of periplasmic and cytoplasmic expression showed that the dimer was associated with cytoplasmic expression. The disulfide bond was partially reduced in the WT protein purified from the cytoplasm and the protein irreversibly unfolded. Periplasmic expression produced monomeric protein with a fully formed disulfide bond and mostly reversible refolding. Crystallization of a disulfide-bond free variant, C22A/C99V, purified from the periplasm yielded a structure of a monomeric form, while crystallization of C22A/C99V from the cytoplasm produced an asymmetric dimer. In the dimer, a significant conformational asymmetry was found in the loop residues of the edge β-strands (S50-Y60) containing the highly variable complementarity determining region, CDR2. Two dimeric assemblies were predicted from the crystal packing. Mutation of a residue at one of the interfaces, Y98A, disrupted the dimer in solution. The pleomorphic homodimer may yield insight into the stability of misfolded states and the importance of the conserved disulfide bond in preventing their formation.