Electrically Detectable Photoinduced Polarization Switching in a Molecular Prussian Blue Analogue.

Light, a nondestructive and remotely controllable external stimulus, effectively triggers a variety of electron-transfer phenomena in metal complexes. One prime example includes using light in molecular cyanide-bridged [FeCo] bimetallic Prussian blue analogues, where it switches the system between the electron-transferred metastable state and the system's ground state. If this process is coupled to a ferroelectric-type phase transition, the generation and disappearance of macroscopic polarization, entirely under light control, become possible. In this research, we successfully executed a nonpolar-to-polar phase transition in a trinuclear cyanide-bridged [Fe2Co] complex crystal via directional electron transfer. Intriguingly, by exposing the crystal to the wavelength of light─785 nm─without any electric field─we can drive this ferroelectric phase transition to completely depolarize the crystal, during which a measurable electric current response can be detected. These discoveries signify an important step toward the realization of fully light-controlled ferroelectric memory devices.

Diffusion of a disordered protein on its folded ligand.

Intrinsically disordered proteins often form dynamic complexes with their ligands. Yet, the speed and amplitude of these motions are hidden in classical binding kinetics. Here, we directly measure the dynamics in an exceptionally mobile, high-affinity complex. We show that the disordered tail of the cell adhesion protein E-cadherin dynamically samples a large surface area of the protooncogene ß-catenin. Single-molecule experiments and molecular simulations resolve these motions with high resolution in space and time. Contacts break and form within hundreds of microseconds without a dissociation of the complex. The energy landscape of this complex is rugged with many small barriers (3 to 4 kBT) and reconciles specificity, high affinity, and extreme disorder. A few persistent contacts provide specificity, whereas unspecific interactions boost affinity.

Development of an FeII Complex Exhibiting Intermolecular Proton Shifting Coupled Spin Transition.

In this study, we investigated reversible intermolecular proton shifting (IPS) coupled with spin transition (ST) in a novel FeII complex. The host FeII complex and the guest carboxylic acid anion were connected by intermolecular hydrogen bonds (IHBs). We extended the intramolecular proton transfer coupled ST phenomenon to the intermolecular system. The dynamic phenomenon was confirmed by variable-temperature single-crystal X-ray diffraction, neutron crystallography, and infrared spectroscopy. The mechanism of IPS was further validated using density functional theory calculations. The discovery of IPS-coupled ST in crystalline molecular materials provides good insights into fundamental processes and promotes the design of novel multifunctional materials with tunable properties for various applications, such as optoelectronics, information storage, and molecular devices.

Fully Optical Control of Polarization Current Direction in a Cyanide-Bridged Trinuclear Complex.

As a remote and non-contact stimulus, light offers the potential for manipulating the polarization of ferroelectric materials without physical contact. However, in current research, the non-contact write-read (erase) process lacks direct observation through the stable current as output signal. To address this limitation, we investigated the photoinduced polarization switching capabilities of the cyanide-bridged compound [Fe2Co] using visible light, leading to the achievement of rewritable polarization. By subjecting [Fe2Co] crystals to alternating irradiation with 785 nm and 532 nm light, the polarization changes exhibited a distinct square wave pattern, confirming the reliability of the writing and erasing processes. Initialization involved exposing specific crystal units to 532 nm light for storing "1" or "0" information, while reading was accomplished by scanning the units with 785 nm light, resulting in brief current pulses for "1" states and no current signal for "0" states. This research unveils new possibilities for optical storage systems, paving the way for efficient and rewritable data storage and retrieval technologies, such as the next-generation memories.

Magnetoelectricity Enhanced by Electron Redistribution in a Spin Crossover [FeCo] Complex.

Molecular-based magnetoelectric materials are among the most promising materials for next-generation magnetoelectric memory devices. However, practical application of existing molecular systems has proven difficult largely because the polarization change is far lower than the practical threshold of the ME memory devices. Herein, we successfully obtained an [FeCo] dinuclear complex that exhibits a magnetic field-induced spin crossover process, resulting in a significant polarization change of 0.45 µC cm-2. Mössbauer spectroscopy and theoretical calculations suggest that the asymmetric structural change, coupled with electron redistribution, leads to the observed polarization change. Our approach provides a new strategy toward rationally enhancing the polarization change.

