Fast collapse but slow formation of secondary structure elements in the refolding transition of E. coli adenylate kinase

The various models proposed for protein folding transition differ in their order of appearance of the basic steps during this process. In this study, steady state and time-resolved dynamic non-radiative excitation energy transfer (FRET and trFRET) combined with site specific labeling experiments were applied in order to characterize the initial transient ensemble of Escherichia coli adenylate kinase (AK) molecules upon shifting conditions from those favoring denaturation to refolding and from folding to denaturing. Three sets of labeled AK mutants were prepared, which were designed to probe the equilibrium and transient distributions of intramolecular segmental end-to-end distances. A 176 residue section (residues 28-203), which spans most of the 214 residue molecule, and two short secondary structure chain segments including an alpha-helix (residues 169-188) and a predominantly beta-strand region (residues 188-203), were labeled. Upon fast change of conditions from denaturing to folding, the end-to-end distance of the 176 residue chain section showed an immediate collapse to a mean value of 26 A. Under the same conditions, the two short secondary structure elements did not respond to this shift within the first ten milliseconds, and retained the characteristics of a fully unfolded state. Within the first 10 ms after changes of the solvent from folding to denaturing, only minor changes were observed at the local environments of residues 203 and 169. The response of these same local environments to the shift of conditions from denaturing to folding occurred within the dead time of the mixing device. Thus, the response of the CORE domain of AK to fast transfer from folding to unfolding conditions is slow at all three conformational levels that were probed, and for at least a few milliseconds the ensemble of folded molecules is maintained under unfolding conditions. A different order of the changes was observed upon initiation of refolding. The AK molecules undergo fast collapse to an ensemble of compact structures where the local environment of surface probes seems to be native-like but the two labeled secondary structure elements remain unfolded.     

The effect of distribution of adjacent amino acids in primary structures on the possibility of mutation substitutions in globular proteins

Frequency distributions of adjacent ARs in the primary structures of 320 globular proteins out of different superfamilies were investigated. ARs of every type were compared with the occurrence frequencies of 20 canonic residues at the distances of 1-20 residues according to their primary structure. Amino acid residues were found to be divided into groups of interchangeable residues in the course of globular protein evolution according to the distribution kinds and in terms of Euclidean distances. The use of pancreatic RNases of mammals showed that the approximate preservation of frequency adjacent (1-4 residues according to their primary structure) and characteristics in 5-15 residues mid-interactions may be used in studying the supposed amino acid substitutions in globular protein.     

Fast Subdomain Folding Prior to the Global Refolding Transition of E. coli Adenylate Kinase: A Double Kinetics Study

The rate of folding of globular proteins depends on specific local and nonlocal intramolecular interactions. What is the relative role of these two types of interaction at the initiation of refolding? We address this question by application of a “double kinetics” method based on fast initiation of refolding of site specifically labeled protein samples and detection of the transient distributions of selected intramolecular distances by means of fast measurements of time‐resolved fluorescence resonance energy transfer. We determined the distribution of the distance between the ends of a 44‐chain segment that includes the AMP_bind domain, by labeling residues 28 and 71, in Escherichia coli adenylate kinase (AK) and the distribution of the distance between residues 18 and 203, which depends on the overall order of the molecule. That distribution shows two-state transition to the native intramolecular distance at the same rate as that of the cooperative refolding transition of the AK molecule. In sharp contrast, the distance distribution between residues 28 and 71 is already native like at the end of the dead-time of the mixing device. This fast formation of native short distance between two widely separated chain sections can be either dependent on fast folding of the AMP_bind domain or a result of a very effective nonlocal interaction between specific short clusters of hydrophobic residues. Further experiments on studying the kinetics of folding of selected structural elements in the protein will help determination of the driving force of this early folding event.     

Fluctuations in ion pairs and their stabilities in proteins

