The high-resolution NMR structure of the R21A Spc-SH3:P41 complex: understanding the determinants of binding affinity by comparison with Abl-SH3

Investigating vulnerability of the conserved SARS-CoV-2 spike's heptad repeat 2 as target for fusion inhibitors using chimeric miniproteins

Inhibition of SARS-CoV-2 membrane fusion is a highly desired target to combat COVID-19. The interaction between the spike's heptad repeat (HR) regions 1 (HR1) and 2 (HR2) is a crucial step during the fusion process and these highly conserved HR regions constitute attractive targets for fusion inhibitors. However, the relative importance of each subregion of the long HR1-HR2 interface for viral inhibition remains unclear. Here, we designed, produced, and characterized a series of chimeric miniproteins that mimic two different half subdomains of HR1. The proteins were designed as single polypeptide chains that spontaneously fold into antiparallel trimeric helical bundles aimed at structurally imitate the molecular surface of each HR1 half subregion. All the miniproteins folded stably as helical structures and could bind complementary HR2 peptides with moderate affinity. However, only the miniproteins mimicking the N-terminal HR1 half subdomain, but not those imitating C-terminal one, could inhibit cell infection by SARS-COV-2 real viruses in cell cultures. Most interestingly, the inhibitory activity of the miniproteins correlated with their structural stability, but not with their relative binding affinity for HR2 peptides. These results are highly relevant for designing more focused and active fusion inhibitors targeting the highly conserved HR2 region of the Spike.

Exploring Highly Conserved Regions of SARS-CoV-2 Spike S2 Subunit as Targets for Fusion Inhibition Using Chimeric Proteins

Since the beginning of the COVID-19 pandemic, considerable efforts have been made to develop protective vaccines against SARS-CoV-2 infection. However, immunity tends to decline within a few months, and new virus variants are emerging with increased transmissibility and capacity to evade natural or vaccine-acquired immunity. Therefore, new robust strategies are needed to combat SARS-CoV-2 infection. The viral spike composed of S1 and S2 subunits mediates viral attachment and membrane fusion to infect the host cell. In this process, interaction between the highly conserved heptad repeat 1 and 2 regions (HR1 and HR2) of S2 is crucial and for this reason; these regions are promising targets to fight SARS-CoV-2. Here, we describe the design and characterization of chimeric proteins that structurally imitate the S2 HR1 region in a trimeric coiled-coil conformation. We biophysically characterized the proteins and determined their capacity to bind the HR2 region, as well as their inhibitory activity of SARS-CoV-2 infection in vitro. HR1 mimetic proteins showed conformational heterogeneity and a propensity to form oligomers. Moreover, their structure is composed of subdomains with varied stability. Interestingly, the full HR1 proteins showed high affinity for HR2-derived peptides and SARS-CoV-2 inhibitory activity, whereas smaller proteins mimicking HR1 subdomains had a decreased affinity for their complementary HR2 region and did not inhibit the virus. The results provide insight into effective strategies to create mimetic proteins with broad inhibitory activity and therapeutic potential against SARS-CoV-2.

Non-Stationary Complementary Non-Uniform Sampling (NOSCO NUS) for Fast Acquisition of Serial 2D NMR Titration Data

NMR spectroscopy offers unique benefits for ligand binding studies on isotopically labelled target proteins. These benefits include atomic resolution, direct distinction of binding sites and modes, a lowest detectable affinity limit, and function independent setup. Yet, retracing protein signal assignments from apo to holo states to derive exact dissociation constants and chemical shift perturbation amplitudes (for ligand docking and structure‐based optimization) requires lengthy titration series of 2D heteronuclear correlation spectra at variable ligand concentration that may exceed the protein's lifetime and available spectrometer time. We present a novel method to overcome this critical limitation, based on non‐stationary complementary non‐uniform sampling (NOSCO NUS) combined with a robust particle swarm optimization algorithm. We illustrate its potential in two challenging studies with very distinct protein sizes and binding affinities, showing that NOSCO NUS can reduce measurement times by an order of magnitude to make such highly informative NMR titration studies more broadly feasible.

