Spatial chemical distance based on atomic property fields.

Similarity of compound chemical structures often leads to close pharmacological profiles, including binding to the same protein targets. The opposite, however, is not always true, as distinct chemical scaffolds can exhibit similar pharmacology as well. Therefore, relying on chemical similarity to known binders in search for novel chemicals targeting the same protein artificially narrows down the results and makes lead hopping impossible. In this study we attempt to design a compound similarity/distance measure that better captures structural aspects of their pharmacology and molecular interactions. The measure is based on our recently published method for compound spatial alignment with atomic property fields as a generalized 3D pharmacophoric potential. We optimized contributions of different atomic properties for better discrimination of compound pairs with the same pharmacology from those with different pharmacology using Partial Least Squares regression. Our proposed similarity measure was then tested for its ability to discriminate pharmacologically similar pairs from decoys on a large diverse dataset of 115 protein-ligand complexes. Compared to 2D Tanimoto and Shape Tanimoto approaches, our new approach led to improvement in the area under the receiver operating characteristic curve values in 66 and 58% of domains respectively. The improvement was particularly high for the previously problematic cases (weak performance of the 2D Tanimoto and Shape Tanimoto measures) with original AUC values below 0.8. In fact for these cases we obtained improvement in 86% of domains compare to 2D Tanimoto measure and 85% compare to Shape Tanimoto measure. The proposed spatial chemical distance measure can be used in virtual ligand screening.

Large-scale prediction of protein geometry and stability changes for arbitrary single point mutations.

We have developed a method to both predict the geometry and the relative stability of point mutants that may be used for arbitrary mutations. The geometry optimization procedure was first tested on a new benchmark of 2141 ordered pairs of X-ray crystal structures of proteins that differ by a single point mutation, the largest data set to date. An empirical energy function, which includes terms representing the energy contributions of the folded and denatured proteins and uses the predicted mutant side chain conformation, was fit to a training set consisting of half of a diverse set of 1816 experimental stability values for single point mutations in 81 different proteins. The data included a substantial number of small to large residue mutations not considered by previous prediction studies. After removing 22 (approximately 2%) outliers, the stability calculation gave a standard deviation of 1.08 kcal/mol with a correlation coefficient of 0.82. The prediction method was then tested on the remaining half of the experimental data, giving a standard deviation of 1.10 kcal/mol and covariance of 0.66 for 97% of the test set. A regression fit of the energy function to a subset of 137 mutants, for which both native and mutant structures were available, gave a prediction error comparable to that for the complete training set with predicted side chain conformations. We found that about half of the variation is due to conformation-independent residue contributions. Finally, a fit to the experimental stability data using these residue parameters exclusively suggests guidelines for improving protein stability in the absence of detailed structure information.

Modeling, mutagenesis, and structural studies on the fully conserved phosphate-binding loop (loop 8) of triosephosphate isomerase: toward a new substrate specificity.

Loop 8 (residues 232-242) in triosephosphate isomerase (TIM) is a highly conserved loop that forms a tight binding pocket for the phosphate moiety of the substrate. Its sequence includes the fully conserved, solvent-exposed Leu238. The tight phosphate-binding pocket explains the high substrate specificity of TIM being limited to the in vivo substrates dihydroxyacetone-phosphate and D-glyceraldehyde-3-phosphate. Here we use the monomeric variant of trypanosomal TIM for exploring the structural consequences of shortening this loop. The mutagenesis, guided by extensive modeling calculations and followed up by crystallographic characterization, is aimed at widening the phosphate-binding pocket and, consequently, changing the substrate specificity. Two new variants were characterized. The crystal structures of these variants indicate that in monomeric forms of TIM, the Leu238 side-chain is nicely buried in a hydrophobic cluster. Monomeric forms of wild-type dimeric TIM are known to exist transiently as folding intermediates; our structural analysis suggests that in this monomeric form, Leu238 of loop 8 also adopts this completely buried conformation, which explains its full conservation across the evolution. The much wider phosphate-binding pocket of the new variant allows for the development of a new TIM variant with a different substrate specificity.

Identification of ligands for RNA targets via structure-based virtual screening: HIV-1 TAR.

