Influence of Amino Acid Mutations and Small Molecules on Targeted Inhibition of Proteins Involved in Cancer.

BACKGROUND: Protein-protein interactions (PPIs) are of crucial importance in regulating the biological processes of cells both in normal and diseased conditions. Significant progress has been made in targeting PPIs using small molecules and achieved promising results. However, PPI drug discovery should be further accelerated with better understanding of chemical space along with various functional aspects. OBJECTIVE: In this review, we focus on the advancements in computational research for targeted inhibition of protein-protein interactions involved in cancer. METHODS: Here, we mainly focused on two aspects: (i) understanding the key roles of amino acid mutations in epidermal growth factor receptor (EGFR) as well as mutation-specific inhibitors and (ii) design of small molecule inhibitors for Bcl-2 to disrupt PPIs. RESULTS: The paradigm of PPI inhibition to date reflect the certainty that inclination towards novel and versatile strategies enormously dictate the success of PPI inhibition. As the chemical space highly differs from the normal drug like compounds the lead optimization process has to be given the utmost priority to ensure the clinical success. Here, we provided a broader perspective on effect of mutations in oncogene EGFR connected to Bcl-2 PPIs and focused on the potential challenges. CONCLUSION: Understanding and bridging mutations and altered PPIs will provide insights into the alarming signals leading to massive malfunctioning of a biological system in various diseases. Finding rational elucidations from a pharmaceutical stand point will presumably broaden the horizons in future.

Prediction of protein stability upon point mutations.

Prediction of protein stability upon amino acid substitution is a challenging problem and it will be helpful for designing stable mutants. We have developed a thermodynamic database for proteins and mutants (ProTherm), which has more than 20000 thermodynamic data along with sequence and structure information, experimental conditions and literature information. It is freely accessible at http://gibk26.bse.kyutech.ac.jp/jouhou/protherm/protherm.html. Utilizing the database, we have analysed the relationship between amino acid properties and protein stability and developed different methods, such as average assignment method, distance and torsion potentials and decision tree models to discriminate the stabilizing and destabilizing mutants, and to predict the stability change upon mutation. Our method could distinguish the stabilizing and destabilizing mutants with an accuracy of 82 and 85% respectively from amino acid sequence and protein three-dimensional structure. We obtained the correlation of 0.70 and 0.87, between the experimental and predicted stability changes upon mutations, from sequence and structure respectively. Furthermore, we have developed different web servers for discrimination and prediction and they are freely accessible at http://bioinformatics.myweb.hinet.net/iptree.htm and http://cupsat.tu-bs.de/.

Thermodynamic database for protein-nucleic acid interactions (ProNIT).

MOTIVATION: Protein-nucleic acid interactions are fundamental to the regulation of gene expression. In order to elucidate the molecular mechanism of protein-nucleic acid recognition and analyze the gene regulation network, not only structural data but also quantitative binding data are necessary. Although there are structural databases for proteins and nucleic acids, there exists no database for their experimental binding data. Thus, we have developed a Thermodynamic Database for Protein-Nucleic Acid Interactions (ProNIT). RESULTS: We have collected experimentally observed binding data from the literature. ProNIT contains several important thermodynamic data for protein-nucleic acid binding, such as dissociation constant (K(d)), association constant (K(a)), Gibbs free energy change (DeltaG), enthalpy change (DeltaH), heat capacity change (DeltaC(p)), experimental conditions, structural information of proteins, nucleic acids and the complex, and literature information. These data are integrated into a relational database system together with structural and functional information to provide flexible searching facilities by using combinations of various terms and parameters. A www interface allows users to search for data based on various conditions, with different display and sorting options, and to visualize molecular structures and their interactions. AVAILABILITY: ProNIT is freely accessible at the URL http://www.rtc.riken.go.jp/jouhou/pronit/pronit.html.

Role of medium- and long-range interactions to the stability of the mutants of T4 lysozyme.

