Search

[#17491] article suite à une conférence orale: 13th EG Euroworkshop on 3D object retrieval, 3DOR 2020, Graz, Austria, september 4-5, 2020; International audience; Proteins are natural modular objects usually composed of several domains, each domain bearing a specific function that is mediated through its surface, which is accessible to vicinal molecules. This draws attention to an understudied characteristic of protein structures: surface, that is mostly unexploited by protein structure comparison methods. In the present work, we evaluated the performance of six shape comparison methods, among which three are based on machine learning, to distinguish between 588 multi-domain proteins and to recreate the evolutionary relationships at the protein and species levels of the SCOPe database. The six groups that participated in the challenge submitted a total of 15 sets of results. We observed that the performance of all the methods significantly decreases at the species level, suggesting that shape-only protein comparison is challenging for closely related proteins. Even if the dataset is limited in size (only 588 proteins are considered whereas more than 160,000 protein structures are experimentally solved), we think that this work provides useful insights into the current shape comparison methods performance, and highlights possible limitations to large-scale applications due to the computational cost. (C) 2020 The Author(s). Published by Elsevier Ltd.

Many important functions in a cell are carried out by protein complexes with more than two subunits. Similar to the folding of a single protein, multimeric protein complexes in general follow an energetically favored assembly path. Knowing the assembly path not only provides critical information about the molecular mechanism of the assembly but also serves as a foundation for artificial design of protein complexes, as well as development of drugs that interfere with complex formation. There are experimental approaches for determining the assembly path of a complex; however, such methods are resource intensive. We have recently developed a computational method, Path-LZerD, which predicts the assembly path of a complex by simulating the docking process of the complex. Here, we explain how to use the Path-LZerD software with examples.

We report the performance of the protein docking prediction pipeline of our group and the results for Critical Assessment of Prediction of Interactions (CAPRI) rounds 38-46. The pipeline integrates programs developed in our group as well as other existing scoring functions. The core of the pipeline is the LZerD protein-protein docking algorithm. If templates of the target complex are not found in PDB, the first step of our docking prediction pipeline is to run LZerD for a query protein pair. Meanwhile, in the case of human group prediction, we survey the literature to find information that can guide the modeling, such as protein-protein interface information. In addition to any literature information and binding residue prediction, generated docking decoys were selected by a rank aggregation of statistical scoring functions. The top 10 decoys were relaxed by a short molecular dynamics simulation before submission to remove atom clashes and improve side-chain conformations. In these CAPRI rounds, our group, particularly the LZerD server, showed robust performance. On the other hand, there are failed cases where some other groups were successful. To understand weaknesses of our pipeline, we analyzed sources of errors for failed targets. Since we noted that structure refinement is a step that needs improvement, we newly performed a comparative study of several refinement approaches. Finally, we show several examples that illustrate successful and unsuccessful cases by our group.

Cryo-EM is a powerful method for determining protein structures. But it requires computational assistance. Physics-based computations have the power to give low-free-energy structures and ensembles of populations, but have been computationally limited to only small soluble proteins. Here, we introduce CryoFold. By integrating data of varying sparsity from electron density maps of 3–5 Å resolution with coarse-grained physical knowledge of secondary and tertiary interactions, CryoFold determines ensembles of protein structures directly from sequence. We give six examples showing its broad capabilities, over proteins ranging from 72 to 2000 residues, including membrane and multi-domain proteins, and including results from two EMDB competitions. The ensembles CryoFold predicts starting from the density data of a single known protein conformation encompass multiple low-energy conformations, all of which are experimentally validated and biologically relevant.

Motivation Many important cellular processes involve physical interactions of proteins. Therefore, determining protein quaternary structures provide critical insights for understanding molecular mechanisms of functions of the complexes. To complement experimental methods, many computational methods have been developed to predict structures of protein complexes. One of the challenges in computational protein complex structure prediction is to identify near-native models from a large pool of generated models. Results We developed a convolutional deep neural network-based approach named DOcking decoy selection with Voxel-based deep neural nEtwork (DOVE) for evaluating protein docking models. To evaluate a protein docking model, DOVE scans the protein–protein interface of the model with a 3D voxel and considers atomic interaction types and their energetic contributions as input features applied to the neural network. The deep learning models were trained and validated on docking models available in the ZDock and DockGround databases. Among the different combinations of features tested, almost all outperformed existing scoring functions. Availability and implementation Codes available at http://github.com/kiharalab/DOVE, http://kiharalab.org/dove/. Supplementary information Supplementary data are available at Bioinformatics online.

An increasing number of biological macromolecules have been solved with cryo-electron microscopy (cryo-EM). Over the past few years, the resolutions of density maps determined by cryo-EM have largely improved in general. However, there are still many cases where the resolution is not high enough to model molecular structures with standard computational tools. If the resolution obtained is near the empirical borderline (3-4 Angstroms), a small improvement of resolution will significantly facilitate structure modeling. Here, we report SuperEM, a novel deep learning-based method that uses a three-dimensional generative adversarial network for generating an improved-resolution EM map from an experimental EM map. SuperEM is designed to work with EM maps in the resolution range of 3 Angstroms to 6 Angstroms and has shown an average resolution improvement of 1.0 Angstrom on a test dataset of 36 experimental maps. The generated super-resolution maps are shown to result in better structure modelling of proteins.

We report the performance of protein-protein docking predictions by our group for recent rounds of the Critical Assessment of Prediction of Interactions (CAPRI), a community-wide assessment of state-of-the-art docking methods. Our prediction procedure uses a protein-protein docking program named LZerD developed in our group. LZerD represents a protein surface with 3D Zernike descriptors (3DZD), which are based on a mathematical series expansion of a 3D function. The appropriate soft representation of protein surface with 3DZD makes the method more tolerant to conformational change of proteins upon docking, which adds an advantage for unbound docking. Docking was guided by interface residue prediction performed with BindML and cons-PPISP as well as literature information when available. The generated docking models were ranked by a combination of scoring functions, including PRESCO, which evaluates the native-likeness of residues' spatial environments in structure models. First, we discuss the overall performance of our group in the CAPRI prediction rounds and investigate the reasons for unsuccessful cases. Then, we examine the performance of several knowledge-based scoring functions and their combinations for ranking docking models. It was found that the quality of a pool of docking models generated by LZerD, that is whether or not the pool includes near-native models, can be predicted by the correlation of multiple scores. Although the current analysis used docking models generated by LZerD, findings on scoring functions are expected to be universally applicable to other docking methods. Proteins 2017; 85:513-527. © 2016 Wiley Periodicals, Inc.

Although structures determined at near-atomic resolution are now routinely reported by cryo-electron microscopy (cryo-EM), many density maps are determined at an intermediate resolution, and extracting structure information from these maps is still a challenge. We report a computational method, Emap2sec, that identifies the secondary structures of proteins (α-helices, β-sheets and other structures) in EM maps at resolutions of between 5 and 10 A. Emap2sec uses a three-dimensional deep convolutional neural network to assign secondary structure to each grid point in an EM map. We tested Emap2sec on EM maps simulated from 34 structures at resolutions of 6.0 and 10.0 A, as well as on 43 maps determined experimentally at resolutions of between 5.0 and 9.5 A. Emap2sec was able to clearly identify the secondary structures in many maps tested, and showed substantially better performance than existing methods. Emap2sec uses a deep convolutional neural network to assign protein secondary structures in intermediate-resolution (5–10 A) cryo-EM maps.

Protein-protein interactions are the cornerstone of numerous biological processes. Although an increasing number of protein complex structures have been determined using experimental methods, relatively fewer studies have been performed to determine the assembly order of complexes. In addition to the insights into the molecular mechanisms of biological function provided by the structure of a complex, knowing the assembly order is important for understanding the process of complex formation. Assembly order is also practically useful for constructing subcomplexes as a step toward solving the entire complex experimentally, designing artificial protein complexes, and developing drugs that interrupt a critical step in the complex assembly. There are several experimental methods for determining the assembly order of complexes; however, these techniques are resource-intensive. Here, we present a computational method that predicts the assembly order of protein complexes by building the complex structure. The method, named Path-LzerD, uses a multimeric protein docking algorithm that assembles a protein complex structure from individual subunit structures and predicts assembly order by observing the simulated assembly process of the complex. Benchmarked on a dataset of complexes with experimental evidence of assembly order, Path-LZerD was successful in predicting the assembly pathway for the majority of the cases. Moreover, when compared with a simple approach that infers the assembly path from the buried surface area of subunits in the native complex, Path-LZerD has the strong advantage that it can be used for cases where the complex structure is not known. The path prediction accuracy decreased when starting from unbound monomers, particularly for larger complexes of five or more subunits, for which only a part of the assembly path was correctly identified. As the first method of its kind, Path-LZerD opens a new area of computational protein structure modeling and will be an indispensable approach for studying protein complexes., Author summary Protein-protein interactions, particularly those involving multiple proteins, are the cornerstone of numerous biological processes. Although an increasing number of multi-chain protein complex structures have been determined, fewer studies have been performed to determine the assembly order of complexes. Knowing the assembly order of a complex provides insights into the process of complex formation. Assembly order is also practically useful for reconstructing and determining the structure of a subcomplex of a large protein complex. It also has important applications including designing artificial protein complexes and drugs that prevent the assembly of protein complexes. We present a computational method, Path-LZerD, which predicts the assembly order of a protein complex by simulating its assembly process. This is the first method of this kind. A strong advantage of Path-LZerD is that the assembly order can be predicted even when the overall complex structure is not known. Path-LZerD opens a new area of computational protein structure modeling and will be an indispensable approach for studying protein complexes.

Protein structure prediction has matured over years, particularly those which uses structure templates for building a model. It can build a model with correct overall conformation in cases where appropriate templates are available. Models with the correct topology can be practically useful for limited purposes that need residue-level accuracy, but further improvement of the models can allow the models to be used in tasks that need detailed structures, such as molecular replacement in X-ray crystallography or structure-based drug screening. Thus, model refinement is an important final step in protein structure prediction to bridge predictions to real-life applications. Model refinement is one of the categories in recent rounds of Critical Assessment of techniques in protein Structure Prediction (CASP) and has recently been drawing more attention due to its realized importance. Here we report our group’s performance in the refinement category in CASP12. Our method is based on inexpensive short molecular dynamics (MD) simulations in implicit solvent. Our performance in CASP12 was among the top, which was consistent with the previous round, CASP11. Our method with short MD runs achieved comparable performance with other methods that used longer simulations. Detailed analyses found that improvements typically occurred in entire regions of a structure rather than only in flexible loop regions. The remaining challenge in the structure refinement includes large conformational refinement which involves substantial motions of secondary structure elements or domains.

Disordered protein-protein interactions (PPIs), those involving a folded protein and an intrinsically disordered protein (IDP), are prevalent in the cell, including important signaling and regulatory pathways. IDPs do not adopt a single dominant structure in isolation but often become ordered upon binding. To aid understanding of the molecular mechanisms of disordered PPIs, it is crucial to obtain the tertiary structure of the PPIs. However, experimental methods have difficulty in solving disordered PPIs and existing protein-protein and protein-peptide docking methods are not able to model them. Here we present a novel computational method, IDP-LZerD, which models the conformation of a disordered PPI by considering the biophysical binding mechanism of an IDP to a structured protein, whereby a local segment of the IDP initiates the interaction and subsequently the remaining IDP regions explore and coalesce around the initial binding site. On a dataset of 22 disordered PPIs with IDPs up to 69 amino acids, successful predictions were made for 21 bound and 18 unbound receptors. The successful modeling provides additional support for biophysical principles. Moreover, the new technique significantly expands the capability of protein structure modeling and provides crucial insights into the molecular mechanisms of disordered PPIs., Author summary A substantial fraction of the proteins encoded in genomes are intrinsically disordered proteins (IDPs), which lack a single stable structure in the native state. IDPs serve many functions including mediating protein-protein interactions (PPIs). Such disordered PPIs are prevalent in important regulatory pathways, including many interactions of the tumor suppressor protein p53. To elucidate the molecular mechanisms of disordered PPIs, obtaining tertiary structure information is essential; however, they are difficult to study with experimental techniques and existing computational protein-protein and protein-peptide modeling methods are unable to model disordered PPIs. Here we present a novel computational method for modeling the structure of disordered PPIs, which is the first of this sort. The method, IDP-LZerD, is designed to follow a known biophysical picture of the mechanism of how IDPs interact with structured proteins. IDP-LZerD successfully modeled the majority of disordered PPIs tested. This technique opens up new possibilities for structural studies of IDPs and their interactions.

Community-wide blind prediction experiments such as CAPRI and CASP provide an objective measure of the current state of predictive methodology. Here we describe a community-wide assessment of methods to predict the effects of mutations on protein-protein interactions. Twenty-two groups predicted the effects of comprehensive saturation mutagenesis for two designed influenza hemagglutinin binders and the results were compared with experimental yeast display enrichment data obtained using deep sequencing. The most successful methods explicitly considered the effects of mutation on monomer stability in addition to binding affinity, carried out explicit side-chain sampling and backbone relaxation, evaluated packing, electrostatic, and solvation effects, and correctly identified around a third of the beneficial mutations. Much room for improvement remains for even the best techniques, and large-scale fitness landscapes should continue to provide an excellent test bed for continued evaluation of both existing and new prediction methodologies.

Background We comprehensively analyzed X-ray cocrystal structures of dipeptidyl peptidase IV (DPP-4) and its inhibitor to clarify whether DPP-4 alters its general or partial structure according to the inhibitor used and whether DPP-4 has a common rule for inhibitor binding. Results All the main and side chains in the inhibitor binding area were minimally altered, except for a few side chains, despite binding to inhibitors of various shapes. Some residues (Arg125, Glu205, Glu206, Tyr662 and Asn710) in the area had binding modes to fix a specific atom of inhibitor to a particular spatial position in DPP-4. We found two specific water molecules that were common to 92 DPP-4 structures. The two water molecules were close to many inhibitors, and seemed to play two roles: maintaining the orientation of the Glu205 and Glu206 side chains through a network via the water molecules, and arranging the inhibitor appropriately at the S2 subsite. Conclusions Our study based on high-quality resources may provide a necessary minimum consensus to help in the discovery of a novel DPP-4 inhibitor that is commercially useful. Electronic supplementary material The online version of this article (doi:10.1186/s12900-016-0062-8) contains supplementary material, which is available to authorized users.

Summary An increasing number of biomolecular structures are solved by electron microscopy (EM). However, the quality of structure models determined from EM maps vary substantially. To understand to what extent structure models are supported by information embedded in EM maps, we used two computational structure refinement methods to examine how much structures can be refined using a dataset of 49 maps with accompanying structure models. The extent of structure modification as well as the disagreement between refinement models produced by the two computational methods scaled inversely with the global and the local map resolutions. A general quantitative estimation of deviations of structures for particular map resolutions are provided. Our results indicate that the observed discrepancy between the deposited map and the refined models is due to the lack of structural information present in EM maps and thus these annotations must be used with caution for further applications.

In protein structure prediction, such as template-based modeling and free modeling (ab initio modeling), the step that assesses the quality of protein models is very important. We have developed a model quality assessment (QA) program United3D that uses an optimized clustering method and a simple Cα atom contact-based potential. United3D automatically estimates the quality scores (Qscore) of predicted protein models that are highly correlated with the actual quality (GDT_TS). The performance of United3D was tested in the ninth Critical Assessment of protein Structure Prediction (CASP9) experiment. In CASP9, United3D showed the lowest average loss of GDT_TS (5.3) among the QA methods participated in CASP9. This result indicates that the performance of United3D to identify the high quality models from the models predicted by CASP9 servers on 116 targets was best among the QA methods that were tested in CASP9. United3D also produced high average Pearson correlation coefficients (0.93) and acceptable Kendall rank correlation coefficients (0.68) between the Qscore and GDT_TS. This performance was competitive with the other top ranked QA methods that were tested in CASP9. These results indicate that United3D is a useful tool for selecting high quality models from many candidate model structures provided by various modeling methods. United3D will improve the accuracy of protein structure prediction.

We have devised a power function (PF) that can predict the accuracy of a three-dimensional (3D) structure model of a protein using only amino acid sequence alignments. This Power Function (PF) consists of three parts; (1) the length of a model, (2) a homology identity percent value and (3) the agreement rate between PSI-PRED secondary structure prediction and the secondary structure judgment of a reference protein. The PF value is mathematically computed from the execution process of homology search tools, such as FASTA or various BLAST programs, to obtain the amino acid sequence alignments. There is a high correlation between the global distance test-total score (GDT_TS) value of the protein model based upon the PF score and the GDT_TS(MAX) value used as an index of protein modeling accuracy in the international contest Critical Assessment of Techniques for Protein Structure Prediction (CASP). Accordingly, the PF method is valuable for constructing a highly accurate model without wasteful calculations of homology modeling that is normally performed by an iterative method to move the main chain and side chains in the modeling process. Moreover, a model with higher accuracy can be obtained by combining the models ordered by the PF score with models sorted by the size of the CIRCLE score. The CIRCLE software is a 3D-1D program, in which energetic stabilization is estimated based upon the experimental environment of each amino acid residue in the protein solution or protein crystals.

Almost all proteins express their biological functions through the structural conformation of their specific amino acid sequences. Therefore, acquiring the three-dimensional structures of proteins is very important to elucidate the role of a particular protein. We had built protein structure model databases, which is called RIKEN FAMSBASE (http://famshelp.gsc.riken.jp/famsbase/). The RIKEN FAMSBASE is a genome-wide protein structure model database that contains a large number of protein models from many organisms. The HUMAN FAMSBASE that is one part of the RIKEN FAMSBASE contains many protein models for human genes, which are significant in the pharmaceutical and medicinal fields. We have now implemented an update of the human protein modeling database consisting of 242918 constructed models against the number of 20743 human protein sequences with an improved modeling method called Full Automatic protein Modeling System Developed (FAMSD). The results of our benchmark test of the FAMSD method indicated that it has an excellent capability to pack amino acid side-chains with correct torsion angles in addition to the main-chain, while avoiding the formation of atom-atom collisions that are not found in experimental structures. This new protein structure model database for human genes, which is named HUMAN FAMSD-BASE, is open to the public as a component part of the RIKEN FAMSBASE at http://mammalia.gsc.riken.jp/human_famsd/. A significant improvement of the HUMAN FAMSD-BASE in comparison with the preceding HUMAN FAMSBASE was verified in the benchmark test of this paper. The HUMAN FAMSD-BASE will have an important impact on the progress of biological science.

Selecting the best quality model from a set of predicted structures is one of the most important aspects of protein structure prediction. We have developed model quality assessment programs that select high quality models which account for both the Calpha backbone and side-chain atom positions. The new methods are based on the consensus method with consideration of the side-chain environment of a protein structure and the secondary structure agreement. This Side-chain Environment Consensus (SEC) method is compared with the conventional consensus method, 3D-Jury (Ginalski K. et al., Bioinformatics, 19, 1015-1018 (2003)), which takes into account only the Calpha backbone atoms of the protein model. As the result, it was found that the SEC method selects the models with more accurate positioning of the side-chain atoms than the 3D-Jury method. When the SEC method was used in combination with the 3D-Jury method (3DJ+SEC), models were selected with improved quality both in the Calpha backbone and side-chain atom positions. Moreover, the CIRCLE (CCL) method (Terashi G. et al., Proteins, 69 (Suppl. 8), 98-107 (2007)) based on the 3D-1D profile score has been shown to select the best possible models that are the closest to the native structures from candidate models. Accordingly, the 3DJ+SEC+CCL method, in which CIRCLE is used after reducing the number of candidates by the 3DJ+SEC consensus method, was found to be very effective in selecting high quality models. Thus, the best method (the 3DJ+SEC+CCL method) includes the consensus approaches of the Calpha backbone and the side-chains, the secondary structure agreement and the 3D-1D profile score which corresponds to the free energy-like score in the residues of the protein model. In short, new algorithms are introduced in protein structure evaluation methods that are based on a side-chain consensus score. Additionally, in order to apply the 3DJ+SEC+CCL method and indicate the usefulness of this method, a model of human Cabin1, a protein associated with p53 function and cancer, is created using various internet modeling and alignment servers.

We participated in rounds 6-12 of the critical assessment of predicted interaction (CAPRI) contest as the SKE-DOCK server and human teams. The SKE-DOCK server is based on simple geometry docking and a knowledge base scoring function. The procedure is summarized in the following three steps: (1) protein docking according to shape complementarity, (2) evaluating complex models, and (3) repacking side-chain of models. The SKE-DOCK server did not make use of biological information. On the other hand, the human team tried various intervention approaches. In this article, we describe in detail the processes of the SKE-DOCK server, together with results and reasons for success and failure. Good predicted models were obtained for target 25 by both the SKE-DOCK server and human teams. When the modeled receptor proteins were superimposed on the experimental structures, the smallest Ligand-rmsd values corresponding to the rmsd between the model and experimental structures were 3.307 and 3.324 A, respectively. Moreover, the two teams obtained 4 and 2 acceptable models for target 25. The overall result for both the SKE-DOCK server and human teams was medium accuracy for one (Target 25) out of nine targets.

During Critical Assessment of Protein Structure Prediction (CASP7, Pacific Grove, CA, 2006), fams-ace was entered in the 3D coordinate prediction category as a human expert group. The procedure can be summarized by the following three steps. (1) All the server models were refined and rebuilt utilizing our homology modeling method. (2) Representative structures were selected from each server, according to a model quality evaluation, based on a 3D1D profile score (like Verify3D). (3) The top five models were selected and submitted in the order of the consensus-based score (like 3D-Jury). Fams-ace is a fully automated server and does not require human intervention. In this article, we introduce the methodology of fams-ace and discuss the successes and failures of this approach during CASP7. In addition, we discuss possible improvements for the next CASP.

The formation of a protein-protein complex is responsible for many biological functions; therefore, three-dimensional structures of protein complexes are essential for deeper understandings of protein functions and the mechanisms of diseases at the atomic level. However, compared with individual proteins, complex structures are difficult to solve experimentally because of technical limitations. Thus a method that can predict protein complex structures would be invaluable. In this study, we developed new software, FAMS Complex; a fully automated homology modeling system for protein complex structures consisting of two or more molecules. FAMS Complex requires only sequences and alignments of the target protein as input and constructs all molecules simultaneously and automatically. FAMS Complex is likely to become an essential tool for structure-based drug design, such as in silico screening to accelerate drug discovery before an experimental structure is solved. Moreover, in this post-genomic era when huge amounts of protein sequence information are available, a major goal is the determination of protein-protein interaction networks on a genomic scale. FAMS Complex will contribute to this goal, because its procedure is fully automated and so is suited for large-scale genome wide modeling.

In CASP6, the CHIMERA-group predicted full-atom models of all targets using SKE-CHIMERA, a Web-user interface system for protein structure prediction that allows human intervention at necessary stages; we used a lot of information from our own data and from publicly available data. Using SKE-CHIMERA, we iterated manual step (template selection and alignment by the in-house program CHIMERA) and automatic step (three-dimensional model building by the in-house program FAMS). The official CASP6 assessment showed that CHIMERA-group was one of the most successful predictors in homology modeling, especially for FR/H (Fold Recognition/Homologous). In this article, we introduce the method of CHIMERA-group and discuss its successes and failures in CASP6.

We report a novel method, ChooseLD (CHOOse biological information Semi-Empirically on the Ligand Docking), which uses simulated annealing (SA) based on bioinformatics for protein-ligand flexible docking. The fingerprint alignment score (FPAScore) value is used to determine the docking conformation of the ligand. This method includes the matching of chemical descriptors such as fingerprints (FPs) and the root mean square deviation (rmsd) calculation of the coordinates of atoms of the chemical descriptors. Here, the FPAScore optimization for the translation and rotation of a rigid body is performed using the Metropolis Monte Carlo method. Our ChooseLD method will find wide application in the field of biochemistry and medicine to improve the search for new drugs targeting various proteins implicated in diseases.

International audience; The CAPRI (Critical Assessment of Predicted Interactions) and CASP (Critical Assessment of protein Structure Prediction) experiments have demonstrated the power of community-wide tests of methodology in assessing the current state of the art and spurring progress in the very challenging areas of protein docking and structure prediction. We sought to bring the power of community-wide experiments to bear on a very challenging protein design problem that provides a complementary but equally fundamental test of current understanding of protein-binding thermodynamics. We have generated a number of designed protein-protein interfaces with very favorable computed binding energies but which do not appear to be formed in experiments, suggesting that there may be important physical chemistry missing in the energy calculations. A total of 28 research groups took up the challenge of determining what is missing: we provided structures of 87 designed complexes and 120 naturally occurring complexes and asked participants to identify energetic contributions and/or structural features that distinguish between the two sets. The community found that electrostatics and solvation terms partially distinguish the designs from the natural complexes, largely due to the nonpolar character of the designed interactions. Beyond this polarity difference, the community found that the designed binding surfaces were, on average, structurally less embedded in the designed monomers, suggesting that backbone conformational rigidity at the designed surface is important for realization of the designed function. These results can be used to improve computational design strategies, but there is still much to be learned; for example, one designed complex, which does form in experiments, was classified by all metrics as a nonbinder.

The prediction of a protein three-dimensional (3D) structure is one of the most important challenges in computational structural biology. We have developed an automatic protein 3D structure prediction method called FAMSD. FAMSD is based on a comparative modeling method which consists of the following four steps: (1) generating and selecting sequence alignments between target and template proteins; (2) constructing 3D structure models based on each selected alignment; (3) selecting the best 3D structure model and (4) refining the selected model. In the FAMSD method, sequence alignment programs such as a series of BLAST programs, SP3 and SPARKS2 programs, the homology modeling program FAMS (Full Automatic Modeling System), the model quality estimation program CIRCLE and the molecular dynamics program APRICOT were used in combination to construct high quality protein models. To assess the FAMSD method we have participated in the 8th Critical Assessment of Techniques for Protein Structure Prediction (CASP8) experiment. The results of our original assessment indicate that the FAMSD method offers excellent capability in packing side-chains with the correct torsion angles while avoiding the formation of atom-atom collisions. Since side-chain packing plays a significant role in defining the biological function of proteins, this method is a valuable resource in biological, pharmaceutical and medicinal research efforts.

Many disease‐protein targets have been found from biochemical experiments. Druggable compounds which inhibit or activate those significant protein targets must be researched rapidly. Many researches are using in‐sillico screening program such as DOCK, AutoDock and GOLD using classical mechanical potentials. We report the new method using a new operator which uses simulated annealing based on bioinformatics on the protein‐ligand flexible docking.

In CAPRI Rounds 1 and 2, we assumed that because there are many ionic charges that weaken electrostatic interaction forces in living cells, the hydrophobic interaction force might be important entropically. As a result of Rounds 1 and 2, the predictions for binding sites and geometric centers were acceptable, but those of the binding axes were poor, because only the largest benzene cluster was used for generating the initial docking structures. These were generated by fitting of benzene clusters formed on the surface of receptor and ligand. In CAPRI Rounds 3-5, the grid-scoring sum on the protein-protein interaction surface and the pairwise potential of the amino acid residues, which were indicated as coming easily into the protein-protein interaction regions, were used as the calculation methods, along with the smaller benzene clusters that participated in benzene cluster fitting. Good predicted models were obtained for Targets 11 and 12. When the modeled receptor proteins were superimposed on the experimental structures, the smallest ligand root-mean-square deviation (RMSD) values corresponding to the RMSD between the model and experimental structures were 6.2 A and 7.3 A, respectively.