C9orf72 polyPR directly binds to various nuclear transport components.

The disruption of nucleocytoplasmic transport (NCT) is an important mechanism in neurodegenerative diseases. In the case of C9orf72-ALS, trafficking of macromolecules through the nuclear pore complex (NPC) might get frustrated by the binding of C9orf72-translated arginine-containing dipeptide repeat proteins (R-DPRs) to the Kapß family of nuclear transport receptors. Besides Kapßs, several other types of transport components have been linked to NCT impairments in R-DPR-expressed cells, but the molecular origin of these observations has not been clarified. Here, we adopt a coarse-grained molecular dynamics model at amino acid resolution to study the direct interaction between polyPR, the most toxic DPR, and various nuclear transport components to elucidate the binding mechanisms and provide a complete picture of potential polyPR-mediated NCT defects. We found polyPR to directly bind to several isoforms of the Impα family, CAS (the specific exporter of Impα) and RanGAP. We observe no binding between polyPR and Ran. Longer polyPRs at lower salt concentrations also make contact with RanGEF and NTF2. Analyzing the polyPR contact sites on the transport components reveals that polyPR potentially interferes with RanGTP/RanGDP binding, with nuclear localization signal (NLS)-containing cargoes (cargo-NLS) binding to Impα, with cargo-NLS release from Impα, and with Impα export from the nucleus. The abundance of polyPR-binding sites on multiple transport components combined with the inherent polyPR length dependence makes direct polyPR interference of NCT a potential mechanistic pathway of C9orf72 toxicity.

A One-Bead-Per-Saccharide (1BPS) Model for Glycosaminoglycans.

Glycosaminoglycans (GAGs) are polysaccharide compounds that play key roles in various biological processes. GAGs are important structural components of cartilage and the extracellular matrix of the brain. Due to the large size of these polysaccharides, coarse-grained approaches are indispensable for modeling these biopolymers. We develop a one-bead-per-saccharide model of chondroitin sulfates and hyaluronic acid based on an existing three-bead-per-saccharide coarse-grained model. Our coarse graining is carried out by using iterative Boltzmann inversion (IBI), including an additional coupling potential to incorporate the correlation between dihedral angles. The predictions of the model are verified against those of the existing three-bead-per-saccharin model and the experimental radius of gyration for hyaluronic acid.

Phase separation of intrinsically disordered FG-Nups is driven by highly dynamic FG motifs.

The intrinsically disordered FG-Nups in the central channel of the nuclear pore complex (NPC) form a selective permeability barrier, allowing small molecules to traverse by passive diffusion, while large molecules can only translocate with the help of nuclear transport receptors. The exact phase state of the permeability barrier remains elusive. In vitro experiments have shown that some FG-Nups can undergo phase separation into condensates that display NPC-like permeability barrier properties. Here, we use molecular dynamics simulations at amino acid resolution to study the phase separation characteristics of each of the disordered FG-Nups of the yeast NPC. We find that GLFG-Nups undergo phase separation and reveal that the FG motifs act as highly dynamic hydrophobic stickers that are essential for the formation of FG-Nup condensates featuring droplet-spanning percolated networks. Additionally, we study phase separation in an FG-Nup mixture that resembles the NPC stoichiometry and observe that an NPC condensate is formed containing multiple GLFG-Nups. We find that the phase separation of this NPC condensate is also driven by FG-FG interactions, similar to the homotypic FG-Nup condensates. Based on the observed phase separation behavior, the different FG-Nups of the yeast NPC can be divided into two classes: The FG-Nups (mostly GLFG-type) located in the central channel of the NPC form a highly dynamic percolated network formed by many short-lived FG-FG interactions, while the peripheral FG-Nups (mostly FxFG-type) at the entry and exit of the NPC channel likely form an entropic brush.

Molecular basis of C9orf72 poly-PR interference with the ß-karyopherin family of nuclear transport receptors.

