Search

Intercellular gap junction (GJ) contacts formed by the coupling of connexin (Cx) hemichannels (HCs) embedded into the plasma membranes of neighboring cells play significant role in the development, signaling and malfunctions of mammalian tissues. Understanding and targeting GJ functions, however, calls for finding valid Cx subtype-specific inhibitors. We conjecture the lack of information about binding interactions between the GJ interface forming extracellular EL1 and EL2 loops and peptide mimetics designed to specifically inhibit Cx43HC coupling to Cx43GJ. Here, we explore active spots at the GJ interface using known peptide inhibitors that mimic various segments of EL1 and EL2. Binding interactions of these peptide inhibitors and the non-peptide inhibitor quinine has been modelled in combination with the use of blind docking molecular mechanics (MM). The neuron-specific Cx36HC and astrocyte-specific Cx43HC subtypes were modelled with a template derived from the high-resolution structure of Cx26GJ. GJ-coupled and free Cx36HC and Cx43HC models were obtained by dissection of GJs (GJ-coupled) followed by 50 ns molecular dynamics (free). Molecular mechanics (MM) calculations were performed by the docking of inhibitors, explicitly the designed Cx43 EL1 or EL2 loop sequence mimetics (GAP26, P5 or P180–195, GAP27, Peptide5, respectively) and the Cx36 subtype-specific quinine into the model structures. In order to explore specific binding interactions between inhibitors and CxHC subtypes, MM/Generalized Born Surface Area (MM/GBSA) ΔGbind values for representative conformers of peptide mimetics and quinine were evaluated by mapping the binding surface of Cx36HC and Cx43HC for all inhibitors. Quinine specifically contacts Cx36 EL1 residues V54-C55-N56-T57-L58, P60 and N63. Blocking the vestibule by the side of Cx36HC entry, quinine explicitly interacts with the non-conserved V54, L58, N63 residues of Cx36 EL1. In addition, our work challenges the predicted specificity of peptide mimetics, showing that the docking site of peptides is unrelated to the location of the sequence they mimic. Binding features, such as unaffected EL2 residues and the lack of Cx43 subtype-specificity of peptide mimetics, suggest critical roles for peptide stringency and dimension, possibly pertaining to the Cx subtype-specificity of peptide inhibitors.

Recent decades have brought significant changes to the protein structure research field [...]

Mutual synergistic folding (MSF) proteins belong to a recently emerged subclass of disordered proteins, which are disordered in their monomeric forms but become ordered in their oligomeric forms. They can be identified by experimental methods following their unfolding, which happens in a single-step cooperative process, without the presence of stable monomeric intermediates. Only a limited number of experimentally validated MSF proteins are accessible. The amino acid composition of MSF proteins shows high similarity to globular ordered proteins, rather than to disordered ones. However, they have some special structural features, which makes it possible to distinguish them from globular proteins. Even in the possession of their oligomeric three-dimensional structure, classification can only be performed based on unfolding experiments, which are frequently absent. In this work, we demonstrate a simple protocol using molecular dynamics simulations, which is able to indicate that a protein structure belongs to the MSF subclass. The presumption of the known atomic resolution quaternary structure is an obvious limitation of the method, and because of its high computational time requirements, it is not suitable for screening large databases; still, it is a valuable in silico tool for identification of MSF proteins.