Mechanistic Insights into Microsecond Time-Scale Motion of Solid Proteins Using Complementary 15N and 1H Relaxation Dispersion Techniques

NMR relaxation dispersion methods provide a holistic way to observe microsecond time-scale protein backbone motion both in solution and in the solid state. Different nuclei (¹H and ¹⁵N) and different relaxation dispersion techniques (Bloch-McConnell and near-rotary-resonance) give complementary information about the amplitudes and time scales of the conformational dynamics and provide comprehensive insights into the mechanistic details of the structural rearrangements. In this paper, we exemplify the benefits of the combination of various solution- and solid-state relaxation dispersion methods on a microcrystalline protein (α-spectrin SH3 domain), for which we are able to identify and model the functionally relevant conformational rearrangements around the ligand recognition loop occurring on multiple microsecond time scales. The observed loop motions suggest that the SH3 domain exists in a binding-competent conformation in dynamic equilibrium with a sterically impaired ground-state conformation both in solution and in crystalline form. This inherent plasticity between the interconverting macrostates is compatible with a conformational-preselection model and provides new insights into the recognition mechanisms of SH3 domains.

Exploring Multiple Binding Modes Using Confined Replica Exchange Molecular Dynamics

Molecular docking is extensively applied to determine the position of a ligand on its receptor despite the rather poor correspondence between docking scores and experimental binding affinities found in several studies, especially for systems structurally unrelated with those used in the scoring functions' training sets. Here, we present a method for the prediction of binding modes and binding free energies, which uses replica exchange molecular dynamics in combination with a receptor-shaped piecewise potential, confining the ligand in the proximity of the receptor surface and limiting the accessible conformational space of interest. We assess our methodology with a set of protein receptor-ligand test cases. In every case studied, the method is able to locate the ligand on the experimentally known receptor binding site, and it gives as output the binding free energy. The added value of our approach with respect to other available methods is that it quickly performs a conformational space search, providing a set of bound (or unbound) configurations, which can be used to determine phenomenological structural and energetic properties of an experimental binding state as a result of contributions provided by diversified multiple binding poses.

From Binding-Induced Dynamic Effects in SH3 Structures to Evolutionary Conserved Sectors

Src Homology 3 domains are ubiquitous small interaction modules known to act as docking sites and regulatory elements in a wide range of proteins. Prior experimental NMR work on the SH3 domain of Src showed that ligand binding induces long-range dynamic changes consistent with an induced fit mechanism. The identification of the residues that participate in this mechanism produces a chart that allows for the exploration of the regulatory role of such domains in the activity of the encompassing protein. Here we show that a computational approach focusing on the changes in side chain dynamics through ligand binding identifies equivalent long-range effects in the Src SH3 domain. Mutation of a subset of the predicted residues elicits long-range effects on the binding energetics, emphasizing the relevance of these positions in the definition of intramolecular cooperative networks of signal transduction in this domain. We find further support for this mechanism through the analysis of seven other publically available SH3 domain structures of which the sequences represent diverse SH3 classes. By comparing the eight predictions, we find that, in addition to a dynamic pathway that is relatively conserved throughout all SH3 domains, there are dynamic aspects specific to each domain and homologous subgroups. Our work shows for the first time from a structural perspective, which transduction mechanisms are common between a subset of closely related and distal SH3 domains, while at the same time highlighting the differences in signal transduction that make each family member unique. These results resolve the missing link between structural predictions of dynamic changes and the domain sectors recently identified for SH3 domains through sequence analysis.

Small protein domains as Src Homology 3 often act as docking sites and serve as regulatory elements. To understand their role in the regulation of a protein’s activity, one needs to understand how their backbone and sidechain dynamics are affected when binding to peptides. We have therefore computationally analyzed eight different SH3 domain structures, predicting dynamical effects induced by binding through our MCIT approach that has been shown to correlate well with experimental data. We show first that binding the Src SH3 domain triggers a particular cascade of dynamic effects, which are compatible with an induced fit mechanism reported before. We then combined the predictions for the eight SH3 domains into different consensus models, with the aim of analyzing, for the first time from a structural perspective, commonalities and differences in the transduction mechanisms among these SH3 domains. These consensus results are, on one hand, in agreement with the domain sectors recently identified for the entire family of SH3 domains. On the other hand, they reveal also that differences exist between the different subgroups that were studied here, requiring extensive experimental investigations of the importance of these differences for the proteins wherein these SH3 domains can be found.