This report investigates the effect of systemic protein conformational flexibility on the contribution of ion pairs to protein stability. Toward this goal, we use all NMR conformer ensembles in the Protein Data Bank (1) that contain at least 40 conformers, (2) whose functional form is monomeric, (3) that are nonredundant, and (4) that are large enough. We find 11 proteins adhering to these criteria. Within these proteins, we identify 22 ion pairs that are close enough to be classified as salt bridges. These are identified in the high-resolution crystal structures of the respective proteins or in the minimized average structures (if the crystal structures are unavailable) or, if both are unavailable, in the "most representative" conformer of each of the ensembles. We next calculate the electrostatic contribution of each such ion pair in each of the conformers in the ensembles. This results in a comprehensive study of 1,201 ion pairs, which allows us to look for consistent trends in their electrostatic contributions to protein stability in large sets of conformers. We find that the contributions of ion pairs vary considerably among the conformers of each protein. The vast majority of the ion pairs interconvert between being stabilizing and destabilizing to the structure at least once in the ensembles. These fluctuations reflect the variabilities in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other, and with respect to other charged groups in the remainder of the protein. The higher crystallographic B-factors for the respective side-chains are consistent with these fluctuations. The major conclusion from this study is that salt bridges observed in crystal structure may break, and new salt bridges may be formed. Hence, the overall stabilizing (or, destabilizing) contribution of an ion pair is conformer population dependent.     

Distribution of amino acid residues in the primary protein structure

The amino acid composition and sequence in primary structure of 180 proteins have been studied. It is shown that the distribution of amino acid residues is near to a random one, i.e. it is determined by the amino acid composition. The ratio between statistical and unique character of protein primary structures has been discussed. The amino acid sequence is suggested to be unique in fibrous proteins. In contrast the amino acid sequence in globular proteins is a statistical one. The statistical character of amino acids distribution in globular proteins explains the possibility of sensible text generation under the frame shift mutations, deletions and insertions.     

Calorimetric and structural investigation of the interaction between bovine serum albumin and high molecular weight dextran in water

This work studies specific interactions between a small globular protein and a highly flexible, branched polysaccharide using differential scanning calorimetry (DSC), circular dichroism (CD), fluorescence, and turbidimetry measurements. It uses the system water/bovine serum albumin (BSA)/dextran (D 2000) as a model. Dextran molecules are able to form interpolymeric complexes with BSA in water at both low and high temperatures if the polysaccharide is in excess and if the protein exists in its associated state. It leads to a partial destabilization of the secondary and tertiary structures of the protein and an additional exposure of the hydrophobic tryptophan residues to the surface of globule. If the total concentration of biopolymers in the mixture is high enough, the stability of the protein molecules with respect to unfolding and thermoaggregation is significantly decreased as a result of an increase in the protein hydrophobicity.     

Heme as an optical probe of a conformational transition of ovine recPrP

The prion protein occurs as a globular domain and a leading fragment whose structure is not well-defined. For the ovine species, all of the tryptophan residues are in the initial fragment, while the globular domain is rich in tyrosine residues. Using heme as a spectroscopic probe, we have studied the recombinant prion protein before and after a temperature-induced conformational change. As for most heme proteins, the absorption spectrum of heme-CO displays a red shift upon binding to the protein, and both the Y and W fluorescence are highly quenched. Flash photolysis kinetics of the PrP-heme-CO complex shows a low yield for the bimolecular phase, indicating a pocket around the hemes. By comparing the holoprotein and the truncated sequence corresponding to the globular domain, the stoichiometry was determined to be five hemes for the globular domain and two hemes for the leading fragment. At high temperature, the hemes are released; upon cooling, only two hemes bind, and only the tryptophan fluorescence is quenched; this would indicate that the globular domain has formed a more compact structure, which is inert with respect to the hydrophobic probe. The final state of polymerization is perturbed if the synthetic peptide "N3" (PrP residues 142-166, which include the first helix) is added to the prion protein solution; the temperature cycle no longer reduces the number of heme binding sites. This would indicate that the peptide may alter or inhibit the polymer formation.     

The natively helical chain segment 169-188 of Escherichia coli adenylate kinase is formed in the latest phase of the refolding transition

