Structural analysis of an equilibrium folding intermediate in the apoflavodoxin native ensemble by small-angle X-ray scattering

Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Despite advances in artificial intelligence methods, protein folding remains in many ways an enigma to be solved. Accurate computation of protein folding energetics could help drive fields such as protein and drug design and genetic interpretation. However, the challenge of calculating the state functions governing protein folding from first-principles remains unaddressed. We present here a simple approach that allows us to accurately calculate the energetics of protein folding. It is based on computing the energy of the folded and unfolded states at different temperatures using molecular dynamics simulations. From this, two essential quantities (ΔH and ΔCp) are obtained and used to calculate the conformational stability of the protein (ΔG). With this approach, we have successfully calculated the energetics of two- and three-state proteins, representatives of the major structural classes, as well as small stability differences (ΔΔG) due to changes in solution conditions or variations in an amino acid residue.

EFAMIX, a tool to decompose inline chromatography SAXS data from partially overlapping components

Small‐angle X‐ray scattering (SAXS) is an established technique for structural analysis of biological macromolecules in solution. During the last decade, inline chromatography setups coupling SAXS with size exclusion (SEC‐SAXS) or ion exchange (IEC‐SAXS) have become popular in the community. These setups allow one to separate individual components in the sample and to record SAXS data from isolated fractions, which is extremely important for subsequent data interpretation, analysis, and structural modeling. However, in case of partially overlapping elution peaks, inline chromatography SAXS may still yield scattering profiles from mixtures of components. The deconvolution of these scattering data into the individual fractions is nontrivial and potentially ambiguous. We describe a cross‐platform computer program, EFAMIX, for restoring the scattering and concentration profiles of the components based on the evolving factor analysis (EFA). The efficiency of the program is demonstrated in a number of simulated and experimental SEC‐SAXS data sets. Sensitivity and limitations of the method are explored, and its applicability to IEC‐SAXS data is discussed. EFAMIX requires minimal user intervention and is available to academic users through the program package ATSAS as from release 3.1.

Structure and thermodynamics of transient protein-protein complexes by chemometric decomposition of SAXS datasets

Transient biomolecular interactions play crucial roles in many cellular signaling and regulation processes. However, deciphering the structure of these assemblies is challenging owing to the difficulties in isolating complexes from the individual partners. The additive nature of small-angle X-ray scattering (SAXS) data allows for probing the species present in these mixtures, but decomposition into structural and thermodynamic information is difficult. We present a chemometric approach enabling the decomposition of titration SAXS data into species-specific information. Using extensive synthetic SAXS data, we demonstrate that robust decomposition can be achieved for titrations with a maximum fraction of complex of 0.5 that can be extended to 0.3 when two orthogonal titrations are simultaneously analyzed. The effect of the structural features, titration points, relative concentrations, and noise are thoroughly analyzed. The validation of the strategy with experimental data highlights the power of the approach to provide unique insights into this family of biomolecular assemblies.

REGALS: a general method to deconvolve X-ray scattering data from evolving mixtures

Mixtures of biological macromolecules are inherently difficult to study using structural methods, as increasing complexity presents new challenges for data analysis. Recently, there has been growing interest in studying evolving mixtures using small-angle X-ray scattering (SAXS) in conjunction with time-resolved, high-throughput or chromatography-coupled setups. Deconvolution and interpretation of the resulting datasets, however, are nontrivial when neither the scattering components nor the way in which they evolve are known a priori. To address this issue, the REGALS method (regularized alternating least squares) is introduced, which incorporates simple expectations about the data as prior knowledge, and utilizes parameterization and regularization to provide robust deconvolution solutions. The restraints used by REGALS are general properties such as smoothness of profiles and maximum dimensions of species, making it well suited for exploring datasets with unknown species. Here, REGALS is applied to the analysis of experimental data from four types of SAXS experiment: anion-exchange (AEX) coupled SAXS, ligand titration, time-resolved mixing and time-resolved temperature jump. Based on its performance with these challenging datasets, it is anticipated that REGALS will be a valuable addition to the SAXS analysis toolkit and enable new experiments. The software is implemented in both MATLAB and Python and is available freely as an open-source software package.

Structurally Tunable pH-responsive Phosphine Oxide Based Gels by Facile Synthesis Strategy

