Biomedical subjects

M M Gromiha

Publications and source records attributed to M M Gromiha.

35 records · Page 2Linked to original sources

PubMed1999

Role of structural and sequence information in the prediction of protein stability changes: comparison between buried and partially buried mutations.

Predicting mutation-induced changes in protein stability is one of the greatest challenges in molecular biology. In this work, we analyzed the correlation between stability changes caused by buried and partially buried mutations and changes in 48 physicochemical, energetic and conformational properties. We found that properties reflecting hydrophobicity strongly correlated with stability of buried mutations, and there was a direct relation between the property values and the number of carbon atoms. Classification of mutations based on their location within helix, strand, turn or coil segments improved the correlation of mutations with stability. Buried mutations within beta-strand segments correlated better than did those in alpha-helical segments, suggesting stronger hydrophobicity of the beta-strands. The stability changes caused by partially buried mutations in ordered structures (helix, strand and turn) correlated most strongly and were mainly governed by hydrophobicity. Due to the disordered nature of coils, the mechanism underlying their stability differed from that of the other secondary structures: the stability changes due to mutations within the coil were mainly influenced by the effects of entropy. Further classification of mutations within coils, based on their hydrogen-bond forming capability, led to much stronger correlations. Hydrophobicity was the major factor in determining the stability of buried mutations, whereas hydrogen bonds, other polar interactions and hydrophobic interactions were all important determinants of the stability of partially buried mutations. Information about local sequence and structural effects were more important for the prediction of stability changes caused by partially buried mutations than for buried mutations; they strengthened correlations by an average of 27% among all data sets.

Amino Acids↗

PubMed1999

A simple method for predicting transmembrane alpha helices with better accuracy.

The prediction of a protein's structure from its amino acid sequence has been a long-standing goal of molecular biology. In this work, a new set of conformational parameters for membrane spanning alpha helices was developed using the information from the topology of 70 membrane proteins. Based on these conformational parameters, a simple algorithm has been formulated to predict the transmembrane alpha helices in membrane proteins. A FORTRAN program has been developed which takes the amino acid sequence as input and gives the predicted transmembrane alpha-helices as output. The present method correctly identifies 295 transmembrane helical segments in 70 membrane proteins with only two overpredictions. Furthermore, this method predicts all 45 transmembrane helices in the photosynthetic reaction center, bacteriorhodopsin and cytochrome c oxidase to an 86% level of accuracy and so is better than all other methods published to date.

Algorithms↗

PubMed1998

Analysis of hydrophobic and charged patches and influence of medium- and long-range interactions in molecular chaperones.

The amino acid composition of the aromatic residues Phe, Tyr and Trp are much less significant in chaperones and the residues Cys, Glu, His, Met and Pro vary significantly in chaperones compared to normal globular proteins. In the present work, we have analysed the hydrophobic and charged patches in molecular chaperones which provide more insight for a better understanding of chaperone folding. Also, we have investigated the role of medium- and long-range contacts in chaperones and the preference of amino acid residues influenced by these interactions. Furthermore, the role of hydrophobic and helix-forming residues and disulfide bonding in these interactions have been discussed.

Amino Acids↗

PubMed1998

Solvent accessibility studies on glycosaminoglycans.

Theoretical calculations on the solvent accessible surface area of the commonly occurring glycosaminoglycans (GAGs) have been carried out by employing the solvent accessibility technique of Lee and Richards. The results indicate that the average variation of solvent accessible surface area (ASA) for the different atoms is in the range of 2-28 A2. The average ASA for carbon (9.02 A2) and oxygen (19.3 A2) are relatively very high compared to that for nitrogen and sulfur. The total contribution of ASA by the buried atoms in all the molecules has been found to be small. Heparin complexed with the protein is less accessible to the solvent than the uncomplexed. The iduronic acid is also found to have a varied ASA due to its conformational flexibility. The sulfate groups in heparin have 50% of the total ASA of the molecule. Residue analysis indicates that in general D-configuration residues have higher ASA than L-configuration residues.

Carbohydrate Sequence↗

PubMed1998

Solvent accessibility analysis on the mutants of Hsc70 ATPase fragment.

Molecular chaperones are the cellular proteins which mediate the correct folding of other polypeptides. The concept of 'solvent accessibility' is one of the most powerful tools to understand the structure and stability of protein molecules. The hydrophobic variation of amino acid residues due to point mutations at many active sites of chaperone protein Hsc70 using solvent accessibility analysis is carried out. The numerical indices for several properties of amino acid residues, such as, reduction in accessibility, preference of amino acid residues in interior and surface parts, transfer free energy and the preference of amino acid residues to change their positions (buried/exposed) due to amino acid substitutions for Hsc70 and its mutants were set up. The accessibility of amino acid residues varies much between native and mutant proteins whereas there is no major changes on their conformations. The conformational stability for Hsc70 and its mutants were established and the computed hydrophobic free energy change is around 10 kcal/mol due to single amino acid substitution.

HSP70 Heat-Shock Proteins↗

PubMed1998

Importance of long-range interactions in (alpha/beta)8 barrel fold.

