On the electrostatic and steric similarity of lactam compounds and the natural substrate for bacterial cell-wall biosynthesis

Molecular dynamics investigations on the effect of D amino acid substitution in a triple-helix structure and the stability of collagen

Studies on the structure and stability of peptides and proteins during l-->d configurational change are certainly important for the designing of peptides with new biological activity and protein engineering. The l-->d amino acid (d AA) changes have been observed in aged proteins such as collagen. Hence, in this study, an attempt has been made to explore the effect of the replacement of l amino acid (l AA) in the model collagen-like peptides with d AA and the origin of structural stability (destability) has been traced using the molecular dynamics (MD) method employing the AMBER force field. Our results reveal that the substitution of d AA produces a large local disruption to the triple-helical structure. Formation of a kink (bulge) at the site of substitution is observed from the detailed analysis of MD trajectory. However, this local perturbation of kinked helix changes the direction of the helices and affects the relative orientation of the respective AA residues for helix-helix interaction, enough to affect the overall stability of the model collagen-like peptide. The destabilization energy per d Ala substitution is 7.87 kcal/mol, which is similar to the value for the Gly-->Ala mutation in collagen. Since the Gly-->Ala mutation is involved in genetic disorders such as osteogenesis imperfecta (OI), the l-->d configurational change may produce a similar effect on collagen.

Is glycine a surrogate for a D-amino acid in the collagen triple helix?

Collagen is the most abundant protein in animals. Every third residue in a collagen strand is a glycine with phi, psi = -70 degrees, 175 degrees. A recent computational study suggested that replacing these glycine residues with D-alanine or D-serine would stabilize the collagen triple helix. This hypothesis is of substantial importance, as the glycine residues in collagen constitute nearly 10% of the amino acid residues in humans. To test this hypothesis, we synthesized a series of collagen mimic peptides that contain one or more D-alanine or D-serine residues replacing the canonical glycine residues. Circular dichroism spectroscopy and thermal denaturation experiments indicated clearly that the substitution of glycine with D-alanine or D-serine greatly disfavors the formation of a triple helix. Host-guest studies also revealed that replacing a single glycine residue with D-alanine is more destabilizing than is its replacement with L-alanine, a substitution that results from a common mutation in patients with collagen-related diseases. These data indicate that the glycine residues in collagen are not a surrogate for a D-amino acid and support the notion that the main-chain torsion angles of a glycine residue in the native structure (especially, phi > 0 degrees ) are critical determinants for its beneficial substitution with a D-amino acid in a protein.

D-amino acid residues in peptides and proteins

We have investigated the D-amino acid residues present in Protein Data Bank (PDB) entries, categorizing them into "real" D-residues and artifacts. In polypeptide chains of more than 20 residues, only a single instance of a "real" D-residue, other than those deliberately designed or engineered, was found. This example was the result of a slow chemical epimerization process. Another 12 designed D-residues were found in these longer polypeptide chains. Smaller peptides of 20 or fewer residues contained 479 "real" D-residues, the majority in various gramicidin, actinomycin, or cyclosporin structures. We found 148 PDB entries with "real" D-residues and a further 186, in which all apparent D-residues are artifacts. Investigating the (phi, psi) preferences of the "real" D-residues, we found that the region around (-60 degrees, -45 degrees ) was almost completely unoccupied, even though it is not formally disallowed. We link the low propensity to occupy this region with the alpha-helix destabilizing properties of D-residues.

Conformational analysis of bacterial cell wall peptides indicates how particular conformations have influenced the evolution of penicillin-binding proteins, beta-lactam antibiotics and antibiotic resistance mechanisms

Our aim was to use a conformational analysis technique developed for peptides to identify structural relationships between bacterial cell wall peptides and beta-lactam antibiotics that might help to explain their different actions as substrates and inhibitors of penicillin binding proteins (PBPs). The conformational forms of the model cell wall peptide Ac-L-Lys(Ac)-D-Ala-D-Ala are described by just a few backbone torsion combinations: three C-terminal carboxylate regions, with Tor8 (psi(i+1)) ranges of D3 region (50 degrees to 70 degrees ), D6 region (140 degrees to 170 degrees ) and D9 region (-50 degrees to -70 degrees ) are combined with either of two Tor6 (phi(i))-Tor4 (psi(i)) combinations, C4 region (-50 degrees to -80 degrees ) with B8 region (-40 degrees to -70 degrees ) or C11 region (30 degrees to 50 degrees ) with B2 region (30 degrees to 70 degrees ). From these results, and comparisons with conformational analyses of various beta-lactams and Ac-L-Lys(Ac)-D-Ala-D-Lac, it is concluded that molecular recognition of cell wall peptide substrates by PBPs requires conformers with backbone torsion angles of D3C4B8. beta-Lactam antibiotics are constrained compounds with fewer conformational forms; these match well the backbone torsions of cell wall peptides at D3C4, allowing their recognition and acylation by PBPs, whereas their unique Tor4 produces differently orientated CO and N atoms that appear to prevent subsequent deacylation, leading to their action as suicide substrates. The results are also related to the selective pressures involved in evolution of beta-lactamases from PBPs. From analysis of conformers of Ac-L-Lys(Ac)-D-Ala-D-Ala and the vancomycin-resistant analogue Ac-L-Lys(Ac)-D-Ala-D-Lac, it is concluded that vancomycin may recognise D6C11B2 conformers, giving it complementary substrate specificity to PBPs. This approach could have applications in the rational design of antibiotics targeted against PBPs and their substrates.

The orientation of N-H...O=C and N-H...N hydrogen bonds in biological systems: how good is a point charge as a model for a hydrogen bonding atom?

In order to design new ligands for protein-binding sites of unknown structure, it would be useful to predict the likely sites of hydrogen bonding of an unknown protein fragment to a known molecule. The positions of maxima and minima in the electrostatic potential at appropriate distances from the van der Waals surface were calculated for various small molecules, nucleic acid bases, peptide units and amino acid side chains containing groups which can form the biologically important N-H...O=C and N-H...N hydrogen bonds. Their ability to predict the positions of H and O/N in hydrogen bonded complexes, as predicted by optimising the electrostatic interactions of pairs of such molecules constrained by the molecular shapes, was assessed. It is shown that extrema in the electrostatic potential around the isolated molecules give worthwhile predictions for the locations of hydrogen binding partners. For molecules bound by a single N-H...O=C hydrogen bond, the electrostatic maximum associated with the H is usually less than 1 A from an acceptor atom, while a C=O electrostatic minimum is generally less than 1.5 A from the hydrogen bond proton. However, a significant number of hydrogen bonds form to the opposite lone pair from the electrostatic minimum, in which case the separation is up to 3.3 A. This reflects the broad electrostatic potential well around a carbonyl oxygen between the lone pair directions. The model predicts when neighbouring atoms drastically change the hydrogen bonding characteristics of an N-H or C=O group. Although the geometries of hydrogen bonded complexes are influenced by the other van der Waals contacts between the molecules, particularly multiple hydrogen bonds, these influences are constant when considering hydrogen bonding to a specific uncharacterised binding site. Hence, the consideration of sterically accessible electrostatic extrema will be useful in the design of new ligands.