Interaction code for polar and nonpolar amino acids: "ice-breaker" model]

A novel model for the study of recognition and interaction code of amino acids in peptides, proteins and their complexes has been proposed. The model is designed on the modern notions on the structure and properties of water and hydrophobic bonds. It is assumed that the polar side chains of amino acids during the formation of the hydrophobic bonds act as "ice-breaker", thus destroying the organized structure of water (clusters or "icebergs") around the hydrophobic radicals of amino acids.     

Fuzzy oil drop model to interpret the structure of antifreeze proteins and their mutants

Mutations in proteins introduce structural changes and influence biological activity: the specific effects depend on the location of the mutation. The simple method proposed in the present paper is based on a two-step model of in silico protein folding. The structure of the first intermediate is assumed to be determined solely by backbone conformation. The structure of the second one is assumed to be determined by the presence of a hydrophobic center. The comparable structural analysis of the set of mutants is performed to identify the mutant-induced structural changes. The changes of the hydrophobic core organization measured by the divergence entropy allows quantitative comparison estimating the relative structural changes upon mutation. The set of antifreeze proteins, which appeared to represent the hydrophobic core structure accordant with “fuzzy oil drop” model was selected for analysis.     

Studies on protein folding, unfolding and fluctuations by computer simulation. IV. Hydrophobic interactions.

The theoretical model of proteins on the two-dimensional square lattice, introduced previously, is extended to include the hydrophobic interactions. Two proteins, whose native conformations have different folded patterns, are studied. Units in the protein chains are classified into polar units and nonpolar units. If there is a vacant lattice point next to a nonpolar unit, it is interpreted as being occupied by solvent water and the entropy of the system is assumed to decrease by a certain amount. Besides these hydrophobic free energies, the specific long-range interactions studied in previous papers are assumed to be operative in a protein chain. Equilibrium properties of the folding and unfolding transitions of the two proteins are found to be similar, even though one of them was predicted, based on the one globule model of the transitions, to unfold through a significant intermediate state (or at least to show a tendency toward such a behavior), when the hydrophobic interactions are strongly weighted. The failure of this prediction led to the development of a more refined model of transitions; a non-interacting local structure model. The hydrophobic interactions assumed here have a character of non-specific long-range interactions. Because of this character the hydrophobic interactions have the effect of decelerating the folding kinetics. The deceleration effect is less pronounced in one of the two proteins, whose native conformation is stabilized by many pairs of medium-range interactions. It is therefore inferred that the medium-range interactions have the power to cope with the decelerating effect of the non-specific hydrophobic interactions.     

Redesigning symmetry-related "mini-core" regions of FGF-1 to increase primary structure symmetry: thermodynamic and functional consequences of structural symmetry

Previous reports detailing mutational effects within the hydrophobic core of human acidic fibroblast growth factor (FGF-1) have shown that a symmetric primary structure constraint is compatible with a stably folded protein. In the present report, we investigate symmetrically related pairs of buried hydrophobic residues in FGF-1 (termed "mini-cores") that are not part of the central core. The effect upon the stability and function of FGF-1 mutations designed to increase primary structure symmetry within these "mini-core" regions was evaluated. At symmetry-related positions 22, 64, and 108, the wild-type protein contains either Tyr or Phe side chains. The results show that either residue can be readily accommodated at these positions. At symmetry-related positions 42, 83, and 130, the wild-type protein contains either Cys or Ile side chains. While positions 42 and 130 can readily accommodate either Cys or Ile side chains, position 83 is substantially destabilized by substitution by Ile. Tertiary structure asymmetry in the vicinity of position 83 appears responsible for the inability to accommodate an Ile side chain at this position, and is known to contribute to functional half-life. A mutant form of FGF-1 with enforced primary structure symmetry at positions 22, 64, and 108 (all Tyr) and 42, 83, and 130 (all Cys) is shown to be more stable than the reference FGF-1 protein. The results support the hypothesis that a symmetric primary structure within a symmetric protein superfold represents a solution to achieving a foldable, stable polypeptide, and highlight the role that function may play in the evolution of asymmetry within symmetric superfolds.     

Conformation of a purified "spontaneously" inserting thylakoid membrane protein precursor in aqueous solvent and detergent micelles

