Plasticity of an ultrafast interaction between nucleoporins and nuclear transport receptors.

The mechanisms by which intrinsically disordered proteins engage in rapid and highly selective binding is a subject of considerable interest and represents a central paradigm to nuclear pore complex (NPC) function, where nuclear transport receptors (NTRs) move through the NPC by binding disordered phenylalanine-glycine-rich nucleoporins (FG-Nups). Combining single-molecule fluorescence, molecular simulations, and nuclear magnetic resonance, we show that a rapidly fluctuating FG-Nup populates an ensemble of conformations that are prone to bind NTRs with near diffusion-limited on rates, as shown by stopped-flow kinetic measurements. This is achieved using multiple, minimalistic, low-affinity binding motifs that are in rapid exchange when engaging with the NTR, allowing the FG-Nup to maintain an unexpectedly high plasticity in its bound state. We propose that these exceptional physical characteristics enable a rapid and specific transport mechanism in the physiological context, a notion supported by single molecule in-cell assays on intact NPCs.

Plant polygalacturonase structures specify enzyme dynamics and processivities to fine-tune cell wall pectins.

Polygalacturonases (PGs) fine-tune pectins to modulate cell wall chemistry and mechanics, impacting plant development. The large number of PGs encoded in plant genomes leads to questions on the diversity and specificity of distinct isozymes. Herein, we report the crystal structures of 2 Arabidopsis thaliana PGs, POLYGALACTURONASE LATERAL ROOT (PGLR), and ARABIDOPSIS DEHISCENCE ZONE POLYGALACTURONASE2 (ADPG2), which are coexpressed during root development. We first determined the amino acid variations and steric clashes that explain the absence of inhibition of the plant PGs by endogenous PG-inhibiting proteins (PGIPs). Although their beta helix folds are highly similar, PGLR and ADPG2 subsites in the substrate binding groove are occupied by divergent amino acids. By combining molecular dynamic simulations, analysis of enzyme kinetics, and hydrolysis products, we showed that these structural differences translated into distinct enzyme-substrate dynamics and enzyme processivities: ADPG2 showed greater substrate fluctuations with hydrolysis products, oligogalacturonides (OGs), with a degree of polymerization (DP) of ≤4, while the DP of OGs generated by PGLR was between 5 and 9. Using the Arabidopsis root as a developmental model, exogenous application of purified enzymes showed that the highly processive ADPG2 had major effects on both root cell elongation and cell adhesion. This work highlights the importance of PG processivity on pectin degradation regulating plant development.

Mutation of AtPME2, a pH-Dependent Pectin Methylesterase, Affects Cell Wall Structure and Hypocotyl Elongation.

Pectin methylesterases (PMEs) modify homogalacturonan's chemistry and play a key role in regulating primary cell wall mechanical properties. Here, we report on Arabidopsis AtPME2, which we found to be highly expressed during lateral root emergence and dark-grown hypocotyl elongation. We showed that dark-grown hypocotyl elongation was reduced in knock-out mutant lines as compared to the control. The latter was related to the decreased total PME activity as well as increased stiffness of the cell wall in the apical part of the hypocotyl. To relate phenotypic analyses to the biochemical specificity of the enzyme, we produced the mature active enzyme using heterologous expression in Pichia pastoris and characterized it through the use of a generic plant PME antiserum. AtPME2 is more active at neutral compared to acidic pH, on pectins with a degree of 55-70% methylesterification. We further showed that the mode of action of AtPME2 can vary according to pH, from high processivity (at pH8) to low processivity (at pH5), and relate these observations to the differences in electrostatic potential of the protein. Our study brings insights into how the pH-dependent regulation by PME activity could affect the pectin structure and associated cell wall mechanical properties.

Charge Relaying within a Phospho-Motif Rescue Binding Competency of a Disordered Transcription Factor.