Protein structures are stabilized by both local and long-range interactions. In this work, we analyzed the importance of long-range interactions in (alpha/beta)8 barrel proteins in terms of residue distances. We found that the residues occurring in the range of 21-30 residues apart contribute more toward long-range contacts. Indeed, about 50% of successive strands in these proteins are found to occur at a sequential distance of 21-30 residues. The aromatic amino acid residues Phe, Trp, and Tyr prefer the 4-10 range and all other residues prefer the 21-30 range. Hydrophobic-hydrophobic residue pairs are the most preferred ones for long-range interactions and they may play a key role in the folding and stabilization of (alpha/beta)8 barrel proteins.

Amino Acids↗

PubMed1998

An analysis of the amino acid clustering pattern in (alpha/beta)8 barrel proteins.

The (alpha/beta)8 barrel proteins, in spite of having a common fold, do not show any sequence similarity. In order to understand the factors which are responsible for maintaining the common fold, the three-dimensional structures of 36 (alpha/beta)8 barrel proteins are analyzed for the presence of identical amino acid clusters or physicochemically similar clusters. The results reveal 14 identical amino acid clusters and a large number of physicochemically similar clusters. Further analysis of the similar clusters points to the conservation of secondary structures, the presence of pairs of residues occupying topologically equivalent secondary structures, and the presence of certain key residues which may play a vital role in directing and stabilizing the (alpha/beta)8 barrel fold.

Amino Acids↗

PubMed1998

Protein secondary structure prediction in different structural classes.

Information about the secondary structure of a protein can be helpful in understanding its native folded state. In previous work, it was shown that the medium-range interactions predominate in all-alpha class and the long-range interactions predominate in all-beta class proteins. Based on this, in this work the performance of several structure prediction methods in different structural classes of globular proteins was analyzed. It was found that all the methods predict the secondary structures of all-alpha proteins more accurately than other classes.

Protein Structure, Secondary↗

PubMed1997

The role of DNA bending in Cro protein-DNA interactions.

Binding energy of DNA-Cro protein complexes is analyzed in terms of DNA elasticity, using a sequence-dependent anisotropic bendability (SDAB) model of DNA, developed recently [M.M. Gromiha, M.G. Munteanu, A. Gabrielian and S. Pongor, J. Biol. Phys. 22(1996) 227-243.]. The protein is considered to bind aspecifically to DNA that reduces the freedom of movement in the DNA molecule. In cognate DNA, the Cro protein moves on to form specific interactions and bends DNA. A comparison of the experimental data [Y. Takeda, A. Sarai and V.M. Rivera, Proc. Natl. Acad. Sci. U.S.A. 86 (1989) 439-443.] with the calculated DNA stiffness data shows that delta G of the complex formation increases with stiffness of the ligand when the interactions are nonspecific ones, while an opposite trend is observed for specific binding. Both of these trends are in agreement with our approach using the SDAB model. A decomposition of the energy terms suggests that binding energy in the nonspecific case is used maily to compensate the free energy changes due to entropy lost by DNA, while the energy of specific interactions provide enough energy both to bend the DNA molecule and to change the conformation of the Cro protein upon ligand binding.

Base Sequence↗

PubMed1997

Identification of membrane spanning beta strands in bacterial porins.

The membrane assembly of outer membrane proteins is more complex than that of transmembrane helical proteins owing to the intervention of many charged and polar residues in the membrane. Accordingly, the predictive accuracy of transmembrane beta strands is considerably lower than that of transmembrane alpha helices. In this paper we develop a set of conformational parameters for membrane spanning beta strands. We formulate an algorithm to predict the transmembrane beta strands in the family of bacterial porins based on the conformational parameters and surrounding hydrophobicities of amino acid residues. A Fortran program has been developed which takes the amino acid sequence as the input file and gives the predicted transmembrane beta strand as output. The present method predicts at an accuracy level of 82% for all the bacterial porins considered.

Amino Acid Sequence↗

PubMed1996

Hydrophobic distribution and spatial arrangement of amino acid residues in membrane proteins.

The analysis of known three-dimensional structures of membrane proteins provides an opportunity to understand their structure and stability. In this article we analyse the hydrophobic variation of amino acid residues at various ranges in membrane and aqueous parts of membrane proteins. The numerical indices for several properties of amino acid residues in membrane proteins, such as surrounding hydrophobicity, gain in surrounding hydrophobicity, hydrophobic gain ratio, accessible surface area, preference of amino acid residues in the interior and surface parts, solvent accessible reduction ratio and buriedness, were set up. The relative preference of amino acid residues at various positions of membrane proteins were obtained in a very realistic approach.

Amino Acids↗

PubMed1995

Prediction of protein secondary structures from their hydrophobic characteristics.

Deciphering the native conformation of proteins from their amino acid sequences is one of the greatest challenges in the field of molecular biology. The successful prediction of structural class may help to improve the accuracy levels of structure (secondary and tertiary) predictive schemes in globular proteins. In our earlier works we developed a new surrounding hydrophobicity scale for the 20 amino acid residues applicable for both globular and membrane proteins and used it successfully to predict the transmembrane helical and strand segments in membrane proteins. In this article we propose (i) rules to predict the structural class of proteins and (ii) a new predictive scheme for forecasting secondary structures of globular proteins, with the use of the new hydrophobicity scale. This scheme predicts the structural class and secondary structures of globular proteins to 92 and 82% levels of accuracy, respectively, far better than the levels from other existing methods.

