Many molecular chaperones exist as oligomeric complexes in their functional states, yet the physical determinants underlying such self-assembly behavior, as well as the role of oligomerization in the activity of molecular chaperones in inhibiting protein aggregation, have proven to be difficult to define. Here, we demonstrate direct measurements under native conditions of the changes in the average oligomer populations of a chaperone system as a function of concentration and time and thus determine the thermodynamic and kinetic parameters governing the self-assembly process. We access this self-assembly behavior in real time under native-like conditions by monitoring the changes in the micrometer-scale diffusion of the different complexes in time and space using a microfluidic platform. Using this approach, we find that the oligomerization mechanism of the Hsp70 subdomain occurs in a cooperative manner and involves structural constraints that limit the size of the species formed beyond the limits imposed by mass balance. These results illustrate the ability of microfluidic methods to probe polydisperse protein self-assembly in real time in solution and to shed light on the nature and dynamics of oligomerization processes.

CRM1/XPO1 (CRM1) is a nuclear export chaperone that mediates the export of proteins essential to growth regulation and tumor suppression. Its overexpression in tumors was found to be associated with poor prognosis. Selective inhibitors of nuclear export are in phase I and II clinical trials for several tumor types. Our aim was to investigate CRM1 expression in pancreatic adenocarcinoma (PAC) and its relationship to survivin expression and the proliferative activity. Sections of tissue microarray containing 76 formalin fixed and paraffin embedded PAC were stained by immunohistochemistry (IHC) for CRM1, survivin, and Cyclin A. Expression levels of CRM1 and survivin and the proliferative activity, the S-phase fraction (SPF) in tumor cells, were determined using a quantitative digital image analysis solution (OTMIAS). Sixty-six of the 76 (86%) PAC showed positive staining for CRM1, and 10 (14%) were completely negative. The mean CRM1 expression levels ranged from 0.3 to 53 units and the median from 0.3 to 45 units. There was significant positive correlation between the mean and median expression levels of CRM1 in tumor cells and the mean and median levels of survivin (p

Despite converging epidemiological evidence for the inverse relationship of regular caffeine consumption and risk of developing Parkinson's disease (PD) with animal studies demonstrating protective effect of caffeine in various neurotoxin models of PD, whether caffeine can protect against mutant α-synuclein (α-Syn) A53T-induced neurotoxicity in intact animals has not been examined. Here, we determined the effect of chronic caffeine treatment using the α-Syn fibril model of PD by intra-striatal injection of preformed A53T α-Syn fibrils. We demonstrated that chronic caffeine treatment blunted a cascade of pathological events leading to α-synucleinopathy, including pSer129α-Syn-rich aggregates, apoptotic neuronal cell death, microglia, and astroglia reactivation. Importantly, chronic caffeine treatment did not affect autophagy processes in the normal striatum, but selectively reversed α-Syn-induced defects in macroautophagy (by enhancing microtubule-associated protein 1 light chain 3, and reducing the receptor protein sequestosome 1, SQSTM1/p62) and chaperone-mediated autophagy (CMA, by enhancing LAMP2A). These findings support that caffeine-a strongly protective environment factor as suggested by epidemiological evidence-may represent a novel pharmacological therapy for PD by targeting autophagy pathway.

Potassium-sensitive hypokalaemic and normokalaemic periodic paralysis are inherited skeletal muscle diseases characterized by episodes of flaccid muscle weakness1,2. They are caused by single mutations in positively charged residues ('gating charges') in the S4 transmembrane segment of the voltage sensor of the voltage-gated sodium channel Nav1.4 or the calcium channel Cav1.11,2. Mutations of the outermost gating charges (R1 and R2) cause hypokalaemic periodic paralysis1,2 by creating a pathogenic gating pore in the voltage sensor through which cations leak in the resting state3,4. Mutations of the third gating charge (R3) cause normokalaemic periodic paralysis 5 owing to cation leak in both activated and inactivated states 6 . Here we present high-resolution structures of the model bacterial sodium channel NavAb with the analogous gating-charge mutations7,8, which have similar functional effects as in the human channels. The R2G and R3G mutations have no effect on the backbone structures of the voltage sensor, but they create an aqueous cavity near the hydrophobic constriction site that controls gating charge movement through the voltage sensor. The R3G mutation extends the extracellular aqueous cleft through the entire length of the activated voltage sensor, creating an aqueous path through the membrane. Conversely, molecular modelling shows that the R2G mutation creates a continuous aqueous path through the membrane only in the resting state. Crystal structures of NavAb(R2G) in complex with guanidinium define a potential drug target site. Molecular dynamics simulations illustrate the mechanism of Na+ permeation through the mutant gating pore in concert with conformational fluctuations of the gating charge R4. Our results reveal pathogenic mechanisms of periodic paralysis at the atomic level and suggest designs of drugs that may prevent ionic leak and provide symptomatic relief from hypokalaemic and normokalaemic periodic paralysis.

In order to elucidate features of the denatured state ensembles that exist in equilibrium with the native state under physiological conditions, we performed 1.4-μs molecular dynamics (MD) simulations at 400 K and 450 K using the monomer subunits of three CutA1 mutants from Escherichia coli: an SH-free mutant (Ec0SH) with denaturation temperature (Td) = 85.6 °C, a hydrophobic mutant (Ec0VV) with Td = 113.3 °C, and an ionic mutant (Ec0VV_6) with Td = 136.8 °C. The occupancy of salt bridges by the six substituted charged residues in Ec0VV_6 was 140.1% at 300 K and 89.5% at 450 K, indicating that even in the denatured state, salt bridge occupancy was high, approximately 60% of that at 300 K. From these results, we can infer that proteins from hyperthermophiles with a high ratio of charged residues are stabilized by a decrease in conformational entropy due to ion-ion interactions in the denatured state. The mechanism must be comparable to the stabilization conferred by disulfide bonds within a protein. This suggests that introduction of charged residues, to promote formation of salt bridges in the denatured state, would be a simple way to rationally design stability-enhanced mutants.

Various experimental and computational techniques have been employed over the past decade to provide structural and thermodynamic insights into G Protein-Coupled Receptor (GPCR) dimerization. Here, we use multiple microsecond-long, coarse-grained, biased and unbiased molecular dynamics simulations (a total of ~4 milliseconds) combined with multi-ensemble Markov state models to elucidate the kinetics of homodimerization of a prototypic GPCR, the µ-opioid receptor (MOR), embedded in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/cholesterol lipid bilayer. Analysis of these computations identifies kinetically distinct macrostates comprising several different short-lived dimeric configurations of either inactive or activated MOR. Calculated kinetic rates and fractions of dimers at different MOR concentrations suggest a negligible population of MOR homodimers at physiological concentrations, which is supported by acceptor photobleaching fluorescence resonance energy transfer (FRET) experiments. This study provides a rigorous, quantitative explanation for some conflicting experimental data on GPCR oligomerization.

Reactive oxygen and nitrogen species released by cold physical plasma are being proposed as effectors in various clinical conditions connected to inflammatory processes. As these plasmas can be tailored in a wide range, models to compare and control their biochemical footprint are desired to infer on the molecular mechanisms underlying the observed effects and to enable the discrimination between different plasma sources. Here, an improved model to trace short-lived reactive species is presented. Using FTIR, high-resolution mass spectrometry, and molecular dynamics computational simulation, covalent modifications of cysteine treated with different plasmas were deciphered and the respective product pattern used to generate a fingerprint of each plasma source. Such, our experimental model allows a fast and reliable grading of the chemical potential of plasmas used for medical purposes. Major reaction products were identified to be cysteine sulfonic acid, cystine, and cysteine fragments. Less-abundant products, such as oxidized cystine derivatives or S-nitrosylated cysteines, were unique to different plasma sources or operating conditions. The data collected point at hydroxyl radicals, atomic O, and singlet oxygen as major contributing species that enable an impact on cellular thiol groups when applying cold plasma in vitro or in vivo.

Previous researches have revealed the importance of shear and the orientation dependence in the structural transition of iron. In this work, we introduce a series of shear deformations by adjusting the strain ratio between the longitudinal ([001]) and transversal ([010] and [100]) directions, and then investigate this structural transition under different anisotropic compressions with molecular dynamics simulations. It is found that the shear deformation can lower the transition pressure notably, and even change the nucleation structure and morphology. Under 1D-dominated compression (along (001) direction), there only appears hcp nucleation with a few fcc stacking faults. For other cases, more equivalent planes will be activated and fcc structure begins to nucleate. Under 2D-dominated compression (along (010) and (001) directions), the fcc mass fraction is already over the hcp phase. At last, we compare the variations of shear stress and potential energy for different phases, and present the sliding mechanism under typical anisotropic compressions.

Unc-93 homolog B1 (UNC93B1) is a key regulator of nucleic acid (NA)-sensing Toll-like receptors (TLRs). Loss of NA-sensing TLR responses in UNC93B1-deficient patients facilitates Herpes simplex virus type 1 (HSV-1) encephalitis. UNC93B1 is thought to guide NA-sensing TLRs from the endoplasmic reticulum (ER) to their respective endosomal signaling compartments and to guide the flagellin receptor TLR5 to the cell surface, raising the question of how UNC93B1 mediates differential TLR trafficking. Here, we report that UNC93B1 regulates a step upstream of the differential TLR trafficking process. We discovered that UNC93B1 deficiency resulted in near-complete loss of TLR3 and TLR7 proteins in primary splenic mouse dendritic cells and macrophages, showing that UNC93B1 is critical for maintaining TLR expression. Notably, expression of an ER-retained UNC93B1 version was sufficient to stabilize TLRs and largely restore endosomal TLR trafficking and activity. These data are critical for an understanding of how UNC93B1 can regulate the function of a broad subset of TLRs.

We employ site-specific, linear and nonlinear infrared spectroscopic techniques as well as fluorescence spectroscopy and molecular dynamics simulations to investigate the binding interactions of urea and three of its derivatives, methylurea, 1,3-dimethylurea, and tetramethylurea, with protein aromatic and polar side chains. We find that (1) urea methylation leads to preferential interactions between the cosolvent molecules and aromatic side chains with an affinity that increases with the number of methyl groups; (2) interactions with tetramethylurea cause significant dehydration of aromatic side chains and the effect is most pronounced for tryptophan; and (3) while neither urea nor tetramethylurea shows preferential accumulation around a polar side chain, the number of hydrogen-bond donors around this side chain is significantly decreased in the presence of tetramethylurea. Taken together, our findings suggest that these urea derivatives, especially tetramethylurea, can effectively disrupt hydrophobic interactions in proteins. Additionally, tetramethylurea can promote intramolecular hydrogen-bond formation and hence induce α-helix folding in peptides, as observed.

We used nearest-neighbor searches in chemical space to improve the activity of the antimicrobial peptide dendrimer (AMPD) G3KL and identified dendrimer T7, which has an expanded activity range against Gram-negative pathogenic bacteria including Klebsiellae pneumoniae, increased serum stability, and promising activity in an in vivo infection model against a multidrug-resistant strain of Acinetobacter baumannii. Imaging, spectroscopic studies, and a structural model from molecular dynamics simulations suggest that T7 acts through membrane disruption. These experiments provide the first example of using virtual screening in the field of dendrimers and show that dendrimer size does not limit the activity of AMPDs.

Misfolding mutations in the DNA-binding domain of p53 alter its conformation, affecting the efficiency with which it binds to chromatin to regulate target gene expression and cell cycle checkpoint functions in many cancers, including glioblastoma. Small molecule drugs that recover misfolded p53 structure and function may improve chemotherapy by activating p53-mediated senescence. We constructed and optimized a split Renilla luciferase (RLUC) complementation molecular biosensor (NRLUC-p53-CRLUC) to determine small molecule-meditated folding changes in p53 protein. After initial evaluation of the biosensor in three different cells lines, we engineered endogenously p53P98L mutant (i.e. not affecting the DNA-binding domain) Ln229 glioblastoma cells, to express the biosensor containing one of four different p53 proteins: p53wt, p53Y220C, p53G245S and p53R282W. We evaluated the consequent phenotypic changes in these four variant cells as well as the parental cells after exposure to PhiKan083 and SCH529074, drugs previously reported to activate mutant p53 folding. Specifically, we measured induced RLUC complementation and consequent therapeutic response. Upon stable transduction with the p53 biosensors, we demonstrated that these originally p53P98L Ln229 cells had acquired p53 cellular phenotypes representative of each p53 protein expressed within the biosensor fusion protein. In these engineered variants we found a differential drug response when treated with doxorubicin and temozolomide, either independently or in combination with PhiKan083 or SCH529074. We thus developed a molecular imaging complementation biosensor that mimics endogenous p53 function for use in future applications to screen novel or repurposed drugs that counter the effects of misfolding mutations responsible for oncogenic structural changes in p53.

The dynamics of liquid water feature a variety of time scales, ranging from extremely fast ballistic-like thermal motion, to slower molecular diffusion and hydrogen-bond rearrangements. Here, we utilize coherent X-ray pulses to investigate the sub-100 fs equilibrium dynamics of water from ambient conditions down to supercooled temperatures. This novel approach utilizes the inherent capability of X-ray speckle visibility spectroscopy to measure equilibrium intermolecular dynamics with lengthscale selectivity, by measuring oxygen motion in momentum space. The observed decay of the speckle contrast at the first diffraction peak, which reflects tetrahedral coordination, is attributed to motion on a molecular scale within the first 120 fs. Through comparison with molecular dynamics simulations, we conclude that the slowing down upon cooling from 328 K down to 253 K is not due to simple thermal ballistic-like motion, but that cage effects play an important role even on timescales over 25 fs due to hydrogen-bonding.

Cleavage of collagen by collagenases such as matrix metalloproteinase 1 (MMP-1) is a key step in development, tissue remodeling, and tumor proliferation. The abundant heterotrimeric type I collagen composed of two α1(I) chains and one α2(I) chain is efficiently cleaved by MMP-1 at a unique site in the triple helix, a process which may be initiated by local unfolding within the peptide chains. Atypical homotrimers of the α1(I) chain, found in embryonic and cancer tissues, are very resistant to MMP cleavage. To investigate MMP-1 cleavage, recombinant homotrimers were constructed with sequences from the MMP cleavage regions of human collagen chains inserted into a host bacterial collagen protein system. All triple-helical constructs were cleaved by MMP-1, with α2(I) homotrimers cleaved efficiently at a rate similar to that seen for α1(II) and α1(III) homotrimers, while α1(I) homotrimers were cleaved at a much slower rate. The introduction of destabilizing Gly to Ser mutations within the human collagenase susceptible region of the α2(I) chain did not interfere with MMP-1 cleavage. Molecular dynamics simulations indicated a greater degree of transient hydrogen bond breaking in α2(I) homotrimers compared with α1(I) homotrimers at the MMP-1 cleavage site, and showed an extensive disruption of hydrogen bonding in the presence of a Gly to Ser mutation, consistent with chymotrypsin digestion results. This study indicates that α2(I) homotrimers are susceptible to MMP-1, proves that the presence of an α1(I) chain is not a requirement for α2(I) cleavage, and supports the importance of local unfolding of α2(I) in collagenase cleavage.

Hsp70 chaperones bind to various protein substrates for folding, trafficking, and degradation. Considerable structural information is available about how prokaryotic Hsp70 (DnaK) binds substrates, but less is known about mammalian Hsp70s, of which there are 13 isoforms encoded in the human genome. Here, we report the interaction between the human Hsp70 isoform heat shock cognate 71 KDa protein (Hsc70 or HSPA8) and peptides derived from the microtubule-associated protein tau, which is linked to Alzheimer's disease. For structural studies, we used an Hsc70 construct (called BETA) comprising the substrate-binding domain, but lacking the lid. Importantly, we found that truncating the lid does not significantly impair Hsc70's chaperone activity or allostery in vitro. Using NMR, we show that BETA is partially dynamically disordered in the absence of substrate and that binding of the tau sequence GKVQIINKKG (with a KD = 500 nM) causes dramatic rigidification of BETA. Nuclear Overhauser effect distance measurements revealed that tau binds to the canonical substrate-binding cleft, similar to the binding observed with DnaK. To further develop BETA as a tool for studying Hsc70 interactions, we also measured BETA binding in NMR and fluorescent competition assays to peptides derived from huntingtin, insulin, a second tau-recognition sequence, and a KFERQ-like sequence linked to chaperone-mediated autophagy. We found that the insulin C-peptide binds BETA with high affinity (KD < 100 nM), whereas the others do not (KD > 100 μM). Together, our findings reveal several similarities and differences in how prokaryotic and mammalian Hsp70 isoforms interact with different substrate peptides.

The histone chaperone complex facilitates chromatin transcription (FACT) plays important roles in DNA repair, replication, and transcription. In the formation of this complex, structure-specific recognition protein-1 (SSRP1) heterodimerizes with suppressor of Ty 16 (SPT16). SSRP1 also has SPT16-independent functions, but how SSRP1 functions alone remains elusive. Here, using analytical ultracentrifugation (AUC) and small-angle X-ray scattering (SAXS) techniques, we characterized human SSRP1 and that from the amoeba Dictyostelium discoideum and show that both orthologs form an elongated homodimer in solution. We found that substitutions in the SSRP1 pleckstrin homology domain known to bind SPT16 also disrupt SSRP1 homodimerization. Moreover, AUC and SAXS analyses revealed that SSRP1 homodimerization and heterodimerization with SPT16 (resulting in FACT) involve the same SSRP1 surface, namely, the PH2 region, and that the FACT complex contains only one molecule of SSRP1. These observations suggest that SSRP1 homo- and heterodimerization might be mutually exclusive. Moreover, isothermal titration calorimetry analyses disclosed that SSRP1 binds both histones H2A-H2B and H3-H4 and that disruption of SSRP1 homodimerization decreases its histone-binding affinity. Together, our results provide evidence for regulation of SSRP1 by homodimerization and suggest a potential role for homodimerization in facilitating SPT16-independent functions of SSRP1.

Castration-resistant prostate cancer (CRPC) is characterized by reactivation of androgen receptor (AR) signaling in part by elevated expression of AR splice variants (ARv) including ARv7, a constitutively active, ligand binding domain (LBD)-deficient variant whose expression has been correlated with therapeutic resistance and poor prognosis. In a screen to identify small molecule dual inhibitors of both androgen-dependent and androgen-independent AR gene signatures, we identified the chalcone C86. Binding studies using purified proteins and CRPC cell lysates revealed C86 to interact with heat shock protein 40 (Hsp40). Pulldown studies using biotinylated-C86 found Hsp40 present in a multi-protein complex with full-length (FL-) AR, ARv7 and Hsp70 in CRPC cells. Treatment of CRPC cells with C86 or the allosteric heat shock protein 70 (Hsp70) inhibitor JG98 resulted in rapid protein destabilization of both FL-AR and ARv, including ARv7, concomitant with reduced FL-AR- and ARv7-mediated transcriptional activity. The glucocorticoid receptor (GR), whose elevated expression in a subset of CRPC also leads to androgen-independent AR target gene transcription, was also destabilized by inhibition of Hsp40 or Hsp70. In vivo, Hsp40 or Hsp70 inhibition demonstrated single agent and combinatorial activity in a 22Rv1 CRPC xenograft model. These data reveal that, in addition to recognized roles of Hsp40 and Hsp70 in FL-AR LBD remodeling, ARv lacking the LBD remain dependent on molecular chaperones for stability and function. Our findings highlight the feasibility and potential benefit of targeting the Hsp40/Hsp70 chaperone axis to treat prostate cancer that has become resistant to standard anti-androgen therapy.

BACKGROUND: Influenza viruses are undergoing continuous and rapid evolution. The fatal influenza A/H7N9 has drawn attention since the first wave of infections in March 2013, and raised more grave concerns with its increased potential to spread among humans. Experimental studies have revealed several host and virulence markers, indicating differential host binding preferences which can help estimate the potential of causing a pandemic. Here we systematically investigate the sequence pattern and structural characteristics of novel influenza A/H7N9 using computational approaches. RESULTS: The sequence analysis highlighted mutations in protein functional domains of influenza viruses. Molecular docking and molecular dynamics simulation revealed that the hemagglutinin (HA) of A/Taiwan/1/2017(H7N9) strain enhanced the binding with both avian and human receptor analogs, compared with the previous A/Shanghai/02/2013(H7N9) strain. The Molecular Mechanics - Poisson Boltzmann Surface Area (MM-PBSA) calculation revealed the change of residue-ligand interaction energy and detected the residues with conspicuous binding preference. CONCLUSION: The results are novel and specific to the emerging influenza A/Taiwan/1/2017(H7N9) strain compared with A/Shanghai/02/2013(H7N9). Its enhanced ability to bind human receptor analogs, which are abundant in the human upper respiratory tract, may be responsible for the recent outbreak. Residues showing binding preference were detected, which could facilitate monitoring the circulating influenza viruses.

The hypothalamic-pituitary-adrenal (HPA) axis is the major neuroendocrine axis regulating homeostasis in mammals. Glucocorticoid hormones are rapidly synthesized and secreted from the adrenal gland in response to stress. In addition, under basal conditions glucocorticoids are released rhythmically with both a circadian and an ultradian (pulsatile) pattern. These rhythms are important not only for normal function of glucocorticoid target organs, but also for the HPA axis responses to stress. Several studies have shown that disruption of glucocorticoid rhythms is associated with disease both in humans and in rodents. In this review, we will discuss our knowledge of the negative feedback mechanisms that regulate basal ultradian synthesis and secretion of glucocorticoids, including the role of glucocorticoid and mineralocorticoid receptors and their chaperone protein FKBP51. Moreover, in light of recent findings, we will also discuss the importance of intra-adrenal glucocorticoid receptor signaling in regulating glucocorticoid synthesis.

The development of effective centroid potentials (ECPs) is explored with both the constrained-centroid and quasi-adiabatic force matching using liquid water as a test system. A trajectory integrated with the ECP is free of statistical noises that would be introduced when the centroid potential is approximated on the fly with a finite number of beads. With the reduced cost of ECP, challenging experimental properties can be studied in the spirit of centroid molecular dynamics. The experimental number density of H2O is 0.38% higher than that of D2O. With the ECP, the H2O number density is predicted to be 0.42% higher, when the dispersion term is not refit. After correction of finite size effects, the diffusion constant of H2O is found to be 21% higher than that of D2O, which is in good agreement with the 29.9% higher diffusivity for H2O observed experimentally. Although the ECP is also able to capture the redshifts of both the OH and OD stretching modes in liquid water, there are a number of properties that a classical simulation with the ECP will not be able to recover. For example, the heat capacities of H2O and D2O are predicted to be almost identical and higher than the experimental values. Such a failure is simply a result of not properly treating quantized vibrational energy levels when the trajectory is propagated with classical mechanics. Several limitations of the ECP based approach without bead population reconstruction are discussed.

Understanding the microscopic mechanism of water molecular translational diffusion is a challenging topic in both physics and chemistry. Here, we report an investigation on the interplay between the translational diffusion and the large-amplitude angular jumps of water molecules in bulk water using molecular dynamics simulations. We found that large-amplitude angular jumps are tightly coupled to the translational diffusions. Particularly, we revealed that concurrent rotational jumps of spatially neighboring water molecules induce inter-basin translational jumps, which contributes to the fast component of the water translational diffusion. Consequently, the translational diffusion shows positional heterogeneity; i.e., the neighbors of the water molecules with inter-basin translational jumps have larger probability to diffuse by inter-basin translational jumps. Our control simulations showed that a model water molecule with moderate hydrogen bond strength can diffuse much faster than a simple Lennard-Jones particle in bulk water due to the capability of disturbing the hydrogen bond network of the surrounding water molecules. Our results added to the understanding of the microscopic picture of the water translational diffusion and demonstrated the unique features of water diffusion arising from their hydrogen bond network structure compared with those of the simple liquids.

To develop prediction methods of three-phase equilibrium (coexistence) conditions of methane hydrate by molecular simulations, we examined the use of NVT (isometric-isothermal) molecular dynamics (MD) simulations. NVT MD simulations of coexisting solid hydrate, liquid water, and vapor methane phases were performed at four different temperatures, namely, 285, 290, 295, and 300 K. NVT simulations do not require complex pressure control schemes in multi-phase systems, and the growth or dissociation of the hydrate phase can lead to significant pressure changes in the approach toward equilibrium conditions. We found that the calculated equilibrium pressures tended to be higher than those reported by previous NPT (isobaric-isothermal) simulation studies using the same water model. The deviations of equilibrium conditions from previous simulation studies are mainly attributable to the employed calculation methods of pressure and Lennard-Jones interactions. We monitored the pressure in the methane phase, far from the interfaces with other phases, and confirmed that it was higher than the total pressure of the system calculated by previous studies. This fact clearly highlights the difficulties associated with the pressure calculation and control for multi-phase systems. The treatment of Lennard-Jones interactions without tail corrections in MD simulations also contributes to the overestimation of equilibrium pressure. Although improvements are still required to obtain accurate equilibrium conditions, NVT MD simulations exhibit potential for the prediction of equilibrium conditions of multi-phase systems.

Inspired by the protein folding and DNA packing, we have systematically studied the thermodynamic and kinetic behaviors of single semiflexible homopolymers by Langevin dynamics simulations. In line with experiments, a rich variety of folding products, such as rod-like bundles, hairpins, toroids, and a mixture of them, are observed in the complete diagram of states. Moreover, knotted structures with a significant population are found in a certain range of bending stiffness in thermal equilibrium. As the solvent quality becomes poorer, the population of the intermediate occurring in the folding process increases, which leads to a severe chevron rollover for the folding arm. However, the population of the intermediates in the unfolding process is very low, insufficient to induce unfolding arm rollover. The total types of folding pathways from the coil state to the toroidal state for a semiflexible polymer chain remain unchanged by varying the solvent quality or temperature, whereas the kinetic partitioning into different folding events can be tuned significantly. In the process of knotting, three types of mechanisms, namely, plugging, slipknotting, and sliding, are discovered. Along the folding evolution, a semiflexible homopolymer chain can knot at any stage of folding upon leaving the extended coil state, and the probability to find a knot increases with chain compactness. In addition, we find rich types of knotted topologies during the folding of a semiflexible homopolymer chain. This study should be helpful in gaining insight into the general principles of biopolymer folding.

Chaperone-mediated autophagy (CMA), a form of selective autophagy, maintains cellular proteostasis in response to diverse stress conditions. Whether and how endoplasmic reticulum (ER) stress triggers CMA remains elusive. In our recent study, we demonstrate that various types of ER stress activate the CMA pathway via an EIF2AK3/PERK-MAP2K4/MKK4-MAPK14/p38-dependent manner. We term this process ERICA for ER stress-induced chaperone-mediated autophagy. This pathway is activated in response to stress associated with Parkinson disease and is required for the viability of the SNc dopaminergic neurons in an animal model of Parkinson disease.