Methacarn (methanol-Carnoy) fixation

Summary According to chemical data, methanol raises the shrinkage temperature of collagen significantly more than ethanol (86° C versus 70° C). Since increase of shrinkage temperature appears desirable in tissues to be embedded in paraffin, methanol was substituted for ethanol in Carnoy's fluid. This methanol-Carnoy mixture is referred to as methacarn solution. The fixation-embedding procedure was similar to that described in the study of Carnoy fixation. Methacarn-fixed sections showed little or no shrinkage and compared well with material fixed in Carnoy's or Zenker's fluid. Myofibrils, especially in endothelial and epithelial cells, were more prominent in methacarn- than in Carnoy-fixed tissues.A review of the chemical literature showed that methanol, ethanol and chloroform stabilize or even enhance helical conformations of proteins, presumably by strengthening of hydrogen bonds. Interference with hydrophobic bonds causes unfolding and/or structural rearrangements in globular proteins. The twin-helical structure of DNA collapses in alcoholic solutions. Hence, methacarn fixation can be expected to preserve the helical proteins in myofibrils and collagen, but the conformations of globular proteins and DNA will be significantly altered. Literature on conformational effects produced by fixatives used in electron microscopy was also reviewed. Glutaraldehyde and OsO₄ cause considerable loss of helix (22–29% and 39–66% respectively). KMnO₄ and glutaraldehyde followed by OsO₄ produce extensive transitions from helical to random-coil conformations similar to those seen in powerful denaturants such as 8 M urea. Evidently these fixatives are unsuitable for studies of helical proteins. In contrast ethylene glycol preserves helical conformations.     

Local and nonlocal interactions in globular proteins and mechanisms of alcohol denaturation.

How important are helical propensities in determining the conformations of globular proteins? Using the two-dimensional lattice model and two monomer types, H (hydrophobic) and P (polar), we explore both nonlocal interactions, through an HH contact energy, as developed in earlier work, and local interactions, through a helix energy, σ. By computer enumeration, the partition functions for short chains are obtained without approximation for the full range of both types of energy. When nonlocal interactions dominate, some sequences undergo coil-globule collapse to a unique native structure. When local interactions dominate, all sequences undergo helix–coil transitions. For two different conformational properties, the closest correspondence between the lattice model and proteins in the Protein Data Bank is obtained if the model local interactions are made small compared to the HH contact interaction, suggesting that helical propensities may be only weak determinants of globular protein structures in water. For some HP sequences, varying σ/ leads to additional sharp transitions (sometimes several) and to “conformational switching” between unique conformations. This behavior resembles the transitions of globular proteins in water to helical states in alcohols. In particular, comparison with experiments shows that whereas urea as a denaturant is best modeled as weakening both local and nonlocal interactions, trifluoroethanol is best modeled as mainly weakening HH interactions and slightly enhancing local helical interactions.     

Excluded volume in protein side-chain packing

The excluded volume occupied by protein side-chains and the requirement of high packing density in the protein interior should severely limit the number of side-chain conformations compatible with a given native backbone. To examine the relationship between side-chain geometry and side-chain packing, we use an all-atom Monte Carlo simulation to sample the large space of side-chain conformations. We study three models of excluded volume and use umbrella sampling to effectively explore the entire space. We find that while excluded volume constraints reduce the size of conformational space by many orders of magnitude, the number of allowed conformations is still large. An average repacked conformation has 20 % of its chi angles in a non-native state, a marked reduction from the expected 67 % in the absence of excluded volume. Interestingly, well-packed conformations with up to 50 % non-native chi angles exist. The repacked conformations have native packing density as measured by a standard Voronoi procedure. Entropy is distributed non-uniformly over positions, and we partially explain the observed distribution using rotamer probabilities derived from the Protein Data Bank database. In several cases, native rotamers that occur infrequently in the database are seen with high probability in our simulation, indicating that sequence-specific excluded volume interactions can stabilize rotamers that are rare for a given backbone. In spite of our finding that 65 % of the native rotamers and 85 % of chi(1) angles can be predicted correctly on the basis of excluded volume only, 95 % of positions can accommodate more than one rotamer in simulation. We estimate that, in order to quench the side-chain entropy observed in the presence of excluded volume interactions, other interactions (hydrophobic, polar, electrostatic) must provide an additional stabilization of at least 0.6 kT per residue in order to single out the native state.     

