Application of a directed conformational search for generating 3-D coordinates for protein structures from alpha-carbon coordinates.

A directed conformational search algorithm using the program CONGEN (ref. 3), which samples backbone conformers, is described. The search technique uses information from the partially built structures to direct the search process and is tested on the problem of generating a full set of backbone Cartesian coordinates given only alpha-carbon coordinates. The method has been tested on six proteins of known structure, varying in size and classification, and was able to generate the original backbone coordinates with RMSs ranging from 0.30-0.87A for the alpha-carbons and 0.5-0.99A RMSs for the backbone atoms. Cis peptide linkages were also correctly identified. The procedure was also applied to two proteins available with only alpha-carbon coordinates in the Brookhaven Protein Data Bank; thioredoxin (SRX) and triacylglycerol acylhydrolase (TGL). All-atom models are proposed for the backbone of both these proteins. In addition, the technique was applied to randomized coordinates of flavodoxin to assess the effects of irregularities in the data on the final RMS. This study represents the first time a deterministic conformational search was used on such a large scale.     

A molecular dynamics approach for the generation of complete protein structures from limited coordinate data

Generation of full protein coordinates from limited information, e.g., the Cα coordinates, is an important step in protein homology modeling and structure determination, and molecular dynamics (MD) simulations may prove to be important in this task. We describe a new method, in which the protein backbone is built quickly in a rather crude way and then refined by minimization techniques. Subsequently, the side chains are positioned using extensive MD calculations. The method is tested on two proteins, and results compared to proteins constructed using two other MD-based methods. In the first method, we supplemented an existing backbone building method with a new procedure to add side chains. The second one largely consists of available methodology. The constructed proteins are compared to the corresponding X-ray structures, which became available during this study, and they are in good agreement (backbone RMS values of 0.5–0.7 Å, and all-atom RMS values of 1.5–1.9 Å). This comparative study indicates that extensive MD simulations are able, to some extent, to generate details of the native protein structure, and may contribute to the development of a standardized methodology to predict reliably (parts of) protein structures when only partial coordinate data are available. © 1994 John Wiley & Sons, Inc.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

Stereochemical quality of protein structure coordinates.

Methods have been developed to assess the stereochemical quality of any protein structure both globally and locally using various criteria. Several parameters can be derived from the coordinates of a given structure. Global parameters include the distribution of phi, psi and chi 1 torsion angles, and hydrogen bond energies. There are clear correlations between these parameters and resolution; as the resolution improves, the distribution of the parameters becomes more clustered. These features show a broad distribution about ideal values derived from high-resolution structures. Some structures have tightly clustered distributions even at relatively low resolutions, while others show abnormal scatter though the data go to high resolution. Additional indicators of local irregularity include proline phi angles, peptide bond planarities, disulfide bond lengths, and their chi 3 torsion angles. These stereochemical parameters have been used to generate measures of stereochemical quality which provide a simple guide as to the reliability of a structure, in addition to the most important measures, resolution and R-factor. The parameters used in this evaluation are not novel, and are easily calculated from structure coordinates. A program suite is currently being developed which will quickly check a given structure, highlighting unusual stereochemistry and possible errors.     

Designed molecules that fold to mimic protein secondary structures

Molecules that fold to mimic protein secondary structures have emerged as important targets of bioorganic chemistry. Recently, a variety of compounds that mimic helices, turns, and sheets have been developed, with notable advances in the design of beta-peptides that mimic each of these structures. These compounds hold promise as a step toward synthetic molecules with protein-like properties and as drugs that block protein-protein interactions.     

Three-dimensional model and quaternary structure of the human eye lens protein gamma S-crystallin based on beta- and gamma-crystallin X-ray coordinates and ultracentrifugation.

