Investigations of Ramachandran disallowed conformations in protein domain families

In peptide and protein structures, occurrence of (phi,psi.) 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 (A) over circle). 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 (A) over circle or better, we have analyzed the conservation of 221 disallowed conformations. It is observed that the disallowed conformation is only moderately conservedin 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. (C) 2013 Elsevier B.V. All rights reserved.

PocketMatch (version 2.0): A parallel algorithm for the detection of structural similarities between protein ligand binding-sites

Knowledge of protein-ligand interactions is essential to understand several biological processes and important for applications ranging from understanding protein function to drug discovery and protein engineering. Here, we describe an algorithm for the comparison of three-dimensional ligand-binding sites in protein structures. A previously described algorithm, PocketMatch (version 1.0) is optimised, expanded, and MPI-enabled for parallel execution. PocketMatch (version 2.0) rapidly quantifies binding-site similarity based on structural descriptors such as residue nature and interatomic distances. Atomic-scale alignments may also be obtained from amino acid residue pairings generated. It allows an end-user to compute database-wide, all-to-all comparisons in a matter of hours. The use of our algorithm on a sample dataset, performance-analysis, and annotated source code is also included.

A model for the interaction of trifluoroethanol with peptides and proteins

The structural stabilizing property of 2,2,2-trifluoroethanol (TFE) in peptides has been widely demonstrated, More recently, TFE has been shown to enhance secondary structure content in globular proteins, and to influence quaternary interactions in protein multimers. The molecular mechanisms by which TFE exerts its Influence on peptide and protein structures remain poorly understood. The present analysis integrates the known physical properties of TFE with a variety of experimental observations on the interaction of TFE with peptides and proteins and on the properties of fluorocarbons. Two features of TFE, namely the hydrophobicity of the trifluoromethyl group and the hydrogen bonding character (strong donor and poor acceptor), emerge as the most important factors for rationalising the observed effects of TFE. A model is proposed for TFE interaction with peptides which involves an initial replacement of the hydration shell by fluoroalcohol molecules, a process driven by apolar interactions and favourable entropy of dehydration. Subsequent bifurcated hydrogen-bond formation with peptide carbonyl groups, which leave intramolecular interactions unaffected, promotes secondary structure formation.

Hydrogen bonding in peptide helices - Analysis of two independent helices in the crystal structure of a peptide Boc-Val-Ala-Leu-Aib-Val-Ala-Phe-OMe

The crystal structure determination of the heptapeptide Boc-Val-Ala-Leu-Aib-Val-Ala-Phe-OMe reveals two peptide helices in the asymmetric unit, Crystal parameters are: space group P2(1), a = 10.356(2) Angstrom, b = 19.488(5) Angstrom, c = 23.756(6) Angstrom, beta = 102.25(2)degrees), V = 4685.4 Angstrom(3), Z = 4 and R = 5.7% for 7615 reflections [I>3 sigma(I)]. Both molecules adopt largely alpha-helical conformations with variations at the C-terminus, Helix type Is determined by analysing both 4-->1 and 5-->1 hydrogen-bond interactions and comparison with the results of analysis of protein structures. The presence of two 4-->1 hydrogen-bond interactions, besides four 5-->1 interact ions in both the conformations provides an opportunity to characterize bifurcated hydrogen bonds at high resolution, Comparison of the two helical conformations with related peptide structures suggests that distortions at the C-terminus are more facile than at the N-terminus.

Identification and analysis of extended strands and beta-sheets in globular proteins

A method to identify β-sheets in globular proteins from extended strands, using only α-carbon positions, has been developed. The strands that form β-sheets are picked up by means of simple distance criteria. The method has been tested by applying it to three proteins with accurately known secondary structures. It has also been applied to ten other proteins wherein only α-carbon coordinates are available, and the list of β-sheets obtained. The following points are worth noting: (i) The sheets identified by the algorithm are found to agree satisfactorily with the reported ones based on backbone hydrogen bonding, wherever this information is available. (ii) β-Strands that do not form parts of any sheet are a common feature of protein structures. (iii) Such isolated β-strands tend to be short. (iv) The conformation corresponding to the preferred right-handed twist of the sheet is overwhelmingly observed in both the sheet-forming and isolated β-strands.

