Structure Of Globular Proteins Research Articles

Evaluation of the structural quality of modeled proteins by using globularity criteria

BackgroundThe knowledge of the three-dimensional structure of globular proteins is fundamental for a detailed investigation of their functional properties. Experimental methods are too slow for structure investigation on a large scale, while computational prediction methods offer alternatives that are continuously being improved. The international Comparative Assessment of Structure Prediction (CASP), an "a posteriori" evaluation of the quality of theoretical models when the experimental structure becomes available, demonstrates that predictions can be successful as well as unsuccessful, and this suggests the necessity for evaluations able to discard "a priori" the wrong models.ResultsWe analyzed different structural properties of globular proteins for experimentally solved proteins belonging to the four different structural classes: "mainly alpha", "mainly beta", "alpha/beta" and "alpha+beta". The properties were found to be linearly correlated to protein molecular weight, but with some differences among the four classes. These results were applied to develop an evaluation test of theoretical models based on the expected globular properties of proteins. To verify the success of our test, we applied it to several protein models submitted to the sixth edition of CASP. The best theoretical models, as judged by CASP assessors, were in agreement with the expected properties, while most of the low-quality models had not passed our evaluations.ConclusionThis study supports the need for careful checks to avoid the diffusion of incorrect structural models. Our test allows the evaluation of models in the absence of experimental reference structures, thereby preventing the diffusion of incorrect structural models and the formulation of incorrect functional hypotheses. It can be used to check the globularity of predicted models, and to supplement other methods already used to evaluate their quality.

Molecular Lipophilicity in Protein Modeling and Drug Design

Hydrophobic interactions play a key role in the folding and maintenance of the 3-dimensional structure of proteins, as well as in the binding of ligands (e.g. drugs) to protein targets. Therefore, quantitative assessment of spatial hydrophobic (lipophilic) properties of these molecules is indispensable for the development of efficient computational methods in drug design. One possible solution to the problem lies in application of a concept of the 3-dimensional molecular hydrophobicity potential (MHP). The formalism of MHP utilizes a set of atomic physicochemical parameters evaluated from octanol-water partition coefficients (log P) of numerous chemical compounds. It permits detailed assessment of the hydrophobic and/or hydrophilic properties of various parts of molecules and may be useful in analysis of protein-protein and protein-ligand interactions. This review surveys recent applications of MHP-based techniques to a number of biologically relevant tasks. Among them are: (i) Detailed assessment of hydrophobic/hydrophilic organization of proteins; (ii) Application of this data to the modeling of structure, dynamics, and function of globular and membrane proteins, membrane-active peptides, etc. (iii) Employment of the MHP-based criteria in docking simulations for ligands binding to receptors. It is demonstrated that the application of the MHP-based techniques in combination with other molecular modeling tools (e.g. Monte Carlo and molecular dynamics simulations, docking, etc.) permits significant improvement to the standard computational approaches, provides additional important insights into the intimate molecular mechanisms driving protein assembling in water and in biological membranes, and helps in the computer-aided drug discovery process.

Prediction of Local Structural Stabilities of Proteins from Their Amino Acid Sequences

Hydrogen exchange experiments provide detailed information about the local stability and the solvent accessibility of different regions of the structures of folded proteins, protein complexes, and amyloid fibrils. We introduce an approach to predict protection factors from hydrogen exchange in proteins based on the knowledge of their amino acid sequences without the inclusion of any additional structural information. These results suggest that the propensity of different regions of the structures of globular proteins to undergo local unfolding events can be predicted from their amino acid sequences with an accuracy of 80% or better.

Stability of globular proteins in H2O and D2O

