Uncovering the molecular interactions underlying MBD2 and MBD3 phase separation

Chromatin organization controls DNA's accessibility to regulatory factors to influence gene expression. Heterochromatin, or transcriptionally silent chromatin enriched in methylated DNA and methylated histone tails, self-assembles through multivalent interactions with its associated proteins into a condensed, but dynamic state. Liquid-liquid phase separation (LLPS) of key heterochromatin regulators, such as heterochromatin protein 1 (HP1), plays an essential role in heterochromatin assembly and function. Methyl-CpG-binding protein 2 (MeCP2), the most studied member of the methyl-CpG-binding domain (MBD) family of proteins, has been recently shown to undergo LLPS in the absence and presence of methylated DNA. These studies provide a new mechanistic framework for understanding the role of methylated DNA and its readers in heterochromatin formation. However, the details of the molecular interactions by which other MBD family members undergo LLPS to mediate genome organization and transcriptional regulation are not fully understood. Here, we focus on two MBD proteins, MBD2 and MBD3, that have distinct but interdependent roles in gene regulation. Using an integrated computational and experimental approach, we uncover the homotypic and heterotypic interactions governing MBD2 and MBD3 phase separation and DNA's influence on this process. We show that despite sharing the highest sequence identity and structural homology among all the MBD protein family members, MBD2 and MBD3 exhibit differing residue patterns resulting in distinct phase separation mechanisms. Understanding the molecular underpinnings of MBD protein condensation offers insights into the higher-order, LLPS-mediated organization of heterochromatin.

Interplay between charge distribution and DNA in shaping HP1 paralog phase separation and localization

The heterochromatin protein 1 (HP1) family is a crucial component of heterochromatin with diverse functions in gene regulation, cell cycle control, and cell differentiation. In humans, there are three paralogs, HP1α, HP1β, and HP1γ, which exhibit remarkable similarities in their domain architecture and sequence properties. Nevertheless, these paralogs display distinct behaviors in liquid-liquid phase separation (LLPS), a process linked to heterochromatin formation. Here, we employ a coarse-grained simulation framework to uncover the sequence features responsible for the observed differences in LLPS. We highlight the significance of the net charge and charge patterning along the sequence in governing paralog LLPS propensities. We also show that both highly conserved folded and less-conserved disordered domains contribute to the observed differences. Furthermore, we explore the potential co-localization of different HP1 paralogs in multicomponent assemblies and the impact of DNA on this process. Importantly, our study reveals that DNA can significantly reshape the stability of a minimal condensate formed by HP1 paralogs due to competitive interactions of HP1α with HP1β and HP1γ versus DNA. In conclusion, our work highlights the physicochemical nature of interactions that govern the distinct phase-separation behaviors of HP1 paralogs and provides a molecular framework for understanding their role in chromatin organization.

Sequence-dependent material properties of biomolecular condensates and their relation to dilute phase conformations

Material properties of phase-separated biomolecular condensates, enriched with disordered proteins, dictate many cellular functions. Contrary to the progress made in understanding the sequence-dependent phase separation of proteins, little is known about the sequence determinants of condensate material properties. Using the hydropathy scale and Martini models, we computationally decipher these relationships for charge-rich disordered protein condensates. Our computations yield dynamical, rheological, and interfacial properties of condensates that are quantitatively comparable with experimentally characterized condensates. Interestingly, we find that the material properties of model and natural proteins respond similarly to charge segregation, despite different sequence compositions. Molecular interactions within the condensates closely resemble those within the single-chain ensembles. Consequently, the material properties strongly correlate with molecular contact dynamics and single-chain structural properties. We demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Coarse-Grained Models to Study Protein-DNA Interactions and Liquid-Liquid Phase Separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define nonbonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), which is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond time scale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micrometer-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

A complex network of interdomain interactions underlies the conformational ensemble of monomeric TDP-43 and modulates its phase behavior

TAR DNA-binding protein 43 (TDP-43) is a multidomain protein involved in the regulation of RNA metabolism, and its aggregates have been observed in neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Numerous studies indicate TDP-43 can undergo liquid-liquid phase separation (LLPS) in vitro and is a component of biological condensates. Homo-oligomerization via the folded N-terminal domain (aa:1-77) and the conserved helical region (aa:319-341) of the disordered, C-terminal domain is found to be an important driver of TDP-43 phase separation. However, a comprehensive molecular view of TDP-43 phase separation, particularly regarding the nature of heterodomain interactions, is lacking due to the challenges associated with its stability and purification. Here, we utilize all-atom and coarse-grained (CG) molecular dynamics (MD) simulations to uncover the network of interdomain interactions implicated in TDP-43 phase separation. All-atom simulations uncovered the presence of transient, interdomain interactions involving flexible linkers, RNA-recognition motif (RRM) domains and a charged segment of disordered C-terminal domain (CTD). CG simulations indicate these inter-domain interactions which affect the conformational landscape of TDP-43 in the dilute phase are also prevalent in the condensed phase. Finally, sequence and surface charge distribution analysis coupled with all-atom simulations (at high salt) confirmed that the transient interdomain contacts are predominantly electrostatic in nature. Overall, our findings from multiscale simulations lead to a greater appreciation of the complex interaction network underlying the structural landscape and phase separation of TDP-43.

Sequence-Dependent Material Properties of Biomolecular Codensates and their Relation to Dilute Phase Conformations

Material properties of phase-separated biomolecular assemblies, enriched with disordered proteins, dictate their ability to participate in many cellular functions. Despite the significant effort dedicated to understanding how the sequence of the disordered protein drives its phase separation to form condensates, little is known about the sequence determinants of condensate material properties. Here, we computationally decipher these relationships for charged disordered proteins using model sequences comprised of glutamic acid and lysine residues as well as naturally occurring sequences of LAF1's RGG domain and DDX4's N-terminal domain. We do so by delineating how the arrangement of oppositely charged residues within these sequences influences the dynamical, rheological, and interfacial properties of the condensed phase through equilibrium and non-equilibrium molecular simulations using the hydropathy scale and Martini models. Our computations yield material properties that are quantitatively comparable with experimentally characterized condensate systems. Interestingly, we find that the material properties of both the model and natural proteins respond similarly to the segregation of charges, despite their very different sequence compositions. Condensates of the highly charge-segregated sequences exhibit slower dynamics than the uniformly charge-patterned sequences, because of their comparatively long-lived molecular contacts between oppositely charged residues. Surprisingly, the molecular interactions within the condensate are highly similar to those within a single-chain for all sequences. Consequently, the condensate material properties of charged disordered proteins are strongly correlated with their dense phase contact dynamics and their single-chain structural properties. Our findings demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Sequence-Dependent Conformational Transitions of Disordered Proteins During Condensation

