The Role of Bridging Water and Hydrogen Bonding as Key Determinants of Noncovalent Protein-Carbohydrate Recognition.

Mechanisms of protein-carbohydrate recognition attract a lot of interest due to their roles in various cellular processes and metabolism disorders. We have performed a large-scale analysis of protein structures solved in complex with glucose, galactose and their substituted analogues. We found that, on average, sugar molecules establish five hydrogen bonds (HBs) in the binding site, including one to three HBs with bridging water molecules. The free energy contribution of bridging and direct HBs was estimated using the free energy perturbation (FEP+) methodology for mono- and disaccharides that bind to l-ABP, ttGBP, TrmB, hGalectin-1 and hGalectin-3. We show that removing hydroxy groups that are engaged in direct HBs with the charged groups of Asp, Arg and Glu residues, protein backbone amide or buried water dramatically decreases binding affinity. In contrast, all solvent-exposed hydroxy groups and hydroxy groups engaged in HBs with the solvent-exposed bridging water molecules contribute weakly to binding affinity and so can be replaced to optimize ligand potency. Finally, we rationalize an effect of binding site water replacement on the binding affinity to l-ABP.

PDBFlex: exploring flexibility in protein structures.

The PDBFlex database, available freely and with no login requirements at http://pdbflex.org, provides information on flexibility of protein structures as revealed by the analysis of variations between depositions of different structural models of the same protein in the Protein Data Bank (PDB). PDBFlex collects information on all instances of such depositions, identifying them by a 95% sequence identity threshold, performs analysis of their structural differences and clusters them according to their structural similarities for easy analysis. The PDBFlex contains tools and viewers enabling in-depth examination of structural variability including: 2D-scaling visualization of RMSD distances between structures of the same protein, graphs of average local RMSD in the aligned structures of protein chains, graphical presentation of differences in secondary structure and observed structural disorder (unresolved residues), difference distance maps between all sets of coordinates and 3D views of individual structures and simulated transitions between different conformations, the latter displayed using JSMol visualization software.

MORPH-PRO: a novel algorithm and web server for protein morphing.

BACKGROUND: Proteins are known to be dynamic in nature, changing from one conformation to another while performing vital cellular tasks. It is important to understand these movements in order to better understand protein function. At the same time, experimental techniques provide us with only single snapshots of the whole ensemble of available conformations. Computational protein morphing provides a visualization of a protein structure transitioning from one conformation to another by producing a series of intermediate conformations. RESULTS: We present a novel, efficient morphing algorithm, Morph-Pro based on linear interpolation. We also show that apart from visualization, morphing can be used to provide plausible intermediate structures. We test this by using the intermediate structures of a c-Jun N-terminal kinase (JNK1) conformational change in a virtual docking experiment. The structures are shown to dock with higher score to known JNK1-binding ligands than structures solved using X-Ray crystallography. This experiment demonstrates the potential applications of the intermediate structures in modeling or virtual screening efforts. CONCLUSIONS: Visualization of protein conformational changes is important for characterization of protein function. Furthermore, the intermediate structures produced by our algorithm are good approximations to true structures. We believe there is great potential for these computationally predicted structures in protein-ligand docking experiments and virtual screening. The Morph-Pro web server can be accessed at http://morph-pro.bioinf.spbau.ru.

Structural and functional analysis of the globular head domain of p115 provides insight into membrane tethering.

Molecular tethers have a central role in the organization of the complex membrane architecture of eukaryotic cells. p115 is a ubiquitous, essential tether involved in vesicle transport and the structural organization of the exocytic pathway. We describe two crystal structures of the N-terminal domain of p115 at 2.0 A resolution. The p115 structures show a novel alpha-solenoid architecture constructed of 12 armadillo-like, tether-repeat, alpha-helical tripod motifs. We find that the H1 TR binds the Rab1 GTPase involved in endoplasmic reticulum to Golgi transport. Mutation of the H1 motif results in the dominant negative inhibition of endoplasmic reticulum to Golgi trafficking. We propose that the H1 helical tripod contributes to the assembly of Rab-dependent complexes responsible for the tether and SNARE-dependent fusion of membranes.

Structural changes of vitamin D receptor induced by 20-epi-1alpha,25-(OH)2D3: an insight from a computational analysis.

