MD simulations and FRET reveal an environment-sensitive conformational plasticity of importin-β

Influence of point mutations on PR65 conformational adaptability: Insights from molecular simulations and nanoaperture optical tweezers

PR65 is the HEAT repeat scaffold subunit of the heterotrimeric protein phosphatase 2A (PP2A) and an archetypal tandem repeat protein. Its conformational mechanics plays a crucial role in PP2A function by opening/closing substrate binding/catalysis interface. Using in silico saturation mutagenesis, we identified PR65 "hinge" residues whose substitutions could alter its conformational adaptability and thereby PP2A function, and selected six mutations that were verified to be expressed and soluble. Molecular simulations and nanoaperture optical tweezers revealed consistent results on the specific effects of the mutations on the structure and dynamics of PR65. Two mutants observed in simulations to stabilize extended/open conformations exhibited higher corner frequencies and lower translational scattering in experiments, indicating a shift toward extended conformations, whereas another displayed the opposite features, confirmed by both simulations and experiments. The study highlights the power of single-molecule nanoaperture-based tweezers integrated with in silico approaches for exploring the effect of mutations on protein structure and dynamics.

Tandem-repeat proteins conformational mechanics are optimized to facilitate functional interactions and complexations

The architectures of tandem-repeat proteins are distinct from those of globular proteins. Individual modules, each comprising small structural motifs of 20-40 residues, are arrayed in a quasi one-dimensional fashion to form striking, elongated, horseshoe-like, and superhelical architectures, stabilized solely by short-range interaction. The spring-like shapes of repeat arrays point to elastic modes of action, and these proteins function as adapter molecules or 'hubs,' propagating signals within multi-subunit assemblies in diverse biological contexts. This flexibility is apparent in the dramatic variability observed in the structures of tandem-repeat proteins in different complexes. Here, using computational analysis, we demonstrate the striking ability of just one or a few global motions to recapitulate these structures. These findings show how the mechanics of repeat arrays are robustly enabled by their unique architecture. Thus, the repeating architecture has been optimized by evolution to favor functional modes of motions. The global motions enabling functional transitions can be fully visualized at http://bahargroup.org/tr_web.

Influence of Point Mutations on PR65 Conformational Adaptability: Insights from Nanoaperture Optical Tweezer Experiments and Molecular Simulations

PR65 is the HEAT-repeat scaffold subunit of the heterotrimeric protein phosphatase 2A (PP2A) and an archetypal tandem-repeat protein, forming a spring-like architecture. PR65 conformational mechanics play a crucial role in PP2A function by opening/closing the substrate-binding/catalysis interface. Using in-silico saturation mutagenesis we identified "hinge" residues of PR65, whose substitutions are predicted to restrict its conformational adaptability and thereby disrupt PP2A function. Molecular simulations revealed that a subset of hinge mutations stabilized the extended/open conformation, whereas another had the opposite effect. By trapping in nanoaperture optical tweezer, we characterized PR65 motion and showed that the former mutants exhibited higher corner frequencies and lower translational scattering, indicating a shift towards extended conformations, whereas the latter showed the opposite behavior. Thus, experiments confirm the conformations predicted computationally. The study highlights the utility of nanoaperture-based tweezers for exploring structure and dynamics, and the power of integrating this single-molecule method with in silico approaches.

Single-photon smFRET. I: Theory and conceptual basis

We present a unified conceptual framework and the associated software package for single-molecule Förster resonance energy transfer (smFRET) analysis from single-photon arrivals leveraging Bayesian nonparametrics, BNP-FRET. This unified framework addresses the following key physical complexities of a single-photon smFRET experiment, including: 1) fluorophore photophysics; 2) continuous time kinetics of the labeled system with large timescale separations between photophysical phenomena such as excited photophysical state lifetimes and events such as transition between system states; 3) unavoidable detector artefacts; 4) background emissions; 5) unknown number of system states; and 6) both continuous and pulsed illumination. These physical features necessarily demand a novel framework that extends beyond existing tools. In particular, the theory naturally brings us to a hidden Markov model with a second-order structure and Bayesian nonparametrics on account of items 1, 2, and 5 on the list. In the second and third companion articles, we discuss the direct effects of these key complexities on the inference of parameters for continuous and pulsed illumination, respectively.

Binding stoichiometry and structural model of the HIV-1 Rev/importin β complex

