Search PubMed⌕ Search

Biomedical subjects

M F Sanner

Publications and source records attributed to M F Sanner.

7 recordsLinked to original sources
PubMed1999

Integrating computation and visualization for biomolecular analysis: an example using python and AVS.

One of the challenges in biocomputing is to enable the efficient use of a wide variety of fast-evolving computational methods to simulate, analyze, and understand the complex properties and interactions of molecular systems. Our laboratory investigates several areas including molecular visualization, protein-ligand docking, protein-protein docking, molecular surfaces, and the derivation of phenomenological potentials. In this paper we present an approach based on the Python programming language to achieve a high level of integration between these different computational methods and our primary visualization system AVS. This approach removes many limitations of AVS while increasing dramatically the inter-operability of our computational tools. Several examples are shown to illustrate how this approach enables a high level of integration and inter-operability between different tools, while retaining modularity and avoiding the creation of a large monolithic package that is difficult to extend and maintain.

Capsid↗
PubMed1997

Residue-residue mean-force potentials for protein structure recognition.

We present two new sets of energy functions for protein structure recognition, given the primary sequence of amino acids along the polypeptide chain. The first set of potentials is based on the positions of alpha- and the second on positions of beta- and alpha-carbon atoms of amino acid residues. The potentials are derived using a theory of Boltzmann-like statistics of protein structure. The energy terms incorporate both long-range interactions between residues remote along a chain and short-range interactions between near neighbors. Distance dependence is approximated by a piecewise constant function defined on intervals of equal size. The size of the interval is optimized to preserve as much detail as possible without introducing excessive error due to limited statistics. A database of 214 non-homologous proteins was used both for the derivation of the potentials, and for the 'threading' test originally suggested by Hendlich et al. (1990) J. Mol. Biol., 216, 167-180. Special care is taken to avoid systematic error in this test. For threading, we used 100 non-homologous protein chains of 60-205 residues. The energy of each of the native structures was compared with the energy of 43,000 to 19,000 alternative structures generated by threading. Of these 100 native structures, 92 have the lowest energy with alpha-carbon-based potentials and, even more, 98 of these 100 structures, have the lowest energy with the beta- and alpha-carbon based potentials.

Databases, Factual↗
PubMed1997

Accurate mean-force pairwise-residue potentials for discrimination of protein folds.

We present two new sets of energy functions for protein structure recognition. The first set of potentials is based on the positions of alpha- and the second on positions of beta-carbon atoms of amino acid residues. The potentials are derived using a theory of Boltzmann-like statistics of protein structure by Finkelstein et al. The energy terms incorporate both long-range interactions between residues remote along a chain and short-range interactions between near neighbors. Distance-dependence is approximated by a piecewise constant function defined on intervals of equal size. The size of this interval is optimized. A database of 222 non-homologous proteins was used both for the derivation of the potentials, and for the "threading" test originally suggested by Hendlich et al. For threading, we used 102 non-homologous protein chains of 60 to 200 residues. The energy of each of the native structures was compared with the energy of 45 to 20 thousand alternative structures generated by threading. Of these 102 native structures 94 have the lowest energy with alpha-carbon-based potentials, and even more, 100 of these 102 structures, have the lowest energy with the beta-carbon-based potentials.

Computer Simulation↗
PubMed1997

Real time surface reconstruction for moving molecular fragments.

Recently we introduced the Reduced Surface as an efficient tool to built molecular surfaces. We describe here how this geometric construct can be used to efficiently reconstruct the solvent excluded surface of a protein for which the coordinates of a subset of atoms are changing. We show that, the complexity of that operation is not dependent upon the size of the molecule and is in O[tlog(t)] where t is the maximum of the number of probes and atoms involved in the reconstruction of the surface. The algorithms described here have been implemented and tested on several proteins. The triangulation of the solvent excluded surface of proteins in which a side chain was changing conformation could be updated at rates ranging from 7 to 22 frames per second. We also applied this method to compute the surface area fluctuation of the FIV protease undergoing a constrained molecular dynamics simulation (16 mobile residues). Rate of 6 frames per second were obtained in this case.

Algorithms↗
PubMed1996

Adjusting potential energy functions for lattice models of chain molecules.

Lattice models of proteins can approximate off-lattice structure to arbitrary precision with RMS (root mean squared) deviations roughly equal to half the lattice spacing (Rykunov et al., Proteins 22:100-109, 1995; Reva et al., J. Comp. Biol., 1996). However, even small distortions in the positions of chain links lead to significant errors in lattice-based energy calculations (Reva et al., J. Comp. Chem., 1996). These errors arise mainly from rigid interactions (such as steric repulsion) which change their energies considerably at a range which is much smaller than the usual accuracy of lattice modeling (> 1.0 A). To reduce this error, we suggest a procedure of adjusting energy functions to a given lattice. The general approach is illustrated with energy calculations based on pairwise potentials by Kolinski et al. (J. Chem. Phys. 98:1-14, 1993). At all the lattice spacings, from 0.5-3.8 A, the lattice-adjusted potentials improve the accuracy of lattice-based energy calculations and increase the correlations between off-lattice and lattice energies.

Models, Chemical↗
PubMed1996

Reduced surface: an efficient way to compute molecular surfaces.

Because of their wide use in molecular modeling, methods to compute molecular surfaces have received a lot of interest in recent years. However, most of the proposed algorithms compute the analytical representation of only the solvent-accessible surface. There are a few programs that compute the analytical representation of the solvent-excluded surface, but they often have problems handling singular cases of self-intersecting surfaces and tend to fail on large molecules (more than 10,000 atoms). We describe here a program called MSMS, which is shown to be fast and reliable in computing molecular surfaces. It relies on the use of the reduced surface that is briefly defined here and from which the solvent-accessible and solvent-excluded surfaces are computed. The four algorithms composing MSMS are described and their complexity is analyzed. Special attention is given to the handling of self-intersecting parts of the solvent-excluded surface called singularities. The program has been compared with Connolly's program PQMS [M.L. Connolly (1993) Journal of Molecular Graphics, Vol. 11, pp. 139-141] on a set of 709 molecules taken from the Brookhaven Data Base. MSMS was able to compute topologically correct surfaces for each molecule in the set. Moreover, the actual time spent to compute surfaces is in agreement with the theoretical complexity of the program, which is shown to be O[n log(n)] for n atoms. On a Hewlett-Packard 9000/735 workstation, MSMS takes 0.73 s to produce a triangulated solvent-excluded surface for crambin (1 crn, 46 residues, 327 atoms, 4772 triangles), 4.6 s for thermolysin (3tln, 316 residues, 2437 atoms, 26462 triangles), and 104.53 s for glutamine synthetase (2gls, 5676 residues, 43632 atoms, 476665 triangles).

Chemistry, Physical↗