We employ a new computational tool CCOMP for the comparison of side chain (SC) conformations between crystal structures of homologous protein complexes. The program is applied to the vitamin D receptor (VDR) liganded with 1alpha,25-(OH)(2)D(3) (in 1DB1) or its 20-epi (in 1IE9) analog with an inverted C-20 configuration. This modification yields no detectable changes in the backbone configuration or ligand topology in the receptor binding cavity, yet it dramatically increases transcription, differentiation and antiproliferation activity of the VDR. We applied very stringent criteria during the comparison process. To eliminate errors arising from the different packing of investigated crystals and the thermal flexibility of atoms, we studied complexes belonging to the same space group, having a low R value (0.2) and a B-factor below 40 for compared residues. We find that 20-epi-1alpha,25-(OH)(2)D(3) changes side chain conformation of amino acids residing far away from direct ligand-VDR contacts. We further verify that a number of the reoriented residues were identified in mutational experiments as important for interaction with SRC-1, GRIP, TAFs co-activators and VDR-RXR heterodimerization. Thus, CCOMP analysis of protein complexes may be used for identifying amino acids that could serve as targets for genetic engineering, such as mutagenesis.

Fast procedure for reconstruction of full-atom protein models from reduced representations.

We introduce PULCHRA, a fast and robust method for the reconstruction of full-atom protein models starting from a reduced protein representation. The algorithm is particularly suitable as an intermediate step between coarse-grained model-based structure prediction and applications requiring an all-atom structure, such as molecular dynamics, protein-ligand docking, structure-based function prediction, or assessment of quality of the predicted structure. The accuracy of the method was tested on a set of high-resolution crystallographic structures as well as on a set of low-resolution protein decoys generated by a protein structure prediction algorithm TASSER. The method is implemented as a standalone program that is available for download from http://cssb.biology.gatech.edu/skolnick/files/PULCHRA.

Computational analysis of the active sites in binary and ternary complexes of the vitamin D receptor.

We have developed a program CCOMP that compares overlapping fragments of two protein complexes and identifies differently oriented amino acids. CCOMP initially performs a sequence alignment of the analyzed receptors, then superimposes the corresponding aligned residues, and finally calculates the root mean square deviation (RMSD) of individual atoms, every amino acid and the entire complex. Thus, amino acids important for functional differences between both complexes can be detected. Application of CCOMP to 1alpha,25-(OH)(2)D(3)-hVDR (1DB1) [Proc. Natl. Acad. Sci. U.S.A. 98 (2001) 5491] and 1alpha,25-(OH)(2)D(3)-rVDR-peptide (1RK3) [Biochemistry 43 (2004) 4101] complexes revealed that the peptide (KNHPMLMNLLKDN) mimicking a co-activator sequence significantly changes the side chain conformation of 35 amino acids. Four of these residues (K242, I256, K260, E416) actually contact the peptide, but all of them are essential for biological activity. Only two (L309 and L400) of the 35 differently oriented amino acids contact the ligand. Interestingly, when the peptide is present (1RK3) leucine 400 shifts closer (0.7A) to the vitamin D 26-methyl group. Applying the CCOMP and DSSP programs to binary and ternary VDR complexes also resulted in establishing that seven amino acids (I238, S252, I256, L413, L415, E416, V417) exhibit significant differences in solvent accessibility and are capable of interacting with co-activators.

DNA vaccine expressing the mimotope of GD2 ganglioside induces protective GD2 cross-reactive antibody responses.

The GD2 ganglioside expressed on neuroectodermally derived tumors, including neuroblastoma and melanoma, is weakly immunogenic in tumor-bearing patients and induces predominantly immunoglobulin (Ig)-M antibody responses in the immunized host. Here, we investigated whether interconversion of GD2 into a peptide mimetic form would induce GD2 cross-reactive IgG antibody responses in mice. Screening of the X(15) phage display peptide library with the anti-GD2 monoclonal antibody (mAb) 14G2a led to isolation of mimetic peptide 47, which inhibited the binding of 14G2a antibody to GD2-positive tumor cells. The peptide was also recognized by GD2-specific serum antibodies from a patient with neuroblastoma, suggesting that it bears an internal image of GD2 ganglioside expressed on the tumor cells. The molecular basis for antigenicity of the GD2 mimetic peptide, established by molecular modeling and mutagenesis studies, led to the generation of a 47-LDA mutant with an increased mimicry to GD2. Immunization of mice with peptide 47-LDA-encoded plasmid DNA elicited GD2 cross-reactive IgG antibody responses, which were increased on subsequent boost with GD2 ganglioside. The vaccine-induced antibodies recognized GD2-positive tumor cells, mediated complement-dependent cytotoxicity, and exhibited protection against s.c. human GD2-positive melanoma growth in the severe combined immunodeficient mouse xenograft model. The results from our studies provide insights into approaches for boosting GD2 cross-reactive IgG antibody responses by minigene vaccination with a protective epitope of GD2 ganglioside.

