Low frequency vibrational modes in proteins: Changes induced by point-mutations in the protein-cofactor matrix of bacterial reaction centers

As a step toward understanding their functional role, the low frequency vibrational motions (<300 cm−1) that are coupled to optical excitation of the primary donor bacteriochlorophyll cofactors in the reaction center from Rhodobacter sphaeroides were investigated. The pattern of hydrogen-bonding interaction between these bacteriochlorophylls and the surrounding protein was altered in several ways by mutation of single amino acids. The spectrum of low frequency vibrational modes identified by femtosecond coherence spectroscopy varied strongly between the different reaction center complexes, including between different mutants where the pattern of hydrogen bonds was the same. It is argued that these variations are primarily due to changes in the nature of the individual modes, rather than to changes in the charge distribution in the electronic states involved in the optical excitation. Pronounced effects of point mutations on the low frequency vibrational modes active in a protein-cofactor system have not been reported previously. The changes in frequency observed indicate a strong involvement of the protein in these nuclear motions and demonstrate that the protein matrix can increase or decrease the fluctuations of the cofactor along specific directions.

The use of site-directed fluorophore labeling and donor–donor energy migration to investigate solution structure and dynamics in proteins

The use of molecular genetics for introducing fluorescent molecules enables the use of donor–donor energy migration to determine intramolecular distances in a variety of proteins. This approach can be applied to examine the overall molecular dimensions of proteins and to investigate structural changes upon interactions with specific target molecules. In this report, the donor–donor energy migration method is demonstrated by experiments with the latent form of plasminogen activator inhibitor type 1. Based on the known x-ray structure of plasminogen activator inhibitor type 1, three positions forming the corners of a triangle were chosen. Double Cys substitution mutants (V106C-H185C, H185C-M266C, and M266C-V106C) and corresponding single substitution mutants (V106C, H185C, and M266C) were created and labeled with a sulfhydryl specific derivative of BODIPY (=the D molecule). The side lengths of this triangle were obtained from analyses of the experimental data. The analyses account for the local anisotropic order and rotational motions of the D molecules, as well as for the influence of a partial DD-labeling. The distances, as determined from x-ray diffraction, between the Cα-atoms of the positions V106C–H185C, H185C–M266C, and M266C–V106C were 60.9, 30.8, and 55.1 Å, respectively. These are in good agreement with the distances of 54 ± 4, 38 ± 3, and 55 ± 3 Å, as determined between the BODIPY groups attached via linkers to the same residues. Although the positions of the D-molecules and the Cα-atoms physically cannot coincide, there is a reasonable agreement between the methods.

Methionine residues as endogenous antioxidants in proteins

Cysteine and methionine are the two sulfur-containing residues normally found in proteins. Cysteine residues function in the catalytic cycle of many enzymes, and they can form disulfide bonds that contribute to protein structure. In contrast, the specific functions of methionine residues are not known. We propose that methionine residues constitute an important antioxidant defense mechanism. A variety of oxidants react readily with methionine to form methionine sulfoxide, and surface exposed methionine residues create an extremely high concentration of reactant, available as an efficient oxidant scavenger. Reduction back to methionine by methionine sulfoxide reductases would allow the antioxidant system to function catalytically. The effect of hydrogen peroxide exposure upon glutamine synthetase from Escherichia coli was studied as an in vitro model system. Eight of the 16 methionine residues could be oxidized with little effect on catalytic activity of the enzyme. The oxidizable methionine residues were found to be relatively surface exposed, whereas the intact residues were generally buried within the core of the protein. Furthermore, the susceptible residues were physically arranged in an array that guarded the entrance to the active site.

Selenium redox biochemistry of zinc–sulfur coordination sites in proteins and enzymes

