Gábor Pohl

Secondary structure of short β-peptides as the chiral expression of monomeric building units: A rational and predictive model

Chirality of the monomeric residues controls and determines the prevalent folding of small oligopeptides (from di- to tetramers) composed of 2-aminocyclobutane-1-carboxylic acid (ACBA) derivatives with the same or different absolute and relative configuration. The cis-form of the monomeric ACBA gives rise to two conformers, namely, Z6 and Z8, while the trans-form manifests uniquely as an H8 structure. By combining these subunits in oligo- and polypeptides, their local structural preference remains, thus allowing the rational design of new short foldamers. A lego-type molecular architecture evolves; the overall look depends only on the conformational properties of the structural building units. A versatile and efficient method to predict the backbone folds of designed cyclobutane β-peptides is based on QM calculations. Predictions are corroborated by high-resolution NMR studies on selected stereoisomers, most of them being new foldamers that have been synthesized and characterized for the first time. Thus, the chiral expression of monomeric building units results in the defined secondary structures of small oligomers. As a result of this study, a new set of chirality controlled foldamers is provided to probe as biocompatible biopolymers.

Extended Apolar β-Peptide Foldamers: The Role of Axis Chirality on β-Peptide Sheet Stability

This study is on structure and stability of sheetlike conformers of beta-peptides; never seen new foldamers are reported here for the first time. Single- and double-stranded structures are analyzed, and the seeds of large beta-layers and biocompatible nanomaterials are described here. Both the monomeric, HCO-[NH-CH(2)-CH(2)CO](n)-NH(2), and dimeric forms, [HCO-(beta-Ala)(n)-NH(2)](2) n = 3 and 4, of oligo-beta-alanine supramolecular complexes are evaluated by using an adequate level of theory M052X/6-31G(d) for peptides of this size. Polymers composed of backbone foldamers with the central mu torsion angle set to an anti orientation were all probed. Sheet structures built up of strands with carbonyl groups monotonically facing the same spatial direction, polar strands, were previously assigned and synthesized ( Seebach , D. Chem. Biodiversity 2004 , 1 , 1111 - 1239 ). Now we are presenting a novel beta-peptide sheet structure of alternating carbonyl group orientations, called as apolar strands. These novel secondary structural elements of beta-peptides are structural analogs of beta-pleated sheets of proteins. Interestingly enough, the latter type of apolar strands are foreseen as very stable supramolecular complexes and are more firm by approximately 10 kcal.mol(-1) than the aforementioned polar strands. Furthermore, apolar strands lack the inherent twisting of beta-layers, present in polar strands resulting in the tubular shape. Once the effect of substitution of Hbeta1 and/or Hbeta2 atoms are revealed on foldamer stability, short peptide sequence could be designed and synthesized. These new, conformationally optimized beta-sheetlike nanostructures of increased stability with little or no twisting could be used as enzymatically resistant ( Frackenpohl , J. , Arvidsson , P. I. , Schreiber , J. V. , and Seebach , D. ChemBioChem 2001 , 2 , 445 - 455 ) biomaterials. These newly designed models systems could enlarge the arsenal of durable polyesters of similar chemical constitution (e.g., -[O-CH(CH(3))-CH(2)CO](n)- and -[O-CH(COOH)-CH(2)CO](n)-) already used as artificial heart valves, for example.

The interactions of phenylalanines in β-sheet-like structures from molecular orbital calculations using density functional theory (DFT), MP2, and CCSD(T) methods

We present density functional theory calculations designed to evaluate the importance of π-stacking interactions to the stability of in-register Phe residues within parallel β-sheets, such as amyloids. We have used a model of a parallel H-bonded tetramer of acetylPheNH2 as a model and both functionals that were specifically designed to incorporate dispersion effects (DFs), as well as, several traditional functionals which have not been so designed. None of the functionals finds a global minimum for the π-stacked conformation, although two of the DFs find this to be a local minimum. The stacked phenyls taken from the optimized geometries calculated for each functional have been evaluated using MP2 and CCSD(T) calculations for comparison. The results suggest that π-stacking does not make an important contribution to the stability of this system and (by implication) to amyloid formation.

