Hydrodynamic radii of native and denatured proteins measured by pulse field gradient NMR techniques

Identifying the region responsible for Brucella abortus MucR higher-order oligomer formation and examining its role in gene regulation

MucR is a member of the Ros/MucR family of prokaryotic zinc-finger proteins found in the α-proteobacteria which regulate the expression of genes required for the successful pathogenic and symbiotic interactions of these bacteria with the eukaryotic hosts. The structure and function of their distinctive zinc-finger domain has been well-studied, but only recently the quaternary structure of the full length proteins was investigated demonstrating their ability to form higher-order oligomers. The aim of this study was to identify the region of MucR involved in higher-order oligomer formation by analysing deletion and point mutants of this protein by Light Scattering, and to determine the role that MucR oligomerization plays in the regulatory function of this protein. Here we demonstrate that a conserved hydrophobic region at the N-terminus of MucR is responsible for higher-order oligomer formation and that MucR oligomerization is essential for its regulatory function in Brucella. All these features of MucR are shared by the histone-like nucleoid structuring protein, (H-NS), leading us to propose that the prokaryotic zinc-finger proteins in the MucR/Ros family control gene expression employing a mechanism similar to that used by the H-NS proteins, rather than working as classical transcriptional regulators.

The 2.1 Å structure of protein F9 and its comparison to L1, two components of the conserved poxvirus entry-fusion complex

The poxvirus F9 protein is a component of the vaccinia virus entry fusion complex (EFC) which consists of 11 proteins. The EFC forms a unique apparatus among viral fusion proteins and complexes. We solved the atomic structure of the F9 ectodomain at 2.10 Å. A structural comparison to the ectodomain of the EFC protein L1 indicated a similar fold and organization, in which a bundle of five α-helices is packed against two pairs of β-strands. However, instead of the L1 myristoylation site and hydrophobic cavity, F9 possesses a protruding loop between α-helices α3 and α4 starting at Gly90. Gly90 is conserved in all poxviruses except Salmon gill poxvirus (SGPV) and Diachasmimorpha longicaudata entomopoxvirus. Phylogenetic sequence analysis of all Poxviridae F9 and L1 orthologs revealed the SGPV genome to contain the most distantly related F9 and L1 sequences compared to the vaccinia proteins studied here. The structural differences between F9 and L1 suggest functional adaptations during evolution from a common precursor that underlie the present requirement for each protein.

Probing Polymer Chain Conformation and Fibril Formation of Peptide Conjugates

Covalent conjugates between a synthetic polymer and a peptide hormone were used to probe the molecular extension of these macromolecules and how the polymer modifies the fibril formation of the hormone. NMR spectroscopy of ¹⁵ N labeled parathyroid hormone (PTH) was employed to visualize the conformation of the conjugated synthetic polymer, triggered by small temperature changes via its lower critical solution temperature. A shroud-like polymer conformation dominated the molecular architecture of the conjugated chimeras. PTH readily forms amyloid fibrils, which is probably the physiological storage form of the hormone. The polyacrylate based polymers stimulated the nucleation processes of the peptide.

Characterization of Dynamic IDP Complexes by NMR Spectroscopy

NMR spectroscopy has proven to be a key method for studying intrinsically disordered proteins (IDPs). Nonetheless, traditional NMR methods developed for solving structures of ordered protein complexes are insufficient for the full characterization of dynamic IDP complexes, where the energy landscape is broader and more rugged. Furthermore, due to their high sensitivity to environmental changes, NMR studies of IDP complexes must be conducted with extra care and the observed NMR parameters thoroughly evaluated to enable disentanglement of binding events from ensemble distribution changes. In this chapter, written for the non-NMR expert, we start out by outlining sample preparation for IDP complexes, guide through the recording and evaluation of diagnostic ¹H,¹⁵N-HSQC spectra, and delineate more sophisticated NMR strategies to follow for the particular type of complex. The most relevant experiments are then described in terms of aims, needs, pitfalls, analysis, and expected outcomes, with references to recent examples.

Hierarchical spidroin micellar nanoparticles as the fundamental precursors of spider silks

Many natural silks produced by spiders and insects are unique materials in their exceptional toughness and tensile strength, while being lightweight and biodegradable-properties that are currently unparalleled in synthetic materials. Myriad approaches have been attempted to prepare artificial silks from recombinant spider silk spidroins but have each failed to achieve the advantageous properties of the natural material. This is because of an incomplete understanding of the in vivo spidroin-to-fiber spinning process and, particularly, because of a lack of knowledge of the true morphological nature of spidroin nanostructures in the precursor dope solution and the mechanisms by which these nanostructures transform into micrometer-scale silk fibers. Herein we determine the physical form of the natural spidroin precursor nanostructures stored within spider glands that seed the formation of their silks and reveal the fundamental structural transformations that occur during the initial stages of extrusion en route to fiber formation. Using a combination of solution phase diffusion NMR and cryogenic transmission electron microscopy (cryo-TEM), we reveal direct evidence that the concentrated spidroin proteins are stored in the silk glands of black widow spiders as complex, hierarchical nanoassemblies (∼300 nm diameter) that are composed of micellar subdomains, substructures that themselves are engaged in the initial nanoscale transformations that occur in response to shear. We find that the established micelle theory of silk fiber precursor storage is incomplete and that the first steps toward liquid crystalline organization during silk spinning involve the fibrillization of nanoscale hierarchical micelle subdomains.

Local and non-local topological information in the denatured state ensemble of a β-barrel protein

