Based on results from both equilibrium and kinetic hydrogen exchange studies of Escherichia coli ribonuclease HI (RNase H), a fragment of RNase H (eABCD) was designed. The sequence of eABCD contains less than half of the protein's primary sequence and includes the regions that were shown to be the most protected from hydrogen exchange in all previous studies of RNase H. This core fragment of RNase H encodes a well-ordered protein with native-like properties. When isolated from the full-length monomeric protein, the eABCD fragment forms a stable dimer. However, we show indirectly that the monomeric form of eABCD is folded and has an overall secondary structure similar to the dimeric form.

Desulforedoxin (Dx), isolated from the sulfate reducing bacterium Desulfovibrio gigas, is a small homodimeric (2 × 36 amino acids) protein. Each subunit contains a high-spin iron atom tetrahedrally bound to four cysteinyl sulfur atoms, a metal center similar to that found in rubredoxin (Rd) type proteins. The simplicity of the active center in Dx and the possibility of replacing the iron by other metals make this protein an attractive case for the crystallographic analysis of metal-substituted derivatives. This study extends the relevance of Dx to the bioinorganic chemistry field and is important to obtain model compounds that can mimic the four sulfur coordination of metals in biology. Metal replacement experiments were carried out by reconstituting the apoprotein with In3+, Ga3+, Cd2+, Hg2+, and Ni2+ salts. The In3+ and Ga3+ derivatives are isomorphous with the iron native protein; whereas Cd2+, Hg2+, and Ni2+ substituted Dx crystallized under different experimental conditions, yielding two additional crystal morphologies; their structures were determined by the molecular replacement method. A comparison of the three-dimensional structures for all metal derivatives shows that the overall secondary and tertiary structures are maintained, while some differences in metal coordination geometry occur, namely, bond lengths and angles of the metal with the sulfur ligands. These data are discussed in terms of the entatic state theory.

We describe an extensive test of Geocore, an ab initio peptide folding algorithm. We studied 18 short molecules for which there are structures in the Protein Data Bank; chains are up to 31 monomers long. Except for the very shortest peptides, an extremely simple energy function is sufficient to discriminate the true native state from more than 108 lowest energy conformations that are searched explicitly for each peptide. A high incidence of native-like structures is found within the best few hundred conformations generated by Geocore for each amino acid sequence. Predictions improve when the number of discrete φ/ψ choices is increased.

An analysis of the pairwise side-chain packing geometries of cysteine residues observed in high-resolution protein crystal structures indicates that cysteine pairs have pronounced orientational preferences due to the geometric constraints of disulfide bond formation. A potential function was generated from these observations and used to evaluate models for novel disulfide bonds in human carbonic anhydrase II (HCAII). Three double-cysteine variants of HCAII were purified and the effective concentrations of their thiol groups were determined by titrations with glutathione and dithiothreitol. The effects of the cysteine mutations on the native state structure and stability were characterized by circular dichroism, enzymatic activity, sulfonamide binding, and guanidine hydrochloride titration. These analyses indicate that the PAIRWISE potential is a good predictor of the strength of the disulfide bond itself, but the overall structural and thermodynamic effects on the protein are complicated by additional factors. In particular, the effects of cysteine substitutions on the native state and the stabilization of compact nonnative states by the disulfide can override any stabilizing effect of the cross-link.

The standard molecular model of the fiber of the sickle hemoglobin (HbS: β6 Glu → Val) has been revised to allow both β6 mutation sites to participate in intermolecular contacts, rather than only one β6 site as previously thought, for four molecules per 14-molecule fiber cross section. This structure accurately predicts the copolymerization of hybridized mixtures of HbS with HbA or HbC (β6 Glu → Lys), which could not be reconciled with prior models in which only half the β6 sites were required for assembly. This model suggests new contacts within the fiber and raises the question of whether these cross-linked double strands could possess added stability important in such processes as nucleation.

Nanoflow electrospray mass spectrometry was used to monitor the formation of protein heterodimers of HU proteins from Bacillus stearothermophilus and Bacillus subtilis. This has enabled us to analyze both thermodynamic and kinetic features associated with the dissociation of homodimeric HU proteins. The results obtained correlate well with the kinetics of the protein dissociation process and the free energy difference between homo- and heterodimeric species anticipated from other studies. We suggest that this approach will have general applicability in studying protein association and dissociation under near-equilibrium conditions and will be relevant to a wide range of biological systems.

