Predominant torsional forms adopted by dipeptide conformers in solution: parameters for molecular recognition.

The present paper describes the predominant conformational forms adopted by dipeptides in aqueous solution. More than 50 dipeptides were subjected to conformational analysis using SYBYL Random Search. The resultant collections of conformers for individual dipeptides, for small groups with related side chain residues and for large groups of about 50 dipeptides were visualized graphically and analysed using a novel three-dimensional pseudo-Ramachandran plot. The distribution of conformers, weighted according to the percentage of each in the total conformer pool, was found to be restricted to nine main combinations of backbone psi (psi) and phi (phi) torsion angles. The preferred psi values were in sectors A7 (+150 degrees to +/-180 degrees), A10 (+60 degrees to +90 degrees) and A4 (-60 degrees to -90 degrees), and these were combined with preferred phi values in sectors B12 (-150 degrees to +/-180 degrees), B9 (-60 degrees to -90 degrees) and B2 (+30 degrees to +60 degrees). These combinations of psi and phi values are distinct from those found in common secondary structures of proteins. These results show that although dipeptides can each adopt many conformations in solution, each possesses a profile of common conformers that is quantifiable. A similarly weighted distribution of dipeptide conformers according to distance between amino-terminal nitrogen and carboxyl-terminal carbon shows how the preferred combinations of backbone torsional angles result in particular N-C geometries for the conformers. This approach gives insight into the important conformational parameters of dipeptides that provide the basis for their molecular recognition as substrates by widely distributed peptide transporters. It offers a basis for the rational design of peptide-based bioactive compounds able to exploit these transporters for targeting and delivery.     

Polynucleotide folding in yeast tRNAPhe: elucidation of short-, medium-, and long-range interactions of sugar-phosphate-sugar backbone and base using a "blocked" nucleotide probe

It has recently been proposed that the repeating backbone nucleotide may be regarded as consisting of two blocks of equal magnitude representable by two virtual bonds. Implicit consideration of the nucleotide (ψ,ψ) and internucleotide (ω′,ω) geometry that generate variety in polynucleotide conformations, and of the constancy of the repeating structural moieties (P-C4′ and C4′-P) independent of the above rotations, has enabled us to utilize this scheme in the study of ordered structures such as di-, oligonucleotides and, most significantly, tRNA. The polynucleotide folding dictated by short-, intermediate-, and long-range interactions in the monoclinic and orthorhombic forms is described and compared through circular plots depicting the virtual bond torsions and distance plots constructed independently for backbone as well as bases. The torsions and the bond angles associated with the virtual bonds afford a clear distinction between ordered helical segments from loops and bends of tRNA. Lower virtual bond torsions (?60° to 60°) concomitant with higher values of virtual bond angles characterize various bend regions, while torsions around 160°–210° typify ordered helical strands. The distance plot elucidates the type of interaction associated with various sub-structures (helix–helix, helix–loop, and loop–loop) that form the constituents of different structural domains. Several other features such as the manifestation of the P₁₀ loop and the approximate twofold symmetry in the tRNA molecule are conspicuous on the distance plot.     

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.     

A comprehensive library of blocked dipeptides reveals intrinsic backbone conformational propensities of unfolded proteins

Despite prolonged scientific efforts to elucidate the intrinsic peptide backbone preferences of amino-acids based on understanding of intermolecular forces, many open questions remain, particularly concerning neighboring peptide interaction effects on the backbone conformational distribution of short peptides and unfolded proteins. Here, we show that spectroscopic studies of a complete library of 400 dipeptides reveal that, irrespective of side-chain properties, the backbone conformation distribution is narrow and they adopt polyproline II and β-strand, indicating the importance of backbone peptide solvation and electronic effects. By directly comparing the dipeptide circular dichroism and NMR results with those of unfolded proteins, the comprehensive dipeptides form a complete set of structural motifs of unfolded proteins. We thus anticipate that the present dipeptide library with spectroscopic data can serve as a useful database for understanding the nature of unfolded protein structures and for further refinements of molecular mechanical parameters.     

Conformational limitations of glycylsarcosine as a prototypic substrate for peptide transporters