The folding of predominantly β-sheet proteins is complicated by the presence of a large number of non-local interactions in their native states, which increase the ruggedness of their folding energy landscapes. However, forming non-local contacts early in folding or even in the unfolded state can smooth the energy landscape and facilitate productive folding. We report that several sequence regions of a β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), populate native-like secondary structure to a significant extent in the denatured state in 8 M urea. In addition, we provide evidence for both local and non-local interactions in the denatured state of CRABP1. NMR chemical shift perturbations (CSPs) under denaturing conditions upon substitution of single residues by mutation support the presence of several non-local interactions in topologically key sites, arguing that the denatured state is conformationally restricted and contains topological information for the native fold. Among the most striking non-local interactions are those between the N- and C-terminal regions, which are involved in closure of the native β-barrel. In addition, CSPs support the presence of two features in the denatured state: a major hydrophobic cluster involving residues from various parts of the sequence and a native-like interaction similar to one identified in previous studies as forming early in folding (Budyak et al., Structure 21, 476 [2013]). Taken together, our data support a model in which transient structures involving nonlocal interactions prime early folding interactions in CRABP1, determine its barrel topology, and may protect this predominantly β-sheet protein against aggregation.

Overview of intranasally delivered peptides: key considerations for pharmaceutical development

INTRODUCTION

Intranasal (IN) delivery for peptides provides unique advantages compared to other invasive systemic delivery routes. However, there still lacks a clear understanding on how to evaluate the potential of the peptides for nasal delivery and key considerations for the nasal formulation development.

AREAS COVERED

A retrospective analysis of intranasally delivered peptides was conducted. The goals of this undertaking were 1) to build a database of the key physicochemical and pharmacokinetic properties of peptides delivered by the nasal route, 2) to evaluate formulation attributes applied to IN peptide delivery systems, and 3) to provide key considerations for IN delivery of peptides.

EXPERT OPINION/COMMENTARY

Extensive data mining showed that peptides with molecular weights up to 6000 Da have been delivered intranasally. The high solubility of some peptides highlighted the possibility of delivering sufficient amounts of peptide in the limited volume available for nasal sprays. Permeation enhancers and mucoadhesives have shown promise in improving the IN bioavailability of peptides. Other formulation considerations, such as the type of formulation, pH, osmolality, as well as drug deposition, are reviewed herein. Based on this retrospective analysis, key considerations for nasal peptides formulations were proposed to guide drug discovery and development for IN delivery of peptides.

The lineage-specific, intrinsically disordered N-terminal extension of monothiol glutaredoxin 1 from trypanosomes contains a regulatory region

Glutaredoxins (Grx) are small proteins conserved throughout all the kingdoms of life that are engaged in a wide variety of biological processes and share a common thioredoxin-fold. Among them, class II Grx are redox-inactive proteins involved in iron-sulfur (FeS) metabolism. They contain a single thiol group in their active site and use low molecular mass thiols such as glutathione as ligand for binding FeS-clusters. In this study, we investigated molecular aspects of 1CGrx1 from the pathogenic parasite Trypanosoma brucei brucei, a mitochondrial class II Grx that fulfills an indispensable role in vivo. Mitochondrial 1CGrx1 from trypanosomes differs from orthologues in several features including the presence of a parasite-specific N-terminal extension (NTE) whose role has yet to be elucidated. Previously we have solved the structure of a truncated form of 1CGrx1 containing only the conserved glutaredoxin domain but lacking the NTE. Our aim here is to investigate the effect of the NTE on the conformation of the protein. We therefore solved the NMR structure of the full-length protein, which reveals subtle but significant differences with the structure of the NTE-less form. By means of different experimental approaches, the NTE proved to be intrinsically disordered and not involved in the non-redox dependent protein dimerization, as previously suggested. Interestingly, the portion comprising residues 65–76 of the NTE modulates the conformational dynamics of the glutathione-binding pocket, which may play a role in iron-sulfur cluster assembly and delivery. Furthermore, we disclosed that the class II-strictly conserved loop that precedes the active site is critical for stabilizing the protein structure. So far, this represents the first communication of a Grx containing an intrinsically disordered region that defines a new protein subgroup within class II Grx.

GRPY: An Accurate Bead Method for Calculation of Hydrodynamic Properties of Rigid Biomacromolecules

Two main problems that arise in the context of hydrodynamic bead modeling are an inaccurate treatment of bead overlaps and the necessity of using volume corrections when calculating intrinsic viscosity. We present a formalism based on the generalized Rotne-Prager-Yamakawa approximation that successfully addresses both of these issues. The generalized Rotne-Prager-Yamakawa method is shown to be highly effective for the calculation of transport properties of rigid biomolecules represented as assemblies of spherical beads of different sizes, both overlapping and nonoverlapping. We test the method on simple molecular shapes as well as real protein structures and compare its performance with other computational approaches.

A Method for Determining Structure Ensemble of Large Disordered Protein: Application to a Mechanosensing Protein

Structure characterization of intrinsically disordered proteins (IDPs) remains a key obstacle in understanding their functional mechanisms. Due to the highly dynamic feature of IDPs, structure ensembles instead of static unique structures are often derived from experimental data. Several state-of-the-art computational methods have been developed to select an optimal ensemble from a pregenerated structure pool, but they suffer from low efficiency for large IDPs. Here we present a matching pursuit genetic algorithm (MPGA) for structure ensemble determination, which takes advantages from both matching pursuit (MP) to reduce the search space and genetic algorithm (GA) to reduce the restriction on constraint types. The MPGA method is validated using a reference ensemble with predefined structures. In comparison with the conventional GA, MPGA takes much less computational time for large IDPs. The utility of the method is demonstrated by application to structure ensemble determination of a mechanosensing protein domain with 306 amino acids. The structure ensemble determined reveals that the N-terminal region 1-240 is more compact than the C-terminal region 240-306. The unique structural feature explains why only a small portion of YXXP tyrosine residues can be phosphorylated easily by kinases in the absence of extension force and why the phosphorylation is force-dependent.

A Luminal Loop of Wilson Disease Protein Binds Copper and Is Required for Protein Activity

