Homology modeling by the ICM method

Five models have been built by the ICM method for the Comparative Modeling section of the Meeting on the Critical Assessment of Techniques for Protein Structure Prediction. The targets have homologous proteins with known three-dimensional structure with sequence identity ranging from 25 to 77%. After alignment of the target sequence with the related three-dimensional structure, the modeling procedure consists of two subproblems: side-chain prediction and loop prediction. The ICM method approaches these problems with the following steps: (1) a starting model is created based on the homologous structure with the conserved portion fixed and the noncon-served portion having standard covalent geometry and free torsion angles; (2) the Biased Probability Monte Carlo (BPMC) procedure is applied to search the subspaces of either all the nonconservative side-chain torsion angles or torsion angles in a loop backbone and surrounding side chains. A special algorithm was designed to generate low-energy loop deformations. The BPMC procedure globally optimizes the energy function consisting of ECEPP/3 and solvation energy terms. Comparison of the predictions with the NMR or crystallographic solutions reveals a high proportion of correctly predicted side chains. The loops were not correctly predicted because imprinted distortions of the backbone increased the energy of the near-native conformation and thus made the solution unrecognizable. Interestingly, the energy terms were found to be reliable and the sampling of conformational space sufficient. The implications of this finding for the strategies of future comparative modeling are discussed. © 1995 Wiley-Liss, Inc.     

Structure of the hydration shells of oligo(dA-dT).oligo(dA-dT) and oligo(dA).oligo(dT) tracts in B-type conformation on the basis of Monte Carlo calculations

