De novo design of high-affinity binders of bioactive helical peptides

Many peptide hormones form an α-helix on binding their receptors1-4, and sensitive methods for their detection could contribute to better clinical management of disease5. De novo protein design can now generate binders with high affinity and specificity to structured proteins6,7. However, the design of interactions between proteins and short peptides with helical propensity is an unmet challenge. Here we describe parametric generation and deep learning-based methods for designing proteins to address this challenge. We show that by extending RFdiffusion8 to enable binder design to flexible targets, and to refining input structure models by successive noising and denoising (partial diffusion), picomolar-affinity binders can be generated to helical peptide targets by either refining designs generated with other methods, or completely de novo starting from random noise distributions without any subsequent experimental optimization. The RFdiffusion designs enable the enrichment and subsequent detection of parathyroid hormone and glucagon by mass spectrometry, and the construction of bioluminescence-based protein biosensors. The ability to design binders to conformationally variable targets, and to optimize by partial diffusion both natural and designed proteins, should be broadly useful.

Interplay of Electrostatics and Hydrophobic Effects in the Metamorphic Protein Human Lymphotactin

The human lymphotactin (hLtn) is a protein that features two native states both of which are physiologically relevant: it is a monomer (hLtn10) at 10 °C with 200 mM salt and a dimer (hLtn40) at 40 °C and without salt. Here we focus on the networks of electrostatic and hydrophobic interactions that display substantial changes upon the conversion from hLtn10 to hLtn40 since they are expected to modulate the relative stability of the two folds. In addition to the Arg 23-Arg 43 interaction discussed in previous work, we find several other like-charge pairs that are likely important to the stability of hLtn10. Free energy perturbation calculations are carried out to explicitly evaluate the contribution of the Arg 23-Arg 43 interaction to the hLtn10 stability. hLtn40 features a larger number of salt bridges, and a set of hydrophobic residues undergo major changes in the solvent accessible surface area between hLtn10 and hLtn40, pointing to their importance to the relative stability of the two folds. We also discuss the use of explicit and implicit solvent simulations for characterizing the conformational ensembles under different solution conditions.

Identification of the Ah-receptor structural determinants for ligand preferences

The aryl hydrocarbon receptor (AHR) is a transcription factor that responds to diverse ligands and plays a critical role in toxicology, immune function, and cardiovascular physiology. The structural basis of the AHR for ligand promiscuity and preferences is critical for understanding AHR function. Based on the structure of a closely related protein HIF2α, we modeled the AHR ligand binding domain (LBD) bound to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and benzo(a)pyrene (BaP) and identified residues that control ligand preferences by shape and H-bond potential. Mutations to these residues, particularly Q377 and G298, resulted in robust and opposite changes in the potency of TCDD and BaP and up to a 20-fold change in the ratio of TCDD/BaP efficacy. The model also revealed a flexible "belt" structure; molecular dynamic (MD) simulation suggested that the "belt" and several other structural elements in the AHR-LBD are more flexible than HIF2α and likely contribute to ligand promiscuity. Molecular docking of TCDD congeners to a model of human AHR-LBD ranks their binding affinity similar to experimental ranking of their toxicity. Our study reveals key structural basis for prediction of toxicity and understanding the AHR signaling through diverse ligands.

Conformational control of the binding of the transactivation domain of the MLL protein and c-Myb to the KIX domain of CREB

The KIX domain of CBP is a transcriptional coactivator. Concomitant binding to the activation domain of proto-oncogene protein c-Myb and the transactivation domain of the trithorax group protein mixed lineage leukemia (MLL) transcription factor lead to the biologically active ternary MLL∶KIX∶c-Myb complex which plays a role in Pol II-mediated transcription. The binding of the activation domain of MLL to KIX enhances c-Myb binding. Here we carried out molecular dynamics (MD) simulations for the MLL∶KIX∶c-Myb ternary complex, its binary components and KIX with the goal of providing a mechanistic explanation for the experimental observations. The dynamic behavior revealed that the MLL binding site is allosterically coupled to the c-Myb binding site. MLL binding redistributes the conformational ensemble of KIX, leading to higher populations of states which favor c-Myb binding. The key element in the allosteric communication pathways is the KIX loop, which acts as a control mechanism to enhance subsequent binding events. We tested this conclusion by in silico mutations of loop residues in the KIX∶MLL complex and by comparing wild type and mutant dynamics through MD simulations. The loop assumed MLL binding conformation similar to that observed in the KIX∶c-Myb state which disfavors the allosteric network. The coupling with c-Myb binding site faded, abolishing the positive cooperativity observed in the presence of MLL. Our major conclusion is that by eliciting a loop-mediated allosteric switch between the different states following the binding events, transcriptional activation can be regulated. The KIX system presents an example how nature makes use of conformational control in higher level regulation of transcriptional activity and thus cellular events.

CBP (CREB-binding protein) is a transcriptional regulator of RNA polymerase II-mediated transcription. KIX is a domain of CBP. KIX binding to the activation domain of proto-oncogene protein c-Myb and the transactivation domain of the trithorax group protein mixed lineage leukemia (MLL) transcription factor forms the biologically active ternary MLL∶KIX∶c-Myb complex. This complex has been shown to play a key role in Pol II-mediated transcription. Experimental data show that the binding of the activation domain of MLL to KIX enhances c-Myb binding. Here, we studied the c-Myb∶KIX∶MLL ternary structure and models based on this structure, KIX-only, KIX-MLL, and c-Myb-KIX, by explicit solvent molecular dynamics (MD) simulations. MD simulations can help in figuring out allosteric events and thus functional mechanisms. Our in silico analysis of the dynamic behavior indicated that when MLL is bound, there is allosteric communication between the two KIX binding sites. We observed a shift in the conformational ensemble of KIX upon the binding of the MLL activation domain, leading to a pre-organization of the KIX∶c-Myb binding site and explaining the enhanced c-Myb binding observed experimentally. On the other hand, this is not the case if c-Myb binds to KIX, which suggests that KIX regulation is under conformational control.