Selenium has been increasingly recognized as an essential element in biology and medicine. Its biochemistry resembles that of sulfur, yet differs from it by virtue of both redox potentials and stabilities of its oxidation states. Selenium can substitute for the more ubiquitous sulfur of cysteine and as such plays an important role in more than a dozen selenoproteins. We have chosen to examine zinc–sulfur centers as possible targets of selenium redox biochemistry. Selenium compounds release zinc from zinc/thiolate-coordination environments, thereby affecting the cellular thiol redox state and the distribution of zinc and likely of other metal ions. Aromatic selenium compounds are excellent spectroscopic probes of the otherwise relatively unstable functional selenium groups. Zinc-coordinated thiolates, e.g., metallothionein (MT), and uncoordinated thiolates, e.g., glutathione, react with benzeneseleninic acid (oxidation state +2), benzeneselenenyl chloride (oxidation state 0) and selenocystamine (oxidation state −1). Benzeneseleninic acid and benzeneselenenyl chloride react very rapidly with MT and titrate substoichiometrically and with a 1:1 stoichiometry, respectively. Selenium compounds also catalyze the release of zinc from MT in peroxidation and thiol/disulfide-interchange reactions. The selenoenzyme glutathione peroxidase catalytically oxidizes MT and releases zinc in the presence of t-butyl hydroperoxide, suggesting that this type of redox chemistry may be employed in biology for the control of metal metabolism. Moreover, selenium compounds are likely targets for zinc/thiolate coordination centers in vivo, because the reactions are only partially suppressed by excess glutathione. This specificity and the potential to undergo catalytic reactions at low concentrations suggests that zinc release is a significant aspect of the therapeutic antioxidant actions of selenium compounds in antiinflammatory and anticarcinogenic agents.

Dynamical principles in biological processes: A model of charge migration in proteins and DNA

The generalized master equations (GMEs) that contain multiple time scales have been derived quantum mechanically. The GME method has then been applied to a model of charge migration in proteins that invokes the hole hopping between local amino acid sites driven by the torsional motions of the floppy backbones. This model is then applied to analyze the experimental results for sequence-dependent long-range hole transport in DNA reported by Meggers et al. [Meggers, E., Michel-Beyerle, M. E., & Giese, B. (1998) J. Am. Chem. Soc. 120, 12950–12955]. The model has also been applied to analyze the experimental results of femtosecond dynamics of DNA-mediated electron transfer reported by Zewail and co-workers [Wan, C., Fiebig, T., Kelley, S. O., Treadway, C. R., Barton, J. K. & Zewail, A. H. (1999) Proc. Natl. Acad. Sci. USA 96, 6014–6019]. The initial events in the dynamics of protein folding have begun to attract attention. The GME obtained in this paper will be applicable to this problem.

Structural determinants of oxidative folding in proteins

A method for determining the kinetic fate of structured disulfide species (i.e., whether they are preferentially oxidized or reshuffle back to an unstructured disulfide species) is introduced. The method relies on the sensitivity of unstructured disulfide species to low concentrations of reducing agents. Because a structured des species that preferentially reshuffles generally first rearranges to an unstructured species, a small concentration of reduced DTT (e.g., 260 μM) suffices to distinguish on-pathway intermediates from dead-end species. We apply this method to the oxidative folding of bovine pancreatic ribonuclease A (RNase A) and show that des[40–95] and des[65–72] are productive intermediates, whereas des[26–84] and des[58–110] are metastable dead-end species that preferentially reshuffle. The key factor in determining the kinetic fate of these des species is the relative accessibility of both their thiol groups and disulfide bonds. Productive intermediates tend to be disulfide-secure, meaning that their structural fluctuations preferentially expose their thiol groups, while keeping their disulfide bonds buried. By contrast, dead-end species tend to be disulfide-insecure, in that their structural fluctuations expose their disulfide bonds in concert with their thiol groups. This distinction leads to four generic types of oxidative folding pathways. We combine these results with those of earlier studies to suggest a general three-stage model of oxidative folding of RNase A and other single-domain proteins with multiple disulfide bonds.

Point mutations and sequence variability in proteins: Redistributions of preexisting populations

Here we study the effect of point mutations in proteins on the redistributions of the conformational substates. We show that regardless of the location of a mutation in the protein structure and of its type, the observed movements of the backbone recur largely at the same positions in the structures. Despite the different interactions that are disrupted and formed by the residue substitution, not only are the conformations very similar, but the regions that move are also the same, regardless of their sequential or spatial distance from the mutation. This observation leads us to conclude that, apart from some extreme cases, the details of the interactions are not critically important in determining the protein conformation or in specifying which parts of the protein would be more prone to take on different local conformations in response to changes in the sequence. This finding further illustrates why proteins manifest a robustness toward many mutational events. This nonuniform distribution of the conformer population is consistently observed in a variety of protein structural types. Topology is critically important in determining folding pathways, kinetics, building block cutting, and anatomy trees. Here we show that topology is also very important in determining which regions of the protein structure will respond to sequence changes, regardless of the sequential or spatial location of the mutation.

BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} acts as a phosphate analog in proteins phosphorylated on aspartate: Structure of a BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} complex with phosphoserine phosphatase

Protein phosphoaspartate bonds play a variety of roles. In response regulator proteins of two-component signal transduction systems, phosphorylation of an aspartate residue is coupled to a change from an inactive to an active conformation. In phosphatases and mutases of the haloacid dehalogenase (HAD) superfamily, phosphoaspartate serves as an intermediate in phosphotransfer reactions, and in P-type ATPases, also members of the HAD family, it serves in the conversion of chemical energy to ion gradients. In each case, lability of the phosphoaspartate linkage has hampered a detailed study of the phosphorylated form. For response regulators, this difficulty was recently overcome with a phosphate analog, BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}, which yields persistent complexes with the active site aspartate of their receiver domains. We now extend the application of this analog to a HAD superfamily member by solving at 1.5-Å resolution the x-ray crystal structure of the complex of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} with phosphoserine phosphatase (PSP) from Methanococcus jannaschii. The structure is comparable to that of a phosphoenzyme intermediate: BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} is bound to Asp-11 with the tetrahedral geometry of a phosphoryl group, is coordinated to Mg2+, and is bound to residues surrounding the active site that are conserved in the HAD superfamily. Comparison of the active sites of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}⋅PSP and BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}⋅CeY, a receiver domain/response regulator, reveals striking similarities that provide insights into the function not only of PSP but also of P-type ATPases. Our results indicate that use of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} for structural studies of proteins that form phosphoaspartate linkages will extend well beyond response regulators.

Global properties of the mapping between local amino acid sequence and local structure in proteins.

Local protein structure prediction efforts have consistently failed to exceed approximately 70% accuracy. We characterize the degeneracy of the mapping from local sequence to local structure responsible for this failure by investigating the extent to which similar sequence segments found in different proteins adopt similar three-dimensional structures. Sequence segments 3-15 residues in length from 154 different protein families are partitioned into neighborhoods containing segments with similar sequences using cluster analysis. The consistency of the sequence-to-structure mapping is assessed by comparing the local structures adopted by sequence segments in the same neighborhood in proteins of known structure. In the 154 families, 45% and 28% of the positions occur in neighborhoods in which one and two local structures predominate, respectively. The sequence patterns that characterize the neighborhoods in the first class probably include virtually all of the short sequence motifs in proteins that consistently occur in a particular local structure. These patterns, many of which occur in transitions between secondary structural elements, are an interesting combination of previously studied and novel motifs. The identification of sequence patterns that consistently occur in one or a small number of local structures in proteins should contribute to the prediction of protein structure from sequence.

Cloning the expression of a mammalian gene involved in the reduction of methionine sulfoxide residues in proteins.

An enzyme that reduces methionine sulfoxide [Met(O)] residues in proteins [peptide Met(O) reductase (MsrA), EC 1.8.4.6; originally identified in Escherichia coli] was purified from bovine liver, and the cDNA encoding this enzyme was cloned and sequenced. The mammalian homologue of E. coli msrA (also called pmsR) cDNA encodes a protein of 255 amino acids with a calculated molecular mass of 25,846 Da. This protein has 61% identity with the E. coli MsrA throughout a region encompassing a 199-amino acid overlap. The protein has been overexpressed in E. coli and purified to homogeneity. The mammalian recombinant MsrA can use as substrate, proteins containing Met(O) as well as other organic compounds that contain an alkyl sulfoxide group such as N-acetylMet(O), Met(O), and dimethyl sulfoxide. Northern analysis of rat tissue extracts showed that rat msrA mRNA is present in a variety of organs with the highest level found in kidney. This is consistent with the observation that kidney extracts also contained the highest level of enzyme activity.

A method for distance determination in proteins using a designed metal ion binding site and site-directed spin labeling: evaluation with T4 lysozyme.

