Compartmentalization of the SUMO/RNF4 pathway by SLX4 drives DNA repair.

SLX4, disabled in the Fanconi anemia group P, is a scaffolding protein that coordinates the action of structure-specific endonucleases and other proteins involved in the replication-coupled repair of DNA interstrand cross-links. Here, we show that SLX4 dimerization and SUMO-SIM interactions drive the assembly of SLX4 membraneless compartments in the nucleus called condensates. Super-resolution microscopy reveals that SLX4 forms chromatin-bound clusters of nanocondensates. We report that SLX4 compartmentalizes the SUMO-RNF4 signaling pathway. SENP6 and RNF4 regulate the assembly and disassembly of SLX4 condensates, respectively. SLX4 condensation per se triggers the selective modification of proteins by SUMO and ubiquitin. Specifically, SLX4 condensation induces ubiquitylation and chromatin extraction of topoisomerase 1 DNA-protein cross-links. SLX4 condensation also induces the nucleolytic degradation of newly replicated DNA. We propose that the compartmentalization of proteins by SLX4 through site-specific interactions ensures the spatiotemporal control of protein modifications and nucleolytic reactions during DNA repair.

Molecular dissection of amyloid disaggregation by human HSP70.

The deposition of highly ordered fibrillar-type aggregates into inclusion bodies is a hallmark of neurodegenerative diseases such as Parkinson's disease. The high stability of such amyloid fibril aggregates makes them challenging substrates for the cellular protein quality-control machinery1,2. However, the human HSP70 chaperone and its co-chaperones DNAJB1 and HSP110 can dissolve preformed fibrils of the Parkinson's disease-linked presynaptic protein α-synuclein in vitro3,4. The underlying mechanisms of this unique activity remain poorly understood. Here we use biochemical tools and nuclear magnetic resonance spectroscopy to determine the crucial steps of the disaggregation process of amyloid fibrils. We find that DNAJB1 specifically recognizes the oligomeric form of α-synuclein via multivalent interactions, and selectively targets HSP70 to fibrils. HSP70 and DNAJB1 interact with the fibril through exposed, flexible amino and carboxy termini of α-synuclein rather than the amyloid core itself. The synergistic action of DNAJB1 and HSP110 strongly accelerates disaggregation by facilitating the loading of several HSP70 molecules in a densely packed arrangement at the fibril surface, which is ideal for the generation of 'entropic pulling' forces. The cooperation of DNAJB1 and HSP110 in amyloid disaggregation goes beyond the classical substrate targeting and recycling functions that are attributed to these HSP70 co-chaperones and constitutes an active and essential contribution to the remodelling of the amyloid substrate. These mechanistic insights into the essential prerequisites for amyloid disaggregation may provide a basis for new therapeutic interventions in neurodegeneration.

Data-driven large-scale genomic analysis reveals an intricate phylogenetic and functional landscape in J-domain proteins.

The 70-kD heat shock protein (Hsp70) chaperone system is a central hub of the proteostasis network that helps maintain protein homeostasis in all organisms. The recruitment of Hsp70 to perform different and specific cellular functions is regulated by the J-domain protein (JDP) co-chaperone family carrying the small namesake J-domain, required to interact and drive the ATPase cycle of Hsp70s. Besides the J-domain, prokaryotic and eukaryotic JDPs display a staggering diversity in domain architecture, function, and cellular localization. Very little is known about the overall JDP family, despite their essential role in cellular proteostasis, development, and its link to a broad range of human diseases. In this work, we leverage the exponentially increasing number of JDP gene sequences identified across all kingdoms owing to the advancements in sequencing technology and provide a broad overview of the JDP repertoire. Using an automated classification scheme based on artificial neural networks (ANNs), we demonstrate that the sequences of J-domains carry sufficient discriminatory information to reliably recover the phylogeny, localization, and domain composition of the corresponding full-length JDP. By harnessing the interpretability of the ANNs, we find that many of the discriminatory sequence positions match residues that form the interaction interface between the J-domain and Hsp70. This reveals that key residues within the J-domains have coevolved with their obligatory Hsp70 partners to build chaperone circuits for specific functions in cells.

DNA-segment-capture model for loop extrusion by structural maintenance of chromosome (SMC) protein complexes.