A 3-dimensional model of the human eye lens protein gamma S-crystallin has been constructed using comparative modeling approaches encoded in the program COMPOSER on the basis of the 3-dimensional structure of gamma-crystallin and beta-crystallin. The model is biased toward the monomeric gamma B-crystallin, which is more similar in sequence. Bovine gamma S-crystallin was shown to be monomeric by analytical ultracentrifugation without any tendency to form assemblies up to concentrations in the millimolar range. The connecting peptide between domains was therefore built assuming an intramolecular association as in the monomeric gamma-crystallins. Because the linker has 1 extra residue compared with gamma B and beta B2, the conformation of the connecting peptide was constructed by using a fragment from a protein database. gamma S-crystallin differs from gamma B-crystallin mainly in the interface region between domains. The charged residues are generally paired, although in a different way from both beta- and gamma-crystallins, and may contribute to the different roles of these proteins in the lens.     

The importance of surface loops for stabilizing an eightfold beta alpha barrel protein.

An important step in understanding how a protein folds is to determine those regions of the sequence that are critical to both its stability and its folding pathway. We chose phosphoribosyl anthranilate isomerase from Escherichia coli, which is a monomeric representative of the (beta alpha)8 barrel family of proteins, to construct a variant that carries an internal tandem duplication of the fifth beta alpha module. This (beta alpha)9 variant was enzymically active and therefore must have a wild-type (beta alpha)8 core. It had a choice a priori to fold to three different folding frames, which are distinguished by carrying the duplicated segment as an insert into one out of three different loops. Steady-state kinetic constants, the fluorescence properties of a crucial tryptophan residue, and limited proteolysis showed that the stable (beta alpha)9 variant carries the insertion between beta-strand 5 and alpha-helix 5. This preference can be explained by the important role of loops between alpha helices and beta strands in stabilizing the structure of the enzyme.     

An optimization approach to predicting protein structural class from amino acid composition.

Proteins are generally classified into four structural classes: all-alpha proteins, all-beta proteins, alpha + beta proteins, and alpha/beta proteins. In this article, a protein is expressed as a vector of 20-dimensional space, in which its 20 components are defined by the composition of its 20 amino acids. Based on this, a new method, the so-called maximum component coefficient method, is proposed for predicting the structural class of a protein according to its amino acid composition. In comparison with the existing methods, the new method yields a higher general accuracy of prediction. Especially for the all-alpha proteins, the rate of correct prediction obtained by the new method is much higher than that by any of the existing methods. For instance, for the 19 all-alpha proteins investigated previously by P.Y. Chou, the rate of correct prediction by means of his method was 84.2%, but the correct rate when predicted with the new method would be 100%! Furthermore, the new method is characterized by an explicable physical picture. This is reflected by the process in which the vector representing a protein to be predicted is decomposed into four component vectors, each of which corresponds to one of the norms of the four protein structural classes.     

Calculation of protein backbone geometry from alpha-carbon coordinates based on peptide-group dipole alignment.

An algorithm is proposed for the conversion of a virtual-bond polypeptide chain (connected C alpha atoms) to an all-atom backbone, based on determining the most extensive hydrogen-bond network between the peptide groups of the backbone, while maintaining all of the backbone atoms in energetically feasible conformations. Hydrogen bonding is represented by aligning the peptide-group dipoles. These peptide groups are not contiguous in the amino acid sequence. The first dipoles to be aligned are those that are both sufficiently close in space to be arranged in approximately linear arrays termed dipole paths. The criteria used in the construction of dipole paths are: to assure good alignment of the greatest possible number of dipoles that are close in space; to optimize the electrostatic interactions between the dipoles that belong to different paths close in space; and to avoid locally unfavorable amino acid residue conformations. The equations for dipole alignment are solved separately for each path, and then the remaining single dipoles are aligned optimally with the electrostatic field from the dipoles that belong to the dipole-path network. A least-squares minimizer is used to keep the geometry of the alpha-carbon trace of the resulting backbone close to that of the input virtual-bond chain. This procedure is sufficient to convert the virtual-bond chain to a real chain; in applications to real systems, however, the final structure is obtained by minimizing the total ECEPP/2 (empirical conformational energy program for peptides) energy of the system, starting from the geometry resulting from the solution of the alignment equations. When applied to model alpha-helical and beta-sheet structures, the algorithm, followed by the ECEPP/2 energy minimization, resulted in an energy and backbone geometry characteristic of these alpha-helical and beta-sheet structures. Application to the alpha-carbon trace of the backbone of the crystallographic 5PTI structure of bovine pancreatic trypsin inhibitor, followed by ECEPP/2 energy minimization with C alpha-distance constraints, led to a structure with almost as low energy and root mean square deviation as the ECEPP/2 geometry analog of 5PTI, the best agreement between the crystal and reconstructed backbone being observed for the residues involved in the dipole-path network.     

