Search PubMed⌕ Search

Biomedical subjects

T F Havel

Publications and source records attributed to T F Havel.

27 records · Page 2Linked to original sources

Computational experience with an algorithm for tetrangle inequality bound smoothing.

An important component of computer programs for determining the solution conformation of proteins and other flexible molecules from nuclear magnetic resonance data are the so-called "bound smoothing algorithms", which compute lower and upper limits on the values of all the interatomic distances from the relatively sparse set which can usually be measured experimentally. To date, the only methods efficient enough for use in large problems take account of only the triangle inequality, but an appreciable improvement in the precision of the limits is possible if the algebraic relations between the distances among each quadruple of atoms are also considered. The goal of this paper is to use a recently improved algorithm for computing these "tetrangle inequality limits" to determine just how much improvement really is possible, given the types of experimental data that are usually available.

Algorithms↗

Three-dimensional solution structure of plastocyanin from the green alga Scenedesmus obliquus.

The solution conformation of plastocyanin from the green alga Scenedesmus obliquus has been determined from distance and dihedral angle constraints derived by nuclear magnetic resonance (NMR) spectroscopy. Structures were generated with distance geometry and restrained molecular dynamics calculations. A novel molecular replacement method was also used with the same NMR constraints to generate solution structures of S. obliquus plastocyanin from the x-ray structure of the homologous poplar protein. Scenedesmus obliquus plastocyanin in solution adopts a beta-barrel structure. The backbone conformation is well defined and is similar overall to that of poplar plastocyanin in the crystalline state. The distinctive acidic region of the higher plant plastocyanins, which functions as a binding site for electron transfer proteins and inorganic complexes, differs in both shape and charge in S. obliquus plastocyanin.

Calorimetry↗

Protein structures in solution by nuclear magnetic resonance and distance geometry. The polypeptide fold of the basic pancreatic trypsin inhibitor determined using two different algorithms, DISGEO and DISMAN.

A set of conformational restraints derived from nuclear magnetic resonance (n.m.r.) measurements on solutions of the basic pancreatic trypsin inhibitor (BPTI) was used as input for distance geometry calculations with the programs DISGEO and DISMAN. Five structures obtained with each of these algorithms were systematically compared among themselves and with the crystal structure of BPTI. It is clear that the protein architecture observed in single crystals of BPTI is largely preserved in aqueous solution, with local structural differences mainly confined to the protein surface. The results confirm that protein conformations determined in solution by combined use of n.m.r. and distance geometry are a consequence of the experimental data and do not depend significantly on the algorithm used for the structure determination. The data obtained further provide an illustration that long intramolecular distances in proteins, which are comparable with the radius of gyration, are defined with high precision by relatively imprecise nuclear Overhauser enhancement measurements of a large number of much shorter distances.

Algorithms↗

An evaluation of the combined use of nuclear magnetic resonance and distance geometry for the determination of protein conformations in solution.

An evaluation of the potential of nuclear magnetic resonance (n.m.r.) as a means of determining polypeptide conformation in solution is performed with the aid of a new distance geometry program which is capable of computing complete spatial structures for small proteins from n.m.r. data. Ten sets of geometric constraints which simulate the results available from n.m.r. experiments of varying precision and completeness were extracted from the crystal structure of the basic pancreatic trypsin inhibitor, and conformers consistent with these constraints were computed. Comparison of these computed structures with each other and with the original crystal structure shows that it is possible to determine the global conformation of a polypeptide chain from the distance constraints which are available from n.m.r. experiments. The results obtained with the different data sets also provide a standard by which the quality of protein structures computed from n.m.r. data can be evaluated when no crystal structure is available, and indicate directions in which n.m.r. experiments for protein structure determination could be further improved.

Animals↗

Solution conformation of proteinase inhibitor IIA from bull seminal plasma by 1H nuclear magnetic resonance and distance geometry.

A determination of the solution conformation of the proteinase inhibitor IIA from bull seminal plasma (BUSI IIA) is described. Two-dimensional nuclear Overhauser enhancement spectroscopy (NOESY) was used to obtain a list of 202 distance constraints between individually assigned hydrogen atoms of the polypeptide chain, to identify the positions of the three disulfide bridges, and to locate the single cis peptide bond. Supplementary geometric constraints were derived from the vicinal spin-spin couplings and the locations of certain hydrogen bonds, as determined by nuclear magnetic resonance (n.m.r.). Using a new distance geometry program (DISGEO) which is capable of computing all-atom structures for proteins the size of BUSI IIA, five conformers were computed from the NOE distance constraints alone, and another five were computed with the supplementary constraints included. Comparison of the different structures computed from the n.m.r. data among themselves and with the crystal structures of two homologous proteins shows that the global features of the conformation of BUSI IIA (i.e. the overall dimensions of the molecule and the threading of the polypeptide chain) were well-defined by the available n.m.r. data. In the Appendix, we describe a preliminary energy refinement of the structure, which showed that the constraints derived from the n.m.r. data are compatible with a low energy spatial structure.

Amino Acid Sequence↗

The combinatorial distance geometry method for the calculation of molecular conformation. I. A new approach to an old problem.

A new approach to the long-standing local minimum problem of molecular energy minimization is proposed. The approach relies upon a field of computer mathematics known as combinatorial optimization, together with methods of conformational analysis derived from distance geometry. The advantages over the usual numerical techniques of optimization are, first, that the algorithms derived are globally convergent, and second, that the mathematical problems involved are well-posed and suitable for study within the modern theory of computational complexity. In this paper we introduce the approach, and describe a computer program based on it.

Computers↗

The combinatorial distance geometry method for the calculation of molecular conformation. II. Sample problems and computational statistics.

The performance of a branch and bound algorithm for molecular energy minimization is evaluated on a variety of test problems. Although not at present efficient enough for use in most practical situations, we show that it has distinct advantages over more conventional methods of global minimization. In addition, this study illustrates the technique on which the present algorithm is based, and the problems which must be overcome in developing an efficient algorithm based on similar principles.

Computers↗

The sampling properties of some distance geometry algorithms applied to unconstrained polypeptide chains: a study of 1830 independently computed conformations.

In this paper we study the statistical geometry of ensembles of poly (L-alanine) conformations computed by several different distance geometry algorithms. Since basic theory only permits us to predict the statistical properties of such ensembles a priori when the distance constraints have a very simple form, the only constraints used for these calculations are those necessary to obtain reasonable bond lengths and angles, together with a lack of short- and long-range atomic overlaps. The geometric properties studied include the squared end-to-end distance and radius of gyration of the computed conformations, in addition to the usual rms coordinate and phi/psi angle deviations among these conformations. The distance geometry algorithms evaluated include several variations of the well-known embed algorithm, together with optimizations of the torsion angles using the ellipsoid and variable target function algorithms. The conclusions may be summarized as follows: First, the distribution with which the trial distances are chosen in most implementations of the embed algorithm is not appropriate when no long-range upper bounds on the distances are present, because it leads to unjustifiably expanded conformations. Second, chosing the trial distances independently of one another leads to a lack of variation in the degree of expansion, which in turn produces a relatively low rms square coordinate difference among the members of the ensemble. Third, when short-range steric constraints are present, torsion angle optimizations that start from conformations obtained by choosing their phi/psi angles randomly with a uniform distribution between -180 degrees and +180 degrees do not converge to conformations whose angles are uniformly distributed over the sterically allowed regions of the phi/psi plane. Finally, in an appendix we show how the sampling obtained with the embed algorithm can be substantially improved upon by the proper application of existing methodology.

Algorithms↗