Macroscopic Polarization Change of Mononuclear Valence Tautomeric Cobalt Complexes Through the Use of Enantiopure Ligand.

The crystallization of a complex having electron transfer properties in a polar space group can induce the polarization switching of a crystal in a specific direction, which is attractive for the development of sensors, memory devices, and capacitors. Unfortunately, the probability of crystallization in a polar space group is usually low. Noticing that enantiopure compounds crystallize in Sohncke space groups, this paper reports a strategy for the molecular design of non-ferroelectric polarization switching crystals based on the use of intramolecular electron transfer and chirality. In addition, this paper describes the synthesis of a mononuclear valence tautomeric (VT) cobalt complex bearing an enantiopure ligand. The introduction of enantiomer enables the crystallization of the complex in the polar space group (P21 ). The polarization of the crystals along the b-axis direction is not canceled out and the VT transition is accompanied by a change in the macroscopic polarization of the polar crystal. Polarization switching via electron transfer is realized at around room temperature.

Sequence dependent phase separation of protein-polynucleotide mixtures elucidated using molecular simulations.

Ribonucleoprotein (RNP) granules are membraneless organelles (MLOs), which majorly consist of RNA and RNA-binding proteins and are formed via liquid-liquid phase separation (LLPS). Experimental studies investigating the drivers of LLPS have shown that intrinsically disordered proteins (IDPs) and nucleic acids like RNA and other polynucleotides play a key role in modulating protein phase separation. There is currently a dearth of modelling techniques which allow one to delve deeper into how polynucleotides play the role of a modulator/promoter of LLPS in cells using computational methods. Here, we present a coarse-grained polynucleotide model developed to fill this gap, which together with our recently developed HPS model for protein LLPS, allows us to capture the factors driving protein-polynucleotide phase separation. We explore the capabilities of the modelling framework with the LAF-1 RGG system which has been well studied in experiments and also with the HPS model previously. Further taking advantage of the fact that the HPS model maintains sequence specificity we explore the role of charge patterning on controlling polynucleotide incorporation into condensates. With increased charge patterning we observe formation of structured or patterned condensates which suggests the possible roles of polynucleotides in not only shifting the phase boundaries but also introducing microscopic organization in MLOs.

Polymer effects modulate binding affinities in disordered proteins.

Structural disorder is widespread in regulatory protein networks. Weak and transient interactions render disordered proteins particularly sensitive to fluctuations in solution conditions such as ion and crowder concentrations. How this sensitivity alters folding coupled binding reactions, however, has not been fully understood. Here, we demonstrate that salt jointly modulates polymer properties and binding affinities of 5 disordered proteins from a transcription factor network. A combination of single-molecule Förster resonance energy transfer experiments, polymer theory, and molecular simulations shows that all 5 proteins expand with increasing ionic strengths due to Debye-Hückel charge screening. Simultaneously, pairwise affinities between the proteins increase by an order of magnitude within physiological salt limits. A quantitative analysis shows that 50% of the affinity increase can be explained by changes in the disordered state. Disordered state properties therefore have a functional relevance even if these states are not directly involved in biological functions. Numerical solutions of coupled binding equilibria with our results show that networks of homologous disordered proteins can function surprisingly robustly in fluctuating cellular environments, despite the sensitivity of its individual proteins.

Quantitative NMR Study of Insulin-Degrading Enzyme Using Amyloid-ß and HIV-1 p6 Elucidates Its Chaperone Activity.

