Understanding Drug Resistance of Wild-Type and L38HL Insertion Mutant of HIV-1 C Protease to Saquinavir.

Acquired immunodeficiency syndrome (AIDS) is one of the most challenging infectious diseases to treat on a global scale. Understanding the mechanisms underlying the development of drug resistance is necessary for novel therapeutics. HIV subtype C is known to harbor mutations at critical positions of HIV aspartic protease compared to HIV subtype B, which affects the binding affinity. Recently, a novel double-insertion mutation at codon 38 (L38HL) was characterized in HIV subtype C protease, whose effects on the interaction with protease inhibitors are hitherto unknown. In this study, the potential of L38HL double-insertion in HIV subtype C protease to induce a drug resistance phenotype towards the protease inhibitor, Saquinavir (SQV), was probed using various computational techniques, such as molecular dynamics simulations, binding free energy calculations, local conformational changes and principal component analysis. The results indicate that the L38HL mutation exhibits an increase in flexibility at the hinge and flap regions with a decrease in the binding affinity of SQV in comparison with wild-type HIV protease C. Further, we observed a wide opening at the binding site in the L38HL variant due to an alteration in flap dynamics, leading to a decrease in interactions with the binding site of the mutant protease. It is supported by an altered direction of motion of flap residues in the L38HL variant compared with the wild-type. These results provide deep insights into understanding the potential drug resistance phenotype in infected individuals.

Molecular dynamics simulations of cognate and non-cognate AspRS-tRNAAsp complexes.

Aspartyl tRNA synthetase (AspRS), one of the 20 aminoacyl-tRNA synthetases, plays an important role in protein synthesis by catalyzing the aminoacylation reaction and synthesises Aspartyl-tRNA (tRNAAsp). A typical three-dimensional structure of AspRS comprises three distinct domains for the recognition of cognate tRNA and catalysis, namely, anti-codon binding domain/N-terminal domain, hinge domain and catalytic domain through their interactions with anti-codon loop, D-stem and acceptor arm of cognate tRNA, respectively. In this work, we have studied the structural characteristics of each domain of AspRS to understand the recognition mechanism of tRNAAsp using molecular dynamics simulations. The dynamics of AspRS-tRNAAsp complexes from E.coli (cognate and non-cognate), S.cerevisiae (cognate) and T.thermophilus (non-cognate) were compared to understand the differences in recognition of cognate and non-cognate tRNAs. Our results explain that the conformational changes associated with the recognition of tRNA occur only in the cognate complexes. Among the cognate complexes, the conformational changes in yeast AspRS are highly controlled during tRNAAsp recognition than that of in the E. coli AspRS. Moreover, the functional motions required for the tRNA recognition are observed only in the cognate complexes, and the conformational changes in AspRS and their recognition of tRNAAsp are organism specific.Communicated by Ramaswamy H. Sarma.

Computational approaches for identifying potential inhibitors on targeting protein interactions in drug discovery.

In the era of big data, the interplay of artificial and human intelligence is the demanding job to address the concerns involving exchange of decisions between both sides. Drug discovery is one of the key sources of the big data, which involves synergy among various computational methods to achieve a clinical success. Rightful acquisition, mining and analysis of the data related to ligand and targets are crucial to accomplish reliable outcomes in the entire process. Novel designing and screening tactics are necessary to substantiate a potent and efficient lead compounds. Such methods are emphasized and portrayed in the current review targeting protein-ligand and protein-protein interactions involved in various diseases with potential applications.

Molecular dynamics simulation studies on influenza A virus H5N1 complexed with sialic acid and fluorinated sialic acid.

Abbreviations HA Hemagglutinin MD Molecular Dynamics MM-PBSA Molecular Mechanics Poisson-Boltzmann Surface Area NA Neuraminidase NAMD Nanoscale Molecular Dynamic Simulation PMEMD Particle Mesh Ewald Molecular Dynamics RMSD Root-Mean-Square Deviation RMSF Root-Mean-Square Fluctuation SIA sialic acid VMD Visual Molecular Dynamics Communicated by Ramaswamy H. Sarma.

Identifying a potential receptor for the antibacterial peptide of sponge Axinella donnani endosymbiont.

