Analysis of the conformation of polypeptides : the combined use of energy computations and nuclear magnetic resonance studies.

This review surveys current approaches to the problem of determining the solution conformation of polypeptides. The basic principles of energy computations are described. The utility, problems, and limitations of various theoretical methods are summarized: conformational energy mapping, energy minimization, scanning of selected local conformations, statistical predictive schemes. The need for combining the calculations with experimental studies is pointed out. The information content of various physico-chemical methods is compared for this purpose. The analysis of nmr coupling constants is discussed in more detail. In combination with energy computations, it can furnish specific information on local aspects of the conformation. Examples of such combined studies on small peptides are summarized.     

Prediction of the conformation of ice-nucleation protein by conformational energy calculation

The active conformation of an ice-nucleation protein, whose major portion consists of a long polypeptide segment of nearly repetitive octapeptides, is predicted by the analyses of conformational energy and the mechanism of crystal growth. The protein ideally has an exact octapeptide repetition and is assumed to have a helical conformation. The present study searched for low-energy helical conformations and each of the obtained low-energy conformations examined as to whether it has a surface structure that can promote crystal formation. Two conformations obtained were good candidates for an ice nucleus. Both were found to have on their surfaces an arrangement of hydrogen-bonding sites, which fits well with those of hydrogen bonds in hexagonal ice crystal. Further, one of the two conformations had a hexagonal conformational symmetry consistent with the hexagonal ice crystal structure. The other conformation had a pentagonal conformational symmetry that could enable the growth of an ice crystal--dendritic polycrystalline snow crystal--which grows on metastable cubic ice.     

Theoretical studies on the conformation of monosialogangliosides and disialogangliosides

The preferred conformation of gangliosides GM3, GM2, GM1, GD1a and GD1b have been studied by computing their potential energies. The conformation of NeuNAc in GM3 differs from that expected for the same residue in GM2 and GM1. The NeuNAc residues in GM2 and GM1 exhibit identical conformations. Theory predicts that the terminal NeuNAc of GD1a is conformationally similar to that of GM3 and that the internal one is similar in conformation to those present in GM2 and GM1 in agreement with NMR studies. The differences in chemical shifts of the C2 and C3 carbons of the internal and terminal NeuNAc of GD1a have been attributed to differences in orientation. The present studies suggest that the binding site of cholera toxin is much smaller than that of tetanus toxin. The preferred shape of these gangliosides correlate well with their biological properties.     

Gel electrophoresis in studies of protein conformation and folding

Electrophoresis through polyacrylamide gels is a useful method for distinguishing conformational states of proteins and analyzing the thermodynamic and kinetic properties of transitions between conformations. Although the relationship between protein conformation and electrophoretic mobility is quite complex, relative mobilities provide qualitative estimates of compactness. Conformational states which interconvert slowly on the time scale of the electrophoretic separation can often be resolved, and the rates of interconversion can be estimated. If the transitions are more rapid, then the electrophoretic mobility represents the equilibrium distribution of conformations. Protein unfolding transitions induced by urea are readily studied using slab gels containing a gradient of urea concentration perpendicular to the direction of electrophoresis. Protein applied across the top of such a gel migrates in the presence of continuously varying urea concentrations, and a profile of the unfolding transition is generated directly. Transitions induced by other agents could be studied using analogous gradient gels. Electrophoretic methods are especially suited for studying small quantities of protein, and complex mixtures, since the different components can be separated during the electrophoresis.     

Alteration of the DNA double helix conformation upon incorporation of mispairs as revealed by energy computations and pathways of point mutations.

To explain biochemical and genetic data on spontaneous nucleotide replacements in nucleic acid biosynthesis all the 8 mispairs in normal tautomeric forms have been considered. Possible B-conformations of DNA fragments containing each of such mispairs incorporated between Watson-Crick pairs have been found using computations of the energy of non-bonded interactions via classical potential functions. These conformations have no reduced interatomic contacts. The values of each dihedral angle of the sugar-phosphate backbone fall within the limits of those of double-helical fragments of B-DNA in crystals. These values differ from those of the corresponding angles for the low-energy polynucleotide conformations consisting of canonical pairs by no more than 30 degrees (except for the fragment with the U:U pair for which the C4'-C3'-O-P angle differs by about 50 degrees). The difference in experimentally observed frequencies of various nucleotide replacements in DNA biosynthesis correlates with the difference in the energy of non-bonded interactions and with the extent of the sugar-phosphate backbone distortion for the fragments containing the mispairs which serve as intermediates for the replacements.     

Linkage studies of Friedreich ataxia by means of blood-group and protein markers

Friedreich ataxia (FA) is an autosomal recessive, neuro-degenerative disorder in which the pathogenetic mechanism remains unidentified despite extensive biochemical studies. Genetic-linkage studies provide an alternative approach to determining the basic defect. Linkage analysis between FA and 36 polymorphic-blood-group and protein markers has been carried out on three separate patient populations--16 families from the inbred Acadian population of Louisiana, 21 French-Canadian families from Quebec, and nine apparently unrelated British families--in an attempt to determine the chromosomal location of the disease mutation. Neither evidence of linkage to any of the markers investigated nor heterogeneity among the populations was found for any of the comparisons. The negative lod scores exclude the locus for FA from greater than 20% of the genome.     

Kinetic studies of protein farnesyltransferase mutants establish active substrate conformation