Insulin-degrading enzyme (IDE) hydrolyzes monomeric polypeptides, including amyloid-ß (Aß) and HIV-1 p6. It also acts as a nonproteolytic chaperone to prevent Aß polymerization. Here we compare interactions of Aß and non-amyloidogenic p6 with IDE. Although both exhibited similar proteolysis rates, the binding kinetics to an inactive IDE characterized using relaxation-based NMR were remarkably different. IDE and Aß formed a sparsely populated complex with a lifetime of milliseconds in which a short hydrophobic cleavage segment of Aß was anchored to IDE. Strikingly, a second and more stable complex was significantly populated with a subsecond lifetime owing to multiple intermolecular contacts between Aß and IDE. By selectively sequestering Aß in this nonproductive complex, IDE likely increases the critical concentration required for fibrillization. In contrast, IDE and p6 formed a transient, submillisecond complex involving a single anchoring p6 motif. Modulation of intermolecular interactions, thus, allows IDE to differentiate between non-amyloidogenic and amyloidogenic substrates.

Phosphorylation of the FUS low-complexity domain disrupts phase separation, aggregation, and toxicity.

Neuronal inclusions of aggregated RNA-binding protein fused in sarcoma (FUS) are hallmarks of ALS and frontotemporal dementia subtypes. Intriguingly, FUS's nearly uncharged, aggregation-prone, yeast prion-like, low sequence-complexity domain (LC) is known to be targeted for phosphorylation. Here we map in vitro and in-cell phosphorylation sites across FUS LC We show that both phosphorylation and phosphomimetic variants reduce its aggregation-prone/prion-like character, disrupting FUS phase separation in the presence of RNA or salt and reducing FUS propensity to aggregate. Nuclear magnetic resonance spectroscopy demonstrates the intrinsically disordered structure of FUS LC is preserved after phosphorylation; however, transient domain collapse and self-interaction are reduced by phosphomimetics. Moreover, we show that phosphomimetic FUS reduces aggregation in human and yeast cell models, and can ameliorate FUS-associated cytotoxicity. Hence, post-translational modification may be a mechanism by which cells control physiological assembly and prevent pathological protein aggregation, suggesting a potential treatment pathway amenable to pharmacologic modulation.

Relation between single-molecule properties and phase behavior of intrinsically disordered proteins.

Proteins that undergo liquid-liquid phase separation (LLPS) have been shown to play a critical role in many physiological functions through formation of condensed liquid-like assemblies that function as membraneless organelles within biological systems. To understand how different proteins may contribute differently to these assemblies and their functions, it is important to understand the molecular driving forces of phase separation and characterize their phase boundaries and material properties. Experimental studies have shown that intrinsically disordered regions of these proteins are a major driving force, as many of them undergo LLPS in isolation. Previous work on polymer solution phase behavior suggests a potential correspondence between intramolecular and intermolecular interactions that can be leveraged to discover relationships between single-molecule properties and phase boundaries. Here, we take advantage of a recently developed coarse-grained framework to calculate the θ temperature [Formula: see text], the Boyle temperature [Formula: see text], and the critical temperature [Formula: see text] for 20 diverse protein sequences, and we show that these three properties are highly correlated. We also highlight that these correlations are not specific to our model or simulation methodology by comparing between different pairwise potentials and with data from other work. We, therefore, suggest that smaller simulations or experiments to determine [Formula: see text] or [Formula: see text] can provide useful insights into the corresponding phase behavior.

Sequence determinants of protein phase behavior from a coarse-grained model.

Membraneless organelles important to intracellular compartmentalization have recently been shown to comprise assemblies of proteins which undergo liquid-liquid phase separation (LLPS). However, many proteins involved in this phase separation are at least partially disordered. The molecular mechanism and the sequence determinants of this process are challenging to determine experimentally owing to the disordered nature of the assemblies, motivating the use of theoretical and simulation methods. This work advances a computational framework for conducting simulations of LLPS with residue-level detail, and allows for the determination of phase diagrams and coexistence densities of proteins in the two phases. The model includes a short-range contact potential as well as a simplified treatment of electrostatic energy. Interaction parameters are optimized against experimentally determined radius of gyration data for multiple unfolded or intrinsically disordered proteins (IDPs). These models are applied to two systems which undergo LLPS: the low complexity domain of the RNA-binding protein FUS and the DEAD-box helicase protein LAF-1. We develop a novel simulation method to determine thermodynamic phase diagrams as a function of the total protein concentration and temperature. We show that the model is capable of capturing qualitative changes in the phase diagram due to phosphomimetic mutations of FUS and to the presence or absence of the large folded domain in LAF-1. We also explore the effects of chain-length, or multivalency, on the phase diagram, and obtain results consistent with Flory-Huggins theory for polymers. Most importantly, the methodology presented here is flexible so that it can be easily extended to other pair potentials, be used with other enhanced sampling methods, and may incorporate additional features for biological systems of interest.