Peptide transporters are present in all species to absorb the small peptides that occur ubiquitously as products of proteolysis. The broad substrate specificities of these systems allow them to be exploited therapeutically for delivery of peptidomimetic drugs in microbes and man. To this end, glycylsarcosine is currently used as a standard substrate for assaying peptidomimetic transport by peptide transporters. However, in this study we find it is unsuitable as a general substrate, based on assays of its transport by model bacterial peptide transporters and computer-based conformational analysis of its structure. Of the two generic transporters for di- and tripeptides, exemplified by Dpp and Tpp in Escherichia coli, only Dpp can transport glycylsarcosine. The explanation for this finding came from molecular modelling, which indicated that glycylsarcosine can adopt only a restricted range of conformers compared with typical dipeptides, and that of the conformers with a trans peptide bond, the majority have the specific psi and phi backbone torsion angles needed for molecular recognition and transport by Dpp but none possessed psi and phi torsions required for recognition by Tpp; moreover, 38% of its conformers have cis peptide bonds that are not substrates for any peptide transporter. Thus, using glycylsarcosine as substrate in competition assays with compounds that typically form conformers recognised by both types of peptide transporter will underestimate their transport. These findings have implications for assays of oral availability of peptidomimetic drugs such as beta-lactams, ACE inhibitors and anti-viral compounds, for which glycylsarcosine is routinely used.     

Molecular dynamics simulations of the ligand-binding domain of an N-methyl-D-aspartate receptor

The mechanism of partial agonism at N-methyl-D-aspartate receptors is an unresolved issue, especially with respect to the role of protein dynamics. We have performed multiple molecular dynamics simulations (7 x 20 ns) to examine the behavior of the ligand-binding core of the NR1 subunit with a series of ligands. Our results show that water plays an important role in stabilizing different conformations of the core and how a closed cleft conformation of the protein might be stabilized in the absence of ligands. In the case of ligand-bound simulations with both full and partial agonists, we observed that ligands within the binding cleft may undergo distinct conformational changes, without grossly influencing the degree of cleft closure within the ligand-binding domain. In agreement with recently published crystallographic data, we also observe similar changes in backbone torsions corresponding to the hinge region between the two lobes for the partial agonist, D-cycloserine. This observation rationalizes the classification of D-cycloserine as a partial agonist and should provide a basis with which to predict partial agonism in this class of receptor by analyzing the behavior of these torsions with other potential ligands.     

Cyclic retro-inverso dipeptides with two aromatic side chains. II. Conformational analysis.

The conformations of cis and trans cyclic retro-inverso dipeptides--2-[(4-hydroxy)benzyl]-5-benzyl-4,6(1H,2H,3H,5H)-pyrimidinedi one (c[mTyr-gPhe]), and 2-benzyl-5-amino-5-[(4-hydroxy)benzyl]-4,6(1H,2H,3H,5H)-pyrimidinedione (c[mTyr-gPhe]), and 2-benzyl-5-amino-5-[(4-hydroxy)benzyl]-4,6(1H,2H,3H,5H)-pyrimidinedione (c[(alpha-amino)mTyr-gPhe])--and the parent cyclic dipeptides--c[tyrosyl-phenylalanine] (cis-c[L-Tyr-L-Phe]) and c[tyrosyl-D-phenylalanine] (trans-c[L-Tyr-D-Phe])--were studied by using 1H-nmr spectroscopy and semiempirical energy calculations. In the cis compounds of all the cyclic retro-inverso and parent dipeptides, the most stable conformer has both aromatic side chains sharing the space over the backbone ring in a "face-to-face" fashion. All the trans compounds predominantly assume a "sandwich" conformation in which the two aromatic rings are folded back over the backbone ring on opposite sides. However, different conformational preferences were observed for the backbones between the retro-inverso and parent cyclic dipeptides. The parent cyclic dipeptide trans-c[L-Tyr-D-Phe] adopts two types of boat structures with different side-chain orientations in almost equal amounts: one with the Tyr side chain in a pseudoaxial position and the Phe side chain in a pseudoequatorial position, the other with the Tyr side chain in a pseudoequatorial position and the Phe side chain in a pseudoaxial position. On the other hand, the cyclic retro-inverso dipeptides trans-c[mPhe-gTyr] and trans c[mTyr-gPhe] assume only one type of boat structure in which the malonyl side chain is in a pseudoequatorial and the gem-diamino side chain is in a pseudoaxial position. In addition to the preferred conformations, the conformational energies of the C alpha--C beta bonds in the malonyl and gem-diamino residues were estimated from the temperature variation of vicinal 1H--1H coupling constants for the H--C alpha--C beta--H groupings observed for the trans isomers of cyclic retro-inverso dipeptides. The energies were evaluated to be 1.1 and 1.8 kcal mol-1 for the malonyl and gem-diamino residues, respectively. Applying these energies to the parent cyclic dipeptide trans-c[L-Tyr-D-Phe], the observed fractions of three side-chain conformations are reasonably reproduced. The conformational energies as well as conformational properties of the molecules estimated in this investigation may be useful to refine force constants for both parent and retro-inverso peptides with aromatic side chains.     