Cells possess remarkable control of the folding and entanglement topology of long and flexible chromosomal DNA molecules. It is thought that structural maintenance of chromosome (SMC) protein complexes play a crucial role in this, by organizing long DNAs into series of loops. Experimental data suggest that SMC complexes are able to translocate on DNA, as well as pull out lengths of DNA via a 'loop extrusion' process. We describe a Brownian loop-capture-ratchet model for translocation and loop extrusion based on known structural, catalytic, and DNA-binding properties of the Bacillus subtilis SMC complex. Our model provides an example of a new class of molecular motor where large conformational fluctuations of the motor 'track'-in this case DNA-are involved in the basic translocation process. Quantitative analysis of our model leads to a series of predictions for the motor properties of SMC complexes, most strikingly a strong dependence of SMC translocation velocity and step size on tension in the DNA track that it is moving along, with 'stalling' occuring at subpiconewton tensions. We discuss how the same mechanism might be used by structurally related SMC complexes (Escherichia coli MukBEF and eukaryote condensin, cohesin and SMC5/6) to organize genomic DNA.

Chaperones convert the energy from ATP into the nonequilibrium stabilization of native proteins.

During and after protein translation, molecular chaperones require ATP hydrolysis to favor the native folding of their substrates and, under stress, to avoid aggregation and revert misfolding. Why do some chaperones need ATP, and what are the consequences of the energy contributed by the ATPase cycle? Here, we used biochemical assays and physical modeling to show that the bacterial chaperones GroEL (Hsp60) and DnaK (Hsp70) both use part of the energy from ATP hydrolysis to restore the native state of their substrates, even under denaturing conditions in which the native state is thermodynamically unstable. Consistently with thermodynamics, upon exhaustion of ATP, the metastable native chaperone products spontaneously revert to their equilibrium non-native states. In the presence of ATPase chaperones, some proteins may thus behave as open ATP-driven, nonequilibrium systems whose fate is only partially determined by equilibrium thermodynamics.

Insights into the activation mechanism of human estrogen-related receptor γ by environmental endocrine disruptors.

The estrogen-related receptor γ (ERRγ, NR3B3) is a constitutively active nuclear receptor which has been proposed to act as a mediator of the low-dose effects of a number of environmental endocrine-disrupting chemicals (EDCs) such as the xenoestrogen bisphenol-A (BPA). To better characterize the ability of exogenous compounds to bind and activate ERRγ, we used a combination of cell-based, biochemical, structural and computational approaches. A purposely created stable cell line allowed for the determination of the EC50s for over 30 environmental ERRγ ligands, including previously unknown ones. Interestingly, affinity constants (Kds) of the most potent compounds measured by isothermal titration calorimetry were in the 50-500 nM range, in agreement with their receptor activation potencies. Crystallographic analysis of the interaction between the ERRγ ligand-binding domain (LBD) and compounds of the bisphenol, alkylphenol and naphthol families revealed a partially shared binding mode and minimal alterations of the receptor conformation upon ligand binding. Further biophysical characterizations coupled to molecular dynamics simulations suggested a mechanism through which ERRγ ligands would exhibit their agonistic properties by preserving the transcriptionally active form of the receptor while rigidifying some loop regions with associated functions. This unique mechanism contrasts with the classical one involving a ligand-induced repositioning and stabilization of the C-terminal activation helix H12.

Prion protein ß2-α2 loop conformational landscape.

In transmissible spongiform encephalopathies (TSEs), which are lethal neurodegenerative diseases that affect humans and a wide range of other mammalian species, the normal "cellular" prion protein ([Formula: see text]) is transformed into amyloid aggregates representing the "scrapie form" of the protein ([Formula: see text]). Continued research on this system is of keen interest, since new information on the physiological function of [Formula: see text] in healthy organisms is emerging, as well as new data on the mechanism of the transformation of [Formula: see text] to [Formula: see text] In this paper we used two different approaches: a combination of the well-tempered ensemble (WTE) and parallel tempering (PT) schemes and metadynamics (MetaD) to characterize the conformational free-energy surface of [Formula: see text] The focus of the data analysis was on an 11-residue polypeptide segment in mouse [Formula: see text](121-231) that includes the [Formula: see text]2-[Formula: see text]2 loop of residues 167-170, for which a correlation between structure and susceptibility to prion disease has previously been described. This study includes wild-type mouse [Formula: see text] and a variant with the single-residue replacement Y169A. The resulting detailed conformational landscapes complement in an integrative manner the available experimental data on [Formula: see text], providing quantitative insights into the nature of the structural transition-related function of the [Formula: see text]2-[Formula: see text]2 loop.

