Robust Prediction of Relative Binding Energies for Protein-Protein Complex Mutations Using Free Energy Perturbation Calculations.

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how Protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Cis inhibition of co-expressed DIPs and Dprs shapes neural development.

In Drosophila , two interacting adhesion protein families, Dprs and DIPs, coordinate the assembly of neural networks. While intercellular DIP/Dpr interactions have been well characterized, DIPs and Dprs are often co-expressed within the same cells, raising the question as to whether they also interact in cis . We show, in cultured cells and in vivo, that DIP-α and DIP-δ can interact in cis with their ligands, Dpr6/10 and Dpr12, respectively. When co-expressed in cis with their cognate partners, these Dprs regulate the extent of trans binding, presumably through competitive cis interactions. We demonstrate the neurodevelopmental effects of cis inhibition in fly motor neurons and in the mushroom body. We further show that a long disordered region of DIP-α at the C-terminus is required for cis but not trans interactions, likely because it alleviates geometric constraints on cis binding. Thus, the balance between cis and trans interactions plays a role in controlling neural development.

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations.

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Free Energy Perturbation Calculations of Mutation Effects on SARS-CoV-2 RBD::ACE2 Binding Affinity.

The strength of binding between human angiotensin converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of viral spike protein plays a role in the transmissibility of the SARS-CoV-2 virus. In this study we focus on a subset of RBD mutations that have been frequently observed in infected individuals and probe binding affinity changes to ACE2 using surface plasmon resonance (SPR) measurements and free energy perturbation (FEP) calculations. Our SPR results are largely in accord with previous studies but discrepancies do arise due to differences in experimental methods and to protocol differences even when a single method is used. Overall, we find that FEP performance is superior to that of other computational approaches examined as determined by agreement with experiment and, in particular, by its ability to identify stabilizing mutations. Moreover, the calculations successfully predict the observed cooperative stabilization of binding by the Q498R N501Y double mutant present in Omicron variants and offer a physical explanation for the underlying mechanism. Overall, our results suggest that despite the significant computational cost, FEP calculations may offer an effective strategy to understand the effects of interfacial mutations on protein-protein binding affinities and, hence, in a variety of practical applications such as the optimization of neutralizing antibodies.

Structures of LRP2 reveal a molecular machine for endocytosis.

The low-density lipoprotein (LDL) receptor-related protein 2 (LRP2 or megalin) is representative of the phylogenetically conserved subfamily of giant LDL receptor-related proteins, which function in endocytosis and are implicated in diseases of the kidney and brain. Here, we report high-resolution cryoelectron microscopy structures of LRP2 isolated from mouse kidney, at extracellular and endosomal pH. The structures reveal LRP2 to be a molecular machine that adopts a conformation for ligand binding at the cell surface and for ligand shedding in the endosome. LRP2 forms a homodimer, the conformational transformation of which is governed by pH-sensitive sites at both homodimer and intra-protomer interfaces. A subset of LRP2 deleterious missense variants in humans appears to impair homodimer assembly. These observations lay the foundation for further understanding the function and mechanism of LDL receptors and implicate homodimerization as a conserved feature of the LRP receptor subfamily.

Structural insights into broadly neutralizing antibodies elicited by hybrid immunity against SARS-CoV-2.

ABSTRACTIncreasing spread by SARS-CoV-2 Omicron variants challenges existing vaccines and broadly reactive neutralizing antibodies (bNAbs) against COVID-19. Here we determine the diversity, potency, breadth and structural insights of bNAbs derived from memory B cells of BNT162b2-vaccinee after homogeneous Omicron BA.1 breakthrough infection. The infection activates diverse memory B cell clonotypes for generating potent class I/II and III bNAbs with new epitopes mapped to the receptor-binding domain (RBD). The top eight bNAbs neutralize wildtype and BA.1 potently but display divergent IgH/IgL sequences and neuralization profiles against other variants of concern (VOCs). Two of them (P2D9 and P3E6) belonging to class III NAbs display comparable potency against BA.4/BA.5, although structural analysis reveals distinct modes of action. P3E6 neutralizes all variants tested through a unique bivalent interaction with two RBDs. Our findings provide new insights into hybrid immunity on BNT162b2-induced diverse memory B cells in response to Omicron breakthrough infection for generating diverse bNAbs with distinct structural basis.