Structural disorder in proteins is central to cellular signaling, where conformational plasticity equips molecules to promiscuously interact with different partners. By engaging with multiple binding partners via the rearrangement of its three helices, the nuclear coactivator binding domain (NCBD) of the CBP/p300 transcription factor is a paradigmatic example of promiscuity. Recently, molecular simulations and experiments revealed that, through the establishment of long-range electrostatic interactions, intended as salt-bridges formed between the post-translationally inserted phosphate and positively charged residues in helix H3 of NCBD, phosphorylation triggers NCBD compaction, lowering its affinity for binding partners. By means of extensive molecular simulations, we here investigated the effect of short-range electrostatics on the conformational ensemble of NCBD, by monitoring the interactions between a phosphorylated serine and conserved positively charged residues within the NCBD phospho-motif. We found that empowering proximal electrostatic interactions, as opposed to long-range electrostatics, can reshape the NCBD ensemble rescuing the binding competency of phosphorylated NCBD. Given the conservation of positive charges in phospho-motifs, proximal electrostatic interactions might dampen the effects of phosphorylation and act as a relay to regulate phosphorylated intrinsically disordered proteins, ultimately tuning the binding affinity for different cellular partners.

Mechanical force can enhance c-Src kinase activity by impairing autoinhibition.

Cellular mechanosensing is pivotal for virtually all biological processes, and many molecular mechano-sensors and their way of function are being uncovered. In this work, we suggest that c-Src kinase acts as a direct mechano-sensor. c-Src is responsible for, among others, cell proliferation, and shows increased activity in stretched cells. In its native state, c-Src has little basal activity, because its kinase domain binds to an SH2 and SH3 domain. However, it is known that c-Src can bind to p130Cas, through which force can be transmitted to the membrane. Using molecular dynamics simulations, we show that force acting between the membrane-bound N-terminus of the SH3 domain and p130Cas induces partial SH3 unfolding, thereby impeding rebinding of the kinase domain onto SH2/SH3 and effectively enhancing kinase activity. Forces involved in this process are slightly lower or similar to the forces required to pull out c-Src from the membrane through the myristoyl linker, and key interactions involved in this anchoring are salt bridges between negative lipids and nearby basic residues in c-Src. Thus, c-Src appears to be a candidate for an intriguing mechanosensing mechanism of impaired kinase inhibition, which can be potentially tuned by membrane composition and other environmental factors.

Multilayered allosteric modulation of coupled folding and binding by phosphorylation, peptidyl-prolyl cis/trans isomerization, and diversity of interaction partners.

Intrinsically disordered proteins (IDPs) play key roles in cellular regulation, including signal transduction, transcription, and cell-cycle control. Accordingly, IDPs can commonly interact with numerous different target proteins, and their interaction networks are expected to be highly regulated. However, many of the underlying regulatory mechanisms have remained unclear. Here, we examine the representative case of the nuclear coactivator binding domain (NCBD) of the large multidomain protein CBP, a hub in transcriptional regulation, and the interaction with several of its binding partners. Single-molecule Förster resonance energy transfer measurements show that phosphorylation of NCBD reduces its binding affinity, with effects that vary depending on the binding partner and the site and number of modifications. The complexity of the interaction is further increased by the dependence of the affinities on peptidyl-prolyl cis/trans isomerization in NCBD. Overall, our results reveal the potential for allosteric regulation on at least three levels: the different affinities of NCBD for its different binding partners, the differential modulation of these affinities by phosphorylation, and the effect of peptidyl-prolyl cis/trans isomerization on binding.

Structural and mechanistic insights into mechanoactivation of focal adhesion kinase.

Focal adhesion kinase (FAK) is a key signaling molecule regulating cell adhesion, migration, and survival. FAK localizes into focal adhesion complexes formed at the cytoplasmic side of cell attachment to the ECM and is activated after force generation via actomyosin fibers attached to this complex. The mechanism of translating mechanical force into a biochemical signal is not understood, and it is not clear whether FAK is activated directly by force or downstream to the force signal. We use experimental and computational single-molecule force spectroscopy to probe the mechanical properties of FAK and examine whether force can trigger activation by inducing conformational changes in FAK. By comparison with an open and active mutant of FAK, we are able to assign mechanoactivation to an initial rupture event in the low-force range. This activation event occurs before FAK unfolding at forces within the native range in focal adhesions. We are also able to assign all subsequent peaks in the force landscape to partial unfolding of FAK modules. We show that binding of ATP stabilizes the kinase domain, thereby altering the unfolding hierarchy. Using all-atom molecular dynamics simulations, we identify intermediates along the unfolding pathway, which provide buffering to allow extension of FAK in focal adhesions without compromising functionality. Our findings strongly support that forces in focal adhesions applied to FAK via known interactions can induce conformational changes, which in turn, trigger focal adhesion signaling.