Effects of guanidine hydrochloride (GdnHCl) on the structure and dynamics of wild-type Humicola lanuginosa lipase (HLL) and its two mutants were studied. The latter were S146A (with the active site Ser replaced by Ala) and the single Trp mutant W89m, with substitutions W117F, W221H, and W260H. Steady-state, stopped-flow, and time-resolved laser-induced fluorescence spectroscopy were carried out as a function of [GdnHCl]. The maximum emission wavelength and fluorescence lifetimes revealed the microenvironment of the tryptophan(s) in these lipases to become more polar upon increasing [GdnHCl]. However, significant extent of tertiary structure in GdnHCl is suggested by the observation that both wild-type HLL and W89m remain catalytically active at rather high GdnHCl concentrations of >6 and 4.0 M, respectively. Changes in steady-state emission anisotropy, as well as variation in rotational correlation times and residual anisotropy values, demonstrate that upon increasing [GdnHCl] the structure of the lipases became more loose, with increasing amplitude of structural fluctuations. Finally, intermediate states in the course of exposure of the proteins to GdnHCl were revealed by stopped-flow fluorescence measurements.

We describe here an algorithm for distinguishing sequential from nonsequentially folding proteins. Several experiments have recently suggested that most of the proteins that are synthesized in the eukaryotic cell may fold sequentially. This proposed folding mechanism in vivo is particularly advantageous to the organism. In the absence of chaperones, the probability that a sequentially folding protein will misfold is reduced significantly. The problem we address here is devising a procedure that would differentiate between the two types of folding patterns. Footprints of sequential folding may be found in structures where consecutive fragments of the chain interact with each other. In such cases, the folding complexity may be viewed as being lower. On the other hand, higher folding complexity suggests that at least a portion of the polypeptide backbone folds back upon itself to form three-dimensional (3D) interactions with noncontiguous portion(s) of the chain. Hence, we look at the mechanism of folding of the molecule via analysis of its complexity, that is, through the 3D interactions formed by contiguous segments on the polypeptide chain. To computationally splice the structure into consecutively interacting fragments, we either cut it into compact hydrophobic folding units or into a set of hypothetical, transient, highly populated, contiguous fragments (“building blocks” of the structure). In sequential folding, successive building blocks interact with each other from the amino to the carboxy terminus of the polypeptide chain. Consequently, the results of the parsing differentiate between sequentially vs. nonsequentially folded chains. The automated assessment of the folding complexity provides insight into both the likelihood of misfolding and the kinetic folding rate of the given protein. In terms of the funnel free energy landscape theory, a protein that truly follows the mechanism of sequential folding, in principle, encounters smoother free energy barriers. A simple sequentially folded protein should, therefore, be less error prone and fold faster than a protein with a complex folding pattern.

Six helix surface positions of protein G (Gβ1) were redesigned using a computational protein design algorithm, resulting in the five fold mutant Gβ1m2. Gβ1m2 is well folded with a circular dichroism spectrum nearly identical to that of Gβ1, and a melting temperature of 91 °C, ∼6 °C higher than that of Gβ1. The crystal structure of Gβ1m2 was solved to 2.0 Å resolution by molecular replacement. The absence of hydrogen bond or salt bridge interactions between the designed residues in Gβ1m2 suggests that the increased stability of Gβ1m2 is due to increased helix propensity and more favorable helix dipole interactions.

Conserved deamidation of PKA catalytic subunit isozymes Cα and Cβ—more than 25% at Asn2 in vivo in both cases—has been shown to yield Asp2- and isoAsp2-containing isozymes (Jedrzejewski PT, Girod A, Tholey A, König N, Thullner S, Kinzel V, Bossemeyer D, 1998, Protein Sci 7:457–469). Isoaspartate formation in proteins in vivo is indicative of succinimide intermediates involved in both the initial deamidation reaction as well as the “repair” of isoAsp to Asp by the action of protein l-isoaspartyl (d-aspartyl) O-methyl transferase (PIMT). l-Succinimide is prone to racemization to d-succinimide, which may hydrolyze to d-isoAsp- and d-Asp-containing diastereomers with, respectively, no and poor substrate character for PIMT. To analyze native PKA catalytic subunit from cardiac muscle for these isomers the N-terminal tryptic peptides (T1) of the enzyme were analyzed following procedures refined specifically with a set of corresponding synthetic peptides. The methods combined high resolution high-performance liquid chromatography and a new mass spectrometric procedure for the discrimination between Asp- and isoAsp-residues in peptides (Lehmann et al., 2000). The results demonstrate the occurrence of d-isoAsp- and d-Asp-containing T1 fragments in addition to the l-isomers. The small amount of the l-isoAsp isomer, representing only part of the d-isoAsp isomer, and the relatively large amounts of the l-Asp and d-Asp isomers argues for an effective action of PIMT present in cardiac tissue.

