A novel exhaustive search algorithm for predicting the conformation of polypeptide segments in proteins

We present a fast ab initio method for the prediction of local conformations in proteins. The program, PETRA, selects polypeptide fragments from a computer-generated database (APD) encoding all possible peptide fragments up to twelve amino acids long. Each fragment is defined by a representative set of eight straight phi/psi pairs, obtained iteratively from a trial set by calculating how fragments generated from them represent the protein databank (PDB). Ninety-six percent (96%) of length five fragments in crystal structures, with a resolution better than 1.5 A and less than 25% identity, have a conformer in the database with less than 1 A root-mean-square deviation (rmsd). In order to select segments from APD, PETRA uses a set of simple rule-based filters, thus reducing the number of potential conformations to a manageable total. This reduced set is scored and sorted using rmsd fit to the anchor regions and a knowledge-based energy function dependent on the sequence to be modelled. The best scoring fragments can then be optimized by minimization of contact potentials and rmsd fit to the core model. The quality of the prediction made by PETRA is evaluated by calculating both the differences in rmsd and backbone torsion angles between the final model and the native fragment. The average rmsd ranges from 1.4 A for three residue loops to 3.9 A for eight residue loops.     

Conformational analysis of polypeptides and proteins for the study of protein folding, molecular recognition, and molecular design

Conformational energy calculations provide an understanding as to how interatomic interactions lead to the three-dimensional structures of polypeptides and proteins, and how these molecules interact with other molecules. Illustrative results of such calculations pertain to model systems (-helices and -sheets, and interactions between them), to various open-chain and cyclic peptides, to fibrous proteins, to globular proteins, and to enzyme-substrate complexes. In most cases, the validity of the computations is established by experimental tests of the predicted structures.This article was presented during the proceedings of the International Conference on Macromolecular Structure and Function, held at the National Defence Medical College, Tokorozawa, Japan, December 1985. This paper first appeared in the Israel Journal of Chemistry, Vol. 27, 1986.     

Benchmarking of programs that predict the position of transmembrane segments in beta-barrel proteins