Design of a Protein with Improved Thermal Stability by an Evolution-Based Generative Model.

Efficient design of functional proteins with higher thermal stability remains challenging especially for highly diverse sequence variants. Considering the evolutionary pressure on protein folds, sequence design optimizing evolutionary fitness could help designing folds with higher stability. Using a generative evolution fitness model trained to capture variation patterns in natural sequences, we designed artificial sequences of a proteinaceous inhibitor of pectin methylesterase enzymes. These inhibitors have considerable industrial interest to avoid phase separation in fruit juice manufacturing or reduce methanol in distillates, averting chromatographic passages triggering unwanted aroma loss. Six out of seven designs with up to 30 % divergence to other inhibitor sequences are functional and two have improved thermal stability. This method can improve protein stability expanding functional protein sequence space, with traits valuable for industrial applications and scientific research.

Probing the Paradigm of Promiscuity for N-Heterocyclic Carbene Complexes and their Protein Adduct Formation.

Metal complexes can be considered a "paradigm of promiscuity" when it comes to their interactions with proteins. They often form adducts with a variety of donor atoms in an unselective manner. We have characterized the adducts formed between a series of isostructural N-heterocyclic carbene (NHC) complexes with Ru, Os, Rh, and Ir centers and the model protein hen egg white lysozyme by X-ray crystallography and mass spectrometry. Distinctive behavior for the metal compounds was observed with the more labile Ru and Rh complexes targeting mainly a surface l-histidine moiety through cleavage of p-cymene or NHC co-ligands, respectively. In contrast, the more inert Os and Ir derivatives were detected abundantly in an electronegative binding pocket after undergoing ligand exchange of a chlorido ligand for an amino acid side chain. Computational studies supported the binding profiles and hinted at the role of the protein microenvironment for metal complexes eliciting selectivity for specific binding sites on the protein.

Decoupling of size and shape fluctuations in heteropolymeric sequences reconciles discrepancies in SAXS vs. FRET measurements.

Unfolded states of proteins and native states of intrinsically disordered proteins (IDPs) populate heterogeneous conformational ensembles in solution. The average sizes of these heterogeneous systems, quantified by the radius of gyration (RG ), can be measured by small-angle X-ray scattering (SAXS). Another parameter, the mean dye-to-dye distance (RE ) for proteins with fluorescently labeled termini, can be estimated using single-molecule Förster resonance energy transfer (smFRET). A number of studies have reported inconsistencies in inferences drawn from the two sets of measurements for the dimensions of unfolded proteins and IDPs in the absence of chemical denaturants. These differences are typically attributed to the influence of fluorescent labels used in smFRET and to the impact of high concentrations and averaging features of SAXS. By measuring the dimensions of a collection of labeled and unlabeled polypeptides using smFRET and SAXS, we directly assessed the contributions of dyes to the experimental values RG and RE For chemically denatured proteins we obtain mutual consistency in our inferences based on RG and RE , whereas for IDPs under native conditions, we find substantial deviations. Using computations, we show that discrepant inferences are neither due to methodological shortcomings of specific measurements nor due to artifacts of dyes. Instead, our analysis suggests that chemical heterogeneity in heteropolymeric systems leads to a decoupling between RE and RG that is amplified in the absence of denaturants. Therefore, joint assessments of RG and RE combined with measurements of polymer shapes should provide a consistent and complete picture of the underlying ensembles.

CONAN: A Tool to Decode Dynamical Information from Molecular Interaction Maps.