Intrinsically disordered proteins (IDPs) can form biomolecular condensates through phase separation. It is recognized that the conformation of IDPs in the dense and dilute phases as well as at the interfaces of condensates can critically impact the resulting properties associated with their functionality. However, a comprehensive understanding of the conformational transitions of IDPs during condensation remains elusive. In this study, we employ a coarse-grained polyampholyte model, comprising an equal number of oppositely charged residues-glutamic acid and lysine-whereby conformations and phase behavior can be readily tuned by altering the protein sequence. By manipulating the sequence patterns from perfectly alternating to block-like, we obtain chains with ideal-like conformations to semi-compact structures in the dilute phase, while in the dense phase, the chain conformation is approximately that of an ideal chain, irrespective of the protein sequence. By performing simulations at different concentrations, we find that the chains assemble from the dilute phase through small oligomeric clusters to the dense phase, accompanied by a gradual swelling of the individual chains. We further demonstrate that these findings are applicable to several naturally occurring proteins involved in the formation of biological condensates. Concurrently, we delve deeper into the chain conformations within the condensate, revealing that chains at the interface show a strong sequence dependence, but remain more collapsed than those in the bulk-like dense phase. This study addresses critical gaps in our knowledge of IDP conformations within condensates as a function of protein sequence.

Interplay between charge distribution and DNA in shaping HP1 paralog phase separation and localization

The heterochromatin protein 1 (HP1) family is a crucial component of heterochromatin with diverse functions in gene regulation, cell cycle control, and cell differentiation. In humans, there are three paralogs, HP1α, HP1β, and HP1γ, which exhibit remarkable similarities in their domain architecture and sequence properties. Nevertheless, these paralogs display distinct behaviors in liquid-liquid phase separation (LLPS), a process linked to heterochromatin formation. Here, we employ a coarse-grained simulation framework to uncover the sequence features responsible for the observed differences in LLPS. We highlight the significance of the net charge and charge patterning along the sequence in governing paralog LLPS propensities. We also show that both highly conserved folded and less-conserved disordered domains contribute to the observed differences. Furthermore, we explore the potential co-localization of different HP1 paralogs in multicomponent assemblies and the impact of DNA on this process. Importantly, our study reveals that DNA can significantly reshape the stability of a minimal condensate formed by HP1 paralogs due to competitive interactions of HP1α with HP1β and HP1γ versus DNA. In conclusion, our work highlights the physicochemical nature of interactions that govern the distinct phase-separation behaviors of HP1 paralogs and provides a molecular framework for understanding their role in chromatin organization.

How a disordered linker in the Polycomb protein Polyhomeotic tunes phase separation and oligomerization

The Polycomb Group (PcG) complex PRC1 represses transcription, forms condensates in cells, and modifies chromatin architecture. These processes are connected through the essential, polymerizing Sterile Alpha Motif (SAM) present in the PRC1 subunit Polyhomeotic (Ph). In vitro, Ph SAM drives formation of short oligomers and phase separation with DNA or chromatin in the context of a Ph truncation ("mini-Ph"). Oligomer length is controlled by the long disordered linker (L) that connects the SAM to the rest of Ph--replacing Drosophila PhL with the evolutionarily diverged human PHC3L strongly increases oligomerization. How the linker controls SAM polymerization, and how polymerization and the linker affect condensate formation are not know. We analyzed PhL and PHC3L using biochemical assays and molecular dynamics (MD) simulations. PHC3L promotes mini-Ph phase separation and makes it relatively independent of DNA. In MD simulations, basic amino acids in PHC3L form contacts with acidic amino acids in the SAM. Engineering the SAM to make analogous charge-based contacts with PhL increased polymerization and phase separation, partially recapitulating the effects of the PHC3L. Ph to PHC3 linker swaps and SAM surface mutations alter Ph condensate formation in cells, and Ph function in Drosophila imaginal discs. Thus, SAM-driven phase separation and polymerization are conserved between flies and mammals, but the underlying mechanisms have diverged through changes to the disordered linker.

A coarse-grained DNA model to study protein-DNA interactions and liquid-liquid phase separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define non-bonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), that is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond timescale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micron-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

Role of Strong Localized vs Weak Distributed Interactions in Disordered Protein Phase Separation

Interaction strength and localization are critical parameters controlling the single-chain and condensed-state properties of intrinsically disordered proteins (IDPs). Here, we decipher these relationships using coarse-grained heteropolymers comprised of hydrophobic (H) and polar (P) monomers as model IDPs. We systematically vary the fraction of P monomers XP and employ two distinct particle-based models that include either strong localized attractions between only H-H pairs (HP model) or weak distributed attractions between both H-H and H-P pairs (HP+ model). To compare different sequences and models, we first carefully tune the attraction strength for all sequences to match the single-chain radius of gyration. Interestingly, we find that this procedure produces similar conformational ensembles, nonbonded potential energies, and chain-level dynamics for single chains of almost all sequences in both models, with some deviations for the HP model at large XP. However, we observe a surprisingly rich phase behavior for the sequences in both models that deviates from the expectation that similarity at the single-chain level will translate to a similar phase-separation propensity. Coexistence between dilute and dense phases is only observed up to a model-dependent XP, despite the presence of favorable interchain interactions, which we quantify using the second virial coefficient. Instead, the limited number of attractive sites (H monomers) leads to the self-assembly of finite-sized clusters of different sizes depending on XP. Our findings strongly suggest that models with distributed interactions favor the formation of liquid-like condensates over a much larger range of sequence compositions compared to models with localized interactions.

Determining interchromophore effects for energy transport in molecular networks using machine-learning algorithms

Nature uses chromophore networks, with highly optimized structural and energetic characteristics, to perform important chemical functions. Due to its modularity, predictable aggregation characteristics, and established synthetic protocols, structural DNA nanotechnology is a promising medium for arranging chromophore networks with analogous structural and energetic controls. However, this high level of control creates a greater need to know how to optimize the systems precisely. This study uses the system's modularity to produce variations of a coupled 14-Site chromophore network. It uses machine-learning algorithms and spectroscopy measurements to reveal the energy-transport roles of these Sites, paying particular attention to the cooperative and inhibitive effects they impose on each other for transport across the network. The physical significance of these patterns is contextualized, using molecular dynamics simulations and energy-transport modeling. This analysis yields insights about how energy transfers across the Donor-Relay and Relay-Acceptor interfaces, as well as the energy-transport pathways through the homogeneous Relay segment. Overall, this report establishes an approach that uses machine-learning methods to understand, in fine detail, the role that each Site plays in an optoelectronic molecular network.