Three easy pieces

BACKGROUND

Differential scanning calorimetry is a powerful method that provides a complete thermodynamic characterization of the stability of a protein as a function of temperature. There are, however, circumstances that preclude a complete analysis of DSC data. The most common ones are irreversible denaturation transitions or transitions that take place at temperatures that are beyond the temperature limit of the instrument. Even for a protein that undergoes reversible thermal denaturation, the extrapolation of the thermodynamic data to lower temperatures, usually 25°C, may become unreliable due to difficulties in the determination of ΔCp.

METHODS

The combination of differential scanning calorimetry and isothermal chemical denaturation allows reliable thermodynamic analysis of protein stability under less than ideal conditions.

RESULTS AND CONCLUSIONS

This paper demonstrates how DSC can be used in combination with chemical denaturation to address three different scenarios: 1) estimation of an accurate ΔCp value for a reversible denaturation using as a test system the envelope HIV-1 glycoprotein gp120; 2) determination of the Gibbs energy of stability in the region in which thermal denaturation is irreversible using HEW lysozyme at different pH values; and, 3) determination of Gibbs energy of stability for a thermostable protein, thermolysin.

Insights into Substrate Specificity of NlpC/P60 Cell Wall Hydrolases Containing Bacterial SH3 Domains

Bacterial SH3 (SH3b) domains are commonly fused with papain-like Nlp/P60 cell wall hydrolase domains. To understand how the modular architecture of SH3b and NlpC/P60 affects the activity of the catalytic domain, three putative NlpC/P60 cell wall hydrolases were biochemically and structurally characterized. These enzymes all have γ-d-Glu-A₂pm (A₂pm is diaminopimelic acid) cysteine amidase (or dl-endopeptidase) activities but with different substrate specificities. One enzyme is a cell wall lysin that cleaves peptidoglycan (PG), while the other two are cell wall recycling enzymes that only cleave stem peptides with an N-terminal l-Ala. Their crystal structures revealed a highly conserved structure consisting of two SH3b domains and a C-terminal NlpC/P60 catalytic domain, despite very low sequence identity. Interestingly, loops from the first SH3b domain dock into the ends of the active site groove of the catalytic domain, remodel the substrate binding site, and modulate substrate specificity. Two amino acid differences at the domain interface alter the substrate binding specificity in favor of stem peptides in recycling enzymes, whereas the SH3b domain may extend the peptidoglycan binding surface in the cell wall lysins. Remarkably, the cell wall lysin can be converted into a recycling enzyme with a single mutation.

Peptidoglycan is a meshlike polymer that envelops the bacterial plasma membrane and bestows structural integrity. Cell wall lysins and recycling enzymes are part of a set of lytic enzymes that target covalent bonds connecting the amino acid and amino sugar building blocks of the PG network. These hydrolases are involved in processes such as cell growth and division, autolysis, invasion, and PG turnover and recycling. To avoid cleavage of unintended substrates, these enzymes have very selective substrate specificities. Our biochemical and structural analysis of three modular NlpC/P60 hydrolases, one lysin, and two recycling enzymes, show that they may have evolved from a common molecular architecture, where the substrate preference is modulated by local changes. These results also suggest that new pathways for recycling PG turnover products, such as tracheal cytotoxin, may have evolved in bacteria in the human gut microbiome that involve NlpC/P60 cell wall hydrolases.

Multimeric and differential binding of CIN85/CD2AP with two atypical proline-rich sequences from CD2 and Cbl-b*