The analysis of contacts is a powerful tool to understand biomolecular function in a series of contexts, from the investigation of dynamical behavior at equilibrium to the study of nonequilibrium dynamics in which the system moves between multiple states. We thus propose a tool called CONtact ANalysis (CONAN) that, from molecular dynamics (MD) trajectories, analyzes interresidue contacts, creates videos of time-resolved contact maps, and performs correlation, principal component, and cluster analysis, revealing how specific contacts relate to functionally relevant states sampled by MD. We present how CONAN can identify features describing the dynamics of ubiquitin both at equilibrium and during mechanical unfolding. Additionally, we show the analysis of MD trajectories of an α-synuclein mutant peptide that undergoes an α-ß conformational transition that can be easily monitored using CONAN, which identifies the multiple states that the peptide explores along its conformational dynamics. The high versatility and ease of use of the software make CONAN a tool that can significantly facilitate the understanding of the complex dynamical behavior of proteins or other biomolecules. CONAN and its documentation are freely available for download on GitHub.

Structural and dynamical characterization of the pH-dependence of the pectin methylesterase-pectin methylesterase inhibitor complex.

Pectin methylesterases (PMEs) catalyze the demethylesterification of pectin, one of the main polysaccharides in the plant cell wall, and are of critical importance in plant development. PME activity generates highly negatively charged pectin and mutates the physiochemical properties of the plant cell wall such that remodeling of the plant cell can occur. PMEs are therefore tightly regulated by proteinaceous inhibitors (PMEIs), some of which become active upon changes in cellular pH. Nevertheless, a detailed picture of how this pH-dependent inhibition of PME occurs at the molecular level is missing. Herein, using an interdisciplinary approach that included homology modeling, MD simulations, and biophysical and biochemical characterizations, we investigated the molecular basis of PME3 inhibition by PMEI7 in Arabidopsis thaliana Our complementary approach uncovered how changes in the protonation of amino acids at the complex interface shift the network of interacting residues between intermolecular and intramolecular. These shifts ultimately regulate the stability of the PME3-PMEI7 complex and the inhibition of the PME as a function of the pH. These findings suggest a general model of how pH-dependent proteinaceous inhibitors function. Moreover, they enhance our understanding of how PMEs may be regulated by pH and provide new insights into how this regulation may control the physical properties and structure of the plant cell wall.

Combined Experimental and Computational Approaches Reveal Distinct pH Dependence of Pectin Methylesterase Inhibitors.

The fine-tuning of the degree of methylesterification of cell wall pectin is a key to regulating cell elongation and ultimately the shape of the plant body. Pectin methylesterification is spatiotemporally controlled by pectin methylesterases (PMEs; 66 members in Arabidopsis [Arabidopsis thaliana]). The comparably large number of proteinaceous pectin methylesterase inhibitors (PMEIs; 76 members in Arabidopsis) questions the specificity of the PME-PMEI interaction and the functional role of such abundance. To understand the difference, or redundancy, between PMEIs, we used molecular dynamics (MD) simulations to predict the behavior of two PMEIs that are coexpressed and have distinct effects on plant development: AtPMEI4 and AtPMEI9. Simulations revealed the structural determinants of the pH dependence for the interaction of these inhibitors with AtPME3, a major PME expressed in roots. Key residues that are likely to play a role in the pH dependence were identified. The predictions obtained from MD simulations were confirmed in vitro, showing that AtPMEI9 is a stronger, less pH-independent inhibitor compared with AtPMEI4. Using pollen tubes as a developmental model, we showed that these biochemical differences have a biological significance. Application of purified proteins at pH ranges in which PMEI inhibition differed between AtPMEI4 and AtPMEI9 had distinct consequences on pollen tube elongation. Therefore, MD simulations have proven to be a powerful tool to predict functional diversity between PMEIs, allowing the discovery of a strategy that may be used by PMEIs to inhibit PMEs in different microenvironmental conditions and paving the way to identify the specific role of PMEI diversity in muro.

On the electrophoretic mobilities of partially charged oligosaccharides as a function of charge patterning and degree of polymerization.