Statistical and conformational analysis of the electron density of protein side chains

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.     

Rotamer strain as a determinant of protein structural specificity

We present direct evidence for a change in protein structural specificity due to hydrophobic core packing. High resolution structural analysis of a designed core variant of ubiquitin reveals that the protein is in slow exchange between two conformations. Examination of side-chain rotamers indicates that this dynamic response and the lower stability of the protein are coupled to greater strain and mobility in the core. The results suggest that manipulating the level of side-chain strain may be one way of fine tuning the stability and specificity of proteins.     

Side-chain rotamers in protein α-α hairpins and a mechanism of their selection

Several regularities were observed for the distribution of side-chain rotamers in α-α hairpins of globular proteins. In left-turned α-α hairpins, most side chains adopt t rotamers in d-positions and g^? rotamers in g-positions. In right-turned α-α hairpins, most side-chains adopt g^? rotamers in a-positions and t rotamers in e-positions. Analysis of these regularities suggested that selection of the side-chain conformation in α-α hairpins depends on two main factors, the mode of α-helix packing and the positions of side chains in α-helices. The regularities were explained by the squeezing mechanism: interhelical interactions bring the α-helices close together so that the side chains are squeezed out of the helix-helix interface and adopt unique conformations.     

A designed four-alpha-helix bundle that binds the volatile general anesthetic halothane with high affinity

The structural features of volatile anesthetic binding sites on proteins are being examined with the use of a defined model system consisting of a four-alpha-helix bundle scaffold with a hydrophobic core. Previous work has suggested that introducing a cavity into the hydrophobic core improves anesthetic binding affinity. The more polarizable methionine side chain was substituted for a leucine, in an attempt to enhance the dispersion forces between the ligand and the protein. The resulting bundle variant has an improved affinity (K(d) = 0.20 +/- 0.01 mM) for halothane binding, compared with the leucine-containing bundle (K(d) = 0.69 +/- 0.06 mM). Photoaffinity labeling with (14)C-halothane reveals preferential labeling of the W15 residue in both peptides, supporting the view that fluorescence quenching by bound anesthetic reports both the binding energetics and the location of the ligand in the hydrophobic core. The rates of amide hydrogen exchange were similar for the two bundles, suggesting that differences in binding affinity were not due to changes in protein stability. Binding of halothane to both four-alpha-helix bundle proteins stabilized the native folded conformations. Molecular dynamics simulations of the bundles illustrate the existence of the hydrophobic core, containing both W15 residues. These results suggest that in addition to packing defects, enhanced dispersion forces may be important in providing higher affinity anesthetic binding sites. Alternatively, the effect of the methionine substitution on halothane binding energetics may reflect either improved access to the binding site or allosteric optimization of the dimensions of the binding pocket. Finally, preferential stabilization of folded protein conformations may represent a fundamental mechanism of inhaled anesthetic action.     

Conformational stability of a thrombin-binding peptide derived from the hirudin C-terminus.

The COOH-terminal region of hirudin represents an independent functional domain that binds to an anion-binding exosite of thrombin and inhibits the interaction of thrombin with fibrinogen and regulatory proteins in blood coagulation. The thrombin-bound structure of the peptide fragment, hirudin 55-65, has been determined by use of transferred NOE spectroscopy [Ni, F., Konishi, Y., & Scheraga, H. A. (1990) Biochemistry 29, 4479-4489]. The stability of the thrombin-bound conformation has been characterized further by a combined NMR and theoretical analysis of the conformational ensemble accessible by the hirudin peptide. Medium- and long-range NOE's were found for the free hirudin peptide in aqueous solution and in a mixture of dimethyl sulfoxide and water at both ambient (25 degrees C) and low (0 degrees C) temperatures, suggesting that ordered conformations are highly populated in solution. The global folding of these conformations is similar to that in the thrombin-bound state, as indicated by NOE's involving the side-chain protons of residues Phe(56), Ile(59), Pro(60), Tyr(63), and Leu(64). Residues Glu(61), Glu(62), Tyr(63), and Leu(64) all contain approximately 50% of helical conformations calculated from the ratio of the sequential dNN and d alpha N NOE's. Among the helical ensemble, active 3(10)-helical conformations were found by an analysis of the medium-range [(i,i+2) and (i,i+3)] NOE's involving the last six residues of the peptide. An analysis of the side-chain rotamers revealed that, upon binding to thrombin, there may be a rotation around the alpha CH-beta CH bond of Ile(59) such that Ile(59) adopts a gauche- (chi 1 = +60) conformation in contrast to the highly populated trans (chi 1 = -60) found for Ile(59) in the free peptide. However, the thrombin-bound conformation of the hirudin peptide is still an intrinsically stable conformer, and the preferred conformational ensemble of the peptide contains a large population of the active conformation. The apparent preference for a gauche- (chi 1 = +60) side-chain conformation of Ile(59) in the bound state may be explained by the existence of a positively charged arginine residue among the hydrophobic residues in the thrombin exosite.     