The CD2AP (CD2-associated protein) and CIN85 (Cbl-interacting protein of 85 kDa) adaptor proteins each employ three Src homology 3 (SH3) domains to cluster protein partners and ensure efficient signal transduction and down-regulation of tyrosine kinase receptors. Using NMR, isothermal titration calorimetry and small-angle X-ray scattering methods, we have characterized several binding modes of the N-terminal SH3 domain (SH3A) of CD2AP and CIN85 with two natural atypical proline-rich regions in CD2 (cluster of differentiation 2) and Cbl-b (Casitas B-lineage lymphoma), and compared these data with previous studies and published crystal structures. Our experiments show that the CD2AP-SH3A domain forms a type II dimer with CD2 and both type I and type II dimeric complexes with Cbl-b. Like CD2AP, the CIN85-SH3A domain forms a type II complex with CD2, but a trimeric complex with Cbl-b, whereby the type I and II interactions take place at the same time. Together, these results explain how multiple interactions among similar SH3 domains and ligands produce a high degree of diversity in tyrosine kinase, cell adhesion or T-cell signaling pathways.

A unified conformational selection and induced fit approach to protein-peptide docking

Protein-peptide interactions are vital for the cell. They mediate, inhibit or serve as structural components in nearly 40% of all macromolecular interactions, and are often associated with diseases, making them interesting leads for protein drug design. In recent years, large-scale technologies have enabled exhaustive studies on the peptide recognition preferences for a number of peptide-binding domain families. Yet, the paucity of data regarding their molecular binding mechanisms together with their inherent flexibility makes the structural prediction of protein-peptide interactions very challenging. This leaves flexible docking as one of the few amenable computational techniques to model these complexes. We present here an ensemble, flexible protein-peptide docking protocol that combines conformational selection and induced fit mechanisms. Starting from an ensemble of three peptide conformations (extended, a-helix, polyproline-II), flexible docking with HADDOCK generates 79.4% of high quality models for bound/unbound and 69.4% for unbound/unbound docking when tested against the largest protein-peptide complexes benchmark dataset available to date. Conformational selection at the rigid-body docking stage successfully recovers the most relevant conformation for a given protein-peptide complex and the subsequent flexible refinement further improves the interface by up to 4.5 Å interface RMSD. Cluster-based scoring of the models results in a selection of near-native solutions in the top three for ∼75% of the successfully predicted cases. This unified conformational selection and induced fit approach to protein-peptide docking should open the route to the modeling of challenging systems such as disorder-order transitions taking place upon binding, significantly expanding the applicability limit of biomolecular interaction modeling by docking.

High-resolution crystal structure of spectrin SH3 domain fused with a proline-rich peptide

A new chimeric protein, named WT-CIIA, was designed by connecting the proline-rich decapeptide PPPVPPYSAG to the C-terminus of the alpha-spectrin SH3 domain through a natural twelve-residue linker to obtain a single-chain model that would imitate intramolecular SH3-ligand interaction. The crystal structure of this fusion protein was determined at 1.7 Å resolution. The asymmetric unit of the crystal contained two SH3 globules contacting with one PPPVPPY fragment located between them. The domains are related by the two-fold non-crystallographic axis and the ligand lies in two opposite orientations with respect to the conservative binding sites of SH3 domains.

Conformational and functional analysis of molecular dynamics trajectories by self-organising maps

BACKGROUND

Molecular dynamics (MD) simulations are powerful tools to investigate the conformational dynamics of proteins that is often a critical element of their function. Identification of functionally relevant conformations is generally done clustering the large ensemble of structures that are generated. Recently, Self-Organising Maps (SOMs) were reported performing more accurately and providing more consistent results than traditional clustering algorithms in various data mining problems. We present a novel strategy to analyse and compare conformational ensembles of protein domains using a two-level approach that combines SOMs and hierarchical clustering.

RESULTS

The conformational dynamics of the α-spectrin SH3 protein domain and six single mutants were analysed by MD simulations. The Cα's Cartesian coordinates of conformations sampled in the essential space were used as input data vectors for SOM training, then complete linkage clustering was performed on the SOM prototype vectors. A specific protocol to optimize a SOM for structural ensembles was proposed: the optimal SOM was selected by means of a Taguchi experimental design plan applied to different data sets, and the optimal sampling rate of the MD trajectory was selected. The proposed two-level approach was applied to single trajectories of the SH3 domain independently as well as to groups of them at the same time. The results demonstrated the potential of this approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space in a simple 2D visualisation, as well as of effectively highlighting differences in the conformational dynamics directly related to biological functions.