Binding of the Tat protein to TAR RNA is necessary for viral replication of HIV-1. We screened the Available Chemicals Directory (ACD) to identify ligands to bind to a TAR RNA structure using a four-step docking procedure: rigid docking first, followed by three steps of flexible docking using a pseudobrownian Monte Carlo minimization in torsion angle space with progressively more detailed conformational sampling on a progressively smaller list of top-ranking compounds. To validate the procedure, we successfully docked ligands for five RNA complexes of known structure. For ranking ligands according to binding avidity, an empirical binding free energy function was developed which accounts, in particular, for solvation, isomerization free energy, and changes in conformational entropy. System-specific parameters for the function were derived on a training set of RNA/ligand complexes with known structure and affinity. To validate the free energy function, we screened the entire ACD for ligands for an RNA aptamer which binds L-arginine tightly. The native ligand ranked 17 out of ca. 153,000 compounds screened, i.e., the procedure is able to filter out >99.98% of the database and still retain the native ligand. Screening of the ACD for TAR ligands yielded a high rank for all known TAR ligands contained in the ACD and suggested several other potential TAR ligands. Eight of the highest ranking compounds not previously known to be ligands were assayed for inhibition of the Tat-TAR interaction, and two exhibited a CD50 of ca. 1 microM.

TNF-alpha convertase enzyme from human arthritis-affected cartilage: isolation of cDNA by differential display, expression of the active enzyme, and regulation of TNF-alpha.

A snake venom-like protease isolated by a differential display screen between normal and osteoarthritis (OA)-affected cartilage (designated as cSVP) has a cDNA sequence identical to TNF-alpha convertase enzyme (TACE). TACE shows the presence of an unknown prodomain, a cysteine switch, a catalytic domain, a zinc binding region, a disintegrin region, an EGF-like domain, a transmembrane domain, and a unique cytoplasmic region. A TACE construct harboring the signal + prodomain + catalytic region (TACE-SPCdeltaDETCy), expressed in baculovirus could cleave preferentially (approximately 12-fold) the TNF-specific peptide over the matrix metalloproteases peptide in vitro. This recombinant protein also cleaved the natural substrate GST-ProTNF-alpha to TNF-alpha (17 kDa) in vitro. The mRNA for TACE, which is broadly distributed and differentially expressed in a variety of human tissues, is up-regulated in arthritis-affected cartilage, but not normal cartilage. OA-affected cartilage also expressed TNF-alpha mRNA that was not detected in normal cartilage. The OA-affected cartilage (in explant assays) spontaneously released TNF-alpha and IL-8 in ex vivo conditions. Addition of TNF-alphaR fused to IgG Fc fragment (TNF-alphaR:Fc) in the presence or absence of soluble IL-1R (with which it acted additively) significantly attenuated the spontaneous/autocrine release of articular IL-8 in this assay. These experiments demonstrate a functional paracrine/autocrine role of TNF-alpha in OA-affected cartilage that may depend, in part, on up-regulated levels of chondrocyte-derived TACE.

The 1.8 A crystal structure of the dimeric peroxisomal 3-ketoacyl-CoA thiolase of Saccharomyces cerevisiae: implications for substrate binding and reaction mechanism.

The dimeric, peroxisomal 3-ketoacyl-CoA thiolase catalyses the conversion of 3-ketoacyl-CoA into acyl-CoA, which is shorter by two carbon atoms. This reaction is the last step of the beta-oxidation pathway. The crystal structure of unliganded peroxisomal thiolase of the yeast Saccharomyces cerevisiae has been refined at 1.8 A resolution. An unusual feature of this structure is the presence of two helices, completely buried in the dimer and sandwiched between two beta-sheets. The analysis of the structure shows that the sequences of these helices are not hydrophobic, but generate two amphipathic helices. The helix in the N-terminal domain exposes the polar side-chains to a cavity at the dimer interface, filled with structured water molecules. The central helix in the C-terminal domain exposes its polar residues to an interior polar pocket. The refined structure has also been used to predict the mode of binding of the substrate molecule acetoacetyl-CoA, as well as the reaction mechanism. From previous studies it is known that Cys125, His375 and Cys403 are important catalytic residues. In the proposed model the acetoacetyl group fits near the two catalytic cysteine residues, such that the oxygen atoms point towards the protein interior. The distance between SG(Cys125) and C3(acetoacetyl-CoA) is 3.7 A. The O2 atom of the docked acetoacetyl group makes a hydrogen bond to N(Gly405), which would favour the formation of the covalent bond between SG(Cys125) and C3(acetoacetyl-CoA) of the intermediate complex of the two-step reaction. The CoA moiety is proposed to bind in a groove on the surface of the protein molecule. Most of the interactions of the CoA molecule are with atoms of the loop domain. The three phosphate groups of the CoA moiety are predicted to interact with side-chains of lysine and arginine residues, which are conserved in the dimeric thiolases.