Bayesian statistical analysis of protein side-chain rotamer preferences.

We present a Bayesian statistical analysis of the conformations of side chains in proteins from the Protein Data Bank. This is an extension of the backbone-dependent rotamer library, and includes rotamer populations and average chi angles for a full range of phi, psi values. The Bayesian analysis used here provides a rigorous statistical method for taking account of varying amounts of data. Bayesian statistics requires the assumption of a prior distribution for parameters over their range of possible values. This prior distribution can be derived from previous data or from pooling some of the present data. The prior distribution is combined with the data to form the posterior distribution, which is a compromise between the prior distribution and the data. For the chi 2, chi 3, and chi 4 rotamer prior distributions, we assume that the probability of each rotamer type is dependent only on the previous chi rotamer in the chain. For the backbone-dependence of the chi 1 rotamers, we derive prior distributions from the product of the phi-dependent and psi-dependent probabilities. Molecular mechanics calculations with the CHARMM22 potential show a strong similarity with the experimental distributions, indicating that proteins attain their lowest energy rotamers with respect to local backbone-side-chain interactions. The new library is suitable for use in homology modeling, protein folding simulations, and the refinement of X-ray and NMR structures.     

Hydrated and Dehydrated Tertiary Interactions-Opening and Closing-of a Four-Helix Bundle Peptide

The structural heterogeneity and thermal denaturation of a dansyl-labeled four-helix bundle homodimeric peptide was studied with steady-state and time-resolved fluorescence spectroscopy and with circular dichroism (CD). At room temperature the fluorescence decay of the polarity-sensitive dansyl, located in the hydrophobic core region, can be described by a broad distribution of fluorescence lifetimes, reflecting the heterogeneous microenvironment. However, the lifetime distribution is nearly bimodal, which we ascribe to the presence of two major conformational subgroups. Since the fluorescence lifetime reflects the water content of the four-helix bundle conformations, we can use the lifetime analysis to monitor the change in hydration state of the hydrophobic core of the four-helix bundle. Increasing the temperature from 9°C to 23°C leads to an increased population of molten-globule-like conformations with a less ordered helical backbone structure. The fluorescence emission maximum remains constant in this temperature interval, and the hydrophobic core is not strongly affected. Above 30°C the structural dynamics involve transient openings of the four-helix bundle structure, as evidenced by the emergence of a water-quenched component and less negative CD. Above 60°C the homodimer starts to dissociate, as shown by the increasing loss of CD and narrow, short-lived fluorescence lifetime distributions.     

Using quantum mechanics to improve estimates of amino acid side chain rotamer energies

Amino acid side chains adopt a discrete set of favorable conformations typically referred to as rotamers. The relative energies of rotamers partially determine which side chain conformations are more often observed in protein structures and accurate estimates of these energies are important for predicting protein structure and designing new proteins. Protein modelers typically calculate side chain rotamer energies by using molecular mechanics (MM) potentials or by converting rotamer probabilities from the protein database (PDB) into relative free energies. One limitation of the knowledge‐based energies is that rotamer preferences observed in the PDB can reflect internal side chain energies as well as longer‐range interactions with the rest of the protein. Here, we test an alternative approach for calculating rotamer energies. We use three different quantum mechanics (QM) methods (second order Møller‐Plesset (MP2), density functional theory (DFT) energy calculation using the B3LYP functional, and Hartree‐Fock) to calculate the energy of amino acid rotamers in a dipeptide model system, and then use these pre‐calculated values in side chain placement simulations. Energies were calculated for over 36,000 different conformations of leucine, isoleucine, and valine dipeptides with backbone torsion angles from the helical and strand regions of the Ramachandran plot. In a subset of cases these energies differ significantly from those calculated with standard molecular mechanics potentials or those derived from PDB statistics. We find that in these cases the energies from the QM methods result in more accurate placement of amino acid side chains in structure prediction tests. Proteins 2008. © 2007 Wiley‐Liss, Inc.     

