FuzPred: a web server for the sequence-based prediction of the context-dependent binding modes of proteins.

Proteins form complex interactions in the cellular environment to carry out their functions. They exhibit a wide range of binding modes depending on the cellular conditions, which result in a variety of ordered or disordered assemblies. To help rationalise the binding behavior of proteins, the FuzPred server predicts their sequence-based binding modes without specifying their binding partners. The binding mode defines whether the bound state is formed through a disorder-to-order transition resulting in a well-defined conformation, or through a disorder-to-disorder transition where the binding partners remain conformationally heterogeneous. To account for the context-dependent nature of the binding modes, the FuzPred method also estimates the multiplicity of binding modes, the likelihood of sampling multiple binding modes. Protein regions with a high multiplicity of binding modes may serve as regulatory sites or hot-spots for structural transitions in the assembly. To facilitate the interpretation of the predictions, protein regions with different interaction behaviors can be visualised on protein structures generated by AlphaFold. The FuzPred web server (https://fuzpred.bio.unipd.it) thus offers insights into the structural and dynamical changes of proteins upon interactions and contributes to development of structure-function relationships under a variety of cellular conditions.

Alternatively spliced exon regulates context-dependent MEF2D higher-order assembly during myogenesis.

During muscle cell differentiation, the alternatively spliced, acidic ß-domain potentiates transcription of Myocyte-specific Enhancer Factor 2 (Mef2D). Sequence analysis by the FuzDrop method indicates that the ß-domain can serve as an interaction element for Mef2D higher-order assembly. In accord, we observed Mef2D mobile nuclear condensates in C2C12 cells, similar to those formed through liquid-liquid phase separation. In addition, we found Mef2D solid-like aggregates in the cytosol, the presence of which correlated with higher transcriptional activity. In parallel, we observed a progress in the early phase of myotube development, and higher MyoD and desmin expression. In accord with our predictions, the formation of aggregates was promoted by rigid ß-domain variants, as well as by a disordered ß-domain variant, capable of switching between liquid-like and solid-like higher-order states. Along these lines, NMR and molecular dynamics simulations corroborated that the ß-domain can sample both ordered and disordered interactions leading to compact and extended conformations. These results suggest that ß-domain fine-tunes Mef2D higher-order assembly to the cellular context, which provides a platform for myogenic regulatory factors and the transcriptional apparatus during the developmental process.

Protein interactions: anything new?

How do proteins interact in the cellular environment? Which interactions stabilize liquid-liquid phase separated condensates? Are the concepts, which have been developed for specific protein complexes also applicable to higher-order assemblies? Recent discoveries prompt for a universal framework for protein interactions, which can be applied across the scales of protein communities. Here, we discuss how our views on protein interactions have evolved from rigid structures to conformational ensembles of proteins and discuss the open problems, in particular related to biomolecular condensates. Protein interactions have evolved to follow changes in the cellular environment, which manifests in multiple modes of interactions between the same partners. Such cellular context-dependence requires multiplicity of binding modes (MBM) by sampling multiple minima of the interaction energy landscape. We demonstrate that the energy landscape framework of protein folding can be applied to explain this phenomenon, opening a perspective toward a physics-based, universal model for cellular protein behaviors.

Conformational buffering underlies functional selection in intrinsically disordered protein regions.

Many disordered proteins conserve essential functions in the face of extensive sequence variation, making it challenging to identify the mechanisms responsible for functional selection. Here we identify the molecular mechanism of functional selection for the disordered adenovirus early gene 1A (E1A) protein. E1A competes with host factors to bind the retinoblastoma (Rb) protein, subverting cell cycle regulation. We show that two binding motifs tethered by a hypervariable disordered linker drive picomolar affinity Rb binding and host factor displacement. Compensatory changes in amino acid sequence composition and sequence length lead to conservation of optimal tethering across a large family of E1A linkers. We refer to this compensatory mechanism as conformational buffering. We also detect coevolution of the motifs and linker, which can preserve or eliminate the tethering mechanism. Conformational buffering and motif-linker coevolution explain robust functional encoding within hypervariable disordered linkers and could underlie functional selection of many disordered protein regions.

Insights on cross-species transmission of SARS-CoV-2 from structural modeling.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global pandemic that has infected more than 31 million people in more than 180 countries worldwide. Like other coronaviruses, SARS-CoV-2 is thought to have been transmitted to humans from wild animals. Given the scale and widespread geographical distribution of the current pandemic and confirmed cases of cross-species transmission, the question of the extent to which this transmission is possible emerges, as well as what molecular features distinguish susceptible from non-susceptible animal species. Here, we investigated the structural properties of several ACE2 orthologs bound to the SARS-CoV-2 spike protein. We found that species known not to be susceptible to SARS-CoV-2 infection have non-conservative mutations in several ACE2 amino acid residues that disrupt key polar and charged contacts with the viral spike protein. Our models also allow us to predict affinity-enhancing mutations that could be used to design ACE2 variants for therapeutic purposes. Finally, our study provides a blueprint for modeling viral-host protein interactions and highlights several important considerations when designing these computational studies and analyzing their results.

Structural models of human ACE2 variants with SARS-CoV-2 Spike protein for structure-based drug design.

Emergence of coronaviruses poses a threat to global health and economy. The current outbreak of SARS-CoV-2 has infected more than 28,000,000 people and killed more than 915,000. To date, there is no treatment for coronavirus infections, making the development of therapies to prevent future epidemics of paramount importance. To this end, we collected information regarding naturally-occurring variants of the Angiotensin-converting enzyme 2 (ACE2), an epithelial receptor that both SARS-CoV and SARS-CoV-2 use to enter the host cells. We built 242 structural models of variants of human ACE2 bound to the receptor binding domain (RBD) of the SARS-CoV-2 surface spike glycoprotein (S protein) and refined their interfaces with HADDOCK. Our dataset includes 140 variants of human ACE2 representing missense mutations found in genome-wide studies, 39 mutants with reported effects on the recognition of the RBD, and 63 predictions after computational alanine scanning mutagenesis of ACE2-RBD interface residues. This dataset will help accelerate the design of therapeutics against SARS-CoV-2, as well as contribute to prevention of possible future coronaviruses outbreaks.

Insights on cross-species transmission of SARS-CoV-2 from structural modeling.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global pandemic that has infected more than 14 million people in more than 180 countries worldwide. Like other coronaviruses, SARS-CoV-2 is thought to have been transmitted to humans from wild animals. Given the scale and widespread geographical distribution of the current pandemic, the question emerges whether human-to-animal transmission is possible and if so, which animal species are most at risk. Here, we investigated the structural properties of several ACE2 orthologs bound to the SARS-CoV-2 spike protein. We found that species known not to be susceptible to SARS-CoV-2 infection have non-conservative mutations in several ACE2 amino acid residues that disrupt key polar and charged contacts with the viral spike protein. Our models also predict affinity-enhancing mutations that could be used to design ACE2 variants for therapeutic purposes. Finally, our study provides a blueprint for modeling viral-host protein interactions and highlights several important considerations when designing these computational studies and analyzing their results.

p15PAF binding to PCNA modulates the DNA sliding surface.

p15PAF is an oncogenic intrinsically disordered protein that regulates DNA replication and lesion bypass by interacting with the human sliding clamp PCNA. In the absence of DNA, p15PAF traverses the PCNA ring via an extended PIP-box that contacts the sliding surface. Here, we probed the atomic-scale structure of p15PAF-PCNA-DNA ternary complexes. Crystallography and MD simulations show that, when p15PAF occupies two subunits of the PCNA homotrimer, DNA within the ring channel binds the unoccupied subunit. The structure of PCNA-bound p15PAF in the absence and presence of DNA is invariant, and solution NMR confirms that DNA does not displace p15PAF from the ring wall. Thus, p15PAF reduces the available sliding surfaces of PCNA, and may function as a belt that fastens the DNA to the clamp during synthesis by the replicative polymerase (pol δ). This constraint, however, may need to be released for efficient DNA lesion bypass by the translesion synthesis polymerase (pol η). Accordingly, our biochemical data show that p15PAF impairs primer synthesis by pol η-PCNA holoenzyme against both damaged and normal DNA templates. In light of our findings, we discuss the possible mechanistic roles of p15PAF in DNA replication and suppression of DNA lesion bypass.

Structural basis of human PCNA sliding on DNA.

Sliding clamps encircle DNA and tether polymerases and other factors to the genomic template. However, the molecular mechanism of clamp sliding on DNA is unknown. Using crystallography, NMR and molecular dynamics simulations, here we show that the human clamp PCNA recognizes DNA through a double patch of basic residues within the ring channel, arranged in a right-hand spiral that matches the pitch of B-DNA. We propose that PCNA slides by tracking the DNA backbone via a 'cogwheel' mechanism based on short-lived polar interactions, which keep the orientation of the clamp invariant relative to DNA. Mutation of residues at the PCNA-DNA interface has been shown to impair the initiation of DNA synthesis by polymerase δ (pol δ). Therefore, our findings suggest that a clamp correctly oriented on DNA is necessary for the assembly of a replication-competent PCNA-pol δ holoenzyme.

ASC Pyrin Domain Self-associates and Binds NLRP3 Protein Using Equivalent Binding Interfaces.

Death domain superfamily members typically act as adaptors mediating in the assembly of supramolecular complexes with critical apoptosis and inflammation functions. These modular proteins consist of death domains, death effector domains, caspase recruitment domains, and pyrin domains (PYD). Despite the high structural similarity among them, only homotypic interactions participate in complex formation, suggesting that subtle factors differentiate each interaction type. It is thus critical to identify these factors as an essential step toward the understanding of the molecular basis of apoptosis and inflammation. The proteins apoptosis-associated speck-like protein containing a CARD (ASC) and NLRP3 play key roles in the regulation of apoptosis and inflammation through self-association and protein-protein interactions mediated by their PYDs. To better understand the molecular basis of their function, we have characterized ASC and NLRP3 PYD self-association and their intermolecular interaction by solution NMR spectroscopy and analytical ultracentrifugation. We found that ASC self-associates and binds NLRP3 PYD through equivalent protein regions, with higher binding affinity for the latter. These regions are located at opposite sides of the protein allowing multimeric complex formation previously shown in ASC PYD fibril assemblies. We show that NLRP3 PYD coexists in solution as a monomer and highly populated large-order oligomerized species. Despite this, we determined its monomeric three-dimensional solution structure by NMR and characterized its binding to ASC PYD. Using our novel structural data, we propose molecular models of ASC·ASC and ASC·NLRP3 PYD early supramolecular complexes, providing new insights into the molecular mechanisms of inflammasome and apoptosis signaling.

Intrinsic order and disorder in the bcl-2 member harakiri: insights into its proapoptotic activity.

Harakiri is a BH3-only member of the Bcl-2 family that localizes in membranes and induces cell death by binding to prosurvival Bcl-x(L) and Bcl-2. The cytosolic domain of Harakiri is largely disorder with residual α-helical conformation according to previous structural studies. As these helical structures could play an important role in Harakiri's function, we have used NMR and circular dichroism to fully characterize them at the residue-atomic level. In addition, we report structural studies on a peptide fragment spanning Harakiri's C-terminal hydrophobic sequence, which potentially operates as a transmembrane domain. We initially checked by enzyme immunoassays and NMR that peptides encompassing different lengths of the cytosolic domain are functional as they bind Bcl-x(L) and Bcl-2. The structural data in water indicate that the α-helical conformation is restricted to a 25-residue segment comprising the BH3 domain. However, structure calculation was precluded because of insufficient NMR restraints. To bypass this problem we used alcohol-water mixture to increase structure population and confirmed by NMR that the conformation in both milieus is equivalent. The resulting three-dimensional structure closely resembles that of peptides encompassing the BH3 domain of BH3-only members in complex with their prosurvival partners, suggesting that preformed structural elements in the disordered protein are central to binding. In contrast, the transmembrane domain forms in micelles a monomeric α-helix with a population close to 100%. Its three-dimensional structure here reported reveals features that explain its function as membrane anchor. Altogether these results are used to propose a tentative structural model of how Harakiri works.

Characterization of a novel interaction between Bcl-2 members Diva and Harakiri.

Interactions within proteins of the Bcl-2 family are key in the regulation of apoptosis. The death-inducing members control apoptotic mechanisms partly by antagonizing the prosurvival proteins through heterodimer formation. Structural and biophysical studies on these complexes are providing important clues to understand their function. To help improve our knowledge on protein-protein interactions within the Bcl-2 family we have studied the binding between two of its members: mouse Diva and human Harakiri. Diva has been shown to perform both prosurvival and killing activity. In contrast, Harakiri induces cell death by interacting with antiapoptotic Bcl-2 members. Here we show using ELISA and NMR that Diva and Harakiri can interact in vitro. Combining the NMR data with the previously reported three-dimensional structure of Diva we find that Harakiri binds to a specific region in Diva. This interacting surface is equivalent to the known binding area of prosurvival Bcl-2 members from the reported structures of the complexes, suggesting that Diva could function at the structural level similarly to the antiapoptotic proteins of the Bcl-2 family. We illustrate this result by building a structural model of the heterodimer using molecular docking and the NMR data as restraints. Moreover, combining circular dichroism and NMR we also show that Harakiri is largely unstructured with residual (13%) α-helical conformation. This result agrees with intrinsic disorder previously observed in other Bcl-2 members. In addition, Harakiri constructs of different length were studied to identify the region critical for the interaction. Differential affinity for Diva of these constructs suggests that the amino acid sequence flanking the interacting region could play an important role in binding.