The refolding transition of Escherichia coli adenylate kinase (AK) was investigated by monitoring the refolding kinetics of a selected 20 residue helical segment in the CORE domain of the protein. Residues 169 and 188 were labeled by 1-acetamido-methyl-pyrene, and by bimane, respectively. The experiment combines double-jump stopped-flow fast mixing initiation of refolding and time-resolved F?rster energy transfer spectroscopy for monitoring the conformational transitions (double-kinetics experiment). Two kinetic phases were found in the denaturant-induced unfolding of AK. In the first phase, the fluorescence quantum yields of both probes decreased. The distribution of the distances between them transformed from the native state's narrow distribution with the mean distance corresponding to the distance in the crystal structure, to a distribution compatible with an unordered structure. In the second, slow step of denaturation, neither the fluorescence parameters of the probes nor the distance distribution between them changed. This step appeared to be a transformation of the fast-folding species formed in the first phase, to the slow-folding species. Refolding of the fast-folding species of the denatured state of AK was also a two-phase process. During the first fast phase, within less than 5ms, the fluorescence emission of both probes increased, but the distance distribution between the labeled sites was unchanged. Only during the second slow refolding step did the intramolecular distance distribution change from the characteristic of the denatured state to the narrow distribution of the native state. This experiment shows that for the case of the CORE domain of AK, the large helical segment of residues 169-188 was not formed in the first compaction step of refolding. The helical conformation of this segment is established only in the second, much slower, refolding phase, simultaneously with the completion of the native structure.     

Structure Prediction of the Second Extracellular Loop in G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) play key roles in living organisms. Therefore, it is important to determine their functional structures. The second extracellular loop (ECL2) is a functionally important region of GPCRs, which poses significant challenge for computational structure prediction methods. In this work, we evaluated CABS, a well-established protein modeling tool for predicting ECL2 structure in 13 GPCRs. The ECL2s (with between 13 and 34 residues) are predicted in an environment of other extracellular loops being fully flexible and the transmembrane domain fixed in its x-ray conformation. The modeling procedure used theoretical predictions of ECL2 secondary structure and experimental constraints on disulfide bridges. Our approach yielded ensembles of low-energy conformers and the most populated conformers that contained models close to the available x-ray structures. The level of similarity between the predicted models and x-ray structures is comparable to that of other state-of-the-art computational methods. Our results extend other studies by including newly crystallized GPCRs.     

Comparison of the transition state ensembles for folding of Im7 and Im9 determined using all-atom molecular dynamics simulations with phi value restraints

Delineation of the structural properties of transition states is key to deriving models for protein folding. Here we describe the structures of the transition states of the bacterial immunity proteins Im7 and Im9 obtained by all-atom molecular dynamics simulations with phi value restraints derived from protein engineering experiments. This pair of proteins is of special interest because, at pH 7 and 10 degrees C, Im7 folds via an intermediate while Im9 folds with a two-state transition. The structures of the transition states for Im7 and Im9, together with their radii of gyration and distances from the native state, are similar. The typical distance between any two members of the transition state ensemble of both proteins is large, with that of Im9 nearly twice that of Im7. Thus, a broad range of structures make up the transition state ensembles of these proteins. The ensembles satisfy the set of rather low phi values and yet are consistent with high beta(T) values (> 0.85 for both proteins). For both Im7 and Im9 the inter-helical angles are highly variable in the transition state ensembles, although the native contacts between helices I and IV are well conserved. By measuring the distribution of the accessible surface area for each residue we show that the hydrophobic residues that are buried in the native state remain buried in the transition state, corresponding to a hydrophobic collapse to a relatively ordered globule. The data provide new insights into the structural properties of the transition states of these proteins at an atomic level of detail and show that molecular dynamics simulations with phi value restraints can significantly enhance the knowledge of the transition state ensembles (TSE) provided by the experimental phi values alone.     

Many Local Motions Cooperate to Produce the Adenylate Kinase Conformational Transition