Monte Carlo simulations [(N, V, T)-ensemble] were performed for the hydration shell of poly(dA-dT).poly(dA-dT) in canonical B form and for the hydration shell of poly(dA).poly(dT) in canonical B conformation and in a conformation with narrow minor groove, highly inclined bases, but with a nearly zero-inclined base pair plane (B' conformation). We introduced helical periodic boundary conditions with a rather small unit cell and a limited number of water molecules to reduce the dimensionality of the configuration space. The coordinates of local maxima of water density and the properties of one- and two-membered water bridges between polar groups of the DNA were obtained. The AT-alternating duplex hydration mirrors the dyad symmetry of polar group distribution. At the dApdT step, a water bridge between the two carbonyl oxygens O2 of thymines is formed as in the central base-pair step of Dickerson's dodecamer. In the major groove, 5-membered water chains along the tetranucleotide pattern d(TATA).d(TATA) are observed. The hydration geometry of poly(dA).poly(dT) in canonical B conformation is distinguished by autonomous primary hydration of the base-pair edges in both grooves. When this polymer adopts a conformation with highly inclined bases and narrow minor groove, the water density distribution in the minor groove is in excellent agreement with Dickerson's spine model. One local maximum per base pair of the first layer is located near the dyad axis between adjacent base pairs, and one local maximum per base pair in the second shell lies near the dyad axis of the base pair itself. The water bridge between the two strands formed within the first layer was observed with high probability. But the water molecules of the second layer do not have a statistically favored orientation necessary for bridging first layer waters. In the major groove, the hydration geometry of the (A.T) base-pair edge resembles the main features of the AT-pair hydration derived from other sequences for the canonical B form. The preference of the B' conformation for oligo(dA).oligo(dT) tracts may express the tendency to common hydration of base-pair edges of successive base pairs in the grooves of B-type DNA. The mean potential energy of hydration of canonical B-DNA was estimated to be -60 to -80 kJ/mole nucleotides in dependence on the (G.C) contents. Because of the small system size, this estimation is preliminary.     

ADP_EM: fast exhaustive multi-resolution docking for high-throughput coverage

MOTIVATION: Efficient fitting tools are needed to take advantage of a fast growth of atomic models of protein domains from crystallography or comparative modeling, and low-resolution density maps of larger molecular assemblies. Here, we report a novel fitting algorithm for the exhaustive and fast overlay of partial high-resolution models into a low-resolution density map. The method incorporates a fast rotational search based on spherical harmonics (SH) combined with a simple translational scanning. RESULTS: This novel combination makes it possible to accurately dock atomic structures into low-resolution electron-density maps in times ranging from seconds to a few minutes. The high-efficiency achieved with simulated and experimental test cases preserves the exhaustiveness needed in these heterogeneous-resolution merging tools. The results demonstrate its efficiency, robustness and high-throughput coverage. AVAILABILITY: http://sbg.cib.csic.es/Software/ADP_EM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

REVCOM: a robust Bayesian method for evolutionary rate estimation

MOTIVATION: Evolutionary conservation estimated from a multiple sequence alignment is a powerful indicator of the functional significance of a residue and helps to predict active sites, ligand binding sites, and protein interaction interfaces. Many algorithms that calculate conservation work well, provided an accurate and balanced alignment is used. However, such a strong dependence on the alignment makes the results highly variable. We attempted to improve the conservation prediction algorithm by making it more robust and less sensitive to (1) local alignment errors, (2) overrepresentation of sequences in some branches and (3) occasional presence of unrelated sequences. RESULTS: A novel method is presented for robust constrained Bayesian estimation of evolutionary rates that avoids overfitting independent rates and satisfies the above requirements. The method is evaluated and compared with an entropy-based conservation measure on a set of 1494 protein interfaces. We demonstrated that approximately 62% of the analyzed protein interfaces are more conserved than the remaining surface at the 5% significance level. A consistent method to incorporate alignment reliability is proposed and demonstrated to reduce arbitrary variation of calculated rates upon inclusion of distantly related or unrelated sequences into the alignment.     

Ab initio prediction of peptide-MHC binding geometry for diverse class I MHC allotypes

Since determining the crystallographic structure of all peptide-MHC complexes is infeasible, an accurate prediction of the conformation is a critical computational problem. These models can be useful for determining binding energetics, predicting the structures of specific ternary complexes with T-cell receptors, and designing new molecules interacting with these complexes. The main difficulties are (1) adequate sampling of the large number of conformational degrees of freedom for the flexible peptide, (2) predicting subtle changes in the MHC interface geometry upon binding, and (3) building models for numerous MHC allotypes without known structures. Whereas previous studies have approached the sampling problem by dividing the conformational variables into different sets and predicting them separately, we have refined the Biased-Probability Monte Carlo docking protocol in internal coordinates to optimize a physical energy function for all peptide variables simultaneously. We also imitated the induced fit by docking into a more permissive smooth grid representation of the MHC followed by refinement and reranking using an all-atom MHC model. Our method was tested by a comparison of the results of cross-docking 14 peptides into HLA-A*0201 and 9 peptides into H-2K(b) as well as docking peptides into homology models for five different HLA allotypes with a comprehensive set of experimental structures. The surprisingly accurate prediction (0.75 A backbone RMSD) for cross-docking of a highly flexible decapeptide, dissimilar to the original bound peptide, as well as docking predictions using homology models for two allotypes with low average backbone RMSDs of less than 1.0 A illustrate the method's effectiveness. Finally, energy terms calculated using the predicted structures were combined with supervised learning on a large data set to classify peptides as either HLA-A*0201 binders or nonbinders. In contrast with sequence-based prediction methods, this model was also able to predict the binding affinity for peptides to a different MHC allotype (H-2K(b)), not used for training, with comparable prediction accuracy.     

Damped-dynamics flexible fitting

In fitting atomic structures into EM maps, it often happens that the map corresponds to a different conformation of the structure. We have developed a new methodology to handle these situations that preserves the covalent geometry of the structure and allows the modeling of large deformations. The first goal is achieved by working in generalized coordinates (positional and internal coordinates), and the second by avoiding harmonic potentials. Instead, we use dampers (shock absorbers) between every pair of atoms, combined with a force field that attracts the atomic structure toward incompletely occupied regions of the EM map. The trajectory obtained by integrating the resulting equations of motion converges to a conformation that, in our validation cases, was very close to the target atomic structure. Compared to current methods, our approach is more efficient and robust against wrong solutions and to overfitting, and does not require user intervention or subjective decisions. Applications to the computation of transition pathways between known conformers, homology and loop modeling, as well as protein docking, are also discussed.     

Flexible ligand docking to multiple receptor conformations: a practical alternative

State of the art docking algorithms predict an incorrect binding pose for about 50-70% of all ligands when only a single fixed receptor conformation is considered. In many more cases, lack of receptor flexibility results in meaningless ligand binding scores, even when the correct pose is obtained. Incorporating conformational rearrangements of the receptor binding pocket into predictions of both ligand binding pose and binding score is crucial for improving structure-based drug design and virtual ligand screening methodologies. However, direct modeling of protein binding site flexibility remains challenging because of the large conformational space that must be sampled, and difficulties remain in constructing a suitably accurate energy function. Here we show that using multiple fixed receptor conformations, either experimentally determined by crystallography or NMR, or computationally generated, is a practical shortcut that may improve docking calculations. In several cases, such an approach has led to experimentally validated predictions.     

The transthyretin-related protein family.

A number of proteins related to the homotetrameric transport protein transthyretin (TTR) forms a highly conserved protein family, which we present in an integrated analysis of data from different sources combined with an initial biochemical characterization. Homologues of the transthyretin-related protein (TRP) can be found in a wide range of species including bacteria, plants and animals, whereas transthyretins have so far only been identified in vertebrates. A multiple sequence alignment of 49 TRP sequences from 47 species to TTR suggests that the tertiary and quaternary features of the three-dimensional structure are most likely preserved. Interestingly, while some of the TRP orthologues show as little as 30% identity, the residues at the putative ligand-binding site are almost entirely conserved. RT/PCR analysis in Caenorhabditis elegans confirms that one TRP gene is transcribed, spliced and predominantly expressed in the worm, which suggests that at least one of the two C. elegans TRP genes encodes a functional protein. We used double-stranded RNA-mediated interference techniques in order to determine the loss-of-function phenotype for the two TRP genes in C. elegans but detected no apparent phenotype. The cloning and initial characterization of purified TRP from Escherichia coli reveals that, while still forming a homotetramer, this protein does not recognize thyroid hormones that are the natural ligands of TTR. The ligand for TRP is not known; however, genomic data support a functional role involving purine catabolism especially linked to urate oxidase (uricase) activity.     

Role of lysine187 within the second extracellular loop of the type A cholecystokinin receptor in agonist-induced activation. Use of complementary charge-reversal mutagenesis to define a functionally important interdomain interaction

Activation of guanine nucleotide-binding protein (G protein)-coupled receptors is believed to involve conformational change that exposes a domain for G protein coupling at the cytosolic surface of the helical confluence, although the mechanisms for achieving this are not well understood. This conformational change can be achieved by docking a diverse variety of agonist ligands, known to occur by interacting with different regions of these receptors. In this study, we focus on the importance of a specific basic residue (Lys187) within the second extracellular loop of the receptor for the peptide hormone, cholecystokinin. Alanine-replacement and charge-reversal mutagenesis of this residue showed that it had no effect on the binding of natural peptide and nonpeptidyl ligands of this receptor but markedly interfered with agonist-stimulated signaling. It was demonstrated that this negative effect on biological activity could be eliminated with the truncation of the first 30 residues of the amino-terminal tail of this receptor. Complementary charge-reversal mutagenesis of each of the five conserved acidic residues within this region of the receptor in the presence of the charge-reversed Lys187 revealed that only the Asp5 mutant fully reversed the negative functional impact of the Lys187 charge reversal. Thus, we have demonstrated that a basic residue within the second extracellular loop of the cholecystokinin receptor interacts with a specific acidic residue within the amino terminus of this receptor. This residue-residue interaction is nicely accommodated within a new molecular model of the agonist-occupied cholecystokinin receptor.     

A new method for modeling large-scale rearrangements of protein domains

A method for modeling large-scale rearrangements of protein domains connected by a single- or a double-stranded linker is proposed. Multidomain proteins may undergo substantial domain displacements, while their intradomain structure remains essentially unchanged. The method allows automatic identification of an interdomain linker and builds an all-atom model of a protein structure in internal coordinates. Torsion angles belonging to the interdomain linkers and side chains potentially able to form domain interfaces are set free while all remaining torsions, bond lengths, and bond angles are fixed. Large-scale sampling of the reduced torsion conformational subspace is effected with the “biased probability Monte Carlo-minimization” method [Abagyan, R.A., Totrov, M.M. (1994): J. Mol. Biol. 235, 983–1002]. Solvation and side-chain entropic contributions are added to the energy function. A special procedure has been developed to generate concerted deformations of a double-stranded interdomain linker in such a way that the polypeptide chain continuity is preserved. The method was tested on Bence-Jones protein with a single-stranded linker and lysine/arginine/ornithine-binding (LAO) protein with a double-stranded linker. For each protein, structurally diverse low-energy conformations with ideal covalent geometry were generated, and an overlap between two sets of conformations generated starting from the crystallographically determined “closed” and “open” forms was found. One of the low-energy conformations generated in a run starting from the LAO “closed” form was only 2.2 Å away from the structure of the “open” form. The method can be useful in predicting the scope of possible domain rearrangements of a multidomain protein. Proteins 27:410–424, 1997. © 1997 Wiley-Liss, Inc.