Organic anion-transporting polypeptide 3A1 (OATP3A1) is a membrane transporter mediating the cellular uptake of various hormones such as estrone-3-sulfate, prostaglandins E1 and E2 and thyroxine. OATP3A1 is widely expressed in the human body and its presence in tissue-blood barriers, neurons and muscle cells marks it as a potential pharmacological target. Herein we demonstrate that an otherwise membrane impermeant, zwitterionic fluorescent coumarin probe, bearing a sulfonate function is a potent substrate of human OATP3A1, thus readily transported into HEK-293-OATP3A1 cells allowing functional investigation and the screen of drug interactions of the OATP3A1 transporter. At the same time, dyes lacking either the sulfonate motif or the coumarin scaffold showed a dramatic decrease in affinity or even a complete loss of transport. Furthermore, we observed a distinct inhibition/activation pattern in the OATP3A1-mediated uptake of closely related fluorescent coumarin derivatives differing only in the presence of the sulfonate moiety. Additionally, we detected a synergistic effect between one of the probes tested and the endogenous OATP substrate estrone-3-sulfate. These data, together with docking results indicate the presence of at least two cooperative substrate binding sites in OATP3A1. Besides providing the first sensitive probe for testing OATP3A1 substrate/inhibitor interactions, our results also help to understand substrate recognition and transport mechanism of the poorly characterized OATP3A1. Moreover, coumarins are good candidates for OATP3A1-targeted drug delivery and as pharmacological modulators of OATP3A1.

Several intrinsically disordered proteins (IDPs) are capable to adopt stable structures without interacting with a folded partner. When the folding of all interacting partners happens at the same time, coupled with the interaction in a synergistic manner, the process is called Mutual Synergistic Folding (MSF). These complexes represent a discrete subset of IDPs. Recently, we collected information on their complexes and created the MFIB (Mutual Folding Induced by Binding) database. In a previous study, we compared homodimeric MSF complexes with homodimeric and monomeric globular proteins with similar amino acid sequence lengths. We concluded that MSF homodimers, compared to globular homodimeric proteins, have a greater solvent accessible main-chain surface area on the contact surface of the subunits, which becomes buried during dimerization. The main driving force of the folding is the mutual shielding of the water-accessible backbones, but the formation of further intermolecular interactions can also be relevant. In this paper, we will report analyses of heterodimeric MSF complexes. Our results indicate that the amino acid composition of the heterodimeric MSF monomer subunits slightly diverges from globular monomer proteins, while after dimerization, the amino acid composition of the overall MSF complexes becomes more similar to overall amino acid compositions of globular complexes. We found that inter-subunit interactions are strengthened, and additionally to the shielding of the solvent accessible backbone, other factors might play an important role in the stabilization of the heterodimeric structures, likewise energy gain resulting from the interaction of the two subunits with different amino acid compositions. We suggest that the shielding of the &beta, sheet backbones and the formation of a buried structural core along with the general strengthening of inter-subunit interactions together could be the driving forces of MSF protein structural ordering upon dimerization.

Intrinsically Unstructured/Disordered Proteins (IUPs/IDPs) exist as highly flexible conformational ensembles without adopting a stable three-dimensional structure. Experimental and bioinformatical studies in the past two decades have shown that these proteins play a central role in various signaling and regulatory processes. Accordingly, their frequency in higher eukaryotes reaches high proportions and their malfunction can be connected to a wide variety of diseases. Recognizing the biological importance of these proteins motivated researchers to understand various aspects of disordered proteins and protein segments from the viewpoint of biochemistry, molecular biology and pharmacology. In general, IDPs are difficult to study experimentally because of the lack of a unique structure in their isolated form. Nevertheless, taking advantage of ongoing efforts in the collection, cataloguing, and annotation of known IDPs in publicly available databases, various bioinformatics tools were developed over the last few years. These methods enable the further identification and characterization of IDPs using only the amino acid sequence. In this chapter—after a brief introduction to IDPs in general—we present a small survey of current methods aimed at identifying disordered proteins or protein segments, focusing on those that are publicly available as web servers. We also discuss in more detail approaches that predict disordered regions and specific regions involved in protein binding by modeling the physical background of protein disorder. Furthermore, we argue that the heterogeneity of disordered segments needs to be taken into account for a better understanding of protein disorder and the correct use and interpretation of the output of disorder prediction algorithms.

