Amino acid residue: is it structural or functional?

Specific Amino Acid Residues in the Three Loops of Snake Cytotoxins Determine Their Membrane Activity and Provide a Rationale for a New Classification of These Toxins

Cytotoxins (CTs) are three-finger membrane-active toxins present mainly in cobra venom. Our analysis of the available CT amino acid sequences, literature data on their membrane activity, and conformational equilibria in aqueous solution and detergent micelles allowed us to identify specific amino acid residues which interfere with CT incorporation into membranes. They include Pro9, Ser28, and Asn/Asp45 within the N-terminal, central, and C-terminal loops, respectively. There is a hierarchy in the effect of these residues on membrane activity: Pro9 > Ser28 > Asn/Asp45. Taking into account all the possible combinations of special residues, we propose to divide CTs into eight groups. Group 1 includes toxins containing all of the above residues. Their representatives demonstrated the lowest membrane activity. Group 8 combines CTs that lack these residues. For the toxins from this group, the greatest membrane activity was observed. We predict that when solely membrane activity determines the cytotoxic effects, the activity of CTs from a group with a higher number should exceed that of CTs from a group with a lower number. This classification is supported by the available data on the cytotoxicity and membranotropic properties of CTs. We hypothesize that the special amino acid residues within the loops of the CT molecule may indicate their involvement in the interaction with non-lipid targets.

Functionally and structurally relevant residues of enzymes: are they segregated or overlapping?

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.

Noncovalent Cross-links in Context with Other Structural and Functional Elements of Proteins

Proteins are heteropolymers with evolutionary selected native sequences of residues. These native sequences code for unique and stable 3D structures indispensable for biochemical activity and for proteolysis resistance, the latter which guarantees an appropriate lifetime for the protein in the protease rich cellular environment. Cross-links between residues close in space but far in the primary structure are required to maintain the folded structure of proteins. Some of these cross-links are covalent, most frequently disulfide bonds, but the majority of the cross-links are sets of cooperative noncovalent long-range interactions. In this paper we focus on special clusters of noncovalent long-range interactions: the Stabilization Centers (SCs). The relation between the SCs and secondary structural elements as well as the relation between SCs and functionally important regions of proteins are presented to show a detailed picture of these clusters, which are believed to be primarily responsible for major aspects of protein stability.

Interaction of cardiotoxins with membranes: a molecular modeling study

Incorporation of beta-sheet proteins into membrane is studied theoretically for the first time, and the results are validated by the direct experimental data. Using Monte Carlo simulations with implicit membrane, we explore spatial structure, energetics, polarity, and mode of insertion of two cardiotoxins with different membrane-destabilizing activity. Both proteins, classified as P- and S-type cardiotoxins, are found to retain the overall "three-finger" fold interacting with membrane core and lipid/water interface by the tips of the "fingers" (loops). The insertion critically depends upon the structure, hydrophobicity, and electrostatics of certain regions. The simulations reveal apparently distinct binding modes for S- and P-type cardiotoxins via the first loop or through all three loops, respectively. This rationalizes an earlier empirical classification of cardiotoxins into S- and P-type, and provides a basis for the analysis of experimental data on their membrane affinities. Accomplished with our previous simulations of membrane alpha-helices, the computational method may be used to study partitioning of proteins with diverse folds into lipid bilayers.

Membrane binding motif of the P-type cardiotoxin

Carditoxins (CTXs) from cobra snake venoms, the basic 60-62 residue all-beta sheet polypeptides, are known to bind to and impair the function of cell membranes. To assess the membrane induced conformation and orientation of CTXs, the interaction of the P-type cardiotoxin II from Naja oxiana snake venom (CTII) with perdeuterated dodecylphosphocholine (DPC) was studied using ( 1 )H-NMR spectroscopy and diffusion measurements. Under conditions where the toxin formed a well-defined complex with DPC, the spatial structure of CTII with respect to the presence of tightly bound water molecules in loop II, was calculated using the torsion angle dynamics program DYANA. The structure was found to be similar, except for subtle changes in the tips of all three loops, to the previously described "major" form of CTII in aqueous solution illustrated by the "trans" configuration of the Val7-Pro8 peptide bond. No "minor" form with the "cis" configuration of the above bond was found in the micelle-bound state. The broadening of the CTII backbone proton signals by 5, 16-doxylstearate relaxation probes, together with modeling based on the spatial structure of CTII, indicated a periphery mode of binding of the toxin molecule to the micelle and revealed its micelle interacting domain. The latter includes a hydrophobic region of CTII within the extremities of loops I and III (residues 5-11, 46-50), the basement of loop II (residues 24-29,31-37) and the belt of polar residues encircling these loops (lysines 4,5,12,23,50, serines 11,46, histidine 31, arginine 36). It is suggested that this structural motif and the mode of binding can be realized during interaction of CTXs with lipid and biological membranes.

Stabilization of proteins by enhancement of inter-residue hydrophobic contacts: lessons of T4 lysozyme and barnase