Fully or partially charged oligosaccharide molecules play a key role in many areas of biology, where their fine structures are crucial in determining their functionality. However, the separation of specific charged oligosaccharides from similar moieties that typically coexist in extracted samples, even for those that are unbranched, and in cases where each saccharide moiety can only carry a single charge or not, is far from trivial. Typically such molecules are characterized by a degree of polymerization n and a number m (and distribution) of charged residues, and must be separated from a plethora of similar species possessing different combinations of n and m. Furthermore, the separation of the possible n!/m!(n-m)! isomers of each species of fixed n and m is a formidable challenge to analytical chemists. Herein, we report the results of molecular dynamics simulations that have been performed in order to calculate the free solution electrophoretic mobilities of galacturonides and charged oligosaccharides derived from digests of the important plant cell-wall polysaccharide pectin. The simulations are compared with an experiment and are found to correctly predict the loss of resolution of fully charged species above a critical degree of polymerization n and the ionic strength dependence of the electrophoretic mobilities of different partially charged oligosaccharides. It is expected that having a predictive tool for the calculation of the electrophoretic mobilities of differently charged oligosaccharide species in hand will allow experimental conditions that optimize the resolution of particular species to be ascertained and understood.

Structure and Properties of a Non-processive, Salt-requiring, and Acidophilic Pectin Methylesterase from Aspergillus niger Provide Insights into the Key Determinants of Processivity Control.

Many pectin methylesterases (PMEs) are expressed in plants to modify plant cell-wall pectins for various physiological roles. These pectins are also attacked by PMEs from phytopathogens and phytophagous insects. The de-methylesterification by PMEs of the O6-methyl ester groups of the homogalacturonan component of pectin, exposing galacturonic acids, can occur processively or non-processively, respectively, describing sequential versus single de-methylesterification events occurring before enzyme-substrate dissociation. The high resolution x-ray structures of a PME from Aspergillus niger in deglycosylated and Asn-linked N-acetylglucosamine-stub forms reveal a 10⅔-turn parallel ß-helix (similar to but with less extensive loops than bacterial, plant, and insect PMEs). Capillary electrophoresis shows that this PME is non-processive, halophilic, and acidophilic. Molecular dynamics simulations and electrostatic potential calculations reveal very different behavior and properties compared with processive PMEs. Specifically, uncorrelated rotations are observed about the glycosidic bonds of a partially de-methyl-esterified decasaccharide model substrate, in sharp contrast to the correlated rotations of processive PMEs, and the substrate-binding groove is negatively not positively charged.

The Cell Wall Arabinose-Deficient Arabidopsis thaliana Mutant murus5 Encodes a Defective Allele of REVERSIBLY GLYCOSYLATED POLYPEPTIDE2.

Traditional marker-based mapping and next-generation sequencing was used to determine that the Arabidopsis (Arabidopsis thaliana) low cell wall arabinose mutant murus5 (mur5) encodes a defective allele of REVERSIBLY GLYCOSYLATED POLYPEPTIDE2 (RGP2). Marker analysis of 13 F2 confirmed mutant progeny from a recombinant mapping population gave a rough map position on the upper arm of chromosome 5, and deep sequencing of DNA from these 13 lines gave five candidate genes with G→A (C→T) transitions predicted to result in amino acid changes. Of these five, only insertional mutant alleles of RGP2, a gene that encodes a UDP-arabinose mutase that interconverts UDP-arabinopyranose and UDP-arabinofuranose, exhibited the low cell wall arabinose phenotype. The identities of mur5 and two SALK insertional alleles were confirmed by allelism tests and overexpression of wild-type RGP2 complementary DNA placed under the control of the 35S promoter in the three alleles. The mur5 mutation results in the conversion of cysteine-257 to tyrosine-257 within a conserved hydrophobic cluster predicted to be distal to the active site and essential for protein stability and possible heterodimerization with other isoforms of RGP.

Molecular Dynamics Simulations Illuminate the Role of Counterion Condensation in the Electrophoretic Transport of Homogalacturonans.

