Deep-UV resonance Raman analysis of the Rhodobacter capsulatus cytochrome bc₁complex reveals a potential marker for the transmembrane peptide backbone

Classical strategies for structure analysis of proteins interacting with a lipid phase typically correlate ensemble secondary structure content measurements with changes in the spectroscopic responses of localized aromatic residues or reporter molecules to map regional solvent environments. Deep-UV resonance Raman (DUVRR) spectroscopy probes the vibrational modes of the peptide backbone itself, is very sensitive to the ensemble secondary structures of a protein, and has been shown to be sensitive to the extent of solvent interaction with the peptide backbone [ Wang , Y. , Purrello , R. , Georgiou , S. , and Spiro , T. G. ( 1991 ) J. Am. Chem. Soc. 113 , 6368 - 6377 ]. Here we show that a large detergent solubilized membrane protein, the Rhodobacter capsulatus cytochrome bc(1) complex, has a distinct DUVRR spectrum versus that of an aqueous soluble protein with similar overall secondary structure content. Cross-section calculations of the amide vibrational modes indicate that the peptide backbone carbonyl stretching modes differ dramatically between these two proteins. Deuterium exchange experiments probing solvent accessibility confirm that the contribution of the backbone vibrational mode differences are derived from the lipid solubilized or transmembrane α-helical portion of the protein complex. These findings indicate that DUVRR is sensitive to both the hydration status of a protein's peptide backbone, regardless of primary sequence, and its secondary structure content. Therefore, DUVRR may be capable of simultaneously measuring protein dynamics and relative water/lipid solvation of the protein.     

Evolutionary conservation of protein vibrational dynamics

The aim of the present work is to study the evolutionary divergence of vibrational protein dynamics. To this end, we used the Gaussian Network Model to perform a systematic analysis of normal mode conservation on a large dataset of proteins classified into homologous sets of family pairs and superfamily pairs. We found that the lowest most collective normal modes are the most conserved ones. More precisely, there is, on average, a linear correlation between normal mode conservation and mode collectivity. These results imply that the previously observed conservation of backbone flexibility (B-factor) profiles is due to the conservation of the most collective modes, which contribute the most to such profiles. We discuss the possible roles of normal mode robustness and natural selection in the determination of the observed behavior. Finally, we draw some practical implications for dynamics-based protein alignment and classification and discuss possible caveats of the present approach.     

Bending of the calmodulin central helix: a theoretical study.

The crystal structure of calcium-calmodulin (CaM) reveals a protein with a typical dumbbell structure. Various spectroscopic studies have suggested that the central linker region of CaM, which is alpha-helical in the crystal structure, is flexible in solution. In particular, NMR studies have indicated the presence of a flexible backbone between residues Lys 77 and Asp 80. This flexibility is related directly to the function of the protein because it enables the N- and C-terminal domains of the protein to move toward each other and bind to the CaM-binding domain of a target protein. We have investigated the flexibility of the CaM central helix by a variety of computational techniques: molecular dynamics (MD) simulations, normal mode analysis (NMA), and essential dynamics (ED) analysis. Our MD results reproduce the experimentally determined location of the bend in a simulation of only the CaM central helix, indicating that the bending point is an intrinsic property of the alpha-helix, for which the remainder of the protein is not important. Interestingly, the modes found by the ED analysis of the MD trajectory are very similar to the lowest frequency modes from the NM analysis and to modes found by an ED analysis of different structures in a set of NMR structures. Electrostatic interactions involving residues Arg 74 and Asp 80 seem to be important for these bending motions and unfolding, which is in line with pH-dependent NMR and CD studies.     

Normal mode analysis of macromolecular motions in a database framework: developing mode concentration as a useful classifying statistic