A right-handed parallel β-helix of 400 residues in 13 tightly packed coils is a major motif of the chains forming the trimeric P22 tailspike adhesin. The β-helix domains of three identical subunits are side-by-side in the trimer and make predominantly hydrophilic inter-subunit contacts (Steinbacher S et al., 1994, Science 265:383–386). After the 13th coil the three individual β-helices terminate and the chains wrap around each other to form three interdigitated β-sheets organized into the walls of a triangular prism. The β-strands then separate and form antiparallel β-sheets, but still defining a triangular prism in which each side is a β-sheet from a different subunit (Seckler R, 1998, J Struct Biol 122:216–222). The subunit interfaces are buried in the triangular core of the prism, which is densely packed with hydrophobic side chains from the three β-sheets. Examination of this structure reveals that its packed core maintains the same pattern of interior packing found in the left-handed β-helix, a single-chain structure. This packing is maintained in both the interdigitated parallel region of the prism and the following antiparallel sheet section. This oligomerization motif for the tailspike β-helices presumably contributes to the very high thermal and detergent stability that is a property of the native tailspike adhesin.

The effects of a number of cryoprotectants on the kinetic and structural properties of glycogen phosphorylase b have been investigated. Kinetic studies showed that glycerol, one of the most commonly used cryoprotectants in X-ray crystallographic studies, is a competitive inhibitor with respect to substrate glucose-1-P with an apparent Ki value of 3.8% (v/v). Cryogenic experiments, with the enzyme, have shown that glycerol binds at the catalytic site and competes with glucose analogues that bind at the catalytic site, thus preventing the formation of complexes. This necessitated a change in the conditions for cryoprotection in crystallographic binding experiments with glycogen phosphorylase. It was found that 2-methyl-2,4-pentanediol (MPD), polyethylene glycols (PEGs) of various molecular weights, and dimethyl sulfoxide (DMSO) activated glycogen phosphorylase b to different extents, by stabilizing its most active conformation, while sucrose acted as a noncompetitive inhibitor and elethylene glycol as an uncompetitve inhibitor with respect to glucose-1-P. A parallel experimental investigation by X-ray crystallography showed that, at 100 K, both MPD and DMSO do not bind at the catalytic site, do not induce any significant conformational change on the enzyme molecule, and hence, are more suitable cryoprotectants than glycerol for binding studies with glycogen phosphorylase.

Pure rat brain tubulin is readily palmitoylated in vitro using [3H]palmitoyl CoA but no added enzymes. A maximum of approximately six palmitic acids are added per dimer in 2–3 h at 36–37 °C under native conditions. Both α and β tubulin are labeled, and 63–73% of the label was hydroxylamine-labile, presumed thioesters. Labeling increases with increasing pH and temperature, and with low concentrations of guanidine HCl or KCl (but not with urea) to a maximum of ∼13 palmitates/dimer. High SDS and guanidine HCl concentrations are inhibitory. At no time could all 20 cysteine residues of the dimer be palmitoylated. Polymerization to microtubules, or use of tubulin S, markedly decreases the accessibility of the palmitoylation sites. Palmitoylation increases the electrophoretic mobility of a portion of α tubulin toward the β band. Palmitoylated tubulin binds a colchicine analogue normally, but during three warm/cold polymerization/depolymerization cycles there is a progressive loss of palmitoylated tubulin, indicating decreased polymerization competence. We postulate that local electrostatic factors are major regulators of reactivity of tubulin cysteine residues toward palmitoyl CoA, and that the negative charges surrounding a number of the cysteines are sensitive to negative charges on palmitoyl CoA.

High mobility group (HMG) proteins 14 and 17 are nonhistone nuclear proteins that have been implicated in control of transcription and chromatin structure. To examine the posttranslational modifications of HMG-14 and -17 in vivo, HMG proteins were prepared from nuclear vs. cytosolic fractions of human K562 cells treated with 12-O-tetradecanoylphorbol 13-acetate (TPA) or okadaic acid (OA) and examined by electrospray mass spectrometry. Analysis of full-length masses demonstrated mono-, di-, and triphosphorylation of HMG-14 and mono- and diphosphorylation of HMG-17 from OA treated cells, whereas HMG-14 and -17 from TPA treated cells were monophosphorylated. Peptide mass and sequence analysis showed major and minor phosphorylation sites, respectively, at Ser24 and Ser28 in HMG-17, and Ser20 and Ser24 in HMG-14. These sites were found in the consensus sequence RRSARLSAK, within the nucleosomal binding domain of each protein. A third phosphorylation site in HMG-14 was located at either Ser6 or Ser7. Interestingly, the proportion of HMG-14 and -17 found in cytosolic pools increased significantly after 1 h of treatment compared to control cells and showed preferential phosphorylation compared with proteins from nuclear fractions. These results suggest that phosphorylation of HMG-14 and -7 interferes with nuclear localization mechanisms in a manner favoring release from nuclei.