Intrinsically disordered proteins (IDPs) lack a well-defined 3D structure. Their disordered nature enables them to interact with several other proteins and to fulfil their vital biological roles, in most cases after coupled folding and binding. In this paper, we analyze IDPs involved in a new mechanism, mutual synergistic folding (MSF). These proteins define a new subset of IDPs. Recently we collected information on these complexes and created the Mutual Folding Induced by Binding (MFIB) database. These protein complexes exhibit considerable structural variation, and almost half of them are homodimers, but there is a significant amount of heterodimers and various kinds of oligomers. In order to understand the basic background of the disordered character of the monomers found in MSF complexes, the simplest part of the MFIB database, the homodimers are analyzed here. We conclude that MFIB homodimeric proteins have a larger solvent-accessible main-chain surface area on the contact surface of the subunits, when compared to globular homodimeric proteins. The main driving force of the dimerization is the mutual shielding of the water-accessible backbones and the formation of extra intermolecular interactions.

Enterprise Resource Planning Systems provide various core business functions in Financials, Human Resources and Logistics. Oracle and SAP, popular and successful ERP vendors worldwide, continuously improve their products. Standard Query Languages play important roles either during standard or customer developments. Because the execution of Oracle e-Business Suite and SAP Business Suite applications are still based on their relational databases, queries and reports are mainly based on SQL codes designed according to old style, but reliable and robust approaches. It is reasonable to explore how new database technologies and SQL options make changes in design methodologies consequently. This paper discusses some new features of innovative SQL development in SAP, reflecting the effectiveness of innovative options in development. It is found out that developments should designed in a proper way, when old style and new methodologies are balanced for the business functions and also for the company where the developments are introduced.

Implementation of Enterprise Resource Planning Systems usually requires long preparation, brainstorming, design, including detailed re-engineering business processes at large companies. After new business functions have been implemented at enterprise level, it takes some time when the end-users are getting familiar with the concepts, the technology and the daily usage of them. Built-in functions, therefore, make the system reliable, robust and well-known environment. Nevertheless, according to customer and business requirements, it is sometimes necessary to change the scope of business functionalities, the user interfaces. Innovation can have a great impact on business operation at company level. SAP Screen Personas demonstrates how an innovative technology makes easier and comfortable to operate business functions in a simplified user interface. It is also discussed that the role of such new toolset should be considered during a post-implementation ERP process.

The definition of stabilization centers was introduced almost two decades ago. They are centers of noncovalent long range interaction clusters, believed to have a role in maintaining the three-dimensional structure of proteins by preventing their decay due to their cooperative long range interactions. Here, this hypothesis is investigated from the viewpoint of thermal stability for the first time, using a large protein thermodynamics database. The positions of amino acids belonging to stabilization centers are correlated with available experimental thermodynamic data on protein thermal stability. Our analysis suggests that stabilization centers, especially solvent exposed ones, do contribute to the thermal stabilization of proteins.

Therapeutics designed to increase synaptic neurotransmitter levels by inhibiting neurotransmitter sodium symporters (NSSs) classify a strategic approach to treat brain disorders such as depression or epilepsy, however, the critical elementary steps that couple downhill flux of sodium to uphill transport of neurotransmitter are not distinguished as yet. Here we present modelling of NSS member neuronal GAT1 with the substrate#947;-aminobutyric acid (GABA), the major inhibitory neurotransmitter. GABA binding is simulated with the occluded conformation of GAT1 homodimer in an explicit lipid/water environment. Simulations performed in the 1-10 ns range of time elucidated persistent formation of halfextended minor and H-bridged major GABA conformations, referred to as binding and traverse conformations, respectively. The traverse GABA conformation was further stabilized by GAT1-bound Na(+)(1). We also observed Na(+)(1) translocation to GAT1-bound Cl(-) as well as the appearance of water molecules at GABA and GAT1-bound Na(+)(2), conjecturing causality. Scaling dynamics suggest that the traverse GABA conformation may be valid for developing substrate inhibitors with high efficacy. The potential for this finding is significant with impact not only in pharmacology but wherever understanding of the mechanism of neurotransmitter uptake is valuable.