BSDD-2.0: Biomolecules Segment Display Device - An updated web accessible computing resource

The existing internet computing resource, Biomolecules Segment Display Device (BSDD), has been updated with several additional useful features. An advanced option is provided to superpose the structural motifs obtained from a search on the Protein Data Bank (PDB) in order to see if the three-dimensional structures adopted by identical or similar sequence motifs are the same. Furthermore, the options to display structural aspects like inter- and intra-molecular interactions, ion-pairs, disulphide bonds, etc. have been provided.The updated resource is interfaced with an up-to-date copy of the public domain PDB as well as 25 and 90% non-redundant protein structures. Further, users can upload the three-dimensional atomic coordinates (PDB format) from the client machine. A free molecular graphics program, JMol, is interfaced with it to display the three-dimensional structures.

The nonplanar peptide unit III. Quantum chemical calculations for related compounds and experimental X-ray diffraction data

The possible nonplanar distortions of the amide group in formamide, acetamide, N-methylacetamide, and N-ethylacetamide have been examined using CNDO/2 and INDO methods. The predictions from these methods are compared with the results obtained from X-ray and neutron diffraction studies on crystals of small open peptides, cyclic peptides, and amides. It is shown that the INDO results are in good agreement with observations, and that the dihedral angles N and defining the nonplanarity of the amide unit are correlated approximately by the relation N = -2, while C is small and uncorrelated with . The present study indicates that the nonplanar distortions at the nitrogen atom of the peptide unit may have to be taken into consideration, in addition to the variation in the dihedral angles (,), in working out polypeptide and protein structures.

Parallel implementation of AutoDock

Computational docking of ligands to protein structures is a key step in structure-based drug design. Currently, the time required for each docking run is high and thus limits the use of docking in a high-throughput manner, warranting parallelization of docking algorithms. AutoDock, a widely used tool, has been chosen for parallelization. Near-linear increases in speed were observed with 96 processors, reducing the time required for docking ligands to HIV-protease from 81 min, as an example, on a single IBM Power-5 processor ( 1.65 GHz), to about 1 min on an IBM cluster, with 96 such processors. This implementation would make it feasible to perform virtual ligand screening using AutoDock.

Conformations of disulfide bridges in proteins

The conformational characteristics of disulfide bridges in proteins have been analyzed using a dataset of 22 protein structures, available at a resolution of 2.0 Å, containing a total of 72 disulfide crosslinks. The parameters used in the analysis include (φ, Ψ) values at Cys residues, bridge dihedral angles χss, χ1i, χ1j, χ2i and χ2j the distances Cαi-Cαj and Cβi-Cβj between the Cα and Cβ atoms of Cys(i) and Cys(j). Eight families of bridge conformations with three or more occurrences have been identified on the basis of these stereochemical parameters. The most populated family corresponds to the "left handed spiral" identified earlier by Richardson ((1981) Adv. Protein Chem. 34, 167–330). Disulfide bridging across antiparallel extended strands is observed in α-lytic protease, crambin, and β-trypsin and this structure is shown to be very similar to those obtained in small cystine peptides. Solvent accessible surface area calculations show that the overwhelming majority of disulfide bridges are inaccessible to solvent.

Recognition of active and inactive catalytic triads: A template based approach

It Is well established that a sequence template along with the database is a powerful tool for identifying the biological function of proteins. Here, we describe a method for predicting the catalytic nature of certain proteins among the several protein structures deposited in the Protein Data Bank (PDB) For the present study, we considered a catalytic triad template (Ser-His-Asp) found in serine proteases We found that a geometrically optimized active site template can be used as a highly selective tool for differentiating an active protein among several inactive proteins, based on their Ser-His-Asp interactions. For any protein to be proteolytic in nature, the bond angle between Ser O-gamma-Ser H-gamma His N-epsilon 2 in the catalytic triad needs to be between 115 degrees and 140 degrees The hydrogen bond distance between Ser H-gamma His N-epsilon 2 is more flexible in nature and it varies from 2 0 angstrom to 27 angstrom while in the case of His H-delta 1 Asp O-delta 1, it is from 1.6 angstrom to 2.0 angstrom In terms of solvent accessibility, most of the active proteins lie in the range of 10-16 angstrom(2), which enables easy accessibility to the substrate These observations hold good for most catalytic triads and they can be employed to predict proteolytic nature of these catalytic triads (C) 2010 Elsevier B V All rights reserved.

