Simulating and analysing configurational landscapes of protein-protein contact formation

Rescaling protein-protein interactions improves Martini 3 for flexible proteins in solution

Multidomain proteins with flexible linkers and disordered regions play important roles in many cellular processes, but characterizing their conformational ensembles is difficult. We have previously shown that the coarse-grained model, Martini 3, produces too compact ensembles in solution, that may in part be remedied by strengthening protein-water interactions. Here, we show that decreasing the strength of protein-protein interactions leads to improved agreement with experimental data on a wide set of systems. We show that the 'symmetry' between rescaling protein-water and protein-protein interactions breaks down when studying interactions with or within membranes; rescaling protein-protein interactions better preserves the binding specificity of proteins with lipid membranes, whereas rescaling protein-water interactions preserves oligomerization of transmembrane helices. We conclude that decreasing the strength of protein-protein interactions improves the accuracy of Martini 3 for IDPs and multidomain proteins, both in solution and in the presence of a lipid membrane.

Unraveling the folding-assisted unbinding mechanism of TCF with its binding partner β-catenin

This study utilizes molecular dynamics simulations aided with multiple walker parallel bias metadynamics to investigate the TCF unbinding mechanism from the β-catenin interface. The results, consistent with experimental binding affinity calculations, unveil a folding-assisted unbinding mechanism.

Reconciling ASPP-p53 binding mode discrepancies through an ensemble binding framework that bridges crystallography and NMR data

ASPP2 and iASPP bind to p53 through their conserved ANK-SH3 domains to respectively promote and inhibit p53-dependent cell apoptosis. While crystallography has indicated that these two proteins employ distinct surfaces of their ANK-SH3 domains to bind to p53, solution NMR data has suggested similar surfaces. In this study, we employed multi-scale molecular dynamics (MD) simulations combined with free energy calculations to reconcile the discrepancy in the binding modes. We demonstrated that the binding mode based solely on a single crystal structure does not enable iASPP's RT loop to engage with p53's C-terminal linker-a verified interaction. Instead, an ensemble of simulated iASPP-p53 complexes facilitates this interaction. We showed that the ensemble-average inter-protein contacting residues and NMR-detected interfacial residues qualitatively overlap on ASPP proteins, and the ensemble-average binding free energies better match experimental KD values compared to single crystallgarphy-determined binding mode. For iASPP, the sampled ensemble complexes can be grouped into two classes, resembling the binding modes determined by crystallography and solution NMR. We thus propose that crystal packing shifts the equilibrium of binding modes towards the crystallography-determined one. Lastly, we showed that the ensemble binding complexes are sensitive to p53's intrinsically disordered regions (IDRs), attesting to experimental observations that these IDRs contribute to biological functions. Our results provide a dynamic and ensemble perspective for scrutinizing these important cancer-related protein-protein interactions (PPIs).

Specifying conformational heterogeneity of multi-domain proteins at atomic resolution

The conformational landscape of multi-domain proteins is inherently linked to their specific functions. This also holds for polyubiquitin chains that are assembled by two or more ubiquitin domains connected by a flexible linker thus showing a large interdomain mobility. However, molecular recognition and signal transduction are associated with particular conformational substates that are populated in solution. Here, we apply high-resolution NMR spectroscopy in combination with dual-scale MD simulations to explore the conformational space of K6-, K29-, and K33-linked diubiquitin molecules. The conformational ensembles are evaluated utilizing a paramagnetic cosolute reporting on solvent exposure plus a set of complementary NMR parameters. This approach unravels a conformational heterogeneity of diubiquitins and explains the diversity of structural models that have been determined for K6-, K29-, and K33-linked diubiquitins in free and ligand-bound states so far. We propose a general application of the approach developed here to demystify multi-domain proteins occurring in nature.

Combining molecular dynamics simulations and scoring method to computationally model ubiquitylated linker histones in chromatosomes

The chromatin in eukaryotic cells plays a fundamental role in all processes during a cell's life cycle. This nucleoprotein is normally tightly packed but needs to be unpacked for expression and division. The linker histones are critical for such packaging processes and while most experimental and simulation works recognize their crucial importance, the focus is nearly always set on the nucleosome as the basic chromatin building block. Linker histones can undergo several modifications, but only few studies on their ubiquitylation have been conducted. Mono-ubiquitylated linker histones (HUb), while poorly understood, are expected to influence DNA compaction. The size of ubiquitin and the globular domain of the linker histone are comparable and one would expect an increased disorder upon ubiquitylation of the linker histone. However, the formation of higher order chromatin is not hindered and ubiquitylation of the linker histone may even promote gene expression. Structural data on chromatosomes is rare and HUb has never been modeled in a chromatosome so far. Descriptions of the chromatin complex with HUb would greatly benefit from computational structural data. In this study we generate molecular dynamics simulation data for six differently linked HUb variants with the help of a sampling scheme tailored to drive the exploration of phase space. We identify conformational sub-states of the six HUb variants using the sketch-map algorithm for dimensionality reduction and iterative HDBSCAN for clustering on the excessively sampled, shallow free energy landscapes. We present a highly efficient geometric scoring method to identify sub-states of HUb that fit into the nucleosome. We predict HUb conformations inside a nucleosome using on-dyad and off-dyad chromatosome structures as reference and show that unbiased simulations of HUb produce significantly more fitting than non-fitting HUb conformations. A tetranucleosome array is used to show that ubiquitylation can even occur in chromatin without too much steric clashes.