Effect of Charge Distribution on the Dynamics of Polyampholytic Disordered Proteins

The stability and physiological function of many biomolecular coacervates depend on the structure and dynamics of intrinsically disordered proteins (IDPs) that typically contain a significant fraction of charged residues. Although the effect of relative arrangement of charged residues on IDP conformation is a well-studied problem, the associated changes in dynamics are far less understood. In this work, we systematically interrogate the effects of charge distribution on the chain-level and segmental dynamics of polyampholytic IDPs in dilute solutions. We study a coarse-grained model polyampholyte consisting of an equal fraction of two oppositely charged residues (glutamic acid and lysine) that undergoes a transition from an ideal chain-like conformation for uniformly charge-patterned sequences to a semi-compact conformation for highly charge-segregated sequences. Changes in the chain-level dynamics with increasing charge segregation correlate with changes in conformation. The chain-level and segmental dynamics conform to simple homopolymer models for uniformly charge-patterned sequences but deviate with increasing charge segregation, both in the presence and absence of hydrodynamic interactions. We discuss the significance of these findings, obtained for a model polyampholyte, in the context of a charge-rich intrinsically disordered region of the naturally occurring protein LAF-1. Our findings have important implications for understanding the effects of charge patterning on the dynamics of polyampholytic IDPs in dilute conditions using polymer scaling theories.

Molecular dynamics simulations of cyanine dimers attached to DNA Holliday junctions

Dye aggregates and their excitonic properties are of interest for their applications to organic photovoltaics, non-linear optics, and quantum information systems. DNA scaffolding has been shown to be effective at promoting the aggregation of dyes in a controllable manner. Specifically, isolated DNA Holliday junctions have been used to achieve strongly coupled cyanine dye dimers. However, the structural properties of the dimers and the DNA, as well as the role of Holliday junction isomerization are not fully understood. To study the dynamics of cyanine dimers in DNA, molecular dynamics simulations were carried out for adjacent and transverse dimers attached to Holliday junctions in two different isomers. It was found that dyes attached to adjacent strands in the junction exhibit stronger dye-DNA interactions and larger inter-dye separations compared to transversely attached dimers, as well as end-to-end arrangements. Transverse dimers exhibit lower inter-dye separations and more stacked configurations. Furthermore, differences in Holliday junction isomer are analyzed and compared to dye orientations. For transverse dyes exhibiting the smaller inter-dye separations, excitonic couplings were calculated and shown to be in agreement with experiment. Our results suggested that dye attachment locations on DNA Holliday junctions affect dye-DNA interactions, dye dynamics, and resultant dye orientations which can guide the design of DNA-templated cyanine dimers with desired properties.

Molecular dynamics simulations reveal dye attachment and DNA Holliday junction isomer effects on dye dimer orientations and excitonic couplings. These simulations can guide synthesis and experiments of dye-DNA structures for excitonic applications.

Developing Bonded Potentials for a Coarse-Grained Model of Intrinsically Disordered Proteins

Recent advances in residue-level coarse-grained (CG) computational models have enabled molecular-level insights into biological condensates of intrinsically disordered proteins (IDPs), shedding light on the sequence determinants of their phase separation. The existing CG models that treat protein chains as flexible molecules connected via harmonic bonds cannot populate common secondary-structure elements. Here, we present a CG dihedral angle potential between four neighboring beads centered at Cα atoms to faithfully capture the transient helical structures of IDPs. In order to parameterize and validate our new model, we propose Cα-based helix assignment rules based on dihedral angles that succeed in reproducing the atomistic helicity results of a polyalanine peptide and folded proteins. We then introduce sequence-dependent dihedral angle potential parameters (εd) and use experimentally available helical propensities of naturally occurring 20 amino acids to find their optimal values. The single-chain helical propensities from the CG simulations for commonly studied prion-like IDPs are in excellent agreement with the NMR-based α-helix fraction, demonstrating that the new HPS-SS model can accurately produce structural features of IDPs. Furthermore, this model can be easily implemented for large-scale assembly simulations due to its simplicity.

Strategies for Controlling the Spatial Orientation of Single Molecules Tethered on DNA Origami Templates Physisorbed on Glass Substrates: Intercalation and Stretching

Nanoarchitectural control of matter is crucial for next-generation technologies. DNA origami templates are harnessed to accurately position single molecules; however, direct single molecule evidence is lacking regarding how well DNA origami can control the orientation of such molecules in three-dimensional space, as well as the factors affecting control. Here, we present two strategies for controlling the polar (θ) and in-plane azimuthal (ϕ) angular orientations of cyanine Cy5 single molecules tethered on rationally-designed DNA origami templates that are physically adsorbed (physisorbed) on glass substrates. By using dipolar imaging to evaluate Cy5′s orientation and super-resolution microscopy, the absolute spatial orientation of Cy5 is calculated relative to the DNA template. The sequence-dependent partial intercalation of Cy5 is discovered and supported theoretically using density functional theory and molecular dynamics simulations, and it is harnessed as our first strategy to achieve θ control for a full revolution with dispersion as small as ±4.5°. In our second strategy, ϕ control is achieved by mechanically stretching the Cy5 from its two tethers, being the dispersion ±10.3° for full stretching. These results can in principle be applied to any single molecule, expanding in this way the capabilities of DNA as a functional templating material for single-molecule orientation control. The experimental and modeling insights provided herein will help engineer similar self-assembling molecular systems based on polymers, such as RNA and proteins.

Understanding Disorder, Vibronic Structure, and Delocalization in Electronically Coupled Dimers on DNA Duplexes

Structural DNA nanotechnology is a promising approach to create chromophore networks with modular structures and Hamiltonians to control the material's functions. The functional behaviors of these systems depend on the interactions of the chromophores' vibronic states, as well as interactions with their environment. To optimize their functions, it is necessary to characterize the chromophore network's structural and energetic properties, including the electronic delocalization in some cases. In this study, parameters of interest are deduced in DNA-scaffolded Cyanine 3 and Cyanine 5 dimers. The methods include steady-state optical measurements, physical modeling, and a genetic algorithm approach. The parameters include the chromophore network's vibronic Hamiltonian, molecular positions, transition dipole orientations, and environmentally induced energy broadening. Additionally, the study uses temperature-dependent optical measurements to characterize the spectral broadening further. These combined results reveal the quantum mechanical delocalization, which is important for functions like coherent energy transport and quantum information applications.