Model of three-dimensional structure of VDR bound with Vitamin D3 analogs substituted at carbon-2.

All Vitamin D analogs possessing the A ring modified at C-2 and showing calcemic activities nest themselves in the VDR binding pocket, oriented towards Tyr 143. Such topology resembles the position of the Vitamin D hormone in hVDRmt [Proc. Natl. Acad. Sci. U.S.A. 98 (2001) 5491]. Conversely, inactive 2beta-methyl-19-nor-analogs anchor the receptor cavity in a distinguishably different manner, namely by their side chain. Moreover, these inactive vitamins have a different conformation around C(6)-C(7) bond. Topology of modeled complexes suggests that a Vitamin D analog will be biologically active if its intercyclic 5,7-diene moiety assumes parallel position to tryptophan aromatic rings; such orientation allows for creating pi-pi interactions. The broad comparison of calcemic activities of the analogs, and their interactions with VDR, revealed that specific hydrophobic contacts are involved in bone calcium mobilization (BCM). These contacts occur between 21-methyl group and a few amino acids (V296, L305 and L309), conserved in the nuclear receptor superfamily. In the inactive 2beta-methyl-19-nor analogs such contacts do not exist. We speculate that two hydrophobic receptor patches, being in close contact with ligand methyl groups, might influence interaction with co-modulators involved in calcium homeostasis.

Protein fragment reconstruction using various modeling techniques.

Recently developed reduced models of proteins with knowledge-based force fields have been applied to a specific case of comparative modeling. From twenty high resolution protein structures of various structural classes, significant fragments of their chains have been removed and treated as unknown. The remaining portions of the structures were treated as fixed - i.e., as templates with an exact alignment. Then, the missed fragments were reconstructed using several modeling tools. These included three reduced types of protein models: the lattice SICHO (Side Chain Only) model, the lattice CABS (Calpha + Cbeta + Side group) model and an off-lattice model similar to the CABS model and called REFINER. The obtained reduced models were compared with more standard comparative modeling tools such as MODELLER and the SWISS-MODEL server. The reduced model results are qualitatively better for the higher resolution lattice models, clearly suggesting that these are now mature, competitive and complementary (in the range of sparse alignments) to the classical tools of comparative modeling. Comparison between the various reduced models strongly suggests that the essential ingredient for the sucessful and accurate modeling of protein structures is not the representation of conformational space (lattice, off-lattice, all-atom) but, rather, the specificity of the force fields used and, perhaps, the sampling techniques employed. These conclusions are encouraging for the future application of the fast reduced models in comparative modeling on a genomic scale.

2-Ethyl and 2-ethylidene analogues of 1alpha,25-dihydroxy-19-norvitamin D(3): synthesis, conformational analysis, biological activities, and docking to the modeled rVDR ligand binding domain.

Novel 19-nor analogues of 1alpha,25-dihydroxyvitamin D(3) were prepared and substituted at C-2 with an ethylidene group. The synthetic pathway was via Wittig-Horner coupling of the corresponding A-ring phosphine oxides with the protected 25-hydroxy Grundmann's ketones. Selective catalytic hydrogenation of 2-ethylidene analogues provided the 2alpha- and 2beta-ethyl compounds. The 2-ethylidene-19-nor compounds with a methyl group from the ethylidene moiety in a trans relationship to the C(6)-C(7) bond (E-isomers) were more potent than the corresponding Z-isomers and the natural hormone in binding to the vitamin D receptor. Both geometrical isomers (E and Z) of (20S)-2-ethylidene-19-norvitamin D(3) and both 2alpha-ethyl-19-norvitamins (in the 20R- and 20S-series) have much higher HL-60 differentiation activity than does 1alpha,25-(OH)(2)D(3). Both E-isomers (20R and 20S) of 2-ethylidene vitamins are characterized by very high calcemic activity in rats. The three-dimensional structure model of the rat vitamin D receptor and the computational docking of four synthesized (20R)-19-norvitamin D(3) analogues into its binding pocket are also reported.