Conformational characteristics of mixed sugar puckered deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp from energy minimization studies

The deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp were subjected to a rigorous theoretical investigation with a view to describing their distinctive conformational characteristics. For each unit 216 probable three-dimensional forms defined by the backbone-base dihedral angles and sugar pucker modes were considered for conformational energy minimization process and scrutinized with reference to properties, such as base-stacking, hydrogen-bonding, internal flexibility and base sequence-phosphate influence. The P-O bond torsions and the phosphate groups were treated with special attention. The results reveal a number of preferred conformational states other than the known helical forms, such as, A-, B-, C-, Z-, and Watson-Crick conformation. Many interesting one-step (change in only one of the dihedral angles or sugar puckers) conformational transitions which involve just about a kcal/mol of energy came to light. The two base sequences CG and GC were noted to differ strikingly in many of their conformational characteristics.     

Predictions of protein backbone bond distances and angles from first principles

A procedure is described, based on a spline-function representation of ab initio peptide conformational geometry maps, that allows one to predict backbone bond distances and angles of proteins as functions of the peptide ?(N-C^α)/Ψ(C^α-C′) torsions with an accuracy comparable to that of high-resolution protein crystallography. For example, for the more than 40 residues of crambin, the rms deviation between predicted and crystallographic values of N-C^α-C′ is 1.9° for the 1.5 Å resolution structure and 1.9° for the 0.83 Å resolution structure, compared with angle variations of < 10°. Accurate information on protein backbone geometries is important for establishing dictionaries of flexible geometry functions for use in empirical peptide and protein modeling. © 1995 John Wiley & Sons, Inc.     

Molecular-dynamics simulations of [Met5]- and [D-Ala2,Met5]-enkephalins. Biological implication of monomeric folded and dimeric unfolded conformations.

To investigate the biologically active conformation of enkephalin, molecular-dynamics simulations were applied to [Met5]- and [D-Ala2,Met5]-enkephalins. The dynamic trajectory of monomeric extended [Met5]-enkephalin was analysed in terms of relative mobility between respective torsions of backbone chain. After 10 ps of the dynamics simulation, the conformational transition was converged into a stationary state among the beta-bend folded forms, where they are stabilized by several intramolecular hydrogen-bond formations. Similar conformational transition was also observed in the dynamics simulation of [D-Ala2,Met5]enkephalin, which is a more mu-receptor-specific peptide than [Met5]enkephalin. The geometrical correspondence between the monomeric enkephalin conformation in the stationary state and morphine molecule (a mu-specific rigid opiate) was surveyed by virtue of the triangular substructures generated by choosing three functional atoms in each molecule, and good resemblances were observed. On the other hand, the dynamics simulation of the antiparallel extended [Met5]enkephalin dimer showed a trajectory different from that of the monomeric one. Two intermolecular hydrogen bonds at Tyr1 (NH3+)...Met5(CO2-) end residues were held throughout the 100 ps simulation, the dimeric structure being consequently kept. The conformational transition of the backbone chains from the antiparallel extended form to the twisted one took place via an intermediate state. Many conformations revealed during the dynamics simulation showed that the relative orientations of each two Tyr1, Gly3, Phe4 and Met5 residues in the dimer are nearly related by a pseudo-C2-symmetry respectively, and both halves of the dimer structure could be further fitted to the monomeric folded enkephalin conformation. The monomeric and dimeric conformations of enkephalin at their stationary states are discussed in relation to the substrate-specificity for mu- and delta-opioid receptors.     

Conformational Characteristics of Mixed Sugar Puckered Deoxydinucleoside Triphosphate Units d-pCpGp and d-pGpCp from Energy Minimization Studies

Abstract

The deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp were subjected to a rigorous theoretical investigation with a view to describing their distinctive conformational characteristics. For each unit 216 probable three-dimensional forms defined by the backbone-base dihedral angles and sugar pucker modes were considered for conformational energy minimization process and scrutinized with reference to properties, such as base-stacking, hydrogen-bonding, internal flexibility and base sequence-phosphate influence. The P-O bond torsions and the phosphate groups were treated with special attention. The results reveal a number of preferred conformational states other than the known helical forms, such as, A-, B-, C-, Z-, and Watson-Crick conformation. Many interesting one-step (change in only one of the dihedral angles or sugar puckers) conformational transitions which involve just about a kcal/mol of energy came to light. The two base sequences CG and GC were noted to differ strikingly in many of their conformational characteristics.     