Catch bond drives stator mechanosensitivity in the bacterial flagellar motor.

The bacterial flagellar motor (BFM) is the rotary motor that rotates each bacterial flagellum, powering the swimming and swarming of many motile bacteria. The torque is provided by stator units, ion motive force-powered ion channels known to assemble and disassemble dynamically in the BFM. This turnover is mechanosensitive, with the number of engaged units dependent on the viscous load experienced by the motor through the flagellum. However, the molecular mechanism driving BFM mechanosensitivity is unknown. Here, we directly measure the kinetics of arrival and departure of the stator units in individual motors via analysis of high-resolution recordings of motor speed, while dynamically varying the load on the motor via external magnetic torque. The kinetic rates obtained, robust with respect to the details of the applied adsorption model, indicate that the lifetime of an assembled stator unit increases when a higher force is applied to its anchoring point in the cell wall. This provides strong evidence that a catch bond (a bond strengthened instead of weakened by force) drives mechanosensitivity of the flagellar motor complex. These results add the BFM to a short, but growing, list of systems demonstrating catch bonds, suggesting that this "molecular strategy" is a widespread mechanism to sense and respond to mechanical stress. We propose that force-enhanced stator adhesion allows the cell to adapt to a heterogeneous environmental viscosity and may ultimately play a role in surface-sensing during swarming and biofilm formation.

Coevolutionary Analysis of Protein Subfamilies by Sequence Reweighting.

Extracting structural information from sequence co-variation has become a common computational biology practice in the recent years, mainly due to the availability of large sequence alignments of protein families. However, identifying features that are specific to sub-classes and not shared by all members of the family using sequence-based approaches has remained an elusive problem. We here present a coevolutionary-based method to differentially analyze subfamily specific structural features by a continuous sequence reweighting (SR) approach. We introduce the underlying principles and test its predictive capabilities on the Response Regulator family, whose subfamilies have been previously shown to display distinct, specific homo-dimerization patterns. Our results show that this reweighting scheme is effective in assigning structural features known a priori to subfamilies, even when sequence data is relatively scarce. Furthermore, sequence reweighting allows assessing if individual structural contacts pertain to specific subfamilies and it thus paves the way for the identification specificity-determining contacts from sequence variation data.

Aromatic Ring Dynamics, Thermal Activation, and Transient Conformations of a 468 kDa Enzyme by Specific 1H-13C Labeling and Fast Magic-Angle Spinning NMR.

Aromatic residues are located at structurally important sites of many proteins. Probing their interactions and dynamics can provide important functional insight but is challenging in large proteins. Here, we introduce approaches to characterize the dynamics of phenylalanine residues using 1H-detected fast magic-angle spinning (MAS) NMR combined with a tailored isotope-labeling scheme. Our approach yields isolated two-spin systems that are ideally suited for artifact-free dynamics measurements, and allows probing motions effectively without molecular weight limitations. The application to the TET2 enzyme assembly of ∼0.5 MDa size, the currently largest protein assigned by MAS NMR, provides insights into motions occurring on a wide range of time scales (picoseconds to milliseconds). We quantitatively probe ring-flip motions and show the temperature dependence by MAS NMR measurements down to 100 K. Interestingly, favorable line widths are observed down to 100 K, with potential implications for DNP NMR. Furthermore, we report the first 13C R1ρ MAS NMR relaxation-dispersion measurements and detect structural excursions occurring on a microsecond time scale in the entry pore to the catalytic chamber and at a trimer interface that was proposed as the exit pore. We show that the labeling scheme with deuteration at ca. 50 kHz MAS provides superior resolution compared to 100 kHz MAS experiments with protonated, uniformly 13C-labeled samples.

Lipid tempering simulation of model biological membranes on parallel platforms.

In this report we have tested a parallel implementation for the simulation of lipid bilayers at the atomistic level, based on a generalized ensemble protocol where only the torsional degrees of freedom of the alkyl chains of the lipids are heated. The results in terms of configurational sampling enhancement have been compared with a conventional simulation produced with a widespread molecular dynamics code. Results show that the proposed thermodynamic-based multiple trajectories parallel protocol for membrane simulations allows for an efficient use of CPU resources with respect to the conventional single trajectory, providing accurate results for area and volume per lipid, membrane thickness, undulation spectra and boosting significantly diffusion and mixing in lipid bilayers due to the sampling enhancement of gauche/trans ratios of the alkyl chain dihedral angles.