Do aligned sequences share the same fold?

Sequence comparison remains a powerful tool to assess the structural relatedness of two proteins. To develop a sensitive sequence-based procedure for fold recognition, we performed an exhaustive global alignment (with zero end gap penalties) between sequences of protein domains with known three-dimensional folds. The subset of 1.3 million alignments between sequences of structurally unrelated domains was used to derive a set of analytical functions that represent the probability of structural significance for any sequence alignment at a given sequence identity, sequence similarity and alignment score. Analysis of overlap between structurally significant and insignificant alignments shows that sequence identity and sequence similarity measures are poor indicators of structural relatedness in the "twilight zone", while the alignment score allows much better discrimination between alignments of structurally related and unrelated sequences for a wide variety of alignment settings. A fold recognition benchmark was used to compare eight different substitution matrices with eight sets of gap penalties. The best performing matrices were Gonnet and Blosum50 with normalized gap penalties of 2.4/0.15 and 2.0/0.15, respectively, while the positive matrices were the worst performers. The derived functions and parameters can be used for fold recognition via a multilink chain of probability weighted pairwise sequence alignments.

Contact area difference (CAD): a robust measure to evaluate accuracy of protein models.

A simple unified measure to evaluate the accuracy of three-dimensional atomic protein models is proposed. This measure is a normalized sum of absolute differences of residue-residue contact surface areas calculated for a reference structure and a model. It employs more rigorous quantitative evaluation of a contact than previously used contact measures. We argue that the contact area difference (CAD) number is a robust single measure to evaluate protein structure predictions in a wide range of model accuracies, from ab initio and threading models to models by homology, since it reflects both backbone topology and side-chain packing, is smooth, continuous and threshold-free, is not sensitive to typical crystallographic errors and ambiguities, adequately penalizes domain and/or secondary structure rearrangements and protein plasticity, and has consistent linear and matrix representations for more detailed analysis. The CAD quality of crystallographic structures, NMR structures, models by homology, and unfolded and misfolded structures is evaluated. It is shown that the CAD number discriminates between models better than Cartesian root-mean-square deviation (cRMSD). Structural variability of the NMR structures was found to be three times larger than deformations of crystallographic structures in different packing environments.

Protein engineering with monomeric triosephosphate isomerase (monoTIM): the modelling and structure verification of a seven-residue loop.

Protein engineering experiments have been carried out with loop-1 of monomeric triosephosphate isomerase (monoTIM). Loop-1 of monoTIM is disordered in every crystal structure of liganded monoTIM, but in the wild-type TIM it is a very rigid dimer interface loop. This loop connects the first beta-strand with the first alpha-helix of the TIM-barrel scaffold. The first residue of this loop, Lys13, is a conserved catalytic residue. The protein design studies with loop-1 were aimed at rigidifying this loop such that the Lys13 side chain points in the same direction as seen in wild type. The modelling suggested that the loop should be made one residue shorter. With the modelling package ICM the optimal sequence of a new seven-residue loop-1 was determined and its structure was predicted. The new variant could be expressed and purified and has been characterized. The catalytic activity and stability are very similar to those of monoTIM. The crystal structure (at 2.6 A resolution) shows that the experimental loop-1 structure agrees well with the modelled loop-1 structure. The direct superposition of the seven loop residues of the modelled and experimental structures results in an r.m.s. difference of 0.5 A for the 28 main chain atoms. The good agreement between the predicted structure and the crystal structure shows that the described modelling protocol can be used successfully for the reliable prediction of loop structures.