Parvalbumins constitute a class of calcium-binding proteins characterized by the presence of several helix-loop-helix (EF-hand) motifs. In a previous study (Revett SP, King G, Shabanowitz J, Hunt DF, Hartman KL, Laue TM, Nelson DJ, 1997, Protein Sci 7:2397–2408), we presented the sequence of the major parvalbumin isoform from the silver hake (Merluccius bilinearis) and presented spectroscopic and structural information on the excised “EF-hand” portion of the protein. In this study, the X-ray crystal structure of the silver hake major parvalbumin has been determined to high resolution, in the frozen state, using the molecular replacement method with the carp parvalbumin structure as a starting model. The crystals are orthorhombic, space group C2221, with a = 75.7 Å, b = 80.7 Å, and c = 42.1 Å. Data were collected from a single crystal grown in 15% glycerol, which served as a cryoprotectant for flash freezing at −188 °C. The structure refined to a conventional R-value of 21% (free R 25%) for observed reflections in the range 8 to 1.65 Å [I > 2σ(I)]. The refined model includes an acetylated amino terminus, 108 residues (characteristic of a β parvalbumin lineage), 2 calcium ions, and 114 water molecules per protein molecule. The resulting structure was used in molecular dynamics (MD) simulations focused primarily on the dynamics of the ligands coordinating the Ca2+ ions in the CD and EF sites. MD simulations were performed on both the fully Ca2+ loaded protein and on a Ca2+ deficient variant, with Ca2+ only in the CD site. There was substantial agreement between the MD and X-ray results in addressing the issue of mobility of key residues in the calcium-binding sites, especially with regard to the side chain of Ser55 in the CD site and Asp92 in the EF site.

The ETS family members display specific DNA binding site preferences. As an example, PU.1 and ETS-1 recognize different DNA sequences with a core element centered over 5′-GGAA-3′ and 5′-GGAA/T-3′, respectively. To understand the molecular basis of this recognition, we carried out site-directed mutagenesis experiments followed by DNA binding studies that use electrophoretic mobility shift assay (EMSA) and surface plasmon resonance methods. EMSA experiments identified amino acid changes A231S and/or N236Y as being important for PU.1 recognition of both 5′-GGAA-3′ and 5′-GGAT-3′ containing oligonucleotides. To confirm these data and obtain accurate binding parameters, we performed kinetic studies using surface plasmon resonance on these mutants. The N236Y substitution revealed a weak protein-DNA interaction with the 5′-GGAA-3′ containing oligonucleotide caused by a faster release of the protein from the DNA (koff tenfold higher than the wild-type protein). With the double mutant A231S-N236Y, we obtained an increase in binding affinity and stability toward both 5′-GGAA-3′ and 5′-GGAT-3′ containing oligonucleotides. We propose that substitution of alanine for serine introduces an oxygen atom that can accept hydrogen and interact with potential water molecules or other atoms to make an energetically favorable hydrogen bond with both 5′-GGAA-3′ and 5′-GGAT-3′ oligonucleotides. The free energy of dissociation for the double mutant A231S-N236Y with 5′-GGAA-3′ (ΔΔG((A231S-N236Y) − (N236Y)) = −1.2 kcal mol−1) confirm the stabilizing effect of this mutant in the protein-DNA complex formation. We conclude that N236Y mutation relaxes the specificity toward 5′-GGAA-3′ and 5′-GGAT-3′ sequences, while A231S mutation modulates the degree of specificity toward 5′-GGAA-3′ and 5′GGAT-3′ sequences. This study explains why wild-type PU.1 does not recognize 5′-GGAT-3′ sequences and in addition broadens our understanding of 5′-GGAA/T-3′ recognition by ETS protein family members.