Homogalacturonans (HGs) are polysaccharide copolymers of galacturonic acid and its methylesterified counterpart. The inter- and intramolecular distributions of the methylesterifed residues are vital behavior-determining characteristics of a sample's structure, and much experimental effort has been directed to their measurement. While many techniques are able to measure the sample-averaged degree of methylesterification (DM), the measurement of inter- and intramolecular charge distributions are challenging. Here, molecular dynamics (MD) simulations are used to calculate the electrophoretic mobilities of HGs that have different amounts and distributions of charges placed along the backbone. The simulations are shown to capture experimental results well, even for low-DM samples that possess high charge densities. In addition, they illuminate the role that local counterion condensation can play in the determination of the electrophoretic mobility of heterogeneous blocky polyelectrolytes that cannot be adequately described by a single chain-averaged charge spacing.

A fast recoiling silk-like elastomer facilitates nanosecond nematocyst discharge.

BACKGROUND: The discharge of the Cnidarian stinging organelle, the nematocyst, is one of the fastest processes in biology and involves volume changes of the highly pressurised (150 bar) capsule of up to 50%. Hitherto, the molecular basis for the unusual biomechanical properties of nematocysts has been elusive, as their structure was mainly defined as a stress-resistant collagenous matrix. RESULTS: Here, we characterise Cnidoin, a novel elastic protein identified as a structural component of Hydra nematocysts. Cnidoin is expressed in nematocytes of all types and immunostainings revealed incorporation into capsule walls and tubules concomitant with minicollagens. Similar to spider silk proteins, to which it is related at sequence level, Cnidoin possesses high elasticity and fast coiling propensity as predicted by molecular dynamics simulations and quantified by force spectroscopy. Recombinant Cnidoin showed a high tendency for spontaneous aggregation to bundles of fibrillar structures. CONCLUSIONS: Cnidoin represents the molecular factor involved in kinetic energy storage and release during the ultra-fast nematocyst discharge. Furthermore, it implies an early evolutionary origin of protein elastomers in basal metazoans.

Origin of Orthogonality of Strain-Promoted Click Reactions.

Site-specific labeling of biomolecules is rapidly advancing due to the discovery of novel mutually orthogonal reactions. Quantum chemistry studies have also increased our understanding of their relative rates, although these have until now been based on highly simplified reactants. Here we examine a set of strain-promoted click-type cycloaddition reactions of n-propyl azide, 3-benzyl tetrazine and 3-benzyl-6-methyl tetrazine with cyclooctenes/ynes, in which we aim to address all relevant structural details of the reactants. Our calculations have included the obligatory handles used to attach the label and biomolecule as these can critically influence the stereochemistry and electron demand of the reaction. We systematically computed orbital gaps, activation and distortion energies using density functional theory and determined experimental rates for validation. Our results challenge the current paradigm of the inverse electron demand for this class of reactions. We found that the ubiquitous handles, when next to the triple bond of cyclooctynes, can switch the Diels-Alder type ligations to normal electron demand, a class we term as SPINEDAC reactions. Electron donating substituents on tetrazine can enhance normal demand but also increase distortion penalties. The presence and isomeric configuration of handles thus determine the reaction speed and regioselectivity. Our findings can be directly utilized in engineering genuine cycloaddition click chemistries for biological labeling.

Disordered regions of human eIF4B orchestrate a dynamic self-association landscape.

Eukaryotic translation initiation factor eIF4B is required for efficient cap-dependent translation, it is overexpressed in cancer cells, and may influence stress granule formation. Due to the high degree of intrinsic disorder, eIF4B is rarely observed in cryo-EM structures of translation complexes and only ever by its single structured RNA recognition motif domain, leaving the molecular details of its large intrinsically disordered region (IDR) unknown. By integrating experiments and simulations we demonstrate that eIF4B IDR orchestrates and fine-tunes an intricate transition from monomers to a condensed phase, in which large-size dynamic oligomers form before mesoscopic phase separation. Single-molecule spectroscopy combined with molecular simulations enabled us to characterize the conformational ensembles and underlying intra- and intermolecular dynamics across the oligomerization transition. The observed sensitivity to ionic strength and molecular crowding in the self-association landscape suggests potential regulation of eIF4B nanoscopic and mesoscopic behaviors such as driven by protein modifications, binding partners or changes to the cellular environment.