A method to detect nonproline cis peptide bonds in proteins

Detection of trans-cis flips and peptide-plane flips in protein structures

A coordinate-based method is presented to detect peptide bonds that need correction either by a peptide-plane flip or by a trans-cis inversion of the peptide bond. When applied to the whole Protein Data Bank, the method predicts 4617 trans-cis flips and many thousands of hitherto unknown peptide-plane flips. A few examples are highlighted for which a correction of the peptide-plane geometry leads to a correction of the understanding of the structure-function relation. All data, including 1088 manually validated cases, are freely available and the method is available from a web server, a web-service interface and through WHAT_CHECK.

Family 42 carbohydrate-binding modules display multiple arabinoxylan-binding interfaces presenting different ligand affinities

Enzymes that degrade plant cell wall polysaccharides display a modular architecture comprising a catalytic domain bound to one or more non-catalytic carbohydrate-binding modules (CBMs). CBMs display considerable variation in primary structure and are grouped into 59 sequence-based families organized in the Carbohydrate-Active enZYme (CAZy) database. Here we report the crystal structure of CtCBM42A together with the biochemical characterization of two other members of family 42 CBMs from Clostridium thermocellum. CtCBM42A, CtCBM42B and CtCBM42C bind specifically to the arabinose side-chains of arabinoxylans and arabinan, suggesting that various cellulosomal components are targeted to these regions of the plant cell wall. The structure of CtCBM42A displays a beta-trefoil fold, which comprises 3 sub-domains designated as alpha, beta and gamma. Each one of the three sub-domains presents a putative carbohydrate-binding pocket where an aspartate residue located in a central position dominates ligand recognition. Intriguingly, the gamma sub-domain of CtCBM42A is pivotal for arabinoxylan binding, while the concerted action of beta and gamma sub-domains of CtCBM42B and CtCBM42C is apparently required for ligand sequestration. Thus, this work reveals that the binding mechanism of CBM42 members is in contrast with that of homologous CBM13s where recognition of complex polysaccharides results from the cooperative action of three protein sub-domains presenting similar affinities.

Estimating the "steric clash" at cis peptide bonds

To account for the scarcity of cis peptide bonds in proteins, especially in nonproline (or secondary amide) cases, a steric-clash argument is often put forward, in a scheme where the R lateral chains are facing parallel one another, and the backbone is kept in an "all- trans"-like arrangement. Although such a steric conflict can be partly relieved through proper adjustment of the backbone dihedral angles, one can try to estimate its associated energy cost. To this end, quantum-chemistry approaches using a differential-torsion protocol and bond-separation-energy analyses are applied to N-ethyl propionamide CH3-CH2-CO-NH-CH2-CH3, regarded as a model capable of exhibiting C beta...C beta interaction as in alanine succession. The calculations provide an increment of 9 kcal/mol, quite close to that obtained in the nearly isostere (gsg) rotamer of n-hexane (10 kcal/mol), suggesting the local effects induced by methyl-methyl contact are similar in both cases. Analogous treatments on larger radicals as encountered in leucine or phenylalanine dimers do not change this increment much, which therefore defines the basic reference per-plaque quota to be overcome along all- cis chains. Explicit modeling indicated it can be reduced by up to a factor of 4.

A DOQSY approach for the elucidation of torsion angle distributions in biopolymers: application to silk

Silk from the wild silkworm Samia cynthia ricini with a molecular mass of about 300kDa consists of alternating repeats of nano-crystalline poly-(Ala) and non-crystalline glycine-rich domains. The backbone torsion angles between pairs of these two amino acids is determined by DOQSY solid-state NMR spectroscopy: the alanine-rich domains are predominantly in a beta-sheet conformation, whereas the glycine-rich domains are found to be partially in an extended beta-sheet conformation and partially in an approximately 3(1)-helical conformation. In the cast film from liquid silk significantly different secondary structures were found: the alanine-rich domains are alpha-helical conformation, whereas the results for glycine-labelled sample are explained by a random-coil state. A detailed error analysis of the technique is presented.

COPS--Cis/trans peptide bond conformation prediction of amino acids on the basis of secondary structure information