The zinc metalloenzyme protein farnesyltransferase (FTase) catalyzes the transfer of a 15-carbon farnesyl moiety from farnesyl diphosphate (FPP) to a cysteine residue near the C-terminus of a protein substrate. Several crystal structures of inactive FTase.FPP.peptide complexes indicate that K164alpha interacts with the alpha-phosphate and that H248beta and Y300beta form hydrogen bonds with the beta-phosphate of FPP [Strickland, C. L., et al. (1998) Biochemistry 37, 16601-16611]. Mutations K164Aalpha, H248Abeta, and Y300Fbeta were prepared and analyzed by single turnover kinetics and ligand binding studies. These mutations do not significantly affect the enzyme affinity for FPP but do decrease the farnesylation rate constant by 30-, 10-, and 500-fold, respectively. These mutations have little effect on the pH and magnesium dependence of the farnesylation rate constant, demonstrating that the side chains of K164alpha, Y300beta, and H248beta do not function either as general acid-base catalysts or as magnesium ligands. Mutation of H248beta and Y300beta, but not K164alpha, decreases the farnesylation rate constant using farnesyl monophosphate (FMP). These data suggest that, contrary to the conclusions derived from analysis of the static crystal structures, the transition state for farnesylation is stabilized by interactions between the alpha-phosphate of the isoprenoid substrate and the side chains of Y300beta and H248beta. These results suggest an active substrate conformation for FTase wherein the C1 carbon of the FPP substrate moves toward the zinc-bound thiolate of the protein substrate to react, resulting in a rearrangement of the diphosphate group relative to its ground state position in the binding pocket.     

Morphogenesis studies by means of X-rays

Development Genes and Evolution -     

Theoretical conformation analysis of MCD peptide

The spatial structure of the MCD-peptide from bee venom has been calculated basing on the known sequence of 22 amino acid. The a priori calculations produce a system of two disulfide bonds, identical to that observed in the native structure. The calculated structure of MCD-peptide is close to that proposed earlier for the homologues peptide tertiapin and is confirmed by NMR and CD data.     

Predicting protein conformation by statistical methods.

The unique folded structure makes a polypeptide a functional protein. The number of known sequences is about a hundred times larger than the number of known structures and the gap is increasing rapidly. The primary goal of all structure prediction methods is to obtain structure-related information on proteins, whose structures have not been determined experimentally. Besides this goal, the development of accurate prediction methods helps to reveal principles of protein folding. Here we present a brief survey of protein structure predictions based on statistical analyses of known sequence and structure data. We discuss the background of these methods and attempt to elucidate principles, which govern structure formation of soluble and membrane proteins.     

Theoretical study of the conformation of the lipoamide arm in a mutant H protein

The lipoamide arm of the H protein plays a pivotal role in the catalytic cycle of the glycine decarboxylase complex (GDC) by being successively methylamine loaded (Hmet), reduced (Hred), and oxidized (Hox). In a previous study, we calculated free-energy surfaces as a function of the lipoamide arm position of the three forms of the wild-type protein and found close agreement with the available experimental data. Our simulations, together with crystallographic and NMR data, showed that the methylamine-loaded arm is locked in a cavity by interaction with Ser12, Glu14, and Asp67. In this work, we investigate the behavior of the methylamine-loaded form of a mutant H protein (HEA) where Glu14 has been replaced by Ala. We find that the arm can still be held in the cavity but that the energy barrier to release of the arm is halved from approximately 40 kcal mol(-1) for Hmet to approximately 12 kcal mol(-1) for HEA. To compensate for the loss of Glu14, the methylamine group shifts toward Ser66 in the mutant form. These results provide a structural basis for the equilibrium between the loaded and the unloaded forms of the arm observed by Gueguen et al. (Gueguen et al., J Biol Chem 1999;274:26344-26352) in HEA.     

Prediction of protein side-chain conformation by packing optimization

We have developed a rapid and completely automatic method for prediction of protein side-chain conformation, applying the simulated annealing algorithm to optimization of side-chain packing (van der Waals) interactions. The method directly attacks the combinatorial problem of simultaneously predicting many residues' conformation, solving in 8 to 12 hours problems for which the systematic search would require over 10(300) central processing unit years. Over a test set of nine proteins ranging in size from 46 to 323 residues, the program's predictions for side-chain atoms had a root-mean-square (r.m.s.) deviation of 1.77 A overall versus the native structures. More importantly, the predictions for core residues were especially accurate, with an r.m.s. value of 1.25 A overall: 80 to 90% of the large hydrophobic side-chains dominating the internal core were correctly predicted, versus 30 to 40% for most current methods. The predictions' main errors were in surface residues poorly constrained by packing and small residues with greater steric freedom and hydrogen bonding interactions, which were not included in the program's potential function. van der Waals interactions appear to be the supreme determinant of the arrangement of side-chains in the core, enforcing a unique allowed packing that in every case so far examined matches the native structure.     

Methods of peptide conformation studies

In solution most of the peptides assume multiple flexible conformations. Determination of the dominant conformers and evaluation of their populations is the aim of peptide conformation studies, in which theoretical and experimental methods play complementary roles. Molecular dynamics or Monte Carlo methods are quite effective in searching the conformational space accessible to a peptide but they are not able to estimate, precisely enough, the populations of various conformations. Therefore, they must be supplemented by experimental data. In this paper, a short review of the experimental methods, most widely used in peptide conformational studies, is presented. Among them NMR plays the leading role. Valuable information is also obtained from hydrogen exchange, fluorescence resonance energy transfer, and circular dichroism measurements. The advantages and shortcomings of these methods are discussed.