Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci. , 10 March 2025

Sec. Structural Biology

Volume 12 - 2025 | https://doi.org/10.3389/fmolb.2025.1512788

Computational electrostatic engineering of nanobodies for enhanced SARS−CoV−2 receptor binding domain recognition

  • 1Central Laboratories, King Faisal University, Al Hofuf, Saudi Arabia
  • 2Centre of Agricultural Biochemistry and Biotechnology (CABB), University of Agriculture, Faisalabad, Pakistan
  • 3Medical and Allied Department, Faisalabad Medical University, Faisalabad, Pakistan
  • 4Department of Biochemistry, Government College University Faisalabad, Faisalabad, Pakistan

This study presents a novel computational approach for engineering nanobodies (Nbs) for improved interaction with receptor-binding domain (RBD) of the SARS-CoV-2 spike protein. Using Protein Structure Reliability reports, RBD (7VYR_R) was selected and refined for subsequent Nb-RBD interactions. By leveraging electrostatic complementarity (EC) analysis, we engineered and characterized five Electrostatically Complementary Nbs (ECSb1-ECSb5) based on the CeVICA library’s SR6c3 Nb. Through targeted modifications in the complementarity-determining regions (CDR) and framework regions (FR), we optimized electrostatic interactions to improve binding affinity and specificity. The engineered Nbs (ECSb3, ECSb4, and ECSb5) demonstrated high binding specificity for AS3, CA1, and CA2 epitopes. Interestingly, ECSb1 and ECSb2 selectively engaged with AS3 and CA1 instead of AS1 and AS2, respectively, due to a preference for residues that conferred superior binding complementarities. Furthermore, ECSbs significantly outperformed SR6c3 Nb in MM/GBSA results, notably, ECSb4 and ECSb3 exhibited superior binding free energies of −182.58 kcal.mol-1 and −119.07 kcal.mol-1, respectively, compared to SR6c3 (−105.50 kcal.mol-1). ECSbs exhibited significantly higher thermostability (100.4–148.3 kcal·mol⁻1) compared to SR6c3 (62.6 kcal·mol⁻1). Similarly, enhanced electrostatic complementarity was also observed for ECSb4-RBD and ECSb3-RBD (0.305 and 0.390, respectively) relative to SR6c3-RBD (0.233). Surface analyses confirmed optimized electrostatic patches and reduced aggregation propensity in the engineered Nb. This integrated EC and structural engineering approach successfully developed engineered Nbs with enhanced binding specificity, increased thermostability, and reduced aggregation, laying the groundwork for novel therapeutic applications targeting the SARS-CoV-2 spike protein.

1 Introduction

The COVID−19 pandemic, caused by SARS-CoV-2, created a global health crisis of unprecedented magnitude, overwhelming healthcare systems (Ranney et al., 2020) and leading to substantial loss of life (Wu and McGoogan, 2020), financial strain (Nicola et al., 2020), unemployment (Coibion et al., 2021), and widespread supply chain disruptions (Golan et al., 2020). Additionally, the pandemic highlighted social and economic disparities, with marginalized communities often bearing a disproportionate burden of the virus’s impact (Egede and Walker, 2020; Tai et al., 2021). While posing immense challenges, the pandemic also spurred innovation in vaccine development (Krammer, 2020) and restructuring of work environments (Kniffin et al., 2021), reshaping the future of global health and society. Notwithstanding this progress, the virus’s swiftly evolving nature necessitated a prompt and pressing need for research on efficacious therapies to curb emerging/evolving strains and mitigate the pandemic’s ongoing impact (Graham, 2020). Developing such therapies entails an understanding of the molecular mechanisms underlying SARS−CoV−2 infection, particularly the role of the spike (S) protein in facilitating viral entry into host cells.

The S protein comprises S1 and S2 subunits, with S2 containing highly conserved regions (Premkumar et al., 2020). A successful SARS-CoV-2 infection hinges on the specific binding of the viral spike (S) protein to cellular entry receptors, primarily the angiotensin−converting enzyme 2 (ACE2) (Lan et al., 2020; Shang et al., 2020; Walls et al., 2020; Wang et al., 2020; Zhou et al., 2020; V’kovski et al., 2021). This receptor engagement has been a key factor in the virus’s adaptation from an animal reservoir and its ongoing evolution (Letko et al., 2020). The immune system combats SARS−CoV−2 by generating a diverse array of antibodies that target various viral epitopes, thereby inhibiting cellular entry, and by activating T cells that eliminate infected host cells. Monoclonal antibodies (mAbs) targeting conserved regions of the S1 or S2 subunits, such as hmAb 1249A8, S2P6, and CC40.8, have demonstrated broad efficacy against SARS-CoV-2 variants (Tian et al., 2020; Cameroni et al., 2021; Chen P. et al., 2021; Starr et al., 2021; Park et al., 2022). These mAbs can directly interfere with viral pathogenesis by neutralizing virions, opsonizing infected cells, or inducing cell death. While mAbs offer broad efficacy, their widespread use is hindered by high production costs (Flaxman et al., 2022; Popping et al., 2023) and the risk of repeated viral evasion (Liu et al., 2021; Wang Q. et al., 2022). Therefore, alternative strategies are needed to address these challenges and develop more affordable and sustainable countermeasures against SARS−CoV−2 infection.

Nanobodies (Nbs), as single−domain antibodies (Hamers-Casterman et al., 1993), offer several distinct advantages over mAbs, including reduced size, enhanced stability, lower immunogenicity, and streamlined manufacturing process (Muyldermans, 2013; Steeland et al., 2016). Their therapeutic potential has garnered increasing attention, particularly following the FDA’s approval of the first Nb therapeutic in 2019 (Morrison, 2019). The superior stability of a Nb arises from their compact, single−domain structure and their highly polar surface. Nevertheless, achieving optimal Nb performance requires careful optimization of their binding affinity and stability (Kunz et al., 2018; Esparza et al., 2020) which are influenced by a sophisticated interplay of factors, including hydrogen bonding, hydrophobic interactions, van der Waals forces, and electrostatic interactions (Chakrabarti and Janin, 2002; Dellisanti et al., 2007; Boehr et al., 2009; Kastritis and Bonvin, 2013). This intricate balance of forces orchestrates robust and highly selective protein complex formation, underscoring their profound significance given that approximately 80% of proteins function within such assemblies (Gavin et al., 2002; de Las Rivas and Fontanillo, 2010).

Collectively, these forces define binding properties by establishing multiple layers of complementarity—surface (Sc), hydrophobic, and electrostatic (EC)—between the interacting entities (Korn and Burnett, 1991; Braden and Poljak, 1995; Meyer et al., 1996; McCoy et al., 1997; Sinha and Smith-Gill, 2002; Basu and Biswas, 2018). Sc arises from the precise geometric fit of molecules facilitated by van der Waals forces and hydrogen bonding (Chothia and Janin, 1975), while hydrophobic complementarity is driven by the association of nonpolar side chains through hydrophobic interactions (Richards, 1977). Conversely, the alignment of opposite surface charges establishes EC, which plays a pivotal role in protein−protein interactions (PPIs) (McCoy et al., 1997; Sheinerman et al., 2000). EC facilitates the initial encounter between binding partners through long−range electrostatic attractions, effectively steering them toward the correct orientation for binding (Vijayakumar et al., 1998; Sheinerman et al., 2000). This directed guidance reduces the entropic barrier associated with random collisions, thereby increasing the likelihood of productive interactions (Schreiber and Fersht, 1995; Selzer et al., 2000). Furthermore, electrostatic interactions enhance binding strength by forming salt bridges and hydrogen bonds at the interface, augmenting both affinity and specificity (Baker and Hubbard, 1984; Jones and Thornton, 1996; Xu et al., 1997). By stabilizing the transition state and lowering the activation energy required for complex formation, EC significantly contributes to the robustness and efficacy of PPIs (Vijayakumar et al., 1998; Sheinerman et al., 2000). Recognizing its critical role, it becomes essential to explore innovative strategies that leverage EC to enhance the binding properties of therapeutic molecules such as Nbs.

However, despite the significant EC role, previous computational approaches for affinity maturation of Nbs targeting the SARS−CoV−2 RBD have primarily relied on two strategies: (1) grafting complementarity−determining regions (CDRs) from neutralizing Abs onto stable scaffold frameworks (Cao et al., 2020; Wu et al., 2020a; Yang et al., 2021; 2023; Gopal et al., 2022; Rangel et al., 2022; Ferraz et al., 2024), and (2) structure−based site−mutagenesis of CDR amino acids within Nb-RBD docked complexes (Laroche et al., 2022; Longsompurana et al., 2023; Zhu et al., 2023; Hannula et al., 2024; Singh et al., 2024). To the best of our knowledge, the simultaneous engineering of CDRs and FRs has not been previously explored, as prior research has typically focused on modifying either CDRs or FRs independently to enhance binding affinity through electrostatic attractions (Sinha et al., 2002; Fukunaga and Tsumoto, 2013; Kiyoshi et al., 2014; Du Q. et al., 2019; Yoshida et al., 2019; Nguyen et al., 2021). This novel approach allows for a more comprehensive optimization of Nb interactions with the SARS−CoV−2 RBD, potentially can lead to the development of more effective therapeutic agents. Leveraging computational techniques, we have predicted the Nbs with enhanced paratope−epitope EC employing this methodology. Unlike earlier EC enhancing studies that primarily targeted binding affinity, our approach systematically evaluates the impact of this strategy on key parameters—including binding interfaces, binding free energy, thermostability, aggregation propensity, and post−translational modification reactivity—with the objective of predicting Nbs with superior binding affinity and stability. To substantiate our approach, we conducted a comparative analysis with the well−characterized Nb SR6c3 from the CeVICA library, which is renowned for its robust neutralizing ability and high thermal stability. Our results on Nb engineering with predicted improvement provide a computational framework for developing effective Nbs against SARS−CoV−2. This approach can offer a promising avenue for engineering of novel and effective therapeutic options by optimizing Sc and hydrophobic complementarity.

2 Materials and methods

2.1 Sequence retrieval and RBD structure analysis

A dataset of ten RBD structures, including 7EAM_A (Li et al., 2021), 7EAM_B (Li et al., 2021), 7B0B_F (Korenkov et al., 2023), 7BZ5_A (Wu et al., 2020b), 6M0J_E (Lan et al., 2020), 7VYR_C (Jeong et al., 2022), 7VYR_R (Jeong et al., 2022), 7BNV_A (Bullen et al., 2021), 8EL2_A (Hermet et al., 2023), and 8EL2_B (Hermet et al., 2023), was retrieved from the Protein Data Bank (PDB). These structures pre−processed employing the Protein Preparation Wizard (PPW) (Madhavi Sastry et al., 2013; Schrödinger, 2023e) with default settings in BioLuminate (Beard et al., 2013; Salam et al., 2014; Zhu et al., 2014; Schrödinger, 2023a) to ensure optimal quality. Subsequently, Protein Structure Reliability (PSR) analysis was performed to assess most suitable structure for the subsequent energy minimization and loop refinement using Prime module (Jacobson et al., 2002; Jacobson et al., 004; Schrödinger, 2023d) employing the VSGB solvation model (Li et al., 2011) and the OPLS3e force field (Roos et al., 2019). This was achieved by specifically addressing the amino acids that displayed high−temperature factors, peptide planarity problems, or steric clashes while maintaining the original sequence as provided. Finally, the structures were further validated through PSR, Structure Quality Reports (SQR), and Ramachandran plot analysis (Ramachandron et al., 1963).

2.2 Epitope selection and their biophysical analysis

We focused on three previously reported neutralizing epitopes, AS1 (Hansen et al., 2020; Pinto et al., 2020), AS2 (Hansen et al., 2020; Robbiani et al., 2020), and AS3 (Yuan et al., 2020), which are established binding sites for single−chain variable fragment (scFv) antibodies (Figure 1A). Additionally, we expanded our repertoire by incorporating two more recently reported immunoglobulin G (IgG)−targeted epitopes, CA1 and CA2 (Contreras et al., 2023). This comprehensive approach leverages both well−studied and newly discovered epitopes to maximize the potential for effective neutralization. To comprehensively characterize the selected epitopes, a detailed analysis of their titration curves and ionizable residues was conducted in BioLuminate. The protonation states of these ionizable residues were determined by employing PROPKA (Li et al., 2005; Bas et al., 2008; Olsson et al., 2011) with default parameters, ensuring accurate representation of the electrostatic properties within the molecular environment.

Figure 1
www.frontiersin.org

Figure 1. (A) Structure of SARS−CoV−2 spike protein (ribbon representation) and neutralizing epitopes on receptor binding domain (RBD), visualized by ChimeraX 1.8. Spike protein structure labelled with sub−domains S1 (Blue) and S2 (Golden) with particular indication of RBD (Red) (A). Enlarged view of the RBD showing the locations of three neutralizing epitopes: AS1, AS2, and AS3 (B). Additional view of RBD highlighting two more recently identified epitopes: CA1 and CA2 (C). (B) Integrated approach of CDR designing and FR engineering to develop a greater electrostatic complementarity. Initial nanobody (blue) and antigen (red) interaction (left) with less EC against epitope with fewer counterpart charged amino acids in CDRs (orange) and FRs (blue). The circular inset shows the amino acid sequence at the binding interface, with the antigen sequence in purple and the nanobody sequence in blue and orange. Modified nanobody with optimized amino acid sequence (Centre) with complete randomization of amino acids in CDRs (green) with engineered FRs (lavender) generating a greater EC against epitope with greater electrostatic potential (blue lines on CDRs).

2.3 Complementarity−determining regions (CDR) designing

SR6c3, a CeVICA Nb with IC₅₀ of 62.7 nM and melting temperature of 72°C, was selected as the reference Nb for the current study due to its robust neutralizing capability and highest thermal stability within the Nb library (Chen X. et al., 2021). Although the monomeric form of SR6v15 exhibits a superior inhibitory potency with an IC₅₀ of 2.18 nM, SR6c3 was preferred for its exceptional thermal resilience coupled with substantial neutralizing efficacy. SR6v15 demonstrated lower refolding capacity, with a heated/non−heated ratio of 0.72, in contrast to SR6c3, which achieved a ratio greater than one indicating an enhanced binding affinity following complete thermal denaturation at 98°C and refolding. Therefore, SR6c3 was chosen based on its superior performance.

Using standard frame sequences of SR6c3 as scaffolds, we designed CDRs to enhance EC between the paratope and the epitope, by prioritizing residues with opposing charges and significant sidechain pKa. For neutral epitope residues, polar amino acids with high free energies of association were chosen (Du H. et al., 2019). CDR1 and CDR2 lengths were constrained to 5-7 amino acids to prevent the potential structural instability, while CDR3 length was allowed to vary randomly. The incorporation of cysteine residues, a defining characteristic of the CeVICA library, was also considered during CDR design (Figure 1B). ColabFold v1.5.2 (Mirdita et al., 2022) was used, with default settings, to predict the structures of five Nbs involving CDR design, alongside the reference Nb SR6c3.

2.4 Framework regions (FR) engineering

The initial Nb structures were subjected to pre-processing using the PPW program with the default settings (Madhavi Sastry et al., 2013; Schrödinger, 2023e) in BioLuminate (Beard et al., 2013; Salam et al., 2014; Zhu et al., 2014; Schrödinger, 2023a). Unlike traditional sequence-based CDR designing approaches, our strategy for engineering the FRs of SR6c3 was essentially structure−based. Pre−processed Nbs were carefully examined using BioLuminate to identify specific residues that are suitable for mutation, as proposed by Fukunaga et al., (Fukunaga and Tsumoto, 2013). A selection range of five residues was defined in the vicinity of the CDRs from FR1, FR2, and FR4, while for longer FR3, this range was extended to eight. It is noteworthy that FR2 hallmark residues (F37, E44, R45, and V48) were deliberately left un-mutated due to their crucial role in the independence of the light chain (Muyldermans et al., 1994; Padlan, 1994; Vu et al., 1997; Muyldermans, 2013).

We refined our residue selection process based on Fukunaga et al., (Fukunaga and Tsumoto, 2013), targeting residues adjacent to CDRs, exposed to the molecular surface, and with sidechains at least 20% accessible to the solvent to ensure the relevance of our mutation targets. Following these specified criteria, we analyzed the electrostatic topography of the RBD and epitopes, then formulated a mutation strategy to introduce charged amino acids into each Nb structure. The initial phase of the strategy involved mutating the neighboring residues around the CDRs to amplify EC between the CDRs region and the epitope, while mutating outlying residues to improve the overall EC between Nb and RBD. Furthermore, we substituted specific hydrophobic residues with polar ones to reduce aggregation and enhance solubility (Figure 1B).