We investigated protein motions using normal modes within a database framework, determining on a large sample the degree to which normal modes anticipate the direction of the observed motion and were useful for motions classification. As a starting point for our analysis, we identified a large number of examples of protein flexibility from a comprehensive set of structural alignments of the proteins in the PDB. Each example consisted of a pair of proteins that were considerably different in structure given their sequence similarity. On each pair, we performed geometric comparisons and adiabatic-mapping interpolations in a high-throughput pipeline, arriving at a final list of 3,814 putative motions and standardized statistics for each. We then computed the normal modes of each motion in this list, determining the linear combination of modes that best approximated the direction of the observed motion. We integrated our new motions and normal mode calculations in the Macromolecular Motions Database, through a new ranking interface at http://molmovdb.org. Based on the normal mode calculations and the interpolations, we identified a new statistic, mode concentration, related to the mathematical concept of information content, which describes the degree to which the direction of the observed motion can be summarized by a few modes. Using this statistic, we were able to determine the fraction of the 3,814 motions where one could anticipate the direction of the actual motion from only a few modes. We also investigated mode concentration in comparison to related statistics on combinations of normal modes and correlated it with quantities characterizing protein flexibility (e.g., maximum backbone displacement or number of mobile atoms). Finally, we evaluated the ability of mode concentration to automatically classify motions into a variety of simple categories (e.g., whether or not they are "fragment-like"), in comparison to motion statistics. This involved the application of decision trees and feature selection (particular machine-learning techniques) to training and testing sets derived from merging the "list" of motions with manually classified ones.     

Application of the local regression method interval partial least-squares to the elucidation of protein secondary structure

The infrared amide bands are sensitive to the conformation of the polypeptide backbone of proteins. Since the backbone of proteins folds in complex spatial arrangements, the amide bands of these proteins result from the superimposition of vibration modes corresponding to the different types of structural motifs (alpha helices, beta sheets, etc.). Initially, band deconvolution techniques were applied to determine the secondary structure of proteins, i.e., the abundance of each structural motif in the polypeptide chain was directly related to the area of the suitable deconvolved vibration modes under the amide I band (1700-1600 cm(-1)). Recently, several multivariate regression methods have been used to predict the secondary structure of proteins as an alternative to the previous methods. They are based on establishing a relationship between a matrix of infrared protein spectra and another that includes their secondary structure, expressed as the fractions of the different structural motifs, determined from X-ray analysis. In this study, we investigated the use of the local regression method interval partial least-squares (iPLS) to seek improvements to the full-spectrum PLS and other regression methods. The local character of iPLS avoids the use of spectral regions that can introduce noise or that can be irrelevant for prediction and focuses on finding specific spectral ranges related to each secondary structure motif in all the proteins. This study has been applied to a representative protein data set with infrared spectra covering a large wavenumber range, including amides I-III bands (1700-1200 cm(-1)). iPLS has revealed new structural mode assignments related to less explored amide bands and has offered a satisfactory predictive ability using a small amount of selected specific spectral information.     

Models with energy penalty on interresidue rotation address insufficiencies of conventional elastic network models

In this study, I present a new elastic network model, to our knowledge, that addresses insufficiencies of two conventional models—the Gaussian network model (GNM) and the anisotropic network model (ANM). It has been shown previously that the GNM is not rotation-invariant due to its energy, which penalizes rigid-body rotation (external rotation). As a result, GNM models are found contaminated with rigid-body rotation, especially in the most collective ones. A new model (EPIRM) is proposed to remove such external component in modes. The extracted internal motions result from a potential that penalizes interresidue stretching and rotation in a protein. The new model is shown to pertinently describe crystallographic temperature factors (B-factors) and protein open↔closed transitions. Also, the capability of separating internal and external motions in GNM slow modes permits reexamining important mechanochemical properties in enzyme active sites. The results suggest that catalytic residues stay closer to rigid-body rotation axes than their immediate backbone neighbors. I show that the cumulative density of states for EPIRM and ANM follow different power laws as functions of low-mode frequencies. When using a cutoff distance of 7.5 Å, The cumulative density of states of EPIRM scales faster than that of all-atom normal mode analysis and slower than that of simple lattices.     