In several experimental techniques D2O rather then H2O is often used as a solvent for proteins. Concerning the influence of the solvent on the stability of the proteins, contradicting results have been reported in literature. In this paper the influence of H2O-D2O solvent substitution on the stability of globular protein structure is determined in a systematic way. The differential scanning calorimetry technique is applied to allow for a thermodynamic analysis of two types of globular proteins: hen's egg lysozyme (LSZ) with relatively strong internal cohesion ("hard" globular protein) and bovine serum albumin (BSA), which is known for its conformational adaptability ("soft" globular protein). Both proteins tend to be more stable in D2O compared to H2O. We explain the increase of protein stability in D2O by the observation that D2O is a poorer solvent for nonpolar amino acids than H2O, implying that the hydrophobic effect is larger in D2O. In case of BSA the transitions between different isomeric forms, at low pH values the Nm and F forms, and at higher pH values Nm and B, were observed by the presence of a supplementary peak in the DSC thermogram. It appears that the pH-range for which the Nm form is the preferred one is wider in D2O than in H2O.

Bhageerath: an energy based web enabled computer software suite for limiting the search space of tertiary structures of small globular proteins

We describe here an energy based computer software suite for narrowing down the search space of tertiary structures of small globular proteins. The protocol comprises eight different computational modules that form an automated pipeline. It combines physics based potentials with biophysical filters to arrive at 10 plausible candidate structures starting from sequence and secondary structure information. The methodology has been validated here on 50 small globular proteins consisting of 2–3 helices and strands with known tertiary structures. For each of these proteins, a structure within 3–6 Å RMSD (root mean square deviation) of the native has been obtained in the 10 lowest energy structures. The protocol has been web enabled and is accessible at .

Protein Modules Conserved Since LUCA

Universal scale of the sequence conservation has been recently introduced based on omnipresence of the protein sequence motifs across species. A large spectrum of short sequences, up to eight residues has been found to reside in all or almost all prokaryotic organisms. By this discovery a principally novel quantitative approach is introduced to the problem of reconstruction of the last universal common ancestor (LUCA). The most conserved elements (protein modules) with defined structures and sequences harboring the omnipresent motifs are outlined in this work, by combining the sequence and protein crystal structure data. The structurally conserved modules involve 25-30 amino acid residues and have appearance of closed loops, loop-n-lock structures. This confirms earlier conclusions on the loop-fold structure of globular proteins. Many of the topmost conserved modules represent the primary closed loop prototypes, that have been derived by whole genome sequence searches. The data presented, thus, make a basis for further developments toward the earliest stages of protein evolution.

Inferring topological features of proteins from amino acid residue networks

Topological properties of native folds are obtained from statistical analysis of 160 low homology proteins covering the four structural classes. This is done analyzing one, two and three-vertex joint distribution of quantities related to the corresponding network of amino acid residues. Emphasis on the amino acid residue hydrophobicity leads to the definition of their center of mass as vertices in this contact network model with interactions represented by edges. The network analysis helps us to interpret experimental results such as hydrophobic scales and fraction of buried accessible surface area in terms of the network connectivity. Moreover, those networks show assortative mixing by degree. To explore the vertex-type dependent correlations, we build a network of hydrophobic and polar vertices. This procedure presents the wiring diagram of the topological structure of globular proteins leading to the following attachment probabilities between hydrophobic–hydrophobic 0.424(5), hydrophobic-polar 0.419(2) and polar–polar 0.157(3) residues.

Prion Protein Helix1 Promotes Aggregation but Is Not Converted into β-Sheet