The copper-transporting ATPase ATP7B is essential for loading of copper ions to copper-dependent enzymes in the secretory pathway; its inactivation results in Wilson disease. In contrast to copper-ion uptake by the cytoplasmic domains, ATP7B-mediated copper-ion release in the Golgi has not been explored yet. We demonstrate here that a luminal loop in ATP7B, rich in histidine/methionine residues, binds reduced copper (Cu(I)) ions, and identified copper-binding residues play an essential role in ATP7B-mediated metal ion release. NMR experiments on short-peptide models demonstrate that three methionine and two histidine residues are specifically involved in Cu(I) ion binding; with these residues replaced by alanines, no Cu(I) ion interaction is detected. Although more than one Cu(I) ion can interact with the wild-type peptide, removing either all histidine or all methionine residues reduces the stoichiometry to one Cu(I) ion binding per peptide. Using a yeast complementation assay, we show that for efficient copper transport by full-length ATP7B, the complete set of histidine and methionine residues in the lumen loop are required. The replacement of histidine/methionine residues by alanines does not perturb overall ATP7B structure, as the localization of ATP7B variants in yeast cells matches that of the wild-type protein. Thus, in similarity to ATP7A, ATP7B also appears to have a luminal “exit” copper ion site.

Conformational response to charge clustering in synthetic intrinsically disordered proteins

BACKGROUND

Recent theoretical and computational studies have shown that the charge content and, most importantly, the linear distribution of opposite charges are major determinants of conformational properties of intrinsically disordered proteins (IDPs). Charge segregation in a sequence can be measured through κ, which represents a normalized measure of charge asymmetry. A strong inverse correlation between κ and radius of gyration has been previously demonstrated for two independent sets of permutated IDP sequences.

METHODS

We used two well-characterized IDPs, namely measles virus N_TAIL and Hendra virus PNT4, sharing a very similar fraction of charged residues and net charge per residue, but differing in proline (Pro) content. For each protein, we have rationally designed a low- and a high-κ variant endowed with the highest and the lowest κ values compatible with their natural amino acid composition. Then, the conformational properties of wild-type and κ-variants have been assessed by biochemical and biophysical techniques.

RESULTS

We confirmed a direct correlation between κ and protein compaction. The analysis of our original data along with those available from the literature suggests that Pro content may affects the responsiveness to charge clustering.

CONCLUSIONS

Charge clustering promotes IDP compaction, but the extent of its effects depends on the sequence context. Proline residues seem to play a role contrasting compaction.

GENERAL SIGNIFICANCE

These results contribute to the identification of sequence determinants of IDP conformational properties. They may also serve as an asset for rational design of non-natural IDPs with tunable degree of compactness.

Structural effects of methylglyoxal glycation, a study on the model protein MNEI

The reaction of free amino groups in proteins with reactive carbonyl species, known as glycation, leads to the formation of mixtures of products, collectively referred to as advanced glycation endproducts (AGEs). These compounds have been implicated in several important diseases, but their role in pathogenesis and clinical symptoms' development is still debated. Particularly, AGEs are often associated to the formation of amyloid deposits in conformational diseases, such as Alzheimer's and Parkinson's disease, and it has been suggested that they might influence the mechanisms and kinetics of protein aggregation. We here present the characterization of the products of glycation of the model protein MNEI with methylglyoxal and their effect on the protein structure. We demonstrate that, despite being an uncontrolled process, glycation occurs only at specific residues of the protein. Moreover, while not affecting the protein fold, it alters its shape and hydrodynamic properties and increases its tendency to fibrillar aggregation. Our study opens the way to in deep structural investigations to shed light on the complex link between protein post-translational modifications, structure, and stability.

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA

We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.

The Cyanobacterial Ribosomal-Associated Protein LrtA from Synechocystis sp. PCC 6803 Is an Oligomeric Protein in Solution with Chameleonic Sequence Properties

The LrtA protein of Synechocystis sp. PCC 6803 intervenes in cyanobacterial post-stress survival and in stabilizing 70S ribosomal particles. It belongs to the hibernating promoting factor (HPF) family of proteins, involved in protein synthesis. In this work, we studied the conformational preferences and stability of isolated LrtA in solution. At physiological conditions, as shown by hydrodynamic techniques, LrtA was involved in a self-association equilibrium. As indicated by Nuclear Magnetic Resonance (NMR), circular dichroism (CD) and fluorescence, the protein acquired a folded, native-like conformation between pH 6.0 and 9.0. However, that conformation was not very stable, as suggested by thermal and chemical denaturations followed by CD and fluorescence. Theoretical studies of its highly-charged sequence suggest that LrtA had a Janus sequence, with a context-dependent fold. Our modelling and molecular dynamics (MD) simulations indicate that the protein adopted the same fold observed in other members of the HPF family (β-α-β-β-β-α) at its N-terminal region (residues 1–100), whereas the C terminus (residues 100–197) appeared disordered and collapsed, supporting the overall percentage of overall secondary structure obtained by CD deconvolution. Then, LrtA has a chameleonic sequence and it is the first member of the HPF family involved in a self-association equilibrium, when isolated in solution.

19F NMR as a Tool for Monitoring Individual Differentially Labeled Proteins in Complex Mixtures

The ability to monitor the behavior of individual proteins in complex mixtures has many potential uses, ranging from analysis of protein interactions in highly concentrated solutions, modeling biological fluids or the intracellular environment, to optimizing biopharmaceutical co-formulations. Differential labeling NMR approaches, which traditionally use ¹⁵N or ¹³C isotope incorporation during recombinant expression, are not always practical in cases when endogenous proteins are obtained from an organism, or where the expression system does not allow for efficient labeling, especially for larger proteins. This study proposes differential labeling of proteins by covalent attachment of ¹⁹F groups with distinct chemical shifts, giving each protein a unique spectral signature which can be monitored by ¹⁹F NMR without signal overlap, even in complex mixtures, and without any interfering signals from the buffer or other unlabeled components. Parameters, such as signal intensities, translational diffusion coefficients, and transverse relaxation rates, which report on the behavior of individual proteins in the mixture, can be recorded even for proteins as large as antibodies at a wide range of concentrations.

Rotational Diffusion Depends on Box Size in Molecular Dynamics Simulations