Inter-residue interactions play an important role to the folding and stability of protein molecules. In this work, we analyze the role of medium- and long-range interactions to the stability of T4 lysozyme mutants. We found that, in buried mutations, the increase in long-range contacts upon mutations destabilizes the protein, whereas, in surface mutations, the increase in long-range contacts increases the stability, indicating the importance of surrounding polar residues to the stability of surface mutations. Further, the increase in medium-range contacts decreases the stability of buried and surface mutations and a direct relationship is observed between the increase of medium-range contacts and increase in stability for partially buried/exposed mutations. Moreover, the relationship between amino acid properties and stability of T4 lysozyme mutants at positions Ile3, Phe53, and Leu99 showed that the effect of medium- and long-range contacts is less for buried mutations and the inter-residue contacts have significant correlation with the stability of partially buried mutations.

Role of medium--and long-range interactions in discriminating globular and membrane proteins.

The analysis of inter-residue interactions in protein structures provides considerable insight to understand their folding and stability. We have previously analyzed the role of medium- and long-range interactions in the folding of globular proteins. In this work, we study the distinct role of such interactions in the three-dimensional structures of membrane proteins. We observed a higher number of long-range contacts in the termini of transmembrane helical (TMH) segments, implying their role in the stabilization of helix-helix interactions. The transmembrane strand (TMS) proteins are having appreciably higher long-range contacts than that in all-beta class of globular proteins, indicating closer packing of the strands in TMS proteins. The residues in membrane spanning segments of TMH proteins have 1.3 times higher medium-range contacts than long-range contacts whereas that of TMS proteins have 14 times higher long-range contacts than medium-range contacts. Residue-wise analysis indicates that in TMH proteins, the residues Cys, Glu, Gly, Pro, Gln, Ser and Tyr have higher long-range contacts than medium-range contacts in contrast with all-alpha class of globular proteins. The charged residue pairs have higher medium-range contacts in all-alpha proteins, whereas hydrophobic residue pairs are dominant in TMH proteins. The information on the preference of residue pairs to form medium-range contacts has been successfully used to discriminate the TMH proteins from all-alpha proteins. The statistical significance of the results obtained from the present study has been verified using randomized structures of TMH and TMS protein templates.

Factors influencing the stability of alpha-helices and beta-strands in thermophilic ribonuclease H.

Understanding the influence of structural parameters is crucial to enhance the thermal stability of proteins. In this work, the stability (deltaG) of residues in different secondary structures of Ribonuclease H (RNase H) has been analyzed with 48 amino acid properties. The properties reflecting hydrophobicity show a good correlation with stability. Further, the linear distribution of surrounding hydrophobicity in alpha-helices, obtained from the three dimensional structure of thermophilic RNase H, agrees well with experimental deltaG values. Moreover, the stability parameters correlate better in alpha-helices than those did in beta-strand segments. Multiple regression analysis, incorporating combinations of three properties from among all possible combinations of the 48 properties, increased the correlation coefficient to 0.77.

Distribution of amino acid residues and residue-residue contacts in molecular chaperones.

The amino acid distribution and residue-residue contacts in molecular chaperones are different when compared to normal globular proteins. The study of molecular chaperones reveals a different surrounding environment to exist for the residues Cys, Trp, and His which may play an important role in determining the chaperone structures. Unlike globular proteins, it has been observed that a one-to-one correspondence between the amino acid distribution in a sequence and the structures of molecular chaperones. The preference of amino acid residues surrounding all 20 types of residues in secondary structures and their accessible surface areas have been analysed.

Important inter-residue contacts for enhancing the thermal stability of thermophilic proteins.

Proteins from thermophilic organisms exhibit high thermal stability, but have structures that are very similar to their mesophilic homologues. In order to gain insight into the basis of thermostability, we have analyzed the medium- and long-range contacts in mesophilic and thermophilic proteins of 16 different families. We found that the thermophiles prefer to have contacts between residues with hydrogen-bond-forming capability. Apart from hydrophobic contacts, more contacts are observed between polar and non-polar residues in thermophiles than mesophiles. Residue-wise analysis showed that Tyr has good contacts with several other residues, and Cys has considerably higher long-range contacts in thermophiles compared with mesophiles. Furthermore, the residues occurring in the range of 31-34 residues apart in the sequence contribute significant long-range contacts to the stability of thermophilic proteins.