A glutaminyl cyclase (QC) fragment library was in silico selected by disconnection of the structure of known QC inhibitors and by lead-like 2D virtual screening of the same set. The resulting fragment library (204 compounds) was acquired from commercial suppliers and pre-screened by differential scanning fluorimetry followed by functional in vitro assays. In this way, 10 fragment hits were identified ( $$\sim $$ 5 % hit rate, best inhibitory activity: 16 $$\upmu \hbox {M}$$ ). The in vitro hits were then docked to the active site of QC, and the best scoring compounds were analyzed for binding interactions. Two fragments bound to different regions in a complementary manner, and thus, linking those fragments offered a rational strategy to generate novel QC inhibitors. Based on the structure of the virtual linked fragment, a 77-membered QC target focused library was selected from vendor databases and docked to the active site of QC. A PubChem search confirmed that the best scoring analogues are novel, potential QC inhibitors.

Understanding of multidrug binding at the atomic level would facilitate drug design and strategies to modulate drug metabolism, including drug transport, oxidation, and conjugation. Therefore we explored the mechanism of promiscuous binding of small molecules by studying the ligand binding domain, the PAS-B domain of the aryl hydrocarbon receptor (AhR). Because of the low sequence identities of PAS domains to be used for homology modeling, structural features of the widely employed HIF-2α and a more recent suitable template, CLOCK were compared. These structures were used to build AhR PAS-B homology models. We performed molecular dynamics simulations to characterize dynamic properties of the PAS-B domain and the generated conformational ensembles were employed in in silico docking. In order to understand structural and ligand binding features we compared the stability and dynamics of the promiscuous AhR PAS-B to other PAS domains exhibiting specific interactions or no ligand binding function. Our exhaustive in silico binding studies, in which we dock a wide spectrum of ligand molecules to the conformational ensembles, suggest that ligand specificity and selection may be determined not only by the PAS-B domain itself, but also by other parts of AhR and its protein interacting partners. We propose that ligand binding pocket and access channels leading to the pocket play equally important roles in discrimination of endogenous molecules and xenobiotics.

There is a delicate balance between stability and flexibility needed for enzyme function. To avoid undesirable alteration of the functional properties during the evolutionary optimization of the structural stability under certain circumstances, and vice versa, to avoid unwanted changes of stability during the optimization of the functional properties of proteins, common sense would suggest that parts of the protein structure responsible for stability and parts responsible for function developed and evolved separately. This study shows that nature did not follow this anthropomorphic logic: the set of residues involved in function and those involved in structural stabilization of enzymes are rather overlapping than segregated.

c ¨¨ ´ ´ ´ ´ ´ ´ Abstract Some structural features underlying the increased thermostability of enzymes from thermophilic organisms relative to their homologues from mesophiles are known from earlier studies.We used cellulase C from Clostridium thermocellum to test whether thermostability can be increased by mutations designed using rules learned from thermophilic proteins.Cellulase C has a TIM barrel fold with an additional helical subdomain.We designed and produced a number of mutants with the aim to increase its thermostability.Five mutants were designed to create new electrostatic interactions.They all retained catalytic activity but exhibited decreased thermostability relative to the wild-type enzyme.Here, the stabilizing contributions are obviously smaller than the destabilization caused by the introduction of the new side chains.In another mutant, the small helical subdomain was deleted.This mutant lost activity but its melting point was only 3 8C lower than that of the wild-type enzyme, which suggests that the subdomain is an independent folding unit and is important for catalytic function.A double mutant was designed to introduce a new disulfide bridge into the enzyme.This mutant is active and has an increased stability (DT s3 8C, m D(DG )s1.73 kcalymol) relative to the wild-type enzyme.Reduction of the disulfide bridge results in destabilization u and an altered thermal denaturation behavior.We conclude that rules learned from thermophilic proteins cannot be used in a straightforward way to increase the thermostability of a protein.Creating a crosslink such as a disulfide bond is a relatively sure-fire method but the stabilization may be smaller than calculated due to coupled destabilizing effects. 2002 Elsevier Science B.V. All rights reserved.

