An energetic correlation of ab initio and NMR studies of the 3'-gauche effect in 3'-substituted thymidines

A straightforward correlation of our experimental NMR findings on 3'-substituted thymidine derivatives with that of the ab initio calculations shows that (i) the delta Go298kNRM of N reversible S equilibrium in nucleoside can be predicted from the ab initio calculated delta ES-N obtained from 6-311++G** level of theory; (ii) the substituent-dependent steric and stereoelectronic effects on the bias of the two-state N reversible S equilibrium in nucleosides can also be predicted from the ab initio calculations with sufficiently large basis functions, and (iii) the necessity of mimicking the solvation behaviour of the experimental NMR measurement condition in the ab initio calculations of biomolecules is also emphasized.     

Calculation and display of electrostatic potentials

A general methodology is developed for incorporating accurate electrostatic information from ab initio molecular orbital calculations into molecular mechanics calculations. Examples are given of the method applied to simple aromatic organic molecules. A program has been developed for displaying the results of the ab initio calculations on a Silicon Graphics workstation. The technique developed here provides an alternative method for including electrostatic interactions in molecular mechanics calculations and is compared with other methods for determining atomic charges.     

High-level ab initio calculations on the energetics of low-lying spin states of biologically relevant transition metal complexes: a first progress report

Although DFT is the unrivaled method of choice for quantum chemical studies of bioinorganic problems, little is known about its ability to predict the energetics of the low-lying electronic states of transition metal complexes. The first high-level ab initio calculations aimed at calibrating DFT vis-a-vis this issue indicate that, despite its many successes, DFT is far from infallible. In the short term, additional calibration of DFT against more elaborate ab initio methods remains an important goal for computational bioinorganic researchers. In the longer term, we are optimistic that high-level ab initio methods such as CASPT2 and CCSD(T) will be regularly used to study realistic molecules of genuine biochemical interest.     

The Hydrolysis Mechanism of Formamide Revisited: Comparison Between ab initio, Semiempirical and DFT Results

The hydrolysis of amides is a model reaction to study peptide hydrolysis. This process has been previously considered in the literature at the ab initio level. In this work, we revisit different reaction mechanisms (water-assisted, non-assisted, neutral and acid-catalyzed) with various theoretical methods : semiempirical, ab initio and Density Functional. The ab initio calculations are carried out at a computational level which is substantially higher than in previous studies. We describe the structure of the transition states and discuss the influence of the catalyst. We also compute the activation free energies for these processes at the Density Functional Theory level. Comparison of the methods allows to outline the main trends of these theoretical approaches which may be useful to design new computational strategies for investigating biological reaction mechanisms through the use of combined Quantum Mechanics/Molecular Mechanics methods.     

Ab initio vibrational calculations on ara-T molecule: application to analysis of IR and Raman spectra

