Beyond the rotamer library: genetic algorithm combined with the disturbing mutation process for upbuilding protein side-chains

The disturbing genetic algorithm, incorporating the disturbing mutation process into the genetic algorithm flow, has been developed to extend the searching space of side-chain conformations and to improve the quality of the rotamer library. Moreover, the growing generation amount idea, simulating the real situation of the natural evolution, is introduced to improve the searching speed. In the calculations using the pseudo energy scoring function of the root mean squared deviation, the disturbing genetic algorithm method has been shown to be highly efficient. With the real energy function based on AMBER force field, the program has been applied to rebuilding side-chain conformations of 25 high-quality crystallographic structures of single-protein and protein-protein complexes. The averaged root mean standard deviation of atom coordinates in side-chains and veracities of the torsion angles of chi(1) and chi(1) + chi(2) are 1.165 A, 88.2 and 72.9% for the buried residues, respectively, and 1.493 A, 79.2 and 64.7% for all residues, showing that the method has equal precision to the program SCWRL, whereas it performs better in the prediction of buried residues and protein-protein interfaces. This method has been successfully used in redesigning the interface of the Basnase-Barstar complex, indicating that it will have extensive application in protein design, protein sequence and structure relationship studies, and research on protein-protein interaction.     

Copper complex of hydroxyl-substituted triazamacrocyclic ligand and its antitumor activity

The protonation constants and the stability constants for the formation of copper (II) complex of the ligand [1,4,7] Triazecan-9-ol (L) were presented. Antitumor activity of CuL complex was reported. Preliminary pharmacological tests showed that it had antitumor activity against HXO-RB44 and BEL-7402 cell lines in vitro. Nuclei of [CuL]-stimulated BEL-7402 cells clearly exhibited condensation and break down into chromatin clumps typical of apoptosis. Also it exhibited perturbation effects to BEL-7402 cell lines cycle and further studies showed that it could cleave supercoiled DNA (pBR 322) to nicked and linear DNA.     

Crystal structures of substrate complexes of malic enzyme and insights into the catalytic mechanism

Malic enzymes catalyze the oxidative decarboxylation of L-malate to pyruvate and CO(2) with the reduction of the NAD(P)(+) cofactor in the presence of divalent cations. We report the crystal structures at up to 2.1 A resolution of human mitochondrial NAD(P)(+)-dependent malic enzyme in different pentary complexes with the natural substrate malate or pyruvate, the dinucleotide cofactor NAD(+) or NADH, the divalent cation Mn(2+), and the allosteric activator fumarate. Malate is bound deep in the active site, providing two ligands for the cation, and its C4 carboxylate group is out of plane with the C1-C2-C3 atoms, facilitating decarboxylation. The divalent cation is positioned optimally to catalyze the entire reaction. Lys183 is the general base for the oxidation step, extracting the proton from the C2 hydroxyl of malate. Tyr112-Lys183 functions as the general acid-base pair to catalyze the tautomerization of the enolpyruvate product from decarboxylation to pyruvate.     

Higher-order interhelical spatial interactions in membrane proteins

Higher-order interactions are important for protein folding and assembly. We introduce the concept of interhelical three-body interactions as derived from Delaunay triangulation and alpha shapes of protein structures. In addition to glycophorin A, where triplets are strongly correlated with protein stability, we found that tight interhelical triplet interactions exist extensively in other membrane proteins, where many types of triplets occur far more frequently than in soluble proteins. We developed a probabilistic model for estimating the value of membrane helical interaction triplet (MHIT) propensity. Because the number of known structures of membrane proteins is limited, we developed a bootstrap method for determining the 95% confidence intervals of estimated MHIT values. We identified triplets that have high propensity for interhelical interactions and are unique to membrane proteins, e.g. AGF, AGG, GLL, GFF and others. A significant fraction (32%) of triplet types contains triplets that may be involved in interhelical hydrogen bond interactions, suggesting the prevalent and important roles of H-bond in the assembly of TM helices. There are several well-defined spatial conformations for triplet interactions on helices with similar parallel or antiparallel orientations and with similar right-handed or left-handed crossing angles. Often, they contain small residues and correspond to the regions of the closest contact between helices. Sequence motifs such as GG4 and AG4 can be part of the three-body interactions that have similar conformations, which in turn can be part of a higher-order cooperative four residue spatial motif observed in helical pairs from different proteins. In many cases, spatial motifs such as serine zipper and polar clamp are part of triplet interactions. On the basis of the analysis of the archaeal rhodopsin family of proteins, tightly packed triplet interactions can be achieved with several different choices of amino acid residues.     

Inferring functional relationships of proteins from local sequence and spatial surface patterns