We show that the rotational dynamics of proteins and nucleic acids determined from molecular dynamics simulations under periodic boundary conditions suffer from significant finite-size effects. We remove the box-size dependence of the rotational diffusion coefficients by adding a hydrodynamic correction k_B T/6 ηV with k_B Boltzmann's constant, T the absolute temperature, η the solvent shear viscosity, and V the box volume. We show that this correction accounts for the finite-size dependence of the rotational diffusion coefficients of horse-heart myoglobin and a B-DNA dodecamer in aqueous solution. The resulting hydrodynamic radii are in excellent agreement with experiment.

Enhancing clot properties through fibrin-specific self-cross-linked PEG side-chain microgels

Excessive bleeding and resulting complications are a major cause of death in both trauma and surgical settings. Recently, there have been a number of investigations into the design of synthetic hemostatic agents with platelet-mimicking activity to effectively treat patients suffering from severe hemorrhage. We developed platelet-like particles from microgels composed of polymers carrying polyethylene glycol (PEG) side-chains and fibrin-targeting single domain variable fragment antibodies (PEG-PLPs). Comparable to natural platelets, PEG-PLPs were found to enhance the fibrin network formation in vitro through strong adhesion to the emerging fibrin clot and physical, non-covalent cross-linking of nascent fibrin fibers. Furthermore, the mechanical reinforcement of the fibrin mesh through the incorporation of particles into the network leads to a ∼three-fold decrease of the overall clot permeability as compared to control clots. However, transport of biomolecules through the fibrin clots, such as peptides and larger proteins is not hindered by the presence of PEG-PLPs and the altered microstructure. Compared to control clots with an elastic modulus of 460+/-260 Pa, PEG-PLP-reinforced fibrin clots exhibit higher degrees of stiffness as demonstrated by the significantly increased average Younǵs modulus of 1770 +/±720 Pa, as measured by AFM force spectroscopy. Furthermore, in vitro degradation studies with plasmin demonstrate that fibrin clots formed in presence of PEG-PLPs withstand hydrolysis for 24 h, indicating enhanced stabilization against exogenous fibrinolysis. The entire set of data suggests that the designed platelet-like particles have high potential for use as hemostatic agents in emergency medicine and surgical settings.

Collapse Transitions of Proteins and the Interplay Among Backbone, Sidechain, and Solvent Interactions

Proteins can collapse into compact globules or form expanded, solvent-accessible, coil-like conformations. Additionally, they can fold into well-defined three-dimensional structures or remain partially or entirely disordered. Recent discoveries have shown that the tendency for proteins to collapse or remain expanded is not intrinsically coupled to their ability to fold. These observations suggest that proteins do not have to form compact globules in aqueous solutions. They can be intrinsically disordered, collapsed, or expanded, and even form well-folded, elongated structures. This ability to decouple collapse from folding is determined by the sequence details of proteins. In this review, we highlight insights gleaned from studies over the past decade. Using a polymer physics framework, we explain how the interplay among sidechains, backbone units, and solvent determines the driving forces for collapsed versus expanded states in aqueous solvents.

Order, Disorder, and Temperature-Driven Compaction in a Designed Elastin Protein

Artificial minielastin constructs have been designed that replicate the structure and function of natural elastins in a simpler context, allowing the NMR observation of structure and dynamics of elastin-like proteins with complete residue-specific resolution. We find that the alanine-rich cross-linking domains of elastin have a partially helical structure, but only when capped by proline-rich hydrophobic domains. We also find that the hydrophobic domains, composed of prominent 6-residue repeats VPGVGG and APGVGV found in natural elastins, appear random coil by both NMR chemical shift analysis and circular dichroism. However, these elastin hydrophobic domains exhibit structural bias for a dynamically disordered conformation that is neither helical nor β sheet with a degree of nonrandom structural bias which is dependent on residue type and position in the sequence. Another nonrandom-coil aspect of hydrophobic domain structure lies in the fact that, in contrast to other intrinsically disordered proteins, these hydrophobic domains retain a relatively condensed conformation whether attached to cross-linking domains or not. Importantly, these domains and the proteins containing them constrict with increasing temperature by up to 30% in volume without becoming more ordered. This property is often observed in nonbiological polymers and suggests that temperature-driven constriction is a new type of protein structural change that is linked to elastin's biological functions of coacervation-driven assembly and elastic recoil.

Improving the Interpretation of Small Molecule Diffusion Coefficients

Diffusion-ordered NMR spectroscopy (DOSY) is increasingly widely used for the analysis of mixtures by NMR spectroscopy, dispersing the signals of different species according to their diffusion coefficients. DOSY is used primarily to distinguish between the signals of different species, with the interpretation of the diffusion coefficients observed usually being purely qualitative, for example to deduce whether one species is bigger or smaller than another. In principle, the actual values of diffusion coefficient obtained carry important information about the sizes of different species and on interactions between species, but the relationship between diffusion coefficient and molecular mass is in general a very complex one. Here a recently proposed analytical relationship between diffusion coefficient and molecular mass for the restricted case of small organic molecules is tested against a wide range of data from the scientific literature and generalized to cover a range of solvents and temperatures.

HullRad: Fast Calculations of Folded and Disordered Protein and Nucleic Acid Hydrodynamic Properties

Hydrodynamic properties are useful parameters for estimating the size and shape of proteins and nucleic acids in solution. The calculation of such properties from structural models informs on the solution properties of these molecules and complements corresponding structural studies. Here we report, to our knowledge, a new method to accurately predict the hydrodynamic properties of molecular structures. This method uses a convex hull model to estimate the hydrodynamic volume of the molecule and is orders of magnitude faster than common methods. It works well for both folded proteins and ensembles of conformationally heterogeneous proteins and for nucleic acids. Because of its simplicity and speed, the method should be useful for the modification of computer-generated, intrinsically disordered protein ensembles and ensembles of flexible, but folded, molecules in which rapid calculation of experimental parameters is needed. The convex hull method is implemented in a Python script called HullRad. The use of the method is facilitated by a web server and the code is freely available for batch applications.

