EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization.

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.

The C-terminal domain of Plasmodium falciparum merozoite surface protein 3 self-assembles into alpha-helical coiled coil tetramer.

Proteins located on the surface of the pathogenic malaria parasite Plasmodium falciparum are objects of intensive studies due to their important role in the invasion of human cells and the accessibility to host antibodies thus making these proteins attractive vaccine candidates. One of these proteins, merozoite surface protein 3 (MSP3) represents a leading component among vaccine candidates; however, little is known about its structure and function. Our biophysical studies suggest that the 40 residue C-terminal domain of MSP3 protein self-assembles into a four-stranded alpha-helical coiled coil structure where alpha-helices are packed "side-by-side". A bioinformatics analysis provides an extended list of known and putative proteins from different species of Plasmodium which have such MSP3-like C-terminal domains. This finding allowed us to extend some conclusions of our studies to a larger group of the malaria surface proteins. Possible structural and functional roles of these highly conserved oligomerization domains in the intact merozoite surface proteins are discussed.

Activation of peroxisome proliferator-activated receptors (PPARs) by their ligands and protein kinase A activators.

The nuclear peroxisome proliferator-activated receptors (PPARs) alpha, beta, and gamma activate the transcription of multiple genes involved in lipid metabolism. Several natural and synthetic ligands have been identified for each PPAR isotype but little is known about the phosphorylation state of these receptors. We show here that activators of protein kinase A (PKA) can enhance mouse PPAR activity in the absence and the presence of exogenous ligands in transient transfection experiments. Activation function 1 (AF-1) of PPARs was dispensable for transcriptional enhancement, whereas activation function 2 (AF-2) was required for this effect. We also show that several domains of PPAR can be phosphorylated by PKA in vitro. Moreover, gel retardation experiments suggest that PKA stabilizes binding of the liganded PPAR to DNA. PKA inhibitors decreased not only the kinase-dependent induction of PPARs but also their ligand-dependent induction, suggesting an interaction between both pathways that leads to maximal transcriptional induction by PPARs. Moreover, comparing PPAR alpha knockout (KO) with PPAR alpha WT mice, we show that the expression of the acyl CoA oxidase (ACO) gene can be regulated by PKA-activated PPAR alpha in liver. These data demonstrate that the PKA pathway is an important modulator of PPAR activity, and we propose a model associating this pathway in the control of fatty acid beta-oxidation under conditions of fasting, stress, and exercise.

Signal peptide and helical bundle domains of virulent canine distemper virus fusion protein restrict fusogenicity.

Persistence in canine distemper virus (CDV) infection is correlated with very limited cell-cell fusion and lack of cytolysis induced by the neurovirulent A75/17-CDV compared to that of the cytolytic Onderstepoort vaccine strain. We have previously shown that this difference was at least in part due to the amino acid sequence of the fusion (F) protein (P. Plattet, J. P. Rivals, B. Zuber, J. M. Brunner, A. Zurbriggen, and R. Wittek, Virology 337:312-326, 2005). Here, we investigated the molecular mechanisms of the neurovirulent CDV F protein underlying limited membrane fusion activity. By exchanging the signal peptide between both F CDV strains or replacing it with an exogenous signal peptide, we demonstrated that this domain controlled intracellular and consequently cell surface protein expression, thus indirectly modulating fusogenicity. In addition, by serially passaging a poorly fusogenic virus and selecting a syncytium-forming variant, we identified the mutation L372W as being responsible for this change of phenotype. Intriguingly, residue L372 potentially is located in the helical bundle domain of the F(1) subunit. We showed that this mutation drastically increased fusion activity of F proteins of both CDV strains in a signal peptide-independent manner. Due to its unique structure even among morbilliviruses, our findings with respect to the signal peptide are likely to be specifically relevant to CDV, whereas the results related to the helical bundle add new insights to our growing understanding of this class of F proteins. We conclude that different mechanisms involving multiple domains of the neurovirulent A75/17-CDV F protein act in concert to limit fusion activity, preventing lysis of infected cells, which ultimately may favor viral persistence.