Inferring properties of disordered chains from FRET transfer efficiencies.

Förster resonance energy transfer (FRET) is a powerful tool for elucidating both structural and dynamic properties of unfolded or disordered biomolecules, especially in single-molecule experiments. However, the key observables, namely, the mean transfer efficiency and fluorescence lifetimes of the donor and acceptor chromophores, are averaged over a broad distribution of donor-acceptor distances. The inferred average properties of the ensemble therefore depend on the form of the model distribution chosen to describe the distance, as has been widely recognized. In addition, while the distribution for one type of polymer model may be appropriate for a chain under a given set of physico-chemical conditions, it may not be suitable for the same chain in a different environment so that even an apparently consistent application of the same model over all conditions may distort the apparent changes in chain dimensions with variation of temperature or solution composition. Here, we present an alternative and straightforward approach to determining ensemble properties from FRET data, in which the polymer scaling exponent is allowed to vary with solution conditions. In its simplest form, it requires either the mean FRET efficiency or fluorescence lifetime information. In order to test the accuracy of the method, we have utilized both synthetic FRET data from implicit and explicit solvent simulations for 30 different protein sequences, and experimental single-molecule FRET data for an intrinsically disordered and a denatured protein. In all cases, we find that the inferred radii of gyration are within 10% of the true values, thus providing higher accuracy than simpler polymer models. In addition, the scaling exponents obtained by our procedure are in good agreement with those determined directly from the molecular ensemble. Our approach can in principle be generalized to treating other ensemble-averaged functions of intramolecular distances from experimental data.

The primary growth of laryngeal squamous cell carcinoma cells in vitro is effectively supported by paired cancer-associated fibroblasts alone.

Most primarily cultured laryngeal squamous cell carcinoma cells are difficult to propagate in vitro and have a low survival rate. However, in our previous work to establish a laryngeal squamous cell carcinoma cell line, we found that laryngeal cancer-associated fibroblasts appeared to strongly inhibit the apoptosis of primarily cultured laryngeal squamous cell carcinoma cells in vitro. In this study, we investigated whether paired laryngeal cancer-associated fibroblasts alone can effectively support the growth of primarily cultured laryngeal squamous cell carcinoma cells in vitro. In all, 29 laryngeal squamous cell carcinoma specimens were collected and primarily cultured. The laryngeal squamous cell carcinoma cells were separated from cancer-associated fibroblasts by differential trypsinization and continuously subcultured. Morphological changes of the cultured laryngeal squamous cell carcinoma cells were observed. Immunocytofluorescence was used to authenticate the identity of the cancer-associated fibroblasts and laryngeal squamous cell carcinoma cells. Flow cytometry was used to quantify the proportion of apoptotic cells. Western blot was used to detect the protein levels of caspase-3. Enzyme-linked immunosorbent assay was used to detect the levels of chemokine (C-X-C motif) ligand 12, chemokine (C-X-C motif) ligand 7, hepatocyte growth factor, and fibroblast growth factor 1 in the supernatants of the laryngeal squamous cell carcinoma and control cells. AMD3100 (a chemokine (C-X-C motif) receptor 4 antagonist) and an anti-chemokine (C-X-C motif) ligand 7 antibody were used to block the tumor-supporting capacity of cancer-associated fibroblasts. Significant apoptotic changes were detected in the morphology of laryngeal squamous cell carcinoma cells detached from cancer-associated fibroblasts. The percentage of apoptotic laryngeal squamous cell carcinoma cells and the protein levels of caspase-3 increased gradually in subsequent subcultures. In contrast, no significant differences in the proliferation capacity of laryngeal squamous cell carcinoma cells cocultured with cancer-associated fibroblasts were detected during subculturing. High level of chemokine (C-X-C motif) ligand 12 was detected in the culture supernatant of cancer-associated fibroblasts. The tumor-supporting effect of cancer-associated fibroblasts was significantly inhibited by AMD3100. Our findings demonstrate that the paired laryngeal cancer-associated fibroblasts alone are sufficient to support the primary growth of laryngeal squamous cell carcinoma cells in vitro and that the chemokine (C-X-C motif) ligand 12/chemokine (C-X-C motif) receptor 4 axis is one of the major contributors.

