Exploring Compactness and Dynamics of Apomyoglobin

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) approach has become a valuable analytical complement to traditional methods. HDX-MS allows the identification of dynamic surfaces in proteins. We have shown that the introduction of various mutations into the amino acid sequence of whale apomyoglobin (apoMb) leads to a change in the number of exchangeable hydrogen atoms, which is associated with a change in its compactness in the native-like condition. Thus, amino acid substitutions V10A, A15S, P120G, and M131A result in an increase in the number of exchangeable hydrogen atoms at the native-like condition, while the mutant form A144S leads to a decrease in the number of exchangeable hydrogen atoms. This may be due to a decrease and increase in the compactness of apoMb structure compared to the wild-type apoMb, respectively. The L9F and L9E mutations did not affect the compactness of the molecule compared to the wild type. We have demonstrated that V10A and M131A substitutions lead to the maximum and large increase correspondently in the average number of exchangeable hydrogen atoms for deuterium, since these substitutions lead to the loss of contacts between important parts of myoglobin structure: helices A, G, and H, which are structured at the early stage of folding.

Computational Modeling of Molecular Structures Guided by Hydrogen-Exchange Data

Data produced by hydrogen-exchange monitoring experiments have been used in structural studies of molecules for several decades. Despite uncertainties about the structural determinants of hydrogen exchange itself, such data have successfully helped guide the structural modeling of challenging molecular systems, such as membrane proteins or large macromolecular complexes. As hydrogen-exchange monitoring provides information on the dynamics of molecules in solution, it can complement other experimental techniques in so-called integrative modeling approaches. However, hydrogen-exchange data have often only been used to qualitatively assess molecular structures produced by computational modeling tools. In this paper, we look beyond qualitative approaches and survey the various paradigms under which hydrogen-exchange data have been used to quantitatively guide the computational modeling of molecular structures. Although numerous prediction models have been proposed to link molecular structure and hydrogen exchange, none of them has been widely accepted by the structural biology community. Here, we present as many hydrogen-exchange prediction models as we could find in the literature, with the aim of providing the first exhaustive list of its kind. From purely structure-based models to so-called fractional-population models or knowledge-based models, the field is quite vast. We aspire for this paper to become a resource for practitioners to gain a broader perspective on the field and guide research toward the definition of better prediction models. This will eventually improve synergies between hydrogen-exchange monitoring and molecular modeling.

Supramolecular organization of Cytochrome-C into quantum-dot decorated macromolecular network under pH and thermal stress

The holo form of Cytochrome-C which is involved in the electron transfer chain of aerobic and anaerobic respiration remains structurally intact by its complex with heme. However, when a prolonged thermal and pH stress was applied, heme was found to abruptly dissociate from the holo protein, resulting in complete collapse of the three-dimensional functional structure. Interestingly, two distinct structures were formed as the consequence of the dissociation event: (i) A macromolecular amyloid-network formed by the collapsed protein fragments, generated by self-oxidation, and (ii) Fe-containing Quantum-Dots (FeQDs) with 2-3 nm diameter formed by heme reorganization. Further adding to intrigue, the FeQDs were re-adsorbed on the surface of the amyloid network leading to FeQD-decorated macromolecular amyloid matrix. The heme-interactant Met80, constituting the amyloidogenic region, initiates the amylogenic cascade, and gradual exposure of Trp59 synergistically emit intrinsic fluorescence alongside FeQDs. The development of the aforementioned events were probed through a multitude of biophysical, chemical and computational analyses like ThT/ANS/intrinsic fluorescence assays, CD-spectroscopy, FETEM/STEM/elemental mapping, Foldamyloid/Foldunfold/Isunstruct/H-protection/LIGplot analyses, etc. The FeQD-decorated amyloid-network was found to exhibit gel-like property, which supported the growth of BHK-21 fibroblast without cytotoxicity. Further studies on FeQD-decorated Cytochrome C amyloid network might open possibilities to design advanced biomaterial for diverse biological applications.

A Review on Computational Approaches for Analyzing Hydrogen- Deuterium (H/D) Exchange of Proteins