Unusual Self-Assembly of the Recombinant Chlamydia trachomatis Major Outer Membrane Protein-Based Fusion Antigen CTH522 Into Protein Nanoparticles

Sexually transmitted Chlamydia trachomatis infects more than 100 million people annually, and untreated chlamydia infections can cause severe complications. Therefore, there is an urgent need for a chlamydia vaccine. The Ctrachomatis major outer membrane protein (MOMP) is highly immunogenic but is a challenging vaccine candidate by being an integral membrane protein, and the immunogenicity depends on a correctly folded structure. We investigated the biophysical properties of the recombinant MOMP-based fusion antigen CTH522, which is tested in early human clinical trials. It consists of a truncated and cysteine-free version of MOMP fused to 4 variable domains from serovars D-G. In the native state, CTH522 did not exist as a monomer but showed an unusual self-assembly into nanoparticles with a negative zeta potential. In contrast to the β-barrel structure of MOMP, native CTH522 contained no well-defined secondary structural elements, and no thermal transitions were measurable. Chemical unfolding resulted in monomers that upon removal of the denaturant self-assembled into higher order structures, comparable to the structure of the native protein. The conformation of CTH522 in nanoparticles is thus not entirely random and contains structural elements stabilized via denaturant-disruptable hydrophobic interactions. In conclusion, CTH522 has an unusual quaternary structure of supramolecular self-assemblies.

On-chip measurements of protein unfolding from direct observations of micron-scale diffusion

The unfolding process of BSA in solution as a function of pH was studied by microfluidic diffusional sizing device.

Investigations of protein folding, unfolding and stability are critical for the understanding of the molecular basis of biological structure and function. We describe here a microfluidic approach to probe the unfolding of unlabelled protein molecules in microliter volumes. We achieve this objective through the use of a microfluidic platform, which allows the changes in molecular diffusivity upon folding and unfolding to be detected directly. We illustrate this approach by monitoring the unfolding of bovine serum albumin in solution as a function of pH. These results show the viability of probing protein stability on chip in small volumes.

Ion pair and solvation dynamics of [Bmim][BF4 ] + water system

In this work, 1-butyl-3-methylimidazolium tetrafluoroborate/water mixtures were analysed over the whole water composition (x_w ) in order to study the rotational and translational behaviour of the ions. We employed a multinuclear NMR approach to determine anion/cation/water diffusion coefficients and longitudinal relaxation rates at different water content. In neat ionic liquids (IL), the cation diffuses faster than the anion, and at low x_w , anions and cations share almost the same diffusion coefficient, but above a critical water concentration, the anion begins to diffuse faster than the cation. We identified this composition as approximately 10% x_w where the ions share the same diffusion coefficient. We found that the water at this composition seems to have a much more dramatic effect in the rotational diffusion of the anion that decreases substantially and approaches that of the anion in the diluted IL. Translational and rotational dynamics of the ions suggest that water is first incorporated in pockets in the nanostructure of the IL allowing the ions to maintain most of the cation/anion interactions present in neat IL but already disrupting some anion/cation interactions due to preferential interaction with the anion. HOESY and NOESY data show that water displays contacts both with the cation and the anion in a positive NOE regime in contrary to the negative regime found for the cation/anion and cation/cation cross-relaxation. This is in accordance with the high relative diffusion coefficient of water and suggests that water molecules can exchange between preferential location sites that allow water to maintain contacts both with the anion and cation.

Concerted millisecond timescale dynamics in the intrinsically disordered carboxyl terminus of γ-tubulin induced by mutation of a conserved tyrosine residue

Tubulins are an ancient family of eukaryotic proteins characterized by an amino-terminal globular domain and disordered carboxyl terminus. These carboxyl termini play important roles in modulating the behavior of microtubules in living cells. However, the atomic-level basis of their function is not well understood. These regions contain multiple acidic residues and their overall charges are modulated in vivo by post-translational modifications, for example, phosphorylation. In this study, we describe an application of NMR and computer Monte Carlo simulations to investigate how the modification of local charge alters the conformational sampling of the γ-tubulin carboxyl terminus. We compared the dynamics of two 39-residue polypeptides corresponding to the carboxyl-terminus of yeast γ-tubulin. One polypeptide comprised the wild-type amino acid sequence while the second contained a Y > D mutation at Y11 in the polypeptide (Y445 in the full protein). This mutation introduces additional negative charge at a site that is phosphorylated in vivo and produces a phenotype with perturbed microtubule function. NMR relaxation measurements show that the Y11D mutation produces dramatic changes in the millisecond-timescale motions of the entire polypeptide. This observation is supported by Monte Carlo simulations that-similar to NMR-predict the WT γ-CT is largely unstructured and that the substitution of Tyr 11 with Asp causes the sampling of extended conformations that are unique to the Y11D polypeptide.

Mechanistic View of hnRNPA2 Low-Complexity Domain Structure, Interactions, and Phase Separation Altered by Mutation and Arginine Methylation

hnRNPA2, a component of RNA-processing membraneless organelles, forms inclusions when mutated in a syndrome characterized by the degeneration of neurons (bearing features of amyotrophic lateral sclerosis [ALS] and frontotemporal dementia), muscle, and bone. Here we provide a unified structural view of hnRNPA2 self-assembly, aggregation, and interaction and the distinct effects of small chemical changes-disease mutations and arginine methylation-on these assemblies. The hnRNPA2 low-complexity (LC) domain is compact and intrinsically disordered as a monomer, retaining predominant disorder in a liquid-liquid phase-separated form. Disease mutations D290V and P298L induce aggregation by enhancing and extending, respectively, the aggregation-prone region. Co-aggregating in disease inclusions, hnRNPA2 LC directly interacts with and induces phase separation of TDP-43. Conversely, arginine methylation reduces hnRNPA2 phase separation, disrupting arginine-mediated contacts. These results highlight the mechanistic role of specific LC domain interactions and modifications conserved across many hnRNP family members but altered by aggregation-causing pathological mutations.