HIV-1 Rev mediates the nuclear export of intron-containing viral RNA transcripts and is essential for viral replication. Rev is imported into the nucleus by the host protein importin β (Impβ), but how Rev associates with Impβ is poorly understood. Here, we report biochemical, mutational, and biophysical studies of the Impβ/Rev complex. We show that Impβ binds two Rev monomers through independent binding sites, in contrast to the 1:1 binding stoichiometry observed for most Impβ cargos. Peptide scanning data and charge-reversal mutations identify the N-terminal tip of Rev helix α2 within Rev's arginine-rich motif (ARM) as a primary Impβ-binding epitope. Cross-linking mass spectrometry and compensatory mutagenesis data combined with molecular docking simulations suggest a structural model in which one Rev monomer binds to the C-terminal half of Impβ with Rev helix α2 roughly parallel to the HEAT-repeat superhelical axis, whereas the other monomer binds to the N-terminal half. These findings shed light on the molecular basis of Rev recognition by Impβ and highlight an atypical binding behavior that distinguishes Rev from canonical cellular Impβ cargos.

A topological data analytic approach for discovering biophysical signatures in protein dynamics

Identifying structural differences among proteins can be a non-trivial task. When contrasting ensembles of protein structures obtained from molecular dynamics simulations, biologically-relevant features can be easily overshadowed by spurious fluctuations. Here, we present SINATRA Pro, a computational pipeline designed to robustly identify topological differences between two sets of protein structures. Algorithmically, SINATRA Pro works by first taking in the 3D atomic coordinates for each protein snapshot and summarizing them according to their underlying topology. Statistically significant topological features are then projected back onto a user-selected representative protein structure, thus facilitating the visual identification of biophysical signatures of different protein ensembles. We assess the ability of SINATRA Pro to detect minute conformational changes in five independent protein systems of varying complexities. In all test cases, SINATRA Pro identifies known structural features that have been validated by previous experimental and computational studies, as well as novel features that are also likely to be biologically-relevant according to the literature. These results highlight SINATRA Pro as a promising method for facilitating the non-trivial task of pattern recognition in trajectories resulting from molecular dynamics simulations, with substantially increased resolution.

Structural features of proteins often serve as signatures of their biological function and molecular binding activity. Elucidating these structural features is essential for a full understanding of underlying biophysical mechanisms. While there are existing methods aimed at identifying structural differences between protein variants, such methods do not have the capability to jointly infer both geometric and dynamic changes, simultaneously. In this paper, we propose SINATRA Pro, a computational framework for extracting key structural features between two sets of proteins. SINATRA Pro robustly outperforms standard techniques in pinpointing the physical locations of both static and dynamic signatures across various types of protein ensembles, and it does so with improved resolution.

Impaired nuclear transport induced by juvenile ALS causing P525L mutation in NLS domain of FUS: A molecular mechanistic study

Amyotrophic lateral sclerosis (ALS) and fronto-temporal lobar degeneration (FTLD) are progressive neurological disorders affecting motor neurons. Cellular aggregates of fused in sarcoma (FUS) protein are found in cytoplasm of ALS and FTLD patients. Nuclear localisation signal (NLS) domain of FUS binds to Karyopherin β2 (Kapβ2), which drives nuclear transport of FUS from cytoplasm. Several pathogenic mutations are reported in FUS NLS, which are associated with its impaired nuclear transport and cytoplasmic mis-localisation. P525L mutation in NLS is most commonly found in cases of juvenile ALS (jALS), which affects individuals below 25 years of age. jALS progresses aggressively causing death within a year of its onset. This study elucidates the molecular mechanism behind jALS-causing P525L mutation hindering nuclear transport of FUS. We perform multiple molecular dynamics simulations in aqueous and hydrophobic solvent to understand the effect of the mutation at molecular level. Dynamics of Kapβ2-FUS complex is better captured in hydrophobic solvent compared to aqueous solvent. P525 and Y526 (PY-motif) of NLS exhibit fine-tuned stereochemical arrangement, which is essential for optimum Kapβ2 binding. P525L causes loss of several native contacts at interface leading to weaker binding, which promotes self-aggregation of FUS in cytoplasm. Native complex samples closed conformation, while mutant complex exhibits open conformation exposing hydrophilic residues of Kapβ2 to hydrophobic solvent. Mutant complex also fails to exhibit spring-like motion essential for its transport through nuclear pore complex. This study provides a mechanistic insight of binding affinity between NLS and Kapβ2 that inhibits self-aggregation of FUS preventing the disease condition.

Comparative study of the interaction of ivermectin with proteins of interest associated with SARS-CoV-2: A computational and biophysical approach