Subunit W of photosystem II (PsbW) is a single-span thylakoid membrane protein that is synthesized with a cleavable hydrophobic signal peptide and integrated into the thylakoid membrane by an apparently spontaneous mechanism. In this study, we have analyzed the secondary structure of the pre-protein at early stages of the insertion pathway, using purified recombinant pre-PsbW. We show that the protein remains soluble in Tris buffer after removal of detergent. Under these conditions pre-PsbW contains no detectable alpha-helix, whereas substantial alpha-helical structure is present in SDS micelles. In aqueous buffer, the tryptophan fluorescence emission characteristics are intermediate between those of solvent-exposed and hydrophobic environments, suggesting the formation of a partially folded structure. If denaturants are excluded from the purification protocol, pre-PsbW purifies instead as a 180-kDa oligomer with substantial alpha-helical structure. Mature-size PsbW was prepared by removal of the presequence, and we show that this protein also contains alpha-helix in detergent but in lower quantities than the pre-protein. We therefore propose that pre-PsbW contains alpha-helical structure in both the mature protein and the signal peptide in nonpolar environments. We propose that pre-PsbW acquires its alpha-helical structure only during the later, membrane-bound stages of the insertion pathway, after which it forms a "helical hairpin"-type loop intermediate in the thylakoid membrane.     

Effects of side-chain characteristics on stability and oligomerization state of a de novo-designed model coiled-coil: 20 amino acid substitutions in position "d"

We describe the de novo design and biophysical characterization of a model coiled-coil protein in which we have systematically substituted 20 different amino acid residues in the central "d" position. The model protein consists of two identical 38 residue polypeptide chains covalently linked at their N termini via a disulfide bridge. The hydrophobic core contained Val and Ile residues at positions "a" and Leu residues at positions "d". This core allowed for the formation of both two-stranded and three-stranded coiled-coils in benign buffer, depending on the substitution at position "d". The structure of each analog was analyzed by CD spectroscopy and their relative stability determined by chemical denaturation using GdnHCI (all analogs denatured from the two-stranded state). The oligomeric state(s) was determined by high-performance size-exclusion chromatography and sedimentation equilibrium analysis in benign medium. Our results showed a thermodynamic stability order (in order of decreasing stability) of: Leu, Met, Ile, Tyr, Phe, Val, Gln, Ala, Trp, Asn, His, Thr, Lys, Ser, Asp, Glu, Arg, Orn, and Gly. The Pro analog prevented coiled-coil formation. The overall stability range was 7.4 kcal/mol from the lowest to the highest analog, indicating the importance of the hydrophobic core and the dramatic effect a single substitution in the core can have upon the stability of the protein fold. In general, the side-chain contribution to the level of stability correlated with side-chain hydrophobicity. Molecular modelling studies, however, showed that packing effects could explain deviations from a direct correlation. In regards to oligomerization state, eight analogs demonstrated the ability to populate exclusively one oligomerization state in benign buffer (0.1 M KCl, 0.05 M K(2)PO(4)(pH 7)). Ile and Val (the beta-branched residues) induced the three-stranded oligomerization state, whereas Tyr, Lys, Arg, Orn, Glu and Asp induced the two-stranded state. Asn, Gln, Ser, Ala, Gly, Phe, Leu, Met and Trp analogs were indiscriminate and populated two-stranded and three-stranded states. Comparison of these results with similar substitutions in position "a" highlights the positional effects of individual residues in defining the stability and numbers of polypeptide chains occurring in a coiled-coil structure. Overall, these results in conjunction with other work now generate a relative thermodynamic stability scale for 19 naturally occurring amino acid residues in either an "a" or "d" position of a two-stranded coiled-coil. Thus, these results will aid in the de novo design of new coiled-coil structures, a better understanding of their structure/function relationships and the design of algorithms to predict the presence of coiled-coils within native protein sequences.     

Study of the "molten globule" intermediate state in protein folding by a hydrophobic fluorescent probe

Binding of the hydrophobic fluorescent probe, 1-anilino-naphthalene-8-sulfonate (ANS), to synthetic polypeptides and proteins with a different structural organization has been studied. It has been shown that ANS has a much stronger affinity to the protein "molten globule" state, with a pronounced secondary structure and compactness, but without a tightly packed tertiary structure as compared with its affinity to the native and coil-like proteins, or to coil-like, alpha-helical, or beta-structural hydrophilic homopolypeptides. The possibility of using ANS for the study of equilibrium and kinetic molten globule intermediates is demonstrated, with carbonic anhydrase, beta-lactamase, and alpha-lactalbumin as examples.     

"De novo" design of peptides with specific lipid-binding properties