Elucidating protein secondary structures using alpha-carbon recurrence quantifications.

Secondary structures of proteins were studied by recurrence quantification analysis (RQA). High-resolution, 3-dimensional coordinates of alpha-carbon atoms comprising a set of 68 proteins were downloaded from the Protein Data Bank. By fine-tuning four recurrence parameters (radius, line, residue, separation), it was possible to establish excellent agreement between percent contribution of alpha-helix and beta-sheet structures determined independently by RQA and that of the DSSP algorithm (Define Secondary Structure of Proteins). These results indicate that there is an equivalency between these two techniques, which are based upon totally different pattern recognition strategies. RQA enhances qualitative contact maps by quantifying the arrangements of recurrent points of alpha carbons close in 3-dimensional space. For example, the radius was systematically increased, moving the analysis beyond local alpha-carbon neighborhoods in order to capture super-secondary and tertiary structures. However, differences between proteins could only be detected within distances up to about 6-11 A, but not higher. This result underscores the complexity of alpha-carbon spacing when super-secondary structures appear at larger distances. Finally, RQA-defined secondary structures were found to be robust against random displacement of alpha carbons upwards of 1 A. This finding has potential import for the dynamic functions of proteins in motion.     

A helix-turn-strand structural motif common in alpha-beta proteins.

By exhaustive structural comparisons, we have found that about one-third of the alpha-helix-turn-beta-strand polypeptides in alpha-beta barrel domains share a common structural motif. The chief characteristics of this motif are that first, the geometry of the turn between the alpha-helix and the beta-strand is somewhat constrained, and second, the beta-strand contains a hydrophobic patch that fits into a hydrophobic pocket on the alpha-helix. The geometry of the turn does not seem to be a major determinant of the alpha-beta unit, because the turns vary in length from four to six residues. However, the motif does not occur when there are few constraints on the geometry of the turn-for instance, when the turns between the alpha-helix and the beta-strands are very long. It also occurs much less frequently in flat-sheet alpha-beta proteins, where the topology is much less regular and the amount of twist on the sheet varies considerably more than in the barrel proteins. The motif may be one of the basic building blocks from which alpha-beta barrels are constructed.     

Thermal isoelectric precipitation of alpha-lactalbumin from a whey protein concentrate: Influence of protein-calcium complexation

The selective precipitation of alpha-lactalbumin (alpha-LA) at a pH around its isoelectric point (4.2) under heat treatment is the basis for a fractionation process of whey proteins. As precipitation is a phenomenon dependent on the protein hydrophobicity, and as the release of the tightly bound calcium occurring at pH around 4 modifies the alpha-LA hydrophobicity, the specific role of calcium on isoelectric precipitation is investigated. A study of the extent of alpha-LA precipitation in a whey protein concentrate under various operating conditions of pH, temperature, protein concentration, and calcium content is presented. We propose a mechanism for this phenomenon as a combination of a complexation equilibrium and of an irreversible precipitation, to account for the influence of temperature, alpha-LA concentration total ionic content, and calcium concentration, and also to estimate the complexation equilibrium constant. (c) 1995 John Wiley & Sons, Inc.     

Properties of integral membrane protein structures: derivation of an implicit membrane potential

