The structure of the gas-phase tyrosine-glycine-glycine tripeptide

An Atomic Counterpoise Method for Estimating Inter- and Intramolecular Basis Set Superposition Errors

An atomic counterpoise method is proposed to calculate estimates of inter- and intramolecular basis set superposition errors. The method estimates the basis set superposition error as a sum of atomic contributions and can be applied for both independent particle and electron correlation models. It is shown that the atomic counterpoise method provides results very similar to the molecular counterpoise method for intermolecular basis set superposition errors at both the HF and MP2 levels of theory with a sequence of increasingly larger basis sets. The advantage of the atomic counterpoise method is that it can be applied with equal ease to estimate intramolecular basis set superposition errors, for which few other methods exist. The atomic counterpoise method is computationally quite efficient, requiring typically double the amount of computer time as required for calculating the uncorrected energy.

Towards an understanding of the nucleation of alpha-para amino benzoic acid from ethanolic solutions: a multi-scale approach

The molecular assembly and subsequent nucleation of para-amino benzoic acid (PABA) from ethanolic solutions is probed using a multi-scale and multi-technique approach. This is applied by examining and interrelating information regarding the molecular, solution-state, cluster, solid-state and surface structures to understand why the alpha form of PABA is crystallised in preference to its low temperature beta form. Calculations suggest that conformational changes within the solute molecule play little or no role in directing the nucleation of either the alpha or beta crystal forms. Combined ab initio and molecular dynamics calculations of the stability of small clusters in solution suggests that the hydrogen-bonded carboxylic acid dimers, present in the alpha structure, are the most stable in solution and play a major role in the self-assembly and polymorphic expression of the alpha form in ethanol in preference to the beta form. These calculations are in good agreement with X-ray small-angle scattering analysis which reveals the presence of PABA clusters in ethanol which are consistent with the size and shape of a carboxylic acid dimer. SAXS studies also reveal the presence of larger cluster structures in a size range 10-40 nm which appear to grow, perhaps reflecting a change in the balance between monomers and dimers within the solution during the nucleation process. The results of crystallisation-kinetics experiments indicate an instantaneous nucleation mechanism where the number of instantaneously nucleated crystallites is calculated to be 1360-660 nuclei per ml and the subsequent growth is found to be only rate limited by diffusion of the growth unit to the crystallite surface. A linear dependence of growth rate with respect to supersaturation is observed for the (0 1 -1) capping face, which is associated with strong π-π stacking interactions. This is consistent with a solid-on-solid mechanism associated with surface roughened growth and concomitant poor lattice-perfection. Conversely, the side (1 0 -1) surface has a growth mechanism consistent with a 2D nucleation birth and spread mechanism. Hence, these mechanisms result in very fast growth along the b-axis and the needle-like morphology that is observed for alpha-PABA.

Comprehensive analysis of energy minima of the 20 natural amino acids

Energy minima of the 20 natural amino acids (capped by a peptide bond at both the N and the C termini, CH3-C(═O)-N(H)-(H)Cα(R)-C(═O)-N(H)-CH3), were obtained by ab initio geometry optimization. Starting with a large number of minima, quickly generated by MarvinView, geometry optimization at the HF/6-31G(d,p) level of theory reduced the number of minima, followed by further optimization at the B3LYP/apc-1 and MP2/cc-pVDZ levels, which caused some minima to disappear and some stable minima to migrate on the Ramachandran map. There is a relation between the number of minima and the size and the flexibility of the side chain. The energy minima of the 20 amino acids are mainly located in the regions of βL, γL, δL, and αL of the Ramachandran map. Multipole moments of atoms occurring in the fragment [-NH-Cα-C(═O)-] common to all 20 amino acids were calculated at the three levels of theory mentioned above. The near parallelism in behavior of these moments between levels of theory is beneficial toward estimating moments with the more expensive B3LYP and MP2 methods from data calculated with the cheaper HF method. Finally, we explored the transferability of properties between different amino acids: the bond length and angles of the common fragment [-NH-Cα(HαCβ)-C'(═O)-] in all amino acids except Gly and Pro. All bond lengths are highly transferable between different amino acids, and the standard deviations are small.

Comprehensive conformational studies of five tripeptides and a deduced method for efficient determinations of peptide structures

Thorough searches on the potential energy surfaces of five tripeptides, GGG, GYG, GWG, TGG, and MGG, were performed by considering all possible combinations of the bond rotational degrees of freedom with a semiempirical and ab initio combined computational approach. Structural characteristics of the obtained stable tripeptide conformers were carefully analyzed. Conformers of the five tripeptides were found to be closely connected with conformers of their constituting dipeptides and amino acids. A method for finding all important tripeptide conformers by optimizing a small number of trial structures generated by suitable superposition of the parent amino acid and dipeptide conformers is thus proposed. Applying the method to another five tripeptides, YGG, FGG, WGG, GFA, and GGF, studied before shows that the new approach is both efficient and reliable by providing the most complete ensembles of tripeptide conformers. The method is further generalized for application to larger peptides by introducing the breeding and mutation concepts in a genetic algorithm way. The generalized method is verified to be capable of finding tetrapeptide conformers with secondary structures of strands, helices, and turns, which are highly populated in larger peptides. This show some promise for the proposed method to be applied for the structural determination of larger peptides.

The role of weakly polar and H-bonding interactions in the stabilization of the conformers of FGG, WGG, and YGG: an aqueous phase computational study

The energetics of intramolecular interactions on the conformational potential energy surface of the terminally protected N-Ac-Phe-Gly-Gly-NHMe (FGG), N-Ac-Trp-Gly-Gly-NHMe (WGG), and N-Ac-Tyr-Gly-Gly-NHMe (YGG) tripeptides was investigated. To identify the representative conformations, simulated annealing molecular dynamics (MD) and density functional theory (DFT) methods were used. The interaction energies were calculated at the BHandHLYP/aug-cc-pVTZ level of theory. In the global minima, 10%, 31%, and 10% of the stabilization energy come from weakly polar interactions, respectively, in FGG, WGG, and YGG. In the prominent cases 46%, 62%, and 46% of the stabilization energy is from the weakly polar interactions, respectively, in FGG, WGG, and YGG. On average, weakly polar interactions account for 15%, 34%, and 9% of the stabilization energies of the FGG, WGG, and YGG conformers, respectively. Thus, weakly polar interactions can make an important energetic contribution to protein structure and function.