In this study, we describe an in silico method to design peptides that can be made of non-natural amino acids and elicit specific membrane-interacting properties. The originality of the method holds in the capacities developed to design peptides from any non-natural amino acids as easily as from natural ones, and to test the structure stability by an angular dynamics rather than the currently-used molecular dynamics. The goal of this study was to design a non-natural tilted peptide. Tilted peptides are short protein fragments able to destabilize lipid membranes and characterized by an asymmetric distribution of hydrophobic residues along their helix structure axis. The method is based on the random generation of peptides and their selection on three main criteria: mean hydrophobicity and the presence of at least one polar residue; tilted insertion at the level of the acyl chains of lipids of a membrane; and conformational stability in that hydrophobic phase. From 10,000,000 randomly-generated peptides, four met all the criteria. One was synthesized and tested for its lipid-destabilizing properties. Biophysical assays showed that the "de novo" peptide made of non-natural amino acids is helical either in solution or into lipids as tested by Fourier transform infrared spectroscopy and is able to induce liposome fusion. These results are in agreement with the calculations and validate the theoretical approach.     

Crystal structure of the "cab"-type beta class carbonic anhydrase from the archaeon Methanobacterium thermoautotrophicum

The structure of the "cab"-type beta class carbonic anhydrase from the archaeon Methanobacterium thermoautotrophicum (Cab) has been determined to 2.1-A resolution using the multiwavelength anomalous diffraction phasing technique. Cab exists as a dimer with a subunit fold similar to that observed in "plant"-type beta class carbonic anhydrases. The active site zinc is coordinated by protein ligands Cys(32), His(87), and Cys(90), with the tetrahedral coordination completed by a water molecule. The major difference between plant- and cab-type beta class carbonic anhydrases is in the organization of the hydrophobic pocket. The structure reveals a Hepes buffer molecule bound 8 A away from the active site zinc, which suggests a possible proton transfer pathway from the active site to the solvent.     

Interaction of N1-, N6- and C8-substituted derivatives of adenosine-5'-triphosphate with the catalytic subunit of cAMP-dependent protein kinase from rabbit skeletal muscles

In order to investigate the structure of the active site of the cAMP-dependent protein kinase catalytic subunit a synthesis of several previously unknown adenosine-5'-triphosphate (ATP) derivatives containing substituents of various nature at N(1), N(C6) and C(8) positions of the purine base was carried out. The interaction of these derivatives with a homogeneous preparation of the catalytic subunit of rabbit skeletal muscle cAMP-dependent protein kinase was investigated. All the nucleotide analogs were found to inhibit the enzyme activity; the inhibition was competitive with respect to ATP. It was assumed that the adenine moiety of the ATP molecule is bound to the active site of protein kinase by the hydrophobic interaction with the aromatic amino acid residues and by formation of the hydrogen bond between the exo-NH2-group of the substrate and a corresponding group of the enzyme. The "correct" binding of ATP to the enzyme active center is defined by the anti-conformation of the nucleotide.     

Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain–side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein’s structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein–protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at . The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation. Figure Fuzzy-oil-drop model applied to identify the ligation site in lysozyme complexed with N-acetylglucosamine (PDB ID:1LMQ) in form of hydrophobicity deficiency (ΔH) profile and three-dimensional distribution of on protein surface     

Fuzzy-oil-drop hydrophobic force field--a model to represent late-stage folding (in silico) of lysozyme

A model of hydrophobic collapse (in silico), which is generally considered to be the driving force for protein folding, is presented in this work. The model introduces the external field in the form of a fuzzy-oil-drop assumed to represent the environment. The drop is expressed in the form of a three-dimensional Gauss function. The usual probability value is assumed to represent the hydrophobicity distribution in the three-dimensional space of the virtual environment. The differences between this idealized hydrophobicity distribution and the one represented by the folded polypeptide chain is the parameter to be minimized in the structure optimization procedure. The size of fuzzy-oil-drop is critical for the folding process. A strong correlation between protein length and the dimension of the native and early-stage molecular form was found on the basis of single-domain proteins analysis. A previously presented early-stage folding (in silico) model was used to create the starting structure for the procedure of late-stage folding of lysozyme. The results of simulation were found to be promising, although additional improvements for the formation of beta-structure and disulfide bonds as well as the participation of natural ligand in folding process seem to be necessary.     

Collapse transitions in protein-like lattice polymers: The effect of sequence patterns