COPS predicts for all 20 naturally occurring amino acids whether the peptide bond in a protein is in cis or trans conformation. The algorithm is based only on secondary structure information of amino acid triplets without considering the amino acid sequence information. Conformation parameters are derived from solved 3D structures deposited in the PDB and led to propensities based on modified Chou-Fasman parameters. COPS analyses amino acid triplets taking only their respective secondary structure into consideration and upon application of a set of rules utilizing the conformation parameters, the N-terminal peptide bond conformation of the middle residue is predicted. COPS was tested on a random selection of protein datasets.

AVAILABILITY

The COPS program and further information are freely available from the FMP website at http://www.fmp-berlin.de/nmr/cops

CONTACT

labudde@fmp-berlin.de.

Computational study of the catalytic domain of human neutrophil collagenase. specific role of the S3 and S'3 subsites in the interaction with a phosphonate inhibitor

Human neutrophil collagenase (HNC, MMP-8) is one of the target enzymes for drug treatment of pathologic extracellular matrix degradation. Peptidomimetic inhibitors bind in the S'-side of the enzyme active site occupying the S'1 primary specificity pocket by their large hydrophobic side-chains. The crystal structure of the complex between the catalytic domain of MMP-8 and Pro-Leu-L-TrpP(OH)2 (PLTP) showed that this phosphonate inhibitor binds in the S side of the active site. This finding was unexpected since it represents the first example of accommodation of the bulky Trp indolyl chain in the S1 rather than in the S'1 subsite. Dynamical and structural factors favouring this uncommon mode of binding were therefore investigated. MD simulations performed on the uncomplexed enzyme show that its structure in aqueous solution is only slightly different from the crystal structure found in the complex with PLTP. ED analysis of the MD simulations, performed on PLTP alternatively interacting with the S- or S'-side of the active site, shows that the enzyme fluctuation increases in both cases. The main contribution to the overall enzyme fluctuation is given by the loop 164-173. The fluctuation of this loop is spread over more degrees of freedom when PLTP interacts with the S-side. This dynamical factor can enhance the preference of PLTP for the S subsites of MMP-8. MD simulations also show that ligation of PLTP in the S subsites is further favoured by better zinc chelation, a cation-pi interaction at the S3 subsite and unstrained binding conformations. The role of the S3, S'3 and S'1 subsites in determining the inhibitor binding is discussed.

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.

C-H...pi-interactions in proteins

A non-redundant set of 1154 protein structures from the Protein Data Bank was examined with respect to close interactions between C-H-donor and pi-acceptor groups. A total of 31,087 interactions were found to satisfy our selection criteria. Their geometric parameters suggest that these interactions can be classified as weak hydrogen bonds.A set of 12 interaction classes were defined based on the division of the donors into three groups and the acceptors into four groups. These classes were examined separately, and the respective interactions described in detail in each class. Most prominent were interactions between aliphatic C-H donors and aromatic pi-acceptors and interactions between aromatic C-H donors and aromatic pi-acceptors. About three-quarters of the Trp-rings, half of all Phe and Tyr-rings and a quarter of all His-rings were found to be involved as acceptors in C-H...pi-interactions. On the donor side, a preference for aromatic C-H groups was observed, but also for the aliphatic side-chains of the long, extended amino acid residues Lys, Arg and Met, and the Pro ring. The average distance between the C-donor and the center-of-mass of the pi-acceptor was observed to be significantly longer in the 174 protein structures determined at >2.5 A resolution. Also, the distribution is significantly wider. This resolution dependence suggests that the force fields commonly used for the refinement of protein structures may not be adequate. C-H...pi-interactions involving aromatic groups either as donor or as acceptor groups are found mostly in the interior of the protein. The more hydrophilic the participating groups are, the closer to the surface are the interactions located. About 40 % of all C-H...pi-interactions occur between amino acid residue side-chains that are separated by nine or less residues in sequence. Dependent on the interaction class, different preferences for secondary structure, residue type and side-chain conformation were observed. It is likely that the C-H...pi-interactions contribute significantly to the overall stability of a protein.