A technique for identifying atoms from a screen image.

Improving the interfaces in molecular graphics applications, making them more natural and easy to use, is an important task, given the current complexity of the displayed objects and of modeling operations. Clicking near an atom center is the usual method of atom selection. However, this method has certain disadvantages when working with images composed of different atomic representations such as sticks, CPK, or dotted surfaces. We propose another technique allowing the user to obtain the correct answer when he or she clicks on any element of the atom image.

Towards protein folding by global energy optimization.

Different components of the theoretical protein folding problem are evaluated critically. It is argued that: (i) as a rule, small- and medium-sized proteins are in the free energy minimum; (ii) long-living metastable states may either appear occasionally with growing protein size, or be selected by evolution for a specific function; (iii) functions discriminating against incorrect folds would fail if they were used directly in the global optimization, unless they approximate the true free energy accurately; (iv) surface and electrostatic free energies should be treated separately; (v) conformational entropy (of side chains in particular) should be taken into account; (vi) Monte Carlo procedures considering all free energy terms and combining global knowledge-based random moves with local optimization have the largest potential for success.

Proposed structure for the DNA-binding domain of the helix-loop-helix family of eukaryotic gene regulatory proteins.

A modelled tertiary structure for the dimeric HLH domain of the E47 protein is presented. Structural information was obtained from the aligned sequences of > 40 members of the HLH family. The information was used to model each monomer as an alpha-helical hairpin, with knobs-into-holes packing of side-chains as found in antiparallel coiled-coil. The dimer forms a four-helix bundle with additional knobs-into-holes packing at the dimer interface. The size and electrostatic properties of core-forming residues are all accounted for in the model. The model does not violate any known properties of protein structure. The monomers are related by two-fold rotational symmetry, in agreement with the observed DNA-binding sites which are imperfect inverted repeats. The N-terminal basic region, in which DNA binding and base specificity reside, forms the first part of helix 1. A prediction based on the model structure is that the HLH domains do not bind to DNA in its B form but require a partially unwound conformation in order to enter the major groove.

Electrophoretic behavior of d(GGAAAAAAGG)n, d(CCAAAAAACC)n, and (CCAAAAAAGG)n and implications for a DNA bending model.

Double stranded multimers (C2A6C2)n, (C2A6G2)n and (G2A6G2)n were prepared from chemically synthesized oligonucleotides to study the influence of sequences flanking the An tract on the curvature of DNA. All these duplexes, including polypurine.polypyrimidine one, exhibit strong retardation in polyacrylamide gel which is indicative of pronounced DNA curvature. It has been proposed previously that among the bends at the boundary with the oligo(A) tract two types should be distinguished: 5'-bends and 3'-bends (Koo et al., 1986) This distinction was deduced from different relative mobilities of two specially designed sequences having phased 5'-bends and 3'-bends, respectively. Our data indicate that the substitutions of nucleotides at both 5' and 3' boundaries of A6 tract result in comparable changes in relative mobility. Therefore, for B-B' bends it is important to take into account not only whether they are at the 5' or 3' end of an oligo(dA) tract, but also the particular sequences at the boundaries of this tract.

Sequence dependent modulating effect of camptothecin on the DNA-cleaving activity of the calf thymus type I topoisomerase.

High-resolution mapping of topol cleavages in the regions of human DNA including the oncogene c-Ha-ras and p53, has revealed three kinds of topol cleavage sites: cleavage sites not affected by camptothecin; cleavage sites reinforced only in the presence of camptothecin, and cleavage sites which weaken in the presence of camptothecin. Statistical analysis of sequences revealed certain nucleotide or dinucleotide preferences for three groups studied. The preferences in camptothecin-reduced sites predominate upstream from the cleavage point, whereas in camptothecin-induced sites the situation is reversed. The influence of camptothecin on cleavage sites induced by two molecular forms of topol has been also studied.

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.

An automatic search for similar spatial arrangements of alpha-helices and beta-strands in globular proteins.