Activation/division of lymphocytes results in increased levels of cytoplasmic activation/proliferation-associated protein-1: prototype of a new family of proteins.

We purified from activated T lymphocytes a novel, highly conserved, 116-kDa, intracellular protein that occurred at high levels in the large, dividing cells of the thymus, was up-regulated when resting T or B lymphocytes or hemopoietic progenitors were activated, and was down-regulated when a monocytic leukemia, M1, was induced to differentiate. Expression of the protein was highest in the thymus and spleen and lowest in tissues with a low proportion of dividing cells such as kidney or muscle, although expression was high in the brain. The protein was localized to the cytosol and was phosphorylated, which is consistent with a previous report that the Xenopus laevis ortholog was phosphorylated by a mitotically activated kinase (1 ). The cDNA was previously mischaracterized as encoding p137, a 137-kDa GPI-linked membrane protein (2 ). We propose that the authentic protein encoded by this cDNA be called cytoplasmic activation/proliferation-associated protein-1 (caprin-1), and show that it is the prototype of a novel family of proteins characterized by two novel protein domains, termed homology regions-1 and -2 (HR-1, HR-2). Although we have found evidence for caprins only in urochordates and vertebrates, two insect proteins exhibit well-conserved HR-1 domains. The HR-1 and HR-2 domains have no known function, although the HR-1 of caprin-1 appeared necessary for formation of multimeric complexes of caprin-1. Overexpression of a fusion protein of enhanced green fluorescent protein and caprin-1 induced a specific, dose-dependent suppression of the proliferation of NIH-3T3 cells, consistent with the notion that caprin-1 plays a role in cellular activation or proliferation.

Characterization of the phosphorylation sites involved in G protein-coupled receptor kinase- and protein kinase C-mediated desensitization of the alpha1B-adrenergic receptor.

Catecholamines as well as phorbol esters can induce the phosphorylation and desensitization of the alpha1B-adrenergic receptor (alpha1BAR). In this study, phosphoamino acid analysis of the phosphorylated alpha1BAR revealed that both epinephrine- and phorbol ester-induced phosphorylation predominantly occurs at serine residues of the receptor. The findings obtained with receptor mutants in which portions of the C-tail were truncated or deleted indicated that a region of 21 amino acids (393-413) of the carboxyl terminus including seven serines contains the main phosphorylation sites involved in agonist- as well as phorbol ester-induced phosphorylation and desensitization of the alpha1BAR. To identify the serines invoved in agonist- versus phorbol ester-dependent regulation of the receptor, two different strategies were adopted, the seven serines were either substituted with alanine or reintroduced into a mutant lacking all of them. Our findings indicate that Ser394 and Ser400 were phosphorylated following phorbol ester-induced activation of protein kinase C, whereas Ser404, Ser408, and Ser410 were phosphorylated upon stimulation of the alpha1BAR with epinephrine. The observation that overexpression of G protein-coupled kinase 2 (GRK2) could increase agonist-induced phosphorylation of Ser404, Ser408, and Ser410, strongly suggests that these serines are the phosphorylation sites of the alpha1BAR for kinases of the GRK family. Phorbol ester-induced phosphorylation of the Ser394 and Ser400 as well as GRK2-mediated phosphorylation of the Ser404, Ser408, and Ser410, resulted in the desensitization of alpha1BAR-mediated inositol phosphate response. This study provides generalities about the biochemical mechanisms underlying homologous and heterologous desensitization of G protein-coupled receptors linked to the activation of phospholipase C.

Aminoterminal amphipathic α-helix AH1 of hepatitis C virus nonstructural protein 4B possesses a dual role in RNA replication and virus production.

