Computer simulations of protein models on lattices have been widely used as an aid in the study of protein folding process. Following the suggestion of Raghunathan and Jernigan (1997, Protein Sci 6:2072–2083) that the cubo-octahedral lattice can allow a more realistic representation of proteins than other lattices, we propose here the use of a new set of internal coordinates θ for the description of a protein model on this lattice. An easy procedure for the conversion of the θ coordinates to the Cartesian coordinates is also described. When the Contact Interaction approach, already proposed by us for simulations on square or cubic lattices, was applied to the cubo-octahedral lattice, the system obeyed the correct thermodynamics derived from the definition of energy. Thus, lattice simulations of protein models, in which secondary structure elements such as α-helices or β-strands can be easily identifiable, can be performed.

Two different recombinant constructs of the N-terminal domain in Escherichia coli DnaJ were uniformly labeled with nitrogen-15 and carbon-13. One, DnaJ(1–78), contains the complete “J-domain,” and the other, DnaJ(1–104), contains both the “J-domain” and a conserved “G/F” extension at the C-terminus. The three-dimensional structures of these proteins have been determined by heteronuclear NMR experiments. In both proteins the “J-domain” adopts a compact structure consisting of a helix-turn-helix-loop-helix-turn-helix motif. In contrast, the “G/F” region in DnaJ(1–104) does not fold into a well-defined structure. Nevertheless, the “G/F” region has been found to have an effect on the packing of the helices in the “J-domain” in DnaJ(1–104). Particularly, the interhelical angles between Helix IV and other helices are significantly different in the two structures. In addition, there are some local conformational changes in the loop region connecting the two central helices. These structural differences in the “J-domain” in the presence of the “G/F” region may be related to the observation that DnaJ(1–78) is incapable of stimulating the ATPase activity of the molecular chaperone protein DnaK despite evidence that sites mediating the binding of DnaJ to DnaK are located in the 1–78 segment.

The lipases from Rhizopus and Rhizomucor are members of the family of Mucorales lipases. Although they display high sequence homology, their stereoselectivity toward triradylglycerols (sn-2 substituted triacylglycerols) varies. Four different triradylglycerols were investigated, which were classified into two groups: flexible substrates with rotatable O′-C1′ ether or ester bonds adjacent to C2 of glycerol and rigid substrates with a rigid N′-C1′ amide bond or a phenyl ring in sn-2. Although Rhizopus lipase shows opposite stereopreference for flexible and rigid substrates (hydrolysis in sn-1 and sn-3, respectively), Rhizomucor lipase hydrolyzes both groups of triradylglycerols preferably in sn-1. To explain these experimental observations, computer-aided molecular modeling was applied to study the molecular basis of stereoselectivity. A generalized model for both lipases of the Mucorales family highlights the residues mediating stereoselectivity: (1) L258, the C-terminal neighbor of the catalytic histidine, and (2) G266, which is located in a loop contacting the glycerol backbone of a bound substrate. Interactions with triradylglycerol substrates are dominated by van der Waals contacts. Stereoselectivity can be predicted by analyzing the value of a single substrate torsion angle that discriminates between sn-1 and sn-3 stereopreference for all substrates and lipases investigated here. This simple model can be easily applied in enzyme and substrate engineering to predict Mucorales lipase variants and synthetic substrates with desired stereoselectivity.

The stability of proteins is known to be affected significantly in the presence of high concentration of salts and is highly pH dependent. Extensive studies have been carried out on the stability of proteins in the presence of simple electrolytes and evaluated in terms of preferential interactions and increase in the surface tension of the medium. We have carried out an in-depth study of the effects of a series of carboxylic acid salts: ethylene diamine tetra acetate, butane tetra carboxylate, propane tricarballylate, citrate, succinate, tartarate, malonate, and gluconate on the thermal stability of five different proteins that vary in their physico-chemical properties: RNase A, cytochrome c, trypsin inhibitor, myoglobin, and lysozyme. Surface tension measurements of aqueous solutions of the salts indicate an increase in the surface tension of the medium that is very strongly correlated with the increase in the thermal stability of proteins. There is also a linear correlation of the increase in thermal stability with the number of carboxylic groups in the salt. Thermal stability has been found to increase by as much as 22 °C at 1 M concentration of salt. Such a high thermal stability at identical concentrations has not been reported before. The differences in the heat capacities of denaturation, ΔCp for RNase A, deduced from the transition curves obtained in the presence of varying concentrations of GdmCl and that of carboxylic acid salts as a function of pH, indicate that the nature of the solvent medium and its interactions with the two end states of the protein control the thermodynamics of protein denaturation. Among the physico-chemical properties of proteins, there seems to be an interplay of the hydrophobic and electrostatic interactions that lead to an overall stabilizing effect. Increase in surface free energy of the solvent medium upon addition of the carboxylic acid salts appears to be the dominant factor in governing the thermal stability of proteins.