NMR solution structure of ImB2, a protein conferring immunity to antimicrobial activity of the type IIa bacteriocin, carnobacteriocin B2

Bacteriocins produced by lactic acid bacteria are potent antimicrobial compounds which are active against closely related bacteria. Producer strains are protected against the effects of their cognate bacteriocins by immunity proteins that are located on the same genetic locus and are coexpressed with the gene encoding the bacteriocin. Several structures are available for class IIa bacteriocins; however, to date, no structures are available for the corresponding immunity proteins. We report here the NMR solution structure of the 111-amino acid immunity protein for carnobacteriocin B2 (ImB2). ImB2 folds into a globular domain in aqueous solution which contains an antiparallel four-helix bundle. Extensive packing by hydrophobic side chains in adjacent helices forms the core of the protein. The C-terminus, containing a fifth helix and an extended strand, is held against the four-helix bundle by hydrophobic interactions with helices 3 and 4. Most of the charged and polar residues in the protein face the solvent. Helix 3 is well-defined to residue 55, and a stretch of nascent helix followed by an unstructured loop joins it to helix 4. No interaction is observed between ImB2 and either carnobacteriocin B2 (CbnB2) or its precursor. Protection from the action of CbnB2 is only observed when ImB2 is expressed within the cell. The loop between helices 3 and 4, and a hydrophobic pocket which it partially masks, may be important for interaction with membrane receptors responsible for sensitivity to class IIa bacteriocins.     

Boltzmann-type distribution of side-chain conformation in proteins

We analyze packing imperfections in globular proteins as reflected in deviations of torsion angles from the equilibrium values for the isolated side chains. The distribution of conformations of methionine and lysine residues in a database of high-resolution structures is compared with energies of model compounds calculated with high-level quantum-mechanics. The distribution of the C-C and C-S torsion angles (chi(3)) correlates well with the Boltzmann factor of the torsion energy, exp(-betaE) of the model compounds C(2)H(5)-C(2)H(5) and C(2)H(5)-S-CH(3). An exponential relation was again found between the relative occurrence of g+, g- and t conformations for C(alpha)-C(beta) bonds in long side chains and the energy differences of rotamers of alpha-amino n-butyric acid, when dependence on backbone conformation was taken into account. The distribution of all 27 rotamers of methionine was correlated with the energy differences between the model's rotamers, corrected for clashes with nearby residues, the correlation being good for a set with backbone in the beta-conformation, but less clear for backbone alpha-conformation. In all correlations, the value of the coefficient beta corresponds to a temperature of circa 300 K. These results can be interpreted with a model that considers the structure of a folded protein as resulting from packing imperfectly complementary parts, with a requirement of an overall low energy. Compromises are required to optimize the fit of nonbonded contacts with surrounding groups, and side chains assume conformations away from the energy minimum. An exponential distribution is a most probable distribution, and this can be established easily under conditions other than thermal equilibrium.     

Configurational-bias sampling technique for predicting side-chain conformations in proteins

Prediction of side-chain conformations is an important component of several biological modeling applications. In this work, we have developed and tested an advanced Monte Carlo sampling strategy for predicting side-chain conformations. Our method is based on a cooperative rearrangement of atoms that belong to a group of neighboring side-chains. This rearrangement is accomplished by deleting groups of atoms from the side-chains in a particular region, and regrowing them with the generation of trial positions that depends on both a rotamer library and a molecular mechanics potential function. This method allows us to incorporate flexibility about the rotamers in the library and explore phase space in a continuous fashion about the primary rotamers. We have tested our algorithm on a set of 76 proteins using the all-atom AMBER99 force field and electrostatics that are governed by a distance-dependent dielectric function. When the tolerance for correct prediction of the dihedral angles is a <20 degrees deviation from the native state, our prediction accuracies for chi1 are 83.3% and for chi1 and chi2 are 65.4%. The accuracies of our predictions are comparable to the best results in the literature that often used Hamiltonians that have been specifically optimized for side-chain packing. We believe that the continuous exploration of phase space enables our method to overcome limitations inherent with using discrete rotamers as trials.     

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.     