Design and synthesis of nanostructured responsive gels have attracted increasing attention, particularly in the biomedical domain. Polymer chain configurations and nanodomain sizes within the network can be used to steer their functions as drug carriers. Here, a catalyst-free facile one-step synthesis strategy is reported for the design of pH-responsive gels and controlled structures in nanoscale. Transparent and impurity free gels were directly synthesized from trivinylphosphine oxide (TVPO) and cyclic secondary diamine monomers via Michael addition polymerization under mild conditions. NMR analysis confirmed the consumption of all TVPO and the absence of side products, thereby eliminating post purification steps. The small-angle X-ray scattering (SAXS) elucidates the nanoscale structural features in gels, that is, it demonstrates the presence of collapsed nanodomains within gel networks and it was possible to tune the size of these domains by varying the amine monomers and the nature of the solvent. The fabricated gels demonstrate structure tunability via solvent-polymer interactions and pH specific drug release behavior. Three different anionic dyes (acid blue 80, acid blue 90, and fluorescein) of varying size and chemistry were incorporated into the hydrogel as model drugs and their release behavior was studied. Compared to acidic pH, a higher and faster release of acid blue 80 and fluorescein was observed at pH 10, possibly because of their increased solubility in alkaline pH. In addition, their release in phosphate buffered saline (PBS) and simulated body fluid (SBF) matrix was positively influenced by the ionic interaction with positively charged metal ions. In the case of hydrogel containing acid blue 90 a very low drug release (<1%) was observed, which is due to the reaction of its accessible free amino group with the vinyl groups of the TVPO. In vitro evaluation of the prepared hydrogel using human dermal fibroblasts indicates no cytotoxic effects, warranting further research for biomedical applications. Our strategy of such gel synthesis lays the basis for the design of other gel-based functional materials.

Direct examination of the relevance for folding, binding and electron transfer of a conserved protein folding intermediate

Near the minimum free energy basin of proteins where the native ensemble resides, partly unfolded conformations of slightly higher energy can be significantly populated under native conditions. It has been speculated that they play roles in molecular recognition and catalysis, but they might represent contemporary features of the evolutionary process without functional relevance. Obtaining conclusive evidence on these alternatives is difficult because it requires comparing the performance of a given protein when populating and when not populating one such intermediate, in otherwise identical conditions. Wild type apoflavodoxin populates under native conditions a partly unfolded conformation (10% of molecules) whose unstructured region includes the binding sites for the FMN cofactor and for redox partner proteins. We recently engineered a thermostable variant where the intermediate is no longer detectable. Using the wild type and variant, we assess the relevance of the intermediate comparing folding kinetics, cofactor binding kinetics, cofactor affinity, X-ray structure, intrinsic dynamics, redox potential of the apoflavodoxin-cofactor complex (Fld), its affinity for partner protein FNR, and electron transfer rate within the Fld/FNR physiological complex. Our data strongly suggest the intermediate state, conserved in long-chain apoflavodoxins, is not required for the correct assembly of flavodoxin nor does it contribute to shape its electron transfer properties. This analysis can be applied to evaluate other native basin intermediates.

Folding of proteins with a flavodoxin-like architecture

The flavodoxin-like fold is a protein architecture that can be traced back to the universal ancestor of the three kingdoms of life. Many proteins share this α-β parallel topology and hence it is highly relevant to illuminate how they fold. Here, we review experiments and simulations concerning the folding of flavodoxins and CheY-like proteins, which share the flavodoxin-like fold. These polypeptides tend to temporarily misfold during unassisted folding to their functionally active forms. This susceptibility to frustration is caused by the more rapid formation of an α-helix compared to a β-sheet, particularly when a parallel β-sheet is involved. As a result, flavodoxin-like proteins form intermediates that are off-pathway to native protein and several of these species are molten globules (MGs). Experiments suggest that the off-pathway species are of helical nature and that flavodoxin-like proteins have a nonconserved transition state that determines the rate of productive folding. Folding of flavodoxin from Azotobacter vinelandii has been investigated extensively, enabling a schematic construction of its folding energy landscape. It is the only flavodoxin-like protein of which cotranslational folding has been probed. New insights that emphasize differences between in vivo and in vitro folding energy landscapes are emerging: the ribosome modulates MG formation in nascent apoflavodoxin and forces this polypeptide toward the native state.