Structuro-functional organization of delta sleep-inducing peptide

Theoretical conformational analysis was carried out for a nonapeptide hormone (delta sleep-inducing peptide). Possible structure of the neuropeptide under physiological conditions may be described by a set of low-energy conformations belonging to nine different forms of the backbone. A solution of the "reverse conformational problem" for delta sleep inducing peptide enables one to predict modified amino acid sequences (D-Ala3-, Pro4-, Pro6-, Pro7, and Tyr7-analogs), which may assume one of the low-energy states of the native hormone. The influence of the solute was not taken into account in our calculations.     

HU protein employs similar mechanisms of minor-groove recognition in binding to different B-DNA sites: demonstration by Raman spectroscopy

The sequence isomers d(CGCAAATTTGCG) and d(TCAAGGCCTTGA) form self-complementary duplexes that present distinct targets for binding of the homodimeric architectural protein HU of Bacillus stearothermophilus (HUBst). Raman spectroscopy shows that although each duplex structure is of the B-DNA type, there are subtle conformational dissimilarities between them, involving torsion angles of the phosphodiester backbone and the arrangements of stacked bases. Each DNA duplex forms a stable stoichiometric (1:1) complex with HUBst, in which the structure of the HUBst dimer is largely conserved. However, the Raman signature of each DNA duplex is perturbed significantly and similarly with HUBst binding, as reflected in marker bands assigned to localized vibrations of the phosphodiester moieties and base residues. The spectral perturbations identify a reorganization of the DNA backbone and partial unstacking of bases with HUBst binding, which is consistent with non-sequence-specific minor-groove recognition. Prominent among the HUBst-induced perturbations of B-DNA are a conversion of approximately one-third of the alpha/beta/gamma torsions from the canonical g(-)/t/g(+) conformation to an alternative conformation, an equivalent conversion of deoxyadenosyl moieties from the C2'-endo/anti to the C3'-endo/anti conformation, and appreciable unstacking of purines. The results imply that each solution complex is characterized by structural perturbations extending throughout the 12-bp sequence. Comparison with previously studied protein/DNA complexes suggests that binding of HUBst bends DNA by approximately 70 degrees.     

Principles of carbopeptoid folding: a molecular dynamics simulation study.

The conformational spaces of five oligomers of tetrahydrofuran-based carbopeptoids in chloroform and dimethyl sulfoxide were investigated through nine molecular dynamics simulations. Prompted by nuclear magnetic resonance experiments that indicated various stable folds for some but not all of these carbopeptoids, their folding behaviour was investigated as a function of stereochemistry, chain length and solvent. The conformational distributions of these molecules were analysed in terms of occurrence of hydrogen bonds, backbone torsional-angle distributions, conformational clustering and solute configurational entropy. While a cis-linkage across the tetrahydrofuran ring favours right-handed helical structures, a trans-linkage results in a larger conformational variability. Intra-solute hydrogen bonding is reduced with increasing chain length and with increasing solvent polarity. Solute configurational entropies confirm the picture obtained: they are smaller for cis- than for trans-linked peptides, for chloroform than for dimethyl sulfoxide as solvent and for shorter peptide chains. The simulations provide an atomic picture of molecular conformational variability that is consistent with the available experimental data.     

Conformational energy studies of linear dipeptides H-X-L-Pro-OH

Empirical conformational energy calculations with the use of ECEPP energy functions have been carried out for linear dipeptides H-X-L -Pro-OH, with X = Gly, L -Ala, D -Ala, L -Leu, D -Leu, L -Phe, and D -Phe, in different states of protonation of the end groups. The results of these calculations are compared with the previously reported experimental equilibrium populations for the cis and trans isomers of the X-Pro bond in the different species. For all the protonation states of the seven dipeptides, the calculated nonbonded interactions and the conformational entropy term lead to a preference of the trans forms over the cis isomers by at least 1 kcal/mol. The electrostatic interactions stabilize the cis conformations in all species except the cationic forms of the D ,L -peptides, and it could further be shown that only the carbonyl group of X and the two end groups contribute significantly to the total electrostatic energy. One of the principal results of the experimental studies, i.e., the occurrence of 5–15% cis-proline in all the peptides with an uncharged C-terminus, was corroborated by our investigation of the cationic species. A detailed assessment of the electrostatic contribution to the total energy of the different conformations of H-Gly-L -Pro-OH indicates that the standard ECEPP parameters tend to overestimate the electrostatic interactions in aqueous solutions of the X-Pro dipeptides.     

