Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments.

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Modeling the instantaneous normal mode spectra of liquids as that of unstable elastic media.

We study the instantaneous normal mode (INM) spectrum of a simulated soft-sphere liquid at different equilibrium temperatures T We find that the spectrum of eigenvalues [Formula: see text] has a sharp maximum near (but not at) [Formula: see text] and decreases monotonically with [Formula: see text] on both the stable and unstable sides of the spectrum. The spectral shape strongly depends on temperature. It is rather asymmetric at low temperatures (close to the dynamical critical temperature) and becomes symmetric at high temperatures. To explain these findings we present a mean-field theory for [Formula: see text], which is based on a heterogeneous elasticity model, in which the local shear moduli exhibit spatial fluctuations, including negative values. We find good agreement between the simulation data and the model calculations, done with the help of the self-consistent Born approximation (SCBA), when we take the variance of the fluctuations to be proportional to the temperature T More importantly, we find an empirical correlation of the positions of the maxima of [Formula: see text] with the low-frequency exponent of the density of the vibrational modes of the glasses obtained by quenching to [Formula: see text] from the temperature T We discuss the present findings in connection to the liquid to glass transformation and its precursor phenomena.

Computational evidences of a misfolding event in an aggregation-prone light chain preceding the formation of the non-native pathogenic dimer.

Antibody light chain amyloidosis is a disorder in which protein aggregates, mainly composed of immunoglobulin light chains, deposit in diverse tissues impairing the correct functioning of organs. Interestingly, due to the high susceptibility of antibodies to mutations, AL amyloidosis appears to be strongly patient-specific. Indeed, every patient will display their own mutations that will make the proteins involved prone to aggregation thus hindering the study of this disease on a wide scale. In this framework, determining the molecular mechanisms that drive the aggregation could pave the way to the development of patient-specific therapeutics. Here, we focus on a particular patient-derived light chain, which has been experimentally characterized. We investigated the early phases of the aggregation pathway through extensive full-atom molecular dynamics simulations, highlighting a structural rearrangement and the exposure of two hydrophobic regions in the aggregation-prone species. Next, we moved to consider the pathological dimerization process through docking and molecular dynamics simulations, proposing a dimeric structure as a candidate pathological first assembly. Overall, our results shed light on the first phases of the aggregation pathway for a light chain at an atomic level detail, offering new structural insights into the corresponding aggregation process.

Dynamical changes of SARS-CoV-2 spike variants in the highly immunogenic regions impact the viral antibodies escaping.

The prolonged circulation of the SARS-CoV-2 virus resulted in the emergence of several viral variants, with different spreading features. Moreover, the increased number of recovered and/or vaccinated people introduced a selective pressure toward variants able to evade the immune system, developed against the former viral versions. This process results in reinfections. Aiming to study the latter process, we first collected a large structural dataset of antibodies in complex with the original version of SARS-CoV-2 Spike protein. We characterized the peculiarities of such antibodies population with respect to a control dataset of antibody-protein complexes, highlighting some statistically significant differences between these two sets of antibodies. Thus, moving our attention to the Spike side of the complexes, we identify the Spike region most prone to interaction with antibodies, describing in detail also the energetic mechanisms used by antibodies to recognize different epitopes. In this framework, fast protocols able to assess the effect of novel mutations on the cohort of developed antibodies would help establish the impact of the variants on the population. Performing a molecular dynamics simulation of the trimeric form of the SARS-CoV-2 Spike protein for the wild type and two variants of concern, that is, the Delta and Omicron variants, we described the physicochemical features and the conformational changes experienced locally by the variants with respect to the original version. Hence, combining the dynamical information with the structural study on the antibody-spike dataset, we quantitatively explain why the Omicron variant has a higher capability of escaping the immune system than the Delta variant, due to the higher conformational variability of the most immunogenic regions. Overall, our results shed light on the molecular mechanism behind the different responses the SARS-CoV-2 variants display against the immune response induced by either vaccines or previous infections. Moreover, our analysis proposes an approach that can be easily extended to both other SARS-CoV-2 variants or different molecular systems.

Thermometer: a webserver to predict protein thermal stability.

