Exploring single-domain antibody thermostability by molecular dynamics simulation

Nanobody engineering: computational modelling and design for biomedical and therapeutic applications

Nanobodies, the smallest functional antibody fragment derived from camelid heavy-chain-only antibodies, have emerged as powerful tools for diverse biomedical applications. In this comprehensive review, we discuss the structural characteristics, functional properties, and computational approaches driving the design and optimisation of synthetic nanobodies. We explore their unique antigen-binding domains, highlighting the critical role of complementarity-determining regions in target recognition and specificity. This review further underscores the advantages of nanobodies over conventional antibodies from a biosynthesis perspective, including their small size, stability, and solubility, which make them ideal candidates for economical antigen capture in diagnostics, therapeutics, and biosensing. We discuss the recent advancements in computational methods for nanobody modelling, epitope prediction, and affinity maturation, shedding light on their intricate antigen-binding mechanisms and conformational dynamics. Finally, we examine a direct example of how computational design strategies were implemented for improving a nanobody-based immunosensor, known as a Quenchbody. Through combining experimental findings and computational insights, this review elucidates the transformative impact of nanobodies in biotechnology and biomedical research, offering a roadmap for future advancements and applications in healthcare and diagnostics.

Novel Polystyrene-Binding Nanobody for Enhancing Immunoassays: Insights into Affinity, Immobilization, and Application Potential

Nanobodies, which represent the next generation of antibodies due to their unique properties, face a significant limitation in their poor physical adsorption on solid supports. In this study, we successfully discovered polystyrene binding nanobodies from a synthetic nanobody library. Notably, bivalent nanobody B2 exhibited high affinity for polystyrene (0.7 nM for ELISA saturation binding analysis and 15.6 nM for isothermal titration calorimetry), displaying a pH-dependent behavior. Remarkably, hydrophobic and electrostatic interactions contribute minimally to the binding process. Molecular modeling provided insights into the interaction between B2 and polystyrene, revealing that the Trp51 residue within the CDR2 loop formed an aromatic H-bond with polystyrene at a distance of 2.74 Å, thus explaining the observed reduction in B2 affinity caused by Trp51 mutations. To explore B2's potential in protein immobilization, we constructed a bispecific nanobody by fusing B2 to an anticarcinoembryonic antigen nanobody 11C12, which cannot be immobilized on polystyrene through passive adsorption. Remarkably, the fusion construct achieved effective immobilization on polystyrene within 5 min by passing the need for periplasmic protein purification despite its low expression level. Moreover, the fusion construct demonstrated excellent linearity in the chemiluminescent enzyme immunoassay. For the first time, this study reports a simplified and seamless platform for the oriented immobilization of nanobody. Importantly, the entire process eliminated the need for protein purification, enabling efficient and rapid immobilization of fusion proteins directly from crude cell extracts, even when the expression level was low. Our developed process dramatically reduced the processing time from 2.5 days to just 5 min.

Exploring the mechanism of the PTP1B inhibitors by molecular dynamics and experimental study

Protein tyrosine phosphatase 1B (PTP1B) has proven to be an attractive target for the treatment of cancer, diabetes and other diseases. Although many PTP1B inhibitors with various scaffolds have been developed, there is still a lack of PTP1B inhibitor with high specificity and acceptable pharmacological properties. Therefore, it is urgent to develop more methods to explore complex action mode of PTP1B and ligands for designing ideal PTP1B modulators. In this work, we developed a potential molecular dynamics (MD) analytic mode to analyze the mechanism of active compounds 6a and 6e against PTP1B from different perspectives, including the stable ability, interactions and binding site of ligand and protein, the binding energy, relative movement between residues and changes in protein internal interactions. The simulated results demonstrated that compound 6a bound more stably to the active pocket of PTP1B than 6e due to its smaller molecular volume (326 Å3), matched electronegativity, and enhanced the positive correlation motion of residues, especially for WPD loop and P loop. Lastly, compound 6a as a competitive inhibitor for PTP1B was verified by enzyme kinetic assay. This work successfully studied the mechanism of compound 6a against PTP1B from various aspects, enriched the analysis of interaction mode between PTP1B and inhibitors. In summary, we hope that this work could provide more theoretical information for designing and developing more novel and ideal PTP1B inhibitors in the future.

Characterization of Foot-and-Mouth Disease Virus Serotype O-Specific Single Domain Antibody Expressed in the pET Expression System

