Bacterial Luciferases from Vibrio harveyi and Photobacterium leiognathi Demonstrate Different Conformational Stability as Detected by Time-Resolved Fluorescence Spectroscopy.

Detecting the folding/unfolding pathways of biological macromolecules is one of the urgent problems of molecular biophysics. The unfolding of bacterial luciferase from Vibrio harveyi is well-studied, unlike that of Photobacterium leiognathi, despite the fact that both of them are actively used as a reporter system. The aim of this study was to compare the conformational transitions of these luciferases from two different protein subfamilies during equilibrium unfolding with urea. Intrinsic steady-state and time-resolved fluorescence spectra and circular dichroism spectra were used to determine the stages of the protein unfolding. Molecular dynamics methods were applied to find the differences in the surroundings of tryptophans in both luciferases. We found that the unfolding pathway is the same for the studied luciferases. However, the results obtained indicate more stable tertiary and secondary structures of P. leiognathi luciferase as compared to enzyme from V. harveyi during the last stage of denaturation, including the unfolding of individual subunits. The distinctions in fluorescence of the two proteins are associated with differences in the structure of the C-terminal domain of α-subunits, which causes different quenching of tryptophan emissions. The time-resolved fluorescence technique proved to be a more effective method for studying protein unfolding than steady-state methods.

Right- and left-handed three-helix proteins. II. Similarity and differences in mechanical unfolding of proteins.

Here, we study mechanical properties of eight 3-helix proteins (four right-handed and four left-handed ones), which are similar in size under stretching at a constant speed and at a constant force on the atomic level using molecular dynamics simulations. The analysis of 256 trajectories from molecular dynamics simulations with explicit water showed that the right-handed three-helix domains are more mechanically resistant than the left-handed domains. Such results are observed at different extension velocities studied (192 trajectories obtained at the following conditions: v = 0.1, 0.05, and 0.01 Å ps(-1) , T = 300 K) and under constant stretching force (64 trajectories, F = 800 pN, T = 300 K). We can explain this by the fact, at least in part, that the right-handed domains have a larger number of contacts per residue and the radius of cross section than the left-handed domains.

All-atom and coarse-grained simulations of the forced unfolding pathways of the SNARE complex.

The SNARE complex, consisting of three proteins (VAMP2, syntaxin, and SNAP-25), is thought to drive membrane fusion by assembling into a four-helix bundle through a zippering process. In support of the above zippering model, a recent single-molecule optical tweezers experiment by Gao et al. revealed a sequential unzipping of SNARE along VAMP2 in the order of the linker domain → the C-terminal domain → the N-terminal domain. To offer detailed structural insights to this unzipping process, we have performed all-atom and coarse-grained steered molecular dynamics (sMD) simulations of the forced unfolding pathways of SNARE using different models and force fields. Our findings are summarized as follows: First, the sMD simulations based on either an all-atom force field (with an implicit solvent model) or a coarse-grained Go model were unable to capture the forced unfolding pathway of SNARE as observed by Gao et al., which may be attributed to insufficient simulation time and inaccurate force fields. Second, the sMD simulations based on a reparameterized coarse-grained model (i.e., modified elastic network model) were able to predict a sequential unzipping of SNARE in good agreement with the findings by Gao et al. The key to this success is to reparameterize the intrahelix and interhelix nonbonded force constants against the pair-wise residue-residue distance fluctuations collected from all-atom MD simulations of SNARE. Therefore, our finding supports the importance of accurately describing the inherent dynamics/flexibility of SNARE (in the absence of force), in order to correctly simulate its unfolding behaviors under force. This study has established a useful computational framework for future studies of the zippering function of SNARE and its perturbations by point mutations with amino-acid level of details, and more generally the forced unfolding pathways of other helix bundle proteins.

Role of the copper ion in pseudoazurin during the mechanical unfolding process.

Metalloproteins require the corresponding metal cofactors to exert their proper function. The presence of metal cofactors in the metalloprotein makes it more difficult to investigate its folding and unfolding process. In this study, we employed atomic-force-microscopy-based single-molecule force spectroscopy to reveal the unfolding process of pseudoazurin (PAZ) that belongs to blue copper proteins. Our study shows that holo-PAZ requires a higher rupture force for mechanical unfolding comparing with the apo-PAZ. This result demonstrates that the copper atom not only enables PAZ access to transfer electron, but should also have an influence on its stability. The results also suggest that the electronic configuration of the metal cofactors has a striking effect on the strength of the organometallic bonds. Moreover, the results also reveal that there is an intermediate state during the unfolding process of PAZ. This study provides insight into the characteristics of metalloproteins and leads to a better knowledge of their interaction at the individual molecule level.

Mechanism for the Unfolding of the TOP7 Protein in Steered Molecular Dynamics Simulations as Revealed by Mutual Information Analysis.