Comparison between long-range interactions and contact order in determining the folding rate of two-state proteins: application of long-range order to folding rate prediction.

The contact order is believed to be an important factor for understanding protein folding mechanisms. In our earlier work, we have shown that the long-range interactions play a vital role in protein folding. In this work, we analyzed the contribution of long-range contacts to determine the folding rate of two-state proteins. We found that the residues that are close in space and are separated by at least ten to 15 residues in sequence are important determinants of folding rates, suggesting the presence of a folding nucleus at an interval of approximately 25 residues. A novel parameter "long-range order" has been proposed to predict protein folding rates. This parameter shows as good a relationship with the folding rate of two-state proteins as contact order. Further, we examined the minimum limit of residue separation to determine the long-range contacts for different structural classes. We observed an excellent correlation between long-range order and folding rate for all classes of globular proteins. We suggest that in mixed-class proteins, a larger number of residues can serve as folding nuclei compared to all-alpha and all-beta proteins. A simple statistical method has been developed to predict the folding rates of two-state proteins using the long-range order that produces an agreement with experimental results that is better or comparable to other methods in the literature.

Clusters in alpha/beta barrel proteins: implications for protein structure, function, and folding: a graph theoretical approach.

The alpha/beta barrel fold is adopted by most enzymes performing a variety of catalytic reactions, but with very low sequence similarity. In order to understand the stabilizing interactions important in maintaining the alpha/beta barrel fold, we have identified residue clusters in a dataset of 36 alpha/beta barrel proteins that have less than 10% sequence identity within themselves. A graph theoretical algorithm is used to identify backbone clusters. This approach uses the global information of the nonbonded interaction in the alpha/beta barrel fold for the clustering procedure. The nonbonded interactions are represented mathematically in the form of an adjacency matrix. On diagonalizing the adjacency matrix, clusters and cluster centers are obtained from the highest eigenvalue and its corresponding vector components. Residue clusters are identified in the strand regions forming the beta barrel and are topologically conserved in all 36 proteins studied. The residues forming the cluster in each of the alpha/beta protein are also conserved among the sequences belonging to the same family. The cluster centers are found to occur in the middle of the strands or in the C-terminal of the strands. In most cases, the residues forming the clusters are part of the active site or are located close to the active site. The folding nucleus of the alpha/beta fold is predicted based on hydrophobicity index evaluation of residues and identification of cluster centers. The predicted nucleation sites are found to occur mostly in the middle of the strands. Proteins 2001;43:103-112.

Solvent accessibility studies on polysaccharides.

The solvent accessible surface area (ASA) of the polysaccharides, namely (i) carrageenan (1CAR); (ii) agarose (1AGA); (iii) guaran (GUR); (iv) capsular polysaccharide (1CAP); and (v) hyaluronan (1HUA), have been computed using the solvent accessibility technique of Lee and Richards. The results show that the average variation of ASA for the various atoms in the molecules lie in the range 1-30 A(2). Irrespective of position of sulfation, either at two or four in the sugar residues in 1CAR, the charged groups interact almost equally with the solvent. The ASA values for the chains A and B in 1AGA and 1CAR indicate that there are not much interchain interactions and the chains in both the molecules interact equally with the solvent. Residue-wise analysis indicates that the ASAs of residues vary alternately, high-low-high value pattern that is similar to that of the hydrophobic behaviour of beta-strands in proteins. The results also suggest that in these polysaccharides D-configuration residues have higher ASA than L-configuration residues.

Thermodynamic databases for proteins and protein-nucleic acid interactions.