Foot and mouth disease (FMD) is an extremely contagious disease of cloven-hoofed domesticated and wild animals, resulting in significant economic losses in many parts of the world. FMD virus (FMDV) serotype O is responsible for approximately 70% of global outbreaks. For detection of FMDV antigen or antibody, ELISAs are used worldwide and have several limitations, such as batch-to-batch variation in generating immunobiologicals, high production cost and ethical concerns over animal sacrifice. The use of single domain antibody (sdAb) or variable N-terminal domain of the heavy chain of heavy chain antibody (VHH) found naturally in camels has proven their effectiveness in diagnostics and therapeutics. In the present study, the anti-FMDV serotype O-specific VHH-C1 gene sequence (Accession no. KJ751546) was retrieved from the NCBI database. The gene was synthesized commercially in the pBluescript KS+ cloning vector and expressed in E. coli BL21 (DE3) cells using the pET303/CT-His expression system with a C-terminal 6X-His tag. The expressed sdAb, verified by SDS‒PAGE and western blotting, was purified by Ni-chelate chromatography and used as a coating antibody in double antibody sandwich (DAS) ELISA for FMDV detection and typing. The sdAb exhibited a high binding affinity for FMDV serotype O, without any cross-reactivity toward serotypes A and Asia-1. It exhibited better thermostability up to 85 °C than conventional rabbit polyclonal anti-FMDV sera. The potential of sdAbs thus produced without sacrificing lab animals could be explored for replacing polyclonal sera in DAS-ELISA as well as for developing biosensors or lateral flow devices for FMDV type O detection.

NANOBODIES®: A Review of Diagnostic and Therapeutic Applications

NANOBODY® (a registered trademark of Ablynx N.V) molecules (Nbs), also referred to as single domain-based VHHs, are antibody fragments derived from heavy-chain only IgG antibodies found in the Camelidae family. Due to their small size, simple structure, high antigen binding affinity, and remarkable stability in extreme conditions, nanobodies possess the potential to overcome several of the limitations of conventional monoclonal antibodies. For many years, nanobodies have been of great interest in a wide variety of research fields, particularly in the diagnosis and treatment of diseases. This culminated in the approval of the world's first nanobody based drug (Caplacizumab) in 2018 with others following soon thereafter. This review will provide an overview, with examples, of (i) the structure and advantages of nanobodies compared to conventional monoclonal antibodies, (ii) methods used to generate and produce antigen-specific nanobodies, (iii) applications for diagnostics, and (iv) ongoing clinical trials for nanobody therapeutics as well as promising candidates for clinical development.

General Trends of the Camelidae Antibody VHHs Domain Dynamics

Conformational flexibility plays an essential role in antibodies' functional and structural stability. They facilitate and determine the strength of antigen-antibody interactions. Camelidae express an interesting subtype of single-chain antibody, named Heavy Chain only Antibody. They have only one N-terminal Variable domain (V_HH) per chain, composed of Frameworks (FRs) and Complementarity Determining regions (CDRs) like their VH and VL counterparts in IgG. Even when expressed independently, V_HH domains display excellent solubility and (thermo)stability, which helps them to retain their impressive interaction capabilities. Sequence and structural features of V_HH domains contributing to these abilities have already been studied compared to classical antibodies. To have the broadest view and understand the changes in dynamics of these macromolecules, large-scale molecular dynamics simulations for a large number of non-redundant V_HH structures have been performed for the first time. This analysis reveals the most prevalent movements in these domains. It reveals the four main classes of V_HHs dynamics. Diverse local changes were observed in CDRs with various intensities. Similarly, different types of constraints were observed in CDRs, while FRs close to CDRs were sometimes primarily impacted. This study sheds light on the changes in flexibility in different regions of V_HH that may impact their in silico design.

VHH Structural Modelling Approaches: A Critical Review

V_HH, i.e., VH domains of camelid single-chain antibodies, are very promising therapeutic agents due to their significant physicochemical advantages compared to classical mammalian antibodies. The number of experimentally solved V_HH structures has significantly improved recently, which is of great help, because it offers the ability to directly work on 3D structures to humanise or improve them. Unfortunately, most V_HHs do not have 3D structures. Thus, it is essential to find alternative ways to get structural information. The methods of structure prediction from the primary amino acid sequence appear essential to bypass this limitation. This review presents the most extensive overview of structure prediction methods applied for the 3D modelling of a given V_HH sequence (a total of 21). Besides the historical overview, it aims at showing how model software programs have been shaping the structural predictions of V_HHs. A brief explanation of each methodology is supplied, and pertinent examples of their usage are provided. Finally, we present a structure prediction case study of a recently solved V_HH structure. According to some recent studies and the present analysis, AlphaFold 2 and NanoNet appear to be the best tools to predict a structural model of V_HH from its sequence.

Engineering of single-domain antibodies for next-generation snakebite antivenoms