Streptococcus pneumoniae TIGR4 Flavodoxin: Structural and Biophysical Characterization of a Novel Drug Target

Streptococcus pneumoniae (Sp) strain TIGR4 is a virulent, encapsulated serotype that causes bacteremia, otitis media, meningitis and pneumonia. Increased bacterial resistance and limited efficacy of the available vaccine to some serotypes complicate the treatment of diseases associated to this microorganism. Flavodoxins are bacterial proteins involved in several important metabolic pathways. The Sp flavodoxin (Spfld) gene was recently reported to be essential for the establishment of meningitis in a rat model, which makes SpFld a potential drug target. To facilitate future pharmacological studies, we have cloned and expressed SpFld in E. coli and we have performed an extensive structural and biochemical characterization of both the apo form and its active complex with the FMN cofactor. SpFld is a short-chain flavodoxin containing 146 residues. Unlike the well-characterized long-chain apoflavodoxins, the Sp apoprotein displays a simple two-state thermal unfolding equilibrium and binds FMN with moderate affinity. The X-ray structures of the apo and holo forms of SpFld differ at the FMN binding site, where substantial rearrangement of residues at the 91-100 loop occurs to permit cofactor binding. This work will set up the basis for future studies aiming at discovering new potential drugs to treat S. pneumoniae diseases through the inhibition of SpFld.

The mechanism of water/ion exchange at a protein surface: a weakly bound chloride in Helicobacter pylori apoflavodoxin

Binding/unbinding of small ligands, such as ions, to/from proteins influences biochemical processes such as protein folding, enzyme catalysis or protein/ligand recognition. We have investigated the mechanism of chloride/water exchange at a protein surface (that of the apoflavodoxin from Helicobacter pylori) using classical all-atom molecular dynamics simulations. They reveal a variety of chloride exit routes and residence times; the latter is related to specific coordination modes of the anion. The role of solvent molecules in the mechanism of chloride unbinding has been studied in detail. We see no temporary increase in chloride coordination along the release process. Instead, the coordination of new water molecules takes place in most cases after the chloride/protein atom release event has begun. Moreover, the distribution function of water entrance events into the first chloride solvation shell peaks after chloride protein atom dissociation events. All these observations together seem to indicate that water molecules simply fill the vacancies left by the previously coordinating protein residues. We thus propose a step-by-step dissociation pathway in which protein/chloride interactions gradually break down before new water molecules progressively fill the vacant positions left by protein atoms. As observed for other systems, water molecules associated with bound chloride or with protein atoms have longer residence times than those bound to the free anion. The implications of the exchange mechanism proposed for the binding of the FMN (Flavin Mononucleotide) protein cofactor are discussed.

Rational stabilization of complex proteins: a divide and combine approach

Increasing the thermostability of proteins is often crucial for their successful use as analytic, synthetic or therapeutic tools. Most rational thermostabilization strategies were developed on small two-state proteins and, unsurprisingly, they tend to fail when applied to the much more abundant, larger, non-fully cooperative proteins. We show that the key to stabilize the latter is to know the regions of lower stability. To prove it, we have engineered apoflavodoxin, a non-fully cooperative protein on which previous thermostabilizing attempts had failed. We use a step-wise combination of structure-based, rationally-designed, stabilizing mutations confined to the less stable structural region, and obtain variants that, according to their van't Hoff to calorimetric enthalpy ratios, exhibit fully-cooperative thermal unfolding with a melting temperature of 75°C, 32 degrees above the lower melting temperature of the non-cooperative wild type protein. The ideas introduced here may also be useful for the thermostabilization of complex proteins through formulation or using specific stabilizing ligands (e.g. pharmacological chaperones).

Reverse micelles as a platform for dynamic nuclear polarization in solution NMR of proteins