Thermodynamic data regarding proteins and their interactions are important for understanding the mechanisms of protein folding, protein stability, and molecular recognition. Although there are several structural databases available for proteins and their complexes with other molecules, databases for experimental thermodynamic data on protein stability and interactions are rather scarce. Thus, we have developed two electronically accessible thermodynamic databases. ProTherm, Thermodynamic Database for Proteins and Mutants, contains numerical data of several thermodynamic parameters of protein stability, experimental methods and conditions, along with structural, functional, and literature information. ProNIT, Thermodynamic Database for Protein-Nucleic Acid Interactions, contains thermodynamic data for protein-nucleic acid binding, experimental conditions, structural information of proteins, nucleic acids and the complex, and literature information. These data have been incorporated into 3DinSight, an integrated database for structure, function, and properties of biomolecules. A WWW interface allows users to search for data based on various conditions, with different display and sorting options, and to visualize molecular structures and their interactions. These thermodynamic databases, together with structural databases, help researchers gain insight into the relationship among structure, function, and thermodynamics of proteins and their interactions, and will become useful resources for studying proteins in the postgenomic era.

Evaluation of free energy landscape for base-amino acid interactions using ab initio force field and extensive sampling.

Structural data of protein-DNA complex show redundancy and flexibility in base-amino acid interactions. To understand the origin of the specificity in protein-DNA recognition, we calculated the interaction free energy, enthalpy, entropy, and minimum energy maps for AT-Asn, GC-Asn, AT-Ser, and GC-Ser by means of a set of ab initio force field with extensive conformational sampling. We found that the most preferable interactions in these pairs are stabilized by hydrogen bonding, and are mainly enthalpy driven. However, minima in the free energy maps are not necessarily the same as those in the minimum energy map or enthalpy maps, due to the entropic effect. The effect of entropy is particularly important in the case of GC-Asn. Experimentally observed structures of base-amino acid interactions are within preferable regions in the calculated free energy maps, where there are many different interaction configurations with similar energy. The full geometry optimization procedure using ab initio molecular orbital method was applied to get the optimal interaction geometries for AT-Asn, GC-Asn, AT-Ser, and GC-Ser. We found that there are various base-amino acid combinations with similar interaction energies. These results suggest that the redundancy and conformational flexibility in the base-amino acid interactions play an important role in the protein-DNA recognition.

Importance of surrounding residues for protein stability of partially buried mutations.

For understanding the factors influencing protein stability, we have analyzed the relationship between changes in protein stability caused by partially buried mutations and changes in 48 physico-chemical, energetic and conformational properties of amino acid residues. Multiple regression equations were derived to predict the stability of protein mutants and the efficiency of the method has been verified with both back-check and jack-knife tests. We observed a good agreement between experimental and computed stabilities. Further, we have analyzed the effect of sequence window length from 1 to 12 residues on each side of the mutated residue to include the sequence information for predicting protein stability and we found that the preferred window length for obtaining the highest correlation is different for each secondary structure; the preferred window length for helical, strand and coil mutations are, respectively, 0, 9 and 4 residues on both sides of the mutant residues. However, all the secondary structures have significant correlation for a window length of one residue on each side of the mutant position, implying the role of short-range interactions. Extraction of surrounding residue information for various distances (3 to 20A) around the mutant position showed the highest correlation at 8A, 6A and 7A, respectively, for mutations in helical, strand and coil segments. Overall, the information about the surrounding residues within the sphere of 7 to 8A, may explain better the stability in all subsets of partially buried mutations implying that this distance is sufficient to accommodate the residues influenced by major intramolecular interactions for the stability of protein structures.

On the thermal unfolding character of globular proteins.

A theoretical model is presented to study the stepwise thermal unfolding of globular proteins using the stabilizing/destabilizing characters of amino acid residues in protein crystals. A multiple regression relation connecting the melting temperature and the amounts of stabilizing and destabilizing groups of residues in a protein, when used for the thermal behavior of peptide segments, provides reliable results on the stepwise unfolding nature of the protein. In ribonuclease A, the shell residues 16-22 are predicted to unfold earlier in the temperature range 30-45 degrees C; the beta-sheet structures undergo thermal denaturation as a single cooperative unit and there is evidence indicating the segment 106-118 as a nucleation site. In ribonuclease S, the S-peptide unfolds earlier than S-protein. The predicted average and the range of melting temperatures, and the folding pathways of a set of globular proteins, agree very well with the experimental results. The results obtained in the present study indicate that (i) most of the nucleation parts possess high relative thermal stability, (ii) the unfolded state retains some residual structure, and (iii) some segments undergo gradual and overlapping thermal denaturation.