CONCLUSIONS

The use of a two-level approach combining SOMs and hierarchical clustering for conformational analysis of structural ensembles of proteins was proposed. It can easily be extended to other study cases and to conformational ensembles from other sources.

High-resolution structure of an alpha-spectrin SH3-domain mutant with a redesigned hydrophobic core

The alpha-spectrin SH3 domain (Spc-SH3) is a small modular domain which has been broadly used as a model protein in folding studies and these studies have sometimes been supported by structural information obtained from the coordinates of Spc-SH3 mutants. The structure of B5/D48G, a multiple mutant designed to improve the hydrophobic core and as a consequence the protein stability, has been solved at 1 A resolution. The crystals belonged to the orthorhombic space group P2(1)2(1)2(1), with unit-cell parameters a=24.79, b=37.23, c=62.95 A. This mutant also bears a D48G substitution in the distal loop and this mutation has also been reported to increase the stability of the protein by itself. The structure of the B5/D48G mutant shows a highly packed hydrophobic core and a more ordered distal loop compared with previous Spc-SH3 structures.

Solution structure of histone chaperone ANP32B: interaction with core histones H3-H4 through its acidic concave domain

Eukaryotic gene expression is regulated by histone deposition onto and eviction from nucleosomes, which are mediated by several chromatin-modulating factors. Among them, histone chaperones are key factors that facilitate nucleosome assembly. Acidic nuclear phosphoprotein 32B (ANP32B) belongs to the ANP32 family, which shares N-terminal leucine-rich repeats (LRRs) and a C-terminal variable anionic region. The C-terminal region functions as an inhibitor of histone acetylation, but the functional roles of the LRR domain in chromatin regulation have remained elusive. Here, we report that the LRR domain of ANP32B possesses histone chaperone activity and forms a curved structure with a parallel beta-sheet on the concave side and mostly helical elements on the convex side. Our analyses revealed that the interaction of ANP32B with the core histones H3-H4 occurs on its concave side, and both the acidic and hydrophobic residues that compose the concave surface are critical for histone binding. These results provide a structural framework for understanding the functional mechanisms of acidic histone chaperones.

Role of interfacial water molecules in proline-rich ligand recognition by the Src homology 3 domain of Abl

The interaction of Abl-Src homology 3 domain (SH3) with the high affinity peptide p41 is the most notable example of the inconsistency existing between the currently accepted description of SH3 complexes and their binding thermodynamic signature. We had previously hypothesized that the presence of interfacial water molecules is partially responsible for this thermodynamic behavior. We present here a thermodynamic, structural, and molecular dynamics simulation study of the interaction of p41 with Abl-SH3 and a set of mutants designed to alter the water-mediated interaction network. Our results provide a detailed description of the dynamic properties of the interfacial water molecules and a molecular interpretation of the thermodynamic effects elicited by the mutations in terms of the modulation of the water-mediated hydrogen bond network. In the light of these results, a new dual binding mechanism is proposed that provides a better description of proline-rich ligand recognition by Abl-SH3 and that has important implications for rational design.

Binding mechanism of an SH3 domain studied by NMR and ITC