Single-molecule spectroscopy reveals chaperone-mediated expansion of substrate protein.

Molecular chaperones are an essential part of the machinery that avoids protein aggregation and misfolding in vivo. However, understanding the molecular basis of how chaperones prevent such undesirable interactions requires the conformational changes within substrate proteins to be probed during chaperone action. Here we use single-molecule fluorescence spectroscopy to investigate how the DnaJ-DnaK chaperone system alters the conformational distribution of the denatured substrate protein rhodanese. We find that in a first step the ATP-independent binding of DnaJ to denatured rhodanese results in a compact denatured ensemble of the substrate protein. The following ATP-dependent binding of multiple DnaK molecules, however, leads to a surprisingly large expansion of denatured rhodanese. Molecular simulations indicate that hard-core repulsion between the multiple DnaK molecules provides the underlying mechanism for disrupting even strong interactions within the substrate protein and preparing it for processing by downstream chaperone systems.

Large-Scale Conformational Transitions and Dimerization Are Encoded in the Amino-Acid Sequences of Hsp70 Chaperones.

Hsp70s are a class of ubiquitous and highly conserved molecular chaperones playing a central role in the regulation of proteostasis in the cell. Hsp70s assist a myriad of cellular processes by binding unfolded or misfolded substrates during a complex biochemical cycle involving large-scale structural rearrangements. Here we show that an analysis of coevolution at the residue level fully captures the characteristic large-scale conformational transitions of this protein family, and predicts an evolutionary conserved-and thus functional-homo-dimeric arrangement. Furthermore, we highlight that the features encoding the Hsp70 dimer are more conserved in bacterial than in eukaryotic sequences, suggesting that the known Hsp70/Hsp110 hetero-dimer is a eukaryotic specialization built on a pre-existing template.

Communication: Role of explicit water models in the helix folding/unfolding processes.

In the last years, it has become evident that computer simulations can assume a relevant role in modelling protein dynamical motions for their ability to provide a full atomistic image of the processes under investigation. The ability of the current protein force-fields in reproducing the correct thermodynamics and kinetics systems behaviour is thus an essential ingredient to improve our understanding of many relevant biological functionalities. In this work, employing the last developments of the metadynamics framework, we compare the ability of state-of-the-art all-atom empirical functions and water models to consistently reproduce the folding and unfolding of a helix turn motif in a model peptide. This theoretical study puts in evidence that the choice of the water models can influence the thermodynamic and the kinetics of the system under investigation, and for this reason cannot be considered trivial.

The allosteric communication pathways in KIX domain of CBP.

Allosteric regulation plays an important role in a myriad of biomacromolecular processes. Specifically, in a protein, the process of allostery refers to the transmission of a local perturbation, such as ligand binding, to a distant site. Decades after the discovery of this phenomenon, models built on static images of proteins are being reconsidered with the knowledge that protein dynamics plays an important role in its function. Molecular dynamics simulations are a valuable tool for studying complex biomolecular systems, providing an atomistic description of their structure and dynamics. Unfortunately, their predictive power has been limited by the complexity of the biomolecule free-energy surface and by the length of the allosteric timescale (in the order of milliseconds). In this work, we are able to probe the origins of the allosteric changes that transcription factor mixed lineage leukemia (MLL) causes to the interactions of KIX domain of CREB-binding protein (CBP) with phosphorylated kinase inducible domain (pKID), by combing all-atom molecular dynamics with enhanced sampling methods recently developed in our group. We discuss our results in relation to previous NMR studies. We also develop a general simulations protocol to study allosteric phenomena and many other biological processes that occur in the micro/milliseconds timescale.

Free-energy landscape of protein oligomerization from atomistic simulations.

In the realm of protein-protein interactions, the assembly process of homooligomers plays a fundamental role because the majority of proteins fall into this category. A comprehensive understanding of this multistep process requires the characterization of the driving molecular interactions and the transient intermediate species. The latter are often short-lived and thus remain elusive to most experimental investigations. Molecular simulations provide a unique tool to shed light onto these complex processes complementing experimental data. Here we combine advanced sampling techniques, such as metadynamics and parallel tempering, to characterize the oligomerization landscape of fibritin foldon domain. This system is an evolutionarily optimized trimerization motif that represents an ideal model for experimental and computational mechanistic studies. Our results are fully consistent with previous experimental nuclear magnetic resonance and kinetic data, but they provide a unique insight into fibritin foldon assembly. In particular, our simulations unveil the role of nonspecific interactions and suggest that an interplay between thermodynamic bias toward native structure and residual conformational disorder may provide a kinetic advantage.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems.

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Structural Insights into the Activation of Human Aryl Hydrocarbon Receptor by the Environmental Contaminant Benzo[a]pyrene and Structurally Related Compounds.