Prion diseases are caused by the aggregation of the native alpha-helical prion protein PrP(C) into its pathological beta-sheet-rich isoform PrP(Sc). In current models of PrP(Sc), helix1 is assumed to be preferentially converted into beta-sheet during aggregation of PrP(C). This was supported by the NMR structure of PrP(C) since, in contrast to the isolated helix1, helix2 and helix3 are connected by a small loop and are additionally stabilized by an interhelical disulfide bond. However, helix1 is extremely hydrophilic and has a high helix propensity. This prompted us to investigate the role of helix1 in prion aggregation using humPrP(23-159) including helix1 (144-156) compared with the C-terminal-truncated isoform humPrP(23-144) corresponding to the pathological human stop mutations Q160Stop and Y145Stop, respectively. Most unexpectedly, humPrP(23-159) aggregated significantly faster compared with the truncated fragment humPrP(23-144), clearly demonstrating that helix1 is involved in the aggregation process. However, helix1 is not resistant to digestion with proteinase K in fibrillar humPrP(23-159), suggesting that helix1 is not converted to beta-sheet. This is confirmed by Fourier transformation infrared spectroscopy since there is almost no difference in beta-sheet content of humPrP(23-159) fibrils compared with humPrP(23-144). In conclusion, we provide strong direct evidence that in contrast to earlier assumptions helix1 is not converted into beta-sheet during aggregation of PrP(C) to PrP(Sc).

Glasslike Structure of Globular Proteins and the Boson Peak

Vibrational spectra of proteins and topologically disordered solids display a common anomaly at low frequencies, known as boson peak. We show that such feature in globular proteins can be deciphered in terms of an energy landscape picture, as it is for glassy systems. Exploiting the tools of Euclidean random matrix theory, we clarify the physical origin of such anomaly in terms of a mechanical instability of the system. As a natural explanation, we argue that such instability is relevant for proteins in order for their molecular functions to be optimally rooted in their structures.

Luminescent spectral properties of rhodamine derivatives while binding to serum albumin

We have studied the effect of blood serum albumin on the absorption and fluorescence spectra of rhodamine C (RC), rhodamine 6G (R6G), and rhodamine 3B (R3B). Interaction of the dye with protein is assessed using the binding parameters: binding constants and concentrations of binding sites. We have studied the effect of temperature on the binding parameters. We have observed that heating a mixture of the dye solution with protein for 30 min leads to an increase in the binding constant for rhodamine 3B with protein by a factor of 2, while the concentration of binding sites increases by a factor of 2.3. This is explained by features of the globular protein structure and a change in its conformation when heated. We have shown that rhodamine 3B at a concentration of 10−5 M is the most effective among the studied rhodamine dyes for application as a fluorescent probe when studying conformational changes in blood serum protein.

Branched-Chain Amino Acids: Enzyme and Substrate Regulation1–3

The three branched-chain amino acids (BCAAs) are the most hydrophobic of the amino acids and play crucial roles in determining the structures of globular proteins as well as the interaction of the transmembrane domains of membranous proteins with phospholipid bilayers. However, the three BCAAs do not behave identically. In terms of protein secondary structure, valine and isoleucine exhibit a definite preference for the beta-structure, whereas leucine has a higher preference for the alpha-helix. Although mutation of one BCAA to another is commonly regarded as conservative, there are well-documented examples of such substitutions that have a significant effect on protein function. The occurrence of BCAA in nature is, therefore, attributable to their primary role in protein structure, not to their secondary metabolic roles. These functions are important for almost all proteins; therefore, BCAA commonly account for approximately 20-25% of most dietary proteins. Dietary BCAA largely escape first-pass splanchnic metabolism. The first steps in their catabolism are common to all three, involving the BCAA aminotransferase (BCAT) and branched-chain alpha-keto acid dehydrogenase (BCKD). Their further metabolism employs distinct pathways to different end-products (glucose and/or ketone bodies). However, the fact that the flux-generating step for the catabolism of the three BCAAs occurs at one of the common steps indicates that the production of these downstream products are not individually regulated and, hence, may not play important individual roles. The catabolism of the BCAAs is highly regulated by both allosteric and covalent mechanisms. BCKD is inhibited by phosphorylation and activated by dephosphorylation. Allosteric inhibition of the kinase by the branched-chain keto acids (BCKA) (particularly by alpha-ketoisocaproate) serves both as a mechanism for promoting the catabolism of excess quantities of these amino acids as well as for conserving low concentrations of these dietary essential amino acids. Cytosolic and mitochondrial isoenzymes of BCAT have been identified. They are thought to play an important role in brain neurotransmitter metabolism.