Functional properties of the spike glycoprotein of the emerging SARS-CoV-2 variant B.1.1.529.

The recently emerged B.1.1.529 (Omicron) severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant has a highly divergent spike (S) glycoprotein. We compared the functional properties of B.1.1.529 BA.1 S with those of previous globally prevalent SARS-CoV-2 variants, D614G and B.1.617.2. Relative to these variants, B.1.1.529 S exhibits decreases in processing, syncytium formation, virion incorporation, and ability to mediate infection of cells with high TMPRSS2 expression. B.1.1.529 and B.1.617.2 S glycoproteins bind ACE2 with higher affinity than D614G S. The unliganded B.1.1.529 S trimer is less stable at low temperatures than the other SARS-CoV-2 Ss, a property related to its more "open" S conformation. Upon ACE2 binding, the B.1.1.529 S trimer sheds S1 at 37°C, but not at 0°C. B.1.1.529 pseudoviruses are relatively resistant to neutralization by sera from patients with coronavirus disease 2019 (COVID-19) and vaccinees. These properties of the B.1.1.529 S glycoprotein likely influence the transmission, cytopathic effects, and immune evasion of this emerging variant.

Affinity requirements for control of synaptic targeting and neuronal cell survival by heterophilic IgSF cell adhesion molecules.

Neurons in the developing brain express many different cell adhesion molecules (CAMs) on their surfaces. CAM-binding affinities can vary by more than 200-fold, but the significance of these variations is unknown. Interactions between the immunoglobulin superfamily CAM DIP-α and its binding partners, Dpr10 and Dpr6, control synaptic targeting and survival of Drosophila optic lobe neurons. We design mutations that systematically change interaction affinity and analyze function in vivo. Reducing affinity causes loss-of-function phenotypes whose severity scales with the magnitude of the change. Synaptic targeting is more sensitive to affinity reduction than is cell survival. Increasing affinity rescues neurons that would normally be culled by apoptosis. By manipulating CAM expression together with affinity, we show that the key parameter controlling circuit assembly is surface avidity, which is the strength of adherence between cell surfaces. We conclude that CAM binding affinities and expression levels are finely tuned for function during development.

How clustered protocadherin binding specificity is tuned for neuronal self-/nonself-recognition.

The stochastic expression of fewer than 60 clustered protocadherin (cPcdh) isoforms provides diverse identities to individual vertebrate neurons and a molecular basis for self-/nonself-discrimination. cPcdhs form chains mediated by alternating cis and trans interactions between apposed membranes, which has been suggested to signal self-recognition. Such a mechanism requires that cPcdh cis dimers form promiscuously to generate diverse recognition units, and that trans interactions have precise specificity so that isoform mismatches terminate chain growth. However, the extent to which cPcdh interactions fulfill these requirements has not been definitively demonstrated. Here, we report biophysical experiments showing that cPcdh cis interactions are promiscuous, but with preferences favoring formation of heterologous cis dimers. Trans homophilic interactions are remarkably precise, with no evidence for heterophilic interactions between different isoforms. A new C-type cPcdh crystal structure and mutagenesis data help to explain these observations. Overall, the interaction characteristics we report for cPcdhs help explain their function in neuronal self-/nonself-discrimination.

Synaptogenic activity of the axon guidance molecule Robo2 underlies hippocampal circuit function.

Synaptic connectivity within adult circuits exhibits a remarkable degree of cellular and subcellular specificity. We report that the axon guidance receptor Robo2 plays a role in establishing synaptic specificity in hippocampal CA1. In vivo, Robo2 is present and required postsynaptically in CA1 pyramidal neurons (PNs) for the formation of excitatory (E) but not inhibitory (I) synapses, specifically in proximal but not distal dendritic compartments. In vitro approaches show that the synaptogenic activity of Robo2 involves a trans-synaptic interaction with presynaptic Neurexins, as well as binding to its canonical extracellular ligand Slit. In vivo 2-photon Ca2+ imaging of CA1 PNs during spatial navigation in awake behaving mice shows that preventing Robo2-dependent excitatory synapse formation cell autonomously during development alters place cell properties of adult CA1 PNs. Our results identify a trans-synaptic complex linking the establishment of synaptic specificity to circuit function.