Nucleocytoplasmic transport (NCT) is affected in several neurodegenerative diseases including C9orf72-ALS. It has recently been found that arginine-containing dipeptide repeat proteins (R-DPRs), translated from C9orf72 repeat expansions, directly bind to several importins. To gain insight into how this can affect nucleocytoplasmic transport, we use coarse-grained molecular dynamics simulations to study the molecular interaction of poly-PR, the most toxic DPR, with several Kapßs (importins and exportins). We show that poly-PR-Kapß binding depends on the net charge per residue (NCPR) of the Kapß, salt concentration of the solvent, and poly-PR length. Poly-PR makes contact with the inner surface of most importins, which strongly interferes with Kapß binding to cargo-NLS, IBB, and RanGTP in a poly-PR length-dependent manner. Longer poly-PRs at higher concentrations are also able to make contact with the outer surface of importins that contain several binding sites to FG-Nups. We also show that poly-PR binds to exportins, especially at lower salt concentrations, interacting with several RanGTP and FG-Nup binding sites. Overall, our results suggest that poly-PR might cause length-dependent defects in cargo loading, cargo release, Kapß transport and Ran gradient across the nuclear envelope.

A designer FG-Nup that reconstitutes the selective transport barrier of the nuclear pore complex.

Nuclear Pore Complexes (NPCs) regulate bidirectional transport between the nucleus and the cytoplasm. Intrinsically disordered FG-Nups line the NPC lumen and form a selective barrier, where transport of most proteins is inhibited whereas specific transporter proteins freely pass. The mechanism underlying selective transport through the NPC is still debated. Here, we reconstitute the selective behaviour of the NPC bottom-up by introducing a rationally designed artificial FG-Nup that mimics natural Nups. Using QCM-D, we measure selective binding of the artificial FG-Nup brushes to the transport receptor Kap95 over cytosolic proteins such as BSA. Solid-state nanopores with the artificial FG-Nups lining their inner walls support fast translocation of Kap95 while blocking BSA, thus demonstrating selectivity. Coarse-grained molecular dynamics simulations highlight the formation of a selective meshwork with densities comparable to native NPCs. Our findings show that simple design rules can recapitulate the selective behaviour of native FG-Nups and demonstrate that no specific spacer sequence nor a spatial segregation of different FG-motif types are needed to create selective NPCs.

Phase Separation of Toxic Dipeptide Repeat Proteins Related to C9orf72 ALS/FTD.

The expansion mutation in the C9orf72 gene is the most common known genetic cause for amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). This mutation can produce five dipeptide repeat proteins (DPRs), of which three are known to be toxic: poly-PR, poly-GR, and poly-GA. The toxicity of poly-GA is attributed to its aggregation in the cytoplasm, whereas for poly-PR and poly-GR, several toxicity pathways have been proposed. The toxicity of the DPRs has been shown to depend on their length, but the underlying molecular mechanism of this length dependence is not well understood. To address the possible role of phase separation in DPR toxicity, a one-bead-per-amino-acid (1BPA) coarse-grained molecular dynamics model is used to study the single-molecule and phase-separation properties of the DPRs. We find a strong dependence of the phase-separation behavior on both DPR length and concentration, with longer DPRs having a higher propensity to phase separate and form condensed phases with higher concentrations. The critical lengths required for phase separation (25 for poly-PR and 50 for poly-GA) are comparable to the toxicity threshold limit of 30 repeats found for the expansion mutation in patient cells, suggesting that phase separation could play an important role in DPR toxicity.

Age-dependent deterioration of nuclear pore assembly in mitotic cells decreases transport dynamics.

Nuclear transport is facilitated by the Nuclear Pore Complex (NPC) and is essential for life in eukaryotes. The NPC is a long-lived and exceptionally large structure. We asked whether NPC quality control is compromised in aging mitotic cells. Our images of single yeast cells during aging, show that the abundance of several NPC components and NPC assembly factors decreases. Additionally, the single-cell life histories reveal that cells that better maintain those components are longer lived. The presence of herniations at the nuclear envelope of aged cells suggests that misassembled NPCs are accumulated in aged cells. Aged cells show decreased dynamics of transcription factor shuttling and increased nuclear compartmentalization. These functional changes are likely caused by the presence of misassembled NPCs, as we find that two NPC assembly mutants show similar transport phenotypes as aged cells. We conclude that NPC interphase assembly is a major challenge for aging mitotic cells.

The Effect of FG-Nup Phosphorylation on NPC Selectivity: A One-Bead-Per-Amino-Acid Molecular Dynamics Study.