The collapse transition of lattice protein-like heteropolymers has been studied by means of the Monte Carlo method. The protein model has been reduced to the α-carbon trace restricted to a high coordination lattice. The sequences of model heteropolymers contain two types of mers: hydrophobic/nonpolar (H) and hydrophilic/polar (P). Interactions of HH and PP pairs were assumed to be negative (weaker attractions of PP pairs) while the contact energy for HP pairs was equal to zero. All sequence-specific short-range interactions have been neglected in the present studies. It has been found that homopolymeric chains undergo a smooth collapse transition to a dense globular state. The globule lacks any signatures of local ordering that could be interpreted as a model of protein secondary structure. Hetero-polymers with the sequences of hydrophilic and hydrophobic residues characteristic for α- and β-type proteins undergo a somewhat sharper (though continuous) collapse transition to a dense globular state with elements of local ordering controlled by the sequence. The helical pattern induces more secondary structure than the β-type pattern. For all examined sequences the level of local ordering was lower than the average secondary structure content of globular proteins. The results are compared with other theoretical work and with known experimental facts. The implications for the reduced modeling of protein systems are briefly discussed. © 1997 John Wiley & Sons, Inc. Biopoly 42: 537–548, 1997     

Mechanical Properties of Thermo-moulded Biofilms in Relation to Proteins/Starch Interactions

The effects of hydrophilic and hydrophobic characteristics of proteins on the interactions with corn starch were investigated in this study. The model system included corn starch and proteins, i.e. zein, gliadin, gluten, soy protein and rapeseed protein. The blend films were prepared by thermo-moulding in gentle conditions at 70 °C in order to avoid starch gelatinization, with respect to water content, and avoid protein denaturation. The effects of different kinds of proteins on structure and mechanical behaviour of blend biomaterials were characterised by scanning electron microscopy (SEM) and tensile test, respectively. The effects of different kinds of proteins on intermolecular interactions between proteins and starch were investigated by dynamical mechanical thermal analysis. Based on the solubility measurement results, almost all protein films showed the similar solubility to the natural protein powders, resulting from the weak influence of mild thermo-moulding treatment on protein inner structure. Different morphologies were observed for different proteins and corresponding blends, which are relatively loose protein architecture that appeared for hydrophobic protein and blend films, and uniform and densely packed architecture for hydrophilic ones. Moreover, different mechanical behaviours were obtained for different proteins and corresponding blends. No significantly increased strength for hydrophilic protein blends with starch added can be explained that there is weak intermolecular interaction between both components based on SEM observation. However, the addition of corn starch granules in hydrophobic protein networks was assumed that starch destroyed or weakened the protein network, resulting in the decrease of mechanical strength.     

Folding patterns of porin and bacteriorhodopsin.

Porin spans the outer membrane of Escherichia coli with most of the protein embedded within the membrane. It lacks pronounced hydrophobic domains and consists predominantly of beta-pleated sheet. These observations require the accommodation of polar and ionizable residues in an environment that has a low dielectric constant. Owing to a currently limited understanding of the constraints governing membrane protein structure, a minimal approach to structure prediction is proposed that identifies segments causing polypeptides to reverse their direction (turn identification). The application of this procedure avoids hydrophobicity parameters and yields a model of porin which is in good agreement with all experimental data available. The presence of polar and ionizable residues within membrane boundaries implies a dense (saturating) network of hydrogen bond donor and acceptor groups. Application to a paradigm of hydrophobic membrane proteins, bacteriorhodopsin, reveals a pattern consistent with its alpha-helical folding. The postulated structure includes significantly more polar residues in the membrane domain than have been assumed previously, suggesting that there are also hydrogen bonding networks in bacteriorhodopsin. Extensive networks permeating protein interior and surfaces would explain the extraordinary stability and the tight interactions between functional units in the formation of crystalline arrays of both proteins.     

Hydrophobic Regions in Myelin Proteins Characterized Through Analysis of "Hydropathic"Profiles

Computer-generated "hydropathic" profiles were constructed for graphic comparison of the amino acid sequences for P2 protein, 18.5 kilodalton (kDa) myelin basic protein (BP), and myelin proteolipid protein (PLP). Profiles were also obtained for cytochrome b5, a membrane protein known to be capable of reversible association with lipid bilayers and of a size comparable to that of the myelin BPs. Analysis of the PLP sequence produced profiles generally compatible with the suggestions that PLP has three transbilayer and two bilayer intercalating segments. Profiles for P2 and 18.5 kDa BP were found to contain hydrophilic segments separated by relatively short hydrophobic regions. Whereas hydropathic indices in hydrophobic regions of P2, 18.5 kDa BP, and PLP fall in the value ranges recently reported for cores of globular proteins and intrabilayer domains of membrane proteins, hydrophobic sections of P2 and 18.5 kDa BP have hydropathic indices similar to those in the hydrophobic core (transprotein) regions of globular proteins. None of them are comparable to the region of cytochrome b5 known to anchor that protein in its membrane or to the segments of PLP sequence proposed as intrabilayer domains. This comparison suggests that neither BP has structural characteristics compatible with insertion into the hydrocarbon core of the myelin lipid bilayer, a conclusion that is consistent with a recently published study that identified the bilayer penetrating proteins of myelin with a hydrophobic probe. The above findings suggest an enhancement for some details of myelin architecture and a cautious approach to interpreting data for BP intercalation into bilayers.     