Secondary Structure of Globular Proteins in Adsorption Layers at the Solution-Air Interface by the Data of Fourier Transform IR Spectroscopy

The secondary structures of globular proteins of different structural types (α, (α + β), and β) are studied in the adsorption layers at the solution-air interface using Fourier transform infrared spectroscopy. Similar investigations are performed for the initial proteins in a powdered state; these data are used as a control. It is shown that, at the water-air interface, either the secondary structure of the studied proteins changes insignificantly or remains intact depending on the structural type of a protein. The adsorption of proteins of the α-helical type (bovine and human serum albumins) and the (α + β) type (lysozyme) virtually do not result in changes in their secondary structure. A protein of β-structural type, α-chymotrypsin, shows a lower stability at the water-air interface.

Sequence swapping does not result in conformation swapping for the β4/β5 and β8/β9 β‐hairpin turns in human acidic fibroblast growth factor

The beta-turn is the most common type of nonrepetitive structure in globular proteins, comprising ~25% of all residues; however, a detailed understanding of effects of specific residues upon beta-turn stability and conformation is lacking. Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold and contains a total of five beta-hairpin structures (antiparallel beta-sheets connected by a reverse turn). beta-Turns related by the characteristic threefold structural symmetry of this superfold exhibit different primary structures, and in some cases, different secondary structures. As such, they represent a useful system with which to study the role that turn sequences play in determining structure, stability, and folding of the protein. Two turns related by the threefold structural symmetry, the beta4/beta5 and beta8/beta9 turns, were subjected to both sequence-swapping and poly-glycine substitution mutations, and the effects upon stability, folding, and structure were investigated. In the wild-type protein these turns are of identical length, but exhibit different conformations. These conformations were observed to be retained during sequence-swapping and glycine substitution mutagenesis. The results indicate that the beta-turn structure at these positions is not determined by the turn sequence. Structural analysis suggests that residues flanking the turn are a primary structural determinant of the conformation within the turn.

A computational pathway for bracketing native-like structures for small alpha helical globular proteins

Impressive advances in the applications of bioinformatics for protein structure prediction coupled with growing structural databases on one hand and the insurmountable time-scale problem with ab initio computational methods on the other continue to raise doubts whether a computational solution to the protein folding problem--categorized as an NP-hard problem--is within reach in the near future. Combining some specially designed biophysical filters and vector algebra tools with ab initio methods, we present here a promising computational pathway for bracketing native-like structures of small alpha helical globular proteins departing from secondary structural information. The automated protocol is initiated by generating multiple structures around the loops between secondary structural elements. A set of knowledge-based biophysical filters namely persistence length and radius of gyration, developed and calibrated on approximately 1000 globular proteins, is introduced to screen the trial structures to filter out improbable candidates for the native and reduce the size of the library of probable structures. The ensemble so generated encompasses a few structures with native-like topology. Monte Carlo optimizations of the loop dihedrals are then carried out to remove steric clashes. The resultant structures are energy minimized and ranked according to a scoring function tested previously on a series of decoy sets vis-a-vis their corresponding natives. We find that the 100 lowest energy structures culled from the ensemble of energy optimized trial structures comprise at least a few to within 3-5 angstroms of the native. Thus the formidable "needle in a haystack" problem is narrowed down to finding an optimal solution amongst a computationally tractable number of alternatives. Encouraging results obtained on twelve small alpha helical globular proteins with the above outlined pathway are presented and discussed.

Universal positions in globular proteins

The description of globular protein structures as an ensemble of contiguous 'closed loops' or 'tightened end fragments' reveals fold elements crucial for the formation of stable structures and for navigating the very process of protein folding. These are the ends of the loops, which are spatially close to each other but are situated apart in the polypeptide chain by 25-30 residues. They also correlate with the locations of highly conserved hydrophobic residues (referred to as topohydrophobic), in a structural alignment of the members of a protein family. This study analysed these positions in 111 representatives of different protein folds, and then carried out dynamic Monte Carlo simulations of the first steps of the folding process, aimed at predicting the origins of the assembling folds. The simulations demonstrated that there is an obvious trend for certain sets of residues, named 'mostly interacting residues', to be buried at the early stages of the folding process. Location of these residues at the loop ends and correlation with topohydrophobic positions are demonstrated, thereby giving a route to simulations of the protein folding process.