Induction of the integrated stress response in the rat cornea

Keratoconus (KC), a progressive, degenerative corneal disease, represents the second leading indication for corneal transplantation globally. We have previously demonstrated that components of the Integrated Stress Response (ISR) are upregulated in human keratoconic donor tissue, and treatment of normal tissue with ISR agonists attenuates collagen production. With no consistently accepted animal models available for translational KC research, we sought to establish an in vivo model based on ISR activation to elucidate its role in the development of the KC phenotype. Four-week-old female SD rats were treated with topical SAL003 formulated as a nanosuspension or vehicle every 48 h for four doses. Animals were subject to monitoring for ocular inflammation and discomfort before being euthanized at 1, 14, or 28 days after treatment was withdrawn. Schirmer's tear test, intraocular pressure, and body weight measurements were obtained at baseline and prior to euthanasia. Globes were subject to routine histopathology, immunohistochemistry for ATF4, and qPCR for Col1a1 expression. ANOVAs and Student's t tests were used to assess statistical significance (α = 0.05). SAL003 treatment did not produce any adverse ocular or systemic phenotype but did result in decreased keratocyte density. Col1a1 transcripts were reduced, corresponding to nuclear ATF4 expression within the axial cornea. In vivo topical treatment with a gel-formulated ISR agonist recapitulates key features of the activated ISR including nuclear ATF4 expression and decreased extracellular matrix (ECM) production. Exogenous ISR agonists may present one approach to establishing a rodent model for keratoconus, a charge essential for future evaluations of pathogenesis and therapeutic interventions.

Use of dual-electron probes reveals the role of ferritin as an iron depot in ex vivo erythropoiesis

In the finely regulated process of mammalian erythropoiesis, the path of the labile iron pool into mitochondria for heme production is not well understood. Existing models for erythropoiesis do not include a central role for the ubiquitous iron storage protein ferritin; one model proposes that incoming endosomal Fe³⁺ bound to transferrin enters the cytoplasm through an ion transporter after reduction to Fe²⁺ and is taken up into mitochondria through mitoferrin-1 transporter. Here, we apply a dual three-dimensional imaging and spectroscopic technique, based on scanned electron probes, to measure Fe³⁺ in ex vivo human hematopoietic stem cells. After seven days in culture, we observe cells displaying a highly specialized architecture with anchored clustering of mitochondria and massive accumulation of nanoparticles containing high iron concentrations localized to lysosomal storage depots, identified as ferritin. We hypothesize that lysosomal ferritin iron depots enable continued heme production after expulsion of most of the cellular machinery.

Improved coarse-grained model for studying sequence dependent phase separation of disordered proteins

We present improvements to the hydropathy scale (HPS) coarse-grained (CG) model for simulating sequence-specific behavior of intrinsically disordered proteins (IDPs), including their liquid-liquid phase separation (LLPS). The previous model based on an atomistic hydropathy scale by Kapcha and Rossky (KR scale) is not able to capture some well-known LLPS trends such as reduced phase separation propensity upon mutations (R-to-K and Y-to-F). Here, we propose to use the Urry hydropathy scale instead, which was derived from the inverse temperature transitions in a model polypeptide with guest residues X. We introduce two free parameters to shift (Δ) and scale (µ) the overall interaction strengths for the new model (HPS-Urry) and use the experimental radius of gyration for a diverse group of IDPs to find their optimal values. Interestingly, many possible (Δ, µ) combinations can be used for typical IDPs, but the phase behavior of a low-complexity (LC) sequence FUS is only well described by one of these models, which highlights the need for a careful validation strategy based on multiple proteins. The CG HPS-Urry model should enable accurate simulations of protein LLPS and provide a microscopically detailed view of molecular interactions.

A predictive coarse-grained model for position-specific effects of post-translational modifications

Biomolecules undergo liquid-liquid phase separation (LLPS), resulting in the formation of multicomponent protein-RNA membraneless organelles in cells. However, the physiological and pathological role of post-translational modifications (PTMs) on the biophysics of phase behavior is only beginning to be probed. To study the effect of PTMs on LLPS in silico, we extend our transferable coarse-grained model of intrinsically disordered proteins to include phosphorylated and acetylated amino acids. Using the parameters for modified amino acids available for fixed-charge atomistic force fields, we parameterize the size and atomistic hydropathy of the coarse-grained-modified amino acid beads and, hence, the interactions between the modified and natural amino acids. We then elucidate how the number and position of phosphorylated and acetylated residues alter the protein's single-chain compactness and its propensity to phase separate. We show that both the number and the position of phosphorylated threonines/serines or acetylated lysines can serve as a molecular on/off switch for phase separation in the well-studied disordered regions of Fused in Sarcoma (FUS) and DDX3X, respectively. We also compare modified residues to their commonly used PTM mimics for their impact on chain properties. Importantly, we show that the model can predict and capture experimentally measured differences in the phase behavior for position-specific modifications, showing that the position of modifications can dictate phase separation. In sum, this model will be useful for studying LLPS of post-translationally modified intrinsically disordered proteins and predicting how modifications control phase behavior with position-specific resolution.

Use of hydrochlorothiazide and risk of nonmelanoma skin cancer in Koreans: a retrospective cohort study using administrative healthcare data

BACKGROUND

The incidence of skin cancer is increasing because of the ageing population and ultraviolet exposure, and previous studies have revealed that long-term use of hydrochlorothiazide (HCTZ), an antihypertensive agent, has been associated with an increased risk of nonmelanoma skin cancer (NMSC). However, the association of NMSC and HCTZ within East Asian populations is unclear.

AIM

To investigate the risk of NMSC in Korean subjects using HCTZ.

METHODS

A retrospective cohort study was conducted using the administrative healthcare data. The study enrolled 62 243 patients exposed to HCTZ with a cumulative dose of ≥ 2500 mg and 62 243 unexposed subjects matched 1 : 1 with the patients for age, sex and income level.

RESULTS

There was a significant difference in the cumulative incidence of NMSC between the two groups (log-rank P < 0.01). Cox regression analysis was conducted after adjusting for potential confounders, and showed the risk for NMSC in the group exposed to HCTZ was significantly higher than that of the unexposed group (hazard ratio = 1.48; 95% CI 1.03-2.13). In the subgroup analysis, the oldest age group (≥ 70 years) showed increased cumulative incidence of NMSC with statistical significance compared with the unexposed control group (log-rank P < 0.01).

CONCLUSIONS