Probing the Action of Chemical Denaturant on an Intrinsically Disordered Protein by Simulation and Experiment.

Chemical denaturants are the most commonly used agents for unfolding proteins and are thought to act by better solvating the unfolded state. Improved solvation is expected to lead to an expansion of unfolded chains with increasing denaturant concentration, providing a sensitive probe of the denaturant action. However, experiments have so far yielded qualitatively different results concerning the effects of chemical denaturation. Studies using Förster resonance energy transfer (FRET) and other methods found an increase in radius of gyration with denaturant concentration, but most small-angle X-ray scattering (SAXS) studies found no change. This discrepancy therefore challenges our understanding of denaturation mechanism and more generally the accuracy of these experiments as applied to unfolded or disordered proteins. Here, we use all-atom molecular simulations to investigate the effect of urea and guanidinium chloride on the structure of the intrinsically disordered protein ACTR, which can be studied by experiment over a wide range of denaturant concentration. Using unbiased molecular simulations with a carefully calibrated denaturant model, we find that the protein chain indeed swells with increasing denaturant concentration. This is due to the favorable association of urea or guanidinium chloride with the backbone of all residues and with the side-chains of almost all residues, with denaturant-water transfer free energies inferred from this association in reasonable accord with experimental estimates. Interactions of the denaturants with the backbone are dominated by hydrogen bonding, while interactions with side-chains include other contributions. By computing FRET efficiencies and SAXS intensities at each denaturant concentration, we show that the simulation trajectories are in accord with both experiments on this protein, demonstrating that there is no fundamental inconsistency between the two types of experiment. Agreement with experiment also supports the picture of chemical denaturation described in our simulations, driven by weak association of denaturant with the protein. Our simulations support some assumptions needed for each experiment to accurately reflect changes in protein size, namely, that the commonly used FRET chromophores do not qualitatively alter the results and that possible effects such as preferential solvent partitioning into the interior of the chain do not interfere with the determination of radius of gyration from the SAXS experiments.

Consistent View of Polypeptide Chain Expansion in Chemical Denaturants from Multiple Experimental Methods.

There has been a long-standing controversy regarding the effect of chemical denaturants on the dimensions of unfolded and intrinsically disordered proteins: A wide range of experimental techniques suggest that polypeptide chains expand with increasing denaturant concentration, but several studies using small-angle X-ray scattering (SAXS) have reported no such increase of the radius of gyration (Rg). This inconsistency challenges our current understanding of the mechanism of chemical denaturants, which are widely employed to investigate protein folding and stability. Here, we use a combination of single-molecule Förster resonance energy transfer (FRET), SAXS, dynamic light scattering (DLS), and two-focus fluorescence correlation spectroscopy (2f-FCS) to characterize the denaturant dependence of the unfolded state of the spectrin domain R17 and the intrinsically disordered protein ACTR in two different denaturants. Standard analysis of the primary data clearly indicates an expansion of the unfolded state with increasing denaturant concentration irrespective of the protein, denaturant, or experimental method used. This is the first case in which SAXS and FRET have yielded even qualitatively consistent results regarding expansion in denaturant when applied to the same proteins. To more directly illustrate this self-consistency, we used both SAXS and FRET data in a Bayesian procedure to refine structural ensembles representative of the observed unfolded state. This analysis demonstrates that both of these experimental probes are compatible with a common ensemble of protein configurations for each denaturant concentration. Furthermore, the resulting ensembles reproduce the trend of increasing hydrodynamic radius with denaturant concentration obtained by 2f-FCS and DLS. We were thus able to reconcile the results from all four experimental techniques quantitatively, to obtain a comprehensive structural picture of denaturant-induced unfolded state expansion, and to identify the most likely sources of earlier discrepancies.