Although the hydrophobic interactions are considered as the main contributors to the protein stability, not much examples of protein stabilization by rational increasing of this type of interactions still can be found in literature. This is partly due to the lack of proper theoretical "measure" of hydrophobic interactions and their changes upon mutations. In the present paper the molecular hydrophobicity potential approach is used to assess how the changes in type and the strength of inter-residue contacts upon single amino acid mutations are correlated with the changes in thermodynamic stability of T4 lysozyme and barnase mutants, and which factors affect these correlations. Mutations changing unfavorable hydrophilic-to-hydrophobic contacts into favorable hydrophobic were found to enhance the thermodynamic stability in more than 81 % of cases, if these mutations do not create steric bumps and do not involve proline residues and hydrogen-bonded side-chains. Mutations increasing hydrophobic contributions (according to molecular hydrophobicity potential formalism) lead to increase of thermodynamic stability in more than 94% of cases for certain type of mutations (i.e., mutations not involving charged residues, Pro and residues with side-chain hydrogen bonds, when these mutations do not introduce steric bumps and do not involve strongly exposed residues and residues situated at helix N- and C-cap positions). For this type of mutations the correlation was found between the change in hydrophobic contributions of mutated residues deltaCphob and thermodynamic parameters deltaTm (change in melting temperature) and deltadeltaG (change in free energy of unfolding). Although the correlation coefficients were larger if the experimental structures of mutants were used for the calculations (correlation coefficients r(exp) deltaC,deltaT = .85 and r(exp) deltaC,deltadeltaG = 0.87) than if the modeled structures were used instead (r(mod) deltaC,deltaT = 0.74 and r(mod)deltaC,deltadeltaG = 0.76), the modelled structures of mutants in the vast majority of cases can be used for qualitative predictition of the protein stabilization. Basing on the analysis of mutations increasing hydrophobic contributions in T4 lysozyme the substitution matrix was derived, which can be used to decide which new residue should be put instead the old one to increase the stability of protein. The estimation shows that the number of potential mutation sites for enhancement of hydrophobic interactions in T4 lysozyme is quite large, and only approximately 10 per cent of them were studied thus far. Basing on the current analysis of T4 lysozyme and barnase mutations the algorithm for increasing of protein stability via increasing of hydrophobic interactions for the proteins with known spatial structure is proposed.

Recognizing misfolded and distorted protein structures by the assumption-based similarity score

A new similarity score (sigma-score) is proposed which is able to find the correct protein structure among the very close alternatives and to distinguish between correct and deliberately misfolded structures. This score is based on the general principle 'similar likes similar', and it favors hydrophobic and hydrophilic contacts, and disfavors hydrophobic-to-hydrophilic contacts in proteins. The values of sigma-scores calculated for the high-resolution protein structures from the representative set are compared with those of alternatives: (i) very close alternatives which are only slightly distorted by conformational energy minimization in vacuo; (ii) alternatives with subsequently growing distortions, generated by molecular dynamics simulations in vacuo; (iii) structures derived by molecular dynamics simulation in solvent at 300 K; (iv) deliberately misfolded protein models. In nearly all tested cases the similarity score can successfully distinguish between experimental structure and its alternatives, even if the root mean square displacement of all heavy atoms is less than 1 A. The confidence interval of the similarity score was estimated using the high-resolution X-ray structures of domain pairs related by non-crystallographic symmetry. The similarity score can be used for the evaluation of the general quality of the protein models, choosing the correct structures among the very close alternatives, characterization of models simulating folding/unfolding, etc.

Molecular modeling of HIV-1 coreceptor CCR5 and exploring of conformational space of its extracellular domain in molecular dynamics simulation

The chemokine receptor CCR5 functions as a major fusion coreceptor for macrophage-tropic human immunodeficiency virus entry into cell. Here we report a three-dimensional model of CCR5 built using molecular modeling approach. Because the virus binds to extracellular domain of the receptor, special attention was given to conformational flexibility, hydrogen bonding, and environmental polarity properties of this protein part. Such data were obtained in the result of molecular dynamics study of the extracellular domain. It was shown that during the simulation the extracellular segments form a compact globular domain with numerous long-range hydrogen bonds between them. First loop of the receptor stays quite rigid while N-terminal region and loops 2, 3 are rather flexible. A number of amino acid residues disposed in unfavourable environment and, therefore, potentially involved in binding of CCR5 to viral glycoproteins and chemokines, was delineated. Comparison of the results with available experimental data permits a proposal that such residues in loop-1 and N-terminal part of the receptor are important for HIV-1 entry, while those in loops 2 and 3 participate in ligand binding. Perspectives of rational alteration of virus-binding activity of CCR5 are discussed.

A new method to characterize hydrophobic organization of proteins: application to rational protein engineering of barnase

We present a new algorithm for characterization of protein spatial structure basing on the molecular hydrophobicity potential approach. The method is illustrated by the analysis of three-dimensional structure of barnase and barnase-barstar complex. Current approach enables identification of amino acid residues situated in unfavorable environment (these residues may be "active" for binding), and to map quantitatively hydrophobic, hydrophilic and unfavorable hydrophobic-hydrophilic intra- and inter-molecular contacts involving backbone and side-chain segments of amino acid residues. Calculation of individual contributions of amino acid residues to such contacts permits identification of structurally-important residues. The contact plots obtained with molecular hydrophobicity potential calculations, provide easy rules to choose sites for mutations, which can increase a strength of intra- or inter-molecular hydrophobic interactions. The unfavorable hydrophobic-hydrophilic contact can be mutated to favorable hydrophobic, and already existing weak hydrophobic contact can be strengthen by increasing hydrophobicity of residues in contact. Basing on the analysis of the contact plots, we suggest several mutations of barnase which are supposed to increase intramolecular hydrophobic interactions, and thus might lead to increased stability of the protein. Part of these mutations was studied previously experimentally, and indeed stabilized barnase. The other of predicted mutations were not studied experimentally yet. Several new mutations of barnase and barstar are also proposed to enhance the hydrophobic interactions on their binding interface.