The use of molecular genetics to introduce both a metal ion binding site and a nitroxide spin label into the same protein opens the use of paramagnetic metalnitroxyl interactions to estimate intramolecular distances in a wide variety of proteins. In this report, a His-Xaa3-His metal ion binding motif was introduced at the N terminus of the long interdomain helix of T4 lysozyme (Lys-65 --> His/Gln-69 --> His) of three mutants, each containing a single nitroxide-labeled cysteine residue at position 71, 76, or 80. The results show that Cu(II)-induced relaxation effects on the nitroxide can be quantitatively analyzed in terms of interspin distance in the range of 10-25 A using Redfield theory, as first suggested by Leigh [Leigh, J.S. (1970) J. Chem. Phys. 52, 2608-2612]. Of particular interest is the observation that distances can be determined both under rigid lattice conditions in frozen solution and in the presence of motion of the spins at room temperature under physiological conditions. The method should be particularly attractive for investigating structure in membrane proteins that are difficult to crystallize. In the accompanying paper, the technique is applied to a polytopic membrane protein, lactose permease.

Distance determination in proteins using designed metal ion binding sites and site-directed spin labeling: application to the lactose permease of Escherichia coli.

As shown in the accompanying paper, the magnetic dipolar interaction between site-directed metal-nitroxide pairs can be exploited to measure distances in T4 lysozyme, a protein of known structure. To evaluate this potentially powerful method for general use, particularly with membrane proteins that are difficult to crystallize, both a paramagnetic metal ion binding site and a nitroxide side chain were introduced at selected positions in the lactose permease of Escherichia coli, a paradigm for polytopic membrane proteins. Thus, three individual cysteine residues were introduced into putative helix IV of a lactose permease mutant devoid of native cysteine residues containing a high-affinity divalent metal ion binding site in the form of six contiguous histidine residues in the periplasmic loop between helices III and IV. In addition, the construct contained a biotin acceptor domain in the middle cytoplasmic loop to facilitate purification. After purification and spin labeling, electron paramagnetic resonance spectra were obtained with the purified proteins in the absence and presence of Cu(II). The results demonstrate that positions 103, 111, and 121 are 8, 14, and > 23 A from the metal binding site. These data are consistent with an alpha-helical conformation of transmembrane domain IV of the permease. Application of the technique to determine helix packing in lactose permease is discussed.

On the existence and implications of an inverse folding code in proteins.

The existence of a code relating the set of possible sequences at a given position in a protein backbone to the local structure at that location is investigated. It is shown that only 73% of 4-C alpha structure fragments in a sample of 114 protein structures exhibit a preference for a particular set of sequences. The remaining structures can accommodate essentially any sequence. The structures that encode specific sequence distributions include the classical "secondary" structures, with the notable exception of planar (beta) bends. It is suggested that this has implications as to the mechanism of folding in proteins with extensive sheet/barrel structure. The possible role of structures that do not encode specific sequences as mutation hot spots is noted.

Efficient folding of proteins with multiple disulfide bonds in the Escherichia coli cytoplasm

Under physiological conditions, the Escherichia coli cytoplasm is maintained in a reduced state that strongly disfavors the formation of stable disulfide bonds in proteins. However, mutants in which the reduction of both thioredoxins and glutathione is impaired (trxB gor mutants) accumulate oxidized, enzymatically active alkaline phosphatase in the cytoplasm. These mutants grow very poorly in the absence of an exogenous reductant and accumulate extragenic suppressors at a high frequency. One such suppressor strain, FA113, grows almost as rapidly as the wild type in the absence of reductant, exhibits slightly faster kinetics of disulfide bond formation, and has fully induced activity of the transcriptional activator, OxyR. FA113 gave substantially higher yields of properly oxidized proteins compared with wild-type or trxB mutant strains. For polypeptides with very complex patterns of disulfide bonds, such as vtPA and the full-length tPA, the amount of active protein was further enhanced up to 15-fold by co-expression of TrxA (thioredoxin 1) mutants with different redox potentials, or 20-fold by the protein disulfide isomerase, DsbC. Remarkably, higher yields of oxidized, biologically active proteins were obtained by expression in the cytoplasm of E. coli FA113 compared with what could be achieved via secretion into the periplasm of a wild-type strain, even under optimized conditions. These results demonstrate that the cytoplasm can be rendered sufficiently oxidizing to allow efficient formation of native disulfide bonds without compromising cell viability.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.