For the computational evaluation of Nbs’ humanization profiles, AbNatiV (Ramon et al., 2024), a deep learning tool was used with default settings. The amino acid sequences of the newly designed Nbs and SR6c3 were submitted to AbNatiV, utilizing the VHH domain mode. However, alignment with CeVICA revealed the adoption of the same VHH hallmark residues (F37, E44, R45, and V48) in lieu of human residues, reflecting their advantageous role in post−minus pre−affinity maturation (Chen X. et al., 2021). ColabFold v1.5.2 (Mirdita et al., 2022) was then employed with default settings to predict the structural conformations of the five Electrostatically Complementary Nbs (ECSbs) and SR6c3 based on their consensus sequences.

2.5 Structural refinement and analyses

Initial structural refinement of SR6c3 and ECSbs was performed using PPW (Madhavi Sastry et al., 2013; Schrödinger, 2023e) with default settings in BioLuminate (Beard et al., 2013; Salam et al., 2014; Zhu et al., 2014; Schrödinger, 2023a), addressing potential issues like steric clashes, peptide planarity, temperature factors, and dihedral angles. Each ECSb then underwent a meticulous PSR assessment, resulting in the selection of SR6c3, ECSb2, ECSb3, and ECSb5 for subsequent surface analyses. ECSb1 and ECSb4 required additional energy minimizations and loop refinement at specific amino acids, which were conducted using Prime employing the VSGB solvation model (Li et al., 2011) and the OPLS3e force field (Roos et al., 2019), followed by PSR reassessment and structure finalization (Jacobson et al., 2002; 2004; Schrödinger, 2023d). All six structures were further evaluated using SQRs and Ramachandran plots (Ramachandran et al., 1963). Post evaluation, an exhaustive analysis of electrostatically charged patches, aggregation scores using AggScore (Sankar et al., 2018), and residues reactivity profiles of each structure was conducted. These comprehensive analyses, performed using BioLuminate’s protein surface analysis tool, were subsequently compared with analogous analyses of SR6c3. Following the surface analyses, the FoldX Suite 4.0 (Schymkowitz et al., 2005), integrated as a YASARA View (Krieger and Vriend, 2014), was employed for the computation of the thermal stability of Nbs. To further investigate the grounds underneath the enhanced stability of ECSbs compared to SR6c3, we utilized ProDy’s InSty module (Bakan et al., 2011; Bakan et al., 2014; Zhang et al., 2021). Using default parameters (distance, angle, and cutoff distance), we analyzed seven types of intramolecular interactions, including hydrogen bonds (HBs), salt bridges (SBs), repulsive ionic interactions (RIB), hydrophobic interactions (HPh), disulfide bonds (DiBs), π−stacking (Pistacking), and π−cation (PiCation) interactions. Additionally, key residues with significant contributions to these interactions were identified, particularly those playing prominent roles in the interaction network of each Nb.

2.6 Protein-protein docking

To identify potential binding sites on the pre-processed RBD for Nbs (ECSbs and SR6c3), we initially employed an FFT-based global docking approach using PIPER within BioLuminate (Kozakov et al., 2006; Chuang et al., 2008; Fernández-Recio and Sternberg, 2010; Desta et al., 2020; Schrödinger, 2023c). This procedure generated 30 unique conformations for each RBD-Nb complex. Due to the inherent rigidity of PIPER’s docking methodology, each docked pose was pre-processed before calculating the binding free energy using the Molecular Mechanics Generalized Born Surface Area (MM/GBSA) method (Genheden and Ryde, 2015). To identify the most optimal pose between the Nbs and RBD, binding free energy calculations were performed, which subsequently provided the foundation for in-depth interaction analysis and subsequent MD simulations.

To visually validate the EC of the binding interfaces in Nb-RBD complexes, we utilized the APBS (Adaptive Poisson−Boltzmann Solver) (Jurrus et al., 2018) plugin within PyMOL version 3.1.1 (Schrodinger, 2024) to generate and visualize electrostatic potential maps. Molecular structures were prepared using PDB2PQR (Dolinsky et al., 2004; 2007) with default parameters. Electrostatic potentials were calculated over a range of ±2.0 kcal.mol−1·e, employing a solvent-excluded surface and a grid spacing of 0.50 Å.

The rigidity of FFT-based approach may limit conformational sampling; therefore, we subsequently employed the Rosetta flexible docking module (Gray et al., 2003; Wang et al., 2005; Chaudhury and Gray, 2008; Marze et al., 2018). This flexible protocol refines both backbone and side-chain conformations. Accordingly, all structures were meticulously prepared to ensure optimal backbone and side-chain integrity. Preparation involved reassigning chain IDs: Chain A was designated for RBD, previously labeled as Chain R, while Chain B was assigned to Nb, formerly designated as Chain A. The residues were then renumbered so that the RBD encompassed residues 1–186 and the Nb commenced from residue 187 onward. Subsequently cleaning each PDB file to preserve only standard atomic records, then assessing initial energetic properties with the Rosetta Energy Function 2015 (REF 2015) (Alford et al., 2017). To relieve steric clashes without disrupting the overall fold, a constrained FastRelax protocol was applied (Nivón et al., 2013). Afterward, interfacial residues were identified using a 6.0 Å distance criterion, enabling selective side−chain repacking confined to the interface via the PackRotamersMover with REF 2015. The prepacked structures were subsequently used for further steps.

Docking was conducted using the Rosetta protein docking module to generate 1000 decoys for each complex, wherein the RBD and Nbs were imported and subjected to FastRelax to achieve local geometric optimization. The optimized structures were subsequently assembled into a single pose by introducing a jump, preserving the native configurations derived from the PIPER-docked conformations. Interfacial residues were identified, and side-chain repacking was selectively confined to these regions to maintain the native-like conformation of the remaining protein structure, thereby minimizing structural perturbations outside the binding interface. A final relaxation phase was applied to further refine both side-chain and backbone orientations, resulting in a set of convergent RBD−Nb complexes. For each docking trajectory, the RMSD was calculated based on the heavy atoms of the interface residues, using the PIPER-docked poses as reference structures. This analysis was conducted utilizing the Rosetta InterfaceAnalyzerMover (Benjamin Stranges and Kuhlman, 2013). To characterize the docking funnel, the interface score (I_score) was plotted against the interface RMSD (I_rmsd), facilitating the identification of energetically favorable conformations.

The top-ranked decoys, exhibiting optimal binding energetics and structural stability, were selected for subsequent post−docking analyses. Post-docking analyses included the calculation of binding free energy (ddG), calculations of solvent-accessible surface area (SASA), and evaluations of Sc and EC. The ddG of each bound complex was recalculated using REF 2015 (Alford et al., 2017), and the change in SASA upon complex formation (dSASA) was determined. Shape and electrostatic complementarity (EC) were assessed using the Shape Complementarity Calculator and Electrostatic Complementarity Calculator, providing insights into the geometric fit and electrostatic compatibility at the interface. Additionally, interface residues were identified using the Interface Analyzer Mover in PyRosetta 4.0 (Chaudhury et al., 2010), which is configured to analyze designated chains by computing packing statistics and dSASA. Residues are classified as interface residues based on distance−based contact criteria, with both the total count and the specific set of residues defining the interface extracted from the analysis. Subsequently, the interaction types between these partnered residues are determined by generating a Pose object from the PDB structure and applying the Interface Analyzer Mover (Benjamin Stranges and Kuhlman, 2013) to evaluate intermolecular interactions. In this step, hydrogen bonding and van der Waals contributions were quantified using the HB and LJ Att/LJ Rep energy terms, respectively; salt bridges are identified from the Elec, while π-π, π-cation, and π-anion interactions are inferred from the side−chain geometry and relative orientation of residues. Finally, an interaction network was constructed by parsing dataset of residue-residue interaction data with Pandas and building a graph using NetworkX, where residues serve as nodes and their interactions as edges.

2.7 Molecular dynamics simulations

MD simulations were conducted to assess the stability of I_score based top ranked decoys of each Nbs-RBD complex using Desmond (Banks et al., 2005) software within the Schrodinger suite. Nbs-RBD complexes were prepared, optimized, and minimized using the PPW (Madhavi Sastry et al., 2013; Schrödinger, 2023e) in Maestro (Schrödinger, 2023b). The System Builder tool was employed to prepare all systems, utilizing the OPLS-2005 (Bowers et al., 2006) force field with the TIP3P (Hornak et al., 2006) solvation model in an orthorhombic box. To achieve the physiological equilibrium, the systems were neutralized by introducing 0.15 M sodium chloride. The simulations were executed at 1 atm and 310.0 K using Martyna-Tobias-Klein barostat (Martyna et al., 1994) and Nose-Hoover thermostat in an NPT ensemble (Nosé, 1984; Hoover, 1985). Prior to the simulation, the models underwent a relaxation process, and the stability of these complexes was analyzed during a 100 ns MD production run. Trajectories were saved at 100 ps intervals, and the simulation’s stability was assessed by comparing the Nb-RBD’s root mean square deviation (RMSD) and root mean square fluctuation (RMSF) over time (Maiorov and Crippen, 1994).

2.8 MM/GBSA analysis

Following MD simulations, binding free energy (ΔG) calculations were performed at intervals of 20 ns using the Prime (Jacobson et al., 2002; 2004; Schrödinger, 2023d) on equilibrated MD trajectories to quantify the thermodynamic drivers of Nb-RBD interactions. The total ΔG was computed as:

ΔG=GComplexGRBD+GNb

where GComplex, GRBD, and GNb represent the free energies of the bound complex, isolated RBD, and isolated Nb, respectively. The MM/GBSA energy decomposition dissects ΔG into six key components:

ΔGCoulomb quantifies electrostatic interactions between charged residues, calculated using Coulomb’s law, capturing the influence of charge−charge interactions on ΔG. ΔGvdW encompasses van der Waals contributions, including both attractive dispersion forces and repulsive Pauli exclusion effects, reflecting short−range atomic interactions. ΔGSolvGB quantifies polar solvation energy, modeled through the Generalized Born approximation to account for solvent−mediated electrostatic shielding and desolvation penalties. ΔGLipo measures hydrophobic burial energy, proportional to the reduction in solvent−accessible surface area (SASA) of non−polar residues, capturing the entropic gain from solvent exclusion. ΔGHbond denotes energy contributions from hydrogen bond formation, derived from geometric and electrostatic criteria, highlighting the role of specific polar interactions. Finally, ΔGCovalent represents internal strain energy arising from bonded interactions (bonds, angles, dihedrals) within the Nb or RBD during binding, reflecting conformational changes induced by complex formation.

To account for entropic effects, a quasi−harmonic entropy analysis was integrated into the workflow. Conformational entropy changes (ΔSbinding) were computed from MD trajectories using Maestro’s Simulation Interaction Analysis tool (Schrödinger, 2023b), which employs a quasi−harmonic approximation to estimate vibrational entropy:

ΔSbinding=ScomplexSRBD+SNb

The total binding free energy was then thermodynamically refined using:

ΔG=ΔGMM/GBSATΔSbinding

where T is the simulation temperature (300 K).

3 Results

3.1 RBD refinement and analysis

A thorough review using Protein Structure Reliability reports identified RBD (7VYR_R) as the most suitable candidate due to its minimal structural defects. This RBD structure was further refined by removing ions and crystallographic water molecules located beyond 5 Å, introducing disulfide bonds, and capping termini to get stable and correct conformation. Hydrogen bond networks were optimized through the prudently incorporation of hydrogens and the rotation of the amide (Asn, Gln, and His) and hydroxyl (Ser, Thr, and Tyr) functional groups. Energy minimization improved peptide planarity for specific bonds (C336−P337, E340−V341, T345−R346, N360−C361, and A388−D389) and reduced temperature factor for L518, H519, and A520. Loop optimization improved backbone dihedral angles (Φ and Ψ) for A372, K386, F392, D442, and Y505 (Supplementary Table S1). These refinements effectively addressed the structural issues present in RBD structure (Figure 2).

Figure 2
www.frontiersin.org