Conformational transitions are functionally important in many proteins. In the enzyme adenylate kinase (AK), two small domains (LID and NMP) close over the larger CORE domain; the reverse (opening) motion limits the rate of catalytic turnover. Here, using double-well Gō simulations of Escherichia coli AK, we elaborate on previous investigations of the AK transition mechanism by characterizing the contributions of rigid-body (Cartesian), backbone dihedral, and contact motions to transition-state (TS) properties. In addition, we compare an apo simulation to a pseudo-ligand-bound simulation to reveal insights into allostery. In Cartesian space, LID closure precedes NMP closure in the bound simulation, consistent with prior coarse-grained models of the AK transition. However, NMP-first closure is preferred in the apo simulation. In backbone dihedral space, we find that, as expected, backbone fluctuations are reduced in the O/C transition in parts of all three domains. Among these “quenching” residues, most in the CORE, especially residues 11–13, are rigidified in the TS of the bound simulation, while residues 42–44 in the NMP are flexible in the TS. In contact space, in both apo and bound simulations, one nucleus of closed-state contacts includes parts of the NMP and CORE; CORE–LID contacts are absent in the TS of the apo simulation but formed in the TS of the bound simulation. From these results, we predict mutations that will perturb the opening and/or closing transition rates by changing the entropy of dihedrals and/or the enthalpy of contacts. Furthermore, regarding allostery, the fully closed structure is populated in the apo simulation, but our contact results imply that ligand binding shifts the preferred O/C transition pathway, thus precluding a simple conformational selection mechanism. Finally, the analytical approach and the insights derived from this work may inform the rational design of flexibility and allostery in proteins.     

Structure Prediction of the Second Extracellular Loop in G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) play key roles in living organisms. Therefore, it is important to determine their functional structures. The second extracellular loop (ECL2) is a functionally important region of GPCRs, which poses significant challenge for computational structure prediction methods. In this work, we evaluated CABS, a well-established protein modeling tool for predicting ECL2 structure in 13 GPCRs. The ECL2s (with between 13 and 34 residues) are predicted in an environment of other extracellular loops being fully flexible and the transmembrane domain fixed in its x-ray conformation. The modeling procedure used theoretical predictions of ECL2 secondary structure and experimental constraints on disulfide bridges. Our approach yielded ensembles of low-energy conformers and the most populated conformers that contained models close to the available x-ray structures. The level of similarity between the predicted models and x-ray structures is comparable to that of other state-of-the-art computational methods. Our results extend other studies by including newly crystallized GPCRs.     

Distribution of amino acid residues and residue-residue contacts in molecular chaperones.

The amino acid distribution and residue-residue contacts in molecular chaperones are different when compared to normal globular proteins. The study of molecular chaperones reveals a different surrounding environment to exist for the residues Cys, Trp, and His which may play an important role in determining the chaperone structures. Unlike globular proteins, it has been observed that a one-to-one correspondence between the amino acid distribution in a sequence and the structures of molecular chaperones. The preference of amino acid residues surrounding all 20 types of residues in secondary structures and their accessible surface areas have been analysed.     

Selective observation of the disordered import signal of a globular protein by in-cell NMR: The example of frataxins

We have exploited the capability of in-cell NMR to selectively observe flexible regions within folded proteins to carry out a comparative study of two members of the highly conserved frataxin family which are found both in prokaryotes and in eukaryotes. They all contain a globular domain which shares more than 50% identity, which in eukaryotes is preceded by an N-terminal tail containing the mitochondrial import signal. We demonstrate that the NMR spectrum of the bacterial ortholog CyaY cannot be observed in the homologous E. coli system, although it becomes fully observable as soon as the cells are lysed. This behavior has been observed for several other compact globular proteins as seems to be the rule rather than the exception. The NMR spectrum of the yeast ortholog Yfh1 contains instead visible signals from the protein. We demonstrate that they correspond to the flexible N-terminal tail indicating that this is flexible and unfolded. This flexibility of the N-terminus agrees with previous studies of human frataxin, despite the extensive sequence diversity of this region in the two proteins. Interestingly, the residues that we observe in in-cell experiments are not visible in the crystal structure of a Yfh1 mutant designed to destabilize the first helix. More importantly, our results show that, in cell, the protein is predominantly present not as an aggregate but as a monomeric species.     

Prediction of the location of structural domains in globular proteins