The consequences of amino acid substitutions at the dimer interface for the strength of the interactions between the monomers and for the catalytic function of the dimeric enzyme alkaline phosphatase from Escherichia coli have been investigated. The altered enzymes R10A, R10K, R24A, R24K, T59A, and R10A/R24A, which have amino acid substitutions at the dimer interface, were characterized using kinetic assays, ultracentrifugation, and transverse urea gradient gel electrophoresis. The kinetic data for the wild-type and altered alkaline phosphatases show comparable catalytic behavior with kcat values between 51.3 and 69.5 s−1 and Km values between 14.8 and 26.3 μM. The ultracentrifugation profiles indicate that the wild-type enzyme is more stable than all the interface-modified enzymes. The wild-type enzyme is dimeric in the pH range of pH 4.0 and above, and disassembled at pH 3.5 and below. All the interface-modified enzymes, however, are apparently monomeric at pH 4.0, begin assembly at pH 5.0, and are not fully assembled into the dimeric form until pH 6.0. The results from transverse urea gradient gel electrophoresis show clear and reproducible differences both in the position and the shape of the unfolding patterns; all these modified enzymes are more sensitive to the denaturant and begin to unfold at urea concentrations between 1.0 and 1.5 M; the wild-type enzyme remains in the folded high mobility form beyond 2.5 M urea. Alkaline phosphatase H370A, modified at the active site and not at the dimer interface, resembles the wild-type enzyme both in ultracentrifugation and electrophoresis studies. The results obtained suggest that substitution of a single amino acid at the interface sacrifices not only the integrity of the assembled dimer, but also the stability of the monomer fold, even though the activity of the enzyme at optimal pH remains unaffected and does not appear to depend on interface stability.

Pokeweed antiviral protein (PAP) is a ribosome-inactivating protein (RIP), which enzymatically removes a single adenine base from a conserved, surface exposed loop sequence of ribosomal rRNA. We now present unprecedented experimental evidence that PAP can release not only adenine but guanine as well from Escherichia coli rRNA, albeit at a rate 20 times slower than for adenine. We also report X-ray structure analysis and supporting modeling studies for the interactions of PAP with guanine. Our modeling studies indicated that PAP can accommodate a guanine base in the active site pocket without large conformational changes. This prediction was experimentally confirmed, since a guanine base was visible in the active site pocket of the crystal structure of the PAP-guanine complex.

Thymidylate synthase (TS) converts dUMP to dTMP by reductive methylation, where 5,10-methylenetetrahydrofolate is the source of both the methylene group and reducing equivalents. The mechanism of this reaction has been extensively studied, mainly using the enzyme from Escherichia coli. Bacillus subtilis contains two genes for TSs, ThyA and ThyB. The ThyB enzyme is very similar to other bacterial TSs, but the ThyA enzyme is quite different, both in sequence and activity. In ThyA TS, the active site histidine is replaced by valine. In addition, the B. subtilis enzyme has a 2.4-fold greater kcat than the E. coli enzyme. The structure of B. subtilis thymidylate synthase in a ternary complex with 5-fluoro-dUMP and 5,10-methylenetetrahydrofolate has been determined to 2.5 Å resolution. Overall, the structure of B. subtilis TS (ThyA) is similar to that of the E. coli enzyme. However, there are significant differences in the structures of two loops, the dimer interface and the details of the active site. The effects of the replacement of histidine by valine and a serine to glutamine substitution in the active site area, and the addition of a loop over the carboxy terminus may account for the differences in kcat found between the two enzymes.

Quantum-chemical methods are used to shed light on the functional role of residues involved in the resistance of HIV-1 reverse transcriptase against nucleoside-analog drugs. Ab initio molecular dynamics simulations are carried out for models representing the adduct between the triphosphate substrate and the nucleoside binding site. The triphosphate is considered either deprotonated or protonated at the γ-position. Although the protonated form already experiences large rearrangements in the ps time scale, the fully deprotonated state exhibits a previously unrecognized low-barrier hydrogen bond between Lys65 and γ-phosphate. Absence of this interaction in Lys65→Arg HIV-1 RT might play a prominent role in the resistance of this mutant for nucleoside analogs (Gu Z et al., 1994b, Antimicrob Agents Chemother 38:275–281; Zhang D et al., 1994, Antimicrob Agents Chemother 38:282–287). Water molecules present in the active site, not detected in the X-ray structure, form a complex H-bond network. Among these waters, one may be crucial for substrate recognition as it bridges Gln151 and Arg72 with the β-phosphate. Absence of this stabilizing interaction in Gln151→Met HIV-1 RT mutant may be a key factor for the known drug resistance of this mutant toward dideoxy-type drugs and AZT (Shirasaka T et al., 1995, Proc Natl Acad Sci USA 92:2398–2402; Iversen AK et al., 1996, J Virol 70:1086–1090).