Algorithms↗

PubMed1994

On the conformational stability of oligonucleotide duplexes and tRNA molecules.

Thermodynamic experiments provide a wealth of data about the conformational stability, viz., the free energy difference (delta G) between folded and unfolded states of DNA/RNA duplexes. However, there is no acceptable view about how the various non-covalent forces contribute individually to the observed stability. In particular, the role of the hydrophobic force is not clearly known. In this paper we quantitatively enumerate the stability factors, hydrogen bonding, base stacking, van der Waals, electrostatic, and hydrophobic interactions from the knowledge of the crystal structures of 15 DNA/RNA duplexes and two tRNA molecules, and translate them into free energy contributions to the stability of nucleic acid systems. Taking the experimental delta G values and computed component free energy terms for a set of duplexes, we set up multiple regression equations to predict their stabilities. After back-check and validity tests, we apply this model to predict delta G values for a large number of duplexes and two tRNA molecules (tRNAphe and tRNAasp). There is excellent agreement between the theoretical predictions and experimental observations. The considered duplexes with four to 16 base-pairs and the tRNA molecules have delta G values in a narrow range, 5-20 kcal mol-1, a range seen in a variety of globular proteins. There is no relationship between delta G and N, the number of nucleotides in the molecule. Base-stacking, hydrogen bonding and van der Waals factors contribute significantly, whereas hydrophobic and electrostatic factors contribute, respectively, marginally and minimally. The major factor which gives sequence specificity is base-stacking. The new set of atomic solvation parameters (ASPs) derived to estimate hydrophobic free energy brings to light the dangers of using already available ASPs, which emphasize the role of the hydrophobic factor unrealistically.

Base Sequence↗

PubMed1993

Prediction of transmembrane beta-strands from hydrophobic characteristics of proteins.

The assembly of outer-membrane proteins consisting of beta-strands as transmembrane segments is somewhat more complex when compared to the assembly of inner membrane proteins having alpha-helices as transmembrane parts. This is probably due to the difference in the amino acid sequences of the transmembrane part strands and helices. Because of this feature, most predictive schemes which are successful in predicting transmembrane helical segments fail to predict transmembrane strand segments. Here we propose a new predictive scheme for forecasting the transmembrane strand segments in outer-membrane proteins with the use of the general surrounding hydrophobicity scale developed both for the globular and membrane proteins in the preceding article. Two major features of the scheme are (i) that it does not solely depend on the amphipathic character of a sequence segment while identifying it as a transmembrane strand, and it is capable of predicting strands in varied lengths, a facility to reflect the variation in the membrane surfaces. This scheme predicts the transmembrane beta-strands in porin from R. capsulatus at 76% accuracy (giving due weights to over- and under-predictions when compared to X-ray results). The predicted beta-structure contents in OmpA, porin from E. coli and maltoporin compared with the Raman spectroscopic results at 95% level. These accuracy levels are far superior than the levels obtained from other existing methods. Apart from the above four proteins for which experimental informations are available, ten other outer-membrane proteins, for which there is no information about their secondary structure, are considered in the forecast. The results are analysed and the common features in the folds of the set of fourteen outer-membrane proteins are deduced.

Amino Acid Sequence↗

PubMed1993

Prediction of transmembrane helices from hydrophobic characteristics of proteins.

Membrane proteins, requiring to be embedded into the lipid bilayers, have evolved to have amino acid sequences that will fold with a hydrophobic surface in contact with the alkane chains of the lipids and polar surface in contact with the aqueous phases on both sides of the membrane and the polar head groups of the lipids. It is generally assumed that the characteristics of the aqueous parts of the membrane proteins are similar to those of normal globular proteins, and the embedded parts are highly hydrophobic. In our earlier works, we introduced the concept of 'surrounding hydrophobicity' and developed a hydrophobicity scale for the 20 amino acid residues, and applied it successfully to the study of the family of globular proteins. In this work we use the concept of surrounding hydrophobicity to indicate quantitatively how the aqueous parts of membrane proteins compare with the normal globular proteins, and how rich the embedded parts are in their hydrophobic activity. We then develop a surrounding hydrophobicity scale applicable to membrane proteins, by mixing judicially the surrounding hydrophobicities observed in the crystals of the membrane protein, photosynthetic reaction center from the bacterium Rhodopseudomonas viridis, porin from Rhodobacter capsulatus and a set of 64 globular proteins. A predictive scheme based on this scale predicts from amino acid sequence, transmembrane segments in PRC and randomly selected 26 membrane proteins to 80% level of accuracy. This is a much higher predictive power when compared to the existing popular methods. A new procedure to measure the amphipathicity of sequence segments is proposed, and it is used to characterize the transmembrane parts of the sample membrane proteins.

Chemical Phenomena↗

PubMedDate not supplied

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.

Adenine↗

PubMedDate not supplied

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.

DNA↗