The location of structural domains in proteins is predicted from the amino acid sequence, based on the analysis of a computed contact map for the protein, the average distance map (ADM). Interactions between residues i and j in a protein are subdivided into several ranges, according to the separation |i-j| in the amino acid sequence. Within each range, average spatial distances between every pair of amino acid residues are computed from a data base of known protein structures. Infrequently occurring pairs are omitted as being statistically insignificant. The average distances are used to construct a predicted ADM. The ADM is analyzed for the occurrence of regions with high densities of contacts (compact regions). Locations of rapid changes of density between various parts of the map are determined by the use of scanning plots of contact densities. These locations serve to pinpoint the distribution of compact regions. This distribution, in turn, is used to predict boundaries of domains in the protein. The technique provides an objective method for the location of domains both on a contact map derived from a known three-dimensional protein structure, the real distance map (RDM), and on an ADM. While most other published methods for the identification of domains locate them in the known three-dimensional structure of a protein, the technique presented here also permits the prediction of domains in proteins of unknown spatial structure, as the construction of the ADM for a given protein requires knowledge of only its amino acid sequence.     

Conformations of folded proteins in restricted spaces

A new method is presented to examine the complete range of folded topologies accessible in the compact state of globular proteins. The procedure is to generate all conformations, with volume exclusion, upon a lattice in a space restricted to the individual protein's known compact conformational space. Using one lattice point per residue, we find 10(2)-10(4) possible compact conformations for the five small globular proteins studied. Subsequently, these conformations are evaluated in terms of residue-specific, pairwise contact energies that favor nonbonded, hydrophobic interactions. Native structures for the five proteins are always found within the best 2% of all conformers generated. This novel method is simple and general and can be used to determine a small group of most favorable overall arrangements for the folding of specific amino acid sequences within a restricted space.     

Polypeptide secondary structure determination by nuclear magnetic resonance observation of short proton-proton distances

The use of proton-proton nuclear Overhauser enhancement (NOE) distance information for identification of polypeptide secondary structures in non-crystalline proteins was investigated by stereochemical studies of standard secondary structures and by statistical analyses of the secondary structures in the crystal conformations of a group of globular proteins. Both regular helix and beta-sheet secondary structures were found to contain a dense network of short 1H-1H distances. The results obtained imply that the combined information on all these distances obtained from visual inspection of the two-dimensional NOE (NOESY) spectra is sufficient for determination of the helical and beta-sheet secondary structures in small globular proteins. Furthermore, cis peptide bonds can be identified from unique, short sequential proton-proton distances. Limitations of this empirical approach are that the exact start or end of a helix may be difficult to define when the adjoining residues form a tight turn, and that unambiguous identification of tight turns can usually be obtained only in the hairpins of antiparallel beta-structures. The short distances between protons in pentapeptide segments of the different secondary structures have been tabulated to provide a generally applicable guide for the analysis of NOESY spectra of proteins.     

Evidence of Conformational Selection Driving the Formation of Ligand Binding Sites in Protein-Protein Interfaces

Many protein-protein interactions (PPIs) are compelling targets for drug discovery, and in a number of cases can be disrupted by small molecules. The main goal of this study is to examine the mechanism of binding site formation in the interface region of proteins that are PPI targets by comparing ligand-free and ligand-bound structures. To avoid any potential bias, we focus on ensembles of ligand-free protein conformations obtained by nuclear magnetic resonance (NMR) techniques and deposited in the Protein Data Bank, rather than on ensembles specifically generated for this study. The measures used for structure comparison are based on detecting binding hot spots, i.e., protein regions that are major contributors to the binding free energy. The main tool of the analysis is computational solvent mapping, which explores the surface of proteins by docking a large number of small “probe” molecules. Although we consider conformational ensembles obtained by NMR techniques, the analysis is independent of the method used for generating the structures. Finding the energetically most important regions, mapping can identify binding site residues using ligand-free models based on NMR data. In addition, the method selects conformations that are similar to some peptide-bound or ligand-bound structure in terms of the properties of the binding site. This agrees with the conformational selection model of molecular recognition, which assumes such pre-existing conformations. The analysis also shows the maximum level of similarity between unbound and bound states that is achieved without any influence from a ligand. Further shift toward the bound structure assumes protein-peptide or protein-ligand interactions, either selecting higher energy conformations that are not part of the NMR ensemble, or leading to induced fit. Thus, forming the sites in protein-protein interfaces that bind peptides and can be targeted by small ligands always includes conformational selection, although other recognition mechanisms may also be involved.