Proteins in Wastewater and Wastewater Sludges

Although there is some information on the total amounts of proteins in wastewater and sludges1"3 and on the amino acids in them,4-11 especially in activated sludge,12 there is almost no evidence on the nature of the proteins in these materials. A knowledge of the nature of proteins in wastewater, sludges, and similar substances would be useful not only for determining the pollutional effects on the environment and the changes in the protein structures during decomposition or treatment, but also for determining the possible usage of the resulting materials in agriculture,13 includ ing animal nutrition.

A theory on the migration of an extraneous electron across hydrogen bonds in polypeptides

The migrating electrons in biological systems normally are extraneous and taking this into account the electron delocalisation across the hydrogen bonds in proteins is re-examined. It is seen that an extraneous electron can travel rapidly via the low-lying virtual orbitals of the hydrogen-bonded π-electronic structure of peptide units in proteins. The frequency of electron transfer decreases slowly with an increase in the path length. However, the coupling of electron and protonic motions enhances this frequency. Transfer of electrons across the hydrogen bonds in accordance with the double-exchange mechanism does not appear to be possible. This theory offers a possibility for an extraneous electron to transfer within protein structures.

A C-H…O hydrogen bond stabilized polypeptide chain reversal motif at the C-terminus helices in proteins.

The serendipitous observation of a C–Hcdots, three dots, centeredO hydrogen bond mediated polypeptide chain reversal in synthetic peptide helices has led to a search for the occurrence of a similar motif in protein structures. From a dataset of 634 proteins, 1304 helices terminating in a Schellman motif have been examined. The C–Hcdots, three dots, centeredO interaction between the T−4 CαH and T+1 C=O group (Ccdots, three dots, centeredO≤3.5 Å) becomes possible only when the T+1 residue adopts an extended β conformation (T is defined as the helix terminating residue adopting an αL conformation). In all, 111 examples of this chain reversal motif have been identified and the compositional and conformational preferences at positions T−4, T, and T+1 determined. A marked preference for residues like Ser, Glu and Gln is observed at T−4 position with the motif being further stabilized by the formation of a side-chain–backbone Ocdots, three dots, centeredH–N hydrogen bond involving the side-chain of residue T−4 and the N–H group of residue T+3. In as many as 57 examples, the segment following the helix was extended with three to four successive residues in β conformation. In a majority of these cases, the succeeding β strand lies approximately antiparallel with the helix, suggesting that the backbone C–Hcdots, three dots, centeredO interactions may provide a means of registering helices and strands in an antiparallel orientation. Two examples were identified in which extended registry was detected with two sets of C–Hcdots, three dots, centeredO hydrogen bonds between (T−4) CαHcdots, three dots, centeredC=O (T+1) and (T−8) CαHcdots, three dots, centeredC=O (T+3).

Empirical torsional potential functions from protiesn structure DATA Phi- and Psi-Potentials for Non-glycyl Amino Acid Residues

The torsional potential functions Vt(φ) and Vt(ψ) around single bonds N–Cα and Cα-C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (φ, ψ)-plane with the value of Vtot(φ, ψ), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in l-configuration, are Vt(φ) = – 1.0 cos (φ + 60°); Vt(ψ) = – 0.5 cos (ψ + 60°) – 1.0 cos (2ψ + 30°) – 0.5 cos (3ψ + 30°). The dipeptide energy maps Vtot(φ, ψ) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line ψ = 0°. These functions, derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembmne helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of ? helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge oil the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. (C) 1996 John Wiley & Sons, Inc.