The well-balanced stability of protein structures allows large-scale fluctuations, which are indispensable in many biochemical functions, ensures the long-term persistence of the equilibrium structure and it regulates the degradation of proteins to provide amino acids for biosynthesis. This balance is studied in the present work with two sets of proteins by analyzing stabilization centers, defined as certain clusters of residues involved in cooperative long-range interactions. One data set contains 56 proteins, which belong to 16 families of homologous proteins, derived from organisms of various physiological temperatures. The other set is composed of 31 major histocompatibility complex (MHC)–peptide complexes, which represent peptide transporters complexed with peptide ligands that apparently contribute to the stabilization of the MHC proteins themselves. We show here that stabilization centers, which had been identified as special clusters of residues that protect the protein structure, evolved to serve also as regulators of function – related degradation of useless protein as part of protein housekeeping.

The comparison of the three-dimensional structures of thermophilic (Thermus thermophilus) and mesophilic (Escherichia coli) 3-isopropylmalate dehydrogenases (IPMDH, EC 1.1.1.85) suggested that the existence of extra ion pairs in the thermophilic enzyme found in the intersubunit region may be an important factor for thermostability. As a test of our assumption, glutamine 200 in the E. coli enzyme was turned into glutamate (Q200E mutant) to mimic the thermophilic enzyme at this site by creating an intersubunit ion pair which can join existing ion clusters. At the same site in the thermophilic enzyme we changed glutamate 190 into glutamine (E190Q), hereby removing the corresponding ion pair. These single amino acid replacements resulted in increased thermostability of the mesophilic and decreased thermostability of the thermophilic enzyme, as measured by spectropolarimetry and differential scanning microcalorimetry.

Intrinsically unstructured/disordered proteins (IUPs/IDPs) exist as highly flexible conformational ensembles without adopting a stable three-dimensional structure. Experimental and bioinformatical studies in the past two decades have shown that these proteins play a central role in various signaling and regulatory processes. Accordingly, their frequency in higher eukaryotes reaches high proportions and their malfunction can be connected to a wide variety of diseases. Recognizing the biological importance of these proteins motivated researchers to understand various aspects of disordered proteins and protein segments from the viewpoint of biochemistry, molecular biology and pharmacology. In general, IDPs are difficult to study experimentally because of the lack of a unique structure in the isolated form. Nevertheless, various bioinformatics tools developed over the last few years enable their identification and characterization using only the amino acid sequence. In this chapter — after a brief introduction to IDPs in general — we present a small survey of current methods aimed at identifying disordered proteins or protein segments, focusing on those that are publicly available as web servers. We also discuss in more detail approaches that predict disordered regions and specific regions involved in protein binding by modeling the physical background of protein disorder. Furthermore, we argue that the heterogeneity of disordered segments needs to be taken into account for a better understanding of protein disorder and the correct use and interpretation of the output of disorder prediction algorithms.

The leuB gene from the psychrotrophic strain Vibrio sp. I5 has been cloned and sequenced. The gene codes for 3-isopropylmalate dehydrogenase, a 360-residue, dimeric enzyme involved in the biosynthesis of leucine. Three recently solved homologous isopropylmalate dehydrogenase (IPMDH) crystal structures from thermophilic and mesophilic organisms have been used to build a homology model for the psychrotrophic IPMDH and to deduce the possible structural reasons for its decreased thermostability. According to our model the psychrotrophic IPMDH contains fewer stabilizing interactions than its mesophilic and thermophilic counterparts. Elements that have been identified as destabilizing in the comparison of the psychrotrophic, mesophilic and thermophilic IPMDHs are a smaller number of salt-bridges, a reduction in aromatic-aromatic interactions, fewer proline residues and longer surface loops. In addition, there are a number of substitutions of otherwise strictly conserved residues that can be linked to thermostability.

