Automation of conformational analysis and other molecular modeling calculations

A software system has been developed for facilitating modeling calculations on large numbers of molecules. Using the system, it is possible to subject one or more molecules to a series of calculations, each requiring use of a different computer program. No user intervention is required: where necessary, output from one program is used automatically as input to the next. Names are assigned to output files automatically and in a systematic manner. As an example, the system can be used to perform a succession of calculations aimed at identifying the major low-energy conformers of each of a set of molecules, starting only from their chemical connectivities. The reliability of the results has been tested by calculations on 40 molecules taken from the Cambridge Structural Database. The observed crystal structure geometry could be found for the majority of these molecules.     

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.     

HAMOG: Molecular graphics program for chemistry 5 biochemistry 5 molecular biology and enzyme research

HAMOG is a computer graphics program written in C for personal computers. Clear menus and a contextsensitive help option make the program easy to operate for occasional users. HAMOG provides a flexible environment for displaying and manipulating molecules and molecular systems. Special functions allow the investigation of structure-activity relationships of biologically active molecules. These include the calculation of molecular electrostatic potentials and fields, the superposition of molecules and the calculation of steric accessibilities. The visualization and manipulation of protein structures immediately readable from the Brookhaven Protein Data Bank files are also possible using HAMOG. The construction of any peptide or protein structure is very simple.     

A simulation of T4 bacteriophage assembly and operation

A model is presented for the self-assembly and operation of a bacteriophage comparable with the T4 bacteriophage that infects Escherichia coli. The model treats protein molecules as simple units obeying the principle free energy minimization, and exhibiting the properties of quasi-equivalence and conformational switching. A computer program incorporating the model has been developed. The results of simulation using this program are presented.     

Automatic construction of restriction site maps.

A computer program is described which constructs maps of restriction endonuclease cleavage sites in DNA molecules, given only the fragment lengths. The program utilizes fragment length data from single and double restriction enzyme digests to generate maps for linear or circular molecules. The search for a map can be limited to the unknown (insert) region of a recombinant phage or plasmid. Typical restriction maps with four or five enzymes which cut at three to five unknown sites can be calculated in a few minutes.     

MOLPACK: molecular graphics,for studying the packing of protein molecules in the crystallographic unit cell

A graphics program, MOLPACK, has been developed on the Silicon Graphics IRIS-4D computer system for displaying the packing of proteins in the crystallographic unit cell. In addition to the normal viewing operations of rotation, translation and scaling, the program has the ability to translate molecules along the cell axes while maintaining their crystallographic equivalent positions within the unit cell. This allows the user to observe the packing of protein molecules generated by molecular replacement, to create a new packing model or to locate an unknown molecule. A special feature of the program is that up to four independent molecules can be manipulated in the asymmetric unit.     

The Search for the Spatial and Electronic Requirements of a Drug

A new program called GAMMA (genetic algorithm for multiple molecule alignment) has been developed for the superimposition of several three-dimensional chemical structures. Superimposition of molecules and evaluation of structural similarity is an important task in drug design and pharmaceutical research. Similarities of compounds are determined by this program either based on their structural or their physicochemical properties by defining different matching criteria. These matching criteria are atomic properties such as atomic number or partial atomic charges. The program is based on a combination of a genetic algorithm with a numerical optimization process. A major goal of this hybrid procedure is to address the conformational flexibility of ligand molecules adequately. Thus, only one conformation per structure is necessary and the program can work even when only one conformation of a compound is stored in a database. The genetic algorithm optimizes in a nondeterministic process the size and the geometric fit of the overlay. The geometric fit of the conformations is further improved by changing torsional angles combining the genetic algorithm and the directed tweak method. The determination of the fitness of a superimposition is based on the Pareto optimization. As an application the superimposition of a set of Cytochrome P450c17 enzyme inhibitors has been performed.Electronic Supplementary Material available.     

Display and interpretation of solvent electron density distributions in insulin crystals