Photoreduction of ferricytochrome c in the presence of potassium ferrocyanide

Ferricytochrome c has been previously shown to photoreduce in the presence of ferrocyanide anions, but the process has been poorly understood.

The histidine phosphocarrier protein, HPr, binds to the highly thermostable regulator of sigma D protein, Rsd, and its isolated helical fragments

The phosphotransferase system (PTS) controls the preferential use of sugars in bacteria and it is also involved in other processes, such as chemotaxis. It is formed by a protein cascade in which the first two proteins are general (namely, EI and HPr) and the others are sugar-specific permeases. The Rsd protein binds specifically to the RNA polymerase (RNAP) σ⁷⁰ factor. We first characterized the conformational stability of Escherichia coli Rsd. And second, we delineated the binding regions of Streptomyces coelicolor, HPr^sc, and E. coli Rsd, by using fragments derived from each protein. To that end, we used several biophysical probes, namely, fluorescence, CD, NMR, ITC and BLI. Rsd had a free energy of unfolding of 15 kcal mol^-1 at 25 °C, and a thermal denaturation midpoint of 103 °C at pH 6.5. The affinity between Rsd and HPr^sc was 2 μM. Interestingly enough, the isolated helical-peptides, comprising the third (RsdH3) and fourth (RsdH4) Rsd helices, also interacted with HPr^sc in a specific manner, and with affinities similar to that of the whole Rsd. Moreover, the isolated peptide of HPr^sc, HPr^9-30, comprising the active site, His15, also was bound to intact Rsd with similar affinity. Therefore, binding between Rsd and HPr^sc was modulated by the two helices H3 and H4 of Rsd, and the regions around the active site of HPr^sc. This implies that specific fragments of Rsd and HPr^sc can be used to interfere with other protein-protein interactions (PPIs) of each other protein.

The energy cost of polypeptide knot formation and its folding consequences

Knots are natural topologies of chains. Yet, little is known about spontaneous knot formation in a polypeptide chain—an event that can potentially impair its folding—and about the effect of a knot on the stability and folding kinetics of a protein. Here we used optical tweezers to show that the free energy cost to form a trefoil knot in the denatured state of a polypeptide chain of 120 residues is 5.8 ± 1 kcal mol⁻¹. Monte Carlo dynamics of random chains predict this value, indicating that the free energy cost of knot formation is of entropic origin. This cost is predicted to remain above 3 kcal mol⁻¹ for denatured proteins as large as 900 residues. Therefore, we conclude that naturally knotted proteins cannot attain their knot randomly in the unfolded state but must pay the cost of knotting through contacts along their folding landscape.

The effect of knots on protein stability and folding kinetics is not well understood. Here the authors combine optical tweezer experiments and calculations to experimentally determine the energy cost for knot formation, which indicates that knotted proteins evolved specific folding pathways because knot formation in unfolded chains is unfavorable.

Structural disorder and induced folding within two cereal, ABA stress and ripening (ASR) proteins

Abscisic acid (ABA), stress and ripening (ASR) proteins are plant-specific proteins involved in plant response to multiple abiotic stresses. We previously isolated the ASR genes and cDNAs from durum wheat (TtASR1) and barley (HvASR1). Here, we show that HvASR1 and TtASR1 are consistently predicted to be disordered and further confirm this experimentally. Addition of glycerol, which mimics dehydration, triggers a gain of structure in both proteins. Limited proteolysis showed that they are highly sensitive to protease degradation. Addition of 2,2,2-trifluoroethanol (TFE) however, results in a decreased susceptibility to proteolysis that is paralleled by a gain of structure. Mass spectrometry analyses (MS) led to the identification of a protein fragment resistant to proteolysis. Addition of zinc also induces a gain of structure and Hydrogen/Deuterium eXchange-Mass Spectrometry (HDX-MS) allowed identification of the region involved in the disorder-to-order transition. This study is the first reported experimental characterization of HvASR1 and TtASR1 proteins, and paves the way for future studies aimed at unveiling the functional impact of the structural transitions that these proteins undergo in the presence of zinc and at achieving atomic-resolution conformational ensemble description of these two plant intrinsically disordered proteins (IDPs).

The dynamics of the monomeric restriction endonuclease BcnI during its interaction with DNA

Endonucleases that generate DNA double strand breaks often employ two independent subunits such that the active site from each subunit cuts either DNA strand. Restriction enzyme BcnI is a remarkable exception. It binds to the 5΄-CC/SGG-3΄ (where S = C or G, ‘/’ designates the cleavage position) target as a monomer forming an asymmetric complex, where a single catalytic center approaches the scissile phosphodiester bond in one of DNA strands. Bulk kinetic measurements have previously shown that the same BcnI molecule cuts both DNA strands at the target site without dissociation from the DNA. Here, we analyse the BcnI DNA binding and target recognition steps at the single molecule level. We find, using FRET, that BcnI adopts either ‘open’ or ‘closed’ conformation in solution. Next, we directly demonstrate that BcnI slides over long distances on DNA using 1D diffusion and show that sliding is accompanied by occasional jumping events, where the enzyme leaves the DNA and rebinds immediately at a distant site. Furthermore, we quantify the dynamics of the BcnI interactions with cognate and non-cognate DNA, and determine the preferred binding orientation of BcnI to the target site. These results provide new insights into the intricate dynamics of BcnI–DNA interactions.

Structural and hydrodynamic properties of an intrinsically disordered region of a germ cell-specific protein on phase separation

