Stability issues of covalently and noncovalently bonded peptide subunits

The Quixote project: Collaborative and Open Quantum Chemistry data management in the Internet age

Computational Quantum Chemistry has developed into a powerful, efficient, reliable and increasingly routine tool for exploring the structure and properties of small to medium sized molecules. Many thousands of calculations are performed every day, some offering results which approach experimental accuracy. However, in contrast to other disciplines, such as crystallography, or bioinformatics, where standard formats and well-known, unified databases exist, this QC data is generally destined to remain locally held in files which are not designed to be machine-readable. Only a very small subset of these results will become accessible to the wider community through publication.

In this paper we describe how the Quixote Project is developing the infrastructure required to convert output from a number of different molecular quantum chemistry packages to a common semantically rich, machine-readable format and to build respositories of QC results. Such an infrastructure offers benefits at many levels. The standardised representation of the results will facilitate software interoperability, for example making it easier for analysis tools to take data from different QC packages, and will also help with archival and deposition of results. The repository infrastructure, which is lightweight and built using Open software components, can be implemented at individual researcher, project, organisation or community level, offering the exciting possibility that in future many of these QC results can be made publically available, to be searched and interpreted just as crystallography and bioinformatics results are today.

Although we believe that quantum chemists will appreciate the contribution the Quixote infrastructure can make to the organisation and and exchange of their results, we anticipate that greater rewards will come from enabling their results to be consumed by a wider community. As the respositories grow they will become a valuable source of chemical data for use by other disciplines in both research and education.

The Quixote project is unconventional in that the infrastructure is being implemented in advance of a full definition of the data model which will eventually underpin it. We believe that a working system which offers real value to researchers based on tools and shared, searchable repositories will encourage early participation from a broader community, including both producers and consumers of data. In the early stages, searching and indexing can be performed on the chemical subject of the calculations, and well defined calculation meta-data. The process of defining more specific quantum chemical definitions, adding them to dictionaries and extracting them consistently from the results of the various software packages can then proceed in an incremental manner, adding additional value at each stage.

Not only will these results help to change the data management model in the field of Quantum Chemistry, but the methodology can be applied to other pressing problems related to data in computational and experimental science.

How stable is a collagen triple helix? An ab initio study on various collagen and beta-sheet forming sequences

Collagen forms the well characterized triple helical secondary structure, stabilized by interchain H-bonds. Here we have investigated the stability of fully optimized collagen triple helices and beta-pleated sheets by using first principles (ab initio and DFT) calculations so as to determine the secondary structure preference depending on the amino acid composition. Models composed of a total of 18 amino acid residues were studied at six different amino acid compositions: (i) L-alanine only, (ii) glycine only, (iii) L-alanines and glycine, (iv) L-alanines and D-alanine, (v) L-prolines with glycine, (vi) L-proline, L-hydroxyproline, and glycine. The last two, v and vi, were designed to mimic the core part of collagen. Furthermore, ii, iii, and iv model the binding and/or recognition sites of collagen. Finally, i models the G-->A replacement, rare in collagen. All calculated structures show great resemblance to those determined by X-ray crystallography. Calculated triple helix formation affinities correlate well with experimentally determined stabilities derived from melting point (T(m)) data of different collagen models. The stabilization energy of a collagen triple helical structure over that of a beta-pleated sheet is 2.1 kcal mol(-1) per triplet for the [(-Pro-Hyp-Gly-)(2)](3) collagen peptide. This changes to 4.8 kcal mol(-1) per triplet of destabilization energy for the [(-Ala-Ala-Gly-)(2)](3) sequence, known to be disfavored in collagen. The present study proves that by using first principles methods for calculating stabilities of supramolecular complexes, such as collagen and beta-pleated sheets, one can obtain stability data in full agreement with experimental observations, which envisage the applicability of QM in molecular design.

Theoretical study on tertiary structural elements of beta-peptides: nanotubes formed from parallel-sheet-derived assemblies of beta-peptides

