GoBean: a Java GUI application for visual exploration of GO term enrichments

We have developed a biologist-friendly, Java GUI application (GoBean) for GO term enrichment analysis. It was designed to be a comprehensive and flexible GUI tool for GO term enrichment analysis, combining the merits of other programs and incorporating extensive graphic exploration of enrichment results. An intuitive user interface with multiple panels allows for extensive visual scrutiny of analysis results. The program includes many essential and useful features, such as enrichment analysis algorithms, multiple test correction methods, and versatile filtering of enriched GO terms for more focused analyses. A unique graphic interface reflecting the GO tree structure was devised to facilitate comparisons of multiple GO analysis results, which can provide valuable insights for biological interpretation. Additional features to enhance user convenience include built in ID conversion, evidence code-based gene-GO association filtering, set operations of gene lists and enriched GO terms, and user -provided data files. It is available at http://neon.gachon.ac.kr/GoBean/.     

Experience in simulating the structural and dynamic features of small proteins using table supercomputers

The results of theoretical studies of the structural and dynamic features of peptides and small proteins have been presented that were carried out by quantum chemical and molecular dynamics methods in high-performance graphic stations, “table supercomputers,” using distributed calculations by the CUDA technology.     

XtalView/Xfit--A versatile program for manipulating atomic coordinates and electron density

Xfit is a model-building and map viewing program in XtalView that is used by the structural biology community including researchers in the fields of crystallography, molecular modeling, and electron microscopy. Among its distinguishing features are built-in fast Fourier transforms that allow users flexibility in map calculations including the creation of OMIT maps and the updating of structure factors to reflect model changes from within the program. Written in C and using the freely available XView toolkit, it is highly portable to almost any X-windows based workstation including Intel-based LINUX systems. Its user interface is designed to aid in facile model-building and contains a semiautomated fitting system that allows the user to interactively and rapidly build chain de novo into an electron density map. The program is highly optimized to allow such features as interactive contour levels and map calculations to be completed within a few seconds. Features in the latest version including phase-combination, solvent-flattening, automated water addition, and small-probe dot contact surfaces, as well as basic design features, are discussed.     

The World Wide Web as a graphical user interface to program macros for molecular graphics,molecular modeling,and structure-based drug design

In this article we describe how the World Wide Web (WWW or Web) has been employed to provide access to computational chemistry software and protein structure data via program macros. We show how the combination of Web technology and macros can automate both the running of chemistry software and the execution of complex operations on protein structures. The current version of the system supports the molecular visualization packages GRASP,¹ RASMOL,² MOLVIEWER-OGL³ and INSIGHT95,⁴ and the ligand design tool GRID⁵ and includes more than 175 in-house protein-ligand complexes. The approach enables in-experienced users to confidently make full use of sophisticated modeling techniques by offering only sensible options, hiding parameter settings, and controlling program invocation and macro excution. Our interface provides both the expert and non-expert alike with powerful tools for protein structure visualization, molecular modeling, and rational drug design.     

Mechanisms for Two-Step Proton Transfer Reactions in the Outward-Facing Form of MATE Transporter

Bacterial pathogens or cancer cells can acquire multidrug resistance, which causes serious clinical problems. In cells with multidrug resistance, various drugs or antibiotics are extruded across the cell membrane by multidrug transporters. The multidrug and toxic compound extrusion (MATE) transporter is one of the five families of multidrug transporters. MATE from Pyrococcus furiosus uses H⁺ to transport a substrate from the cytoplasm to the outside of a cell. Crystal structures of MATE from P. furiosus provide essential information on the relevant H⁺-binding sites (D41 and D184). Hybrid quantum mechanical/molecular mechanical simulations and continuum electrostatic calculations on the crystal structures predict that D41 is protonated in one structure (Straight) and, both D41 and D184 protonated in another (Bent). All-atom molecular dynamics simulations suggest a dynamic equilibrium between the protonation states of the two aspartic acids and that the protonation state affects hydration in the substrate binding cavity and lipid intrusion in the cleft between the N- and C-lobes. This hypothesis is examined in more detail by quantum mechanical/molecular mechanical calculations on snapshots taken from the molecular dynamics trajectories. We find the possibility of two proton transfer (PT) reactions in Straight: the 1st PT takes place between side-chains D41 and D184 through a transient formation of low-barrier hydrogen bonds and the 2nd through another H⁺ from the headgroup of a lipid that intrudes into the cleft resulting in a doubly protonated (both D41 and D184) state. The 1st PT affects the local hydrogen bond network and hydration in the N-lobe cavity, which would impinge on the substrate-binding affinity. The 2nd PT would drive the conformational change from Straight to Bent. This model may be applicable to several prokaryotic H⁺-coupled MATE multidrug transporters with the relevant aspartic acids.     

Construction and optimisation of a computer model for a bacterial membrane