The Hydrophobicity of the H3 Histone Fold differs from the Hydrophobicity of the other three Folds

The eukaryotic histone dimers, H3–H4 and H2A–H2B, are formed in the cytosol prior to being transported into the nucleus and assembled into the nucleosome. Residue side-chain distances from the interior of the histone dimers are obtained with an ellipsoidal spatial metric and structural information provided by X-ray analyses at atomic resolution of the nucleosome core particles. While the spatial hydrophobic moment profiles of the dimers are comparable with profiles obtained previously that characterize the hydrophobic core of single-chain, single-domain globular soluble proteins, correlation coefficients between the side-chain hydrophobicities and distances from the interior of the H3–H4 dimer and H2A–H2B dimer differ significantly. This difference is traced to the H3 histone fold, which segregates fewer hydrophobic residues within the protein interior than the three other folds. Examination of the correlation coefficient between residue hydrophobicity and side-chain distance from the dimer interior over local regions of the fold sequence shows that the region of reduced correlation is associated mainly with the residues at the carboxyl end of the H3 histone fold, the helical region of the fold involved in the H3–H3 binding of the (H3–H4)₂ tetramer of the nucleosome. Hydrophobic interactions apparently contribute to the binding of this fourfold helical bundle and this evolutionary requirement may trade off against the requirement for H3–H4 dimer stability. The present results provide a different view than previously proposed, albeit of similar origin, to account for the reduced stability of the H3–H4 dimer compared with the H2A–H2B dimer.Reviewing Editor: Dr. Martin Kreitman     

The focal adhesion targeting (FAT) region of focal adhesion kinase is a four-helix bundle that binds paxillin.

Focal adhesion kinase (FAK) is a tyrosine kinase found in focal adhesions, intracellular signaling complexes that are formed following engagement of the extracellular matrix by integrins. The C-terminal 'focal adhesion targeting' (FAT) region is necessary and sufficient for localizing FAK to focal adhesions. We have determined the crystal structure of FAT and show that it forms a four-helix bundle that resembles those found in two other proteins involved in cell adhesion, alpha-catenin and vinculin. The binding of FAT to the focal adhesion protein, paxillin, requires the integrity of the helical bundle, whereas binding to another focal adhesion protein, talin, does not. We show by mutagenesis that paxillin binding involves two hydrophobic patches on opposite faces of the bundle and propose a model in which two LD motifs of paxillin adopt amphipathic helices that augment the hydrophobic core of FAT, creating a six-helix bundle.     

Iterative assembly of helical proteins by optimal hydrophobic packing

We present a method for the computer-based iterative assembly of native-like tertiary structures of helical proteins from alpha-helical fragments. For any pair of helices, our method, called MATCHSTIX, first generates an ensemble of possible relative orientations of the helices with various ways to form hydrophobic contacts between them. Those conformations having steric clashes, or a large radius of gyration of hydrophobic residues, or with helices too far separated to be connected by the intervening linking region, are discarded. Then, we attempt to connect the two helical fragments by using a robotics-based loop-closure algorithm. When loop closure is feasible, the algorithm generates an ensemble of viable interconnecting loops. After energy minimization and clustering, we use a representative set of conformations for further assembly with the remaining helices, adding one helix at a time. To efficiently sample the conformational space, the order of assembly generally proceeds from the pair of helices connected by the shortest loop, followed by joining one of its adjacent helices, always proceeding with the shorter connecting loop. We tested MATCHSTIX on 28 helical proteins each containing up to 5 helices and found it to heavily sample native-like conformations. The average rmsd of the best conformations for the 17 helix-bundle proteins that have 2 or 3 helices is less than 2 A; errors increase somewhat for proteins containing more helices. Native-like states are even more densely sampled when disulfide bonds are known and imposed as restraints. We conclude that, at least for helical proteins, if the secondary structures are known, this rapid rigid-body maximization of hydrophobic interactions can lead to small ensembles of highly native-like structures. It may be useful for protein structure prediction.