Nuclear pore complexes (NPCs) are large protein complexes embedded in the nuclear envelope separating the cytoplasm from the nucleoplasm in eukaryotic cells. They function as selective gates for the transport of molecules in and out of the nucleus. The inner wall of the NPC is coated with intrinsically disordered proteins rich in phenylalanine-glycine repeats (FG-repeats), which are responsible for the intriguing selectivity of NPCs. The phosphorylation state of the FG-Nups is controlled by kinases and phosphatases. In the current study, we extended our one-bead-per-amino-acid (1BPA) model for intrinsically disordered proteins to account for phosphorylation. With this, we performed molecular dynamics simulations to probe the effect of phosphorylation on the Stokes radius of isolated FG-Nups, and on the structure and transport properties of the NPC. Our results indicate that phosphorylation causes a reduced attraction between the residues, leading to an extension of the FG-Nups and the formation of a significantly less dense FG-network inside the NPC. Furthermore, our simulations show that upon phosphorylation, the transport rate of inert molecules increases, while that of nuclear transport receptors decreases, which can be rationalized in terms of modified hydrophobic, electrostatic, and steric interactions. Altogether, our models provide a molecular framework to explain how extensive phosphorylation of FG-Nups decreases the selectivity of the NPC.

Hemagglutinin-Mediated Membrane Fusion: A Biophysical Perspective.

Influenza hemagglutinin (HA) is a viral membrane protein responsible for the initial steps of the entry of influenza virus into the host cell. It mediates binding of the virus particle to the host-cell membrane and catalyzes fusion of the viral membrane with that of the host. HA is therefore a major target in the development of antiviral strategies. The fusion of two membranes involves high activation barriers and proceeds through several intermediate states. Here, we provide a biophysical description of the membrane fusion process, relating its kinetic and thermodynamic properties to the large conformational changes taking place in HA and placing these in the context of multiple HA proteins working together to mediate fusion. Furthermore, we highlight the role of novel single-particle experiments and computational approaches in understanding the fusion process and their complementarity with other biophysical approaches.

Spatial structure of disordered proteins dictates conductance and selectivity in nuclear pore complex mimics.

Nuclear pore complexes (NPCs) lined with intrinsically disordered FG-domains act as selective gatekeepers for molecular transport between the nucleus and the cytoplasm in eukaryotic cells. The underlying physical mechanism of the intriguing selectivity is still under debate. Here, we probe the transport of ions and transport receptors through biomimetic NPCs consisting of Nsp1 domains attached to the inner surface of solid-state nanopores. We examine both wildtype FG-domains and hydrophilic SG-mutants. FG-nanopores showed a clear selectivity as transport receptors can translocate across the pore whereas other proteins cannot. SG mutant pores lack such selectivity. To unravel this striking difference, we present coarse-grained molecular dynamics simulations that reveal that FG-pores exhibit a high-density, nonuniform protein distribution, in contrast to a uniform and significantly less-dense protein distribution in the SG-mutant. We conclude that the sequence-dependent density distribution of disordered proteins inside the NPC plays a key role for its conductivity and selective permeability.

Artificial muscle-like function from hierarchical supramolecular assembly of photoresponsive molecular motors.

A striking feature of living systems is their ability to produce motility by amplification of collective molecular motion from the nanoscale up to macroscopic dimensions. Some of nature's protein motors, such as myosin in muscle tissue, consist of a hierarchical supramolecular assembly of very large proteins, in which mechanical stress induces a coordinated movement. However, artificial molecular muscles have often relied on covalent polymer-based actuators. Here, we describe the macroscopic contractile muscle-like motion of a supramolecular system (comprising 95% water) formed by the hierarchical self-assembly of a photoresponsive amphiphilic molecular motor. The molecular motor first assembles into nanofibres, which further assemble into aligned bundles that make up centimetre-long strings. Irradiation induces rotary motion of the molecular motors, and propagation and accumulation of this motion lead to contraction of the fibres towards the light source. This system supports large-amplitude motion, fast response, precise control over shape, as well as weight-lifting experiments in water and air.

Computation of Hemagglutinin Free Energy Difference by the Confinement Method.

Hemagglutinin (HA) mediates membrane fusion, a crucial step during influenza virus cell entry. How many HAs are needed for this process is still subject to debate. To aid in this discussion, the confinement free energy method was used to calculate the conformational free energy difference between the extended intermediate and postfusion state of HA. Special care was taken to comply with the general guidelines for free energy calculations, thereby obtaining convergence and demonstrating reliability of the results. The energy that one HA trimer contributes to fusion was found to be 34.2 ± 3.4kBT, similar to the known contributions from other fusion proteins. Although computationally expensive, the technique used is a promising tool for the further energetic characterization of fusion protein mechanisms. Knowledge of the energetic contributions per protein, and of conserved residues that are crucial for fusion, aids in the development of fusion inhibitors for antiviral drugs.