Given the magnitude of the global snakebite crisis, strategies to ensure the quality of antivenom, as well as the availability and sustainability of its supply are under development by several research groups. Recombinant DNA technology has allowed the engineering of monoclonal antibodies and recombinant fragments as alternatives to conventional antivenoms. Besides having higher therapeutic efficacy, with broad neutralization capacity against local and systemic toxicity, novel antivenoms need to be safe and cost-effective. Due to the biological and physical chemical properties of camelid single-domain antibodies, with high volume of distribution to distal tissue, their modular format, and their versatility, their biotechnological application has grown considerably in recent decades. This article presents the most up-to-date developments concerning camelid single-domain-based antibodies against major toxins from snake venoms, the main venomous animals responsible for reported envenoming cases and related human deaths. A brief discussion on the composition, challenges, and perspectives of antivenoms is presented, as well as the road ahead for next-generation antivenoms based on single-domain antibodies.

Exploring the dynamic mechanism of allosteric drug SHP099 inhibiting SHP2E69K

The E69K mutation is one of the most frequent protein tyrosine phosphatase-2 (SHP2) mutations in leukemia, and it can cause the increase in the protein activity. Recent studies have shown that the E69K mutation was fairly sensitive to the allosteric inhibitor of SHP2 (SHP099). However, the molecular mechanism of the allosteric drug SHP099 inhibiting SHP2^E69K remains unclear. Thus, the molecular dynamic simulations and the post-dynamics analyses (RMSF, PCA, DCCM, RIN and the binding free energies) for SHP2^WT, SHP2^WT-SHP099, SHP2^E69K and SHP2^E69K-SHP099 were carried out, respectively. Owing to the strong binding affinity of SHP099 to residues Thr219 and Arg220, the flexibility of linker region (residues Val209-Arg231) was reduced. Moreover, the presence of SHP099 kept the autoinhibition state of the SHP2 protein through enhancing the interactions between the linker region and Q loop in PTP domain, such as Thr219/Val490, Thr219/Asn491, Arg220/Ile488 and Leu254/Asn491. In addition, it was found that the residues (Thr219, Arg220, Leu254 and Asn491) might be the key residues responsible for the conformational changes of protein. Overall, this study may provide an important basis for understanding how the SHP099 effectively inhibited the SHP2^E69K activity at the molecular level.

Computational Design of Nitrile Hydratase from Pseudonocardia thermophila JCM3095 for Improved Thermostability

High thermostability and catalytic activity are key properties for nitrile hydratase (NHase, EC 4.2.1.84) as a well-industrialized catalyst. In this study, rational design was applied to tailor the thermostability of NHase from Pseudonocardia thermophila JCM3095 (PtNHase) by combining FireProt server prediction and molecular dynamics (MD) simulation. Site-directed mutagenesis of non-catalytic residues provided by the rational design was subsequentially performed. The positive multiple-point mutant, namely, M10 (αI5P/αT18Y/αQ31L/αD92H/βA20P/βP38L/βF118W/βS130Y/βC189N/βC218V), was obtained and further analyzed. The Melting temperature (T_m) of the M10 mutant showed an increase by 3.2 °C and a substantial increase in residual activity of the enzyme at elevated temperatures was also observed. Moreover, the M10 mutant also showed a 2.1-fold increase in catalytic activity compared with the wild-type PtNHase. Molecular docking and MD simulations demonstrated better substrate affinity and improved thermostability for the mutant.

Identification of potential leukocyte antigen-related protein (PTP-LAR) inhibitors through 3D QSAR pharmacophore-based virtual screening and molecular dynamics simulation

Owing to its negative regulatory role in insulin signaling, protein tyrosine phosphatase of leukocyte antigen-related protein (PTP-LAR) was widely thought as a potential drug target for diabetes. Now, it was urgent to search for potential LAR inhibitors targeting diabetes. Initially, the pharmacophore models of LAR inhibitors were established with the application of the HypoGen module. The cost analysis, test set validation, as well as Fischer's test was used to verify the efficiency of pharmacophore model. Then, the best pharmacophore model (Hypo-1-LAR) was applied for the virtual screening of the ZINC database. And 30 compounds met the Lipinski's rule of five. Among them, 10 compounds with better binding affinity than the known LAR inhibitor (BDBM50296375) were discovered by docking studies. Finally, molecular dynamics simulations and post-analysis experiments (RMSD, RMSF, PCA, DCCM and RIN) were conducted to explore the effect of ligands (ZINC97018474 and Compound 1) on LAR and preliminary understand why ZINC97018474 had better inhibitory activity than Compound 1 (BDBM50296375). Communicated by Ramaswamy H. Sarma.