In this study, we revealed that the cumulative use of HCTZ (≥ 2500 mg) could increase the risk of NMSC in Koreans, especially the older age group. Thus, HCTZ could be a risk factor for NMSC in East Asian as well as white populations.

Molecular Details of Protein Condensates Probed by Microsecond Long Atomistic Simulations

The formation of membraneless organelles in cells commonly occurs via liquid-liquid phase separation (LLPS) and is in many cases driven by multivalent interactions between intrinsically disordered proteins (IDPs). Investigating the nature of these interactions, and their effect on dynamics within the condensed phase, is therefore of critical importance but very challenging for either simulation or experiment. Here, we study these interactions and their dynamics by pairing a novel multiscale simulation strategy with microsecond all-atom MD simulations of a condensed, IDP-rich phase. We simulate two IDPs this way, the low complexity domain of FUS and the N-terminal disordered domain of LAF-1, and find good agreement with experimental information about average density, water content, and residue-residue contacts. We go significantly beyond what is known from experiments by showing that ion partitioning within the condensed phase is largely driven by the charge distribution of the proteins and-in the cases considered-shows little evidence of preferential interactions of the ions with the proteins. Furthermore, we can probe the microscopic diffusive dynamics within the condensed phase, showing that water and ions are in dynamic equilibrium between dense and dilute phases, and their diffusion is reduced in the dense phase. Despite their high concentration in the condensate, the protein molecules also remain mobile, explaining the observed liquid-like properties of this phase. We finally show that IDP self-association is driven by a combination of nonspecific hydrophobic interactions as well as hydrogen bonds, salt bridges, and π-π and cation-π interactions. The simulation approach presented here allows the structural and dynamical properties of biomolecular condensates to be studied in microscopic detail and is generally applicable to single- and multicomponent systems of proteins and nucleic acids involved in LLPS.

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.

DNA scaffold supports long-lived vibronic coherence in an indodicarbocyanine (Cy5) dimer

Vibronic coupling between pigment molecules is believed to prolong coherences in photosynthetic pigment–protein complexes. Reproducing long-lived coherences using vibronically coupled chromophores in synthetic DNA constructs presents a biomimetic route to efficient artificial light harvesting. Here, we present two-dimensional (2D) electronic spectra of one monomeric Cy5 construct and two dimeric Cy5 constructs (0 bp and 1 bp between dyes) on a DNA scaffold and perform beating frequency analysis to interpret observed coherences. Power spectra of quantum beating signals of the dimers reveal high frequency oscillations that correspond to coherences between vibronic exciton states. Beating frequency maps confirm that these oscillations, 1270 cm⁻¹ and 1545 cm⁻¹ for the 0-bp dimer and 1100 cm⁻¹ for the 1-bp dimer, are coherences between vibronic exciton states and that these coherences persist for ∼300 fs. Our observations are well described by a vibronic exciton model, which predicts the excitonic coupling strength in the dimers and the resulting molecular exciton states. The energy spacing between those states closely corresponds to the observed beat frequencies. MD simulations indicate that the dyes in our constructs lie largely internal to the DNA base stacking region, similar to the native design of biological light harvesting complexes. Observed coherences persist on the timescale of photosynthetic energy transfer yielding further parallels to observed biological coherences, establishing DNA as an attractive scaffold for synthetic light harvesting applications.

Dyes coupled to DNA display distance-dependent vibronic couplings that prolongs quantum coherences detected with 2D spectroscopy.

Hydropathy Patterning Complements Charge Patterning to Describe Conformational Preferences of Disordered Proteins

Understanding the conformational ensemble of an intrinsically disordered protein (IDP) is of great interest due to its relevance to critical intracellular functions and diseases. It is now well established that the polymer scaling behavior can provide a great deal of information about the conformational properties as well as liquid-liquid phase separation of an IDP. It is, therefore, extremely desirable to be able to predict an IDP's scaling behavior from the protein sequence itself. The work in this direction so far has focused on highly charged proteins and how charge patterning can perturb their structural properties. As naturally occurring IDPs are composed of a significant fraction of uncharged amino acids, the rules based on charge content and patterning are only partially helpful in solving the problem. Here, we propose a new order parameter, sequence hydropathy decoration, which can provide a near-quantitative understanding of scaling and structural properties of IDPs devoid of charged residues. We combine this with a charge patterning parameter, sequence charge decoration, to obtain a general equation, parametrized from extensive coarse-grained simulation data, for predicting protein dimensions from the sequence. We finally test this equation against available experimental data and find a semiquantitative match in predicting the scaling behavior. We also provide guidance on how to extend this approach to experimental data, which should be feasible in the near future.

Upper Critical Field Based on a Width of ΔH = ΔB region in a Superconductor

We studied a method of measuring upper critical field (H_c2) of a superconductor based on a width of ΔH = ΔB region, which appears in a superconductor that volume defects are many and dominant. Here we show basic concepts and details of the method. Although H_c2 of a superconductor is fixed according to a kind of superconductor, it is difficult to measure H_c2 experimentally. Thus, results are different depending on experimental conditions. H_c2 was otained by a theory on a width of ΔH = ΔB region, which is that pinned fluxes at volume defects are picked out and move into an inside of the superconductor when the distance between pinned fluxes is the same as that at H_c2 of the superconductor. H_c2 of MgB₂ obtained by the method was 65.4 Tesla at 0 K, which is quite same as that of Ginzburg-Landau theory. The reason that H_c2 obtained by the method is closer to ultimate H_c2 is based on that F_pinning/F_pickout is more than 4 when pinned fluxes at volume defects of 163 nm radius are depinned, which means that the H_c2 is less sensitive to fluctuation. The method will help to find the ultimate H_c2 of volume defect-dominating superconductors.

Transplant outcomes in kidney recipients with lupus nephritis, and systematic review

Background

Despite improved survival of patients with lupus nephritis (LN), some require kidney transplantation because of progression to end-stage renal disease (ESRD). However, the transplant outcomes of these patients and other recipients have not been thoroughly compared.

Methods

In total, 1848 Korean kidney recipients who underwent transplantation from 1998 to 2017 at two tertiary referral centers were evaluated retrospectively. Among them, 28 recipients with LN, and 50 control recipients matched by age, sex, and donor type, were compared with respect to graft and patient survival. We pooled our data with 17 previous cohort studies in which the graft survival of recipients with LN was described in detail.

Results