We describe a novel approach for inferring functional relationship of proteins by detecting sequence and spatial patterns of protein surfaces. Well-formed concave surface regions in the form of pockets and voids are examined to identify similarity relationship that might be directly related to protein function. We first exhaustively identify and measure analytically all 910,379 surface pockets and interior voids on 12,177 protein structures from the Protein Data Bank. The similarity of patterns of residues forming pockets and voids are then assessed in sequence, in spatial arrangement, and in orientational arrangement. Statistical significance in the form of E and p-values is then estimated for each of the three types of similarity measurements. Our method is fully automated without human intervention and can be used without input of query patterns. It does not assume any prior knowledge of functional residues of a protein, and can detect similarity based on surface patterns small and large. It also tolerates, to some extent, conformational flexibility of functional sites. We show with examples that this method can detect functional relationship with specificity for members of the same protein family and superfamily, as well as remotely related functional surfaces from proteins of different fold structures. We envision that this method can be used for discovering novel functional relationship of protein surfaces, for functional annotation of protein structures with unknown biological roles, and for further inquiries on evolutionary origins of structural elements important for protein function.     

OsSET1, a novel SET-domain-containing gene from rice

A novel SET-domain-containing gene OsSET1 was isolated from rice (Oryza sativa L.). Its deduced protein consists of 895 amino acids. OsSET1 has a high degree of structure similarity to other SET-domain-containing genes such as CLF in higher plants and E(z) in animals. RT-PCR showed that the gene expresses throughout the entire plant. A transient expression assay in onion epidermis revealed that the OsSET1 protein is localized in nuclei. Over-expression of the SET domain of OsSET1 in Arabidopsis resulted in altered shoot development at seedling stages.     

Structure function differences in nonpeptide CCR1 antagonists for human and mouse CCR1

A useful strategy for identifying ligand binding domains of G protein-coupled receptors has been the exploitation of species differences in antagonist potencies. We have used this approach for the CCR1 chemokine receptor with a novel series of antagonists, the 4-hydroxypiperidines, which were discovered by high throughput screening of human CCR1 and subsequently optimized. The structure-activity relationships for a number of different 4-hydroxypiperidine antagonists for human and mouse CCR1 were examined by receptor binding and functional assays. These compounds exhibit major differences in their rank order of potency for the human and mouse chemokine receptor CCR1. For example, the initial lead template, BX 510, which was a highly potent functional antagonist for human CCR1 (K(i) = 21 nM) was >400-fold less active on mouse CCR1 (K(i) = 9150 nM). However, increasing the length of the linker between the piperidine and dibenzothiepine groups by one methylene group generated a compound, BX 511, which was equipotent for both human and mouse CCR1. These and other analogs of the lead template BX 510, which have major differences in potency for human and mouse CCR1, are described, and a model for their interaction with human CCR1 is presented.     

Assessment of requirements for IL-15 and IFN regulatory factors in uterine NK cell differentiation and function during pregnancy

In mouse and human, precursors of NK cell lineage home to decidualizing uteri. To assess the requirement for IL-15, an essential cytokine for NK differentiation in lymphoid tissue, on uterine NK (uNK) cell differentiation, implantation sites from IL-15(-/-) mice were analyzed histologically. IL-15(-/-) implantation sites had no uNK cells, no spiral-artery modification, and lacked the decidual integrity found in normal mice. IL-15(-/-) recipients of C57BL/6 marrow displayed similar pathology. However, implantation sites from recombination-activating gene-2(-/-)gamma(c)(-/-) (alymphoid) recipients of IL-15(-/-) marrow showed normal uNK cells, modified spiral arteries, and well-developed decidua basalis. Deletion of the IFN-regulatory factor (IRF)-1, but not IRF-2 (factors important in peripheral NK cell differentiation) limited but did not prevent uNK cell development. In situ hybridization localized IRF-1 largely to placental trophoblast cells. IRF-1(-/-) marrow transplanted into recombination-activating gene-2(-/-)gamma(c)(-/-) displayed competence for full uNK cell differentiation. IL-15 mRNA expression at implantation sites of IRF-1(-/-) and C57BL/6 was similar, suggesting that, unlike in bone marrow and spleen, IRF-1 does not regulate IL-15 in the pregnant uterus. Terminal differentiation of uNK cells was not promoted in pregnant IRF-1(-/-) mice by 5-day infusion of murine rIL-15, suggesting that IRF-1 deficiency rather than IL-15 deficiency limits uNK cell differentiation in these mice. Further, IRF-1 regulates placental growth, birth weight, and postnatal growth of offspring. These studies indicate that uNK cell development and maturation share some aspects with NK cell development in other tissues, but also display distinctive tissue-specific regulation.