MOTIVATION: Thermal properties of proteins are of great importance for a number of theoretical and practical implications. Predicting the thermal stability of a protein is a difficult and still scarcely addressed task. RESULTS: Here, we introduce Thermometer, a webserver to assess the thermal stability of a protein using structural information. Thermometer is implemented as a publicly available, user-friendly interface. AVAILABILITY AND IMPLEMENTATION: Our server can be found at the following link (all major browser supported): http://service.tartaglialab.com/new_submission/thermometer_file. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Reference-less wavefront shaping in a Hopfield-like rough intensity landscape.

This study introduces a new digital-micromirror based binary-phase wavefront shaping technique, which allows the measurement of the full coupling matrix of a disordered medium without a reference and enables to focusing transmitted light. The coupling matrix takes on a bi-dyadic structure, similar to a Hopfield memory matrix containing two memory patterns. Sequential wavefront optimization in this configuration often stalls due to a rough intensity landscape, resulting in a non-optimal state. To overcome this issue, we propose the Complete Couplings Mapping method, which consistently reaches the theoretically expected maximum intensity.

Physics-informed deep neural network for image denoising.

Image enhancement deep neural networks (DNN) can improve signal to noise ratio or resolution of optically collected visual information. The literature reports a variety of approaches with varying effectiveness. All these algorithms rely on arbitrary data (the pixels' count-rate) normalization, making their performance strngly affected by dataset or user-specific data pre-manipulation. We developed a DNN algorithm capable to enhance images signal-to-noise surpassing previous algorithms. Our model stems from the nature of the photon detection process which is characterized by an inherently Poissonian statistics. Our algorithm is thus driven by distance between probability functions instead than relying on the sole count-rate, producing high performance results especially in high-dynamic-range images. Moreover, it does not require any arbitrary image renormalization other than the transformation of the camera's count-rate into photon-number.

Brillouin microscopy: an emerging tool for mechanobiology.

The role and importance of mechanical properties of cells and tissues in cellular function, development and disease has widely been acknowledged, however standard techniques currently used to assess them exhibit intrinsic limitations. Recently, Brillouin microscopy, a type of optical elastography, has emerged as a non-destructive, label- and contact-free method that can probe the viscoelastic properties of biological samples with diffraction-limited resolution in 3D. This led to increased attention amongst the biological and medical research communities, but it also sparked debates about the interpretation and relevance of the measured physical quantities. Here, we review this emerging technology by describing the underlying biophysical principles and discussing the interpretation of Brillouin spectra arising from heterogeneous biological matter. We further elaborate on the technique's limitations, as well as its potential for gaining insights in biology, in order to guide interested researchers from various fields.

Binding site identification of G protein-coupled receptors through a 3D Zernike polynomials-based method: application to C. elegans olfactory receptors.

Studying the binding processes of G protein-coupled receptors (GPCRs) proteins is of particular interest both to better understand the molecular mechanisms that regulate the signaling between the extracellular and intracellular environment and for drug design purposes. In this study, we propose a new computational approach for the identification of the binding site for a specific ligand on a GPCR. The method is based on the Zernike polynomials and performs the ligand-GPCR association through a shape complementarity analysis of the local molecular surfaces. The method is parameter-free and it can distinguish, working on hundreds of experimentally GPCR-ligand complexes, binding pockets from randomly sampled regions on the receptor surface, obtaining an Area Under ROC curve of 0.77. Given its importance both as a model organism and in terms of applications, we thus investigated the olfactory receptors of the C. elegans, building a list of associations between 21 GPCRs belonging to its olfactory neurons and a set of possible ligands. Thus, we can not only carry out rapid and efficient screenings of drugs proposed for GPCRs, key targets in many pathologies, but also we laid the groundwork for computational mutagenesis processes, aimed at increasing or decreasing the binding affinity between ligands and receptors.

Supervised perceptron learning vs unsupervised Hebbian unlearning: Approaching optimal memory retrieval in Hopfield-like networks.

The Hebbian unlearning algorithm, i.e., an unsupervised local procedure used to improve the retrieval properties in Hopfield-like neural networks, is numerically compared to a supervised algorithm to train a linear symmetric perceptron. We analyze the stability of the stored memories: basins of attraction obtained by the Hebbian unlearning technique are found to be comparable in size to those obtained in the symmetric perceptron, while the two algorithms are found to converge in the same region of Gardner's space of interactions, having followed similar learning paths. A geometric interpretation of Hebbian unlearning is proposed to explain its optimal performances. Because the Hopfield model is also a prototypical model of the disordered magnetic system, it might be possible to translate our results to other models of interest for memory storage in materials.