Figure 2. Structural representation and evaluation of the RBD facilitating ECSb nanobody construction. (A) 3D structure of RBD (red) retrieved from protein data bank underwent structural refinement protocols resulting in a refined RBD structure (blue; (B) with fewer structural issues. (C) A visual comparison of the two structures, highlighting the improvements achieved through the refinement process. Both models fitted well within electron density maps of 7VYR_R.

3.2 Epitopes characterization and electrostatic analysis

The analysis of titration curves over a pH range of 2–12 to determine the pKa values of ionizable residue in each epitope was conducted (Supplementary Figure S1). In addition, PROPKA [59–61] was used to visualize the pKa values of epitope residues at a pH of 7.2 ± 0.2 (Supplementary Table S2). The epitopes were characterized by their number of ionizable residues: CA1 had the most (9), followed by AS3 (6), AS1 (4), AS2 (2), and CA2 (2). This analysis suggested that CA1 would have the strongest and the most favorable EC, followed by AS3 and AS1.

Surface analyses assessed the overall electrostatic charge of the epitopes and RBD. CA1 and AS3 exhibited the highest positive charge, while AS1 demonstrated the most negative charge. These findings established the basis for recommendation for replacing specific charged amino acids in FRs with those selected for mutation from the prior structure.

3.3 ECSbs construction

This study presents a novel approach to FR engineering, building upon the SR6c3 Nb for an integrated strategy (Chen X. et al., 2021). This approach uniquely integrates CDR design with the optimization of the overall electrostatic potential between the epitope and the RBD. Unlike traditional FR engineering approaches that primarily focus on stability, solubility, and immunogenicity (Lawrence et al., 2007; Miklos et al., 2012), this approach empowers FRs to actively contribute to enhanced binding affinity and specificity in conjunction with the CDRs. We strategically engineered specific amino acids in both the FRs and CDRs of the ECSb sequences to minimize aggregation propensity, induce complementary dipole, and optimize electrostatic potential (Figure 3; Supplementary Table S3). Balancing net charges of CDRs and entire ECSbs with their corresponding epitopes and RBD, led to improved ECSb binding to the RBD. To specifically isolate the effect of EC on binding affinity and specificity, a minimal number of hydrophobic residues were introduced, as described earlier (Chothia and Janin, 1975; Rose et al., 1985; Nicholls et al., 1991).

Figure 3
www.frontiersin.org

Figure 3. Framework regions (FRs; upper section) and complementarity−determining regions (CDRs; lower section) sequences displaying mutations in ECSbs sequences generated through integrated engineering strategy. The mutations introduced into the FRs (FR1, FR2, FR3, and FR4) of ECSbs sequences to minimize aggregation propensity are outlined in red and mutations to induce complementary dipole are outlined in blue. The CDR sequences of ECSb variants are aligned and compared with the SR6c3 reference sequence, illustrating the modifications made to optimize binding affinity and electrostatic potential. Asterisks under residues (F37, E44, R45, and V48) represent the VHH hallmark residues.

All Nbs (ECSbs and SR6c3) demonstrated strong nativeness (overall) scores, with values consistently near the threshold level of 0.8, ranging from 0.76 to 0.79 (Figure 4). ECSb2 and SR6c3 exhibited the highest scores (0.79), while ECSb3 had a lower score (0.76). Notably, the FR sequences of all Nbs consistently demonstrated high nativeness, surpassing the threshold with an average score of 0.93, due to their high sequence conservation. In contrast, the CDRs exhibited significantly lower nativeness scores (<0.8), ranging from 0.12 to 0.6, reflecting their completely randomized sequences (Figure 4).

Figure 4
www.frontiersin.org

Figure 4. Humanization profile of nanobodies evaluated by AbNatiV. The bar chart comparing the computational humanization profile of the whole nanobodies (ECSbs and SR6c3) and their overall score, CDR score, and FR scores. The x−axis represents the names, and the y−axis represents the score.

These results demonstrate the potential of our novel approach in Nb engineering, balancing EC with structural integrity. The high nativeness scores of the FRs, coupled with the strategically designed CDRs, suggest that our method successfully integrated FR engineering with CDR design while maintaining overall antibody structure.

3.4 Structure prediction analysis

All predicted structures of Nbs demonstrated high confidence, with pLDDT score consistently exceeding 85 (Supplementary Figure S2). ECSb2 demonstrated the highest average pLDDT scores of 93.6, followed by SR6c3 at 92.9, ECSb5 at 90.5, ECSb1 at 89.6, ECSb4 at 88.8, and ECSb3 at 87.8. Each Nb exhibited distinct variations in pLDDT scores, suggesting distinct structural characteristics within their flexible regions. Consistent with other protein structure prediction tools (Du et al., 2021; Su et al., 2021; Cohen et al., 2022; Lin et al., 2022; 2023; Wang W. et al., 2022; Wu et al., 2022; Ruffolo et al., 2023), pLDDT scores were generally lower in loop regions but remained high (>95) in the non−loop region (Supplementary Figure S2).

3.5 Structural refinement and analyses

Structural refinements of the Nbs were performed using Prime (Jacobson et al., 2002; 2004; Schrödinger, 2023d). However, ECSb3 and ECSb4 exhibited persistent structural anomalies, including distorted backbone and sidechain dihedrals, deviations in bond angles and lengths, non−planer peptides, and the buried unsatisfied hydrogen bond acceptors. These issues were rectified through additional energy minimization and loop refinement, though this process was challenging due to the exacerbating factors such as temperature, bond angle and length deviations, steric clashes, and sidechain dihedrals. The remaining Nb structures demonstrated satisfactory reliability, as validated by PSRs and Ramachandran Plots (Supplementary Table S4; Supplementary Figure S3).

Surface evaluations of the refined structures were conducted using BioLuminate (Beard et al., 2013; Salam et al., 2014; Zhu et al., 2014; Schrödinger, 2023a) to assess their electrostatic characteristics, aggregation tendency, and reactivity. Aggregation score calculations identified SR6c3 as having the highest aggregation propensity (190.548), followed by ECSb5 (70.451), while ECSb3 exhibited the lowest score (0.853) (Table 1). This observation implies that SR6c3 and ECSb5 have a higher propensity for aggregate formation and can lead to structural misfolding during in vitro or in vivo development compared to the remaining Nbs (Lijnzaad et al., 1996; Nissley et al., 2022; Van Gils et al., 2022; Peruzzi et al., 2024). Patch analysis further elucidated the basis of aggregation scores by examining the distribution, intensity, and orientation of hydrophobic patches (Sankar et al., 2018) (Supplementary Figure S4C; Supplementary Figure S5). SR6c3, with its 85% hydrophobic amino acids, displayed five hydrophobic patches with an average intensity of 0.639. In contrast, ECSb1 and ECSb3 had fewer patches, 6 and 4 with average intensity of 0.543 and 0.503, respectively. Regarding distribution and orientation, patches in SR6c3 are clustered majorly around CDRs, driven by the substantial involvement of hydrophobic residues, which facilitate orientation changes due to CDR flexibility. Conversely, ECSb1 and ECSb3 had highly scattered hydrophobic patches along FRs, hindering aggregate formation.

Table 1
www.frontiersin.org

Table 1. Comparison of the sum of aggregation score (AggScore) and thermostability between ECSbs and SR6c3.

Our electrostatic engineering strategy increased the number and intensity of negatively charged patches on ECSbs, enhancing their solubility potential and reducing aggregation propensity (Kramer et al., 2012; Ganesan et al., 2016; Carballo-Amador et al., 2019). Notably, ECSbs exhibited an average of 12 patches with an intensity of 1.26, while SR6c3 exhibited 15 scattered patches with a lower intensity of 1.10 (Supplementary Figures S4A, B; Supplementary Figure S5). The higher average intensity of ECSbs is due to distribution of 64 residues on average in 12 patches compared to SR6c3’s 64 residues distribution in 15 patches (Supplementary Figure S6).

Assimilation of more charged residues; 13 residues in positive patches and 12 residues in negative and four residues in hydrophobic patches in ECSbs’ CDRs result in higher electrostatic steering potential than SR6c3 (Wade et al., 1998; Patil and Nakamura, 2007; Tsutakawa et al., 2017) (Supplementary Figure S4; Supplementary Figure S7). As in contrast, SR6c3 CDRs have 5 positively charged residues, 13 negatively charged, and 10 hydrophobic residues. This finding aligns with the above results of higher aggregation around CDRs in SR6c3 (Supplementary Figure S5; Supplementary Figure S6).

Our integrated approach successfully developed electrostatically complementary patches and potential between the Nbs and specific epitopes, as well as the overall RBD. ECSbs, such as ECSb1, exhibited a balance distribution of positively (16) and negatively charged (14) residues, complementing AS1’s net charge of ‒2e. Across the full ECSb1 structure, 59 residues contribute to positively charged patches and 67 to negatively charged patches, aligning with the overall RBD charge of +2e (Supplementary Figure S6). Similar pattern was consistent across all ECSbs.

BioLuminate (Beard et al., 2013; Salam et al., 2014; Zhu et al., 2014; Schrödinger, 2023a) was also employed to predict reactive hotspots in the Nbs structures, focusing on vulnerabilities to post−translational modifications (PTMs), including cysteine modifications, glycosylation, deamination, oxidation, proteolysis, and Asp isomerization (Farnsworth et al., 1989; Cohen, 2002; Blom et al., 2004; Ohtsubo and Marth, 2006; Aebersold et al., 2018; March and Smith, 2020). These PTMs can significantly impact protein stability, immunogenicity, and function (Tsien, 1998; Hermeling et al., 2004; Walsh, 2010; Kumar et al., 2016; Lee et al., 2023). While no N− or O−linked glycosylation sites were identified in any of the Nbs, several reactive PTM hotspot, including deamidation, oxidation, proteolysis, Aspartate isomerization, and free cysteine sites, were found, with SR6c3 exhibiting the highest number (38) and ECSb3 the fewest (25) (Supplementary Table S5). Proteolysis and oxidation were the predominant PTMs across all Nbs, with SR6c3 showing the highest number of oxidation hotspots (14), indicating a higher susceptibility to oxidative damage in comparison to the other Nbs. In contrast, ECSbs displayed a higher number of proteolytic sites, suggesting increased vulnerability to proteolytic degradation compared to SR6c3.

Thermostability predictions using FoldX (Schymkowitz et al., 2005) revealed significantly higher thermostability in the ECSbs compared to SR6c3, with ECSb4 demonstrating the highest thermostability at 148.3 kcal.mol−1, while SR6c3 had the lowest (62.6 kcal.mol−1) (Table 2). Further analysis of intramolecular interactions, including hydrogen bonds (HBs), salt bridges (SBs), repulsive ionic bonding (RIB), disulfide bonds (DiBs), π−π stacking (Pistacking), π−cation (PiCation) interactions, and hydrophobic interactions (HPh), revealed that ECSbs possess superior characteristics. Electrostatic engineering resulted in a greater number of SBs in ECSbs, which contain an average of seven SBs compared to four in SR6c3 (Supplementary Figure S7). Moreover, SR6c3 possesses four RIBs, while ECSbs typically have only one RIB per structure. SR6c3 also lacks both π−π stacking and π−cation interactions, which significantly contribute to changes in free energy (ΔΔG), another factor of its lower stability. Furthermore, the presence of more DiBs in ECSbs, as compared to SR6c3, enhances their overall stability. In terms of HBs, ECSb5 and ECSb2 exhibit the highest number of interactions, with 33 and 32 HBs respectively, surpassing SR6c3, which has 25 HBs (Supplementary Figure S7). In contrast, SR6c3, due to a greater number of hydrophobic residues, has a higher propensity for aggregation, resulting in 32 HPh interactions compared to an average of 28 HPh interactions in ECSbs.

Table 2
www.frontiersin.org

Table 2. Comparison of the binding free energy (ΔG), Coulombic component of binding free energy (ΔGC), hydrogen bonding component of binding free energy (ΔGHB) between ECSbs and SR6c3.

A total of 10 key residues contributing to Nb stability were identified due to their prominent involvement in the intramolecular interaction landscape (Figure 5). While the precise functions of these residues remain to be fully elucidated, their potential roles can be inferred. Cysteine residues likely stabilize the nanobody through disulfide bond formation. Arginine, a positively charged amino acid, may form salt bridges, contributing to stability and antigen binding. Aromatic residues like tyrosine can participate in various interactions, including pi−pi stacking and hydrogen bonding. Hydrophobic residues like valine often stabilize the protein’s core. However, further research is needed to fully understand the specific roles of these residues in nanobody function.

Figure 5
www.frontiersin.org

Figure 5. 3D structures of the all the nanobodies, (A) ECSb1, (B) ECSb2, (C) ECSb3, (D) ECSb4, (E) ECSb5, and (F) SR6c3 highlighting key residues contributing to their stability. The top 10 critical residues for each nanobody are listed in the boxes. Structural regions are color−coded as follows: purple for CDR1, green for CDR2, orange for CDR3, and red for key stability residues. The background density maps provide a visual context for the Nb structure.

The integrated electrostatic engineering approach effectively increased the number of SBs in ECSbs. Notably, mutations such as S50R, S61D, and W106R in ECSb2; G26D, A50D, N59R, Y95D, and W110D in ECSb3; and A23E, S25R, Q39D, S63D, and W107R in ECSb5 were instrumental in SB formation. Conversely, the two RIB interactions present in SR6c3 (W105−R45 and Y95−R45) are attributed to residues W105 and Y95, which are mutated in ECSbs. The comparison of thermostability results with the intramolecular contact analysis of Nbs elucidated one of the contributing factors to the heightened stability in ECSbs compared to SR6c3 (Vogt et al., 1997; Kumar et al., 2000). These computational analyses indicated that the electrostatic integration strategy was successful in enhancing the thermostability in ECSbs compared to SR6c3.

3.6 Protein−Protein docking

MM/GBSA calculations was used to rank the 30 poses of each Nb-RBD docked system based on their binding free energy. Notably, ECSb5-RBD pose (P) 15 exhibited the highest binding free energy (ΔG) of 90.094 kcal.mol−1, closely followed by P26 of ECSb1-RBD at 88.911 kcal.mol−1. SR6c3-RBD ranked fourth with ΔG of 70.525 kcal.mol−1, while ECSb3-RBD P14 exhibited the lowest ΔG at 57.452 kcal.mol−1 (Table 2).

Each Nb in the docked pose with the highest ΔG exhibited a distinct binding mode relative to the whole spike protein (Figure 6). ECSb1 and ECSb3 bound to a cryptic site within the RBD, whereas ECSb2, ECSb4, and SR6c3 targeted the receptor-binding motif (RBM) region. Detailed analysis of ECSb1 and ECSb2 interactions revealed unexpected binding patterns. Contrary to their intended targets (AS1 and AS2), these Nbs demonstrated a preference for residues that offered a more favourable binding architecture. This unexpected behaviour was characterized by the formation of two hydrogen bonds, resulting from engineering within the FRs. ECSb1 bound to AS3 instead of AS1, whereas ECSb2 bound with CA1 instead of AS2. In contrast, other ECSbs exhibited specific interactions that aligned with their target antigenic sites, highlighting the binding specificity conferred by electrostatic complementation (Figure 6).

Figure 6
www.frontiersin.org

Figure 6. Specific interactions between RBD (A) and the Nbs; ECSb1 (B), ECSb2 (C), ECSb3 (D), ECSb4 (E), ECSb5 (F), and SR6c3 (G) are highlighted in the circles.

Visualization of the electrostatic potential maps for the binding interfaces of each Nb-RBD complex corroborated the results of calculated EC scores (Supplementary Figure S8). The electrostatic maps of the ECSb2-RBD interface revealed similar charge distributions, leading to repulsive interactions that resulted in the lower EC score−. In contrast, ECSb1-RBD and ECSb5-RBD complexes displayed both high EC and Sc, resulting from strong molecular attractions. For the SR6c3-RBD complex, the electrostatic map of SR6c3 exhibited a high negative charge density in the upper region and a high positive charge density in the lower right region, whereas the RBD displayed a slight negative charge density in the upper portion and a high positive charge density in the lower right region (Supplementary Figure S8). This analogous charge distribution reduced the EC score for the SR6c3-RBD complex.

In the second stage of docking, 1,000 decoys were generated for each complex using Rosetta, and the docking funnel was characterized by plotting interface scores (I_score, in Rosetta Energy Units, REU) against interface RMSD (I_rmsd). Among the engineered nanobodies, ECSb4-RBD exhibited the most favorable binding energy, with the lowest I_score of −133.232 REU, followed by ECSb2-RBD (−121.634 REU). In contrast, ECSb1-RBD showed the least favorable I_score (−32.619 REU) among the top-scoring decoys (Figure 7). The reference SR6c3-RBD demonstrated an intermediate I_score of −79.541 REU. For each complex, the decoy with the lowest I_score was extracted from the docking trajectory, representing the energetically optimal conformation for subsequent structural and functional analyses (Figure 9). These results highlight the enhanced interface energetics of ECSb4 and ECSb2 compared to both SR6c3 and other ECSb variants.

Figure 7
www.frontiersin.org

Figure 7. Docking funnels for each Nb-RBD complex. The interface scores (REU) are plotted against I_rmsd (Å) for each complex, (A) ECSb1-RBD, (B) ECSb2-RBD, (C) ECSb3-RBD, (D) ECSb4-RBD, (E) ECSb5-RBD, and (F) SR6c3-RBD. Red colored data points have the lowest interface scores.

Among these structural and functional analyses, the comprehensive scoring elucidated the detailed energetic breakdown of Nb-RBD complexes. The major energetic contributions include Lennard-Jones attractive (LJ Att) and repulsive (LJ Rep) components, in conjunction with the Lazaridis-Karplus solvation energy (LK Sol) and Coulombic electrostatics (Elec), as well as hydrogen bond (HB) energies, collectively informed the calculation of the overall binding free energy (ddG). Among the complexes, ECSb1-RBD displayed the most favorable ddG of −802.83 REU, driven by robust LJ Att. of −1642.61 REU and HB energy of −280.36 REU (Table 3). Coupling with the highest dSASA 15,285.66 Å2, ECSb1--RBD’s scoring analysis yielded a normalized ddG/dSASA ratio of −0.0525. Conversely, although ECSb4-RBD was identified as possessing the lowest I_score (−133.232 REU) in the docking funnel, scoring analysis showed ddG (−592.35 REU) and dSASA (13,852.21 Å2) resulted in a less favorable normalized ddG/dSASA ratio of −0.04276 (Table 3). This lower ddG evidenced by higher LK Sol (973.46 REU) and LJ Rep (372.90 REU) coupled with lower individual energetic contributions of LJ Att; −1617.54 REU, Elec: −390.23 REU, and HB: −246.08 REU.

Table 3
www.frontiersin.org

Table 3. Key scoring components including Lennard-Jones attractive (LJ Att), repulsive (LJ Rep), and solvation (LK Sol) energies, electrostatic contribution (Elec), hydrogen bonding (HB), and the resulting binding free energy change (ddG) per solvent−accessible surface area (dSASA) of Rosetta Energy Function 2015.

While ECSb5-RBD demonstrated robust binding energetics as well (ddG = −718.91 REU) yielding a normalized ddG/dSASA ratio of −0.0507 coupled with dSASA of 14,171.11 Å2 (Table 3). The second highest ddG was due to the major contributions of highest LJ Att (−1680.14 REU) and HB (−288.20 REU) energetics even with the highest LJ Rep of 395.94 REU and moderate score of LK Sol energy (975.81 REU). On the other hand, SR6c3-RBD showed intermediate ddG (−615.84 REU) even with lowest Coulombic (Elec: −357.99 REU) and hydrogen bonding (HB: −241.33 REU) energetics (Table 3). Although, it has lowest yet balanced LJ Att and LJ Rep −1539.38 REU and 363.17 REU, respectively. However, the lowest solvation penalty (LK Sol: 852.80 REU) and second highest buried surface area (15,117.49 Å2) improve ddG by requiring less energetic cost for desolvating the binding interface and increasing buried surface area to amplify the net favorable hydrophobic interactions. These data suggest that while Nb-RBD complexes may offer the most energetically favorable interactions, the interplay between individual energy components and the solvent exposure upon complex formation remains critical for the holistic assessment of binding affinity.

The structure of interfacial interaction networks coupled with Sc and EC further characterized factors underneath these scorings. The ECSb-RBD complexes demonstrated distinct interfacial binding profiles compared to the SR6c3-RBD, as revealed by Sc and EC analyses. ECSb3-RBD achieved the highest Sc (0.7304), attributed to its extensive network of geometrically constrained interactions, including 17 hydrogen bonds (e.g., T10-D236) and five π-π stacking interactions (e.g., R11-F223) (Table 4; Figure 8). These interactions demand precise spatial alignment, enabling exceptional structural fit. In contrast, SR6c3-RBD exhibited moderate Sc (0.6302), limited by narrow interactions network (14 hydrogen bonds, five π−π) and the absence of salt bridges (Supplementary Table S6). Notably, ECSb1-RBD, despite its small interface (18 residues), displayed intermediate Sc (0.6851), supported by six hydrogen bonds and a stabilizing salt bridge (174R-298D), underscoring the impact of targeted polar interactions on structural alignment (Table 4; Figure 8).

Table 4
www.frontiersin.org

Table 4. Interface characteristics including shape complementarity (Sc), electrostatic complementarity (EC) and number of interface residues for each Nb-RBD complex.

Figure 8
www.frontiersin.org

Figure 8. Schematic representation of interfacial interactions in the Nb-RBD complexes. The network of interactions for ECSb1-RBD (A), ECSb2-RBD (B), ECSb3-RBD (C), ECSb4-RBD (D), ECSb5-RBD (E), and SR6c3-RBD (F), with Nb residues shown as blue balls and RBD residues as red balls. The various types of interactions are color−coded as follows: Hydrophobic interactions (dark pink), π-Anion interactions (purple), Hydrogen bonds (green), π-π stacking (blue), Repulsive Ionic Bonding (red), Salt bridges (orange), and π-Cation interactions (grey).

Electrostatic optimization varied markedly across complexes. ECSb5-RBD exhibited the highest EC (0.5295), driven by strategic integration of salt bridges (N8−R236) and π-cation interactions (T10−D249), which enhanced charge complementarity despite its modest interface size (31 residues). Conversely, SR6c3-RBD showed the lowest EC (0.2325), a consequence of its reliance on non-polar interactions (28 hydrophobic, five π-π) and lack of salt bridges (Table 4; Figure 9). ECSb4-RBD, despite possessing the largest interface (67 residues), displayed poor EC (0.3047), as its interaction profile was dominated by hydrophobic clustering (39 instances), which lacks electrostatic specificity (Table 4).

Figure 9
www.frontiersin.org

Figure 9. Structural and stability analysis of Nb-RBD complexes. (A) Solvent−accessible surface area (SASA) over time, highlighting the extent of exposure of hydrophobic and polar regions for the ECSb1-RBD (black), ECSb2-RBD (green), ECSb3-RBD (blue), ECSb4-RBD (cyan), ECSb5-RBD (purple), and SR6c3-RBD (yellow) complexes, (B) Root−mean−square deviation (RMSD) over time, reflecting the overall structural deviation during the simulations, and (C) Root−mean−square fluctuation (RMSF) analysis, indicating the flexibility of residues within each complex across the simulation, with higher RMSF values corresponding to greater flexibility. The data provide insights into the structural stability and conformational changes of the complexes throughout the simulation period.

The interplay between Sc and EC revealed design-driven trade-offs. While ECSb3-RBD achieved structural precision (Sc = 0.7304), its moderate EC (0.3901) suggested that geometric optimization alone insufficiently resolves electrostatic mismatches. In contrast, ECSb5-RBD prioritized charge complementarity, achieving superior EC (0.5295) through targeted electrostatic interactions, even with fewer total interactions (23 vs. SR6c3’s 47). Intermediate performers included ECSb2-RBD, with moderate EC (0.4570) linked to four salt bridges and three π-anion interactions, and ECSb4-RBD, whose low EC (0.3047) reflected over-reliance on hydrophobic contacts.

Collectively, ECSbs outperformed SR6c3-RBD in EC, with ECSb5-RBD achieving a 2.3-fold improvement. The SR6c3-RBD’s suboptimal EC underscored the limitations of natural systems in balancing hydrophobic packing with electrostatic optimization − a challenge mitigated in the ECSb series through rational design of interaction networks. These results highlight the efficacy of tailored electrostatic synergy (ECSb5) and structural precision (ECSb3) in advancing molecular recognition beyond natural benchmarks.

The comprehensive analysis of docking energetics, interface complementarity, and binding affinity reveals distinct mechanistic advantages and trade−offs among the ECSbs compared to the SR6c3. ECSb4-RBD, despite exhibiting the most favorable initial interface energy, demonstrated suboptimal ddG and EC. This disparity arises from its over-reliance on hydrophobic clustering (39 instances) and high solvation penalties, which offset favorable interfacial contacts. Conversely, ECSb1-RBD achieved the most favorable ddG, driven by robust LJ Att and HB, despite its small interface (18 residues). Its superior normalized ddG/dSASA ratio underscores efficient energy utilization per unit of buried surface area, highlighting the importance of packing density in hydrophobic−driven interactions.

ECSb3-RBD emerged as a paradigm of structural precision, with the highest Sc of 0.7304 due to its dense HB network (17) and π-π stacking (5). However, its moderate EC reflects unresolved electrostatic mismatches, illustrating that geometric optimization alone cannot fully compensate for charge misalignment. In contrast, ECSb5-RBD prioritized electrostatic synergy, achieving the highest EC through strategic salt bridges and π−cation interactions, which enabled a competitive ddG and normalized ratio despite fewer total interactions (23 vs. SR6c3’s 47). This underscores the efficacy of targeted charge complementarity in enhancing binding efficiency.

The SR6c3-RBD, while leveraging hydrophobic packing (28 instances) and moderate shape fit, exhibited the lowest EC and intermediate ddG. Its performance highlights inherent limitations in natural systems: a reliance on broad, spatially tolerant interactions limits electrostatic optimization, despite advantages in solvation cost and buried surface area.

3.7 MD simulations

The molecular dynamics (MD) simulations of Nb-RBD complexes revealed critical insights into their structural stability, interfacial dynamics, and solvent accessibility, elucidating the biophysical underpinnings of their binding efficacy. ECSb5-RBD emerged as the most stable complex, exhibiting the lowest average RMSD (3.72 Å) and RMSF (1.41 Å), indicative of a rigid, well−optimized interface (Figures 9B, C). This exceptional stability likely stems from its strategic integration of electrostatic interactions (e.g., salt bridges, π-cation pairs) and a high proportion of secondary structure elements (35.96% SSE, predominantly β-strands), which collectively minimize conformational entropy and resist structural deviations. The low RMSF, particularly at the binding interface, underscores restricted residue-level flexibility—a hallmark of geometrically precise and electrostatically complementary interactions, as previously identified in its high Sc and EC scores.

In stark contrast, ECSb1-RBD displayed the highest RMSD (7.74 Å) and RMSF (2.39 Å), reflecting significant conformational plasticity (Figure 9C). While this dynamic behavior suggests ongoing structural adjustments to optimize binding, the paradoxically low SASA (9,981.11 Å2) indicates a deeply buried interface dominated by hydrophobic clustering (Figure 9A). This duality implies that ECSb1-RBD’s binding mechanism prioritizes strong, enthalpically favorable van der Waals interactions and hydrogen bonds (HB: −280.36 REU) over rigid structural alignment, allowing it to maintain affinity despite flexibility (Table 3). However, the elevated solvation penalty (LK Sol: 973.46 REU) observed in prior scoring analyses likely exacerbates its conformational instability, as desolvation costs destabilize the complex over time.

ECSb2-RBD and ECSb3-RBD exhibited intermediate stability (RMSD: 3.81–3.94 Å) but distinct interfacial characteristics (Figure 9B). ECSb2-RBD’s high SASA (16,594.02 Å2) suggests a solvent−exposed interface, which may enhance adaptability through dynamic interactions but compromises hydrophobic stabilization. Conversely, ECSb3-RBD’s moderate SASA (15,998.14 Å2) aligns with its geometrically constrained network of hydrogen bonds and π-π stacks, balancing rigidity with partial solvent accessibility (Figure 9A). Notably, both variants lack the electrostatic optimization seen in ECSb5-RBD, underscoring the trade-off between structural precision and charge complementarity.

ECSb4-RBD, despite its expansive interface (67 residues), displayed suboptimal stability (RMSD: 4.13 Å) and moderate RMSF (1.59 Å), consistent with its over-reliance on hydrophobic clustering (39 instances) and poor EC (0.3047) (Figure 9). The lack of polar specificity likely permits minor interfacial fluctuations, as non-directional hydrophobic interactions fail to fully constrain residue mobility. This aligns with its elevated solvation penalties (LK Sol: 973.46 REU) and repulsive forces (LJ Rep: 372.90 REU), which collectively destabilize the complex despite favorable initial interfacial energy (Table 3).

The SR6c3-RBD exhibited intermediate metrics across all parameters (RMSD: 4.44 Å; RMSF: 1.94 Å; SASA: 16,048.02 Å2), reflecting its developmental honed balance of hydrophobic packing and moderate structural adaptability (Figures 9A–C). Its relatively low solvation penalty (LK Sol: 852.80 REU) highlights artificial selection’s preference for energetically efficient desolvation and broad interfacial contacts (Table 3). However, its limited EC and sparse salt bridges reveal inherent limitations in natural systems, which prioritize robustness over precision.

3.8 MM/GBSA analysis

The MM/GBSA methodology (Genheden and Ryde, 2015), applied to equilibrated MD trajectories, provides a comprehensive thermodynamic perspective on the binding free energies (ΔG) of ECSbs and SR6c3. This analysis reveals distinct binding mechanisms and highlights the potential of this EC based design of ECSbs.

ECSb1-RBD exhibits moderate binding affinity (ΔG: −51.25 kcal·mol−1), primarily driven by strong van der Waals interactions (ΔGvdW: −56.25 kcal·mol−1) and lipophilic stabilization (ΔGLipo: −20.41 kcal·mol−1) (Table 5; Supplementary Table S7). These contributions reflect its deeply buried interface (low SASA: ∼9,981 Å2) and efficient hydrophobic packing density, which is further evidenced by transiently favorable ΔG values at 80 ns (−74.10 kcal·mol−1), correlating with maximal lipophilic burial (Supplementary Table S7). However, polar solvation penalties (ΔGSolvGB: +10.40 kcal·mol−1) and variable Coulombic contributions (ΔGCoulomb: +12.53 kcal.mol−1) partially offset these gains, aligning with its high desolvation costs (LK Sol: 973.46 REU) noted earlier (Supplementary Table S7; Table 3). The significant conformational flexibility observed in MD simulations (RMSD: 7.74 Å) suggests a substantial entropic penalty, which likely reduces binding affinity (Figure 10). Despite these trade−offs, ECSb1-RBD’s exceptional normalized energy density (ddG/dSASA: −0.0525) highlights its efficiency in leveraging hydrophobic-driven interactions, albeit at the expense of dynamic adaptability (Table 3).

Table 5
www.frontiersin.org

Table 5. Binding free energies (kcal·mol⁻1) for each Nb-RBD complex, measured at 0, 20, 40, 60, 80, and 100 ns, along with the average (Avg) binding energy across the simulations.

In similar, ECSb5-RBD achieves a moderate average ΔG (−64.39 kcal.mol−1), underpinned by ΔGvdW of −93.95 kcal·mol−1 and ΔGLipo of −50.31 kcal.mol−1 contributions (Table 5; Supplementary Table S7). Its structural rigidity (RMSD: 3.72 Å) minimizes entropic penalties, allowing robust hydrophobic network−powered enthalpic gains to dominate. However, polar solvation (ΔGSolvGB: +24.33 kcal·mol−1) and electrostatic interactions (ΔGCoulomb: +16.46 kcal.mol−1) despite highest EC (0.5295) driven by salt bridges and π-cation pairs remains unfavorable (Supplementary Table S7). Notably, its low solvation penalties and high secondary structure content (SSE: 35.96%) minimized entropic costs, reflecting its rigid, pre-organized binding interface. This enthalpic precision underscores the success of rational design in optimizing electrostatic synergy.

ECSb4-RBD yielded highest ΔG (−182.58 kcal.mol−1), yet the decomposition of energy terms reveals critical weaknesses (Table 5). The Coulombic contribution (ΔGCoulomb: +49.02 kcal.mol−1) is destabilizing, aligning with the low EC (0.3047) and underscoring poor electrostatic optimization evidenced by interfacial repulsion (R174-R279 charge clash) (Supplementary Table S7; Table 4). Conversely, van der Waals (ΔGvdW: −116.65 kcal.mol−1) and lipophilic (ΔGLipo: −28.85 kcal·mol−1) terms dominate the stabilization, consistent with the hydrophobic−driven interface (39 instances). However, the solvation energetics (ΔGSolvGB: −29.15 kcal.mol−1) and negligible hydrogen−bond contributions (ΔGHbond: −0.92 kcal.mol−1) further highlight weak polar specificity (Supplementary Table S7). Temporal fluctuations in ΔG (e.g., ΔGCoulomb ranging from −11.72 to +225.25 kcal.mol−1 over 100 ns) reflect transient electrostatic repulsions and interfacial rearrangements, indicative of an electrostatic instability. Contrarily, the progressive decline in covalent binding energy (ΔGCovalent) from −13.22 to −106.33 kcal·mol−1 suggests strengthening interfacial cohesion over time (Supplementary Table S7). This analogous inclination between ΔGCoulomb and ΔGCovalent over simulation trajectory reflects the influence of electrostatic instability on ΔGCovalent. However, the conversed drift in dominate energetics of ΔGvdW (−138.66 to −59.38) and ΔGLipo (−81.03 to +20.97 kcal.mol−1) align with MD observations of rising RMSD (1.96 Å at 0.2 ns to 4.66 Å at 100 ns) (Supplementary Table S7; Figure 10).

ECSb2-RBD and ECSb3-RBD exhibited divergent profiles, reflecting their distinct interfacial strategies. ECSb2-RBD’s average ΔG (−66.94 kcal.mol−1) is driven by strong van der Waals (ΔGvdW: −114.96 kcal.mol−1) and lipophilic contributions (ΔGLipo: −97.91 kcal.mol−1) (Table 5; Supplementary Table S7). At 60 ns of ΔG calculations trajectory, ECSb2-RBD displayed a pronounced outlier of +516.50 kcal.mol−1 driven by anomalously high ΔGCovalent, ΔGCoulomb, and ΔGvdW (+103.81 kcal.mol−1, +207.60 kcal.mol−1, and +331.74 kcal.mol−1, respectively) contributions (Table 5; Supplementary Table S7). This likely reflects transient conformational artifacts or force field limitations. Excluding this outlier, the complex showed high ΔG (−183.63 kcal.mol−1), supported by robust ΔGvdW (−204.30 kcal.mol−1) and ΔGLipo (−100.47 kcal.mol−1) interactions. However, its elevated SASA (∼16,594 Å2) incurred solvation penalties (ΔGSolvGB: + 17.63 kcal·mol−1), while reduced EC (0.4570) limited ΔGCoulomb gains. This aligns with its intermediate stability (RMSD: ∼3.81 Å) and suggests a design compromise between adaptability and interfacial specificity. In contrast, ECSb3-RBD achieves the second highest average ΔG (−119.07 kcal·mol−1), supported by geometric precision (Sc: 0.7304) via extensive hydrophobic bonding (ΔGvdW: −129.73 kcal·mol−1) and π-π stacking (ΔGLipo: −47.67 kcal.mol−1). However, its electrostatic inefficiency (EC: 0.3901) manifested in suboptimal Coulombic contributions (ΔGCoulomb: +17.08 kcal.mol−1). Although its moderate entropy and solvation penalty (ΔGSolvGB: +47.86 kcal.mol−1) tempers affinity, ECSb3-RBD’s structural precision exemplifies the advantages of engineered hydrogen bond networks in enhancing binding efficacy.

Conversely, SR6c3-RBD showed stronger ΔG (−105.50 kcal.mol−1), driven by favorable vdW (ΔGvdW: −80.58 kcal.mol−1) and lipophilic terms (ΔGLipo: −56.59 kcal.mol−1) (Table 5; Supplementary Table S7). ΔGvdW and ΔGLipo partially offset by weak electrostatic contributions (ΔGCoulomb: +175.61 kcal.mol−1) reflecting the lowest EC (0.3047), which were frequently counteracted by polar desolvation (ΔGSolvGB: −96.37 kcal.mol−1) (Table 5; Supplementary Table S7). These calculation trajectories aligned with the development prioritization of hydrophobic packing over charge optimization, a limitation that ECSbs effectively addresses. The highest electrostatic penalty and entropic costs reflect SR6c3-RBD’s natural trade−off between adaptability and stability, further underscoring the advantages of rational design in achieving superior binding profiles.

4 Discussion

This study aimed to explore the potential of electrostatically tailored Nbs as a therapeutic strategy against SARS-CoV-2, focusing on generating high-affinity and stable binding with RBD of the spike protein. Our findings demonstrated the successful development and characterization of ECSbs designed to target the Spike protein’s RBD. By leveraging computational biology and engineering, we created Nbs with enhanced electrostatic profile and binding affinities (Cheloha et al., 2020).

The structural refinement of the RBD effectively addressed various anomalies inherent in the initial model, enhancing its overall accuracy and stability. Disulfide bond formation,− the introduction of hydrogen bonds, loop refinements and energy minimizations of selected residues were pivotal in achieving conformational stability. Key residues (such as C336−P337, T345−R346, and A388−D389) demonstrated improved peptide planarity and reduced thermal fluctuations, particularly around the critical loop regions (A372, K386, and Y505). A prior study has shown that different amino acid residues R403, K/N/T417, L455, F486, Y489, F495, Y501, and Y505 on RBD play a key role in the protein recognition mechanism (Biskupek and Gieldon, 2024). Our methodical approach to ECSb design, focusing on optimizing electrostatic interactions through strategic amino acid modifications, proved effective in enhancing binding affinity and specificity towards the RBD (Khan et al., 2021). This aligns with the previous study (Zhu et al., 2023) that successfully engineered Nbs with higher binding affinity to the S protein using a combination of template selection, mutation analysis, and single−site saturated mutagenesis in CDRs. However, our methodical approach for ECSbs represents a novel approach in the field of Nb engineering and differ principally, as our study focused on ECSb design, through strategic amino acid modification both in CDRs and FRs based on electrostatic complementation (Entzminger et al., 2017; Gąciarz and Ruddock, 2017).

The engineered Nbs in this study exhibited varying binding patterns to the SARS-CoV-2 spike protein. Three of the five Nbs (ECSb3-ECSb5) successfully bound to the targeted antigenic sites, demonstrating their specificity. However, ECSb1 and ECSb2 displayed unexpected binding behavior, binding to alternative sites on the spike protein. ECSb1 bound to AS3 instead of the intended AS1, while ECSb2 bound to CA1 rather than AS2. These findings deviate from recent studies that have shown antibodies binding to the conserved fusion peptide region adjacent to the AS2 epitope (Premkumar et al., 2020; Dacon et al., 2022). The steric constraints imposed by the high density of spike proteins on the viral surface make this region difficult to access for antibodies (Low et al., 2022; Zhu et al., 2023). The unexpected binding patterns observed in ECSb1 and ECSb2 suggest that these Nbs may have evolved to overcome steric hindrances or target alternative regions on the spike protein.

The detailed features of the designed Nbs are presented in the supplemental material, providing a better comprehension of the biophysical principles underlying their binding efficiency. From the pKa−mapped ionizable residues essential for contact specificity (Supplementary Figure S1) to the structural predictions and validations that give our designed structures confidence (Supplementary Figures S2–S4). Supplementary Figures S1–S9 offer a visual compendium that captures everything. Interestingly, the surface patch analyses (Supplementary Figures S4–S6) point to possible aggregation sites, which is an important factor to consider during antibody design. This study demonstrates that our integration strategy simultaneously enhances binding affinity via EC optimization and improves aggregation propensity and thermostability, which should be considered when designing such candidates (Chiu et al., 2019; Muyldermans, 2021).

Electrostatic attraction plays a key factor in antibody-antigen interactions, enhancing binding specificity. Analysis of RBD epitopes based on their ionizable residues identified CA1 as the most favorable epitope for EC. Surface charge analysis of both Nb-RBD complexes and their binding interfaces reinforced the electrostatic differences between epitopes, with CA1 and AS3 being predominantly positively charged, while AS1 was more negative charged. These differences guided epitope-focused Nb design and mutation strategies to improve binding affinity. Engineered CDRs of ECSbs, paired with optimized FRs, enhanced EC, leading to higher binding specificity and stability (Tu et al., 2016). ECSbs demonstrated a better balance between net charge and electrostatic interactions, resulting in improved binding affinity to the RBD. This balance was particularly evident in ECSb2 and ECSb5, which exhibited high native structure scores and outperformed SR6c3 in terms of EC and CDR nativeness. The same idea could be implemented in bispecific antibody design, potentially resulting in more effective and scalable treatments (Fennell et al., 2013).

Rosetta docking, scoring, MD simulations, and MM/GBSA calculations revealed distinct mechanistic advantages and trade-offs for engineered nanobodies (ECSbs) compared to SR6c3. ECSb4-RBD, despite an initially favorable interface energy, exhibited suboptimal ddG and poor EC due to excessive hydrophobic clustering (39 instances) and high solvation penalties. This instability was further corroborated by MD simulations, which revealed moderate RMSF and RMSD values, indicative of limited polar specificity and ongoing interfacial fluctuations (Mangat et al., 2022). Conversely, ECSb3-RBD, with the highest Sc, demonstrated structural precision driven by extensive hydrogen bonding (17 H-bonds) and π-π stacking. While its MM/GBSA ΔG was high, suboptimal Coulombic contributions limited EC. Despite a well-packed interface, moderate EC and a positive ΔGCoulomb suggest unresolved charge misalignments limit binding stability, consistent with previous findings (Zhou and Pang, 2018). ECSb2-RBD showed intermediate performance with a broad SASA, indicating a solvent-exposed interface that enhanced conformational adaptability but reduced hydrophobic stabilization, as seen in its moderate EC and Sc. A transient MM/GBSA outlier, likely from interfacial disruptions, highlighted the flexibility-stability trade-off (Childers and Daggett, 2017). Excluding this outlier revealed a more favorable ΔG, indicating hydrophobic and lipophilic interactions primarily drove binding. ECSb5-RBD’s balance of rigidity (low RMSD/RMSF) and electrostatic optimization (high EC, salt bridges, π-cation interactions) yielded a competitive MM/GBSA ΔG and impressive ddG/dSASA, highlighting the efficacy of targeted electrostatic design for enhancing binding efficiency, even with a moderate interface (Zhang et al., 2011). ECSb1-RBD, relying on a deeply buried hydrophobic interface, exhibited the most favorable ddG but high solvation penalties and entropic costs. SR6c3, in contrast, showed limitations in electrostatic optimization, with moderate ΔG, RMSD, and RMSF, underscoring the advantage of rational design (Derat and Kamerlin, 2022). These analyses demonstrate that ECSbs achieved improved binding affinity through tailored interfacial optimization, highlighting the importance of balancing shape complementarity, electrostatic interactions, and hydrophobic packing for potent nanobody design.

Our approach, which is based on in-depth structural and biophysical investigations, allowed us to precisely tailor Nbs to match the electrostatic environment of the RBD. Specifically, the ECSb1-RBD complex demonstrated a consistent interaction profile across simulated parameters, suggesting a viable treatment target. The dynamic nature of the interactions between ECSb5 and RBD, along with the varying surface area exposure, highlights the ability of modified Nbs to adjust to changing viral structures (Mei et al., 2022; Shi et al., 2022). The higher stability and affinity seen in ECSb1-RBD and ECSb5-RBD complexes are fundamentally derived from the harmonious interplay between surface and electrostatic compatibility at the binding interfaces, as illustrated by the docking poses and complementarity plots (Supplementary Figure S8). This additional data not only supports the conclusions made in the main text but also offers a concrete tool for research into Nb design in the future (Price et al., 2018).

In silico techniques have revolutionized antibody development, making significant contributions across various stages of the process. These approaches have been instrumental in several key areas: identifying antigen epitopes (Chen et al., 2019; Amitai, 2021; Mason et al., 2021; Tubiana et al., 2022; Xu et al., 2022), exploring the vast sequence space of CDRs (Lim et al., 2022; Pooja Mahajan et al., 2022; Prihoda et al., 2022), optimizing CDR sequences, predicting antibody structures (Ruffolo et al., 2023), forecasting binding modes (Schneider et al., 2022), and estimating the binding affinity between antibodies and antigens (Myung et al., 2022). These computational approaches have been particularly valuable in the context of SARS-CoV-2, accelerating the discovery and optimization of therapeutic antibodies (Taft et al., 2021; 2022; Shan et al., 2022; Shaver et al., 2022; Zhang et al., 2022). These advancements are enhancing our ability to design antibodies with improved binding affinity, specificity, and other desirable properties. Our work significantly contributes to the growing body of knowledge in Nb development by providing a comprehensive dataset that correlates residue substitutions in Nb CDRs and their binding affinity to the SARS-CoV-2 spike protein. This dataset not only adds to our understanding of Nb-RBD interactions, but also to support the future in silico driven Nbs engineering efforts. Such engineered Nbs have potential beyond binding affinity. They can be used as scaffolds for targeted drug delivery, conjugated with therapeutic agents or imaging probes (Sakthikumar et al., 2022). Studying Nb-viral protein interactions helps understand resistance mechanisms, informing strategies against emerging variants (Nasution et al., 2018). This multifaceted approach enriches our comprehension of Nb functionality and paves the way for innovative solutions in combating infectious diseases. By leveraging computational technologies, our study can advance the development of potent and targeted Nb−based therapeutics to combat SARS-CoV-2 and other pathogens.

5 Conclusion

Conclusively, our innovative and integrated computational approach to Nb engineering has yielded promising results, with ECSbs demonstrating superior stability, binding affinity, and EC compared to the referenced SR6c3. Notably, ECSbs exhibited significantly higher thermostabilities (100.4–148.3 kcal·mol⁻1) compared to SR6c3 (62.6 kcal·mol⁻1). In terms of binding free energy (ΔG), ECSb4-RBD and ECSb3-RBD achieved superior values of −182.58 and −119.07 kcal·mol⁻1, respectively, versus −105.50 kcal·mol⁻1 for SR6c3-RBD. Additionally, ECSb4-RBD and ECSb3-RBD showed enhanced electrostatic complementarity (EC) values of 0.305 and 0.390, respectively, compared to 0.233 for SR6c3-RBD. These key findings make a strong case for the incorporation of electrostatic principles into Nb design, indicating a significant improvement in our ability to combat viral infections through creative biotechnological developments. further experimental validation and optimization of the ECSb design process are a prerequisite.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

ZI: Formal Analysis, Funding acquisition, Software, Supervision, Writing–original draft, Writing–review and editing. MA: Formal Analysis, Investigation, Methodology, Software, Writing–original draft. UK: Formal Analysis, Writing–original draft. NS: Data curation, Validation, Writing–original draft. IA: Conceptualization, Project administration, Supervision, Validation, Writing–original draft, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The publication of this work was supported by the Deanship of Scientific Research (DSR), Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia KFU250572.

Acknowledgments

Authors acknowledge all the support staff in the lab for their assistance during this work.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2025.1512788/full#supplementary-material

References

Aebersold, R., Agar, J. N., Amster, I. J., Baker, M. S., Bertozzi, C. R., Boja, E. S., et al. (2018). How many human proteoforms are there? Nat. Chem. Biol. 14, 206–214. doi:10.1038/nchembio.2576

PubMed Abstract | CrossRef Full Text | Google Scholar

Alford, R. F., Leaver-Fay, A., Jeliazkov, J. R., O’Meara, M. J., DiMaio, F. P., Park, H., et al. (2017). The Rosetta all-atom energy function for macromolecular modeling and design. J. Chem. Theory Comput. 13, 3031–3048. doi:10.1021/acs.jctc.7b00125

PubMed Abstract | CrossRef Full Text | Google Scholar

Amitai, A. (2021). Viral surface geometry shapes influenza and coronavirus spike evolution through antibody pressure. PLoS Comput. Biol. 17, e1009664. doi:10.1371/journal.pcbi.1009664

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakan, A., Dutta, A., Mao, W., Liu, Y., Chennubhotla, C., Lezon, T. R., et al. (2014). Evol and ProDy for bridging protein sequence evolution and structural dynamics. Bioinformatics 30, 2681–2683. doi:10.1093/bioinformatics/btu336

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakan, A., Meireles, L. M., and Bahar, I. (2011). ProDy: protein dynamics inferred from theory and experiments. Bioinformatics 27, 1575–1577. doi:10.1093/bioinformatics/btr168

PubMed Abstract | CrossRef Full Text | Google Scholar

Baker, E. N., and Hubbard, R. E. (1984). Hydrogen bonding in globular proteins. Prog. Biophys. Mol. Biol. 44, 97–179. doi:10.1016/0079-6107(84)90007-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Banks, J. L., Beard, H. S., Cao, Y., Cho, A. E., Damm, W., Farid, R., et al. (2005). Integrated modeling program, applied chemical theory (IMPACT). J. Comput. Chem. 26, 1752–1780. doi:10.1002/jcc.20292

PubMed Abstract | CrossRef Full Text | Google Scholar

Bas, D. C., Rogers, D. M., and Jensen, J. H. (2008). Very fast prediction and rationalization of pKa values for protein-ligand complexes. Proteins Struct. Funct. Genet. 73, 765–783. doi:10.1002/prot.22102

PubMed Abstract | CrossRef Full Text | Google Scholar

Basu, S., and Biswas, P. (2018). Salt-bridge dynamics in intrinsically disordered proteins: a trade-off between electrostatic interactions and structural flexibility. Biochim. Biophys. Acta Proteins Proteom 1866, 624–641. doi:10.1016/j.bbapap.2018.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Beard, H., Cholleti, A., Pearlman, D., Sherman, W., and Loving, K. A. (2013). Applying physics-based scoring to calculate free energies of binding for single amino acid mutations in protein-protein complexes. PLoS One 8, e82849. doi:10.1371/journal.pone.0082849

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamin Stranges, P., and Kuhlman, B. (2013). A comparison of successful and failed protein interface designs highlights the challenges of designing buried hydrogen bonds. Protein Sci. 22, 74–82. doi:10.1002/pro.2187

PubMed Abstract | CrossRef Full Text | Google Scholar

Biskupek, I., and Gieldon, A. (2024). Two-stage recognition mechanism of the SARS-CoV-2 receptor-binding domain to angiotensin-converting enzyme-2 (ACE2). Int. J. Mol. Sci. 25, 679. doi:10.3390/ijms25010679

PubMed Abstract | CrossRef Full Text | Google Scholar

Blom, N., Sicheritz-Pontén, T., Gupta, R., Gammeltoft, S., and Brunak, S. (2004). Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics 4, 1633–1649. doi:10.1002/pmic.200300771

PubMed Abstract | CrossRef Full Text | Google Scholar

Boehr, D. D., Nussinov, R., and Wright, P. E. (2009). The role of dynamic conformational ensembles in biomolecular recognition. Nat. Chem. Biol. 5, 789–796. doi:10.1038/nchembio.232

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowers, K. J., Chow, E., Xu, H., Dror, R. O., Eastwood, M. P., Gregersen, B. A., et al. (2006). “Scalable algorithms for molecular dynamics simulations on commodity clusters,” in Proceedings of the 2006 ACM/IEEE conference on supercomputing. doi:10.1145/1188455.1188544

CrossRef Full Text | Google Scholar

Braden, B. C., and Poljak, R. J. (1995). Structural features of the reactions between antibodies and protein antigens. FASEB J. 9, 9–16. doi:10.1096/fasebj.9.1.7821765

PubMed Abstract | CrossRef Full Text | Google Scholar

Bullen, G., Galson, J. D., Hall, G., Villar, P., Moreels, L., Ledsgaard, L., et al. (2021). Cross-reactive SARS-CoV-2 neutralizing antibodies from deep mining of early patient responses. Front. Immunol. 12, 678570. doi:10.3389/fimmu.2021.678570

PubMed Abstract | CrossRef Full Text | Google Scholar

Cameroni, E., Bowen, J. E., Rosen, L. E., Saliba, C., Zepeda, S. K., Culap, K., et al. (2021). Broadly neutralizing antibodies overcome SARS-CoV-2 Omicron antigenic shift. Nature 602, 664–670. doi:10.1038/s41586-021-04386-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, L., Goreshnik, I., Coventry, B., Case, J. B., Miller, L., Kozodoy, L., et al. (2020). De novo design of picomolar SARS-CoV-2 miniprotein inhibitors. Science 1979, 426–431. doi:10.1126/science.abd9909

PubMed Abstract | CrossRef Full Text | Google Scholar

Carballo-Amador, M. A., McKenzie, E. A., Dickson, A. J., and Warwicker, J. (2019). Surface patches on recombinant erythropoietin predict protein solubility: engineering proteins to minimise aggregation. BMC Biotechnol. 19, 26. doi:10.1186/s12896-019-0520-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Chakrabarti, P., and Janin, J. (2002). Dissecting protein-protein recognition sites. Proteins Struct. Funct. Genet. 47, 334–343. doi:10.1002/prot.10085

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaudhury, S., and Gray, J. J. (2008). Conformer selection and induced fit in flexible backbone protein-protein docking using computational and NMR ensembles. J. Mol. Biol. 381, 1068–1087. doi:10.1016/j.jmb.2008.05.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaudhury, S., Lyskov, S., and Gray, J. J. (2010). PyRosetta: a script-based interface for implementing molecular modeling algorithms using Rosetta. Bioinformatics 26, 689–691. doi:10.1093/bioinformatics/btq007

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheloha, R. W., Harmand, T. J., Wijne, C., Schwartz, T. U., and Ploegh, H. L. (2020). Exploring cellular biochemistry with nanobodies. J. Biol. Chem. 295, 15307–15327. doi:10.1074/jbc.REV120.012960

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, B., Khodadoust, M. S., Olsson, N., Wagar, L. E., Fast, E., Liu, C. L., et al. (2019). Predicting HLA class II antigen presentation through integrated deep learning. Nat. Biotechnol. 37, 1332–1343. doi:10.1038/s41587-019-0280-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, P., Nirula, A., Heller, B., Gottlieb, R. L., Boscia, J., Morris, J., et al. (2021a). SARS-CoV-2 neutralizing antibody LY-CoV555 in outpatients with covid-19. N. Engl. J. Med. 384, 229–237. doi:10.1056/nejmoa2029849

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Gentili, M., Hacohen, N., and Regev, A. (2021b). A cell-free nanobody engineering platform rapidly generates SARS-CoV-2 neutralizing nanobodies. Nat. Commun. 12, 5506. doi:10.1038/s41467-021-25777-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Childers, M. C., and Daggett, V. (2017). Insights from molecular dynamics simulations for computational protein design. Mol. Syst. Des. Eng. 2, 9–33. doi:10.1039/c6me00083e

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiu, M. L., Goulet, D. R., Teplyakov, A., and Gilliland, G. L. (2019). Antibody structure and function: the basis for engineering therapeutics. Antibodies 8, 55. doi:10.3390/antib8040055

PubMed Abstract | CrossRef Full Text | Google Scholar

Chothia, C., and Janin, J. (1975). Principles of protein-protein recognition. Nature 256, 705–708. doi:10.1038/256705a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chuang, G. Y., Kozakov, D., Brenke, R., Comeau, S. R., and Vajda, S. (2008). DARS (Decoys as the Reference State) potentials for protein-protein docking. Biophys. J. 95, 4217–4227. doi:10.1529/biophysj.108.135814

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, P. (2002). The origins of protein phosphorylation. Nat. Cell. Biol. 4, E127–E130. doi:10.1038/ncb0502-e127

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, T., Halfon, M., and Schneidman-Duhovny, D. (2022). NanoNet: rapid and accurate end-to-end nanobody modeling by deep learning. Front. Immunol. 13, 958584. doi:10.3389/fimmu.2022.958584

PubMed Abstract | CrossRef Full Text | Google Scholar

Coibion, O., Gorodnichenko, Y., and Weber, M. (2021). The cost of the COVID-19 crisis: lockdowns, macroeconomic expectations, and consumer spending. SSRN Electron. J. doi:10.2139/ssrn.3603849

CrossRef Full Text | Google Scholar

Contreras, M., Vicente, J., Cerón, J. J., Martinez Subiela, S., Urra, J. M., Rodríguez-del-Río, F. J., et al. (2023). Antibody isotype epitope mapping of SARS-CoV-2 Spike RBD protein: targets for COVID-19 symptomatology and disease control. Eur. J. Immunol. 53, e2250206. doi:10.1002/eji.202250206

PubMed Abstract | CrossRef Full Text | Google Scholar

Dacon, C., Tucker, C., Peng, L., Lee, C. C. D., Lin, T. H., Yuan, M., et al. (2022). Broadly neutralizing antibodies target the coronavirus fusion peptide. Science 1979, 728–735. doi:10.1126/science.abq3773

PubMed Abstract | CrossRef Full Text | Google Scholar

de Las Rivas, J., and Fontanillo, C. (2010). Protein-protein interactions essentials: key concepts to building and analyzing interactome networks. PLoS Comput. Biol. 6, e1000807. doi:10.1371/journal.pcbi.1000807

PubMed Abstract | CrossRef Full Text | Google Scholar

Dellisanti, C. D., Yao, Y., Stroud, J. C., Wang, Z. Z., and Chen, L. (2007). Crystal structure of the extracellular domain of nAChR alpha1 bound to alpha-bungarotoxin at 1.94 A resolution. Nat. Neurosci. 10, 953–962. doi:10.1038/nn1942

PubMed Abstract | CrossRef Full Text | Google Scholar

Derat, E., and Kamerlin, S. C. L. (2022). Computational advances in protein engineering and enzyme design. J. Phys. Chem. B 126, 2449–2451. doi:10.1021/acs.jpcb.2c01198

PubMed Abstract | CrossRef Full Text | Google Scholar

Desta, I. T., Porter, K. A., Xia, B., Kozakov, D., and Vajda, S. (2020). Performance and its limits in rigid body protein-protein docking. Structure 28, 1071–1081. doi:10.1016/j.str.2020.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Dolinsky, T. J., Czodrowski, P., Li, H., Nielsen, J. E., Jensen, J. H., Klebe, G., et al. (2007). PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations. Nucleic Acids Res. 35, W522–W525. doi:10.1093/nar/gkm276

PubMed Abstract | CrossRef Full Text | Google Scholar

Dolinsky, T. J., Nielsen, J. E., McCammon, J. A., and Baker, N. A. (2004). PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. 32, W665–W667. doi:10.1093/nar/gkh381

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, H., Hu, X., Duan, H., Yu, L., Qu, F., Huang, Q., et al. (2019a). Principles of inter-amino-acid recognition revealed by binding energies between homogeneous oligopeptides. ACS Cent. Sci. 5, 97–108. doi:10.1021/acscentsci.8b00723

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Q., Damschroder, M., Pabst, T. M., Hunter, A. K., Wang, W. K., and Luo, H. (2019b). Process optimization and protein engineering mitigated manufacturing challenges of a monoclonal antibody with liquid-liquid phase separation issue by disrupting inter-molecule electrostatic interactions. MAbs 11, 789–802. doi:10.1080/19420862.2019.1599634

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Z., Su, H., Wang, W., Ye, L., Wei, H., Peng, Z., et al. (2021). The trRosetta server for fast and accurate protein structure prediction. Nat. Protoc. 16, 5634–5651. doi:10.1038/s41596-021-00628-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Egede, L. E., and Walker, R. J. (2020). Structural racism, social risk factors, and covid-19 — a dangerous convergence for black Americans. N. Engl. J. Med. 383, e77. doi:10.1056/nejmp2023616

PubMed Abstract | CrossRef Full Text | Google Scholar

Entzminger, K. C., Hyun, J. M., Pantazes, R. J., Patterson-Orazem, A. C., Qerqez, A. N., Frye, Z. P., et al. (2017). De novo design of antibody complementarity determining regions binding a FLAG tetra-peptide. Sci. Rep. 7, 10295. doi:10.1038/s41598-017-10737-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Esparza, T. J., Martin, N. P., Anderson, G. P., Goldman, E. R., and Brody, D. L. (2020). High affinity nanobodies block SARS-CoV-2 spike receptor binding domain interaction with human angiotensin converting enzyme. Sci. Rep. 10, 22370. doi:10.1038/s41598-020-79036-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Farnsworth, C. C., Wolda, S. L., Gelb, M. H., and Glomset, J. A. (1989). Human lamin B contains a farnesylated cysteine residue. J. Biol. Chem. 264, 20422–20429. doi:10.1016/s0021-9258(19)47079-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Fennell, B. J., McDonnell, B., Tam, A. S. P., Chang, L., Steven, J., Broadbent, I. D., et al. (2013). CDR-restricted engineering of native human scFvs creates highly stable and soluble bifunctional antibodies for subcutaneous delivery. MAbs 5, 882–895. doi:10.4161/mabs.26201

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernández-Recio, J., and Sternberg, M. J. E. (2010). The 4th meeting on the critical assessment of predicted interaction (CAPRI) held at the mare nostrum, barcelona. Proteins Struct. Funct. Bioinforma. 78, 3065–3066. doi:10.1002/prot.22801

CrossRef Full Text | Google Scholar

Ferraz, M., Adan, W., Lima, T., Santos, A., Paula, S., Dhalia, R., et al. (2024). Design of nanobody targeting SARS-CoV-2 spike glycoprotein using CDR-grafting assisted by molecular simulation and machine learning. doi:10.1101/2024.09.30.615772

CrossRef Full Text | Google Scholar

Flaxman, A. D., Issema, R., Barnabas, R. V., and Ross, J. M. (2022). Estimated health outcomes and costs of COVID-19 prophylaxis with monoclonal antibodies among unvaccinated household contacts in the US. JAMA Netw. Open 5, E228632. doi:10.1001/jamanetworkopen.2022.8632

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukunaga, A., and Tsumoto, K. (2013). Improving the affinity of an antibody for its antigen via long-range electrostatic interactions. Protein Eng. Des. Sel. 26, 773–780. doi:10.1093/protein/gzt053

PubMed Abstract | CrossRef Full Text | Google Scholar

Gąciarz, A., and Ruddock, L. W. (2017). Complementarity determining regions and frameworks contribute to the disulfide bond independent folding of intrinsically stable scFv. PLoS One 12, e0189964. doi:10.1371/journal.pone.0189964

PubMed Abstract | CrossRef Full Text | Google Scholar

Ganesan, A., Siekierska, A., Beerten, J., Brams, M., Van Durme, J., De Baets, G., et al. (2016). Structural hot spots for the solubility of globular proteins. Nat. Commun. 7, 10816. doi:10.1038/ncomms10816

PubMed Abstract | CrossRef Full Text | Google Scholar

Gavin, A. C., Bösche, M., Krause, R., Grandi, P., Marzioch, M., Bauer, A., et al. (2002). Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature 415, 141–147. doi:10.1038/415141a

PubMed Abstract | CrossRef Full Text | Google Scholar

Genheden, S., and Ryde, U. (2015). The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 10, 449–461. doi:10.1517/17460441.2015.1032936

PubMed Abstract | CrossRef Full Text | Google Scholar

Golan, M. S., Jernegan, L. H., and Linkov, I. (2020). Trends and applications of resilience analytics in supply chain modeling: systematic literature review in the context of the COVID-19 pandemic. Environ. Syst. Decis. 40, 222–243. doi:10.1007/s10669-020-09777-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Gopal, R., Fitzpatrick, E., Pentakota, N., Jayaraman, A., Tharakaraman, K., and Capila, I. (2022). Optimizing antibody affinity and developability using a framework–CDR shuffling approach—application to an anti-SARS-CoV-2 antibody. Viruses 14, 2694. doi:10.3390/v14122694

PubMed Abstract | CrossRef Full Text | Google Scholar

Graham, B. S. (2020). Rapid COVID-19 vaccine development. Science 1979, 945–946. doi:10.1126/science.abb8923

PubMed Abstract | CrossRef Full Text | Google Scholar

Gray, J. J., Moughon, S., Wang, C., Schueler-Furman, O., Kuhlman, B., Rohl, C. A., et al. (2003). Protein-protein docking with simultaneous optimization of rigid-body displacement and side-chain conformations. J. Mol. Biol. 331, 281–299. doi:10.1016/S0022-2836(03)00670-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamers-Casterman, C., Atarhouch, T., Muyldermans, S., Robinson, G., Hammers, C., Songa, E. B., et al. (1993). Naturally occurring antibodies devoid of light chains. Nature 363, 446–448. doi:10.1038/363446a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Hannula, L., Kuivanen, S., Lasham, J., Kant, R., Kareinen, L., Bogacheva, M., et al. (2024). Nanobody engineering for SARS-CoV-2 neutralization and detection. Microbiol. Spectr. 12, e0419922. doi:10.1128/spectrum.04199-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, J., Baum, A., Pascal, K. E., Russo, V., Giordano, S., Wloga, E., et al. (2020). Studies in humanized mice and convalescent humans yield a SARS-CoV-2 antibody cocktail. Science 369 (369), 1010–1014. doi:10.1126/science.abd0827

PubMed Abstract | CrossRef Full Text | Google Scholar

Hermeling, S., Crommelin, D. J. A., Schellekens, H., and Jiskoot, W. (2004). Structure-immunogenicity relationships of therapeutic proteins. Pharm. Res. 21, 897–903. doi:10.1023/B:PHAM.0000029275.41323.a6

PubMed Abstract | CrossRef Full Text | Google Scholar

Hermet, P., Delache, B., Herate, C., Wolf, E., Kivi, G., Juronen, E., et al. (2023). Broadly neutralizing humanized SARS-CoV-2 antibody binds to a conserved epitope on Spike and provides antiviral protection through inhalation-based delivery in nonhuman primates. PLoS Pathog. 19, e1011532. doi:10.1371/journal.ppat.1011532

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoover, W. G. (1985). Canonical dynamics: equilibrium phase-space distributions. Phys. Rev. A Coll. Park 31, 1695–1697. doi:10.1103/PhysRevA.31.1695

PubMed Abstract | CrossRef Full Text | Google Scholar

Hornak, V., Abel, R., Okur, A., Strockbine, B., Roitberg, A., and Simmerling, C. (2006). Comparison of multiple amber force fields and development of improved protein backbone parameters. Proteins Struct. Funct. Genet. 65, 712–725. doi:10.1002/prot.21123

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacobson, M. P., Friesner, R. A., Xiang, Z., and Honig, B. (2002). On the role of the crystal environment in determining protein side-chain conformations. J. Mol. Biol. 320, 597–608. doi:10.1016/S0022-2836(02)00470-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacobson, M. P., Pincus, D. L., Rapp, C. S., Day, T. J. F., Honig, B., Shaw, D. E., et al. (2004). A hierarchical approach to all-atom protein loop prediction. Proteins Struct. Funct. Genet. 55, 351–367. doi:10.1002/prot.10613

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeong, B. S., Cha, J. S., Hwang, I., Kim, U., Adolf-Bryfogle, J., Coventry, B., et al. (2022). Computational design of a neutralizing antibody with picomolar binding affinity for all concerning SARS-CoV-2 variants. MAbs 14, 2021601. doi:10.1080/19420862.2021.2021601

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, S., and Thornton, J. M. (1996). Principles of protein-protein interactions. Proc. Natl. Acad. Sci. U. S. A. 93, 13–20. doi:10.1073/pnas.93.1.13

PubMed Abstract | CrossRef Full Text | Google Scholar

Jurrus, E., Engel, D., Star, K., Monson, K., Brandi, J., Felberg, L. E., et al. (2018). Improvements to the APBS biomolecular solvation software suite. Protein Sci. 27, 112–128. doi:10.1002/pro.3280

PubMed Abstract | CrossRef Full Text | Google Scholar

Kastritis, P. L., and Bonvin, A. M. J. J. (2013). On the binding affinity of macromolecular interactions: daring to ask why proteins interact. J. R. Soc. Interface 10, 20120835. doi:10.1098/rsif.2012.0835

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, A., Mohammad, A., Haq, I., Nasar, M., Ahmad, W., Yousafi, Q., et al. (2021). Structural-dynamics and binding analysis of rbd from sars-cov-2 variants of concern (Voces) and grp78 receptor revealed basis for higher infectivity. Microorganisms 9, 2331. doi:10.3390/microorganisms9112331

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiyoshi, M., Caaveiro, J. M. M., Miura, E., Nagatoishi, S., Nakakido, M., Soga, S., et al. (2014). Affinity improvement of a therapeutic antibody by structure-based computational design: generation of electrostatic interactions in the transition state stabilizes the antibody-antigen complex. PLoS One 9, e87099. doi:10.1371/journal.pone.0087099

PubMed Abstract | CrossRef Full Text | Google Scholar

Kniffin, K. M., Narayanan, J., Anseel, F., Antonakis, J., Ashford, S. P., Bakker, A. B., et al. (2021). COVID-19 and the workplace: implications, issues, and insights for future research and action. Am. Psychol. 76, 63–77. doi:10.1037/amp0000716

PubMed Abstract | CrossRef Full Text | Google Scholar

Korenkov, M., Zehner, M., Cohen-Dvashi, H., Borenstein-Katz, A., Kottege, L., Janicki, H., et al. (2023). Dissecting the impact of somatic hypermutation on SARS-CoV-2 neutralization and viral escape. bioRxiv, 2023–2025. doi:10.1101/2023.05.09.539943

CrossRef Full Text | Google Scholar

Korn, A. P., and Burnett, R. M. (1991). Distribution and complementarity of hydropathy in multisubunit proteins. Proteins Struct. Funct. Bioinforma. 9, 37–55. doi:10.1002/prot.340090106

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozakov, D., Brenke, R., Comeau, S. R., and Vajda, S. (2006). PIPER: an FFT-based protein docking program with pairwise potentials. Proteins Struct. Funct. Genet. 65, 392–406. doi:10.1002/prot.21117

PubMed Abstract | CrossRef Full Text | Google Scholar

Kramer, R. M., Shende, V. R., Motl, N., Pace, C. N., and Scholtz, J. M. (2012). Toward a molecular understanding of protein solubility: increased negative surface charge correlates with increased solubility. Biophys. J. 102, 1907–1915. doi:10.1016/j.bpj.2012.01.060

PubMed Abstract | CrossRef Full Text | Google Scholar

Krammer, F. (2020). SARS-CoV-2 vaccines in development. Nature 586, 516–527. doi:10.1038/s41586-020-2798-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Krieger, E., and Vriend, G. (2014). YASARA View - molecular graphics for all devices - from smartphones to workstations. Bioinformatics 30, 2981–2982. doi:10.1093/bioinformatics/btu426

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Prakash, S., Gupta, K., Dongre, A., Balaram, P., and Balaram, H. (2016). Unexpected functional implication of a stable succinimide in the structural stability of Methanocaldococcus jannaschii glutaminase. Nat. Commun. 7, 12798. doi:10.1038/ncomms12798

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Tsai, C. J., and Nussinov, R. (2000). Factors enhancing protein thermostability. Protein Eng. 13, 179–191. doi:10.1093/protein/13.3.179

PubMed Abstract | CrossRef Full Text | Google Scholar

Kunz, P., Zinner, K., Mücke, N., Bartoschik, T., Muyldermans, S., and Hoheisel, J. D. (2018). The structural basis of nanobody unfolding reversibility and thermoresistance. Sci. Rep. 8, 7934. doi:10.1038/s41598-018-26338-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Lan, J., Ge, J., Yu, J., Shan, S., Zhou, H., Fan, S., et al. (2020). Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor. Nature 581, 215–220. doi:10.1038/s41586-020-2180-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Laroche, A., Orsini Delgado, M. L., Chalopin, B., Cuniasse, P., Dubois, S., Sierocki, R., et al. (2022). Deep mutational engineering of broadly-neutralizing nanobodies accommodating SARS-CoV-1 and 2 antigenic drift. MAbs 14, 2076775. doi:10.1080/19420862.2022.2076775

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawrence, M. S., Phillips, K. J., and Liu, D. R. (2007). Supercharging proteins can impart unusual resilience. J. Am. Chem. Soc. 129, 10110–10112. doi:10.1021/ja071641y

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. M., Hammarén, H. M., Savitski, M. M., and Baek, S. H. (2023). Control of protein stability by post-translational modifications. Nat. Commun. 14, 201. doi:10.1038/s41467-023-35795-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Letko, M., Marzi, A., and Munster, V. (2020). Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses. Nat. Microbiol. 5, 562–569. doi:10.1038/s41564-020-0688-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Robertson, A. D., and Jensen, J. H. (2005). Very fast empirical prediction and rationalization of protein pK a values. Proteins Struct. Funct. Genet. 61, 704–721. doi:10.1002/prot.20660

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Abel, R., Zhu, K., Cao, Y., Zhao, S., and Friesner, R. A. (2011). The VSGB 2.0 model: a next generation energy model for high resolution protein structure modeling. Proteins Struct. Funct. Bioinforma. 79, 2794–2812. doi:10.1002/prot.23106

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Xue, W., Zheng, Q., Song, S., Yang, C., Xiong, H., et al. (2021). Cross-neutralizing antibodies bind a SARS-CoV-2 cryptic site and resist circulating variants. Nat. Commun. 12, 5652. doi:10.1038/s41467-021-25997-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Lijnzaad, P., Berendsen, H. J. C., and Argos, P. (1996). Hydrophobic patches on the surfaces of protein structures. Proteins Struct. Funct. Genet. 25, 389–397. doi:10.1002/(SICI)1097-0134(199607)25:3<389::AID-PROT10>3.0.CO;2-E

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, Y. W., Adler, A. S., and Johnson, D. S. (2022). Predicting antibody binders and generating synthetic antibodies using deep learning. MAbs 14, 2069075. doi:10.1080/19420862.2022.2069075

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Z., Akin, H., Rao, R., Hie, B., Zhu, Z., Lu, W., et al. (2022). Language models of protein sequences at the scale of evolution enable accurate structure prediction. bioRxiv 07 (20), 500902. Available at: https://www.biorxiv.org/content/10.1101/2022.07.20.500902v1%0Ahttps://www.biorxiv.org/content/10.1101/2022.07.20.500902v1.

Google Scholar

Lin, Z., Akin, H., Rao, R., Hie, B., Zhu, Z., Lu, W., et al. (2023). Evolutionary-scale prediction of atomic-level protein structure with a language model. Science 379, 1123–1130.

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Iketani, S., Guo, Y., Chan, J. F.-W., Wang, M., Liu, L., et al. (2021). Striking antibody evasion manifested by the Omicron variant of SARS-CoV-2. Nature 602, 676–681. doi:10.1038/s41586-021-04388-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Longsompurana, P., Rungrotmongkol, T., Plongthongkum, N., Wangkanont, K., Wolschann, P., and Pooarporn, R. P. (2023). Computational design of novel nanobodies targeting the receptor binding domain of variants of concern of SARS-CoV-2. PLoS One 18, e0293263. doi:10.1371/journal.pone.0293263

PubMed Abstract | CrossRef Full Text | Google Scholar

Low, J. S., Jerak, J., Tortorici, M. A., McCallum, M., Pinto, D., Cassotta, A., et al. (2022). ACE2-binding exposes the SARS-CoV-2 fusion peptide to broadly neutralizing coronavirus antibodies. Science 1979, 735–742. doi:10.1126/science.abq2679

PubMed Abstract | CrossRef Full Text | Google Scholar

Madhavi Sastry, G., Adzhigirey, M., Day, T., Annabhimoju, R., and Sherman, W. (2013). Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments. J. Comput. Aided Mol. Des. 27, 221–234. doi:10.1007/s10822-013-9644-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Maiorov, V. N., and Crippen, G. M. (1994). Significance of root-mean-square deviation in comparing three-dimensional structures of globular proteins. J. Mol. Biol. 235, 625–634. doi:10.1006/jmbi.1994.1017

PubMed Abstract | CrossRef Full Text | Google Scholar

Mangat, H. K., Rani, M., Pathak, R. K., Yadav, I. S., Utreja, D., Chhuneja, P. K., et al. (2022). Virtual screening, molecular dynamics and binding energy-MM-PBSA studies of natural compounds to identify potential EcR inhibitors against Bemisia tabaci Gennadius. PLoS One 17, e0261545. doi:10.1371/journal.pone.0261545

PubMed Abstract | CrossRef Full Text | Google Scholar

March, J., and Smith, M. B. (2020). Advanced organic chemistry: reactions, mechanisms, and structure. Chichester: John Wiley and Sons, 2357.

Google Scholar

Martyna, G. J., Tobias, D. J., and Klein, M. L. (1994). Constant pressure molecular dynamics algorithms. J. Chem. Phys. 101, 4177–4189. doi:10.1063/1.467468

CrossRef Full Text | Google Scholar

Marze, N. A., Roy Burman, S. S., Sheffler, W., and Gray, J. J. (2018). Efficient flexible backbone protein-protein docking for challenging targets. Bioinformatics 34, 3461–3469. doi:10.1093/bioinformatics/bty355

PubMed Abstract | CrossRef Full Text | Google Scholar

Mason, D. M., Friedensohn, S., Weber, C. R., Jordi, C., Wagner, B., Meng, S. M., et al. (2021). Optimization of therapeutic antibodies by predicting antigen specificity from antibody sequence via deep learning. Nat. Biomed. Eng. 5, 600–612. doi:10.1038/s41551-021-00699-9

PubMed Abstract | CrossRef Full Text | Google Scholar

McCoy, A. J., Chandana Epa, V., and Colman, P. M. (1997). Electrostatic complementarity at protein/protein interfaces. J. Mol. Biol. 268, 570–584. doi:10.1006/jmbi.1997.0987

PubMed Abstract | CrossRef Full Text | Google Scholar

Mei, Y., Chen, Y., Sivaccumar, J. P., An, Z., Xia, N., and Luo, W. (2022). Research progress and applications of nanobody in human infectious diseases. Front. Pharmacol. 13, 963978. doi:10.3389/fphar.2022.963978

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, M., Wilson, P., and Schomburg, D. (1996). Hydrogen bonding and molecular surface shape complementarity as a basis for protein docking. J. Mol. Biol. 264, 199–210. doi:10.1006/jmbi.1996.0634

PubMed Abstract | CrossRef Full Text | Google Scholar

Miklos, A. E., Kluwe, C., Der, B. S., Pai, S., Sircar, A., Hughes, R. A., et al. (2012). Structure-based design of supercharged, highly thermoresistant antibodies. Chem. Biol. 19, 449–455. doi:10.1016/j.chembiol.2012.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Mirdita, M., Schütze, K., Moriwaki, Y., Heo, L., Ovchinnikov, S., and Steinegger, M. (2022). ColabFold: making protein folding accessible to all. Nat. Methods 19, 679–682. doi:10.1038/s41592-022-01488-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Morrison, C. (2019). Nanobody approval gives domain antibodies a boost. Nat. Rev. Drug Discov. 18, 485–487. doi:10.1038/d41573-019-00104-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Muyldermans, S. (2013). Nanobodies: natural single-domain antibodies. Annu. Rev. Biochem. 82, 775–797. doi:10.1146/annurev-biochem-063011-092449

PubMed Abstract | CrossRef Full Text | Google Scholar

Muyldermans, S. (2021). A guide to: generation and design of nanobodies. FEBS J. 288, 2084–2102. doi:10.1111/febs.15515

PubMed Abstract | CrossRef Full Text | Google Scholar

Muyldermans, S., Atarhouch, T., Saldanha, J., Barbosa, J. A. R. G., and Hamers, R. (1994). Sequence and structure of vh domain from naturally occurring camel heavy chain immunoglobulins lacking light chains. Protein Eng. Des. Sel. 7, 1129–1135. doi:10.1093/protein/7.9.1129

PubMed Abstract | CrossRef Full Text | Google Scholar

Myung, Y., Pires, D. E. V., and Ascher, D. B. (2022). CSM-AB: graph-based antibody–antigen binding affinity prediction and docking scoring function. Bioinformatics 38, 1141–1143. doi:10.1093/bioinformatics/btab762

PubMed Abstract | CrossRef Full Text | Google Scholar

Nasution, M. A. F., Toepak, E. P., Alkaff, A. H., and Tambunan, U. S. F. (2018). Flexible docking-based molecular dynamics simulation of natural product compounds and Ebola virus Nucleocapsid (EBOV NP): a computational approach to discover new drug for combating Ebola. BMC Bioinforma. 19, 419. doi:10.1186/s12859-018-2387-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, H., Lan, P. D., Nissley, D. A., O’Brien, E. P., and Li, M. S. (2021). Electrostatic interactions explain the higher binding affinity of the CR3022 antibody for SARS-CoV-2 than the 4A8 antibody. J. Phys. Chem. B 125, 7368–7379. doi:10.1021/acs.jpcb.1c03639

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicholls, A., Sharp, K. A., and Honig, B. (1991). Protein folding and association: insights from the interfacial and thermodynamic properties of hydrocarbons. Proteins Struct. Funct. Bioinforma. 11, 281–296. doi:10.1002/prot.340110407

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicola, M., Alsafi, Z., Sohrabi, C., Kerwan, A., Al-Jabir, A., Iosifidis, C., et al. (2020). The socio-economic implications of the coronavirus pandemic (COVID-19): a review. Int. J. Surg. 78, 185–193. doi:10.1016/j.ijsu.2020.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Nissley, D. A., Jiang, Y., Trovato, F., Sitarik, I., Narayan, K. B., To, P., et al. (2022). Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less functional. Nat. Commun. 13, 3081. doi:10.1038/s41467-022-30548-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Nivón, L. G., Moretti, R., and Baker, D. (2013). A pareto-optimal refinement method for protein design scaffolds. PLoS One 8, e59004. doi:10.1371/journal.pone.0059004

PubMed Abstract | CrossRef Full Text | Google Scholar

Nosé, S. (1984). A molecular dynamics method for simulations in the canonical ensemble. Mol. Phys. 52, 255–268. doi:10.1080/00268978400101201

CrossRef Full Text | Google Scholar

Ohtsubo, K., and Marth, J. D. (2006). Glycosylation in cellular mechanisms of health and disease. Cell. 126, 855–867. doi:10.1016/j.cell.2006.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Olsson, M. H. M., Søndergaard, C. R., Rostkowski, M., and Jensen, J. H. (2011). PROPKA3: consistent treatment of internal and surface residues in empirical p K a predictions. J. Chem. Theory Comput. 7, 525–537. doi:10.1021/ct100578z

PubMed Abstract | CrossRef Full Text | Google Scholar

Padlan, E. A. (1994). Anatomy of the antibody molecule. Mol. Immunol. 31, 169–217. doi:10.1016/0161-5890(94)90001-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, Y.-J., De Marco, A., Starr, T. N., Liu, Z., Pinto, D., Walls, A. C., et al. (2022). Antibody-mediated broad sarbecovirus neutralization through ACE2 molecular mimicry. Science 375, 449–454. doi:10.1126/science.abm8143

PubMed Abstract | CrossRef Full Text | Google Scholar

Patil, A., and Nakamura, H. (2007). The role of charged surface residues in the binding ability of small hubs in protein-protein interaction networks. Biophys. (Oxf) 3, 27–35. doi:10.2142/biophysics.3.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Peruzzi, J. A., Steinkühler, J., Vu, T. Q., Gunnels, T. F., Hu, V. T., Lu, P., et al. (2024). Hydrophobic mismatch drives self-organization of designer proteins into synthetic membranes. Nat. Commun. 15, 3162. doi:10.1038/s41467-024-47163-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinto, D., Park, Y. J., Beltramello, M., Walls, A. C., Tortorici, M. A., Bianchi, S., et al. (2020). Cross-neutralization of SARS-CoV-2 by a human monoclonal SARS-CoV antibody. Nature 583, 290–295. doi:10.1038/s41586-020-2349-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Pooja Mahajan, S., Ruffolo, J., Frick, R., and Gray, J. J. (2022). Towards deep learning models for target-specific antibody design. Biophys. J. 121, 528a. doi:10.1016/j.bpj.2021.11.2783

CrossRef Full Text | Google Scholar

Popping, S., Nichols, B. E., Appelman, B., Biemond, J. J., Vergouwe, M., Rosendaal, F. R., et al. (2023). Health outcomes and cost-effectiveness of monoclonal SARS-CoV-2 antibodies as pre-exposure prophylaxis. JAMA Netw. Open 6, E2321985. doi:10.1001/jamanetworkopen.2023.21985

PubMed Abstract | CrossRef Full Text | Google Scholar

Premkumar, L., Segovia-Chumbez, B., Jadi, R., Martinez, D. R., Raut, R., Markmann, A. J., et al. (2020). The receptor-binding domain of the viral spike protein is an immunodominant and highly specific target of antibodies in SARS-CoV-2 patients. Sci. Immunol. 5, eabc8413. doi:10.1126/SCIIMMUNOL.ABC8413

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, A., Schroter, S., Clarke, M., and McAneney, H. (2018). Role of supplementary material in biomedical journal articles: surveys of authors, reviewers and readers. BMJ Open 8, e021753. doi:10.1136/bmjopen-2018-021753

PubMed Abstract | CrossRef Full Text | Google Scholar

Prihoda, D., Maamary, J., Waight, A., Juan, V., Fayadat-Dilman, L., Svozil, D., et al. (2022). BioPhi: a platform for antibody design, humanization, and humanness evaluation based on natural antibody repertoires and deep learning. MAbs 14, 2020203. doi:10.1080/19420862.2021.2020203

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramachandran, G. N., Ramakrishnan, C., and Sasisekharan, V. (1963). Stereochemistry of polypeptide chain configurations. J. Mol. Biol. 7, 95–99. doi:10.1016/s0022-2836(63)80023-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramon, A., Ali, M., Atkinson, M., Saturnino, A., Didi, K., Visentin, C., et al. (2024). Assessing antibody and nanobody nativeness for hit selection and humanization with AbNatiV. Nat. Mach. Intell. 6, 74–91. doi:10.1038/s42256-023-00778-3

CrossRef Full Text | Google Scholar

Rangel, M. A., Bedwell, A., Costanzi, E., Taylor, R. J., Russo, R., Bernardes, G. J. L., et al. (2022). Fragment-based computational design of antibodies targeting structured epitopes. Sci. Adv. 8, eabp9540. doi:10.1126/sciadv.abp9540

PubMed Abstract | CrossRef Full Text | Google Scholar

Ranney, M. L., Griffeth, V., and Jha, A. K. (2020). Critical supply shortages — the need for ventilators and personal protective equipment during the covid-19 pandemic. N. Engl. J. Med. 382, e41. doi:10.1056/nejmp2006141

PubMed Abstract | CrossRef Full Text | Google Scholar

Richards, F. M. (1977). Areas, volumes, packing and protein structure. Annu. Rev. Biophys. Bioeng. 6, 151–176. doi:10.1146/annurev.bb.06.060177.001055

PubMed Abstract | CrossRef Full Text | Google Scholar

Robbiani, D. F., Gaebler, C., Muecksch, F., Lorenzi, J. C. C., Wang, Z., Cho, A., et al. (2020). Convergent antibody responses to SARS-CoV-2 in convalescent individuals. Nature 584, 437–442. doi:10.1038/s41586-020-2456-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Roos, K., Wu, C., Damm, W., Reboul, M., Stevenson, J. M., Lu, C., et al. (2019). OPLS3e: Extending Force Field Coverage for Drug-Like Small Molecules. J. Chem. Theory Comput. 15, 1863–1874. doi:10.1021/acs.jctc.8b01026

PubMed Abstract | CrossRef Full Text | Google Scholar

Rose, G. D., Geselowitz, A. R., Lesser, G. J., Lee, R. H., and Zehfus, M. H. (1985). Hydrophobicity of amino acid residues in globular proteins. Science 229 (229), 834–838. doi:10.1126/science.4023714

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruffolo, J. A., Chu, L. S., Mahajan, S. P., and Gray, J. J. (2023). Fast, accurate antibody structure prediction from deep learning on massive set of natural antibodies. Nat. Commun. 14, 2389. doi:10.1038/s41467-023-38063-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakthikumar, K., Isamura, B. K., and Krause, R. W. M. (2022). An integrated analysis of mechanistic insights into biomolecular interactions and molecular dynamics of bio-inspired Cu(II) and Zn(II) complexes towards DNA/BSA/SARS-CoV-2 3CLpro by molecular docking-based virtual screening and FRET detection. Biomolecules 12, 1883. doi:10.3390/biom12121883

PubMed Abstract | CrossRef Full Text | Google Scholar

Salam, N. K., Adzhigirey, M., Sherman, W., and Pearlman, D. A. (2014). Structure-based approach to the prediction of disulfide bonds in proteins. Protein Eng. Des. Sel. 27, 365–374. doi:10.1093/protein/gzu017

PubMed Abstract | CrossRef Full Text | Google Scholar

Sankar, K., Krystek, S. R., Carl, S. M., Day, T., and Maier, J. K. X. (2018). AggScore: prediction of aggregation-prone regions in proteins based on the distribution of surface patches. Proteins Struct. Funct. Bioinforma. 86, 1147–1156. doi:10.1002/prot.25594

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, C., Buchanan, A., Taddese, B., and Deane, C. M. (2022). DLAB: deep learning methods for structure-based virtual screening of antibodies. Bioinformatics 38, 377–383. doi:10.1093/bioinformatics/btab660

PubMed Abstract | CrossRef Full Text | Google Scholar

Schreiber, G., and Fersht, A. R. (1995). Energetics of protein-protein interactions: analysis ofthe Barnase-Barstar interface by single mutations and double mutant cycles. J. Mol. Biol. 248, 478–486. doi:10.1016/S0022-2836(95)80064-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Schrödinger (2023a). Schrödinger release 2023-3: BioLuminate. Schrödinger, LLC. Available at: https://www.schrodinger.com/platform/products/bioluminate/.

Google Scholar

Schrödinger (2023b). Schrödinger release 2023-3: Maestro. Schrödinger, LLC. Available at: https://www.schrodinger.com/platform/products/maestro/.

Google Scholar

Schrödinger (2023c). Schrödinger release 2023-3: PIPER. Schrödinger, LLC. Available at: https://www.schrodinger.com/platform/products/piper/.

Google Scholar

Schrödinger (2023d). Schrödinger release 2023-3: Prime. Schrödinger, LLC. Available at: https://www.schrodinger.com/platform/products/prime/.

Google Scholar

Schrödinger (2023e). Schrödinger release 2023-3: protein preparation wizard. Schrödinger, LLC. Available at: https://www.schrodinger.com/life-science/learn/white-papers/protein-preparation-workflow/.

Google Scholar

Schrödinger, L. L. C. (2024). The PyMOL molecular graphics system. Schrödinger, LLC. Available at: https://www.pymol.org/.

Google Scholar

Schymkowitz, J., Borg, J., Stricher, F., Nys, R., Rousseau, F., and Serrano, L. (2005). The FoldX web server: an online force field. Nucleic Acids Res. 33, W382–W388. doi:10.1093/nar/gki387

PubMed Abstract | CrossRef Full Text | Google Scholar

Selzer, T., Albeck, S., and Schreiber, G. (2000). Rational design of faster associating and tighter binding protein complexes. Nat. Struct. Biol. 7, 537–541. doi:10.1038/76744

PubMed Abstract | CrossRef Full Text | Google Scholar

Shan, S., Luo, S., Yang, Z., Hong, J., Su, Y., Ding, F., et al. (2022). Deep learning guided optimization of human antibody against SARS-CoV-2 variants with broad neutralization. Proc. Natl. Acad. Sci. U. S. A. 119, e2122954119. doi:10.1073/pnas.2122954119

PubMed Abstract | CrossRef Full Text | Google Scholar

Shang, J., Ye, G., Shi, K., Wan, Y., Luo, C., Aihara, H., et al. (2020). Structural basis of receptor recognition by SARS-CoV-2. Nature 581, 221–224. doi:10.1038/s41586-020-2179-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Shaver, J. M., Smith, J., and Amimeur, T. (2022). “Deep learning in therapeutic antibody development,” in Methods in molecular biology. doi:10.1007/978-1-0716-1787-8_19

CrossRef Full Text | Google Scholar

Sheinerman, F. B., Norel, R., and Honig, B. (2000). Electrostatic aspects of protein–protein interactions. Curr. Opin. Struct. Biol. 10, 153–159. doi:10.1016/S0959-440X(00)00065-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Z., Li, X., Wang, L., Sun, Z., Zhang, H., Chen, X., et al. (2022). Structural basis of nanobodies neutralizing SARS-CoV-2 variants. Structure 30, 707–720.e5. doi:10.1016/j.str.2022.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, V., Bhutkar, M., Choudhary, S., Nehul, S., Kumar, R., Singla, J., et al. (2024). Structure-guided mutations in CDRs for enhancing the affinity of neutralizing SARS-CoV-2 nanobody. Biochem. Biophys. Res. Commun. 734, 150746. doi:10.1016/j.bbrc.2024.150746

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinha, N., Mohan, S., Lipschultz, C. A., and Smith-Gill, S. J. (2002). Differences in electrostatic properties at antibody-antigen binding sites: implications for specificity and cross-reactivity. Biophys. J. 83, 2946–2968. doi:10.1016/S0006-3495(02)75302-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinha, N., and Smith-Gill, S. J. (2002). Electrostatics in protein binding and function. Curr. Protein Pept. Sci. 3, 601–614. doi:10.2174/1389203023380431

PubMed Abstract | CrossRef Full Text | Google Scholar

Starr, T. N., Czudnochowski, N., Liu, Z., Zatta, F., Park, Y. J., Addetia, A., et al. (2021). SARS-CoV-2 RBD antibodies that maximize breadth and resistance to escape. Nature 597, 97–102. doi:10.1038/s41586-021-03807-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Steeland, S., Vandenbroucke, R. E., and Libert, C. (2016). Nanobodies as therapeutics: big opportunities for small antibodies. Drug Discov. Today 21, 1076–1113. doi:10.1016/j.drudis.2016.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, H., Wang, W., Du, Z., Peng, Z., Gao, S. H., Cheng, M. M., et al. (2021). Improved protein structure prediction using a new multi-scale network and homologous templates. Adv. Sci. 8, e2102592. doi:10.1002/advs.202102592

PubMed Abstract | CrossRef Full Text | Google Scholar

Taft, J. M., Weber, C. R., Gao, B., Ehling, R. A., Han, J., Frei, L., et al. (2021). Predictive profiling of SARS-CoV-2 variants by deep mutational learning. SSRN Electron. J. doi:10.2139/ssrn.3982671

CrossRef Full Text | Google Scholar

Taft, J. M., Weber, C. R., Gao, B., Ehling, R. A., Han, J., Frei, L., et al. (2022). Deep mutational learning predicts ACE2 binding and antibody escape to combinatorial mutations in the SARS-CoV-2 receptor-binding domain. Cell. 185, 4008–4022.e14. doi:10.1016/j.cell.2022.08.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Tai, D. B. G., Shah, A., Doubeni, C. A., Sia, I. G., and Wieland, M. L. (2021). The disproportionate impact of COVID-19 on racial and ethnic minorities in the United States. Clin. Infect. Dis. 72, 703–706. doi:10.1093/cid/ciaa815

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, X., Li, C., Huang, A., Xia, S., Lu, S., Shi, Z., et al. (2020). Potent binding of 2019 novel coronavirus spike protein by a SARS coronavirus-specific human monoclonal antibody. Emerg. Microbes Infect. 9, 382–385. doi:10.1080/22221751.2020.1729069

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsien, R. Y. (1998). The green fluorescent protein. Annu. Rev. Biochem. 67, 509–544. doi:10.1146/annurev.biochem.67.1.509

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsutakawa, S. E., Thompson, M. J., Arvai, A. S., Neil, A. J., Shaw, S. J., Algasaier, S. I., et al. (2017). Phosphate steering by Flap Endonuclease 1 promotes 5′-flap specificity and incision to prevent genome instability. Nat. Commun. 8, 15855. doi:10.1038/ncomms15855

PubMed Abstract | CrossRef Full Text | Google Scholar

Tu, C., Terraube, V., Tam, A. S. P., Stochaj, W., Fennell, B. J., Lin, L., et al. (2016). A combination of structural and empirical analyses delineates the key contacts mediating stability and affinity increases in an optimized biotherapeutic single-chain fv (scFv). J. Biol. Chem. 291, 1267–1276. doi:10.1074/jbc.M115.688010

PubMed Abstract | CrossRef Full Text | Google Scholar

Tubiana, J., Schneidman-Duhovny, D., and Wolfson, H. J. (2022). ScanNet: an interpretable geometric deep learning model for structure-based protein binding site prediction. Nat. Methods 19, 730–739. doi:10.1038/s41592-022-01490-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Gils, J. H. M., Gogishvili, D., Van Eck, J., Bouwmeester, R., Van Dijk, E., and Abeln, S. (2022). How sticky are our proteins? Quantifying hydrophobicity of the human proteome. Bioinforma. Adv. 2, vbac002. doi:10.1093/bioadv/vbac002

PubMed Abstract | CrossRef Full Text | Google Scholar

Vijayakumar, M., Wong, K. Y., Schreiber, G., Fersht, A. R., Szabo, A., and Zhou, H. X. (1998). Electrostatic enhancement of diffusion-controlled protein-protein association: comparison of theory and experiment on barnase and barstar. J. Mol. Biol. 278, 1015–1024. doi:10.1006/jmbi.1998.1747

PubMed Abstract | CrossRef Full Text | Google Scholar

V’kovski, P., Kratzel, A., Steiner, S., Stalder, H., and Thiel, V. (2021). Coronavirus biology and replication: implications for SARS-CoV-2. Nat. Rev. Microbiol. 19, 155–170. doi:10.1038/s41579-020-00468-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogt, G., Woell, S., and Argos, P. (1997). Protein thermal stability, hydrogen bonds, and ion pairs. J. Mol. Biol. 269, 631–643. doi:10.1006/jmbi.1997.1042

PubMed Abstract | CrossRef Full Text | Google Scholar

Vu, K. B., Ghahroudi, M. A., Wyns, L., and Muyldermans, S. (1997). Comparison of llama VH sequences from conventional and heavy chain antibodies. Mol. Immunol. 34, 1121–1131. doi:10.1016/S0161-5890(97)00146-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wade, R. C., Gabdoulline, R. R., Lüdemann, S. K., and Lounnas, V. (1998). Electrostatic steering and ionic tethering in enzyme-ligand binding: insights from simulations. Proc. Natl. Acad. Sci. U. S. A. 95, 5942–5949. doi:10.1073/pnas.95.11.5942

PubMed Abstract | CrossRef Full Text | Google Scholar

Walls, A. C., Park, Y. J., Tortorici, M. A., Wall, A., McGuire, A. T., and Veesler, D. (2020). Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell. 181, 281–292. doi:10.1016/j.cell.2020.02.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Walsh, G. (2010). Biopharmaceutical benchmarks. Nat. Biotechnol. 28, 917–924. doi:10.1038/nbt0910-917

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Schueler-Furman, O., and Baker, D. (2005). Improved side-chain modeling for protein–protein docking. Protein Sci. 14, 1328–1339. doi:10.1110/ps.041222905

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Guo, Y., Iketani, S., Nair, M. S., Li, Z., Mohri, H., et al. (2022a). Antibody evasion by SARS-CoV-2 Omicron subvariants BA.2.12.1, BA.4 and BA.5. Nature 608, 603–608. doi:10.1038/s41586-022-05053-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Zhang, Y., Wu, L., Niu, S., Song, C., Zhang, Z., et al. (2020). Structural and functional basis of SARS-CoV-2 entry by using human ACE2. Cell. 181, 894–904. doi:10.1016/j.cell.2020.03.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Peng, Z., and Yang, J. (2022b). Single-sequence protein structure prediction using supervised transformer protein language models. Nat. Comput. Sci. 2, 804–814. doi:10.1038/s43588-022-00373-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, R., Ding, F., Wang, R., Shen, R., Zhang, X., Luo, S., et al. (2022). High-resolution de novo structure prediction from primary sequence. bioRxiv. Available at: https://www.biorxiv.org/content/10.1101/2022.07.21.500999v1 (Accessed November 1, 2024).

Google Scholar

Wu, Y., Li, C., Xia, S., Tian, X., Kong, Y., Wang, Z., et al. (2020a). Identification of human single-domain antibodies against SARS-CoV-2. Cell. Host Microbe 27, 891–898. doi:10.1016/j.chom.2020.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Wang, F., Shen, C., Peng, W., Li, D., Zhao, C., et al. (2020b). A noncompeting pair of human neutralizing antibodies block COVID-19 virus binding to its receptor ACE2. Science 368, 1274–1278. doi:10.1126/science.abc2241

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Z., and McGoogan, J. M. (2020). Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72314 cases from the Chinese center for disease control and prevention. JAMA - J. Am. Med. Assoc. 323, 1239–1242. doi:10.1001/jama.2020.2648

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, D., Lin, S. L., and Nussinov, R. (1997). Protein binding versus protein folding: the role of hydrophilic bridges in protein associations. J. Mol. Biol. 265, 68–84. doi:10.1006/jmbi.1996.0712

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Z., Davila, A., Wilamowski, J., Teraguchi, S., and Standley, D. M. (2022). Improved antibody-specific epitope prediction using AlphaFold and AbAdapt. ChemBioChem 23, e202200303. doi:10.1002/cbic.202200303

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Zhang, Z., Yang, F., Zhang, H., Wu, H., Zhu, F., et al. (2021). Computational design and modeling of nanobodies toward SARS-CoV-2 receptorbinding domain. Chem. Biol. Drug Des. 98, 1–18. doi:10.1111/cbdd.13847

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Duan, H., Liu, X., Zhang, X., Pan, S., Zhang, F., et al. (2023). Broad sarbecovirus neutralizing antibodies obtained by computational design and synthetic library screening. J. Virol. 97, e0061023. doi:10.1128/jvi.00610-23

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, K., Kuroda, D., Kiyoshi, M., Nakakido, M., Nagatoishi, S., Soga, S., et al. (2019). Exploring designability of electrostatic complementarity at an antigen-antibody interface directed by mutagenesis, biophysical analysis, and molecular dynamics simulations. Sci. Rep. 9, 4482. doi:10.1038/s41598-019-40461-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, M., Wu, N. C., Zhu, X., Lee, C. C. D., So, R. T. Y., Lv, H., et al. (2020). A highly conserved cryptic epitope in the receptor binding domains of SARS-CoV-2 and SARS-CoV. Science 368, 630–633. doi:10.1126/science.abb7269

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Ghosh, S., and Pal, R. (2022). Predicting binding affinities of emerging variants of SARS-CoV-2 using spike protein sequencing data: observations, caveats and recommendations. Brief. Bioinform 23, bbac128. doi:10.1093/bib/bbac128

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Krieger, J. M., Zhang, Y., Kaya, C., Kaynak, B., Mikulska-Ruminska, K., et al. (2021). ProDy 2.0: increased scale and scope after 10 years of protein dynamics modelling with Python. Bioinformatics 37, 3657–3659. doi:10.1093/bioinformatics/btab187

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Witham, S., and Alexov, E. (2011). On the role of electrostatics in protein-protein interactions. Phys. Biol. 8, 035001. doi:10.1088/1478-3975/8/3/035001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, H.-X., and Pang, X. (2018). Electrostatic interactions in protein structure, folding, binding, and condensation. Chem. Rev. 118, 1691–1741. doi:10.1021/acs.chemrev.7b00305

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, P., Yang, X. L., Wang, X. G., Hu, B., Zhang, L., Zhang, W., et al. (2020). A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579, 270–273. doi:10.1038/s41586-020-2012-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, K., Day, T., Warshaviak, D., Murrett, C., Friesner, R., and Pearlman, D. (2014). Antibody structure determination using a combination of homology modeling, energy-based refinement, and loop prediction. Proteins Struct. Funct. Bioinforma. 82, 1646–1655. doi:10.1002/prot.24551

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., An, K., Yan, J., Xu, P., and Bai, C. (2023). In silico optimization of SARS-CoV-2 spike specific nanobodies. Front. Biosci. - Landmark 28, 67. doi:10.31083/j.fbl2804067

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: electrostatic potential, nanobodies, ACE-2 receptor, SARS-CoV-2, spike protein, receptor binding domain, molecular docking simulation

Citation: Iqbal Z, Asim M, Khan UA, Sultan N and Ali I (2025) Computational electrostatic engineering of nanobodies for enhanced SARS−CoV−2 receptor binding domain recognition. Front. Mol. Biosci. 12:1512788. doi: 10.3389/fmolb.2025.1512788

Received: 17 October 2024; Accepted: 11 February 2025;
Published: 10 March 2025.

Edited by:

Elena Rosini, University of Insubria, Italy

Reviewed by:

Marco Orlando, University of Milano-Bicocca, Italy
Chang Liu, The University of Texas at Austin, United States

Copyright © 2025 Iqbal, Asim, Khan, Sultan and Ali. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zafar Iqbal, emFmYXJAa2Z1LmVkdS5zYQ==; Irfan Ali , aXJmYW4uY2FiYkB1YWYuZWR1LnBr

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Research integrity at Frontiers

Man ultramarathon runner in the mountains he trains at sunset

94% of researchers rate our articles as excellent or good

Learn more about the work of our research integrity team to safeguard the quality of each article we publish.


Find out more