The role of heat-shock and chaperone proteins in protein folding: possible molecular mechanisms.

Recently some heat-shock proteins have been linked to functions of 'chaperoning' protein folding in vivo. Here current experimental evidence is reviewed and possible requirements for such an activity are discussed. It is proposed that one mode of chaperone action is to actively unfold misfolded or badly aggregated proteins to a conformation from which they could refold spontaneously; that improperly folded proteins are recognized by excessive stretches of solvent-exposed backbone, rather than by exposed hydrophobic patches; and that the molecular mechanism for unfolding is either repeated binding and dissociation ('plucking') or translocation of the protein backbone through a binding cleft ('threading'), allowing the threaded chain to refold spontaneously. The observed hydrolysis of ATP would provide the energy for active unfolding. These hypotheses can be applied to both monomeric folding and oligomeric assembly and are sufficiently detailed to be open to directed experimental verification.     

Modeling backbone flexibility to achieve sequence diversity: the design of novel alpha-helical ligands for Bcl-xL

Computational protein design can be used to select sequences that are compatible with a fixed-backbone template. This strategy has been used in numerous instances to engineer novel proteins. However, the fixed-backbone assumption severely restricts the sequence space that is accessible via design. For challenging problems, such as the design of functional proteins, this may not be acceptable. Here, we present a method for introducing backbone flexibility into protein design calculations and apply it to the design of diverse helical BH3 ligands that bind to the anti-apoptotic protein Bcl-xL, a member of the Bcl-2 protein family. We demonstrate how normal mode analysis can be used to sample different BH3 backbones, and show that this leads to a larger and more diverse set of low-energy solutions than can be achieved using a native high-resolution Bcl-xL complex crystal structure as a template. We tested several of the designed solutions experimentally and found that this approach worked well when normal mode calculations were used to deform a native BH3 helix structure, but less well when they were used to deform an idealized helix. A subsequent round of design and testing identified a likely source of the problem as inadequate sampling of the helix pitch. In all, we tested 17 designed BH3 peptide sequences, including several point mutants. Of these, eight bound well to Bcl-xL and four others showed weak but detectable binding. The successful designs showed a diversity of sequences that would have been difficult or impossible to achieve using only a fixed backbone. Thus, introducing backbone flexibility via normal mode analysis effectively broadened the set of sequences identified by computational design, and provided insight into positions important for binding Bcl-xL.     

Altering the RNA-binding mode of the U1A RBD1 protein

The N-terminal RNA-binding domain (RBD1) of the human U1A protein is evolutionarily designed to bind its RNA targets with great affinity and specificity. The physical mechanisms that modulate the coupling (local cooperativity) among amino acid residues on the extensive binding surface of RBD1 are investigated here, using mutants that replace a highly conserved glycine residue. This glycine residue, at the strand/loop junction of beta3/loop3, is found in U1A RBD1, and in most RBD domains, suggesting it has a specific role in modulation of RNA binding. Here, two RBD1 proteins are constructed in which that residue (Gly53) is replaced by either alanine or valine. These new proteins are shown by NMR methods and molecular dynamics simulations to be very similar to the wild-type RBD1, both in structure and in their backbone dynamics. However, RNA-binding assays show that affinity for the U1 snRNA stem-loop II RNA target is reduced by nearly 200-fold for the RBD1-G53A protein, and by 1.6 x 10(4)-fold for RBD1-G53V. The mode of RNA binding by RBD1-G53A is similar to that of RBD1-WT, displaying its characteristic non-additive free energies of base recognition and its salt-dependence. The binding mode of RBD1-G53V is altered, having lost its salt-dependence and displaying site-independence of base recognition. The molecular basis for this alteration in RNA-binding properties is proposed to result from the inability of the RNA to induce a change in the structure of the free protein to produce a high-affinity complex.     

Quantitative organization of the known protein x-ray structures. I. Methods and short-length-scale results

