Application of an integrated computational antibody engineering platform to design SARS-CoV-2 neutralizers.

As the COVID-19 pandemic continues to spread, hundreds of new initiatives including studies on existing medicines are running to fight the disease. To deliver a potentially immediate and lasting treatment to current and emerging SARS-CoV-2 variants, new collaborations and ways of sharing are required to create as many paths forward as possible. Here, we leverage our expertise in computational antibody engineering to rationally design/engineer three previously reported SARS-CoV neutralizing antibodies and share our proposal towards anti-SARS-CoV-2 biologics therapeutics. SARS-CoV neutralizing antibodies, m396, 80R and CR-3022 were chosen as templates due to their diversified epitopes and confirmed neutralization potency against SARS-CoV (but not SARS-CoV-2 except for CR3022). Structures of variable fragment (Fv) in complex with receptor binding domain (RBD) from SARS-CoV or SARS-CoV-2 were subjected to our established in silico antibody engineering platform to improve their binding affinity to SARS-CoV-2 and developability profiles. The selected top mutations were ensembled into a focused library for each antibody for further screening. In addition, we convert the selected binders with different epitopes into the trispecific format, aiming to increase potency and to prevent mutational escape. Lastly, to avoid antibody-induced virus activation or enhancement, we suggest application of NNAS and DQ mutations to the Fc region to eliminate effector functions and extend half-life.

Improving the analytical toolbox to investigate copurifying host cell proteins presence: N-(4)-(ß-acetylglucosaminyl)- l-asparaginase case study.

Levels of host cell proteins (HCPs) in purification intermediates and drug substances (DS) of monoclonal antibodies (mAbs) must be carefully monitored for the production of safe and efficacious biotherapeutics. During the development of mAb1, an immunoglobulin G1 product, unexpected results generated with HCP Enzyme-Linked Immunosorbent Assay (ELISA) kit triggered an investigation which led to the identification of a copurifying HCP called N-(4)-(ß-acetylglucosaminyl)-l-asparaginase (AGA, EC3.5.1.26) by liquid chromatography-tandem mass spectrometry (LC-MS/MS). The risk assessment performed indicated a low immunogenicity risk for the copurifying HCP and an ad hoc stability study demonstrated no mAb glycan cleavage and thus no impact on product quality. Fractionation studies performed on polishing steps revealed that AGA was coeluted with the mAb. Very interestingly, the native digestion protocol implemented to go deeper in the MS-HCP profiling was found to be incompatible with correct AGA detection in last purification intermediate and DS, further suggesting a hitchhiking behavior of AGA. In silico surface characterization of AGA also supports this hypothesis. Finally, the combined support of HCP ELISA results and MS allowed process optimization and removal of this copurifying HCP.

The impact of proline isomerization on antigen binding and the analytical profile of a trispecific anti-HIV antibody.

Proline cis-trans conformational isomerization is a mechanism that affects different types of protein functions and behaviors. Using analytical characterization, structural analysis, and molecular dynamics simulations, we studied the causes of an aberrant two-peak size-exclusion chromatography profile observed for a trispecific anti-HIV antibody. We found that proline isomerization in the tyrosine-proline-proline (YPP) motif in the heavy chain complementarity-determining region (CDR)3 domain of one of the antibody arms (10e8v4) was a component of this profile. The pH effect on the conformational equilibrium that led to these two populations was presumably caused by a histidine residue (H147) in the light chain that is in direct contact with the YPP motif. Finally, we demonstrated that, due to chemical equilibrium between the cis and trans proline conformers, the antigen-binding potency of the trispecific anti-HIV antibody was not significantly affected in spite of a potential structural clash of 10e8v4 YPtransPtrans conformers with the membrane-proximal ectodomain region epitope in the GP41 antigen. Altogether, these results reveal at mechanistic and molecular levels the effect of proline isomerization in the CDR on the antibody binding and analytical profiles, and support further development of the trispecific anti-HIV antibody.

Intracellular delivery of therapeutic antibodies into specific cells using antibody-peptide fusions.

Because of their favorable properties as macromolecular drugs, antibodies are a very successful therapeutic modality for interfering with disease-relevant targets in the extracellular space or at the cell membrane. However, a large number of diseases involve cytosolic targets and designing antibodies able to efficiently reach intracellular compartments would expand the antibody-tractable conditions. Here, we genetically fused cell penetrating peptides (CPPs) at various positions to an antibody targeting cancer cells, evaluated the developability features of the resulting antibody-peptide fusions and the ability of selected constructs to reach the cytosol. We first determined positions in the IgG structure that were permissive to CPP incorporation without destabilizing the antibody. Fusing CPPs to the C-terminus of the light chain and either before or after the hinge had the least effect on antibody developability features. These constructs were further evaluated for cell penetration efficiency. Two out of five tested CPPs significantly enhanced antibody penetration into the cytosol, in particular when fused before or after the hinge. Finally, we demonstrate that specific antibody binding to the cell surface target is necessary for efficient cell penetration of the CPP-antibody fusions. This study provides a solid basis for further exploration of therapeutic antibodies for intracellular targets.