During the median follow-up period of 9.5 years (maximum 21 years), graft failure (GF) occurred in 10.7% and 16.0% of LN and control recipients, respectively. No differences were found in the rates of GF and death-censored graft failure or patient survival between the two groups. The risks of acute T cell-mediated and antibody-mediated rejection were also similar between the two groups. The pooled analysis showed similar 1- and 5-year graft survival rates between LN and control recipients.

Conclusions

Kidney transplantation is an acceptable option in patients with concurrent LN and ESRD.

Temperature-Controlled Liquid-Liquid Phase Separation of Disordered Proteins

The liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) is a commonly observed phenomenon within the cell, and such condensates are also highly attractive for applications in biomaterials and drug delivery. A better understanding of the sequence-dependent thermoresponsive behavior is of immense interest as it will aid in the design of protein sequences with desirable properties and in the understanding of cellular response to heat stress. In this work, we use a transferable coarse-grained model to directly probe the sequence-dependent thermoresponsive phase behavior of IDPs. To achieve this goal, we develop a unique knowledge-based amino acid potential that accounts for the temperature-dependent effects on solvent-mediated interactions for different types of amino acids. Remarkably, we are able to distinguish between more than 35 IDPs with upper or lower critical solution temperatures at experimental conditions, thus providing direct evidence that incorporating the temperature-dependent solvent-mediated interactions to IDP assemblies can capture the difference in the shape of the resulting phase diagrams. Given the success of the model in predicting experimental behavior, we use it as a high-throughput screening framework to scan through millions of disordered sequences to characterize the composition dependence of protein phase separation.

The Differential Flow of Epidural Local Anaesthetic via Needle or Catheter: A Prospective Randomized Double-blind Study

The extent of epidural anaesthesia and pattern of spread of contrast medium, using different injection techniques, has not been well documented. Therefore, in this prospective, randomized, double-blind study, the extent of anaesthesia and pattern of spread of contrast medium following an epidural bolus injection, via either a Tuohy needle or an epidural catheter, were compared.

The study had two parts. In the first, 59 of 79 patients scheduled for a lower extremity operation under epidural anaesthesia were randomly allocated to one of the two groups. Anaesthesia was achieved with an epidural injection of 10 to 15 ml (including a 3 ml test dose) of 0.75% ropivacaine and fentanyl 25 μg via either a Tuohy needle (Group N, n=31) or a catheter (Group C, n=28). The level of sensory anaesthesia was recorded. In the second part, the remaining 20 patients were randomized to initially receive 5 ml of contrast medium via either a Tuohy needle (Group NE, n=10) or a catheter (Group CE, n=10). The extent of spread was recorded radiologically.

Unilateral or missed blocks and additional dose requirement were absent in Groups N and C. No differences were found in the extent of sensory anaesthesia or the spread of contrast medium. Twenty per cent of catheter tips lay outside the lateral margins of the vertebral bodies.

We found that an epidural bolus injection, via either a Tuohy needle or a catheter, made no difference in regard to spread of local anaesthetic or contrast medium in the epidural space.

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.

A natural infection by the red sea bream iridovirus-type Megalocytivirus in the golden mandarin fish Siniperca scherzeri

An outbreak of a Megalocytivirus infection was found in the golden mandarin fish Siniperca scherzeri during September and October 2016, in Korea. Phylogeny and genetic diversity based on the major capsid protein (MCP) and adenosine triphosphatase (ATPase) genes showed a new strain. Designated as GMIV, this strain derived from the golden mandarin fish was suggested to belong to the red sea bream iridovirus (RSIV)-subgroup I. Additionally, this train clustered with the ehime-1 strain from red sea bream Pagrus major in Japan and was distinguished from circulating isolates (RSIV-type subgroup II and turbot reddish body iridovirus [TRBIV] type) in Korea. The infection level, evaluated by qPCR, ranged from 8.18 × 102 to 7.95 × 106  copies/mg of tissue individually, suggesting that the infected fish were in the disease-transmitting stage. The diseased fish showed degenerative changes associated with cytomegaly in the spleen as general sign of Megalocytivirus infection. The results confirm that the RSIV-type Megalocytivirus might have crossed the environmental and species barriers to cause widespread infection in freshwater fish.

Computational modeling highlights the role of the disordered Formin Homology 1 domain in profilin-actin transfer

Formins accelerate actin polymerization, assumed to occur through flexible Formin Homology 1 (FH1) domain-mediated transfer of profilin-actin to the barbed end. To study FH1 properties and address sequence effects, including varying length/distribution of profilin-binding proline-rich motifs, we performed all-atom simulations of a set of representative FH1 domains of formins: mouse mDia1 and mDia2, budding yeast Bni1 and Bnr1, and fission yeast Cdc12, For3, and Fus1. We find FH1 has flexible regions between high-propensity polyproline helix regions. A coarse-grained model retaining sequence specificity, assuming rigid polyproline segments, describes their size. Multiple bound profilins or profilin-actin complexes expand mDia1-FH1, which may be important in cells. Simulations of the barbed end bound to Bni1-FH1-FH2 dimer show that the leading FH1 can better transfer profilin or profilin-actin, with decreasing probability as the distance from FH2 increases.

A multicentre, open, investigator-initiated phase IV clinical trial to evaluate the efficacy and safety of ingenol mebutate gel, 0·015% on the face and scalp, and 0·05% on the trunk and extremities, in Korean patients with actinic keratosis (PERFECT)

BACKGROUND

Ingenol mebutate gel is a novel, field-directed topical treatment for actinic keratosis (AK). Most pivotal studies have targeted Western populations. No clinical study has been conducted to investigate its efficacy and safety in Asian populations.

OBJECTIVES

To evaluate the efficacy and safety of ingenol mebutate gel for treating AK of face/scalp and trunk/extremities in a large Asian (Korean) population.

PATIENTS AND METHODS

In this multicentre, open-label, interventional, parallel-group, prospective phase IV study (PERFECT, trial registration no.: NCT02716714), the eligible patients were allocated into either the face/scalp or the trunk/extremities group, according to their selected treatment area location. After application of ingenol mebutate gel, the participants were followed up for 6 months. The primary efficacy endpoint was complete clearance (CC) of AK lesions in the selected treatment area at day 57. Quality of life was evaluated using Skindex-29. Safety endpoints included local skin responses, scar, pigmentation, pain and adverse events.

RESULTS

In total, 78·1% [95% confidence interval (CI) 66·86-86·92%] of subjects had CC at day 57, with 76·6% (95% CI 64·31-86·25%) in the face/scalp group and 88·9% (95% CI 51·75-99·72%) in the trunk/extremities group. Among them, CC was sustained in 88·9% (48 of 54, 95% CI 77·37-95·81%) at month 6. The local skin responses significantly increased 1 day after the treatment compared with baseline, and decreased afterwards. Among the total subjects, 7·8% (6 of 77) had hyperpigmentation on the application area. Scars were not reported.