To reveal the structural basis of the increased thermal stability of 3-isopropylmalate dehydrogenase (IPMDH) from Thermus thermophilus, an extreme thermophile, the homology-based structural model of one mesophilic (Escherichia coli) counterpart, was constructed. Both IPMDHs are homodimeric proteins. We built a model of one subunit using the 3-D structures of the Th. thermophilus IPMDH and the homologous E.coli isocitrate dehydrogenase. Energy minimization and molecular dynamics simulated annealing were performed on the dimer, including a surrounding solvation shell. No serious errors were detected in the refined model using the 3-D profile method. The resulting structure was scrutinized and compared with the structure of the Th.thermophilus IPMDH. Significant differences were found in the non-specific interactions including the hydrophobic effect. The model predicts a higher number of ion pairs in the Th.thermophilus than in the E.coli enzyme. An increase was observed in the stabilities of alpha-helical regions in the thermophilic protein. The preliminary X-ray coordinates of the E.coli IPMDH were received after the completion of this work, allowing an assessment of the model in terms of the X-ray structure. The comparison proved that most of the structural features underlying the stability differences between the two enzymes were predicted correctly.

Phosphorylation at Ser-133 of the kinase inducible domain of CREB (KID) triggers its binding to the KIX domain of CBP via a concomitant coil-to-helix transition. The exact role of this key event is still puzzling: it does not switch between disordered and ordered states, nor its direct interactions fully account for selectivity. Hence, we reasoned that phosphorylation may shift the conformational preferences of KID towards a binding-competent state. To this end we investigated the intrinsic structural properties of the unbound KID in phosphorylated and unphosphorylated forms by simulated annealing and molecular dynamics simulations. Although helical populations show subtle differences, phosphorylation reduces the flexibility of the turn segment connecting the two helices in the complexed structure and induces a transient structural element that corresponds to its bound conformation. It is stabilized by the pSer-133–Arg-131 interaction, which is absent from the unphosphorylated KID. Diminishing this coupling decreases the 3.1 kcal/mol contribution of pSer-133 to the binding free energy (ΔGbind) of the phosphorylated KID to KIX by 1.1 kcal/mol, as computed in reference to Ser-133. In a binding competent form of the S133E KID mutant, the contribution of Glu-133 to ΔGbind is by 1.5 kcal/mol smaller than that of pSer, suggesting that altered structural properties due to pSer → Glu replacement impair the binding affinity. Thus, we propose that phoshorylation contributes to selectivity not merely by the direct interactions of the phosphate group with KIX, but also by promoting the formation of a transient structural element in the highly conserved turn segment. Proteins 2006. © 2006 Wiley-Liss, Inc.

The flexibility of prolyl oligopeptidase has been investigated using molecular dynamics (MD) and molecular framework approaches to delineate the route of the substrate to the active site. The selectivity of the enzyme is mediated by a seven-bladed beta-propeller that in the crystal structure does not indicate the possible passage for the substrate to the catalytic center. Its open topology however, could allow the blades to move apart and let the substrate into the large central cavity. Flexibility analysis of prolyl oligopeptidase structure using the FIRST (Floppy Inclusion and Rigid Substructure Topology) approach and the atomic fluctuations derived from MD simulations demonstrated the rigidity of the propeller domain, which does not permit the substrate to approach the active site through this domain. Instead, a smaller tunnel at the inter-domain region comprising the highly flexible N-terminal segment of the peptidase domain and a facing hydrophilic loop from the propeller (residues 192-205) was identified by cross-correlation analysis and essential dynamics as the only potential pathway for the substrate. The functional importance of the flexible loop has been also verified by kinetic analysis of the enzyme with a split loop. Catalytic effect of engineered disulfide bridges was rationalized by characterizing the concerted motions of the two domains.