We employed mutual information (MI) analysis to detect motions affecting the mechanical resistance of the human-engineered protein Top7. The results are based on the MI analysis of pair contact correlations measured in steered molecular dynamics (SMD) trajectories and their statistical dependence on global unfolding. This study is the first application of the MI analysis to SMD forced unfolding, and we furnish specific SMD recommendations for the utility of parameters and options in the TimeScapes package. The MI analysis provided a global overview of the effect of perturbation on the stability of the protein. We also employed a more conventional trajectory analysis for a detailed description of the mechanical resistance of Top7. Specifically, we investigated 1) the hydropathy of the interactions of structural segments, 2) the H2O concentration near residues relevant for unfolding, and 3) the changing hydrogen bonding patterns and main chain dihedral angles. The results show that the application of MI in the study of protein mechanical resistance can be useful for the engineering of more resistant mutants when combined with conventional analysis. We propose a novel mutation design based on the hydropathy of residues that would stabilize the unfolding region by mimicking its more stable symmetry mate. The proposed design process does not involve the introduction of covalent crosslinks, so it has the potential to preserve the conformational space and unfolding pathway of the protein.

Loading-device effects on the protein-unfolding mechanisms using molecular-dynamic simulations.

Experimental force spectroscopy has been effectively utilized for measuring structural characterization of biomolecules and mechanical properties of biomaterials. Specifically, atomic force microscopy (AFM) has been widely used to portray biomolecular characterization in single-molecule experiment by observing the unfolding behavior of the proteins. Not only the experimental techniques enable us to characterize globular protein, but computational methods like molecular dynamics (MD) also gives insight into understanding biomolecular structures. To better comprehend the behavior of biomolecules, conditions such as pulling velocities and loading rates are put to the test, yet there are still limitations in understanding the unfolding behavior of biomolecules with the effect of different loading devices. In this study, we performed an all-atom MD and steered molecular dynamics (SMD) simulations considering different loading device effects such as "soft" and "stiff" to characterize the anisotropic unfolding behavior of ubiquitin protein. We found out the anisotropic unfolding pathways of the protein through the broken number of hydrogen bonds and geometric secondary structures of the biomolecule. Our study provides the importance for usage of various loading-devices on biomolecules when analyzing the structural compositions and the characteristics of globular biomolecules.

Investigating the ion dependence of the first unfolding step of GTPase-Associating Center ribosomal RNA.

The interactions in the tertiary structure of a ribosomal RNA fragment in the GTPase Associating Center (GAC) have been experimentally studied, but the roles of the bound and diffuse cations in its folding pathway have not yet been fully elucidated. Melting experiments have shown that the temperature of the first of the two distinguishable transitions in the unfolding pathway of the GAC RNA can be regulated by altering the magnesium concentration, yet the physical interpretation of such ion-dependent effects on folding have not been clearly understood in spite of the availability of crystal structures that depict many GAC RNA-ion interactions. Here, we use umbrella sampling and molecular dynamics (MD) simulations to provide a physical description for the first transition in this unfolding pathway, with a focus on the role of a chelated magnesium ion. Our results indicate that the presence of cations mediating the local interaction of two loops stabilizes the folded state relative to the unfolded or partially folded states. Also, our findings suggest that a bridging magnesium ion between the two loops improves the stabilizing effect. This is consistent with the multistep unfolding pathway proposed for the GAC RNA and highlights the importance of ions in the first unfolding step. The results suggest how MD simulations can provide insight into RNA unfolding pathways as a complementary approach to experiments.

Thermal stability and unfolding pathways of hyperthermophilic and mesophilic periplasmic binding proteins studied by molecular dynamics simulation.

The ribose binding protein (RBP), a sugar-binding periplasmic protein, is involved in the transport and signaling processes in both prokaryotes and eukaryotes. Although several cellular and structural studies have been reported, a description of the thermostability of RBP at the molecular level remains elusive. Focused on the hyperthermophilic Thermoytoga maritima RBP (tmRBP) and mesophilic Escherichia coli homolog (ecRBP), we applied molecular dynamics simulations at four different temperatures (300, 380, 450, and 500 K) to obtain a deeper insight into the structural features responsible for the reduced thermostability of the ecRBP. The simulations results indicate that there are distinct structural differences in the unfolding pathway between the two homologs and the ecRBP unfolds faster than the hyperthermophilic homologs at certain temperatures in accordance with the lower thermal stability found experimentally. Essential dynamics analysis uncovers that the essential subspaces of ecRBP and tmRBP are non-overlapping and these two proteins show different directions of motion within the simulations trajectories. Such an understanding is required for designing efficient proteins with characteristics for a particular application.

Multiple Unfolding Intermediates Obtained by Molecular Dynamic Simulations under Stretching for Immunoglobulin-Binding Domain of Protein G.

We have studied the mechanical properties of the immunoglobulin-binding domain of protein G at the atomic level under stretching at constant velocity using molecular dynamics simulations. We have found that the unfolding process can occur either in a single step or through intermediate states. Analysis of the trajectories from the molecular dynamic simulations showed that the mechanical unfolding of the immunoglobulin-binding domain of protein G is triggered by the separation of the terminal beta-strands and the order in which the secondary-structure elements break is practically the same in two- and multi-state events and at the different extension velocities studied. It is seen from our analysis of 24 trajectories that the theoretical pathway of mechanical unfolding for the immunoglobulin-binding domain of protein G does not coincide with that proposed in denaturant studies in the absence of force.