Despite tremendous advances in recent years, solution NMR remains fundamentally restricted due to its inherent insensitivity. Dynamic nuclear polarization (DNP) potentially offers significant improvements in this respect. The basic DNP strategy is to irradiate the EPR transitions of a stable radical and transfer this nonequilibrium polarization to the hydrogen spins of water, which will in turn transfer polarization to the hydrogens of the macromolecule. Unfortunately, these EPR transitions lie in the microwave range of the electromagnetic spectrum where bulk water absorbs strongly, often resulting in catastrophic heating. Furthermore, the residence times of water on the surface of the protein in bulk solution are generally too short for efficient transfer of polarization. Here we take advantage of the properties of solutions of encapsulated proteins dissolved in low viscosity solvents to implement DNP in liquids. Such samples are largely transparent to the microwave frequencies required and thereby avoid significant heating. Nitroxide radicals are introduced into the reverse micelle system in three ways: attached to the protein, embedded in the reverse micelle shell, and free in the aqueous core. Significant enhancements of the water resonance ranging up to ∼-93 at 0.35 T were observed. We also find that the hydration properties of encapsulated proteins allow for efficient polarization transfer from water to the protein. These and other observations suggest that merging reverse micelle encapsulation technology with DNP offers a route to a significant increase in the sensitivity of solution NMR spectroscopy of proteins and other biomolecules.

Understanding the relevance of local conformational stability and dynamics to the aggregation propensity of an IgG1 and IgG2 monoclonal antibodies

Aggregation of monoclonal antibodies is often a multi-step process involving structural alterations in monomeric proteins and subsequent formation of soluble or insoluble oligomers. The role of local conformational stability and dynamics of native and/or partially altered structures in determining the aggregation propensity of monoclonal antibodies, however, is not well understood. Here, we investigate the role of conformational stability and dynamics of regions with distinct solvent exposure in determining the aggregation propensity of an IgG1 and IgG2 monoclonal antibody. The temperatures employed span the pre-unfolding range (10-40°C) and the onset temperatures (T onset ) for exposure of apolar residues (≈ 50°C), alterations in secondary structures (≈ 60°C) and initiation of visible aggregate formation (≈ 60°C). Solvent-exposed regions were found to precede solvent-shielded regions in an initiation of aggregation for both proteins. Such a process was observed upon alterations in overall tertiary structure while retaining the secondary structures in both the proteins. In addition, a greater dynamic nature of solvent-shielded regions in potential intermediates of IgG1 and the improved conformational stability increased its resistance to aggregation when compared to IgG2. These results suggest that local conformational stability and fluctuations of partially altered structures can influence the aggregation propensity of immunoglobulins.

Impact and progress in small and wide angle X-ray scattering (SAXS and WAXS)

The advances made in small and wide angle X-ray scattering over the past decades have had a large impact on structural biology. Many new insights into challenging biological probes including large and transient complexes, flexible macromolecules as well as other exciting objects of various sizes were gained with this low resolution technique. Here, we review the recent developments in the experimental setups and in software for data collection and analysis, specifically for hybrid approaches. These progresses have allowed scientists to address a number of intriguing questions which could not be answered with other structural methods alone.

Effect of surfactant hydrophobicity on the pathway for unfolding of ubiquitin

This paper describes the interaction between ubiquitin (UBI) and three sodium n-alkyl sulfates (SC(n)S) that have the same charge (Z = -1) but different hydrophobicity (n = 10, 12, or 14). Increasing the hydrophobicity of the n-alkyl sulfate resulted in (i) an increase in the number of distinct intermediates (that is, complexes of UBI and surfactant) that form along the pathway of unfolding, (ii) a decrease in the minimum concentrations of surfactant at which intermediates begin to form (i.e., a more negative ΔG(binding) of surfactant for UBI), and (iii) an increase in the number of surfactant molecules bound to UBI in each intermediate or complex. These results demonstrate that small changes in the hydrophobicity of a surfactant can significantly alter the binding interactions with a folded or unfolded cytosolic protein.

Protein dynamics governed by interfaces of high polarity and low packing density

The folding pathway, three-dimensional structure and intrinsic dynamics of proteins are governed by their amino acid sequences. Internal protein surfaces with physicochemical properties appropriate to modulate conformational fluctuations could play important roles in folding and dynamics. We show here that proteins contain buried interfaces of high polarity and low packing density, coined as LIPs: Light Interfaces of high Polarity, whose physicochemical properties make them unstable. The structures of well-characterized equilibrium and kinetic folding intermediates indicate that the LIPs of the corresponding native proteins fold late and are involved in local unfolding events. Importantly, LIPs can be identified using very fast and uncomplicated computational analysis of protein three-dimensional structures, which provides an easy way to delineate the protein segments involved in dynamics. Since LIPs can be retained while the sequences of the interacting segments diverge significantly, proteins could in principle evolve new functional features reusing pre-existing encoded dynamics. Large-scale identification of LIPS may contribute to understanding evolutionary constraints of proteins and the way protein intrinsic dynamics are encoded.