In macromolecular crystallography, three-dimensional contour surfaces are useful for interactive computer graphics displays of the protein electron density but are less effective for presenting static images of large volumes of solvent density. A raster-based computer graphics program which displays depth-cued projections of continuous density distributions has been developed to analyze the distribution of solvent atoms in macromolecular crystals. Maps of the water distribution in the cubic insulin crystal show some well-ordered waters, which are bound to surrounding protein atoms by multiple hydrogen bonds, and an ill-defined solvent structure at a greater distance from the protein surface. Molecular dynamics calculations were used to assist in the interpretation of the time-varying solvent structure within two enclosed cavities in the crystal. Two water molecules that ligate a sodium ion were almost immobile during the simulation but the majority of water molecules were found to move rapidly between the density maxima identified from the crystallographic refinement.     

Analysis of structural water and CH···π interactions in HIV-1 protease and PTP1B complexes using a hydrogen bond prediction tool,HBPredicT

A hydrogen bond prediction tool HBPredicT is developed for detecting structural water molecules and CH···π interactions in PDB files of protein-ligand complexes. The program adds the missing hydrogen atoms to the protein, ligands, and oxygen atoms of water molecules and subsequently all the hydrogen bonds in the complex are located using specific geometrical criteria. Hydrogen bonds are classified into various types based on (i) donor and acceptor atoms, and interactions such as (ii) protein-protein, (iii) protein-ligand, (iv) protein-water, (v) ligand-water, (vi) water-water, and (vii) protein-water-ligand. Using the information in category (vii), the water molecules which form hydrogen bonds with the ligand and the protein simultaneously–the structural water–is identified and retrieved along with the associated ligand and protein residues. For CH···π interactions, the relevant portions of the corresponding structures are also extracted in the output. The application potential of this program is tested using 19 HIV-1 protease and 11 PTP1B inhibitor complexes. All the systems showed presence of structural water molecules and in several cases, the CH···π interaction between ligand and protein are detected. A rare occurrence of CH···π interactions emanating from both faces of a phenyl ring of the inhibitor is identified in HIV-1 protease 1D4L.     

Ab initio Molecular Electrostatic Potential Grid Maps for Quantitative Similarity Calculations of Organic Compounds

The program MolSim designed to calculate the similarity of different molecules quantitatively in a fast and easy way is described. The molecular similarity is estimated for the molecular shape as well as for the electrostatic potentials of the molecules derived from ab initio calculations. A grid-based method is used to determine the steric and electrostatic similarities between a lead compound and the corresponding test set by calculating the Spearman correlation coefficient. The superpositioning of the molecules was accomplished with the SEAL algorithm incorporating a Monte Carlo simulated annealing approach while preserving the conformational flexibility of the calculated structures.The ability of the program was tested on a set of Sandalwood odour compounds, a class of substances that is difficult to analyse with respect to its structure-activity relationship because of the structural diversity of Sandalwood odour compounds, in contrast to their high selectivity and pronounced structural specificity. The application of the program on a small test set of these compounds showed that the program is able to explain the Sandalwood odour activity correctly.Electronic Supplementary Material available.     

A simple method to display molecular orbitals with computer graphics

A computer program named LOBE was developed to draw molecular orbitals as lobes on a graphic display. With this program, any molecular orbital of large molecules can be displayed quickly. This program is suitable not only for general-purpose computers but also for microcomputers. A sample application is used to illustrate the program.     

Electron microscopy and computerized evaluation of some partially denatured group P resistance plasmids.

DNA of the R plasmids RP1, RP4 and RP8 was isolated from various hosts. The lengths of these plasmid molecules were determined by electron microscopy: RP1 and RP4 were about 19 micron long, RP8 measured 31 micron. An RP4 plasmid mutant, designated RP4a, was isolated from Escherichai coli; it was about 1 micron shorter than normal RP4 DNA. To investigate the molecular relationship between RP4, RP4a and RP8 DNAs of these plasmids were partially denatured and examined in the electron microscope. Measurements of the length and denaturation pattern of the DNA molecules were used to construct physical maps. A new computer program was devised for the alignment of the circular molecules, and the effect of variations of different parameters on the reliability of the program was tested. A comparison of the denaturation pattern of RP4 and RP8 indicated that RP8 was composed of total RP4 plus an additional DNA fragment. The RP4a mutant plasmid could be defined as a deletion mutant with loss of 1 micron DNA.     

MULTIHYDRO and MONTEHYDRO: conformational search and Monte Carlo calculation of solution properties of rigid or flexible bead models