CONCLUSIONS

Ingenol mebutate is effective for the treatment of AK in Asians, with tolerable safety profiles.

Isolation and initial characterization of new betanodaviruses in shellfish

Betanodaviruses cause the disease viral nervous necrosis (VNN) in finfish. Using a novel approach with two consecutive PCRs, detection semi-nested two-step RT-PCR (DSN-2 RT-PCR) and discriminative multiplex two-step RT-PCR (DMT-2 RT-PCR), we have identified the presence of a new type of betanodavirus in shellfish and called it Korean shellfish nervous necrosis virus (KSNNV). Partial nucleotide sequences of the T4 region in RNA2 fragment of KSNNVs were 73%-75% homologous to those of other reported genotypes and formed a new cluster of betanodavirus in phylogenetic tree analysis. Successful isolation of KSNNV was achieved in two of six shellfish samples containing high concentrations of virus using the blind passage method, and the typical shapes of betanodavirus were confirmed in KSNNV-KOR1 by electron microscopy. In the experimental infection test, seven of 14 fish species showed susceptibility to KSNNV-KOR1 isolate but without clinical signs or death. Although the range of susceptible host species was not significantly different from the RGNNV type, the concentration of KSNNV in the brain of infected fish (10² -10⁵ copies/mg brain) was much lower compared to that found in sevenband grouper (Epinephelus septemfasciatus Thunberg) sampled in the moribund stage with RGNNV infection (10⁶ -10⁷ copies/mg brain). However, histopathological analyses showed the presence of multiple vacuoles in brains of all KSNNV-infected fish at 14 days postinjection. In detection test, as a single or multiple type with the other genotype(s) (RGNNV or BFNNV), the prevalence of KSNNV was 8.4% and 8.7% in domestic (62 of 741 samples) and Chinese samples (12 of 138 samples), respectively, but not in finfish. We propose that KSNNVs obtained from shellfish be classified into a separate and new genotype of betanodavirus.

Optical Properties of Vibronically Coupled Cy3 Dimers on DNA Scaffolds

We examine the effect of electronic coupling on the optical properties of Cy3 dimers attached to DNA duplexes as a function of base pair (bp) separation using steady-state and time-resolved spectroscopy. For close Cy3-Cy3 separations, 0 and 1 bp between dyes, intermediate to strong electronic coupling is revealed by modulation of the absorption and fluorescence properties including spectral band shape, peak wavelength, and excited-state lifetime. Using a vibronic exciton model, we estimate coupling strengths of 150 and 266 cm^-1 for the 1 and 0 bp separations, respectively, which are comparable to those found in natural light-harvesting complexes. For the strongest electronic coupling (0 bp separation), we observe that the absorption band shape is strongly affected by the base pairs that surround the dyes, where more strongly hydrogen-bonded G-C pairs produce a red-shifted absorption spectrum consistent with a J-type dimer. This effect is studied theoretically using molecular dynamics simulation, which predicts an in-line dye configuration that is consistent with the experimental J-type spectrum. When the Cy3 dimers are in a standard aqueous buffer, the presence of relatively strong electronic coupling is accompanied by decreased fluorescence lifetime, suggesting that it promotes nonradiative relaxation in cyanine dyes. However, we show that the use of a viscous solvent can suppress this nonradiative recombination and thereby restore the dimer fluorescent emission. Ultrafast transient absorption measurements of Cy3 dimers in both standard aqueous buffer and viscous glycerol buffer suggest that sufficiently strong electronic coupling increases the probability of excited-state relaxation through a dark state that is related to Cy3 torsional motion.

High prevalence of betanodavirus barfin flounder nervous necrosis virus as well as red-spotted grouper nervous necrosis virus genotype in shellfish

Using two serially executed PCRs, the discriminative multiplex two-step RT-PCR (DMT-2 RT-PCR) following the detection seminested two-step RT-PCR (DSN-2 RT-PCR), we found a high frequency presence of BFNNV genotype as well as RGNNV in various domestic and imported shellfish. This was definitely different from the previous reports of outbreaks and asymptomatic infection only by the RGNNV genotype in cultured finfish in Korea. Cultivation of NNV entrapped in shellfish was performed successfully by a blind passage. Thus, in an attempt to elucidate the epidemiology of betanodavirus, experiments conducted on 969 shellfish samples concluded that (i) distribution of NNV genotype, especially BFNNV, in shellfish is clearly different from that found in finfish of the world; (ii) unlike RGNNV, which showed a high rate in summer, BFNNV showed no seasonal variation and this result suggests BFNNVs in the marine environment remain fairly constant throughout the year; and (iii) the entrapped virus in shellfish was alive and culturable in vitro. These results are the first report of high level prevalence of in vitro culturable NNV in shellfish, for both BFNNV and RGNNV, which may present a potential risk in transmitting nodaviruses to host species in a marine environment.

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.

Liquid liquid phase separation (LLPS) of low-complexity protein sequences has emerged as an important research topic due to its relevance to membraneless organelles and intracellular compartmentalization. However a molecular level understanding of LLPS cannot be easily obtained by experimental methods due to difficulty of determining structural properties of phase separated protein assemblies, and of choosing appropriate mutations. Here we advance a coarse-grained computational framework for accessing the long time scale phase separation process and for obtaining molecular details of LLPS, in conjunction with state of the art enhanced sampling methods. We are able to qualitatively capture the changes of the phase diagram due to specific mutations, inclusion of a folded domain, and variation of chain length. The model is flexible and can be used with different knowledge-based potential energy functions, as we demonstrate. We expect a wide application of the presented framework for advancing our understanding of the formation of liquid-like protein assemblies.

Comparative evaluation of MCP gene in worldwide strains of Megalocytivirus (Iridoviridae family) for early diagnostic marker

Members of the Iridoviridae family have been considered as aetiological agents of iridovirus diseases, causing fish mortalities and economic losses all over the world. Virus identification based on candidate gene sequencing is faster, more accurate and more reliable than other traditional phenotype methodologies. Iridoviridae viruses are covered by a protein shell (capsid) encoded by the important candidate gene, major capsid protein (MCP). In this study, we investigated the potential of the MCP gene for use in the diagnosis and identification of infections caused Megalocytivirus of the Iridoviridae family. We selected data of 66 Iridoviridae family isolates (53 strains of Megalocytivirus, eight strains of iridoviruses and five strains of Ranavirus) infecting various species of fish distributed all over the world. A total of 53 strains of Megalocytivirus were used for designing the complete primer sets for identifying the most hypervariable region of the MCP gene. Further, our in silico analysis of 102 sequences of related and unrelated viruses reconfirms that primer sets could identify strains more specifically and offers a useful and fast alternative for routine clinical laboratory testing. Our findings suggest that phenotype observation along with diagnosis using universal primer sets can help detect infection or carriers at an early stage.