Distributions of each amino acid in the trans-membrane domain were calculated as a function of the membrane normal using all currently available alpha-helical membrane protein structures with resolutions better than 4 A. The results were compared with previous sequence- and structure-based analyses. Calculation of the average hydrophobicity along the membrane normal demonstrated that the protein surface in the membrane domain is in fact much more hydrophobic than the protein core. While hydrophobic residues dominate the membrane domain, the interfacial regions of membrane proteins were found to be abundant in the small residues glycine, alanine, and serine, consistent with previous studies on membrane protein packing. Charged residues displayed nonsymmetric distributions with a preference for the intracellular interface. This effect was more prominent for Arg and Lys resulting in a direct confirmation of the positive inside rule. Potentials of mean force along the membrane normal were derived for each amino acid by fitting Gaussian functions to the residue distributions. The individual potentials agree well with experimental and theoretical considerations. The resulting implicit membrane potential was tested on various membrane proteins as well as single trans-membrane alpha-helices. All membrane proteins were found to be at an energy minimum when correctly inserted into the membrane. For alpha-helices both interfacial (i.e. surface bound) and inserted configurations were found to correspond to energy minima. The results demonstrate that the use of trans-membrane amino acid distributions to derive an implicit membrane representation yields meaningful residue potentials.     

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on Cα atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 Å, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low. © 1995 Wiley-Liss, Inc.     

Structural analysis of the 2.8 A model of Xylose isomerase from Actinoplanes missouriensis

The structure of Xylose isomerase (X.I.) from Actinoplanes missouriensis has been solved to 2.8 Angstroms resolution. Phases were determined from a single Eu3+ derivative and from the noncrystallographic 222 symmetry of the tetrameric molecule. An atomic model was built and subjected to restrained crystallographic refinement. The resulting model is shown to be closely similar to the recently reported X.I.'s structures from three other bacterial sources. Each monomer is found to be composed of an eight-stranded alpha/beta "T.I.M." barrel forming an N-terminal domain of 328 residues followed by a large loop of 66 residues embracing an adjacent subunit. Analysis of intersubunit packing shows that the X.I. tetramer is an assembly of two tight dimers. The beta barrel fits a simple hyperboloid model as other T.I.M. barrels do. The active site, identified as the binding site for the inhibitor xylitol, is located at the carboxyl end of the beta strands in the barrel next to a pair of binding sites for Eu3+ ions, which are assumed to be sites for the divalent ions involved in catalysis. Active sites in the tetramer are oriented towards the interface between dimers. It is suggested that subunit interfaces might stabilize the active site region and this might explain the oligomeric nature of other alpha/beta barrel enzymes.     

Nicotine-evoked transmitter release from synaptosomes: functional association of specific presynaptic acetylcholine receptors and voltage-gated calcium channels.

It has previously been shown that nicotine-evoked dopamine release from rat striatal synaptosomes and nicotine-evoked norepinephrine release from hippocampal synaptosomes are mediated by distinct nicotinic acetylcholine receptor (nAChR) subtypes. In the present study, the functional association of these nicotinic receptors with specific subtypes of voltage-gated calcium channels was examined. Cd(2+) (200 microM), as well as omega-conotoxin MVIIC (5 microM), blocks approximately 85% of nicotine-evoked dopamine release from striatal synaptosomes, indicating a major involvement of calcium channels. Furthermore, the toxin-susceptibility suggests that these calcium channels contain alpha(1A) and/or alpha(1B) subunits. Inhibition of nicotine-evoked dopamine release by conotoxins alpha-MII and omega-GVIA is additive and indicates that presynaptic alpha3beta2 nAChRs are functionally coupled to alpha(1A), but not alpha(1B), calcium channel subtypes. Conversely, insensitivity to alpha-AuIB and sensitivity to omega-MVIIC indicate that non-alpha3beta2/alpha3beta4-containing nAChRs are functionally coupled to alpha(1B)-containing calcium channels. In contrast, Cd(2+) blocks only 65% of nicotine-evoked norepinephrine release from hippocampal synaptosomes, indicating that a substantial fraction of this release occurs through mechanisms not involving calcium channels. This Cd(2+)-insensitive component of release is blocked by alpha-AuIB and therefore appears to be triggered by Ca(2+) flowing directly through the channels of presynaptic alpha3beta4 nAChRs. Thus, these data indicate that different presynaptic termini can have distinctive functional associations of specific nAChRs and voltage-gated calcium channels.