The SARS-CoV-2 pandemic has accelerated the study of existing drugs. The mixture of homologs called ivermectin (avermectin-B1a [HB1a] + avermectin-B1b [HB1b]) has shown antiviral activity against SARS-CoV-2 in vitro. However, there are few reports on the behavior of each homolog. We investigated the interaction of each homolog with promising targets of interest associated with SARS-CoV-2 infection from a biophysical and computational-chemistry perspective using docking and molecular dynamics. We observed a differential behavior for each homolog, with an affinity of HB1b for viral structures, and of HB1a for host structures considered. The induced disturbances were differential and influenced by the hydrophobicity of each homolog and of the binding pockets. We present the first comparative analysis of the potential theoretical inhibitory effect of both avermectins on biomolecules associated with COVID-19, and suggest that ivermectin through its homologs, has a multiobjective behavior.

Structural deformability induced in proteins of potential interest associated with COVID-19 by binding of homologues present in ivermectin: Comparative study based in elastic networks models

The COVID-19 pandemic has accelerated the study of the potential of multi-target drugs (MTDs). The mixture of homologues called ivermectin (avermectin-B1a + avermectin-B1b) has been shown to be a MTD with potential antiviral activity against SARS-CoV-2 in vitro. However, there are few reports on the effect of each homologue on the flexibility and stiffness of proteins associated with COVID-19, described as ivermectin targets. We observed that each homologue was stably bound to the proteins studied and was able to induce detectable changes with Elastic Network Models (ENM). The perturbations induced by each homologue were characteristic of each compound and, in turn, were represented by a disruption of native intramolecular networks (interactions between residues). The homologues were able to slightly modify the conformation and stability of the connection points between the Cα atoms of the residues that make up the structural network of proteins (nodes), compared to free proteins. Each homologue was able to modified differently the distribution of quasi-rigid regions of the proteins, which could theoretically alter their biological activities. These results could provide a biophysical-computational view of the potential MTD mechanism that has been reported for ivermectin.

The Polyglutamine Expansion at the N-Terminal of Huntingtin Protein Modulates the Dynamic Configuration and Phosphorylation of the C-Terminal HEAT Domain

The polyQ expansion in huntingtin protein (HTT) is the prime cause of Huntington's disease (HD). The recent cryoelectron microscopy (cryo-EM) structure of HTT-HAP40 complex provided the structural information on its HEAT-repeat domains. Here, we present analyses of the impact of polyQ length on the structure and function of HTT via an integrative structural and biochemical approach. The cryo-EM analysis of normal (Q23) and disease (Q78) type HTTs shows that the structures of apo HTTs significantly differ from the structure of HTT in a HAP40 complex and that the polyQ expansion induces global structural changes in the relative movements among the HTT domains. In addition, we show that the polyQ expansion alters the phosphorylation pattern across HTT and that Ser2116 phosphorylation in turn affects the global structure and function of HTT. These results provide a molecular basis for the effect of the polyQ segment on HTT structure and activity, which may be important for HTT pathology.

Large-Scale Conformational Changes and Protein Function: Breaking the in silico Barrier

Large-scale conformational changes are essential to link protein structures with their function at the cell and organism scale, but have been elusive both experimentally and computationally. Over the past few years developments in cryo-electron microscopy and crystallography techniques have started to reveal multiple snapshots of increasingly large and flexible systems, deemed impossible only short time ago. As structural information accumulates, theoretical methods become central to understand how different conformers interconvert to mediate biological function. Here we briefly survey current in silico methods to tackle large conformational changes, reviewing recent examples of cross-validation of experiments and computational predictions, which show how the integration of different scale simulations with biological information is already starting to break the barriers between the in silico, in vitro, and in vivo worlds, shedding new light onto complex biological problems inaccessible so far.

Non-canonical amino acid labeling in proteomics and biotechnology

Metabolic labeling of proteins with non-canonical amino acids (ncAAs) provides unique bioorthogonal chemical groups during de novo synthesis by taking advantage of both endogenous and heterologous protein synthesis machineries. Labeled proteins can then be selectively conjugated to fluorophores, affinity reagents, peptides, polymers, nanoparticles or surfaces for a wide variety of downstream applications in proteomics and biotechnology. In this review, we focus on techniques in which proteins are residue- and site-specifically labeled with ncAAs containing bioorthogonal handles. These ncAA-labeled proteins are: readily enriched from cells and tissues for identification via mass spectrometry-based proteomic analysis; selectively purified for downstream biotechnology applications; or labeled with fluorophores for in situ analysis. To facilitate the wider use of these techniques, we provide decision trees to help guide the design of future experiments. It is expected that the use of ncAA labeling will continue to expand into new application areas where spatial and temporal analysis of proteome dynamics and engineering new chemistries and new function into proteins are desired.

Karyopherins regulate nuclear pore complex barrier and transport function

Kapinos et al. show that nuclear pore complex permeability and cargo release functionalities are concomitantly regulated by karyopherin occupancy and turnover in a systematic continuum. This highlights increasingly important roles for the soluble nucleocytoplasmic transport machinery that depart from established views of the nuclear pore complex selectivity mechanism.