The main steps in the construction of a computer model for a bacterial membrane are described. The membrane has been built of 72 lipid molecules, 54 of which being 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylethanolamine (POPE) and 18--1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidyl-rac-glycerol (POPG) molecules (thus in the proportion of 3:1). The membrane was hydrated with 1955 water molecules (approximately 27 water molecules per lipid). To neutralise the electronic charge (-e) on each POPG molecule, 18 sodium ions (Na+) were added to the membrane close to the POPG phosphate groups. The atomic charges on the POPE and POPG headgroups were obtained from ab initio quantum mechanical restrained electrostatic potential fitting (RESP) (Bayly et al., 1993, J. Phys. Chem. 97, 10269) using the GAMESS program at the 6-31G* level (Schmidt et al., 1993, J. Comput. Chem. 14, 1347). The model constructed in this way provided an initial structure for subsequent molecular dynamics simulation studies intended to elucidate the atomic level interactions responsible for the structure and dynamics of the bacterial membrane.     

COBRAMM 2.0 — A software interface for tailoring molecular electronic structure calculations and running nanoscale (QM/MM) simulations

We present a new version of the simulation software COBRAMM, a program package interfacing widely known commercial and academic software for molecular modeling. It allows a problem-driven tailoring of computational chemistry simulations with effortless ground and excited-state electronic structure computations. Calculations can be executed within a pure QM or combined quantum mechanical/molecular mechanical (QM/MM) framework, bridging from the atomistic to the nanoscale. The user can perform all necessary steps to simulate ground state and photoreactions in vacuum, complex biopolymer, or solvent environments. Starting from ground-state optimization, reaction path computations, initial conditions sampling, spectroscopy simulation, and photodynamics with deactivation events, COBRAMM is designed to assist in characterization and analysis of complex molecular materials and their properties. Interpretation of recorded spectra range from steady-state to time-resolved measurements. Various tools help the user to set up the system of interest and analyze the results.     

Electric Field Distribution of an Optical Nanocavity Embedded with a Single Molecule

Using state-of-the-art quantum mechanical calculations, we investigate the spatial profile of the electric field enhancements induced within an optical cavity embedded with a variety of organic molecules. We observed marked differences in the spectral positions, spectral intensities, maximum achievable electric field, and spatial profile of the electric field with a remarkable sensitivity to the embedded molecular type. In a broader perspective, our quantum mechanical calculations provide quantitative access to the molecule-dependent electric field distributions and unveil intricate and rich optical features.

     

EzMol: A Web Server Wizard for the Rapid Visualization and Image Production of Protein and Nucleic Acid Structures

EzMol is a molecular visualization Web server in the form of a software wizard, located at http://www.sbg.bio.ic.ac.uk/ezmol/. It is designed for easy and rapid image manipulation and display of protein molecules, and is intended for users who need to quickly produce high-resolution images of protein molecules but do not have the time or inclination to use a software molecular visualization system. EzMol allows the upload of molecular structure files in PDB format to generate a Web page including a representation of the structure that the user can manipulate. EzMol provides intuitive options for chain display, adjusting the color/transparency of residues, side chains and protein surfaces, and for adding labels to residues. The final adjusted protein image can then be downloaded as a high-resolution image. There are a range of applications for rapid protein display, including the illustration of specific areas of a protein structure and the rapid prototyping of images.     

Activated Complex Theory of Barite Scale Control Processes

A theoretical framework useful for the estimation of scale dissolution rate constants is introduced. The model consists of (a) a dissolution mechanism, (b) a quantum mechanical force field capable of describing the bulk and surface properties of barite, (c) a complete study of the structure of the dissolver molecule complex in solution, (d) a three dimensional periodic system useful for mapping the barite-dissolver interactions, including the localization of the activated complex and (e) a rate expression to estimate the dissolution rate constants from the properties of the activated complex. Our results show that molecular modeling, through a combination of molecular mechanics and high level quantum mechanical calculations, provide a new and insightful information about scale control processes. This revised version was published online in June 2006 with corrections to the Cover Date.     

Approximate quantum trajectory dynamics for reactive processes in condensed phase

A method of molecular dynamics with quantum corrections, practical for studies of large molecular systems, is reviewed. The approach is based on the Bohmian formulation of the time-dependent Schrödinger equation in which a wavefunction is represented by an ensemble of interdependent trajectories. The quantum effects come from the quantum potential acting on trajectories on par with the usual classical potential. The quantum potential is determined from the evolving nuclear wavefunction, i.e. from the quantum trajectory (QT) ensemble itself. For practical and conceptual reasons the quantum potential and corresponding quantum nuclear effect are computed only for the selected light nuclei. For studies of reactive chemical processes, the classical potential is computed on-the-fly using the density functional tight binding method of electronic structure. A massively parallel implementation, based on the message passing interface allows for efficient simulations of ensembles of thousands of trajectories describing systems of up to 200 atoms. As a biochemical application, the approximate QT approach is used to model the tunnelling-dominated proton transfer in soybean-lipoxygenase-1. A materials science application is represented by a study of the nuclear quantum effect on adsorption of hydrogen and deuterium on a C₃₇H₁₅ molecule, which is a model ‘flake’ of graphene.     