Protein Composition Determines the Effect of Crowding on the Properties of Disordered Proteins

Unlike dilute experimental conditions under which biological molecules are typically characterized, the cell interior is crowded by macromolecules, which affects both the thermodynamics and kinetics of in vivo processes. Although the excluded-volume effects of macromolecular crowding are expected to cause compaction of unfolded and disordered proteins, the extent of this effect is uncertain. We use a coarse-grained model to represent proteins with varying sequence content and directly observe changes in chain dimensions in the presence of purely repulsive spherical crowders. We find that the extent of crowding-induced compaction is dependent not only on crowder size and concentration, but also on the properties of the protein itself. In fact, we observe a nonmonotonic trend between the dimensions of the polypeptide chain in bulk and the degree of compaction: the most extended chains experience up to 24% compaction, the most compact chains show virtually no change, and intermediate chains compress by up to 40% in size at a 40% crowder volume fraction. Free-volume theory combined with an impenetrable ellipsoidal representation of the chains predicts the crowding effects only for collapsed protein chains. An additional scaling factor, which can be easily computed from protein-crowder potential of mean force, corrects for the penetrability of extended chains and is sufficient to capture the observed nonmonotonic trend in compaction.

Delta-frequency stimulation of cerebellar projections can compensate for schizophrenia-related medial frontal dysfunction

Schizophrenia involves abnormalities in the medial frontal cortex that lead to cognitive deficits. Here we investigate a novel strategy to normalize medial frontal brain activity by stimulating cerebellar projections. We used an interval timing task to study elementary cognitive processing that requires both frontal and cerebellar networks that are disrupted in patients with schizophrenia. We report three novel findings. First, patients with schizophrenia had dysfunctional delta rhythms between 1-4 Hz in the medial frontal cortex. We explored cerebellar-frontal interactions in animal models and found that both frontal and cerebellar neurons were modulated during interval timing and had delta-frequency interactions. Finally, delta-frequency optogenetic stimulation of thalamic synaptic terminals of lateral cerebellar projection neurons rescued timing performance as well as medial frontal activity in a rodent model of schizophrenia-related frontal dysfunction. These data provide insight into how the cerebellum influences medial frontal networks and the role of the cerebellum in cognitive processing.

Cloning and expression analysis of innate immune genes from red sea bream to assess different susceptibility to megalocytivirus infection

As suggested by the Office International des Epizooties (OIE), fishes belonging to the genus Oplegnathus are more sensitive to megalocytivirus infection than other fish species including red sea bream (Pagrus major). To assess the roles of the innate immune response to these different susceptibilities, we cloned the genes encoding inflammatory factors including IL-8 and COX-2, and the antiviral factor like Mx from red sea bream for the first time and performed phylogenetic and structural analysis. Analysed expression levels of IL-1β, IL-8 and COX-2 and the antiviral factor like Mx genes performed with in vivo challenge experiment showed no difference in inflammatory gene expression or respiratory burst activity between red sea bream and rock bream (Oplegnathus fasciatus). However, the Mx gene expression levels in red sea bream were markedly higher than those in rock bream, suggesting the importance of type I interferon (IFN)-induced proteins, particularly Mx, during megalocytivirus infection, rather than inflammation-related genes. The in vitro challenge experiments using embryonic primary cultures derived from both fish species showed no difference in cytopathic effects (CPE), viral replication profiles, and inflammatory and Mx gene expression pattern between the two fish species.

Cyclic fatigue resistance, torsional resistance, and metallurgical characteristics of V taper 2 and V taper 2H rotary NiTi files

The aim of this study was to compare the cyclic fatigue resistance, torsional resistance, and metallurgical characteristics of conventional NiTi wire (V taper 2, V2) and CM wire (V taper 2H, V2H)-based files. Cyclic fatigue and torsional resistance of V2 and V2H were investigated by measuring the number of cycles to fracture, maximum torque at fracture, and maximum angle at fracture. The typical patterns of fatigue and torsional fractures were investigated using a scanning electron microscope (SEM). The metallurgical characteristics were investigated by differential scanning calorimetry (DSC) from -100 °C to 100 °C. The austenite finishing temperature (Af) of each instrument was also measured. The microstructures of the instruments were investigated by a transmission electron microscope (TEM) along with selected area diffraction pattern analysis. The results were statistically analyzed by Mann-Whitney U-test (p = 0.05). V2H showed significantly higher cyclic fatigue resistance and torsional resistance than V2. SEM images of the fractured surfaces showed typical patterns of fatigue and torsional fracture. The DSC analysis of V2 showed one small peak in both the heating and cooling curves. The Af of V2 was -0.32 °C. V2H showed two remarkable peaks in the heating curve and one remarkable peak in the cooling curve. The Af of V2H was 33.25 °C. The TEM analysis showed that both V2 and V2H are mainly composed of austenite. In conclusion, V2H showed higher cyclic fatigue resistance and torsional resistance than V2. The superior properties of V2H could be attributed to the annealing effect and possibly the martensite phase. SCANNING 38:564-570, 2016. © 2016 Wiley Periodicals, Inc.

Efficacy of Pulsed Electromagnetic Therapy for Chronic Lower Back Pain: A Randomized, Double-blind, Placebo-controlled Study

This randomized, double-blind, placebo-controlled clinical trial studied the effectiveness of pulsed electromagnetic therapy (PEMT) in patients with chronic lower back pain. Active PEMT ( n = 17) or placebo treatment ( n = 19) was performed three times a week for 3 weeks. Patients were assessed using a numerical rating scale (NRS) and revised Oswestry disability scores for 4 weeks after therapy. PEMT produced significant pain reduction throughout the observation period compared with baseline values. The percentage change in the NRS score from baseline was significantly greater in the PEMT group than the placebo group at all three time-points measured. The mean revised Oswestry disability percentage after 4 weeks was significantly improved from the baseline value in the PEMT group, whereas there were no significant differences in the placebo group. In conclusion, PEMT reduced pain and disability and appears to be a potentially useful therapeutic tool for the conservative management of chronic lower back pain.