Generating a conformational landscape of ubiquitin chains at atomistic resolution by back-mapping based sampling

Ubiquitin chains are flexible multidomain proteins that have important biological functions in cellular signalling. Computational studies with all-atom molecular dynamics simulations of the conformational spaces of polyubiquitins can be challenging due to the system size and a multitude of long-lived meta-stable states. Coarse graining is an efficient approach to overcome this problem-at the cost of losing high-resolution details. Recently, we proposed the back-mapping based sampling (BMBS) approach that reintroduces atomistic information into a given coarse grained (CG) sampling based on a two-dimensional (2D) projection of the conformational landscape, produces an atomistic ensemble and allows to systematically compare the ensembles at the two levels of resolution. Here, we apply BMBS to K48-linked tri-ubiquitin, showing its applicability to larger systems than those it was originally introduced on and demonstrating that the algorithm scales very well with system size. In an extension of the original BMBS we test three different seeding strategies, i.e. different approaches from where in the CG landscape atomistic trajectories are initiated. Furthermore, we apply a recently introduced conformational clustering algorithm to the back-mapped atomistic ensemble. Thus, we obtain insight into the structural composition of the 2D landscape and illustrate that the dimensionality reduction algorithm separates different conformational characteristics very well into different regions of the map. This cluster analysis allows us to show how atomistic trajectories sample conformational states, move through the projection space and in sum converge to an atomistic conformational landscape that slightly differs from the original CG map, indicating a correction of flaws in the CG template.

Improving Martini 3 for Disordered and Multidomain Proteins

Coarse-grained molecular dynamics simulations are a useful tool to determine conformational ensembles of proteins. Here, we show that the coarse-grained force field Martini 3 underestimates the global dimensions of intrinsically disordered proteins (IDPs) and multidomain proteins when compared with small-angle X-ray scattering (SAXS) data and that increasing the strength of protein-water interactions favors more expanded conformations. We find that increasing the strength of interactions between protein and water by ca. 10% results in improved agreement with the SAXS data for IDPs and multidomain proteins. We also show that this correction results in a more accurate description of self-association of IDPs and folded proteins and better agreement with paramagnetic relaxation enhancement data for most IDPs. While simulations with this revised force field still show deviations to experiments for some systems, our results suggest that it is overall a substantial improvement for coarse-grained simulations of soluble proteins.

Conformational ensembles of intrinsically disordered proteins and flexible multidomain proteins

Intrinsically disordered proteins (IDPs) and multidomain proteins with flexible linkers show a high level of structural heterogeneity and are best described by ensembles consisting of multiple conformations with associated thermodynamic weights. Determining conformational ensembles usually involves the integration of biophysical experiments and computational models. In this review, we discuss current approaches to determine conformational ensembles of IDPs and multidomain proteins, including the choice of biophysical experiments, computational models used to sample protein conformations, models to calculate experimental observables from protein structure, and methods to refine ensembles against experimental data. We also provide examples of recent applications of integrative conformational ensemble determination to study IDPs and multidomain proteins and suggest future directions for research in the field.

Multiscale simulations of protein and membrane systems

Classical multiscale simulations are perfectly suited to investigate biological soft matter systems. Owing to the bridging between microscopically realistic and lower-resolution models or the integration of a hierarchy of subsystems, one gets access to biologically relevant system sizes and timescales. In recent years, increasingly complex systems and processes have come into focus such as multidomain proteins, phase separation processes in biopolymer solutions, multicomponent biomembranes, or multiprotein complexes up to entire viruses. The review shows factors that have contributed to this progress - from improved models to machine-learning-based analysis and scale-bridging methods.

Non-Canonical Interaction between Calmodulin and Calcineurin Contributes to the Differential Regulation of Plant-Derived Calmodulins on Calcineurin