Dependence of internal friction on folding mechanism.

An outstanding challenge in protein folding is understanding the origin of "internal friction" in folding dynamics, experimentally identified from the dependence of folding rates on solvent viscosity. A possible origin suggested by simulation is the crossing of local torsion barriers. However, it was unclear why internal friction varied from protein to protein or for different folding barriers of the same protein. Using all-atom simulations with variable solvent viscosity, in conjunction with transition-path sampling to obtain reaction rates and analysis via Markov state models, we are able to determine the internal friction in the folding of several peptides and miniproteins. In agreement with experiment, we find that the folding events with greatest internal friction are those that mainly involve helix formation, while hairpin formation exhibits little or no evidence of friction. Via a careful analysis of folding transition paths, we show that internal friction arises when torsion angle changes are an important part of the folding mechanism near the folding free energy barrier. These results suggest an explanation for the variation of internal friction effects from protein to protein and across the energy landscape of the same protein.

Multiscale approach to the determination of the photoactive yellow protein signaling state ensemble.

The nature of the optical cycle of photoactive yellow protein (PYP) makes its elucidation challenging for both experiment and theory. The long transition times render conventional simulation methods ineffective, and yet the short signaling-state lifetime makes experimental data difficult to obtain and interpret. Here, through an innovative combination of computational methods, a prediction and analysis of the biological signaling state of PYP is presented. Coarse-grained modeling and locally scaled diffusion map are first used to obtain a rough bird's-eye view of the free energy landscape of photo-activated PYP. Then all-atom reconstruction, followed by an enhanced sampling scheme; diffusion map-directed-molecular dynamics are used to focus in on the signaling-state region of configuration space and obtain an ensemble of signaling state structures. To the best of our knowledge, this is the first time an all-atom reconstruction from a coarse grained model has been performed in a relatively unexplored region of molecular configuration space. We compare our signaling state prediction with previous computational and more recent experimental results, and the comparison is favorable, which validates the method presented. This approach provides additional insight to understand the PYP photo cycle, and can be applied to other systems for which more direct methods are impractical.

A comparative analysis of clustering algorithms: O2 migration in truncated hemoglobin I from transition networks.

The ligand migration network for O2-diffusion in truncated Hemoglobin N is analyzed based on three different clustering schemes. For coordinate-based clustering, the conventional k-means and the kinetics-based Markov Clustering (MCL) methods are employed, whereas the locally scaled diffusion map (LSDMap) method is a collective-variable-based approach. It is found that all three methods agree well in their geometrical definition of the most important docking site, and all experimentally known docking sites are recovered by all three methods. Also, for most of the states, their population coincides quite favourably, whereas the kinetics of and between the states differs. One of the major differences between k-means and MCL clustering on the one hand and LSDMap on the other is that the latter finds one large primary cluster containing the Xe1a, IS1, and ENT states. This is related to the fact that the motion within the state occurs on similar time scales, whereas structurally the state is found to be quite diverse. In agreement with previous explicit atomistic simulations, the Xe3 pocket is found to be a highly dynamical site which points to its potential role as a hub in the network. This is also highlighted in the fact that LSDMap cannot identify this state. First passage time distributions from MCL clusterings using a one- (ligand-position) and two-dimensional (ligand-position and protein-structure) descriptor suggest that ligand- and protein-motions are coupled. The benefits and drawbacks of the three methods are discussed in a comparative fashion and highlight that depending on the questions at hand the best-performing method for a particular data set may differ.

Discovering mountain passes via torchlight: methods for the definition of reaction coordinates and pathways in complex macromolecular reactions.

The long-timescale dynamics of macromolecular systems can be oftentimes viewed as a reaction connecting metastable states of the system. In the past decade, various approaches have been developed to discover the collective motions associated with this dynamics. The corresponding collective variables are used in many applications, e.g., to understand the reaction mechanism, to quantify the system's free energy landscape, to enhance the sampling of the reaction path, and to determine the reaction rate. In this review we focus on a number of key developments in this field, providing an overview of several methods along with their relative regimes of applicability.