Neutralizing antibody 5-7 defines a distinct site of vulnerability in SARS-CoV-2 spike N-terminal domain.

Antibodies that potently neutralize SARS-CoV-2 target mainly the receptor-binding domain or the N-terminal domain (NTD). Over a dozen potently neutralizing NTD-directed antibodies have been studied structurally, and all target a single antigenic supersite in NTD (site 1). Here, we report the cryo-EM structure of a potent NTD-directed neutralizing antibody 5-7, which recognizes a site distinct from other potently neutralizing antibodies, inserting a binding loop into an exposed hydrophobic pocket between the two sheets of the NTD ß sandwich. Interestingly, this pocket was previously identified as the binding site for hydrophobic molecules, including heme metabolites, but we observe that their presence does not substantially impede 5-7 recognition. Mirroring its distinctive binding, antibody 5-7 retains neutralization potency with many variants of concern (VOCs). Overall, we reveal that a hydrophobic pocket in NTD proposed for immune evasion can be used by the immune system for recognition.

Paired heavy- and light-chain signatures contribute to potent SARS-CoV-2 neutralization in public antibody responses.

Understanding mechanisms of protective antibody recognition can inform vaccine and therapeutic strategies against SARS-CoV-2. We report a monoclonal antibody, 910-30, targeting the SARS-CoV-2 receptor-binding site for ACE2 as a member of a public antibody response encoded by IGHV3-53/IGHV3-66 genes. Sequence and structural analyses of 910-30 and related antibodies explore how class recognition features correlate with SARS-CoV-2 neutralization. Cryo-EM structures of 910-30 bound to the SARS-CoV-2 spike trimer reveal binding interactions and its ability to disassemble spike. Despite heavy-chain sequence similarity, biophysical analyses of IGHV3-53/3-66-encoded antibodies highlight the importance of native heavy:light pairings for ACE2-binding competition and SARS-CoV-2 neutralization. We develop paired heavy:light class sequence signatures and determine antibody precursor prevalence to be ∼1 in 44,000 human B cells, consistent with public antibody identification in several convalescent COVID-19 patients. These class signatures reveal genetic, structural, and functional immune features that are helpful in accelerating antibody-based medical interventions for SARS-CoV-2.

Neutralizing antibody 5-7 defines a distinct site of vulnerability in SARS-CoV-2 spike N-terminal domain.

Antibodies that potently neutralize SARS-CoV-2 target mainly the receptor-binding domain or the N-terminal domain (NTD). Over a dozen potently neutralizing NTD-directed antibodies have been studied structurally, and all target a single antigenic supersite in NTD (site 1). Here we report the 3.7 Å resolution cryo-EM structure of a potent NTD-directed neutralizing antibody 5-7, which recognizes a site distinct from other potently neutralizing antibodies, inserting a binding loop into an exposed hydrophobic pocket between the two sheets of the NTD ß-sandwich. Interestingly, this pocket has been previously identified as the binding site for hydrophobic molecules including heme metabolites, but we observe their presence to not substantially impede 5-7 recognition. Mirroring its distinctive binding, antibody 5-7 retains a distinctive neutralization potency with variants of concern (VOC). Overall, we reveal a hydrophobic pocket in NTD proposed for immune evasion can actually be used by the immune system for recognition. HIGHLIGHTS: Cryo-EM structure of neutralizing antibody 5-7 in complex with SARS CoV-2 spike5-7 recognizes NTD outside of the previously identified antigenic supersite5-7 binds to a site known to accommodate numerous hydrophobic ligandsStructural basis of 5-7 neutralization tolerance to some variants of concern.

Modular basis for potent SARS-CoV-2 neutralization by a prevalent VH1-2-derived antibody class.