Mechanisms of influenza viral membrane fusion.

Influenza viral particles are enveloped by a lipid bilayer. A major step in infection is fusion of the viral and host cellular membranes, a process with large kinetic barriers. Influenza membrane fusion is catalyzed by hemagglutinin (HA), a class I viral fusion protein activated by low pH. The exact nature of the HA conformational changes that deliver the energy required for fusion remains poorly understood. This review summarizes our current knowledge of HA structure and dynamics, describes recent single-particle experiments and modeling studies, and discusses their role in understanding how multiple HAs mediate fusion. These approaches provide a mechanistic picture in which HAs independently and stochastically insert into the target membrane, forming a cluster of HAs that is collectively able to overcome the barrier to membrane fusion. The new experimental and modeling approaches described in this review hold promise for a more complete understanding of other viral fusion systems and the protein systems responsible for cellular fusion.

Energetics of Transport through the Nuclear Pore Complex.

Molecular transport across the nuclear envelope in eukaryotic cells is solely controlled by the nuclear pore complex (NPC). The NPC provides two types of nucleocytoplasmic transport: passive diffusion of small molecules and active chaperon-mediated translocation of large molecules. It has been shown that the interaction between intrinsically disordered proteins that line the central channel of the NPC and the transporting cargoes is the determining factor, but the exact mechanism of transport is yet unknown. Here, we use coarse-grained molecular dynamics simulations to quantify the energy barrier that has to be overcome for molecules to pass through the NPC. We focus on two aspects of transport. First, the passive transport of model cargo molecules with different sizes is studied and the size selectivity feature of the NPC is investigated. Our results show that the transport probability of cargoes is significantly reduced when they are larger than ∼5 nm in diameter. Secondly, we show that incorporating hydrophobic binding spots on the surface of the cargo effectively decreases the energy barrier of the pore. Finally, a simple transport model is proposed which characterizes the energy barrier of the NPC as a function of diameter and hydrophobicity of the transporting particles.

Two fundamental mechanisms govern the stiffening of cross-linked networks.

Biopolymer networks, such as those constituting the cytoskeleton of a cell or biological tissue, exhibit a nonlinear strain-stiffening behavior when subjected to large deformations. Interestingly, rheological experiments on various in vitro biopolymer networks have shown similar strain-stiffening trends regardless of the differences in their microstructure or constituents, suggesting a universal stiffening mechanism. In this article, we use computer simulations of a random network comprised of cross-linked biopolymer-like fibers to substantiate the notion that this universality lies in the existence of two fundamental stiffening mechanisms. After showing that the large strain response is accompanied by the development of a stress path, i.e., a percolating path of axially stressed fibers and cross-links, we demonstrate that the strain stiffening can be caused by two distinctly different mechanisms: 1) the pulling out of stress-path undulations; and 2) reorientation of the stress path. The former mechanism is bending-dominated and can be recognized by a power-law dependence with exponent 3/2 of the shear modulus on stress, whereas the latter mechanism is stretching-dominated and characterized by a power-law exponent 1/2. We demonstrate how material properties of the constituents, as well as the network microstructure, can affect the transition between the two stiffening mechanisms and, as such, control the dominant power-law scaling behavior.

Probing the disordered domain of the nuclear pore complex through coarse-grained molecular dynamics simulations.

The distribution of disordered proteins (FG-nups) that line the transport channel of the nuclear pore complex (NPC) is investigated by means of coarse-grained molecular dynamics simulations. A one-bead-per-amino-acid model is presented that accounts for the hydrophobic/hydrophilic and electrostatic interactions between different amino acids, polarity of the solvent, and screening of free ions. The results indicate that the interaction of the FG-nups forms a high-density, doughnut-like distribution inside the NPC, which is rich in FG-repeats. We show that the obtained distribution is encoded in the amino-acid sequence of the FG-nups and is driven by both electrostatic and hydrophobic interactions. To explore the relation between structure and function, we have systematically removed different combinations of FG-nups from the pore to simulate inviable and viable NPCs that were previously studied experimentally. The obtained density distributions show that the maximum density of the FG-nups inside the pore does not exceed 185 mg/mL in the inviable NPCs, whereas for the wild-type and viable NPCs, this value increases to 300 mg/mL. Interestingly, this maximum density is not correlated to the total mass of the FG-nups, but depends sensitively on the specific combination of essential Nups located in the central plane of the NPC.