Defining the nature of thermal intermediate in 3 state folding proteins: apoflavodoxin, a study case

The early stages of the thermal unfolding of apoflavodoxin have been determined by using atomistic multi microsecond-scale molecular dynamics (MD) simulations complemented with a variety of experimental techniques. Results strongly suggest that the intermediate is reached very early in the thermal unfolding process and that it has the properties of an “activated” form of the native state, where thermal fluctuations in the loops break loop-loop contacts. The unrestrained loops gain then kinetic energy corrupting short secondary structure elements without corrupting the core of the protein. The MD-derived ensembles agree with experimental observables and draw a picture of the intermediate state inconsistent with a well-defined structure and characteristic of a typical partially disordered protein. Our results allow us to speculate that proteins with a well packed core connected by long loops might behave as partially disordered proteins under native conditions, or alternatively behave as three state folders. Small details in the sequence, easily tunable by evolution, can yield to one or the other type of proteins.

A simplistic view of protein structure tends to emphasize the opposition between the native state and the denatured ensemble of unfolded conformations. In addition to these extreme conformations, proteins subjected to a variety of perturbations often populate alternative partly unfolded conformations, some of which are close in energy to the native state and, accordingly, can be populated under native or quasi-native conditions. There is increasing evidence that these “perturbed” conformations participate in protein function or, in some cases, are related to the outcome of folding diseases. We have used the “state of the art” molecular dynamics combined with a variety of experimental techniques to characterize for the first time, to our knowledge, the thermal intermediate of a three-state folding protein (apoflavodoxin). Based on our results we have been able to suggest a general mechanism of thermal unfolding in complex proteins and to determine interesting links between thermal intermediates and partially unfolded proteins.

Structural analysis of intrinsically disordered proteins by small-angle X-ray scattering

Small-angle scattering of X-rays (SAXS) is an established method to study the overall structure and structural transitions of biological macromolecules in solution. For folded proteins, the technique provides three-dimensional low resolution structures ab initio or it can be used to drive rigid-body modeling. SAXS is also a powerful tool for the quantitative analysis of flexible systems, including intrinsically disordered proteins (IDPs), and is highly complementary to the high resolution methods of X-ray crystallography and NMR. Here we present the basic principles of SAXS and review the main approaches to the characterization of IDPs and flexible multidomain proteins using SAXS. Together with the standard approaches based on the analysis of overall parameters, a recently developed Ensemble Optimization Method (EOM) is now available. The latter method allows for the co-existence of multiple protein conformations in solution compatible with the scattering data. Analysis of the selected ensembles provides quantitative information about flexibility and also offers insights into structural features. Examples of the use of SAXS and combined approaches with NMR, X-ray crystallography, and computational methods to characterize completely or partially disordered proteins are presented.

Mechanism of FMN Binding to the Apoflavodoxin from Helicobacter pylori

Flavodoxins are bacterial electron transport proteins whose redox competence is due to the presence of a tightly but noncovalently bound FMN molecule. While the thermodynamics of the complex are understood, the mechanism of association between the apoflavodoxin and the redox cofactor is not so clear. We investigate here the mechanism of FMN binding to the apoflavodoxin from Helicobacter pylori, an essential protein that is being used as a target to develop antimicrobials. This flavodoxin is structurally peculiar as it lacks the typical bulky residue interacting with the FMN re face but bears instead a small alanine. FMN binding is biphasic, regardless of the presence of phosphate molecules in solution, while riboflavin binding takes place in a single step, the rate constant of which coincides with the fast phase of FMN binding. A mutational study at the isoalloxazine and phosphate subsites for FMN binding clearly indicates that FMN association is always limited by interaction with the isoalloxazine subsite because mutating residues that interact with the phosphate moiety of FMN in the native complex hardly changes the observed rate constants and amplitudes. In contrast, replacing tyr92, which interacts with the isoalloxazine, greatly lowers the rate constants. Our analysis indicates that the two FMN binding phases observed are related neither with alternative or sequential interaction with the two binding subsites nor with the presence of bound phosphate. It is possible that they reflect the intrinsic conformational heterogeneity of the apoflavodoxin ensemble.