Antibodies with heavy chains that derive from the VH1-2 gene constitute some of the most potent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-neutralizing antibodies yet identified. To provide insight into whether these genetic similarities inform common modes of recognition, we determine the structures of the SARS-CoV-2 spike in complex with three VH1-2-derived antibodies: 2-15, 2-43, and H4. All three use VH1-2-encoded motifs to recognize the receptor-binding domain (RBD), with heavy-chain N53I-enhancing binding and light-chain tyrosines recognizing F486RBD. Despite these similarities, class members bind both RBD-up and -down conformations of the spike, with a subset of antibodies using elongated CDRH3s to recognize glycan N343 on a neighboring RBD-a quaternary interaction accommodated by an increase in RBD separation of up to 12 Å. The VH1-2 antibody class, thus, uses modular recognition encoded by modular genetic elements to effect potent neutralization, with the VH-gene component specifying recognition of RBD and the CDRH3 component specifying quaternary interactions.

Potent SARS-CoV-2 neutralizing antibodies directed against spike N-terminal domain target a single supersite.

Numerous antibodies that neutralize SARS-CoV-2 have been identified, and these generally target either the receptor-binding domain (RBD) or the N-terminal domain (NTD) of the viral spike. While RBD-directed antibodies have been extensively studied, far less is known about NTD-directed antibodies. Here, we report cryo-EM and crystal structures for seven potent NTD-directed neutralizing antibodies in complex with spike or isolated NTD. These structures defined several antibody classes, with at least one observed in multiple convalescent donors. The structures revealed that all seven antibodies target a common surface, bordered by glycans N17, N74, N122, and N149. This site-formed primarily by a mobile ß-hairpin and several flexible loops-was highly electropositive, located at the periphery of the spike, and the largest glycan-free surface of NTD facing away from the viral membrane. Thus, in contrast to neutralizing RBD-directed antibodies that recognize multiple non-overlapping epitopes, potent NTD-directed neutralizing antibodies appear to target a single supersite.

Paired heavy and light chain signatures contribute to potent SARS-CoV-2 neutralization in public antibody responses.

Understanding protective mechanisms of antibody recognition can inform vaccine and therapeutic strategies against SARS-CoV-2. We discovered a new antibody, 910-30, that targets the SARS-CoV-2 ACE2 receptor binding site as a member of a public antibody response encoded by IGHV3-53/IGHV3-66 genes. We performed sequence and structural analyses to explore how antibody features correlate with SARS-CoV-2 neutralization. Cryo-EM structures of 910-30 bound to the SARS-CoV-2 spike trimer revealed its binding interactions and ability to disassemble spike. Despite heavy chain sequence similarity, biophysical analyses of IGHV3-53/3-66 antibodies highlighted the importance of native heavy:light pairings for ACE2 binding competition and for SARS-CoV-2 neutralization. We defined paired heavy:light sequence signatures and determined antibody precursor prevalence to be ~1 in 44,000 human B cells, consistent with public antibody identification in several convalescent COVID-19 patients. These data reveal key structural and functional neutralization features in the IGHV3-53/3-66 public antibody class to accelerate antibody-based medical interventions against SARS-CoV-2. HIGHLIGHTS: A molecular study of IGHV3-53/3-66 public antibody responses reveals critical heavy and light chain features for potent neutralizationCryo-EM analyses detail the structure of a novel public antibody class member, antibody 910-30, in complex with SARS-CoV-2 spike trimerCryo-EM data reveal that 910-30 can both bind assembled trimer and can disassemble the SARS-CoV-2 spikeSequence-structure-function signatures defined for IGHV3-53/3-66 class antibodies including both heavy and light chainsIGHV3-53/3-66 class precursors have a prevalence of 1:44,000 B cells in healthy human antibody repertoires.

Structural Basis of Antibody Conformation and Stability Modulation by Framework Somatic Hypermutation.