A computer program, MULTIHYDRO, has been constructed for the calculation of hydrodynamic coefficients and other solution properties of multiple possible conformations of a bead model. With minimal additional programming to describe the model under study, this program interfaces efficiently with HYDRO for the calculation of solution properties, including hydrodynamic coefficients, radius of gyration, covolume, etc. A useful application is the conformation search of rigid macromolecules, because many possible conformations can be evaluated in a single run of the program. In this paper we also pay attention to the properties of flexible macromolecules, in the so-called Monte Carlo rigid-body approximation, which is virtually exact for the simpler solution properties. The theoretical aspects of the procedure are described, and we show how MULTIHYDRO can be employed for this calculation. However, for flexible molecules, a more general simulation scheme is importance-sampling Monte Carlo generation. We describe how this procedure is implemented in another computer program, MONTEHYDRO. Examples of the usage of these tools are provided.     

Investigating protein-protein interaction surfaces using a reduced stereochemical and electrostatic model

A method of calculating the electrostatic potential energy between two molecules, using finite difference potential, is presented. A reduced charge set is used so that the interaction energy can be calculated as the two static molecules explore their full six-dimensional configurational space. The energies are contoured over surfaces fixed to each molecule with an interactive computer graphics program. For two crystal structures (trypsin-trypsin inhibitor and anti-lysozyme Fab-lysozyme), it is found that the complex corresponds to highly favourable interacting regions in the contour plots. These matches arise from a small number of protruding basic residues interacting with enhanced negative potential in each case. The redox pair cytochrome c peroxidase-cytochrome c exhibits an extensive favourably interacting surface within which a possible electron transfer complex may be defined by an increased electrostatic complementarity, but a decreased electrostatic energy. A possible substrate transfer configuration for the glycolytic enzyme pair glyceraldehyde phosphate dehydrogenase-phosphoglycerate kinase is presented.     

Harnessing biosynthesis and quantitative NMR for late stage functionalization of lead molecules: Application to the M1 positive allosteric modulator (PAM) program

A facile method for late stage diversification of lead molecules for the M1 PAM program using biosynthesis is described. Liver microsomes from several species are screened to identify a high turnover system. Subsequent incubations using less than 1?mg of substrate generate nanomole quantities of drug metabolites that are purified, characterized by microcryoprobe NMR spectroscopy, and quantified to known concentrations to enable rapid biology testing. The late-stage diversification of lead compounds provides rapid SAR feedback to the medicinal chemistry design cycle.     

Program for the visualization of inorganic crystals

A program has been developed to manipulate images of inorganic structures and organic molecules on ALLIANT VFX/40 using the PHIGS + standard. This article reviews algorithms for representing spheres, ellipsoids and various polyhedrons involved in inorganic chemistry. The program also supports the display and manipulation of animated frames from dynamics simulations. Many graphical facilities have been implemented and we discuss their interest in the field of molecular graphics.     

The calculation of sedimentation coefficients in vertical rotors

A program for the calculation of sedimentation coefficients of molecules centrifuged in sucrose in vertical rotors has been developed. The program has been tested with both protein and RNA of known sedimentation coefficients. The preparation can accept any shape of gradient in the 0–70% sucrose and any temperature in the range of 0–60°C. The program can be used with any vertical rotor for which the dimensions are known.     

The influence of N-dialkyl and other cationic substituents on DNA intercalation and genotoxicity

DNA intercalation by small chemical molecules can result in frameshift mutagenesis and chromosomal breakage. With evidence mounting that broadly diverse structures are capable of intercalating between DNA base pairs, it becomes important to better define those structural features that enhance intercalation strength and those that confer genotoxicity particularly among those intercalators that do not have the classical planar tricyclic fused ring structure. A chemical substituent that is present on many pharmaceutical and other biologically active molecules is the N-dialkyl group. In the present study, we investigate if and how the presence of an aromatic N-dialkyl or other cationic group affects the genotoxicity and DNA intercalation ability of 26 selected acridines, phenothiazines, benzophenones, triphenylethylenes and other classes of molecules. The data were obtained from the literature, from experiments using a cell-based DNA intercalation assay, and from modeling studies using a three-dimensional computational DNA docking program. It is demonstrated that cationic substitution can enhance both genotoxicity and electrostatic interactions within a chemical/DNA intercalation complex.