Structural class prediction: an application of residue distribution along the sequence.

Deciphering the native conformation of proteins from their amino acid sequences is one of the most challenging problems in molecular biology. Information on the secondary structure of a protein can be helpful in understanding its native folded state. In our earlier work on molecular chaperones, we have analyzed the hydrophobic and charged patches, short-, medium- and long-range contacts and residue distributions along the sequence. In this article, we have made an attempt to predict the structural class of globular and chaperone proteins based on the information obtained from residue distributions. This method predicts the structural class with an accuracy of 93 and 96%, respectively, for the four- and three-state models in a training set of 120 globular proteins, and 90 and 96%, respectively, for a test set of 80 proteins. We have used this information and methodology to predict the structural classes of chaperones. Interestingly most of the chaperone proteins are predicted under alpha/beta or mixed folding type.

ProTherm, version 2.0: thermodynamic database for proteins and mutants.

ProTherm 2.0 is the second release of the Thermo-dynamic Database for Proteins and Mutants that includes numerical data for several thermodynamic parameters, structural information, experimental methods and conditions, functional and literature information. The present release contains >5500 entries, an approximately 67% increase over the previous version. In addition, we have included information about reversibility of data, details about buffer and ion concentrations and the surrounding residues in space for all mutants. A WWW interface enables users to search data based on various conditions with different sorting options for outputs. Further, ProTherm has links with other structural and literature databases, and the mutation sites and surrounding residues are automatically mapped on the structures and can be directly viewed through 3DinSight developed in our laboratory. The ProTherm database is freely available through the WWW at http://www.rtc.riken.go.jp/protherm.html

Structure based sequence dependent stiffness scale for trinucleotides: a direct method.

A new set of stiffness parameters for all the 32trinucleotide units has been set up directly from thethree dimensional structures of DNA molecules. It wasobserved that GAC/GTC is the stiffest trinucleotideand ACC/GGT is the most flexible one. The averagestiffness values computed for a set of operatorsequences using the new parameters correlate very wellwith the protein-DNA binding specificity and bindingfree energy change of 434 repressor and Cro repressor,respectively. The new structure based stiffness scalecan explain the protein-DNA binding specificity to thelevel of 0.92.

Important amino acid properties for enhanced thermostability from mesophilic to thermophilic proteins.

Understanding the role of various interactions in enhancing the thermostability of proteins is important not only for clarifying the mechanism of protein stability but also for designing stable proteins. In this work, we have analyzed the thermostability of 16 different families by comparing mesophilic and thermophilic proteins with 48 various physicochemical, energetic and conformational properties. We found that the increase in shape, s (location of branch point in side chain) increases the thermostability, whereas, an opposite trend is observed for Gibbs free energy change of hydration for native proteins, GhN, in 14 families. A good correlation is observed between these two properties and the simultaneous increases of -GhN and s is necessary to enhance the thermostability from mesophile to thermophile. The increase in shape, which tends to increase with increasing number of carbon atoms both for polar and non-polar residues, may generate more packing and compactness, and the position of beta and higher order branches may be important for better packing. On the other hand, the increase in -GhN in thermophilic proteins increases the solubility of the proteins. This tendency counterbalances the increases in insolubility and unfolding heat capacity change due to the increase in the number of carbon atoms. Thus, the present results suggest that the stability of thermophilic proteins may be achieved by a balance between better packing and solubility.

Influence of medium and long range interactions in protein folding.

Protein structures are stabilized by both local and long range interactions. In this work, we analyze the residue-residue contacts and the role of medium- and long-range interactions in globular proteins belonging to different structural classes. The results show that while medium range interactions predominate in all-alpha class proteins, long-range interactions predominate in all-beta class. Based on this, we analyze the performance of several structure prediction methods in different structural classes of globular proteins and found that all the methods predict the secondary structures of all-alpha proteins more accurately than other classes. Also, we observed that the residues occurring in the range of 21-30 residues apart contributes more towards long-range contacts and about 85% of residues are involved in long-range contacts. Further, the preference of residue pairs to the folding and stability of globular proteins is discussed.