Coarse-Grained Potentials for Local Interactions in Unfolded Proteins.

Recent studies have revealed the key role of natively unfolded proteins in many important biological processes. In order to study the conformational changes of these proteins, a one-bead-per-amino-acid coarse grained (CG) model is developed, and a method is proposed to extract the potential functions for the local interactions between CG beads. Experimentally obtained Ramachandran data for the coil regions of proteins are converted into distributions of pseudo-bond and pseudo-dihedral angles between neighboring alpha-carbons in the polypeptide chain. These are then used to derive bending and torsion potentials, which are residue and sequence specific. The validity of the developed model is testified by studying the radius of gyration as well as the hydrodynamic properties of chemically denatured proteins.

Release of content through mechano-sensitive gates in pressurized liposomes.

Mechano-sensitive channels are ubiquitous membrane proteins that activate in response to increasing tension in the lipid membrane. They facilitate a sudden, nonselective release of solutes and water that safeguards the integrity of the cell in hypo- or hyper-osmotic shock conditions. We have simulated the rapid release of content from a pressurized liposome through a particular mechano-sensitive protein channel, MscL, embedded in the liposomal membrane. We show that a single channel is able to relax the liposome, stressed to the point of bursting, in a matter of microseconds. We map the full activation-deactivation cycle of MscL in near-atomic detail and are able to quantify the rapid decrease in liposomal stress as a result of channel activation. This provides a computational tool that opens the way to contribute to the rational design of functional nano-containers.

Forced reptation revealed by chain pull-out simulations.

We report computation results obtained from extensive molecular dynamics simulations of tensile disentanglement of connector chains placed at the interface between two polymer bulks. Each polymer chain (either belonging to the bulks or being a connector) is treated as a sequence of beads interconnected by springs, using a coarse-grained representation based on the Kremer-Grest model, extended to account for stiffness along the chain backbone. Forced reptation of the connectors was observed during their disentanglement from the bulk chains. The extracted chains are clearly seen following an imaginary "tube" inside the bulks as they are pulled out. The entropic and energetic responses to the external deformation are investigated by monitoring the connector conformation tensor and the modifications of the internal parameters (bonds, bending, and torsion angles along the connectors). The work needed to separate the two bulks is computed from the tensile force induced during debonding in the connector chains. The value of the work reached at total separation is considered as the debonding energy G. The most important parameters controlling G are the length (n) of the chains placed at the interface and their areal density. Our in silico experiments are performed at relatively low areal density and are disregarded if chain scission occurs during disentanglement. As predicted by the reptation theory, for this pure pull-out regime, the power exponent from the scaling G proportional, variant n(a) is a approximately 2, irrespective of chain stiffness. Small variations are found when the connectors form different number of stitches at the interface, or when their length is randomly distributed in between the two bulks. Our results show that the effects of the number of stitches and of the randomness of the block lengths have to be considered together, especially when comparing with experiments where they cannot be controlled rigorously. These results may be significant for industrial applications, such reinforcement of polymer-polymer adhesion by connector chains, when incorporated as constitutive laws at higher time/length scales in finite element calculations.

Mechanosensitive membrane channels in action.

The tension-driven gating process of MscL from Mycobacterium tuberculosis, Tb-MscL, has been addressed at near-atomic detail using coarse-grained molecular dynamics simulations. To perform the simulations, a novel coarse-grained peptide model based on a thermodynamic parameterization of the amino-acid side chains has been applied. Both the wild-type Tb-MscL and its gain-of-function mutant V21D embedded in a solvated lipid bilayer have been studied. To mimic hypoosmotic shock conditions, simulations were performed at increasing levels of membrane tension approaching the rupture threshold of the lipid bilayer. Both the wild-type and the mutant channel are found to undergo significant conformational changes in accordance with an irislike expansion mechanism, reaching a conducting state on a microsecond timescale. The most pronounced expansion of the pore has been observed for the V21D mutant, which is consistent with the experimentally shown gain-of-function phenotype of the V21D mutant.