Nonstructural protein 4B (NS4B) is a key organizer of hepatitis C virus (HCV) replication complex formation. In concert with other nonstructural proteins, it induces a specific membrane rearrangement, designated as membranous web, which serves as a scaffold for the HCV replicase. The N-terminal part of NS4B comprises a predicted and a structurally resolved amphipathic α-helix, designated as AH1 and AH2, respectively. Here, we report a detailed structure-function analysis of NS4B AH1. Circular dichroism and nuclear magnetic resonance structural analyses revealed that AH1 folds into an amphipathic α-helix extending from NS4B amino acid 4 to 32, with positively charged residues flanking the helix. These residues are conserved among hepaciviruses. Mutagenesis and selection of pseudorevertants revealed an important role of these residues in RNA replication by affecting the biogenesis of double-membrane vesicles making up the membranous web. Moreover, alanine substitution of conserved acidic residues on the hydrophilic side of the helix reduced infectivity without significantly affecting RNA replication, indicating that AH1 is also involved in virus production. Selective membrane permeabilization and immunofluorescence microscopy analyses of a functional replicon harboring an epitope tag between NS4B AH1 and AH2 revealed a dual membrane topology of the N-terminal part of NS4B during HCV RNA replication. Luminal translocation was unaffected by the mutations introduced into AH1, but was abrogated by mutations introduced into AH2. In conclusion, our study reports the three-dimensional structure of AH1 from HCV NS4B, and highlights the importance of positively charged amino acid residues flanking this amphipathic α-helix in membranous web formation and RNA replication. In addition, we demonstrate that AH1 possesses a dual role in RNA replication and virus production, potentially governed by different topologies of the N-terminal part of NS4B.

Characterisation of the putative effector interaction site of the regulatory HbpR protein from Pseudomonas azelaica by site-directed mutagenesis.

Bacterial transcription activators of the XylR/DmpR subfamily exert their expression control via σ(54)-dependent RNA polymerase upon stimulation by a chemical effector, typically an aromatic compound. Where the chemical effector interacts with the transcription regulator protein to achieve activation is still largely unknown. Here we focus on the HbpR protein from Pseudomonas azelaica, which is a member of the XylR/DmpR subfamily and responds to biaromatic effectors such as 2-hydroxybiphenyl. We use protein structure modeling to predict folding of the effector recognition domain of HbpR and molecular docking to identify the region where 2-hydroxybiphenyl may interact with HbpR. A large number of site-directed HbpR mutants of residues in- and outside the predicted interaction area was created and their potential to induce reporter gene expression in Escherichia coli from the cognate P(C) promoter upon activation with 2-hydroxybiphenyl was studied. Mutant proteins were purified to study their conformation. Critical residues for effector stimulation indeed grouped near the predicted area, some of which are conserved among XylR/DmpR subfamily members in spite of displaying different effector specificities. This suggests that they are important for the process of effector activation, but not necessarily for effector specificity recognition.

Mutations in the TGFβ binding-protein-like domain 5 of FBN1 are responsible for acromicric and geleophysic dysplasias.

Geleophysic (GD) and acromicric dysplasia (AD) belong to the acromelic dysplasia group and are both characterized by severe short stature, short extremities, and stiff joints. Although AD has an unknown molecular basis, we have previously identified ADAMTSL2 mutations in a subset of GD patients. After exome sequencing in GD and AD cases, we selected fibrillin 1 (FBN1) as a candidate gene, even though mutations in this gene have been described in Marfan syndrome, which is characterized by tall stature and arachnodactyly. We identified 16 heterozygous FBN1 mutations that are all located in exons 41 and 42 and encode TGFβ-binding protein-like domain 5 (TB5) of FBN1 in 29 GD and AD cases. Microfibrillar network disorganization and enhanced TGFβ signaling were consistent features in GD and AD fibroblasts. Importantly, a direct interaction between ADAMTSL2 and FBN1 was demonstrated, suggesting a disruption of this interaction as the underlying mechanism of GD and AD phenotypes. Although enhanced TGFβ signaling caused by FBN1 mutations can trigger either Marfan syndrome or GD and AD, our findings support the fact that TB5 mutations in FBN1 are responsible for short stature phenotypes.