A fast search algorithm to reveal similar polypeptide backbone structural motifs in proteins is proposed. It is based on the vector representation of a polypeptide chain fold in which the elements of regular secondary structures are approximated by linear segments (Abagyan and Maiorov, J. Biomol. Struct. Dyn. 5, 1267-1279 (1988)). The algorithm permits insertions and deletions in the polypeptide chain fragments to be compared. The fast search algorithm implemented in FASEAR program is used for collecting beta alpha beta supersecondary structure units in a number of alpha/beta proteins of Brookhaven Data Bank. Variation of geometrical parameters specifying backbone chain fold is estimated. It appears that the conformation of the majority of the fragments, although almost all of them are right-handed, is quite different from that of standard beta alpha beta units. Apart from searching for specific type of secondary structure motif, the algorithm allows automatically to identify new recurrent folding patterns in proteins. It may be of particular interest for the development of tertiary template approach for prediction of protein three-dimensional structure as well for constructing artificial polypeptides with goal-oriented conformation.

New methodology for computer-aided modelling of biomolecular structure and dynamics. 1. Non-cyclic structures.

A general methodology is proposed for the conformational modelling of biomolecular systems. The approach allows one: (i) to describe the system under investigation by an arbitrary set of internal variables, i.e., torsion angles, bond angles, and bond lengths; it offers a possibility to pass from the free structure to a completely fixed one with the number of variables from 3N to zero, respectively, where N is the number of atoms; (ii) to consider both, a single molecule and a complex of many molecules, (e.g., proteins, water, ligands, etc.) in terms of one universal model; (iii) to study the dynamics of the system using explicit analytical Lagrangian equations of motion, thus opening up possibilities for investigations of slow concerted motions such as domain oscillations in proteins etc.; (iv) to calculate the partial derivatives of various functions of conformation, e.g., the conformational energy or external constraints imposed, using a standard efficient procedure regardless of the variables and the structure of the system. The approach is meant to be used in various investigations concerning the conformations and dynamics of biomacromolecules.

New methodology for computer-aided modelling of biomolecular structure and dynamics. 2. Local deformations and cycles.

A new methodology for the conformational modelling of biomolecular systems (1) is extended to local deformations of chain molecules and to flexible molecular rings. It is shown that these two cases may be reduced to considering an equivalent molecular model with a regular tree-like topology. A simple procedure is developed to analyze any flexible rings (the five- and six-membered sugar rings of carbohydrates and nucleic acids, in particular) and local deformation regions by energy minimization. Dynamic equations are also derived for such molecular systems. As a result, a unified approach is proposed for the efficient energy minimization and simulation of dynamic behavior of multimolecular systems having any set of variable internal coordinates, local deformation regions and cycles. Advantages and domains of applicability of the approach are discussed.

Structural basis of stable bending in DNA containing An tracts. Different types of bending.

Structural determinants of DNA bending of different types have been studied by theoretical conformational analysis of duplexes. Their terminal parts were fixed either in an ordinary low-energy B-like conformation or in "anomalous" conformations with a narrowed minor groove typical of An tracts. The anomalous conformations had different negative tilt angles (up to about zero), different propeller twists and minor groove widths. Calculations have been performed for DNA fragments AnTm, TnAm, AnGCTm, AnCGTm, TmGCAn, TmCGAn which are the models of the junction of two anomalous structures on An and Tm tracts. On the AT step of the AnTm fragment the minor groove can be easily narrowed so that a whole unbent fragment of anomalous structure is formed on AnTm. According to our energy estimates, there should not be any reliable bending on AnTm. In contrast, in all other cases there was a pronounced roll-like bending into the major groove in the chemical symmetry region. Calculations of the junction between the anomalous and ordinary B-like structure for GnTm and CnAm have shown that there is an equilibrium bending with a tilt component towards the chain having the anomalous structure at the 5'-end. From our calculations it is impossible to determine precisely the direction of bending, though it can be suggested that the roll component of bending might be directed towards the major groove. The anomalous structure is the main reason of bending; alternations of pyrimidines and purines can modulate the value and the direction of equilibrium bending (only the value in the case of self-complementary fragments).(ABSTRACT TRUNCATED AT 250 WORDS)