Calmodulin (CaM) serves as an important Ca2+ signaling hub that regulates many protein signaling pathways. Recently, it was demonstrated that plant CaM homologues can regulate mammalian targets, often in a manner that opposes the impact of the mammalian CaM (mCaM). However, the molecular basis of how CaM homologue mutations differentially impact target activation is unclear. To understand these mechanisms, we examined two CaM isoforms found in soybean plants that differentially regulate a mammalian target, calcineurin (CaN). These CaM isoforms, sCaM-1 and sCaM-4, share >90 and ∼78% identity with the mCaM, respectively, and activate CaN with comparable or reduced activity relative to mCaM. We used molecular dynamics (MD) simulations and fluorometric assays of CaN-dependent dephosphorylation of MUF-P to probe whether calcium and protein-protein binding interactions are altered by plant CaMs relative to mCaM as a basis for differential CaN regulation. In the presence of CaN, we found that the two sCaMs' Ca2+ binding properties, such as their predicted coordination of Ca2+ and experimentally measured EC50 [Ca2+] values are comparable to mCaM. Furthermore, the binding of CaM to the CaM binding region (CaMBR) in CaN is comparable among the three CaMs, as evidenced by MD-predicted binding energies and experimentally measured EC50 [CaM] values. However, mCaM and sCaM-1 exhibited binding with a secondary region of CaN's regulatory domain that is weakened for sCaM-4. We speculate that this secondary interaction affects the turnover rate (kcat) of CaN based on our modeling of enzyme activity, which is consistent with our experimental data. Together, our data describe how plant-derived CaM variants alter CaN activity through enlisting interactions other than those directly influencing Ca2+ binding and canonical CaMBR binding, which may additionally play a role in the differential regulation of other mammalian targets.

Protein-Protein Interface Topology as a Predictor of Secondary Structure and Molecular Function Using Convolutional Deep Learning

To power the specific recognition and binding of protein partners into functional complexes, a wealth of information about the structure and function of the partners is necessarily encoded into the global shape of protein-protein interfaces and their local topological features. To identify whether this is the case, this study uses convolutional deep learning methods (typically leveraged for 2D image recognition) on 3D voxel representations of protein-protein interfaces colored by burial depth. A novel two-stage network fed with voxelizations of each interface at two distinct resolutions achieves balance between performance and computational cost. From the shape of the interfaces, the network tries to predict the presence of secondary structure motifs at the interface and the molecular function of the corresponding complex. Secondary structure and certain classes of function are found to be very well predicted, validating the hypothesis that interface shape is a conveyor of higher-level information. Interface patterns triggering the recognition of specific classes are also identified and described.

Machine Learning Driven Analysis of Large Scale Simulations Reveals Conformational Characteristics of Ubiquitin Chains

Understanding the conformational characteristics of protein complexes in solution is crucial for a deeper insight in their biological function. Molecular dynamics simulations performed on high performance computing plants and with modern simulation techniques can be used to obtain large data sets that contain conformational and thermodynamic information about biomolecular systems. While this can in principle give a detailed picture of protein-protein interactions in solution and therefore complement experimental data, it also raises the challenge of processing exceedingly large high-dimensional data sets with several million samples. Here we present a novel method for the characterization of protein-protein interactions, which combines a neural network based dimensionality reduction technique to obtain a two-dimensional representation of the conformational space with a density based clustering algorithm for state detection and a metric which assesses the (dis)similarity between different conformational spaces. This method is highly scalable and therefore makes the analysis of massive data sets computationally tractable. We demonstrate the power of this approach to large scale data analysis by characterizing highly dynamic conformational phase spaces of differently linked ubiquitin (Ub) oligomers from coarse-grained simulations. We are able to extract a protein-protein interaction model for two unlinked Ub proteins which is then used to determine how the Ub-Ub interaction pattern is altered in Ub oligomers by the introduction of a covalent linkage. We find that the Ub chain conformational ensemble depends highly on the linkage type and for some cases also on the Ub chain length. By this, we obtain insight into the conformational characteristics of different Ub chains and how this may contribute to linkage type and chain length specific recognition.

Conformational and functional characterization of artificially conjugated non-canonical ubiquitin dimers

Ubiquitylation is an eminent posttranslational modification referring to the covalent attachment of single ubiquitin molecules or polyubiquitin chains to a target protein dictating the fate of such labeled polypeptide chains. Here, we have biochemically produced artificially Lys11-, and Lys27-, and Lys63-linked ubiquitin dimers based on click-chemistry generating milligram quantities in high purity. We show that the artificial linkage used for the conjugation of two ubiquitin moieties represents a fully reliable surrogate of the natural isopeptide bond by acquiring highly resolved nuclear magnetic resonance (NMR) spectroscopic data including ligand binding studies. Extensive coarse grained and atomistic molecular dynamics (MD) simulations allow to extract structures representing the ensemble of domain-domain conformations used to verify the experimental data. Advantageously, this methodology does not require individual isotopic labeling of both ubiquitin moieties as NMR data have been acquired on the isotopically labeled proximal moiety and complementary MD simulations have been used to fully interpret the experimental data in terms of domain-domain conformation. This combined approach intertwining NMR spectroscopy with MD simulations makes it possible to describe the conformational space non-canonically Lys11-, and Lys27-linked ubiquitin dimers occupy in a solution averaged ensemble by taking atomically resolved information representing all residues in ubiquitin dimers into account.

Multi-resolution simulations of intracellular processes

This is an introduction to the Thematic Issue of Interface Focus containing papers by speakers of the Theo Murphy International Scientific Meeting on ‘Multi-resolution simulations of intracellular processes’