We developed a molecular simulation method suitable for estimation of binding free energy, called the filling potential method, based on the concept of the Taboo search, which is a type of self-avoiding random walk consisting of a cycle of local-minimum searches and transition state searches. The filling potential method is an umbrella potential sampling method, and enables the ligand molecule to drift from its local minima automatically. In the case of the filling potential method, the umbrella potential is a combination of Gaussian-type repulsive potentials, which are located on the trajectory of the ligand. Without setting the reaction coordinates a priori, this method searches for and determines the suitable reaction coordinates by successive generation of umbrella potentials based on its trajectory analysis. The weighted histogram analysis for these trajectories gives the binding free energy of the ligand to the receptor protein. It was applied to the complex of thermolysin and its inhibitors in explicit water, and the free energy surfaces with the stable binding state and the energy barrier were examined. The calculated binding free energies agreed well with the experimental results.

No abstract

The RIKEN Mouse Gene Encyclopaedia Project, a systematic approach to determining the full coding potential of the mouse genome, involves collection and sequencing of full-length complementary DNAs and physical mapping of the corresponding genes to the mouse genome. We organized an international functional annotation meeting (FANTOM) to annotate the first 21,076 cDNAs to be analysed in this project. Here we describe the first RIKEN clone collection, which is one of the largest described for any organism. Analysis of these cDNAs extends known gene families and identifies new ones.

We have developed a novel assay system for systematic analysis of protein-protein interactions (PPIs) that is characteristic of a PCR-mediated rapid sample preparation and a high-throughput assay system based on the mammalian two-hybrid method. Using gene-specific primers, we successfully constructed the assay samples by two rounds of PCR with up to 3.6 kb from the first-round PCR fragments. In the assay system, we designed all the steps to be performed by adding only samples, reagents, and cells into 384-well assay plates using two types of semiautomatic multiple dispensers. The system enabled us examine more than 20,000 assay wells per day. We detected 145 interactions in our pilot study using 3500 samples derived from mouse full-length enriched cDNAs. Analysis of the interaction data showed both several significant interaction clusters and predicted functions of a few uncharacterized proteins. In combination with our comprehensive mouse full-length cDNA clone bank covering a large part of the whole genes, our high-throughput assay system will discover many interactions to facilitate understanding of the function of uncharacterized proteins and the molecular mechanism of crucial biological processes, and also enable completion of a rough draft of the entire PPI panel in certain cell types or tissues of mouse within a short time.

A compound's synthetic accessibility (SA) is an important aspect of drug design, since in some cases computer-designed compounds cannot be synthesized. There have been several reports on SA prediction, most of which have focused on the difficulties of synthetic reactions based on retro-synthesis analyses, reaction databases and the availability of starting materials. We developed a new method of predicting SA using commercially available compound databases and molecular descriptors. SA was estimated from the probability of existence of substructures consisting of the compound in question, the number of symmetry atoms, the graph complexity, and the number of chiral centers of the compound. The probabilities of the existence of given substructures were estimated based on a compound library. The predicted SA results reproduced the expert manual assessments with a Pearson correlation coefficient of 0.56. Since our method required a compound database and not a reaction database, it should be easy to customize the prediction for compound vendors. The correlation between the sales price of approved drugs and the SA values was also examined and found to be weak. The price most likely depends on the total cost of development and other factors.

Chemicals and toxins are useful tools to elucidate the structure-function relationship of various proteins including ion channels. The HERG channel is blocked by many compounds and this may cause life-threatening cardiac arrhythmia. Besides block, some chemicals such as the class III anti-arrhythmic agent nifekalant stimulate HERG at low potentials by shifting its activation curve towards hyperpolarizing voltages. This is called "facilitation". Here, we report mutations and simulations analyzing the association between nifekalant and channel pore residues for block and facilitation. Alanine-scanning mutagenesis was performed in the pore region of HERG. The mutations at the base of the pore helix (T623A), the selectivity filter (V625A) and the S6 helix (G648A, Y652A and F656A) abolished and S624A attenuated both block and facilitation induced by the drug. On the other hand, the mutation of other residues caused either an increase or a decrease in nifekalant-induced facilitation without affecting block. An open-state homology model of the HERG pore suggested that T623, S624, Y652 and F656 faced the central cavity, and were positioned within geometrical range for the drug to be able to interact with all of them at the same time. Of these, S649 was the only polar residue located within possible interaction distance from the drug held in its blocking position. Further mutations and flexible-docking simulations suggest that the size, but not the polarity, of the side chain at S649 is critical for drug induced facilitation.

No abstract

A molecular dynamics (MD) simulation program for biological macromolecules was implemented with a non-Ewald scheme for long-ranged electrostatic interactions and run on a general purpose graphics processing unit (GPU). We recently developed several non-Ewald methods to compute the electrostatic energies with high precision. In particular, the zero-dipole summation (ZD) method, which takes into account the neutralities of charges and dipoles in a truncated subset, enables the calculation of electrostatic interactions with high accuracy and low computational cost, and its algorithm is simple enough to be implemented in a GPU. We developed an MD program with the space decomposition algorithm, myPresto/psygene, and applied it to several biological macromolecular systems with GPUs implementing the ZD method. Rapid computing performance with high accuracy was obtained.