Parallel or polar strands of beta-peptides spontaneously form nanotubes of different sizes in a vacuum as determined by ab initio calculations. Stability and conformational features of [CH3CO-(beta-Ala)k-NHCH3]l (1 < or = k < or = 4, 2 < or = l < or = 4) models were computed at different levels of theory (e.g., B3LYP/6-311++G(d,p)// B3LYP/6-31G(d), with consideration of BSSE). For the first time, calculations demonstrate that sheets of beta-peptides display nanotubular characteristics rather than two-dimensional extended beta-layers, as is the case of alpha-peptides. Of the configurations studied, k = l = 4 gave the most stable nanotubular structure, but larger assemblies are expected to produce even more stable nanotubes. As with other nanosystems such as cyclodextrane, these nanotubes can also incorporate small molecules, creating a diverse range of applications for these flexible, biocompatible, and highly stable molecules. The various side chains of beta-peptides can make these nanosystems rather versatile. Energetic and structural features of these tubular model systems are detailed in this paper. It is hoped that the results presented in this paper will stimulate experimental research in the field of nanostructure technology involving beta-peptides.

Determining suitable lego-structures to estimate stability of larger peptide nanostructures using computational methods

Nanofibers, nanofilms and nanotubes constructed of one to four strands of oligo-alpha- and oligo-beta-peptides were obtained by using carefully selected building units. Lego-type approaches based on thermoneutral isodesmic reactions can be used to reconstruct the total energies of both linear and tubular periodic nanostructures with acceptable accuracy. Total energies of several different nanostructures were accurately determined with errors typically falling in the subchemical range. Thus, attention will be focused on the description of suitable isodesmic reactions that have enabled the determination of the total energy of polypeptides and therefore offer a very fast, efficient and accurate method to obtain energetic information on large and even very large nanosystems.

Toward a rational design of beta-peptide structures

Intrinsic conformational characteristics of beta-peptides built up from simple achiral and chiral beta-amino acid residues (i.e., HCO-beta-Ala-NH2, HCO-beta-Abu-NH2) were studied using quantum chemical calculations and 1H-NMR spectroscopy. A conformer-based systematic and uniform nomenclature was introduced to differentiate conformers. Geometry optimizations were performed on all homoconformers of both HCO-(beta-Ala)(k)-NH2 and HCO-(beta-Abu)(k)-NH2 (1 < or = k < or = 6) model systems at the RHF/3-21G and RHF/6-311++G(d, p) levels of theory. To test for accuracy and precision, additional computations were carried out at several levels of theory [e.g., RHF/6-31G(d), and B3LYP/6-311++G(d, p)]. To display the folding preference, the relative stability of selected conformers as function of the length of the polypeptide chain was determined. Ab initio population distribution of hexapeptides and the conformational ensemble of synthetic models composed of beta-Ala and beta-Abu studied using 1H-NMR in different solvents were compared at a range of temperatures. Helical preference induced by various steric effects of nonpolar side chains was tested using higher level ab initio methods for well-known model systems such as: HCO-(beta-HVal-beta-HAla-beta-HLeu)2-NH2, HCO-(ACHC)6-NH2, HCO-(trans-ACPC)6-NH2, and HCO-(cis-ACPC)6-NH2. The relative stabilities determined by theoretical methods agreed well with most experimental data, supporting the theory that the local conformational preference influenced by steric effects is a key determining factor of the global fold both in solution and in the gas phase.

A self-stabilized model of the chymotrypsin catalytic pocket. The energy profile of the overall catalytic cycle

A model of the catalytic triad of chymotrypsin is built assuring the arrangement and properties as they are within the complete enzyme. The model contains 18 amino acid residues of chymotrypsin and its substrate. A total of 135 atoms (including 70 heavy atoms) were subjected to full ab initio geometry optimizations through 127 individual steps along the reaction coordinate of the complete catalytic mechanism. It was shown that the described model of the catalytic apparatus forms a self-stabilized molecule ensemble without the rest of the enzyme and substrate. According to the calculations, the formations of the first and second tetrahedral intermediates in the model have 20.3 and 15.7 kcal/mol activation energy barriers, respectively. Removing elements of the catalytic apparatus such as the (1) catalytic aspartate or (2) the anion hole, as well as (3) inserting a water molecule "early" in the catalytic process, or (4) introducing conformational rigidity of the substrate, results in an increase of the above energy barrier of the first catalytic step in the model by 6.4, 13.7, 3.7, and 4.1 kcal/mol, respectively. Based on the calculated process one can conclude that the catalytic reaction in this model is much more similar to the reaction in the enzyme than to the reference reaction. To our knowledge, this is the first model system that mimics the complete catalytic mechanism.

All-cis Cyclic Peptides

Amide bonds -NH-CO- preferentially exist in trans conformations, the cis conformation being thermodynamically unfavored with respect to the trans by about 2 kcal/mol. Yet, the main reason most proteins or peptides cannot be made from cis-peptide plaques only lies in that connecting them into open chains appears to be sterically impracticable. It is possible, however, to build all-cis cyclic peptides in which all cis-plaques are efficiently locked. The present work examines, through quantum calculations, the structural and energetic issues associated with these peculiar arrangements. Systematic exploration at DFT-B3LYP level of the potential-energy surfaces for all-cis cyclopolyglycines cG(n)(c) (n = 2-10,15), and to a lesser extent, for all-cis cyclopolyalanines and all-cis cyclopolyphenylalanines confirms that all these structures are true minima. Optimal ring size occurs around eight peptide units, resulting in planar cG7(c), cG8(c), and cG9(c). In smaller systems, the ring strain is relieved through nonplanar cup-like distortions, particularly in cG6(c). From 10 peptide units and beyond, the ring framework distorts into a saddle-edge shape. These molecules disclose some molecular flexibility, as combinatorial tilting of the plaques may give sets of minima close in energy. Indexes based on isodesmic reactions are used to estimate the energy for joining all-cis or all-trans plaques into cyclic peptides. One of them, the mean plaque-junction energy (MPJE) suggests that within sensible sizes from six peptide units and beyond, all-cis plaque association is almost equally favorable as all-trans one. The frame of radiating cis-amide bonds can be considered as defining a new kind of peptidic material, endowed with specific self-assembling properties.

Structure and stability of β-pleated sheets

Beside alpha-helices, beta-sheets are the most common secondary structure elements of proteins. In this article, the question of structure and stability of parallel and antiparallel sheets of various lengths is addressed. All data obtained are compared to a selected set of protein structures. In antiparallel beta-sheets, one of the two possible H-bonded structures (containing 14 atoms in the H-bonded pseudoring) is energetically more favored and also more abundant in proteins than the other one (with 10 atoms involved in the pseudoring). Parallel beta-sheets and their subunits are energetically less stable and indeed found to occur more rarely in proteins. Antiparallel hairpins are disfavored compared to beta-sheets formed by sequentially separated strands. Agreement between theory and experimental data indicates that characterization of structural building blocks at an appropriately accurate level of theory is a useful tool to get insight into fundamentals of protein structure.