How flexible are DNA constituents? The quantum-mechanical study

Relaxed force constants (RFCs) and vibrational root-mean-square deviations have been evaluated by the original calculation method for conformational parameters of the DNA structural units and their constituents: nucleic acid bases (uracile, thymine, cytosine, adenine and guanine) and their 'building blocks' (benzene, pyrimidine, imidazole and purine molecules), as well as the DNA backbone structural units: tetrahydrofuran, 1,2-dideoxyribose, methanol and orthophosphoric acid. It has been found that the RFCs for nomenclature torsions beta, gamma, epsilon; and sugar pseudorotation angle P in 1,2-dideoxyribose are sensible to the molecule conformation and their values are in the range of 1-25 kcal/(mole·rad2) obeying the inequality K(γ)> K(ε) > K(ρ) > K(β). The RFCs values for endocyclic torsions of nucleic acid bases six-member rings lie within 15-45?kcal/(mole·rad2) in pyrimidines and within 20-60?kcal/(mole·rad2) in purines. It is shown that the quantum zero-point motion effectively neglects the amino group non-planarity in cytosine, adenine and partially in guanine.     

Nucleic acid model building: the multiple backbone solutions associated with a given base morphology

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (chi) and sugar torsions (both the ring puckering and the exocyclic C5'-C4' (psi) torsion) are treated as independent variables and the resulting O3'...O5' distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles--phi', omega', omega and phi--about the C3'-O3'-P-O5'-C5' bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of the local sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.     

Theoretical studies of DNA-RNA hybrid conformations.

Molecular modelling has been used to probe the conformational preferences of double stranded DNA-RNA hybrids. As might be expected, the sugars of the DNA strand have higher conformational flexibility, but, for the majority of the repetitive sequences studied, these sugars prefer a C2-endo pucker, while ribose sugars uniformly adopt a C3-endo pucker. This gives rise to a strongly heteronomous duplex conformation. One exception to this rule involves the thymidine strand of poly(dT).poly(rA), which marginally prefers a C3-endo pucker. Our study further indicates that the DNA strands of the hybrids favour backbone torsions in the canonical B domain, rather than the modified values proposed on the basis of fibre diffraction studies. Backbone conformational transitions can nevertheless be induced leading to an alpha gamma-flip (alpha:gamma, g-/g(+)-->t/t) or to the alpha beta gamma-flip form proposed from fibre studies (alpha:beta:gamma, g-/t/g(+)-->t/g+/t). The latter transition is also found to be linked to BI-->BII transitions (epsilon:zeta, t/g(-)-->g-/t).     

Characterization of the binding properties of SIRT2 inhibitors with a N-(3-phenylpropenoyl)-glycine tryptamide backbone

SIRT2 inhibitors with a N-(3-phenylpropenoyl)-glycine tryptamide backbone were studied. This backbone has been developed in our group, and it is derived from a compound originally found by virtual screening. In addition, compounds with a smaller 3-phenylpropenoic acid tryptamide backbone were also included in the study. Binding modes for the new compounds and the previously reported compounds were analyzed with molecular modelling methods. The approach, which included a combination of molecular dynamics, molecular docking and cluster analysis, showed that certain docking poses were favourable despite the conformational variation in the target protein. The N-(3-phenylpropenoyl)-glycine tryptamide backbone is also a good backbone for SIRT2 inhibitors, and the series of compounds includes several potent SIRT2 inhibitors.     

Binding mode analysis between membrane dipeptidase and its substrates

Membrane dipeptidase (MDP) is a membrane-bound glycoprotein involved in the hydrolysis of dipeptides, showing specific activity for dipeptides. Recent study showed that membrane dipeptidase was the receptor for a lung-targeting peptide identified by in vivo phage display and the crystal structure of the cilastatin-liganded human renal dipeptidase was determined. We performed a pharmacophore-based virtual screening and molecular docking in order to characterize MDP binding interactions with its substrates. A ligand-based pharmacophore model represented only a slight enrichment because of a lacked variety and centralization of ligand features. Molecular docking study was used to incorporate ligand conformational changes in the binding sites and the performance was much better than pharmacophore model; only 10% of compound library needed to be screened in order to detect all included active compounds. In addition, we found that one of the crystallographically observed water molecules plays an important role in the binding modes between MDP and its substrate.