Molecular cloning and homology modelling of a subtilisin-like serine protease from the marine fungus, <Emphasis Type="Italic">Engyodontium album</Emphasis> BTMFS10

An alkaline protease gene (Eap) was isolated for the first time from a marine fungus, Engyodontium album. Eap consists of an open reading frame of 1,161 bp encoding a prepropeptide consisting of 387 amino acids with a calculated molecular mass of 40.923 kDa. Homology comparison of the deduced amino acid sequence of Eap with other known proteins indicated that Eap encode an extracellular protease that belongs to the subtilase family of serine protease (Family S8). A comparative homology model of the Engyodontium album protease (EAP) was developed using the crystal structure of proteinase K. The model revealed that EAP has broad substrate specificity similar to Proteinase K with preference for bulky hydrophobic residues at P1 and P4. Also, EAP is suggested to have two disulfide bonds and more than two Ca²⁺ binding sites in its 3D structure; both of which are assumed to contribute to the thermostable nature of the protein.     

Kinetic stabilization of Bacillus licheniformis alpha-amylase through introduction of hydrophobic residues at the surface

It is generally assumed that in proteins hydrophobic residues are not favorable at solvent-exposed sites, and that amino acid substitutions on the surface have little effect on protein thermostability. Contrary to these assumptions, we have identified hyperthermostable variants of Bacillus licheniformis alpha-amylase (BLA) that result from the incorporation of hydrophobic residues at the surface. Under highly destabilizing conditions, a variant combining five stabilizing mutations unfolds 32 times more slowly and at a temperature 13 degrees C higher than the wild-type. Crystal structure analysis at 1.7 A resolution suggests that stabilization is achieved through (a) extension of the concept of increased hydrophobic packing, usually applied to cavities, to surface indentations, (b) introduction of favorable aromatic-aromatic interactions on the surface, (c) specific stabilization of intrinsic metal binding sites, and (d) stabilization of a beta-sheet by introducing a residue with high beta-sheet forming propensity. All mutated residues are involved in forming complex, cooperative interaction networks that extend from the interior of the protein to its surface and which may therefore constitute "weak points" where BLA unfolding is initiated. This might explain the unexpectedly large effect induced by some of the substitutions on the kinetic stability of BLA. Our study shows that substantial protein stabilization can be achieved by stabilizing surface positions that participate in underlying cooperatively formed substructures. At such positions, even the apparently thermodynamically unfavorable introduction of hydrophobic residues should be explored.     

The structure of human lipoprotein A-I. Evidence for the "belt" model.

The two main competing models for the structure of discoidal lipoprotein A-I complexes both presume that the protein component is helical and situated around the perimeter of a lipid bilayer disc. However, the more popular "picket fence" model orients the protein helices perpendicular to the surface of the lipid bilayer, while the alternative "belt" model orients them parallel to the bilayer surface. To distinguish between these models, we have investigated the structure of human lipoprotein A-I using a novel form of polarized internal reflection infrared spectroscopy that can characterize the relative orientation of protein and lipid components in the lipoprotein complexes under native conditions. Our results verify lipid bilayer structure in the complexes and point unambiguously to the belt model.     

Role of "anti-restriction" motif in functional activity of anti-restriction protein ArdA pKM101 (incN)

A number of mutant forms of the antirestriction protein ArdA encoded by the ardA gene located in a transmissive IncN plasmid pKM101 have been constructed. Proteins belonging to the Ard family are specific inhibitors of type I restriction--modification enzymes. Single mutational substitutions of negatively charged amino acid residues located in the "antirestriction motif" with hydrophobic alanine, E134A, E137A, D144A, or a double substitution E134A, E137A do not affect the antirestriction activity (Ard) of ArdA but almost completely abolish the antimodification activity (Amd). Mutational substitutions F107D and A110D in the assumed interface ArdA, which determines contact between monomers in the active dimer (Ard)2, cause an approximately 100-fold decrease in the antirestriction protein activity. It is hypothesized that the ArdA protein forms two complexes with the type I restriction--modification enzyme (R2M2S): (1) with a specific region in the S subunit involved in contact with the sK site in DNA; and (2) with a nonspecific region in the R subunit involved in DNA translocation and degradation by restriction endonucleases. The association of ArdA with the specific region inhibits restriction endonuclease and methyltransferase activities simultaneously, whereas the association of ArdA with a nonspecific region inhibits only restriction endonuclease activity of the R2M2S enzyme.