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.

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.

Astrocyte morphogenesis requires self-recognition.

Self-recognition is a fundamental cellular process across evolution and forms the basis of neuronal self-avoidance1-4. Clustered protocadherins (Pcdh), comprising a large family of isoform-specific homophilic recognition molecules, play a pivotal role in neuronal self-avoidance required for mammalian brain development5-7. The probabilistic expression of different Pcdh isoforms confers unique identities upon neurons and forms the basis for neuronal processes to discriminate between self and non-self5,6,8. Whether this self-recognition mechanism exists in astrocytes, the other predominant cell type of the brain, remains unknown. Here, we report that a specific isoform in the Pcdhγ cluster, γC3, is highly enriched in human and murine astrocytes. Through genetic manipulation, we demonstrate that γC3 acts autonomously to regulate astrocyte morphogenesis in the mouse visual cortex. To determine if γC3 proteins act by promoting recognition between processes of the same astrocyte, we generated pairs of γC3 chimeric proteins capable of heterophilic binding to each other, but incapable of homophilic binding. Co-expressing complementary heterophilic binding isoform pairs in the same γC3 null astrocyte restored normal morphology. By contrast, chimeric γC3 proteins individually expressed in single γC3 null mutant astrocytes did not. These data establish that self-recognition is essential for astrocyte development in the mammalian brain and that, by contrast to neuronal self-recognition, a single Pcdh isoform is both necessary and sufficient for this process.

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.

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.

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.

DIP/Dpr interactions and the evolutionary design of specificity in protein families.

Differential binding affinities among closely related protein family members underlie many biological phenomena, including cell-cell recognition. Drosophila DIP and Dpr proteins mediate neuronal targeting in the fly through highly specific protein-protein interactions. We show here that DIPs/Dprs segregate into seven specificity subgroups defined by binding preferences between their DIP and Dpr members. We then describe a sequence-, structure- and energy-based computational approach, combined with experimental binding affinity measurements, to reveal how specificity is coded on the canonical DIP/Dpr interface. We show that binding specificity of DIP/Dpr subgroups is controlled by "negative constraints", which interfere with binding. To achieve specificity, each subgroup utilizes a different combination of negative constraints, which are broadly distributed and cover the majority of the protein-protein interface. We discuss the structural origins of negative constraints, and potential general implications for the evolutionary origins of binding specificity in multi-protein families.

Visualization of clustered protocadherin neuronal self-recognition complexes.

Neurite self-recognition and avoidance are fundamental properties of all nervous systems1. These processes facilitate dendritic arborization2,3, prevent formation of autapses4 and allow free interaction among non-self neurons1,2,4,5. Avoidance among self neurites is mediated by stochastic cell-surface expression of combinations of about 60 isoforms of α-, ß- and γ-clustered protocadherin that provide mammalian neurons with single-cell identities1,2,4-13. Avoidance is observed between neurons that express identical protocadherin repertoires2,5, and single-isoform differences are sufficient to prevent self-recognition10. Protocadherins form isoform-promiscuous cis dimers and isoform-specific homophilic trans dimers10,14-20. Although these interactions have previously been characterized in isolation15,17-20, structures of full-length protocadherin ectodomains have not been determined, and how these two interfaces engage in self-recognition between neuronal surfaces remains unknown. Here we determine the molecular arrangement of full-length clustered protocadherin ectodomains in single-isoform self-recognition complexes, using X-ray crystallography and cryo-electron tomography. We determine the crystal structure of the clustered protocadherin γB4 ectodomain, which reveals a zipper-like lattice that is formed by alternating cis and trans interactions. Using cryo-electron tomography, we show that clustered protocadherin γB6 ectodomains tethered to liposomes spontaneously assemble into linear arrays at membrane contact sites, in a configuration that is consistent with the assembly observed in the crystal structure. These linear assemblies pack against each other as parallel arrays to form larger two-dimensional structures between membranes. Our results suggest that the formation of ordered linear assemblies by clustered protocadherins represents the initial self-recognition step in neuronal avoidance, and thus provide support for the isoform-mismatch chain-termination model of protocadherin-mediated self-recognition, which depends on these linear chains11.

Protocadherin cis-dimer architecture and recognition unit diversity.

Clustered protocadherins (Pcdhs) mediate numerous neural patterning functions, including neuronal self-recognition and non-self-discrimination to direct self-avoidance among vertebrate neurons. Individual neurons stochastically express a subset of Pcdh isoforms, which assemble to form a stochastic repertoire of cis-dimers. We describe the structure of a PcdhγB7 cis-homodimer, which includes the membrane-proximal extracellular cadherin domains EC5 and EC6. The structure is asymmetric with one molecule contributing interface surface from both EC5 and EC6, and the other only from EC6. Structural and sequence analyses suggest that all Pcdh isoforms will dimerize through this interface. Site-directed mutants at this interface interfere with both Pcdh cis-dimerization and cell surface transport. The structure explains the known restrictions of cis-interactions of some Pcdh isoforms, including α-Pcdhs, which cannot form homodimers. The asymmetry of the interface approximately doubles the size of the recognition repertoire, and restrictions on cis-interactions among Pcdh isoforms define the limits of the Pcdh recognition unit repertoire.

γ-Protocadherin structural diversity and functional implications.