Topological aspects of DNA function and protein folding.

The Topological Aspects of DNA Function and Protein Folding international meeting provided an interdisciplinary forum for biological scientists, physicists and mathematicians to discuss recent developments in the application of topology to the study of DNA and protein structure. It had 111 invited participants, 48 talks and 21 posters. The present article discusses the importance of topology and introduces the articles from the meeting's speakers.

The UlaG protein familly defines novel structural and functional motifs grafted on an ancient RNase fold

Background: Bacterial populations are highly successful at colonizing new habitats and adapting to changing environmental conditions, partly due to their capacity to evolve novel virulence and metabolic pathways in response to stress conditions and to shuffle them by horizontal gene transfer (HGT). A common theme in the evolution of new functions consists of gene duplication followed by functional divergence. UlaG, a unique manganese-dependent metallo-b-lactamase (MBL) enzyme involved in L-ascorbate metabolism by commensal and symbiotic enterobacteria, provides a model for the study of the emergence of new catalytic activities from the modification of an ancient fold. Furthermore, UlaG is the founding member of the so-called UlaG-like (UlaGL) protein family, a recently established and poorly characterized family comprising divalent (and perhaps trivalent)metal-binding MBLs that catalyze transformations on phosphorylated sugars and nucleotides. Results: Here we combined protein structure-guided and sequence-only molecular phylogenetic analyses to dissect the molecular evolution of UlaG and to study its phylogenomic distribution, its relatedness with present-day UlaGL protein sequences and functional conservation. Phylogenetic analyses indicate that UlaGL sequences are present in Bacteria and Archaea, with bona fide orthologs found mainly in mammalian and plant-associated Gramnegative and Gram-positive bacteria. The incongruence between the UlaGL tree and known species trees indicates exchange by HGT and suggests that the UlaGL-encoding genes provided a growth advantage under changing conditions. Our search for more distantly related protein sequences aided by structural homology has uncovered that UlaGL sequences have a common evolutionary origin with present-day RNA processing and metabolizing MBL enzymes widespread in Bacteria, Archaea, and Eukarya. This observation suggests an ancient origin for the UlaGL family within the broader trunk of the MBL superfamily by duplication, neofunctionalization and fixation. Conclusions: Our results suggest that the forerunner of UlaG was present as an RNA metabolizing enzyme in the last common ancestor, and that the modern descendants of that ancestral gene have a wide phylogenetic distribution and functional roles. We propose that the UlaGL family evolved new metabolic roles among bacterial and possibly archeal phyla in the setting of a close association with metazoans, such as in the mammalian gastrointestinal tract or in animal and plant pathogens, as well as in environmental settings. Accordingly, the major evolutionary forces shaping the UlaGL family include vertical inheritance and lineage-specific duplication and acquisition of novel metabolic functions, followed by HGT and numerous lineage-specific gene loss events.

Bone morphogenetic proteins: from structure to clinical use

Bone morphogenetic proteins (BMPs) are multi-functional growth factors belonging to the transforming growth factor ß superfamily. Family members are expressed during limb development, endochondral ossification, early fracture, and cartilage repair. The activity of BMPs was first identified in the 1960s but the proteins responsible for bone induction were unknown until the purification and cloning of human BMPs in the 1980s. To date, about 15 BMP family members have been identified and characterized. The signal triggered by BMPs is transduced through serine/threonine kinase receptors, type I and II subtypes. Three type I receptors have been shown to bind BMP ligands, namely: type IA and IB BMP receptors and type IA activin receptors. BMPs seem to be involved in the regulation of cell proliferation, survival, differentiation and apoptosis, but their hallmark is their ability to induce bone, cartilage, ligament, and tendon formation at both heterotopic and orthotopic sites. This suggests that, in the future, they may play a major role in the treatment of bone diseases. Several animal studies have illustrated the potential of BMPs to enhance spinal fusion, repair critical-size defects, accelerate union, and heal articular cartilage lesions. Difficulties in producing and purifying BMPs from bone tissue have prompted the attempts made by several laboratories, including ours, to express these proteins in the recombinant form in heterologous systems. This review focuses on BMP structure, molecular mechanisms of action and significance and potential applications in medical, dental and veterinary practice for the treatment of cartilage and bone-related diseases.