We address herein the problem of delineating the relationships between the known protein structures. In order to study this problem, methods have been developed to represent arbitrarily sized fragments of biopolymer backbone, and to compare distributions of such fragments. These methods are applied to a classification of 123 structures representing the entire set of known x-ray structures. The resulting data are analyzed (on the four-C alpha length scale) to determine both the large-scale organization of the set of known structures (i.e., the relationships between large groups of structures, each comprised of proteins that are structurally related) and its local structure (i.e., the quantitative degree of similarity between any two specific structures). It is shown that the set of structures forms a continuum of structural types, ranging from all-helical to all-sheet/barrel proteins. It is further demonstrated that the density of protein structures is not uniform across this continuum, but rather that structures cluster in certain regions, separated by regions of lower population. The properties of the various regions of the structural space are determined. The existence is demonstrated of strong quantitative correlations between the contents of different types of four-C alpha fragments within protein structures, which imply significant constraints on the types of architecture that can occur in proteins. Analysis of the distribution of structures demonstrates some hitherto unsuspected similarities and suggests that, in some circumstances, neither structural similarity nor sequence homology may be necessary conditions for evolutionary relationship between proteins. It is also suggested that these unsuspected similarities may imply similar folding mechanisms for structures of apparently different global architecture. Cases are also noted in which apparently similar structures may fold by different mechanisms. The connection between structure and dynamic properties is discussed, and a possible role of dynamics in the evolution of protein structures is suggested. The sensitivity of the methods presented herein to anomalies of structure refinement is demonstrated. It is suggested that the present results provide a framework for analyzing experimental results on structural similarity obtained using vibrational circular dichroism spectra, which are sensitive to local backbone structure.     

A method to search for similar protein local structures at ligand-binding sites and its application to adenine recognition

We have developed a method of searching for similar spatial arrangements of atoms around a given chemical moiety in proteins that bind a common ligand. The first step in this method is to consider a set of atoms that closely surround a given chemical moiety. Then, to compare the spatial arrangements of such surrounding atoms in different proteins, they are translated and rotated so that the chemical moieties are superposed on each other. Spatial arrangements of surrounding atoms in a pair of proteins are judged to be similar, when there are many corresponding atoms occupying similar spatial positions. Because the method focuses on the arrangements of surrounding atoms, it can detect structural similarities of binding sites in proteins that are dissimilar in their amino acid sequences or in their chain folds. We have applied this method to identify modes of nucleotide base recognition by proteins. An all-against-all comparison of the arrangements of atoms surrounding adenine moieties revealed an unexpected structural similarity between protein kinases, cAMP-dependent protein kinase (cAPK), and casein kinase-1 (CK1), and D-Ala:D-Ala ligase (DD-ligase) at their adenine-binding sites, despite a lack of similarity in their chain folds. The similar local structure consists of a four-residue segment and three sequentially separated residues. In particular the four-residue segments of these enzymes were found to have nearly identical conformations in their backbone parts, which are involved in the recognition of adenine. This common local structure was also found in substrate-free three-dimensional structures of other proteins that are similar to DD-ligase in the chain fold and of other protein kinases. As the proteins with different folds were found to share a common local structure, these proteins seem to constitute a remarkable example of convergent evolution for the same recognition mechanism. Received: 9 December 1996 / Accepted: 7 February 1997     

Osmolyte perturbation reveals conformational equilibria in spin‐labeled proteins