Nucleocytoplasmic transport is sustained by karyopherins (Kaps) and a Ran guanosine triphosphate (RanGTP) gradient that imports nuclear localization signal (NLS)–specific cargoes (NLS-cargoes) into the nucleus. However, how nuclear pore complex (NPC) barrier selectivity, Kap traffic, and NLS-cargo release are systematically linked and simultaneously regulated remains incoherent. In this study, we show that Kapα facilitates Kapβ1 turnover and occupancy at the NPC in a RanGTP-dependent manner that is directly coupled to NLS-cargo release and NPC barrier function. This is underpinned by the binding affinity of Kapβ1 to phenylalanine–glycine nucleoporins (FG Nups), which is comparable with RanGTP·Kapβ1, but stronger for Kapα·Kapβ1. On this basis, RanGTP is ineffective at releasing standalone Kapβ1 from NPCs. Depleting Kapα·Kapβ1 by RanGTP further abrogates NPC barrier function, whereas adding back Kapβ1 rescues it while Kapβ1 turnover softens it. Therefore, the FG Nups are necessary but insufficient for NPC barrier function. We conclude that Kaps constitute integral constituents of the NPC whose barrier, transport, and cargo release functionalities establish a continuum under a mechanism of Kap-centric control.

In vivo analysis of protein crowding within the nuclear pore complex in interphase and mitosis

The central channel of the nuclear pore complex (NPC) is occupied by non-structured polypeptides with a high content of Phe-Gly (FG) motifs. This protein-rich environment functions as an entropic barrier that prevents the passage of molecules, as well as the binding sites for karyopherins, to regulate macromolecular traffic between the nucleoplasm and the cytoplasm. In this study, we expressed individual Nups fused with a crowding-sensitive probe (GimRET) to determine the spatial distribution of protein-rich domains within the central channel in vivo, and characterize the properties of the entropic barrier. Analyses of the probe signal revealed that the central channel contains two protein-rich domains at both the nucleoplasmic and cytoplasmic peripheries, and a less-crowded central cavity. Karyopherins and other soluble proteins are not the constituents of the protein-rich domains. The time-lapse observation of the post-mitotic reassembly process also revealed how individual protein-rich domains are constructed by a sequential assembly of nucleoporins.

How to operate a nuclear pore complex by Kap-centric control

Nuclear pore complexes (NPCs) mediate molecular transport between the nucleus and cytoplasm in eukaryotic cells. Tethered within each NPC lie numerous intrinsically disordered proteins known as FG nucleoporins (FG Nups) that are central to this process. Over two decades of investigation has converged on a view that a barrier mechanism consisting of FG Nups rejects non-specific macromolecules while promoting the speed and selectivity of karyopherin (Kaps) receptors (and their cargoes). Yet, the number of NPCs in the cell is exceedingly small compared to the number of Kaps, so that in fact there is a high likelihood the pores are always populated by Kaps. Here, we contemplate a view where Kaps actively participate in regulating the selectivity and speed of transport through NPCs. This so-called "Kap-centric" control of the NPC accounts for Kaps as essential barrier reinforcements that play a prerequisite role in facilitating fast transport kinetics. Importantly, Kap-centric control reconciles both mechanistic and kinetic requirements of the NPC, and in so doing potentially resolves incoherent aspects of FG-centric models. On this basis, we surmise that Kaps prime the NPC for nucleocytoplasmic transport by fine-tuning the NPC microenvironment according to the functional needs of the cell.

Impact of the crystallization condition on importin-β conformation

In eukaryotic cells, the exchange of macromolecules between the nucleus and cytoplasm is highly selective and requires specialized soluble transport factors. Many of them belong to the importin-β superfamily, the members of which share an overall superhelical structure owing to the tandem arrangement of a specific motif, the HEAT repeat. This structural organization leads to great intrinsic flexibility, which in turn is a prerequisite for the interaction with a variety of proteins and for its transport function. During the passage from the aqueous cytosol into the nucleus, the receptor passes the gated channel of the nuclear pore complex filled with a protein meshwork of unknown organization, which seems to be highly selective owing to the presence of FG-repeats, which are peptides with hydrophobic patches. Here, the structural changes of free importin-β from a single organism, crystallized in polar (salt) or apolar (PEG) buffer conditions, are reported. This allowed analysis of the structural changes, which are attributable to the surrounding milieu and are not affected by bound interaction partners. The importin-β structures obtained exhibit significant conformational changes and suggest an influence of the polarity of the environment, resulting in an extended conformation in the PEG condition. The significance of this observation is supported by SAXS experiments and the analysis of other crystal structures of importin-β deposited in the Protein Data Bank.