A user-friendly biological workstation

Learning methods developed by artificial intelligence research teams are very efficient for biological sequences analysis but they need running on large computers accessed by terminals. These computers are interfaced with standard displays involving long and unpleasant alphanumerical data handling. The "biological work station" is a personal computer with a color graphic screen providing a user-friendly interface for the artificial intelligence learning programs running on large computers. It provides to biologist a graphical convenient tool for sequence analysis built with efficient man-machine communication methods such as multiwindows, icons and mouse selection. It allows the biologist to edit and display sequences in an efficient and natural way, showing off directly on color pictures the data and the results of learning programs.     

GRASS: a server for the graphical representation and analysis of structures

GRASS (Graphical Representation and Analysis of Structures Server), a new web-based server, is described. GRASS exploits many of the features of the GRASP program and is designed to provide interactive molecular graphics and quantitative analysis tools with a simple interface over the World-Wide Web. Using GRASS, it is now possible to view many surface features of biological macromolecules on either standard workstations used in macromolecular analysis or personal computers. The result is a World-Wide Web-based, platform-independent, easily used tool for macromolecular visualization and structure analysis.     

Homology modeling, molecular dynamics and QM/MM study of the regulatory protein PhoP from Corynebacterium pseudotuberculosis

Corynebacterium pseudotuberculosis is a facultatively intracellular Gram-positive bacterium that causes caseous lymphadenitis, principally in sheep and goats, though sometimes in other species of animals, leading to considerable economic losses. This pathogen has a TCS known as PhoPR, which consists of a sensory histidine kinase protein (PhoR) and an intracellular response regulator protein (PhoP). This system is involved in the regulation of proteins present in various processes, including virulence. The regulation is activated by PhoP protein phosphorylation, an event that requires a magnesium (Mg²⁺) ion. Here we describe the 3D structure of the regulatory response protein (PhoP) of C. pseudotuberculosis through molecular modeling by homology. The model generated provides the first structural information on a full-length member of the OmpR/PhoP subfamily. Classical molecular dynamics was used to investigate the stability of the model. In addition, we used quantum mechanical/molecular mechanical techniques to perform (internal, potential) energy optimizations to determine the interaction energy between the Mg²⁺ ion and the structure of the PhoP protein. Analysis of the interaction energy residue by residue shows that Asp-16 and Asp-59 play an important role in the protein–Mg²⁺ ion interactions. These results may be useful for the future development of a new vaccine against tuberculosis based on genetic attenuation via a point mutation that results in the polar residue Asp-16 and/or Asp-59 being replaced with a nonpolar residue in the DNA-binding domain of PhoP of C. pseudotuberculosis.     

Conjugation in multi-tetrazole derivatives: a new design direction for energetic materials

Multi-tetrazole derivatives with conjugated structures were designed and investigated in this study. Using quantum chemistry methods, the crystal structures, electrostatic potentials (ESPs), multicenter bond orders, HOMO–LUMO energy gaps, and detonation properties of the derivatives were calculated. As expected, these molecules with conjugated structures showed low energies of their crystal structures, molecular layering in their crystals, high average ESPs, high multicenter bond order values, and enhanced detonation properties. The derivative 1,2-di(1H-tetrazol-5-yl)diazene (N2) was predicted to have the best density (1.87 g/cm³), detonation velocity (9006 m/s), and detonation pressure (36.8 GPa) of the designed molecules, while its total crystal energy was low, suggesting that it is relatively stable. Its sensitivity was also low, as the molecular stacking that occurs in its crystal allows external forces to be dissipated into movements of crystal layers. Finally, its multicenter bond order was high, indicating a highly conjugated structure.     

Differences in muscle function during walking and running at the same speed

Individual muscle contributions to body segment mechanical energetics and the functional tasks of body support and forward propulsion in walking and running at the same speed were quantified using forward dynamical simulations to elucidate differences in muscle function between the two different gait modes. Simulations that emulated experimentally measured kinesiological data of young adults walking and running at the preferred walk-to-run transition speed revealed that muscles use similar biomechanical mechanisms to provide support and forward propulsion during the two tasks. The primary exception was a decreased contribution of the soleus to forward propulsion in running, which was previously found to be significant in walking. In addition, the soleus distributed its mechanical power differently to individual body segments between the two gait modes from mid- to late stance. In walking, the soleus transferred mechanical energy from the leg to the trunk to provide support, but in running it delivered energy to both the leg and trunk. In running, earlier soleus excitation resulted in it working in synergy with the hip and knee extensors near mid-stance to provide the vertical acceleration for the subsequent flight phase in running. In addition, greater power output was produced by the soleus and hip and knee extensors in running. All other muscle groups distributed mechanical power among the body segments and provided support and forward propulsion in a qualitatively similar manner in both walking and running.