Native state Hydrogen-Deuterium (H/D) exchange method has been used to study the structures and the unfolding pathways for quite a number of proteins. The H/D exchange method is generally monitored using nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS) techniques. NMR-assisted H/D exchange methods primarily monitor the residue level fluctuation of proteins, whereas MS-assisted H/D exchange methods analyze multifold ensemble conformations of proteins. In this connection, quite a large number of computational tools and algorithms have been developed for processing and analyzing huge amount of the H/D exchange data generated from these techniques. In this review, most of the freely available computational tools associated with the H/D exchange of proteins have been comprehensively reviewed and scopes to improve/ develop novel computational approaches for analyzing the H/D exchange data of proteins have also been brought into fore.

Perdeuteration, large crystal growth and neutron data collection of Leishmania mexicana triose-phosphate isomerase E65Q variant

Triose-phosphate isomerase (TIM) catalyses the interconversion of dihydroxyacetone phosphate and glyceraldehyde 3-phosphate. Two catalytic mechanisms have been proposed based on two reaction-intermediate analogues, 2-phosphoglycolate (2PG) and phosphoglycolohydroxamate (PGH), that have been used as mimics of the cis-enediol(ate) intermediate in several studies of TIM. The protonation states that are critical for the mechanistic interpretation of these structures are generally not visible in the X-ray structures. To resolve these questions, it is necessary to determine the hydrogen positions using neutron crystallography. Neutron crystallography requires large crystals and benefits from replacing all hydrogens with deuterium. Leishmania mexicana triose-phosphate isomerase was therefore perdeuterated and large crystals with 2PG and PGH were produced. Neutron diffraction data collected from two crystals with different volumes highlighted the importance of crystal volume, as smaller crystals required longer exposures and resulted in overall worse statistics.

Exploring the Sequence-based Prediction of Folding Initiation Sites in Proteins

Protein folding is a complex process that can lead to disease when it fails. Especially poorly understood are the very early stages of protein folding, which are likely defined by intrinsic local interactions between amino acids close to each other in the protein sequence. We here present EFoldMine, a method that predicts, from the primary amino acid sequence of a protein, which amino acids are likely involved in early folding events. The method is based on early folding data from hydrogen deuterium exchange (HDX) data from NMR pulsed labelling experiments, and uses backbone and sidechain dynamics as well as secondary structure propensities as features. The EFoldMine predictions give insights into the folding process, as illustrated by a qualitative comparison with independent experimental observations. Furthermore, on a quantitative proteome scale, the predicted early folding residues tend to become the residues that interact the most in the folded structure, and they are often residues that display evolutionary covariation. The connection of the EFoldMine predictions with both folding pathway data and the folded protein structure suggests that the initial statistical behavior of the protein chain with respect to local structure formation has a lasting effect on its subsequent states.

Coarse-Grained Conformational Sampling of Protein Structure Improves the Fit to Experimental Hydrogen-Exchange Data

Monitoring hydrogen/deuterium exchange (HDX) undergone by a protein in solution produces experimental data that translates into valuable information about the protein's structure. Data produced by HDX experiments is often interpreted using a crystal structure of the protein, when available. However, it has been shown that the correspondence between experimental HDX data and crystal structures is often not satisfactory. This creates difficulties when trying to perform a structural analysis of the HDX data. In this paper, we evaluate several strategies to obtain a conformation providing a good fit to the experimental HDX data, which is a premise of an accurate structural analysis. We show that performing molecular dynamics simulations can be inadequate to obtain such conformations, and we propose a novel methodology involving a coarse-grained conformational sampling approach instead. By extensively exploring the intrinsic flexibility of a protein with this approach, we produce a conformational ensemble from which we extract a single conformation providing a good fit to the experimental HDX data. We successfully demonstrate the applicability of our method to four small and medium-sized proteins.

A Residue-Resolved Bayesian Approach to Quantitative Interpretation of Hydrogen-Deuterium Exchange from Mass Spectrometry: Application to Characterizing Protein-Ligand Interactions