Accumulation of somatic hypermutation (SHM) is the primary mechanism to enhance the binding affinity of antibodies to antigens in vivo. However, the structural basis of the effects of many SHMs remains elusive. Here, we integrated atomistic molecular dynamics (MD) simulation and data mining to build a high-throughput structural bioinformatics pipeline to study the effects of individual and combination SHMs on antibody conformation, flexibility, stability, and affinity. By applying this pipeline, we characterized a common mechanism of modulation of heavy-light pairing orientation by frequent SHMs at framework positions 39H, 91H, 38L, and 87L through disruption of a conserved hydrogen-bond network. Q39LH alone and in combination with light chain framework 4 (FWR4L) insertions further modulated the elbow angle between variable and constant domains of many antibodies, resulting in improved binding affinity for a subset of anti-HIV-1 antibodies. Q39LH also alleviated aggregation induced by FWR4L insertion, suggesting remote epistasis between these SHMs. Altogether, this study provides tools and insights for understanding antibody affinity maturation and for engineering functionally improved antibodies.

Cryo-EM Structures of SARS-CoV-2 Spike without and with ACE2 Reveal a pH-Dependent Switch to Mediate Endosomal Positioning of Receptor-Binding Domains.

The SARS-CoV-2 spike employs mobile receptor-binding domains (RBDs) to engage the human ACE2 receptor and to facilitate virus entry, which can occur through low-pH-endosomal pathways. To understand how ACE2 binding and low pH affect spike conformation, we determined cryo-electron microscopy structures-at serological and endosomal pH-delineating spike recognition of up to three ACE2 molecules. RBDs freely adopted "up" conformations required for ACE2 interaction, primarily through RBD movement combined with smaller alterations in neighboring domains. In the absence of ACE2, single-RBD-up conformations dominated at pH 5.5, resolving into a solitary all-down conformation at lower pH. Notably, a pH-dependent refolding region (residues 824-858) at the spike-interdomain interface displayed dramatic structural rearrangements and mediated RBD positioning through coordinated movements of the entire trimer apex. These structures provide a foundation for understanding prefusion-spike mechanics governing endosomal entry; we suggest that the low pH all-down conformation potentially facilitates immune evasion from RBD-up binding antibody.

Ubiquitin-dependent regulation of a conserved DMRT protein controls sexually dimorphic synaptic connectivity and behavior.

Sex-specific synaptic connectivity is beginning to emerge as a remarkable, but little explored feature of animal brains. We describe here a novel mechanism that promotes sexually dimorphic neuronal function and synaptic connectivity in the nervous system of the nematode Caenorhabditis elegans. We demonstrate that a phylogenetically conserved, but previously uncharacterized Doublesex/Mab-3 related transcription factor (DMRT), dmd-4, is expressed in two classes of sex-shared phasmid neurons specifically in hermaphrodites but not in males. We find dmd-4 to promote hermaphrodite-specific synaptic connectivity and neuronal function of phasmid sensory neurons. Sex-specificity of DMD-4 function is conferred by a novel mode of posttranslational regulation that involves sex-specific protein stabilization through ubiquitin binding to a phylogenetically conserved but previously unstudied protein domain, the DMA domain. A human DMRT homolog of DMD-4 is controlled in a similar manner, indicating that our findings may have implications for the control of sexual differentiation in other animals as well.

Structure-Based Design with Tag-Based Purification and In-Process Biotinylation Enable Streamlined Development of SARS-CoV-2 Spike Molecular Probes.

Biotin-labeled molecular probes, comprising specific regions of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike, would be helpful in the isolation and characterization of antibodies targeting this recently emerged pathogen. Here, we design constructs incorporating an N-terminal purification tag, a site-specific protease-cleavage site, the probe region of interest, and a C-terminal sequence targeted by biotin ligase. Probe regions include full-length spike ectodomain as well as various subregions, and we also design mutants that eliminate recognition of the angiotensin-converting enzyme 2 (ACE2) receptor. Yields of biotin-labeled probes from transient transfection range from ∼0.5 mg/L for the complete ectodomain to >5 mg/L for several subregions. Probes are characterized for antigenicity and ACE2 recognition, and the structure of the spike ectodomain probe is determined by cryoelectron microscopy. We also characterize antibody-binding specificities and cell-sorting capabilities of the biotinylated probes. Altogether, structure-based design coupled to efficient purification and biotinylation processes can thus enable streamlined development of SARS-CoV-2 spike ectodomain probes.