Marine sponges and their associated bacteria are rich sources of novel secondary metabolites with therapeutic usefulness. In our earlier work, we have identified a novel antibacterial peptide from the marine sponge Axinella donnani endosymbiotic bacteria. In this work, we have carried out a comparative genomic analysis and identified a set of 60 proteins as probable receptor which is common in all the strains. The analysis on binding substrate showed that ß barrel assembly machinery (BamA) of the outer membrane protein 85 (omp85) superfamily is a potential receptor protein for the antibacterial peptide. It plays a central role in OMP biogenesis, especially in cell viability. Further, the triplet and quartet motifs RGF and YGDG, respectively in L6 loop are conserved over all the strains and these conserved residues interact with antibacterial peptide to inhibit the BamA function, which is essential for OMP biogenesis.

Investigations of Ramachandran disallowed conformations in protein domain families.

In peptide and protein structures, occurrence of (φ,ψ) angles in the disallowed region of the Ramachandran map almost always suggests local regions of error or poor accuracy. However, very rarely genuine disallowed conformations occur as noted in the current study in proteins of known structure available at ultra-high resolution (≤1.2 Å). In the current work, extent of conservation of genuine disallowed conformations in homologous proteins of known structures has been analyzed. From a dataset of 124 protein domain families, with structure of at least one constituent member in each family available at a resolution of 1.2 Å or better, we have analyzed the conservation of 221 disallowed conformations. It is observed that the disallowed conformation is only moderately conserved in protein domain families. In the gross dataset no particular residue type adopting disallowed conformation elicit high conservation of residue type though there are alignment positions in the dataset with complete conservation of both the residue type and the disallowed conformation. Conserved disallowed conformation in protein domain families play biologically significant role in roughly 50% of the cases. The residues with the disallowed conformation or its flanking residues are often located within or around the functional site of the protein.

Sparsely populated residue conformations in protein structures: revisiting "experimental" Ramachandran maps.

The Ramachandran map clearly delineates the regions of accessible conformational (φ-ψ) space for amino acid residues in proteins. Experimental distributions of φ, ψ values in high-resolution protein structures, reveal sparsely populated zones within fully allowed regions and distinct clusters in apparently disallowed regions. Conformational space has been divided into 14 distinct bins. Residues adopting these relatively rare conformations are presented and amino acid propensities for these regions are estimated. Inspection of specific examples in a completely "arid", fully allowed region in the top left quadrant establishes that side-chain and backbone interactions may provide the energetic compensation necessary for populating this region of φ-ψ space. Asn, Asp, and His residues showed the highest propensities in this region. The two distinct clusters in the bottom right quadrant which are formally disallowed on strict steric considerations correspond to the gamma turn (C7 axial) conformation (Bin 12) and the i + 1 position of Type II' ß turns (Bin 13). Of the 516 non-Gly residues in Bin 13, 384 occupied the i + 1 position of Type II' ß turns. Further examination of these turn segments revealed a high propensity to occur at the N-terminus of helices and as a tight turn in ß hairpins. The ß strand-helix motif with the Type II' ß turn as a connecting element was also found in as many as 57 examples.

Molecular dynamics simulations of retinoblastoma protein.

Tumor suppressor proteins play a crucial role in cell cycle regulation. Retinoblastoma protein (pRB) is one among them which regulates G1-S transition by binding with transcription factors. The activity of pRB is deregulated by cyclin dependent kinases-mediated hyper-phosphorylation and also due to cancer-derived mutations. In addition, it is also deactivated by binding of viral onco-proteins such as large T antigen, E1A, and E7. These viral proteins initially recognize pRB through their conserved LxCxE motif and facilitate dissociation of preexisting pRB-E2F complex. Based on these features, molecular dynamics (MD) simulation is performed for four different states of pRB for which the crystal structure is available. The unliganded/apo form and complex forms with E2F and E7 peptides reveal the molecular mechanism behind the activation and inactivation of pRB. In addition, the ternary complex of pRB with both E7 and E2F (for which no crystal structure is available) is modeled and simulated to understand the influence of binding of one ligand on the other. The variations in the three major factors such as conformational changes, inter- and intra-molecular interactions, and binding free energies between the apo and complex forms confirm the possibility for designing a small molecule inhibitor to inhibit pRB-E7 interactions without altering the prebound E2F. The present study deals with the molecular modeling and MD simulations of pRB in free and ligand-bound forms and confirms that pRB could be a valid target for the anticancer drug design when the cancer is induced by the viral onco-proteins and forms a clear base for designing E7 antagonists.

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins.