Isoform-specific interactions of the von Hippel-Lindau tumor suppressor protein.

Deregulation of the von Hippel-Lindau tumor suppressor protein (pVHL) is considered one of the main causes for malignant renal clear-cell carcinoma (ccRCC) insurgence. In human, pVHL exists in two isoforms, pVHL19 and pVHL30 respectively, displaying comparable tumor suppressor abilities. Mutations of the p53 tumor suppressor gene have been also correlated with ccRCC insurgence and ineffectiveness of treatment. A recent proteomic analysis linked full length pVHL30 with p53 pathway regulation through complex formation with the p14ARF oncosuppressor. The alternatively spliced pVHL19, missing the first 53 residues, lacks this interaction and suggests an asymmetric function of the two pVHL isoforms. Here, we present an integrative bioinformatics and experimental characterization of the pVHL oncosuppressor isoforms. Predictions of the pVHL30 N-terminus three-dimensional structure suggest that it may exist as an ensemble of structured and disordered forms. The results were used to guide Yeast two hybrid experiments to highlight isoform-specific binding properties. We observed that the physical pVHL/p14ARF interaction is specifically mediated by the 53 residue long pVHL30 N-terminal region, suggesting that this N-terminus acts as a further pVHL interaction interface. Of note, we also observed that the shorter pVHL19 isoform shows an unexpected high tendency to form homodimers, suggesting an additional isoform-specific binding specialization.

Structural protein reorganization and fold emergence investigated through amino acid sequence permutations.

Correlation between random amino acid sequences and protein folds suggests that proteins autonomously evolved the most stable folds, with stability and function evolving subsequently, suggesting the existence of common protein ancestors from which all modern proteins evolved. To test this hypothesis, we shuffled the sequences of 10 natural proteins and obtained 40 different and apparently unrelated folds. Our results suggest that shuffled sequences are sufficiently stable and may act as a basis to evolve functional proteins. The common secondary structure of modern proteins is well represented by a small set of permuted sequences, which also show the emergence of intrinsic disorder and aggregation-prone stretches of the polypeptide chain.

Evaluation of the steric impact of flavin adenine dinucleotide in Drosophila melanogaster cryptochrome function.

Photoreceptors are crucial components for circadian rhythm entrainment in animals, plants, fungi and cyanobacteria. Cryptochromes (CRYs) are flavin adenine dinucleotide (FAD) containing photoreceptors, and FAD is responsible for signal transduction, in contrast to photolyases where it promotes DNA-damage repair. In this work, we investigated an alternative role for FAD in CRY. We analyzed the Drosophila melanogaster CRY crystal structure by means of molecular dynamics, elucidating how this large co-factor within the receptor could be crucial for CRY structural stability. The co-factor appears indeed to improve receptor motility, providing steric hindrance. Moreover, multiple sequence alignments revealed that conserved motifs in the C-terminal tail could be necessary for functional stability.

Design and analysis of a Petri net model of the Von Hippel-Lindau (VHL) tumor suppressor interaction network.

Von Hippel-Lindau (VHL) syndrome is a hereditary condition predisposing to the development of different cancer forms, related to germline inactivation of the homonymous tumor suppressor pVHL. The best characterized function of pVHL is the ubiquitination dependent degradation of Hypoxia Inducible Factor (HIF) via the proteasome. It is also involved in several cellular pathways acting as a molecular hub and interacting with more than 200 different proteins. Molecular details of pVHL plasticity remain in large part unknown. Here, we present a novel manually curated Petri Net (PN) model of the main pVHL functional pathways. The model was built using functional information derived from the literature. It includes all major pVHL functions and is able to credibly reproduce VHL syndrome at the molecular level. The reliability of the PN model also allowed in silico knockout experiments, driven by previous model analysis. Interestingly, PN analysis suggests that the variability of different VHL manifestations is correlated with the concomitant inactivation of different metabolic pathways.

In silico investigation of PHD-3 specific HIF1-α proline 567 hydroxylation: a new player in the VHL/HIF-1α interaction pathway?

Hypoxia inducible factor 1α (HIF-1α) regulates oxygen homeostasis in the cell through a sensing mechanism involving its hydroxylation and binding to the von Hippel-Lindau (VHL) tumor suppressor. This mechanism is mediated through hydroxylation of HIF-1α proline 564, although in vitro tests have previously shown an alternative hydroxylation at proline 567 by PHD-3. Here, molecular dynamics simulations were used to investigate the structural effect of this alternative hydroxylation. A specific hydrogen bond network rearrangement and improved electrostatic energy for hydroxylated P567 are compatible with an increase in HIF-1α binding affinity. Sequence analysis also confirms P567 to be vastly conserved during evolution, indicating a possible role for this alternative, PHD-3 driven, post translational modification in pVHL-HIF-1α complex formation.