alpha- and 3(10)-helix interconversion: a quantum-chemical study on polyalanine systems in the gas phase and in aqueous solvent.

Helices are among the predominant secondary structures in globular proteins. About 90% of the residues in them are found to be in the alpha-helical conformation, and another 10% in the 3(10) conformation. There is a standing controversy between experimental and some theoretical results, and controversy among theoretical results concerning the predominance of each conformation, in particular, helices. We address this controversy by ab initio Hartree-Fock and density functional theory studies of helices with different lengths in a vacuum and in the aqueous phase. Our results show that (1) in a vacuum, all oligo(Ala) helices of 4-10 residues adopt the 3(10) - conformation; (2) in aqueous solution, the 6-10 residue peptides adopt the alpha-helical conformation; (3) there might be two intermediates between these helical conformers allowing for their interconversion. The relevance of these results to the structure and folding of proteins is discussed.     

Why Are B ions stable species in peptide spectra?

Protonated amino acids and derivatives RCH(NH₂)C(+O)X · H⁺ (X = OH, NH₂, OCH₃) do not form stable acylium ions on loss of HX, but rather the acylium ion eliminates CO to form the immonium ion RCH = NH ₂ ⁺ . By contrast, protonated dipeptide derivatives H₂NCH(R)C(+O)NHCH(R′)C(+O)X · H⁺ [X = OH, OCH₃, NH₂, NHCH(R″)COOH] form stable B₂ ions by elimination of HX. These B₂ ions fragment on the metastable ion time scale by elimination of CO with substantial kinetic energy release (T _1/2 = 0.3–0.5 eV). Similarly, protonated N-acetyl amino acid derivatives CH₃C(+O)NHCH(R′)C(+O)X · H⁺ [X = OH, OCH₃, NH₂, NHCH(R″)COOH] form stable B ions by loss of HX. These B ions also fragment unimolecularly by loss of CO with T _1/2 values of ～ 0.5 eV. These large kinetic energy releases indicate that a stable configuration of the B ions fragments by way of activation to a reacting configuration that is higher in energy than the products, and some of the fragmentation exothermicity of the final step is partitioned into kinetic energy of the separating fragments. We conclude that the stable configuration is a protonated oxazolone, which is formed by interaction of the developing charge (as HX is lost) with the N-terminus carbonyl group and that the reacting configuration is the acyclic acylium ion. This conclusion is supported by the similar fragmentation behavior of protonated 2-phenyl-5-oxazolone and the B ion derived by loss of H-Gly-OH from protonated C₆H₅C(+O)-Gly-Gly-OH. In addition, ab initio calculations on the simplest B ion, nominally HC(+O)NHCH₂CO⁺, show that the lowest energy structure is the protonated oxazolone. The acyclic acylium isomer is 1.49 eV higher in energy than the protonated oxazolone and 0.88 eV higher in energy than the fragmentation products, HC(+O)N⁺H = CH₂ + CO, which is consistent with the kinetic energy releases measured.     

Conformational study of protonated,neutral, and deprotonated formamide

Neutral, protonated, and deprotonated formamide isomers were studied at the 3-21G SCF level with complete geometry optimization. Ten stable structures, ten first-order saddle points, and three second-order saddle points (conformational maxima) are reported. [Total energies are reported in hartrees (1 hartree = 627.51 kcal/mol = 2625.5 kJ/mol) and energy differences are reported in kJ/mol (1 kJ/mol = 0.239 kcal/mol).] Rotational barriers and proton affinities are discussed and compared to isoelectronic amidine species.     

Finiteness of the quadratic primal simplex method when s-monotone index selection rules are applied

Central European Journal of Operations Research - This paper considers the primal quadratic simplex method for linearly constrained convex quadratic programming problems. Finiteness of the...     

Prospects in computational molecular medicine: a millennial mega-project on peptide folding

During the second half of the 20th century, Molecular Computations have reached to a level that can revolutionize chemistry. The next target will be structural biology, which will be followed soon by Molecular Medicine. The present paper outlines where we are at, in this field, at the end of the 20th century, and in what direction the development may take in the new millennium. In view of the gigantic nature of the problem, it is suggested that a suitably designed cooperative Millennial Mega-project might accelerate our schedule.     

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.     

Thermodynamic functions of conformational changes: conformational network of glycine diamide folding, entropy lowering, and informational accumulation

A refined grid of a conformational potential energy surface (PES) and a conformational entropy surface for glycine diamide was generated by ab initio molecular computations. The possible network of reaction paths was recognized in terms of the linear combinations of internal coordinates corresponding to conrotatory and disrotatory modes of motions. Such a Woodward-Hoffmann-like path selection principle was detected for the folding of this peptide from extended to some virtually cyclic structure. It seemed reasonable to assume that this principle (or its generalized form) might be applicable to protein folding. A reaction path network was projected on the potential energy, and a continuous entropy surface was constructed under the condition of reduced dimensionality. The low entropy of the folded conformation indicated an information accumulation between 326% and 1414% with respect to the fully extended or unfolded structure. It is found that the location of existing and 'latent' critical points on the surface is revealed by the extrema and inflection points of the entropy curve.     

Thermodynamic functions of conformational changes. 2. Conformational entropy as a measure of information accumulation

Various folded molecular structures contain different amount of information. The relative amount of information may be related to relative entropy or entropy change. The conformational entropy change for n-butane has been computed as the function of rotation around the central C-C bond. It appears that the g+ or g- conformers contain about 16% more information than the anti-structure. Furthermore, the syn conformation with the two groups eclipsed contained about 42% more information than the fully staggered anti orientation. The conformational entropy function calculated from 3N - 7 internal degrees of freedom was found to be a continuous function.     

Making "real" molecules in virtual space

Predicting "realistic" compounds of given chemical reactions with virtual synthesis tools usually requires the manual intervention of experienced chemists in the enumeration phase for the selection of appropriate reactants, assignment of the corresponding reaction sites, and removal of the unlikely products. To automate the virtual synthesis process, we have moved the expertise intensive parts from the compound library design phase to the reaction library design phase. ChemAxon is building an in silico reaction library containing important preparative transformations, where each reaction definition contains a generic transformation scheme and additional rules to handle the various starting compounds according to the corresponding chemo-, regio-, and stereoselectivity issues. Having well designed reaction definitions in hand, our software tool is able to generate synthetically feasible compound libraries with minimal effort in the enumeration phase.