Complexes between Src-homology 3 domains and proline-rich target peptides can have lifetimes on the order of milliseconds, making them too short-lived for kinetic characterization by conventional methods. Nuclear magnetic resonance (NMR) dynamics experiments are ideally suited to study such rapid binding equilibria, and additionally provide information on partly bound intermediate states. We used NMR together with isothermal titration calorimetry (ITC) to characterize the interaction of the SH3 domain from the Fyn tyrosine kinase with a 12-residue peptide at temperatures between 10 and 50 degrees C. NMR data at all temperatures are consistent with an effectively two-state binding reaction, such that any intermediates are either very weakly populated or exchange extremely rapidly with the free or bound forms. Dissociation rate constants, determined by CPMG and ZZ-exchange NMR experiments, range from k(off)(10 degrees C) = 4.5 s(-1) to k(off)(50 degrees C) = 331 s(-1). ITC data at all temperatures follow a simple two-state interaction model. Binding is favored enthalpically, with a dissociation enthalpy, DeltaH(D)(30 degrees C) = 15.4 kcal mol(-1), and disfavored entropically, with a dissociation entropy, DeltaS(D)(30 degrees C) = 20.0 cal mol(-1) K(-1). The free protein and peptide have significantly higher heat capacity than the bound complex, DeltaC(p) = 352 cal mol(-1) K(-1), which is consistent with the largely hydrophobic character of the binding interface. An Eyring plot of k(off) values gives an activation enthalpy of dissociation, DeltaH(D)(double dagger)(30 degrees C) = 19.3 kcal mol(-1) and exhibits slight curvature consistent with the ITC-derived value of DeltaC(p). The curvature suggests that nonpolar residues of the hydrophobic interface are solvated in the transition state for dissociation. Association rate constants were calculated as k(on) = k(off)/K(D), and range from k(on)(10 degrees C) = 1.03 x 10(8) M(-1) s(-1) to k(on)(50 degrees C) = 2.0 x 10(8) M(-1) s(-1), with an apparent activation enthalpy, DeltaH(A)(double dagger) = 3.4 kcal mol(-1). Both the magnitudes and temperature dependence of k(on) values are consistent with a diffusion-limited association mechanism. The combination of NMR and ITC data sheds light on how the Fyn tyrosine kinase is activated by binding to proline-rich targets, and represents a powerful approach for characterizing transient protein/ligand interactions.

Evaluation of folding co-operativity of a chimeric protein based on the molecular recognition between polyproline ligands and SH3 domains

In previous work, we designed a chimeric protein, named SPCp41, to evaluate the thermodynamics of the interaction between SH3 domains and proline-rich ligands by combining thermal unfolding measurements and mutagenesis. Here, we have investigated the energetic integrity of the chain extension corresponding to the ligand sequence into the native structure, since the opposite will produce changes in the folding mechanism of the SH3 domain that may give rise to undesirable contributions to the thermodynamic parameters. We have analysed the folding-unfolding kinetics under standard conditions (50 mM phosphate pH 7). Kinetic evolutions are well described by a bi-exponential where, on top of the main kinetic phase, a low-populated slower phase appears as a consequence of cis-trans isomerisation of Pro39, as demonstrated by the influence of prolyl isomerases and by mutational analysis. There is also a burst phase possibly due to a productive formation of some helical ensembles. The main evolution, accounting for the true folding kinetics of SPCp41, can be considered as a two-state process, where the folding transition state produces essentially the same picture shown by the circular permutant S19-P20s (the 'nucleus' of the design) and the ligand will dock at the latter stages of the two-state process. Thus, all conclusions argue in favour of the effectiveness of SPCp41 to study energetic, dynamic and structural aspects of SH3-ligand interactions.

A survey of the year 2007 literature on applications of isothermal titration calorimetry

Elucidation of the energetic principles of binding affinity and specificity is a central task in many branches of current sciences: biology, medicine, pharmacology, chemistry, material sciences, etc. In biomedical research, integral approaches combining structural information with in-solution biophysical data have proved to be a powerful way toward understanding the physical basis of vital cellular phenomena. Isothermal titration calorimetry (ITC) is a valuable experimental tool facilitating quantification of the thermodynamic parameters that characterize recognition processes involving biomacromolecules. The method provides access to all relevant thermodynamic information by performing a few experiments. In particular, ITC experiments allow to by-pass tedious and (rarely precise) procedures aimed at determining the changes in enthalpy and entropy upon binding by van't Hoff analysis. Notwithstanding limitations, ITC has now the reputation of being the "gold standard" and ITC data are widely used to validate theoretical predictions of thermodynamic parameters, as well as to benchmark the results of novel binding assays. In this paper, we discuss several publications from 2007 reporting ITC results. The focus is on applications in biologically oriented fields. We do not intend a comprehensive coverage of all newly accumulated information. Rather, we emphasize work which has captured our attention with originality and far-reaching analysis, or else has provided ideas for expanding the potential of the method.