Characterization of interactions between proteins and other molecules is crucial for understanding the mechanisms of action of biological systems and, thus, drug discovery. An increasingly useful approach to mapping these interactions is measurement of hydrogen/deuterium exchange (HDX) using mass spectrometry (HDX-MS), which measures the time-resolved deuterium incorporation of peptides obtained by enzymatic digestion of the protein. Comparison of exchange rates between apo- and ligand-bound conditions results in a mapping of the differential HDX (ΔHDX) of the ligand. Residue-level analysis of these data, however, must account for experimental error, sparseness, and ambiguity due to overlapping peptides. Here, we propose a Bayesian method consisting of a forward model, noise model, prior probabilities, and a Monte Carlo sampling scheme. This method exploits a residue-resolved exponential rate model of HDX-MS data obtained from all peptides simultaneously, and explicitly models experimental error. The result is the best possible estimate of ΔHDX magnitude and significance for each residue given the data. We demonstrate the method by revealing richer structural interpretation of ΔHDX data on two nuclear receptors: vitamin D-receptor (VDR) and retinoic acid receptor gamma (RORγ). The method is implemented in HDX Workbench and as a standalone module of the open source Integrative Modeling Platform.

Characterization of folding cores in the cyclophilin A-cyclosporin A complex

Determining the folding core of a protein yields information about its folding process and dynamics. The experimental procedures for identifying the amino acids that make up the folding core include hydrogen-deuterium exchange and Φ-value analysis and can be expensive and time consuming. Because of this, there is a desire to improve upon existing methods for determining protein folding cores theoretically. We have obtained HDX data for the complex of cyclophilin A with the immunosuppressant cyclosporin A. We compare these data, as well as literature values for uncomplexed cyclophilin A, to theoretical predictions using a combination of rigidity analysis and coarse-grained simulations of protein motion. We find that in this case, the most specific prediction of folding cores comes from a combined approach that models the rigidity of the protein using the first software suite and the dynamics of the protein using the froda tool.

Functional dynamics of hexameric helicase probed by hydrogen exchange and simulation

The biological function of large macromolecular assemblies depends on their structure and their dynamics over a broad range of timescales; for this reason, it is a significant challenge to investigate these assemblies using conventional experimental techniques. One of the most promising experimental techniques is hydrogen-deuterium exchange detected by mass spectrometry. Here, we describe to our knowledge a new computational method for quantitative interpretation of deuterium exchange kinetics and apply it to a hexameric viral helicase P4 that unwinds and translocates RNA into a virus capsid at the expense of ATP hydrolysis. Room-temperature dynamics probed by a hundred nanoseconds of all-atom molecular dynamics simulations is sufficient to predict the exchange kinetics of most sequence fragments and provide a residue-level interpretation of the low-resolution experimental results. The strategy presented here is also a valuable tool to validate experimental data, e.g., assignments, and to probe mechanisms that cannot be observed by x-ray crystallography, or that occur over timescales longer than those that can be realistically simulated, such as the opening of the hexameric ring.

Association between foldability and aggregation propensity in small disulfide-rich proteins

AIMS

Disulfide-rich domains (DRDs) are small proteins whose native structure is stabilized by the presence of covalent disulfide bonds. These domains are versatile and can perform a wide range of functions. Many of these domains readily unfold on disulfide bond reduction, suggesting that in the absence of covalent bonding they might display significant disorder.

RESULTS

Here, we analyzed the degree of disorder in 97 domains representative of the different DRDs families and demonstrate that, in terms of sequence, many of them can be classified as intrinsically disordered proteins (IDPs) or contain predicted disordered regions. The analysis of the aggregation propensity of these domains indicates that, similar to IDPs, their sequences are more soluble and have less aggregating regions than those of other globular domains, suggesting that they might have evolved to avoid aggregation after protein synthesis and before they can attain its compact and covalently linked native structure.

INNOVATION AND CONCLUSION

DRDs, which resemble IDPs in the reduced state and become globular when their disulfide bonds are formed, illustrate the link between protein folding and aggregation propensities and how these two properties cannot be easily dissociated, determining the main traits of the folding routes followed by these small proteins to attain their native oxidized states.

OneG-Vali: a computational tool for detecting, estimating and validating cryptic intermediates of proteins under native conditions

Unfolding pathway of T4 lysozyme under native conditions as predicted by the OneG-Vali has been illustrated. Also, structural contexts of various states (native (N), cryptic intermediates (CIs) and unfolded (U) conformations) of the protein and the population of three CIs are depicted.