Recent evidence suggests that proteins at equilibrium can exist in a manifold of conformational substates, and that these substates play important roles in protein function. Therefore, there is great interest in identifying regions in proteins that are in conformational exchange. Electron paramagnetic resonance spectra of spin‐labeled proteins containing the nitroxide side chain (R1) often consist of two (or more) components that may arise from slow exchange between conformational substates (lifetimes > 100 ns). However, crystal structures of proteins containing R1 have shown that multicomponent spectra can also arise from equilibria between rotamers of the side chain itself. In this report, it is shown that these scenarios can be distinguished by the response of the system to solvent perturbation with stabilizing osmolytes such as sucrose. Thus, site‐directed spin labeling (SDSL) emerges as a new tool to explore slow conformational exchange in proteins of arbitrary size, including membrane proteins in a native‐like environment. Moreover, equilibrium between substates with even modest differences in conformation is revealed, and the simplicity of the method makes it suitable for facile screening of multiple proteins. Together with previously developed strategies for monitoring picosecond to millisecond backbone dynamics, the results presented here expand the timescale over which SDSL can be used to explore protein flexibility.     

Structure, dynamics, and insertion of a chloroplast targeting peptide in mixed micelles

Nuclear-encoded, chloroplast-destined proteins are synthesized with transit sequences that contain all information to get them inside the organelle. Different proteins are imported via a general protein import machinery, but their transit sequences do not share amino acid homology. It has been suggested that interactions between transit sequence and chloroplast envelope membrane lipids give rise to recognizable, structural motifs. In this study a detailed investigation of the structural, dynamical, and topological features of an isolated transit peptide associated with mixed micelles is described. The structure of the preferredoxin transit peptide in these micelles was studied by circular dichroism (CD) and multidimensional NMR techniques. CD experiments indicated that the peptide, which is unstructured in aqueous solution, obtained helical structure in the presence of the micelles. By NMR it is shown that the micelles introduced ill-defined helical structures in the transit peptide. Heteronuclear relaxation experiments showed that the whole peptide backbone is very flexible. The least dynamic segments are two N- and C-terminal helical regions flanking an unstructured proline-rich amino acid stretch. Finally, the insertion of the peptide backbone in the hydrophobic interior of the micelle was investigated by use of hydrophobic spin-labels. The combined data result in a model of the transit peptide structure, backbone dynamics, and insertion upon its interaction with mixed micelles.     

Excited state charge-transfer dynamics study of poplar plastocyanin by ultrafast pump-probe spectroscopy and molecular dynamics simulation

We have applied ultrafast pump-probe spectroscopy to investigate the excited state dynamics of the blue copper protein poplar plastocyanin, by exciting in the blue side of its 600-nm absorption band. The decay of the charge-transfer excited state occurs exponentially with a time constant of approximately 280 fs and is modulated by well visible oscillations. The Fourier transform of the oscillatory component, besides providing most of the vibrational modes found by conventional resonance Raman, presents additional bands in the low frequency region modes, which are reminiscent of collective motions of biological relevance. Notably, a high frequency mode at approximately 508 cm(-1), whose dynamics are consistent with that of the excited state and already observed for other blue copper proteins, is shown to be present also in poplar plastocyanin. This vibrational mode is reproduced by a molecular dynamics simulation involving the excited state of the copper site.     

Comparison of normal mode analyses on a small globular protein in dihedral angle space and Cartesian coordinate space

Normal mode analyses on the protein, bovine pancreatic trypsin inhibitor, in dihedral angle space and Cartesian coordinate space are compared. In Cartesian coordinate space it is found that modes of frequencies lower than 30 cm(-1) contribute 80% of the total mean-square fluctuation and are represented almost completely by motions in the dihedral angles. Bond angle and length fluctuations dominate in modes above 200 cm(-1), but contribute less than 2% to the total mean-square fluctuation. In the low-frequency modes a good correspondence between patterns of atomic displacements was found, but on average the root-mean-square fluctuations of the Cartesian coordinate modes are 13% greater than their dihedral angle counterparts. The main effect of fluctuations in the bond angles and lengths, therefore, is to allow the dihedral angles to become more flexible. As the important subspaces determined from the two methods overlap considerably, dihedral angle space analysis can be applied to proteins too large for Cartesian coordinate space analysis.     

Cooperative binding of polyamines induces the Escherichia coli single-strand binding protein-DNA binding mode transitions.