The FTIR and FT-Raman spectra are reported for the arabinonucleoside ara-T (1-beta-D-arabinofuranosylthymine), which shows antiviral activity. The accurate knowledge of the vibrational modes is a prerequisite for the elucidation of drug-nucleotide and drug-enzyme interactions. The FTIR and FT-Raman spectra of ara-T were recorded from 4000 to 30 cm(-1). A tetradeuterated derivative (deuteration at N3, and hydroxyl groups O'2, O'3, and O'5) was synthesized and the observed isotopic shifts in its spectra were used for the vibrational analysis of ara-T. The theoretical frequencies and the potential energy distribution (PED) of the vibrational modes of ara-T were calculated using the ab initio Hartree-Fock/3-21G method. An assignment of the vibrational spectra of ara-T is proposed considering the scaled PED and the observed band shifts under deuteration. The scaled ab initio frequencies were in reasonable agreement with the experimental data.     

Trypsin revisited: crystallography AT (SUB) atomic resolution and quantum chemistry revealing details of catalysis

A series of crystal structures of trypsin, containing either an autoproteolytic cleaved peptide fragment or a covalently bound inhibitor, were determined at atomic and ultra-high resolution and subjected to ab initio quantum chemical calculations and multipole refinement. Quantum chemical calculations reproduced the observed active site crystal structure with severe deviations from standard stereochemistry and indicated the protonation state of the catalytic residues. Multipole refinement directly revealed the charge distribution in the active site and proved the validity of the ab initio calculations. The combined results confirmed the catalytic function of the active site residues and the two water molecules acting as the nucleophile and the proton donor. The crystal structures represent snapshots from the reaction pathway, close to a tetrahedral intermediate. The de-acylation of trypsin then occurs in true SN2 fashion.     

Progress in protein structure prediction: assessment of CASP3.

The third comparative assessment of techniques of protein structure prediction (CASP3) was held during 1998. This is a blind trial in which structures are predicted prior to having knowledge of the coordinates, which are then revealed to enable the assessment. Three sections at the meeting evaluated different methodologies - comparative modelling, fold recognition and ab initio methods. For some, but not all of the target coordinates, high quality models were submitted in each of these sections. There have been improvements in prediction techniques since CASP2 in 1996, most notably for ab initio methods.     

Study of nonplanarity of peptide bond using theoretical calculations

The conformational dependence of nonplanarity of the peptide bond of formylglycinamide has been studied using ab initio and density functional theory methods. Hartree-Fock self-consistent field theory (HF), M?ller-Plesset perturbation theory (MP2) of ab initio and B3LYP level of theory of dft method have been used employing 6-31++G** basis set. The MP2 method predicts better results than HF and B3LYP levels of theory for conformational stability dependence of nonplanarity. Systematic dependence of planarity deviation has been observed in MP2 theory. The chemical hardness values successfully predict the conformational region, but fail to obey maximum hardness principle. It is concluded that the most reliable dft method could not successfully predict the planarity of peptide bond in comparison with electron correlated method of ab initio method.     

蛋白质结构从头预测方法研究进展

蛋白质结构从头预测是不依赖模板仅从氨基酸序列信息得到天然结构。它的关键是正确定义能量函数、精确选用计算机搜索算法来寻找能量最低值。基于此,本文系统介绍了能量函数和构象搜索策略,并列举了几种比较成功的从头预测方法,通过比较得出结论:基于统计学知识的能量函数是近年来从头预测发展的主要方向,现有从头预测的构象搜索都用到Monte Carlo法。这表明随着蛋白质结构预测研究的深入,能量函数的构建、构象搜索方法的选择、大分子蛋白质结构的从头预测等关键性问题都取得了突破性进展。     

The nature of intermolecular interactions between aromatic amino acid residues

The nature of intermolecular interactions between aromatic amino acid residues has been investigated by a combination of molecular dynamics and ab initio methods. The potential energy surface of various interacting pairs, including tryptophan, phenilalanine, and tyrosine, was scanned for determining all the relevant local minima by a combined molecular dynamics and conjugate gradient methodology with the AMBER force field. For each of these minima, single-point correlated ab initio calculations of the binding energy were performed. The agreement between empirical force field and ab initio binding energies of the minimum energy structures is excellent. Aromatic-aromatic interactions can be rationalized on the basis of electrostatic and van der Waals interactions, whereas charge transfer or polarization phenomena are small for all intermolecular complexes and, particularly, for stacked structures. Proteins 2002;48:117-125.     

The distribution of wearout over evolved reliability structures

A multiple-integral equation, termed the wearout equation, describes the distribution of wearout (or aging) over evolved reliability structures, such as organisms and self-replicating machines, and thus statistically governs virtually all aging properties of the systems. The equation is applied to the computation of ab initio ("from the beginning") life tables for four natural populations of ungulates--wild boar, Dall sheep, African buffalo, and hippopotamus--which represent a broad range of survival characteristics. The good agreement of the ab initio and empirical tables, the best available for testing the theory, demonstrates the basic realism of the wearout equation. If the equation withstands further experimental testing, its analysis may provide insight into fundamental questions in the biology of aging.     

Application of statistical potentials to protein structure refinement from low resolution ab initio models

Recently ab initio protein structure prediction methods have advanced sufficiently so that they often assemble the correct low resolution structure of the protein. To enhance the speed of conformational search, many ab initio prediction programs adopt a reduced protein representation. However, for drug design purposes, better quality structures are probably needed. To achieve this refinement, it is natural to use a more detailed heavy atom representation. Here, as opposed to costly implicit or explicit solvent molecular dynamics simulations, knowledge-based heavy atom pair potentials were employed. By way of illustration, we tried to improve the quality of the predicted structures obtained from the ab initio prediction program TOUCHSTONE by three methods: local constraint refinement, reduced predicted tertiary contact refinement, and statistical pair potential guided molecular dynamics. Sixty-seven predicted structures from 30 small proteins (less than 150 residues in length) representing different structural classes (alpha, beta, alpha;/beta) were examined. In 33 cases, the root mean square deviation (RMSD) from native structures improved by more than 0.3 A; in 19 cases, the improvement was more than 0.5 A, and sometimes as large as 1 A. In only seven (four) cases did the refinement procedure increase the RMSD by more than 0.3 (0.5) A. For the remaining structures, the refinement procedures changed the structures by less than 0.3 A. While modest, the performance of the current refinement methods is better than the published refinement results obtained using standard molecular dynamics.     

Ab initio prediction of the three-dimensional structure of a de novo designed protein: a double-blind case study

Ab initio structure prediction and de novo protein design are two problems at the forefront of research in the fields of structural biology and chemistry. The goal of ab initio structure prediction of proteins is to correctly characterize the 3D structure of a protein using only the amino acid sequence as input. De novo protein design involves the production of novel protein sequences that adopt a desired fold. In this work, the results of a double-blind study are presented in which a new ab initio method was successfully used to predict the 3D structure of a protein designed through an experimental approach using binary patterned combinatorial libraries of de novo sequences. The predicted structure, which was produced before the experimental structure was known and without consideration of the design goals, and the final NMR analysis both characterize this protein as a 4-helix bundle. The similarity of these structures is evidenced by both small RMSD values between the coordinates of the two structures and a detailed analysis of the helical packing.     

Beta(2)-microglobulin and its deamidated variant, N17D form amyloid fibrils with a range of morphologies in vitro.

Amyloid fibrils formed by incubation of recombinant wild-type human beta(2)-microglobulin (beta(2)M) ab initio in vitro at low pH and high ionic strength are short and highly curved. By contrast, fibrils extracted from patients suffering from haemodialysis-related amyloidosis and those formed by seeding growth of the wild-type protein in vitro with fibrils ex vivo are longer and straighter than those previously produced ab initio in vitro. Here we explore the effect of growth conditions on morphology of beta(2)M fibrils formed ab initio in vitro from the wild-type protein, as well as a variant form of beta(2)M in which Asn17 is deamidated to Asp (N17D). We show that deamidation results in significant destabilisation of beta(2)M at neutral pH. Despite this, acidification is still necessary to form amyloid from the mutant protein in vitro. Interestingly, at low pH and low ionic strength long, straight fibrils of recombinant beta(2)M are formed in vitro. The fibrils comprise three distinct morphological types when examined using electron microscopy (EM) and atomic force microscopy (AFM) that vary in periodicity and the number of constituent protofibrils. Using kinetic experiments we suggest that the immature fibrils observed previously do not represent intermediates in the assembly of fully mature amyloid, at least under the conditions studied here.     

Constructing side chains on near-native main chains for ab initio protein structure prediction

Is there value in constructing side chains while searching protein conformational space during an ab initio simulation? If so, what is the most computationally efficient method for constructing these side chains? To answer these questions, four published approaches were used to construct side chain conformations on a range of near-native main chains generated by ab initio protein structure prediction methods. The accuracy of these approaches was compared with a naive approach that selects the most frequently observed rotamer for a given amino acid to construct side chains. An all-atom conditional probability discriminatory function is useful at selecting conformations with overall low all-atom root mean square deviation (r.m.s.d.) and the discrimination improves on sets that are closer to the native conformation. In addition, the naive approach performs as well as more sophisticated methods in terms of the percentage of chi(1) angles built accurately and the all-atom r. m.s.d., between the native and near-native conformations. The results suggest that the naive method would be extremely useful for fast and efficient side chain construction on vast numbers of conformations for ab initio prediction of protein structure.     

Quantum chemistry and enzymes: a next step

In this paper the relevance of ab initio quantum mechanical calculations is briefly reviewed. A method to extend such calculations to the domain of solvent effects and the theoretical investigation of chemical reactions in the active site of a protein is discussed.     

Ab initio conformational maps for disaccharides in gas phase and aqueous solution

Ab initio conformational maps for beta-lactose in both the gas phase and in aqueous solution have been constructed at the HF/6-31G(d,p) level of calculation. The results of the gas-phase ab initio calculations allow us to conclude that a rigid conformational map is able to predict the regions of the minima in the potential energy surface of beta-lactose, in full agreement with those found in the relaxed conformational map. The solvation effects do not give rise to any new local minimum in the potential energy surface of beta-lactose, but just change the relative Boltzmann populations of the conformers found in the gas-phase calculations. The values obtained for heteronuclear spin coupling constant (3J(H,C)), using the seven most stable conformers in solution are in good agreement with the available experimental values. This is a good indication that ab initio rigid conformational maps can be reliably used to sort the most stable conformers of beta-lactose.     

HCPM--program for hierarchical clustering of protein models

HCPM is a tool for clustering protein structures from comparative modeling, ab initio structure prediction, etc. A hierarchical clustering algorithm is designed and tested, and a heuristic is provided for an optimal cluster selection. The method has been successfully tested during the CASP6 experiment.