The aryl hydrocarbon receptor (AHR) is a ligand-dependent transcription factor belonging to the bHLH/PAS protein family and responding to hundreds of natural and chemical substances. It is primarily involved in the defense against chemical insults and bacterial infections or in the adaptive immune response, but also in the development of pathological conditions ranging from inflammatory to neoplastic disorders. Despite its prominent roles in many (patho)physiological processes, the lack of high-resolution structural data has precluded for thirty years an in-depth understanding of the structural mechanisms underlying ligand-binding specificity, promiscuity and activation of AHR. We recently reported a cryogenic electron microscopy (cryo-EM) structure of human AHR bound to the natural ligand indirubin, the chaperone Hsp90 and the co-chaperone XAP2 that provided the first experimental visualization of its ligand-binding PAS-B domain. Here, we report a 2.75 Å resolution structure of the AHR complex bound to the environmental pollutant benzo[a]pyrene (B[a]P). The structure substantiates the existence of a bipartite PAS-B ligand-binding pocket with a geometrically constrained primary binding site controlling ligand binding specificity and affinity, and a secondary binding site contributing to the binding promiscuity of AHR. We also report a docking study of B[a]P congeners that validates the B[a]P-bound PAS-B structure as a suitable model for accurate computational ligand binding assessment. Finally, comparison of our agonist-bound complex with the recently reported structures of mouse and fruit fly AHR PAS-B in different activation states suggests a ligand-induced loop conformational change potentially involved in the regulation of AHR function.

24-h continuous non-invasive multiparameter home monitoring of vitals in patients with Rett syndrome by an innovative wearable technology: evidence of an overlooked chronic fatigue status.

Background: Sleep is disturbed in Rett syndrome (RTT), a rare and progressive neurodevelopmental disorder primarily affecting female patients (prevalence 7.1/100,000 female patients) linked to pathogenic variations in the X-linked methyl-CpG-binding protein 2 (MECP2) gene. Autonomic nervous system dysfunction with a predominance of the sympathetic nervous system (SNS) over the parasympathetic nervous system (PSNS) is reported in RTT, along with exercise fatigue and increased sudden death risk. The aim of the present study was to test the feasibility of a continuous 24 h non-invasive home monitoring of the biological vitals (biovitals) by an innovative wearable sensor device in pediatric and adolescent/adult RTT patients. Methods: A total of 10 female patients (mean age 18.3 ± 9.4 years, range 4.7-35.5 years) with typical RTT and MECP2 pathogenic variations were enrolled. Clinical severity was assessed by validated scales. Heart rate (HR), respiratory rate (RR), and skin temperature (SkT) were monitored by the YouCare Wearable Medical Device (Accyourate Group SpA, L'Aquila, Italy). The average percentage of maximum HR (HRmax%) was calculated. Heart rate variability (HRV) was expressed by consolidated time-domain and frequency-domain parameters. The HR/LF (low frequency) ratio, indicating SNS activation under dynamic exercise, was calculated. Simultaneous continuous measurement of indoor air quality variables was performed and the patients' contributions to the surrounding water vapor partial pressure [PH2O (pt)] and carbon dioxide [PCO2 (pt)] were indirectly estimated. Results: Of the 6,559.79 h of biovital recordings, 5051.03 h (77%) were valid for data interpretation. Sleep and wake hours were 9.0 ± 1.1 h and 14.9 ± 1.1 h, respectively. HRmax % [median: 71.86% (interquartile range 61.03-82%)] and HR/LF [median: 3.75 (interquartile range 3.19-5.05)] were elevated, independent from the wake-sleep cycle. The majority of HRV time- and frequency-domain parameters were significantly higher in the pediatric patients (p ≤ 0.031). The HRV HR/LF ratio was associated with phenotype severity, disease progression, clinical sleep disorder, subclinical hypoxia, and electroencephalographic observations of multifocal epileptic activity and general background slowing. Conclusion: Our findings indicate the feasibility of a continuous 24-h non-invasive home monitoring of biovital parameters in RTT. Moreover, for the first time, HRmax% and the HR/LF ratio were identified as potential objective markers of fatigue, illness severity, and disease progression.