Gelsolin is an actin-binding protein (82 kDa) consisting of six repeated segments (S1–S6), each approximately 120 residues long. It interacts with phospholipids and we previously showed that phosphatidylinositol 4,5-bisphosphate promotes phosphorylation of gelsolin by the tyrosine kinase c-Src. We used a combination of different methods, such as thin-layer chromatography and anti-phosphotyrosine-agarose immunoprecipitation of phosphopeptides combined with matrix assisted laser desorption ionization-mass spectrometry (MALDI-MS) and post source decay (PSD) analysis, to identify the phosphorylation sites in gelsolin. The major phosphorylation site (Tyr438) was located in subdomain 4 (S4). Phosphorylation of gelsolin in the gelsolin-actin2 complex was inhibited by 90%. Gelsolin phosphorylation by c-Src in the presence of lysophosphatidic acid also revealed Tyr438 as the most prominent site. Additional minor sites were found using the anti-phosphotyrosine bead immunoprecipitation method followed by MALDI-MS and PSD analysis. These sites, representing ∼5% of the total phosphate incorporation, were identified as Tyr59, Tyr382, Tyr576, and Tyr624. Based on these results we generated antibodies which specifically recognize Tyr438 phosphorylated gelsolin.

Scolexin is a coagulation-provoking plasma protein induced in response to bacterial or viral infection of larval Manduca sexta, a large lepidopterous insect. Here we report the isolation and sequencing of two cDNA clones that code for scolexin isoforms sharing 80% sequence identity. The scolexin sequences have low but recognizable sequence similarity to members of the chymotrypsin family and represent a new subfamily of chymotrypsin-like serine proteinases. Comparison with known structures reveals the conservation of key catalytic residues and a possible specificity for small nonpolar residues. Most remarkable is the absence of a canonical activation peptide cleavage site. This suggests that the regulation of scolexin activity will involve a novel activation mechanism.

The metallo-β-lactamases require divalent cations such as zinc or cadmium for hydrolyzing the amide bond of β-lactam antibiotics. The crystal structure of the Zn2+-bound enzyme from Bacteroides fragilis contains a binuclear zinc center in the active site. A hydroxide, coordinated to both zinc atoms, is proposed as the moiety that mounts the nucleophilic attack on the carbonyl carbon atom of the β-lactam bond of the substrate. It was previously reported that the replacement of the active site Cys181 by a serine residue severely impaired catalysis while atomic absorption measurements indicated that binding of the two zinc ions remained intact. Contradicting data emerge from recent mass spectrometry results, which show that only a single zinc ion binds to the C181S metallo-β-lactamase. In the current study, the C181S mutant enzyme was examined at the atomic level by determining the crystal structure at 2.6 Å resolution. The overall structure of the mutant enzyme is the same as that of the wild-type enzyme. At the mutation site, the side chain of Ser181 occupies the same position as that of the side chain of Cys181 in the wild-type protein. One zinc ion, Zn1, is present in the crystal structure; however, the site of the second zinc ion, Zn2 is unoccupied. A water molecule is associated with Zn1, reminiscent of the hydroxide seen in the structure of the wild-type enzyme but farther from the metal. The position of the water molecule is off the plane of the carboxylate group of Asp103; therefore, the water molecule may be less nucleophilic than a water molecule which is coplanar with the carboxylate group.

Unlike bovine cationic trypsin, rat anionic trypsin retains activity at high pH. This alkaline stability has been attributed to stabilization of the salt bridge between the N-terminal Ile16 and Asp194 by the surface negative charge (Soman K, Yang A-S, Honig B, Fletterick R., 1989, Biochemistry 28:9918–9926). The formation of this salt bridge controls the conformation of the activation domain in trypsin. In this work we probe the structure of rat trypsinogen to determine the effects of the surface negative charge on the activation domain in the absence of the Ile16–Asp194 salt bridge. We determined the crystal structures of the rat trypsin-BPTI complex and the rat trypsinogen-BPTI complex at 1.8 and 2.2 Å, respectively. The BPTI complex of rat trypsinogen resembles that of rat trypsin. Surprisingly, the side chain of Ile16 is found in a similar position in both the rat trypsin and trypsinogen complexes, although it is not the N-terminal residue and cannot form the salt bridge in trypsinogen. The resulting position of the activation peptide alters the conformation of the adjacent autolysis loop (residues 142–153). While bovine trypsinogen and trypsin have similar CD spectra, the CD spectrum of rat trypsinogen has only 60% of the intensity of rat trypsin. This lower intensity most likely results from increased flexibility around two conserved tryptophans, which are adjacent to the activation domain. The NMR spectrum of rat trypsinogen contains high field methyl signals as observed in bovine trypsinogen. It is concluded that the activation domain of rat trypsinogen is more flexible than that of bovine trypsinogen, but does not extend further into the protein core.