Applications of time-resolved resonance energy transfer measurements in studies of the molecular crowding effect.

The native structures of many globular proteins are only weakly stabilized and form in solution ensembles of multiple conformers. The energy differences between the conformers are assumed to be small. This is the case of flexible multidomain proteins where domain motions were observed. High concentrations of inert macrosolute, which create a crowded or confined environment, can cause shifts of the distribution of the conformers of such proteins towards the more compact structures. This effect may also promote compact structures in partially folded proteins. Time-resolved dynamic non-radiative excitation energy transfer (tr-RET) is suitable for detection of either subtle or major changes in distributions of intramolecular distances in protein molecules in solutions. Two experiments were performed which demonstrated the applicability of tr-RET for detection of the effect of macrosolutes on the conformational ensembles of flexible states of protein molecules. The distribution of distances between residues 203 and 169 in the CORE domain of E. coli adenylate kinase (AK) in the denatured state was determined in the presence of high concentrations of dextran 40. A significant shift of the mean of the distribution was observed without reduction of its width. This was interpreted as a shift to compact structure without change of the degree of disorder of the chain. In a second experiment the distribution of the distance between residues 55 and 169 in AK, which spans the cleft between the CORE and the AMPbind domains, was monitored. No clear effect of high concentrations of dextran 40 was found. These experiments show the strength of the application of tr-RET in investigation of changes in the sub-states of flexible conformations of globular protein. Networks of pairs of labeled sites can be prepared and tr-RET experiments can be performed in order to search for the segments of the protein molecules, which respond to the presence of inert macromolecules in their environment.

Effect of pH on the Stability and Structure of Yeast Hexokinase A: ACIDIC AMINO ACID RESIDUES IN THE CLEFT REGION ARE CRITICAL FOR THE OPENING AND THE CLOSING OF THE STRUCTURE

pH and salts have a marked effect on the stability, structure, and function of many globular proteins due to their ability to influence the electrostatic interactions. In this work, calorimetry, CD, and fluorescence studies have been carried out to understand the pH-dependent conformational changes of the two-domain protein yeast hexokinase A. In conjunction with the crystal structural data available, the present results have enabled the complete characterization and analysis of the pH-dependent conformational changes of the enzyme that have strong implications in understanding its structure-function relationship. The calorimetric profiles show a single thermal transition in the acidic pH range, whereas two independent transitions were observed in the alkaline pH range, suggesting the structural merger of the domains at the acidic pH. Comparison of the thermal transitions at pH 8.5 studied by different techniques suggests that the first transition corresponds to the smaller domain, and the second transition corresponds to the larger domain. The acid-denatured state of hexokinase A has high secondary structure content with little or no tertiary interactions and binds to the hydrophobic dye 8-anilinonaphthalene-1-sulfonic acid, suggesting that it is a molten globule-like state, whereas the alkali-denatured state is less structured than the acid-denatured state but more structured than the urea-denatured state, suggestive of a premolten globule-like state. Structural analysis using the published hexokinase B structure as well as the hexokinase A structure with the revised amino acid sequence in conjunction with the results obtained by us suggests that the ionization state of the acidic residues at the active site could regulate domain movements that are responsible for the opening and the closure of the cleft between the two domains and in turn affect the structure and function of the enzyme.

The occurrence of C--H...O hydrogen bonds in alpha-helices and helix termini in globular proteins.