We propose a method for comparative analysis of the programs that locate transmembrane segments in proteins. On this basis, we have compiled a sample of beta-barrel protein that allows unequivocal assessment of prediction quality. Upon testing of several Internet servers, B2TMR proves the best, with B2TMPRED, HMM-B2TMR, and PROFtmb following.     

Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety‐six of these alignment positions correspond to conserved Gly residue with (φ, ψ) values allowed for non‐glycyl residues. Reasons for this observation were investigated by in‐silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C^β atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (φ, ψ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (φ, ψ) values which are disallowed for Ala. In‐silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C^β atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.     

Surface interactions of gamma-crystallins in the crystal medium in relation to their association in the eye lens

A comparative study of intermolecular interactions in crystals of two homologous low molecular weight proteins, gamma-II and gamma-IIIb crystallins, from calf eye lens was carried out. Crystal packings for these proteins are very different: intermolecular contact areas compose about 33% of the total accessible surface area of gamma-II as compared with 13% in gamma-III. Two key residues seem to be mainly responsible for the differences in protein association in the crystal medium. These are Ser 103 and Leu 155 in gamma-II, which are replaced by Met 103 and His 155 in gamma-IIb. A similar substitution of these residues is observed in different gene products of gamma-crystallins from a number of vertebrates. This is consistent with the existence of a genetically controlled mechanism for determining intermolecular association of gamma-crystallins in the native medium of the lens.     

Modeling of peptides and proteins in a membrane environment: I. A solvation model mimicking a lipid bilayer

A theoretical solvation model of peptides and proteins that mimics the heterogeneous membrane-water system was proposed. Our approach is based on the combined use of atomic parameters of solvation for water and hydrocarbons, which approximates the hydrated polar groups and acyl chains of lipids, respectively. This model was tested in simulations of several peptides: a nonpolar 20-mer polyleucine, a hydrophobic peptide with terminal polar groups, and a strongly amphiphilic peptide. The conformational space of the peptides in the presence of the membrane was studied by the Monte Carlo method. Unlike a polar solvent and vacuum, the membrane-like environment was shown to stabilize the α-helical conformation: low-energy structures have a helicity index of 100% in all cases. At the same time, the energetically most favorable orientations of the peptides relative to the membrane depend on their hydrophobic properties: nonpolar polyleucine is entirely immersed in the bilayer and the hydrophobic peptide with polar groups at the termini adopts a transbilayer orientation, whereas the amphiphilic peptide lies at the interface parallel to the membrane plane. The results of the simulations agree well with the available experimental data for these systems. In the following communications of this series, we plan to describe applications of the solvation model to membrane-bound proteins and peptides with biologically important functional activities.     

Variability of the canonical loop conformations in serine proteinases inhibitors and other proteins

Canonical loops of protein inhibitors of serine proteinases occur in proteins having completely different folds. In this article, conformations of the loops have been analyzed for inhibitors belonging to 10 structurally different families. Using deviation in C_α-C_α distances as a criterion for loop similarity, we found that the P3-P3′ segment defines most properly the length of the loop. When conformational differences among loops of individual inhibitors were compared using root mean square deviation (rmsd) in atomic coordinates for all main chain atoms (Δr method) and rmsd operating in main chain torsion angles (Δt method), differences of up to 2.1 Å and 72.3°, respectively, were observed. Such large values indicate significant conformational differences among individual loops. Nevertheless, the overall geometry of the inhibitor-proteinase interaction is very well preserved, as judged from the similarity of C_α-C_α distances between C_α of catalytic Ser and C_α of P3-P3′ residues in various enzyme-inhibitor complexes. The mode of interaction is very well preserved both in the chymotrypsin and subtilisin families, as distances calculated for subtilisin-inhibitor complexes are almost always within the range of those for chymotrypsin-inhibitor complexes. Complex formation leads to conformational changes of up to 160° for χ₁ angle. Side chains of residue P1 and P2′ adopt in different complexes a similar orientation (χ₁ angle = −60° and −180°, respectively). To check whether the canonical conformation can be found among non–proteinase-inhibitor Brookhaven Protein Data Bank structures, two selection criteria—the allowed main chain dihedral angles and C_α-C_α distances for the P3-P3′ segment—were applied to all these structures. This procedure detected 132 unique hexapeptide segments in 121 structurally and functionally unrelated proteins. Partial preferences for certain amino acids occurring at particular positions in these hexapeptides could be noted. Further restriction of this set to hexapeptides with a highly exposed P1 residue side chain resulted in 40 segments. The possibility of complexes formation between these segments and serine proteinases was ruled out in molecular modeling due to steric clashes. Several structural features that determine the canonical conformation of the loop both in inhibitors and in other proteins can be distinguished. They include main chain hydrogen bonds both within the P3-P3′ segment and with the scaffold region, P3-P4 and P3′-P4′ hydrophobic interactions, and finally either hydrophobic or polar interactions involving the P1′ residue. Proteins 32:459–474, 1998. © 1998 Wiley-Liss, Inc.     

NMR spectroscopy as a tool for the rapid assessment of the conformation of GST-fusion proteins

Glutathione-S-transferase (GST)-fusion proteins are used extensively for structural, biochemical, and functional analyses. Although the conformation of the target protein is of critical importance, confirmation of the folded state of the target is often not undertaken or is cumbersome because of the requirement to first remove the GST tag. Here, we demonstrate that it is possible to record conventional (15)N-HSQC NMR spectra of small GST-fusion proteins and that the observed signals arise almost exclusively from the target protein. This approach constitutes a rapid and straightforward means of assessing the conformation of a GST-fusion protein without having to cleave the GST and should prove valuable, both to biochemists seeking to check the conformation of their proteins prior to functional studies and to structural biologists screening protein constructs for suitability as targets for structural studies.     

Dynamics of the aromatic amino acid residues in the globular conformation of the basic pancreatic trypsin inhibitor (BPTI)

Summary The molecular conformation of the basic pancreatic trypsin inhibitor (BPTI) is known in considerable detail from both X-ray studies in single crystals and NMR studies in solution. The NMR experiments showed that the aromatic rings of the phenylalanyl and tyrosyl residues can undergo rapid rotational motions about the C-C bond. The present paper describes a model investigation of the mechanistic aspects of these intramolecular rotational motions. From calculations of the conformational energies for molecular species derived from the X-ray structure by rotations of individual aromatic rings, it was apparent that the rotational motions of the aromatics could only be understood in a flexible structure. Flexibility was simulated by allowing the protein to relax to an energetically favorable conformation for each of the different rotation states of the aromatic rings. It was then of particular interest to investigate how the perturbations caused by different rotation states of the aromatic rings were propagated in the protein structure. It was found that the rotation axes C-C were only slightly affected ( ¹20°). The most sizeable perturbations are caused by through space interactions with nearby atoms, which move away from the ring center and thus release the steric hindrance opposing the rotational motions. The values for the energy barriers obtained from the energy minimization are of the same order of magnitude as those measured by NMR.     

Detection of homologous proteins by an intermediate sequence search

We developed a variant of the intermediate sequence search method (ISS(new)) for detection and alignment of weakly similar pairs of protein sequences. ISS(new) relates two query sequences by an intermediate sequence that is potentially homologous to both queries. The improvement was achieved by a more robust overlap score for a match between the queries through an intermediate. The approach was benchmarked on a data set of 2369 sequences of known structure with insignificant sequence similarity to each other (BLAST E-value larger than 0.001); 2050 of these sequences had a related structure in the set. ISS(new) performed significantly better than both PSI-BLAST and a previously described intermediate sequence search method. PSI-BLAST could not detect correct homologs for 1619 of the 2369 sequences. In contrast, ISS(new) assigned a correct homolog as the top hit for 121 of these 1619 sequences, while incorrectly assigning homologs for only nine targets; it did not assign homologs for the remainder of the sequences. By estimate, ISS(new) may be able to assign the folds of domains in approximately 29,000 of the approximately 500,000 sequences unassigned by PSI-BLAST, with 90% specificity (1 - false positives fraction). In addition, we show that the 15 alignments with the most significant BLAST E-values include the nearly best alignments constructed by ISS(new).     

The importance of anchorage in determining a strained protein loop conformation.

We examine the role of the conformational restriction imposed by constrained ends of a protein loop on the determination of a strained loop conformation. The Lys 116-Pro 117 peptide bond of staphylococcal nuclease A exists in equilibrium between the cis and trans isomers. The folded protein favors the strained cis isomer with an occupancy of 90%. This peptide bond is contained in a solvent-exposed, flexible loop of residues 112-117 whose ends are anchored by Val 111 and Asn 118. Asn 118 is constrained by 2 side-chain hydrogen bonds. We investigate the importance of this constraint by replacing Asn 118 with aspartate, alanine, and glycine. We found that removing 1 or more of the hydrogen bonds observed in Asn 118 stabilizes the trans configuration over the cis configuration. By protonating the Asp 118 side chain of N118D through decreased pH, the hydrogen bonding character of Asp 118 approached that of Asn 118 in nuclease A, and the cis configuration was stabilized relative to the trans configuration. These data suggest that the rigid anchoring of the loop end is important in establishing the strained cis conformation. The segment of residues 112-117 in nuclease A provides a promising model system for study of the basic principles that determine polypeptide conformations. Such studies could be useful in the rational design or redesign of protein molecules.     

Resonance electroconformational coupling: A proposed mechanism for energy and signal transductions by membrane proteins

Recent experiments show that membrane ATPases are capable of absorbing free energy from an applied oscillating electric field and converting it to chemical bond energy of ATP or chemical potential energy of concentration gradients. Presumably these enzymes would also respond to endogenous transmembrane electric fields of similar intensity and waveform. A mechanism is proposed in which energy coupling is achieved via Coulombic interaction of an electric field and the conformational equilibria of an ATPase. Analysis indicates that only an oscillating or fluctuating electric field can be used by an enzyme to drive a chemical reaction away from equilibrium.In vivo, the stationary transmembrane potential of a cell must be modulated to become locally oscillatory if it is to derive energy and signal transduction processes.     

Factors determining the reconstructive denaturation of proteins in sodium dodecyl sulfate solutions. Further circular dichroism studies on structural reorganization of seven proteins

The factors determining the onset and extent of reconstructive denaturation of proteins were considered by comparing circular dichroism (CD) data of seven proteins and previously published findings. The effects of sodium dodecyl sulfate (SDS) on the conformation of the following proteins were tested: lysozyme, the mitogens fromPhytolacca americana (fractions Pa2 and Pa4), lectin fromWistaria floribunda, ovine lutropin, a Bence Jones protein, and histone H2B. While the helix content of lysozyme was raised by SDS slightly, in the Bence Jones protein andW. floribunda lectin it increased from near zero to about 25–30%. In histone H2B the helix content was raised by SDS even to about 48%. However, no clear indication of helix formation could be observed in the mitogens and lutropin, even at low pH or 2.0–2.5. The tertiary structure of the proteins was perturbed by SDS. It was concluded that the reorganization of secondary structure of the proteins was favored by the following factors: (1) presence of helicogenic amino acid sequences in the protein, (2) availability of positively charged sites of the basic amino acids for interactions with the dodecyl ion, (3) absence of a large surplus of negatively charged sites on the surface of protein, and (4) absence of extensive disulfide cross-linking within the macromolecule. Both hydrophobic and electrostatic interactions occur in reconstructive denaturation, and the newly formed helices are stabilized by hydrophobic shielding by the alkyl chains of the alkyl sulfate.     

Selecting near-native conformations in homology modeling: the role of molecular mechanics and solvation terms

A free energy function, combining molecular mechanics energy with empirical solvation and entropic terms, is used for ranking near-native conformations that occur in the conformational search steps of homology modeling, i.e., side-chain search and loop closure calculations. Correlations between the free energy and RMS deviation from the X-ray structure are established. It is shown that generally both molecular mechanics and solvation/entropic terms should be included in the potential. The identification of near-native backbone conformations is accomplished primarily by the molecular mechanics term that becomes the dominant contribution to the free energy if the backbone is even slightly strained, as frequently occurs in loop closure calculations. Both terms become equally important if a sufficiently accurate backbone conformation is found. Finally, the selection of the best side-chain positions for a fixed backbone is almost completely governed by the solvation term. The discriminatory power of the combined potential is demonstrated by evaluating the free energies of protein models submitted to the first meeting on Critical Assessment of techniques for protein Structure Prediction (CASP1), and comparing them to the free energies of the native conformations.     

Analycys: a database for conservation and conformation of disulphide bonds in homologous protein domains

Disulphide bonds in proteins are known to play diverse roles ranging from folding to structure to function. Thorough knowledge of the conservation status and structural state of the disulphide bonds will help in understanding of the differences in homologous proteins. Here we present a database for the analysis of conservation and conformation of disulphide bonds in SCOP structural families. This database has a wide range of applications including mapping of disulphide bond mutation patterns, identification of disulphide bonds important for folding and stabilization, modeling of protein tertiary structures and in protein engineering. The database can be accessed at: http://bioinformatics.univ-reunion.fr/analycys/.     

Structural effects of point mutations in proteins

A structural database of 11 families of chains differing by a single amino acid substitution has been built. Another structural dataset of 5 families with identical sequences has been used for comparison. The RMSD computed after a global superimposition of the mutated protein on each native one is smaller than the RMSD calculated among proteins of identical sequences. The effect of the perturbation is very local, and not necessarily the highest at the position of the mutation. A RMSD between mutated and native proteins is computed over a 3‐residue or a 7‐residue window at each position. To separate the effects of structural fluctuations due to point mutations from other sources, pair RMSD have been translated into P values which themselves are included in a score called P‐RANK. This score allows highlighting small backbone distortions by comparing these RMSD between mutated and native positions to the RMSD at the same positions in the absence of a mutation. It results from the P‐RANK that 38% of all mutations produce a significant effect on the displacement. When compared with a random distribution of RMSD at un‐mutated positions, we show that, even if the RMSD is greater when the mutation is in loops than in regular secondary structure, the relative effect is more important for regular secondary structures and for buried positions. We confirm the absence of correlation between RMSD and the predicted variation of free energy of folding but we found a small correlation between high RMSD and the error in the prediction of ΔΔG.     

The bioactive conformation of glucose-dependent insulinotropic polypeptide by NMR and CD spectroscopy

Glucose-dependent insulinotropic polypeptide (GIP) is a gastrointestinal incretin hormone, which modulates physiological insulin secretion. Because of its glucose-sensitive insulinotropic activity, there has been a considerable interest in utilizing the hormone as a potential treatment for type 2 diabetes. Structural parameters obtained from NMR spectroscopy combined with molecular modeling techniques play a vital role in the design of new therapeutic drugs. Therefore, to understand the structural requirements for the biological activity of GIP, the solution structure of GIP was investigated by circular dichroism (CD) followed by proton nuclear magnetic resonance (NMR) spectroscopy. CD studies showed an increase in the helical character of the peptide with increasing concentration of trifluoroethanol (TFE) up to 50%. Therefore, the solution structure of GIP in 50% TFE was determined. It was found that there was an alpha-helix between residues 6 and 29, which tends to extend further up to residue 36. The implications of the C-terminal extended helical segment in the inhibitory properties of GIP on gastric acid secretion are discussed. It is shown that the adoption by GIP of an alpha-helical secondary structure is a requirement for its biological activity. Knowledge of the solution structure of GIP will help in the understanding of how the peptide interacts with its receptor and aids in the design of new therapeutic agents useful for the treatment of diabetes.     

Closed conformation of the active site loop of rabbit muscle triosephosphate isomerase in the absence of substrate: evidence of conformational heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.     

Conservation of polyproline II helices in homologous proteins: implications for structure prediction by model building.

Left-handed polyproline II (PPII) helices commonly occur in globular proteins in segments of 4-8 residues. This paper analyzes the structural conservation of PPII-helices in 3 protein families: serine proteinases, aspartic proteinases, and immunoglobulin constant domains. Calculations of the number of conserved segments based on structural alignment of homologous molecules yielded similar results for the PPII-helices, the alpha-helices, and the beta-strands. The PPII-helices are consistently conserved at the level of 100-80% in the proteins with sequence identity above 20% and RMS deviation of structure alignments below 3.0 A. The most structurally important PPII segments are conserved below this level of sequence identity. These results suggest that the PPII-helices, in addition to the other 2 secondary structure classes, should be identified as part of structurally conserved regions in proteins. This is supported by similar values for the local RMS deviations of the aligned segments for the structural classes of PPII-helices, alpha-helices, and beta-strands. The PPII-helices are shown to participate in supersecondary elements such as PPII-helix/alpha-helix. The conservation of PPII-helices depends on the conservation of a supersecondary element as a whole. PPII-helices also form links, possibly flexible, in the interdomain regions. The role of the PPII-helices in model building by homology is 2-fold; they serve as additional conserved elements in the structure allowing improvement of the accuracy of a model and provide correct chain geometry for modeling of the segments equivalenced to them in a target sequence. The improvement in model building is demonstrated in 2 test studies.