Membrane encapsulation is frequently used by the cell to sequester biomolecules and compartmentalize their function. Cells also concentrate molecules into phase-separated protein or protein/nucleic acid "membraneless organelles" that regulate a host of biochemical processes. Here, we use solution NMR spectroscopy to study phase-separated droplets formed from the intrinsically disordered N-terminal 236 residues of the germ-granule protein Ddx4. We show that the protein within the concentrated phase of phase-separated Ddx4, [Formula: see text], diffuses as a particle of 600-nm hydrodynamic radius dissolved in water. However, NMR spectra reveal sharp resonances with chemical shifts showing [Formula: see text] to be intrinsically disordered. Spin relaxation measurements indicate that the backbone amides of [Formula: see text] have significant mobility, explaining why high-resolution spectra are observed, but motion is reduced compared with an equivalently concentrated nonphase-separating control. Observation of a network of interchain interactions, as established by NOE spectroscopy, shows the importance of Phe and Arg interactions in driving the phase separation of Ddx4, while the salt dependence of both low- and high-concentration regions of phase diagrams establishes an important role for electrostatic interactions. The diffusion of a series of small probes and the compact but disordered 4E binding protein 2 (4E-BP2) protein in [Formula: see text] are explained by an excluded volume effect, similar to that found for globular protein solvents. No changes in structural propensities of 4E-BP2 dissolved in [Formula: see text] are observed, while changes to DNA and RNA molecules have been reported, highlighting the diverse roles that proteinaceous solvents play in dictating the properties of dissolved solutes.

Arabidopsis calmodulin-like protein CML36 is a calcium (Ca2+) sensor that interacts with the plasma membrane Ca2+-ATPase isoform ACA8 and stimulates its activity

Calmodulin-like (CML) proteins are major EF-hand-containing, calcium (Ca2+)-binding proteins with crucial roles in plant development and in coordinating plant stress tolerance. Given their abundance in plants, the properties of Ca2+ sensors and identification of novel target proteins of CMLs deserve special attention. To this end, we recombinantly produced and biochemically characterized CML36 from Arabidopsis thaliana We analyzed Ca2+ and Mg2+ binding to the individual EF-hands, observed metal-induced conformational changes, and identified a physiologically relevant target. CML36 possesses two high-affinity Ca2+/Mg2+ mixed binding sites and two low-affinity Ca2+-specific sites. Binding of Ca2+ induced an increase in the α-helical content and a conformational change that lead to the exposure of hydrophobic regions responsible for target protein recognition. Cation binding, either Ca2+ or Mg2+, stabilized the secondary and tertiary structures of CML36, guiding a large structural transition from a molten globule apo-state to a compact holoconformation. Importantly, through in vitro binding and activity assays, we showed that CML36 interacts directly with the regulative N terminus of the Arabidopsis plasma membrane Ca2+-ATPase isoform 8 (ACA8) and that this interaction stimulates ACA8 activity. Gene expression analysis revealed that CML36 and ACA8 are co-expressed mainly in inflorescences. Collectively, our results support a role for CML36 as a Ca2+ sensor that binds to and modulates ACA8, uncovering a possible involvement of the CML protein family in the modulation of plant-autoinhibited Ca2+ pumps.

Functions of intrinsic disorder in transmembrane proteins

Intrinsic disorder is common in integral membrane proteins, particularly in the intracellular domains. Despite this observation, these domains are not always recognized as being disordered. In this review, we will discuss the biological functions of intrinsically disordered regions of membrane proteins, and address why the flexibility afforded by disorder is mechanistically important. Intrinsically disordered regions are present in many common classes of membrane proteins including ion channels and transporters; G-protein coupled receptors (GPCRs), receptor tyrosine kinases and cytokine receptors. The functions of the disordered regions are many and varied. We will discuss selected examples including: (1) Organization of receptors, kinases, phosphatases and second messenger sources into signaling complexes. (2) Modulation of the membrane-embedded domain function by ball-and-chain like mechanisms. (3) Trafficking of membrane proteins. (4) Transient membrane associations. (5) Post-translational modifications most notably phosphorylation and (6) disorder-linked isoform dependent function. We finish the review by discussing the future challenges facing the membrane protein community regarding protein disorder.

An Efficient Method for Estimating the Hydrodynamic Radius of Disordered Protein Conformations

Intrinsically disordered proteins play important roles throughout biology, yet our understanding of the relationship between their sequences, structural properties, and functions remains incomplete. The dynamic nature of these proteins, however, makes them difficult to characterize structurally. Many disordered proteins can attain both compact and expanded conformations, and the level of expansion may be regulated and important for function. Experimentally, the level of compaction and shape is often determined either by small-angle x-ray scattering experiments or pulsed-field-gradient NMR diffusion measurements, which provide ensemble-averaged estimates of the radius of gyration and hydrodynamic radius, respectively. Often, these experiments are interpreted using molecular simulations or are used to validate them. We here provide, to our knowledge, a new and efficient method to calculate the hydrodynamic radius of a disordered protein chain from a model of its structural ensemble. In particular, starting from basic concepts in polymer physics, we derive a relationship between the radius of gyration of a structure and its hydrodynamic ratio, which in turn can be used, for example, to compare a simulated ensemble of conformations to NMR diffusion measurements. The relationship may also be valuable when using NMR diffusion measurements to restrain molecular simulations.

Detection of long-range electrostatic interactions between charged molecules by means of fluorescence correlation spectroscopy

In the present paper, an experimental feasibility study on the detection of long-range intermolecular interactions through three-dimensional molecular diffusion in solution is performed. This follows recent theoretical and numerical analyses reporting that long-range electrodynamic forces between biomolecules could be identified through deviations from Brownian diffusion. The suggested experimental technique was fluorescence correlation spectroscopy (FCS). By considering two oppositely charged molecular species in aqueous solution, namely, lysozymes and fluorescent dye molecules (Alexa488), the diffusion coefficient of the dyes has been measured for different values of the concentration of lysozyme, that is, for different average distances between the oppositely charged molecules. For our model, long-range interactions are of electrostatic origin, suggesting that their action radius can be varied by changing the ionic strength of the solution. The experimental outcomes clearly prove the detectability of long-range intermolecular interactions by means of the FCS technique. Molecular dynamics simulations provide a clear and unambiguous interpretation of the experimental results.