In nature, 1 out of every 10 proteins has an (alpha/beta)(8) (TIM)-barrel fold, and in most cases, pairwise comparisons show no sequence similarity between them. Hence, delineating the key residues that induce very different sequences to share a common fold is important for understanding the folding and stability of TIM-barrel domains. In this work, we propose a new consensus approach for locating these stabilizing residues based on long-range interactions, hydrophobicity, and conservation of amino acid residues. We have identified 957 stabilizing residues in 63 proteins from a nonredundant set of 71 TIM-barrel domains. Most of these residues are located in the 8-stranded beta-sheet, with nearly one half of them oriented toward the interior of the barrel and the other half oriented toward the surrounding alpha-helices. Several stabilizing residues are found in the N- and C-terminal loops, whereas very few appear in the alpha-helices that surround the internal beta-sheet. Further, these 957 residues are placed in 434 stabilizing segments of various sizes, and each domain contains 1-10 of these segments. We found that 8 segments per domain is the most abundant one, and two thirds of the proteins have 7-9 stabilizing segments. Finally, we verified the identified residues with experimental temperature factors and found that these residues are among the ones with less mobility in the considered proteins. We suggest that our new protocol serves as a powerful tool to identify the stabilizing residues in TIM-barrel domains, which can be used as potential candidates for studying protein folding and stability by means of protein engineering experiments.

Summary: SCide is a program to identify stabilization centers from known protein structures. These are residues involved in cooperative long-range contacts, which can be formed between various regions of a single polypeptide chain, or they can belong to different peptides or polypeptides in a complex. The server takes a PDB file as an input, and the result is presented in graphical or text format. Availability: SCide is available on the web at http://www. enzim.hu/scide. The source code can be obtained from the authors on request. Contact: zsuzsa@enzim.hu The analysis of contacts present in protein structures can help to elucidate the difficult relationship between the 3D structure and the sequence encoding it. The three-dimensional structures are strongly related to the thermodynamic properties of proteins (Gromiha et al., 2002). Mutational studies as well as analysis of conservation patterns in proteins suggest that only a small portion of the residues are indispensable for the protein for defining its fold, and ensuring its folding and stability (Reddy et al., 2001). Yet, identifying the relevant positions is a challenging problem even with the knowledge of the structure. Earlier analyses identified some important contacts by locating dense residue clusters, mainly of hydrophobic character (Zehfus, 1995; Kannan and Vishveshwara, 1999). Different aspects of structural requirements can be investigated by focusing on short or long range contacts, since folding is dominated by short-range interactions, while long-range contacts are more likely to play an important role in ensuring stability of the structure and preventing its unfolding (Doszt ´ anyi et al., 1997; Baker, 2000). The SCide program is intended to select stabilization centers (SCs), i.e. amino acids involved in contacts which are important in maintaining the stability of a protein. Two residues are defined as forming a long range contact if they are at least 10 amino acids apart in the sequence, and ∗ To whom correspondence should be addressed. they have at least one heavy atom separation smaller than the sum of their van der Waals radii plus 1 ˚ A. From these contacts, two residues are defined as SC elements when it is possible to select one‐one residues from both flanking tetrapeptides of the two residues such that three residues from each segment make at least seven contacts out of the possible nine. For more details about stabilization centers and their properties in proteins the reader is referred

Calculations of the IR and Raman spectra of molecules containing linear carbon segments are reported. H-C 2 n -H ( n = 1, 2, 3, 4, 5) and H 2 -C 2 n -H 2 ( n = 1, 2, 3, 4, 5) were choosen to model the alternating and the cumulenic structures, respectively. Geometry optimization and force constant calculation were done on the ab initio Hartree-Fock level using Dunning's double-ζ basis set with polarization. The frequencies were scaled by Pulay's method. IR- and Raman-intensities were obtained by calculating dipole- and polarizability derivatives, respectively.