The Escherichia coli single-strand binding (SSB) protein is an essential protein involved in DNA replication, recombination, and repair processes. The tetrameric protein binds to ss nucleic acids in a number of different binding modes in vitro. These modes differ in the number of nucleotides occluded per SSB tetramer and in the type and degree of cooperative complexes that are formed with ss DNA. Although it is not yet known whether only one or all of these modes function in vivo, based on the dramatically different properties of the SSB tetramer in these different ss DNA binding modes, it has been suggested that the different modes may function selectively in replication, recombination, and/or repair. The transitions between these different modes are very sensitive to solution conditions, including salt (concentration, as well as cation and anion type), pH, and temperature. We have examined the effects of multivalent cations, principally the polyamine spermine, on the SSB-ss poly(dT) binding mode transitions and find that the transition from the (SSB)35 to the (SSB)56 binding mode can be induced by micromolar concentrations of polyamines as well as the inorganic cation Co(NH3)6(3+). Furthermore, these multivalent cations, as well as Mg2+, induce the binding mode transition by binding cooperatively to the SSB-poly(dT) complexes. These observations are interesting in light of the fact that polyamines, such as spermidine, are part of the ionic environment in E. coli and hence these cations are likely to affect the distribution of SSB-ss DNA binding modes in vivo.(ABSTRACT TRUNCATED AT 250 WORDS)     

Different mechanistic requirements for prokaryotic and eukaryotic chaperonins: a lattice study

MOTIVATION: The folding of many proteins in vivo and in vitro is assisted by molecular chaperones. A well-characterized molecular chaperone system is the chaperonin GroEL/GroES from Escherichia coli which has a homolog found in the eukaryotic cytosol called CCT. All chaperonins have a ring structure with a cavity in which the substrate protein folds. An interesting difference between prokaryotic and eukaryotic chaperonins is in the nature of the ATP-mediated conformational changes that their ring structures undergo during their reaction cycle. Prokaryotic chaperonins are known to exhibit a highly cooperative concerted change of their cavity surface while in eukaryotic chaperonins the change is sequential. Approximately 70% of proteins in eukaryotic cells are multi-domain whereas in prokaryotes single-domain proteins are more common. Thus, it was suggested that the different modes of action of prokaryotic and eukaryotic chaperonins can be explained by the need of eukaryotic chaperonins to facilitate folding of multi-domain proteins. RESULTS: Using a 2D square lattice model, we generated two large populations of single-domain and double-domain substrate proteins. Chaperonins were modeled as static structures with a cavity wall with which the substrate protein interacts. We simulated both concerted and sequential changes of the cavity surfaces and demonstrated that folding of single-domain proteins benefits from concerted but not sequential changes whereas double-domain proteins benefit also from sequential changes. Thus, our results support the suggestion that the different modes of allosteric switching of prokaryotic and eukaryotic chaperonin rings have functional implications as it enables eukaryotic chaperonins to better assist multi-domain protein folding.     

How does protein architecture facilitate the transduction of ATP chemical-bond energy into mechanical work? The cases of nitrogenase and ATP binding-cassette proteins

Transduction of adenosine triphosphate (ATP) chemical-bond energy into work to drive large-scale conformational changes is common in proteins. Two specific examples of ATP-utilizing proteins are the nitrogenase iron protein and the ATP binding-cassette transporter protein, BtuCD. Nitrogenase catalyzes biological nitrogen fixation whereas BtuCD transports vitamin B(12) across membranes. Both proteins drive their reactions with ATP. To interpret how the mechanical force generated by ATP binding and hydrolysis is propagated in these proteins, a coarse-grained elastic network model is employed. The analysis shows that subunits of the proteins move against each other in a concerted manner. The lowest-frequency modes of the nitrogenase iron protein and of the ATP binding-cassette transporter BtuCD protein are found to link the functionally critical domains, and these modes are suggested to be responsible for (at least the initial stages) large-scale ATP-coupled conformational changes.