Short range side chain-backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution ≤ 1.5 Å). Hydrogen bonds involving residues i - 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (C(n), where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C(10)(CO(δ)(i) …NH(i + 2)), C(13)(CO(δ)(i) …NH(i + 3)) and C(17)(N(δ)H(i) …CO(i - 4)) structures. In contrast, Gln predominantly forms C(16)(CO(ε)(i) …NH(i - 3)), C(12)(N(ε)H(i) …CO(i - 2)), C(15)(N(ε)H(i) …CO(i - 3)) and C(18)(N(ε)H(i) …CO(i - 4)) motifs, with only the C(18) motif being analogous to the Asn C(17) structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone ß-turns and α-turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain-backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered.

Conformational analysis and design of cross-strand disulfides in antiparallel ß-sheets.

Cross-strand disulfides bridge two cysteines in a registered pair of antiparallel ß-strands. A nonredundant data set comprising 5025 polypeptides containing 2311 disulfides was used to study cross-strand disulfides. Seventy-six cross-strand disulfides were found of which 75 and 1 occurred at non-hydrogen-bonded (NHB) and hydrogen-bonded (HB) registered pairs, respectively. Conformational analysis and modeling studies demonstrated that disulfide formation at HB pairs necessarily requires an extremely rare and positive χ¹ value for at least one of the cysteine residues. Disulfides at HB positions also have more unfavorable steric repulsion with the main chain. Thirteen pairs of disulfides were introduced in NHB and HB pairs in four model proteins: leucine binding protein (LBP), leucine, isoleucine, valine binding protein (LIVBP), maltose binding protein (MBP), and Top7. All mutants LIVBP T247C V331C showed disulfide formation either on purification, or on treatment with oxidants. Protein stability in both oxidized and reduced states of all mutants was measured. Relative to wild type, LBP and MBP mutants were destabilized with respect to chemical denaturation, although the sole exposed NHB LBP mutant showed an increase of 3.1°C in T(m). All Top7 mutants were characterized for stability through guanidinium thiocyanate chemical denaturation. Both exposed and two of the three buried NHB mutants were appreciably stabilized. All four HB Top7 mutants were destabilized (ΔΔG° = -3.3 to -6.7 kcal/mol). The data demonstrate that introduction of cross-strand disulfides at exposed NHB pairs is a robust method of improving protein stability. All four exposed Top7 disulfide mutants showed mild redox activity.

Molecular dynamics study of secretory phospholipase A2 of Russell's viper and bovine pancreatic sources.

A comparative molecular dynamics simulation of free and inhibitor-bound form of secretory phospholipase A(2) (sPLA(2)) of Russell's viper discloses the sort of restrictions in active site for inhibitor binding and implies suitable sites for further design of inhibitors based on active site scaffold. This enzyme belongs to group II PLA(2)s and dimerize asymmetrically with difference in orientation of W31 at the gateway of the active site of both the subunits. Hence, the active site of subunit A is open and that of subunit B is inaccessible to monodispersed inhibitors. PLA(2) enzymes are active at solvent-lipid interface and their action could be inhibited at the solvent environment before it reacts with aggregated substrates. Some sPLA(2)s, especially of different venom sources, undergo aggregation in a concentration-dependent manner, associate symmetrically into dimeric or trimeric form, and attain functional monomeric form during their interaction with the aggregated substrate. All sPLA(2)s exhibit catalysis with similar mechanism and show considerable differences in its way of inhibition. This necessitates conformational analysis on asymmetric dimer viper PLA(2) and its comparison with bovine pancreatic sPLA(2) (BPsPLA(2)) which belongs to group IB. BPsPLA(2) exists in monomeric form and does not have W31 at the gateway of hydrophobic pocket. In general, both monomeric and dimeric forms possess conserved active site with six subsites including the residues H48 and D49, and calcium-binding and surface loops. In the PLA(2) inhibitor complexes, the presence of calcium in monomer and W31 in dimer form is the unique feature and it makes the difference only in inhibitory mechanism without altering the catalytic mechanism. With this context, molecular dynamics simulation is performed for monomer and dimer form of sPLA(2)s in both native and complex forms. Comparison of trajectories with respect to fluctuation and deviation discloses the dynamics of surface and calcium-binding loops as well as the difference in dynamics of active site residues of group IB and II sPLA(2). Further, principal component and conformational cluster analyses are performed to substantiate the results.

Disulfide conformation and design at helix N-termini.

To understand structural and thermodynamic features of disulfides within an alpha-helix, a non-redundant dataset comprising of 5025 polypeptide chains containing 2311 disulfides was examined. Thirty-five examples were found of intrahelical disulfides involving a CXXC motif between the N-Cap and third helical positions. GLY and PRO were the most common amino acids at positions 1 and 2, respectively. The N-Cap residue for disulfide bonded CXXC motifs had average (phi,psi) values of (-112 +/- 25.2 degrees , 106 +/- 25.4 degrees ). To further explore conformational requirements for intrahelical disulfides, CYS pairs were introduced at positions N-Cap-3; 1,4; 7,10 in two helices of an Escherichia coli thioredoxin mutant lacking its active site disulfide (nSS Trx). In both helices, disulfides formed spontaneously during purification only at positions N-Cap-3. Mutant stabilities were characterized by chemical denaturation studies (in both oxidized and reduced states) and differential scanning calorimetry (oxidized state only). All oxidized as well as reduced mutants were destabilized relative to nSS Trx. All mutants were redox active, but showed decreased activity relative to wild-type thioredoxin. Such engineered disulfides can be used to probe helix start sites in proteins of unknown structure and to introduce redox activity into proteins. Conversely, a protein with CYS residues at positions N-Cap and 3 of an alpha-helix is likely to have redox activity.

Hybridization between threatened freshwater catfish Mystus gulio (Hamilton & buchanan) and Mystus montanus (Jerdon) by artificial fertilization.

Inter-specific hybrids were produced between the threatened catfish species Mystus gulio x Mystus montanus. The differences in percentage of fertilization and hatching between control and interspecies were significant. The survival of hybrid was significantly lower (24.80 +/- 4.3%) when compared to control (95.1 +/- 3.5%). Time difference in yolk absorption by hybrid (73.30 h) was higher than that of control (72 h). When compared to interspecific fertilized egg the hatching time (24-25 h) and viability of larvae of the control were significantly better. In hybrids more deformed hatchlings (52.7 +/- 4.2 %) were noticed than the control (24.80 +/- 4.3%).

Stereochemical criteria for prediction of the effects of proline mutations on protein stability.

When incorporated into a polypeptide chain, proline (Pro) differs from all other naturally occurring amino acid residues in two important respects. The phi dihedral angle of Pro is constrained to values close to -65 degrees and Pro lacks an amide hydrogen. Consequently, mutations which result in introduction of Pro can significantly affect protein stability. In the present work, we describe a procedure to accurately predict the effect of Pro introduction on protein thermodynamic stability. Seventy-seven of the 97 non-Pro amino acid residues in the model protein, CcdB, were individually mutated to Pro, and the in vivo activity of each mutant was characterized. A decision tree to classify the mutation as perturbing or nonperturbing was created by correlating stereochemical properties of mutants to activity data. The stereochemical properties including main chain dihedral angle phi and main chain amide H-bonds (hydrogen bonds) were determined from 3D models of the mutant proteins built using MODELLER. We assessed the performance of the decision tree on a large dataset of 163 single-site Pro mutations of T4 lysozyme, 74 nsSNPs, and 52 other Pro substitutions from the literature. The overall accuracy of this algorithm was found to be 81% in the case of CcdB, 77% in the case of lysozyme, 76% in the case of nsSNPs, and 71% in the case of other Pro substitution data. The accuracy of Pro scanning mutagenesis for secondary structure assignment was also assessed and found to be at best 69%. Our prediction procedure will be useful in annotating uncharacterized nsSNPs of disease-associated proteins and for protein engineering and design.

Structural compromise of disallowed conformations in peptide and protein structures.

Using a data set of 454 crystal structures of peptides and 80 crystal structures of non-homologous proteins solved at ultra high resolution of 1.2 A or better we have analyzed the occurrence of disallowed Ramachandran (phi, psi) angles. Out of 1492 and 13508 non-glycyl residues in peptides and proteins respectively 12 and 76 residues in the two datasets adopt clearly disallowed combinations of Ramachandran angles. These examples include a number of conformational points which are far away from any of the allowed regions in the Ramachandran map. According to the Ramachandran map a given (phi, psi) combination is considered disallowed when two non-bonded atoms in a system of two-linked peptide units with ideal geometry are prohibitively proximal in space. However, analysis of the disallowed conformations in peptide and protein structures reveals that none of the observations of disallowed conformations in the crystal structures correspond to a short contact between non-bonded atoms. A further analysis of deviations of bond lengths and angles, from the ideal peptide geometry, at the residue positions of disallowed conformations in the crystal structures suggest that individual bond lengths and angles are all within acceptable limits. Thus, it appears that the rare tolerance of disallowed conformations is possible by gentle and acceptable deviations in a number of bond lengths and angles, from ideal geometry, over a series of bonds resulting in a net gross effect of acceptable non-bonded inter-atomic distances.

Design of disulfide-linked thioredoxin dimers and multimers through analysis of crystal contacts.

Disulfide bonds play an important role in protein stability and function. Here, we describe a general procedure for generating disulfide-linked dimers and multimers of proteins of known crystal structures. An algorithm was developed to predict sites in a protein compatible with intermolecular disulfide formation with neighboring molecules in the crystal lattice. A database analysis was carried out on 46 PDB coordinates to verify the general applicability of this algorithm to predict intermolecular disulfide linkages. On the basis of the predictions from this algorithm, mutants were constructed and characterized for a model protein, thioredoxin. Of the five mutants, as predicted, in solution four formed disulfide-linked dimers while one formed polymers. Thermal and chemical denaturation studies on these mutant thioredoxins showed that three of the four dimeric mutants had similar stability to wild-type thioredoxin while one had lower stability. Three of the mutant dimers crystallized readily (in four to seven days) in contrast to the wild-type protein, which is particularly difficult to crystallize and takes more than a month to form diffraction-quality crystals. In two of the three cases, the structure of the dimer was exactly as predicted by the algorithm, while in the third case the relative orientation of the monomers in the dimer was different from the predicted one. This methodology can be used to enhance protein crystallizability, modulate the oligomerization state and to produce linear chains or ordered three-dimensional protein arrays.

Structural assessment of glycyl mutations in invariantly conserved motifs.

Motifs that are evolutionarily conserved in proteins are crucial to their structure and function. In one of our earlier studies, we demonstrated that the conserved motifs occurring invariantly across several organisms could act as structural determinants of the proteins. We observed the abundance of glycyl residues in these invariantly conserved motifs. The role of glycyl residues in highly conserved motifs has not been studied extensively. Thus, it would be interesting to examine the structural perturbations induced by mutation in these conserved glycyl sites. In this work, we selected a representative set of invariant signature (IS) peptides for which both the PDB structure and mutation information was available. We thoroughly analyzed the conformational features of the glycyl sites and their local interactions with the surrounding residues. Using Ramachandran angles, we showed that the glycyl residues occurring in these IS peptides, which have undergone mutation, occurred more often in the L-disallowed as compared with the L-allowed region of the Ramachandran plot. Short range contacts around the mutation site were analyzed to study the steric effects. With the results obtained from our analysis, we hypothesize that any change of activity arising because of such mutations must be attributed to the long-range interaction(s) of the new residue if the glycyl residue in the IS peptide occurred in the L-allowed region of the Ramachandran plot. However, the mutation of those conserved glycyl residues that occurred in the L-disallowed region of the Ramachandran plot might lead to an altered activity of the protein as a result of an altered conformation of the backbone in the immediate vicinity of the glycyl residue, in addition to long range effects arising from the long side chains of the new residue. Thus, the loss of activity because of mutation in the conserved glycyl site might either relate to long range interactions or to local perturbations around the site depending upon the conformational preference of the glycyl residue.

Thermodynamic effects of proline introduction on protein stability.

The amino acid Pro is more rigid than other naturally occurring amino acids and, in proteins, lacks an amide hydrogen. To understand the structural and thermodynamic effects of Pro substitutions, it was introduced at 13 different positions in four different proteins, leucine-isoleucine-valine binding protein, maltose binding protein, ribose binding protein, and thioredoxin. Three of the maltose binding protein mutants were characterized by X-ray crystallography to confirm that no structural changes had occurred upon mutation. In the remaining cases, fluorescence and CD spectroscopy were used to show the absence of structural change. Stabilities of wild type and mutant proteins were characterized by chemical denaturation at neutral pH and by differential scanning calorimetry as a function of pH. The mutants did not show enhanced stability with respect to chemical denaturation at room temperature. However, 6 of the 13 single mutants showed a small but significant increase in the free energy of thermal unfolding in the range of 0.3-2.4 kcal/mol, 2 mutants showed no change, and 5 were destabilized. In five of the six cases, the stabilization was because of reduced entropy of unfolding. However, the magnitude of the reduction in entropy of unfolding was typically several fold larger than the theoretical estimate of -4 cal K(-1) mol(-1) derived from the relative areas in the Ramachandran map accessible to Pro and Ala residues, respectively. Two double mutants were constructed. In both cases, the effects of the single mutations on the free energy of thermal unfolding were nonadditive.

Protein stabilization by introduction of cross-strand disulfides.

Disulfides cross-link residues in a protein that are separated in primary sequence and stabilize the protein through entropic destabilization of the unfolded state. While the removal of naturally occurring disulfides leads to protein destabilization, introduction of engineered disulfides does not always lead to significant stabilization of a protein. We have analyzed naturally occurring disulfides that span adjacent antiparallel strands of beta sheets (cross-strand disulfides). Cross-strand disulfides have recently been implicated as redox-based conformational switches in proteins such as gp120 and CD4. The propensity of these disulfides to act as conformational switches was postulated on the basis of the hypothesis that this class of disulfide is conformationally strained. In the present analysis, there was no evidence to suggest that cross-strand disulfides are more strained compared to other disulfides as assessed by their torsional energy. It was also observed that these disulfides occur solely at non-hydrogen-bonded (NHB) registered pairs of adjacent antiparallel strands and not at hydrogen-bonded (HB) positions as suggested previously. One of the half-cystines involved in cross-strand disulfide formation often occurs at an edge strand. Experimental confirmation of the stabilizing effects of such disulfides was carried out in Escherichia coli thioredoxin. Four pairs of cross-strand cysteines were introduced, two at HB and two at NHB pairs. Disulfides were formed in all four cases. However, as predicted from our analysis, disulfides at NHB positions resulted in an increase in melting temperature of 7-10 degrees C, while at HB positions there was a corresponding decrease of -7 degrees C. The reduced state of all proteins had similar stability.

Conformational analysis of invariant peptide sequences in bacterial genomes.

The functional significance of evolutionarily conserved motifs/patterns of short regions in proteins is well documented. Although a large number of sequences are conserved, only a small fraction of these are invariant across several organisms. Here, we have examined the structural features of the functionally important peptide sequences, which have been found invariant across diverse bacterial genera. Ramachandran angles (phi,psi) have been used to analyze the conformation, folding patterns and geometrical location (buried/exposed) of these invariant peptides in different crystal structures harboring these sequences. The analysis indicates that the peptides preferred a single conformation in different protein structures, with the exception of only a few longer peptides that exhibited some conformational variability. In addition, it is noticed that the variability of conformation occurs mainly due to flipping of peptide units about the virtual C(alpha)...C(alpha) bond. However, for a given invariant peptide, the folding patterns are found to be similar in almost all the cases. Over and above, such peptides are found to be buried in the protein core. Thus, we can safely conclude that these invariant peptides are structurally important for the proteins, since they acquire unique structures across different proteins and can act as structural determinants (SD) of the proteins. The location of these SD peptides on the protein chain indicated that most of them are clustered towards the N-terminal and middle region of the protein with the C-terminal region exhibiting low preference. Another feature that emerges out of this study is that some of these SD peptides can also play the roles of "fold boundaries" or "hinge nucleus" in the protein structure. The study indicates that these SD peptides may act as chain-reversal signatures, guiding the proteins to adopt appropriate folds. In some cases the invariant signature peptides may also act as folding nuclei (FN) of the proteins.