Benchmarking consensus model quality assessment for protein fold recognition

Background: Selecting the highest quality 3D model of a protein structure from a number of alternatives remains an important challenge in the field of structural bioinformatics. Many Model Quality Assessment Programs (MQAPs) have been developed which adopt various strategies in order to tackle this problem, ranging from the so called "true" MQAPs capable of producing a single energy score based on a single model, to methods which rely on structural comparisons of multiple models or additional information from meta-servers. However, it is clear that no current method can separate the highest accuracy models from the lowest consistently. In this paper, a number of the top performing MQAP methods are benchmarked in the context of the potential value that they add to protein fold recognition. Two novel methods are also described: ModSSEA, which based on the alignment of predicted secondary structure elements and ModFOLD which combines several true MQAP methods using an artificial neural network. Results: The ModSSEA method is found to be an effective model quality assessment program for ranking multiple models from many servers, however further accuracy can be gained by using the consensus approach of ModFOLD. The ModFOLD method is shown to significantly outperform the true MQAPs tested and is competitive with methods which make use of clustering or additional information from multiple servers. Several of the true MQAPs are also shown to add value to most individual fold recognition servers by improving model selection, when applied as a post filter in order to re-rank models. Conclusion: MQAPs should be benchmarked appropriately for the practical context in which they are intended to be used. Clustering based methods are the top performing MQAPs where many models are available from many servers; however, they often do not add value to individual fold recognition servers when limited models are available. Conversely, the true MQAP methods tested can often be used as effective post filters for re-ranking few models from individual fold recognition servers and further improvements can be achieved using a consensus of these methods.

FunFOLDQA: a quality assessment tool for protein-ligand binding site residue predictions

The estimation of prediction quality is important because without quality measures, it is difficult to determine the usefulness of a prediction. Currently, methods for ligand binding site residue predictions are assessed in the function prediction category of the biennial Critical Assessment of Techniques for Protein Structure Prediction (CASP) experiment, utilizing the Matthews Correlation Coefficient (MCC) and Binding-site Distance Test (BDT) metrics. However, the assessment of ligand binding site predictions using such metrics requires the availability of solved structures with bound ligands. Thus, we have developed a ligand binding site quality assessment tool, FunFOLDQA, which utilizes protein feature analysis to predict ligand binding site quality prior to the experimental solution of the protein structures and their ligand interactions. The FunFOLDQA feature scores were combined using: simple linear combinations, multiple linear regression and a neural network. The neural network produced significantly better results for correlations to both the MCC and BDT scores, according to Kendall’s τ, Spearman’s ρ and Pearson’s r correlation coefficients, when tested on both the CASP8 and CASP9 datasets. The neural network also produced the largest Area Under the Curve score (AUC) when Receiver Operator Characteristic (ROC) analysis was undertaken for the CASP8 dataset. Furthermore, the FunFOLDQA algorithm incorporating the neural network, is shown to add value to FunFOLD, when both methods are employed in combination. This results in a statistically significant improvement over all of the best server methods, the FunFOLD method (6.43%), and one of the top manual groups (FN293) tested on the CASP8 dataset. The FunFOLDQA method was also found to be competitive with the top server methods when tested on the CASP9 dataset. To the best of our knowledge, FunFOLDQA is the first attempt to develop a method that can be used to assess ligand binding site prediction quality, in the absence of experimental data.

Predicting protein structures and structural annotation of proteomes

Protein structure prediction methods aim to predict the structures of proteins from their amino acid sequences, utilizing various computational algorithms. Structural genome annotation is the process of attaching biological information to every protein encoded within a genome via the production of three-dimensional protein models.