Rapid-Acting and Human Insulins: Hexamer Dissociation Kinetics upon Dilution of the Pharmaceutical Formulation

Purpose

Comparison of the dissociation kinetics of rapid-acting insulins lispro, aspart, glulisine and human insulin under physiologically relevant conditions.

Methods

Dissociation kinetics after dilution were monitored directly in terms of the average molecular mass using combined static and dynamic light scattering. Changes in tertiary structure were detected by near-UV circular dichroism.

Results

Glulisine forms compact hexamers in formulation even in the absence of Zn²⁺. Upon severe dilution, these rapidly dissociate into monomers in less than 10 s. In contrast, in formulations of lispro and aspart, the presence of Zn²⁺ and phenolic compounds is essential for formation of compact R6 hexamers. These slowly dissociate in times ranging from seconds to one hour depending on the concentration of phenolic additives. The disadvantage of the long dissociation times of lispro and aspart can be diminished by a rapid depletion of the concentration of phenolic additives independent of the insulin dilution. This is especially important in conditions similar to those after subcutaneous injection, where only minor dilution of the insulins occurs.

Conclusion

Knowledge of the diverging dissociation mechanisms of lispro and aspart compared to glulisine will be helpful for optimizing formulation conditions of rapid-acting insulins.

Electronic supplementary material

The online version of this article (doi:10.1007/s11095-017-2233-0) contains supplementary material, which is available to authorized users.

Comprehensive structural analysis of designed incomplete polypeptide chains of the replicase nonstructural protein 1 from the severe acute respiratory syndrome coronavirus

The cotranslational folding is recognized as a very cooperative process that occurs after the nearly completion of the polypeptide sequence of a domain. Here we investigated the challenges faced by polypeptide segments of a non-vectorial β-barrel fold. Besides the biological interest behind the SARS coronavirus non-structural protein 1 (nsp1, 117 amino acids), this study model has two structural features that motivated its use in this work: 1- its recombinant production is dependent on the temperature, with greater solubility when expressed at low temperatures. This is an indication of the cotranslational guidance to the native protein conformation. 2- Conversely, nsp1 has a six-stranded, mixed parallel/antiparallel β-barrel with intricate long-range interactions, indicating it will need the full-length protein to fold properly. We used non-denaturing purification conditions that allowed the characterization of polypeptide chains of different lengths, mimicking the landscape of the cotranslational fold of a β-barrel, and avoiding the major technical hindrances of working with the nascent polypeptide bound to the ribosome. Our results showed partially folded states formed as soon as the amino acids of the second β-strand were present (55 amino acids). These partially folded states are different based on the length of polypeptide chain. The native α-helix (amino acids 24-37) was identified as a transient structure (~20-30% propensity). We identified the presence of regular secondary structure after the fourth native β-strand is present (89 amino acids), in parallel to the collapse to a non-native 3D structure. Interestingly the polypeptide sequences of the native strands β2, β3 and β4 have characteristics of α-helices. Our comprehensive analyses support the idea that incomplete polypeptide chains, such as the ones of nascent proteins much earlier than the end of the translation, adopt an abundance of specific transient folds, instead of disordered conformations.

Challenges in Predicting Protein-Protein Interactions from Measurements of Molecular Diffusivity

Dynamic light scattering can be used to measure the diffusivity of a protein within a formulation. The dependence of molecular diffusivity on protein concentration (traditionally expressed in terms of the interaction parameter k_D) is often used to infer whether protein-protein interactions are repulsive or attractive, resulting in solutions that are colloidally stable or unstable, respectively. However, a number of factors unrelated to intermolecular forces can also impact protein diffusion, complicating this interpretation. Here, we investigate the influence of multicomponent diffusion in a ternary protein-salt-water system on protein diffusion and k_D in the context of Nernst-Planck theory. This analysis demonstrates that large changes in protein diffusivity with protein concentration can result even for hard-sphere systems in the absence of protein-protein interactions. In addition, we show that dynamic light scattering measurements of diffusivity made at low ionic strength cannot be reliably used to detect protein conformational changes. We recommend comparing experimentally determined k_D values to theoretically predicted excluded-volume contributions, which will allow a more accurate assessment of protein-protein interactions.

The PipX Protein, When Not Bound to Its Targets, Has Its Signaling C-Terminal Helix in a Flexed Conformation

PipX, an 89-residue protein, acts as a coactivator of the global nitrogen regulator NtcA in cyanobacteria. NtcA-PipX interactions are regulated by 2-oxoglutarate (2-OG), an inverse indicator of the ammonia abundance, and by P_II, a protein that binds to PipX at low 2-OG concentrations. The structure of PipX, when bound to NtcA or P_II, consists of an N-terminal, five-stranded β-sheet (conforming a Tudor-like domain), and two long α-helices. These helices adopt either a flexed conformation, where they are in close contact and in an antiparallel mutual orientation, also packing against the β-sheet, or an open conformation (observed only in the P_II-PipX complex) where the last α-helix moves apart from the rest of the protein. The aim of this work was to study the structure and dynamics of isolated PipX in solution by NMR. The backbone chemical shifts, the hydrogen-exchange, and the NOE patterns indicated that the isolated, monomeric PipX structure was formed by an N-terminal five-stranded β-sheet and two C-terminal α-helices. Furthermore, the observed NOEs between the two helices, and of α-helix2 with β-strand2 suggested that PipX adopted a flexed conformation. The β-strands 1 and 5 were highly flexible, as shown by the lack of interstrand backbone-backbone NOEs; in addition, the ¹⁵N-dynamics indicated that the C terminus of β-strand4 and the following β-turn (Phe42-Thr47), and the C-cap of α-helix1 (Arg70-Asn71) were particularly mobile. These two regions could act as hinges, allowing PipX to interact with its partners, including PlmA in the newly recognized P_II-PipX-PlmA ternary complex.