Capping Amyloid β‑Sheets of the Tau-Amyloid Structure VQIVYK with Hexapeptides Designed To Arrest Growth. An ONIOM and Density Functional Theory Study

We present ONIOM calculations using density functional theory (DFT) as the high and AM1 as the medium level that explore the abilities of different hexapeptide sequences to terminate the growth of a model for the tau-amyloid implicated in Alzheimer’s disease. We delineate and explore several design principles (H-bonding in the side chains, using antiparallel interactions on the growing edge of a parallel sheet, using all-d residues to form rippled interactions at the edge of the sheet, and replacing the H-bond donor N–H’s that inhibit further growth) that can be used individually and in combination to design such peptides that will have a greater affinity for binding to the parallel β-sheet of acetyl-VQIVYK-NHCH3 than the natural sequence and will prevent another strand from binding to the sheet, thus providing a cap to the growing sheet that arrests further growth. We found peptides in which the Q is replaced by an acetyllysine (aK) residue to be particularly promising candidates, particularly if the reverse sequence (KYVIaKV) is used to form an antiparallel interaction with the sheet.

Ab Initio Molecular Dynamics Study of Solvated Electrons in Methanol Clusters

We performed a series of ab initio molecular dynamics simulations to investigate the physical properties of small methanol cluster anions, [(CH3OH)n]-, (n = 8-32). An excess electron was attached to neutral clusters that were prepared to accommodate the electron in interior cavity states or surface bound states. The computed initial binding energies of the electrons to these clusters indicate appealing similarity to the experimentally observed vertical detachment energies. The tendency of the interior state clusters parallels that of the clusters with strong electron binding in the experiments, while the simulated unrelaxed surface state anions are similar to the observed weakly bound species. This assignment is consistent with a previous identification based on hybrid quantum-classical simulations. The time evolution of the cluster anions suggests that interior state electrons slowly move to and relax on the surface, in excess electronic states that appear significantly more stable than the experimentally assigned putative surface states. Based on this result we predict the existence of relaxed surface state isomers of small methanol cluster anions. Due to the kinetic metastability of the experimentally found weakly bound species, we anticipate a serious technical challenge to prepare and identify small methanol cluster anions with relaxed surface states. These more strongly binding surface states are stabilized by dangling hydroxyl hydrogen atoms pointing to the excess electrons charge distribution. In addition, methyl hydrogens also appear to contribute to the stability of these states. During its transition to the surface, the interior excess electron maintains its initial solvent cavity. No signs of non-cavity interior states are observed in the present first principles ab initio molecular dynamics simulations.

Excess electrons in methanol clusters: Beyond the one-electron picture

We performed a series of comparative quantum chemical calculations on various size negatively charged methanol clusters, CH3OHn-. The clusters are examined in their optimized geometries (n = 2-4), and in geometries taken from mixed quantum-classical molecular dynamics simulations at finite temperature (n = 2-128). These latter structures model potential electron binding sites in methanol clusters and in bulk methanol. In particular, we compute the vertical detachment energy (VDE) of an excess electron from increasing size methanol cluster anions using quantum chemical computations at various levels of theory including a one-electron pseudopotential model, several density functional theory (DFT) based methods, MP2 and coupled-cluster CCSD(T) calculations. The results suggest that at least four methanol molecules are needed to bind an excess electron on a hydrogen bonded methanol chain in a dipole bound state. Larger methanol clusters are able to form stronger interactions with an excess electron. The two simulated excess electron binding motifs in methanol clusters, interior and surface states, correlate well with distinct, experimentally found VDE tendencies with size. Interior states in a solvent cavity are stabilized significantly stronger than electron states on cluster surfaces. Although we find that all the examined quantum chemistry methods more or less overestimate the strength of the experimental excess electron stabilization, MP2, LC-BLYP, and BHandHLYP methods with diffuse basis sets provide a significantly better estimate of the VDE than traditional DFT methods (BLYP, B3LYP, X3LYP, PBE0). A comparison to the better performing many electron methods indicates that the examined one-electron pseudopotential can be reasonably used in simulations for systems of larger size.