Stochastic cell-surface expression of α-, ß-, and γ-clustered protocadherins (Pcdhs) provides vertebrate neurons with single-cell identities that underlie neuronal self-recognition. Here we report crystal structures of ectodomain fragments comprising cell-cell recognition regions of mouse γ-Pcdhs γA1, γA8, γB2, and γB7 revealing trans-homodimers, and of C-terminal ectodomain fragments from γ-Pcdhs γA4 and γB2, which depict cis-interacting regions in monomeric form. Together these structures span the entire γ-Pcdh ectodomain. The trans-dimer structures reveal determinants of γ-Pcdh isoform-specific homophilic recognition. We identified and structurally mapped cis-dimerization mutations to the C-terminal ectodomain structures. Biophysical studies showed that Pcdh ectodomains from γB-subfamily isoforms formed cis dimers, whereas γA isoforms did not, but both γA and γB isoforms could interact in cis with α-Pcdhs. Together, these data show how interaction specificity is distributed over all domains of the γ-Pcdh trans interface, and suggest that subfamily- or isoform-specific cis-interactions may play a role in the Pcdh-mediated neuronal self-recognition code.

Molecular basis of sidekick-mediated cell-cell adhesion and specificity.

Sidekick (Sdk) 1 and 2 are related immunoglobulin superfamily cell adhesion proteins required for appropriate synaptic connections between specific subtypes of retinal neurons. Sdks mediate cell-cell adhesion with homophilic specificity that underlies their neuronal targeting function. Here we report crystal structures of Sdk1 and Sdk2 ectodomain regions, revealing similar homodimers mediated by the four N-terminal immunoglobulin domains (Ig1-4), arranged in a horseshoe conformation. These Ig1-4 horseshoes interact in a novel back-to-back orientation in both homodimers through Ig1:Ig2, Ig1:Ig1 and Ig3:Ig4 interactions. Structure-guided mutagenesis results show that this canonical dimer is required for both Sdk-mediated cell aggregation (via trans interactions) and Sdk clustering in isolated cells (via cis interactions). Sdk1/Sdk2 recognition specificity is encoded across Ig1-4, with Ig1-2 conferring the majority of binding affinity and differential specificity. We suggest that competition between cis and trans interactions provides a novel mechanism to sharpen the specificity of cell-cell interactions.

Structural Basis of Diverse Homophilic Recognition by Clustered α- and ß-Protocadherins.

Clustered protocadherin proteins (α-, ß-, and γ-Pcdhs) provide a high level of cell-surface diversity to individual vertebrate neurons, engaging in highly specific homophilic interactions to mediate important roles in mammalian neural circuit development. How Pcdhs bind homophilically through their extracellular cadherin (EC) domains among dozens of highly similar isoforms has not been determined. Here, we report crystal structures for extracellular regions from four mouse Pcdh isoforms (α4, α7, ß6, and ß8), revealing a canonical head-to-tail interaction mode for homophilic trans dimers comprising primary intermolecular EC1:EC4 and EC2:EC3 interactions. A subset of trans interface residues exhibit isoform-specific conservation, suggesting roles in recognition specificity. Mutation of these residues, along with trans-interacting partner residues, altered the specificities of Pcdh interactions. Together, these data show how sequence variation among Pcdh isoforms encodes their diverse strict homophilic recognition specificities, which are required for their key roles in neural circuit assembly.

Molecular logic of neuronal self-recognition through protocadherin domain interactions.

Self-avoidance, a process preventing interactions of axons and dendrites from the same neuron during development, is mediated in vertebrates through the stochastic single-neuron expression of clustered protocadherin protein isoforms. Extracellular cadherin (EC) domains mediate isoform-specific homophilic binding between cells, conferring cell recognition through a poorly understood mechanism. Here, we report crystal structures for the EC1-EC3 domain regions from four protocadherin isoforms representing the α, ß, and γ subfamilies. All are rod shaped and monomeric in solution. Biophysical measurements, cell aggregation assays, and computational docking reveal that trans binding between cells depends on the EC1-EC4 domains, which interact in an antiparallel orientation. We also show that the EC6 domains are required for the formation of cis-dimers. Overall, our results are consistent with a model in which protocadherin cis-dimers engage in a head-to-tail interaction between EC1-EC4 domains from apposed cell surfaces, possibly forming a zipper-like protein assembly, and thus providing a size-dependent self-recognition mechanism.

Crystal structures of Drosophila N-cadherin ectodomain regions reveal a widely used class of Ca²+-free interdomain linkers.

Vertebrate classical cadherins mediate selective calcium-dependent cell adhesion by mechanisms now understood at the atomic level. However, structures and adhesion mechanisms of cadherins from invertebrates, which are highly divergent yet function in similar roles, remain unknown. Here we present crystal structures of three- and four-tandem extracellular cadherin (EC) domain segments from Drosophila N-cadherin (DN-cadherin), each including the predicted N-terminal EC1 domain (denoted EC1') of the mature protein. While the linker regions for the EC1'-EC2' and EC3'-EC4' pairs display binding of three Ca(2+) ions similar to that of vertebrate cadherins, domains EC2' and EC3' are joined in a "kinked" orientation by a previously uncharacterized Ca(2+)-free linker. Biophysical analysis demonstrates that a construct containing the predicted N-terminal nine EC domains of DN-cadherin forms homodimers with affinity similar to vertebrate classical cadherins, whereas deleting the ninth EC domain ablates dimerization. These results suggest that, unlike their vertebrate counterparts, invertebrate cadherins may utilize multiple EC domains to form intercellular adhesive bonds. Sequence analysis reveals that similar Ca(2+)-free linkers are widely distributed in the ectodomains of both vertebrate and invertebrate cadherins.