A comprehensive database analysis of C--H...O hydrogen bonds in 3124 alpha-helices and their corresponding helix termini has been carried out from a nonredundant data set of high-resolution globular protein structures resolved at better than 2.0 A in order to investigate their role in the helix, the important protein secondary structural element. The possible occurrence of 5 --> 1 C--H...O hydrogen bond between the ith residue CH group and (i - 4)th residue C==O with C...O < or = 3.8 A is studied, considering as potential donors the main-chain Calpha and the side-chain carbon atoms Cbeta, Cgamma, Cdelta and Cepsilon. Similar analysis has been carried out for 4 --> 1 C--H...O hydrogen bonds, since the C--H...O hydrogen bonds found in helices are predominantly of type 5 --> 1 or 4 --> 1. A total of 17,367 (9310 of type 5 --> 1 and 8057 of type 4 --> 1) C--H...O hydrogen bonds are found to satisfy the selected criteria. The average stereochemical parameters for the data set suggest that the observed C--H...O hydrogen bonds are attractive interactions. Our analysis reveals that the Cgamma and Cbeta hydrogen atom(s) are frequently involved in such hydrogen bonds. A marked preference is noticed for aliphatic beta-branched residue Ile to participate in 5 --> 1 C--H...O hydrogen bonds involving methylene Cgamma 1 atom as donor in alpha-helices. This may be an enthalpic compensation for the greater loss of side-chain conformational entropy for beta-branched amino acids due to the constraint on side-chain torsion angle, namely, chi1, when they occur in helices. The preference of amino acids for 4 --> 1 C--H...O hydrogen bonds is found to be more for Asp, Cys, and for aromatic residues Trp, Phe, and His. Interestingly, overall propensity for C--H...O hydrogen bonds shows that a majority of the helix favoring residues such as Met, Glu, Arg, Lys, Leu, and Gln, which also have large side-chains, prefer to be involved in such types of weak attractive interactions in helices. The amino acid side-chains that participate in C--H...O interactions are found to shield the acceptor carbonyl oxygen atom from the solvent. In addition, C--H...O hydrogen bonds are present along with helix stabilizing salt bridges. A novel helix terminating interaction motif, X-Gly with Gly at C(cap) position having 5 --> 1 Calpha--H...O, and a chain reversal structural motif having 1 --> 5 Calpha-H...O have been identified and discussed. Our analysis highlights that a multitude of local C--H...O hydrogen bonds formed by a variety of amino acid side-chains and Calpha hydrogen atoms occur in helices and more so at the helix termini. It may be surmised that the main-chain Calpha and the side-chain CH that participate in C--H...O hydrogen bonds collectively augment the cohesive energy and thereby contribute together with the classical N--H...O hydrogen bonds and other interactions to the overall stability of helix and therefore of proteins.

Reconstruction of protein structures from a vectorial representation.

We show that the contact map of the native structure of globular proteins can be reconstructed starting from the sole knowledge of the contact map's principal eigenvector, and present an exact algorithm for this purpose. Our algorithm yields a unique contact map for all 221 globular structures of PDBselect25 of length N</=120. We also show that the reconstructed contact maps allow in turn for the accurate reconstruction of the three-dimensional structure. These results indicate that the reduced vectorial representation provided by the principal eigenvector of the contact map is equivalent to the protein structure itself. This representation is expected to provide a useful tool in bioinformatics algorithms for protein structure comparison and alignment, as well as a promising intermediate step towards protein structure prediction.

A general clustering approach with application to the Miyazawa-Jernigan potentials for amino acids.

In this article, we address the problem of classification of amino acids. Starting from the Miyazawa-Jernigan matrix obtained from the relative positions of amino acids in the crystal structure of globular proteins, we develop a fully unsupervised method of classification for the amino acids. The method is based in the subdominant ultrametric associated to the distance induced by the Miyazawa-Jernigan matrix and the maximum likelihood principle to determine the cluster structure. We obtain a classification consistent with the five groups used in the literature, although with some peculiarities. We also show the stability of our results against changes of the method used to classify the amino acids. Proteins 2004.