Sébastien Fiorucci

Binding site prediction and improved scoring during flexible protein–protein docking with ATTRACT†

The ATTRACT protein–protein docking program combined with a coarse‐grained protein model has been used to predict protein–protein complex structures in CAPRI rounds 13–19. For six targets acceptable or better quality solutions have been submitted (high quality predictions for targets 32, 40, 41, and 42). The improved performance compared to previous rounds can be attributed in part to the inclusion of conformational flexibility during systematic searches and an optimized scoring function. In addition, a recently developed method for the prediction of putative protein binding sites based on the electrostatic penalty to place neutral low dielectric probes on the protein surface was applied to the most recent targets. The approach resulted in useful predictions of putative binding sites that can help to limit the systematic docking searches. Possible improvements of the docking approach in particular at the scoring and refinement steps are discussed. Proteins 2010.

Prediction of Protein-Protein Interaction Sites Using Electrostatic Desolvation Profiles

Protein-protein complex formation involves removal of water from the interface region. Surface regions with a small free energy penalty for water removal or desolvation may correspond to preferred interaction sites. A method to calculate the electrostatic free energy of placing a neutral low-dielectric probe at various protein surface positions has been designed and applied to characterize putative interaction sites. Based on solutions of the finite-difference Poisson equation, this method also includes long-range electrostatic contributions and the protein solvent boundary shape in contrast to accessible-surface-area-based solvation energies. Calculations on a large set of proteins indicate that in many cases (>90%), the known binding site overlaps with one of the six regions of lowest electrostatic desolvation penalty (overlap with the lowest desolvation region for 48% of proteins). Since the onset of electrostatic desolvation occurs even before direct protein-protein contact formation, it may help guide proteins toward the binding region in the final stage of complex formation. It is interesting that the probe desolvation properties associated with residue types were found to depend to some degree on whether the residue was outside of or part of a binding site. The probe desolvation penalty was on average smaller if the residue was part of a binding site compared to other surface locations. Applications to several antigen-antibody complexes demonstrated that the approach might be useful not only to predict protein interaction sites in general but to map potential antigenic epitopes on protein surfaces.

Exploiting Antigenic Diversity for Vaccine Design: THE CHLAMYDIA ArtJ PARADIGM*

We present an interdisciplinary approach that, by incorporating a range of experimental and computational techniques, allows the identification and characterization of functional/immunogenic domains. This approach has been applied to ArtJ, an arginine-binding protein whose orthologs in Chlamydiae trachomatis (CT ArtJ) and pneumoniae (CPn ArtJ) are shown to have different immunogenic properties despite a high sequence similarity (60% identity). We have solved the crystallographic structures of CT ArtJ and CPn ArtJ, which are found to display a type II transporter fold organized in two α-β domains with the arginine-binding region at their interface. Although ArtJ is considered to belong to the periplasm, we found that both domains contain regions exposed on the bacterial surface. Moreover, we show that recombinant ArtJ binds to epithelial cells in vitro, suggesting a role for ArtJ in host-cell adhesion during Chlamydia infection. Experimental epitope mapping and computational analysis of physicochemical determinants of antibody recognition revealed that immunogenic epitopes reside mainly in the terminal (D1) domain of both CPn and CT ArtJ, whereas the surface properties of the respective binding-prone regions appear sufficiently different to assume divergent immunogenic behavior. Neutralization assays revealed that sera raised against CPn ArtJ D1 partially reduce both CPn and CT infectivity in vitro, suggesting that functional antibodies directed against this domain may potentially impair chlamydial infectivity. These findings suggest that the approach presented here, combining functional and structure-based analyses of evolutionary-related antigens can be a valuable tool for the identification of cross-species immunogenic epitopes for vaccine development.

PTools: an opensource molecular docking library.

BackgroundMacromolecular docking is a challenging field of bioinformatics. Developing new algorithms is a slow process generally involving routine tasks that should be found in a robust library and not programmed from scratch for every new software application.ResultsWe present an object-oriented Python/C++ library to help the development of new docking methods. This library contains low-level routines like PDB-format manipulation functions as well as high-level tools for docking and analyzing results. We also illustrate the ease of use of this library with the detailed implementation of a 3-body docking procedure.ConclusionThe PTools library can handle molecules at coarse-grained or atomic resolution and allows users to rapidly develop new software. The library is already in use for protein-protein and protein-DNA docking with the ATTRACT program and for simulation analysis. This library is freely available under the GNU GPL license, together with detailed documentation.

Mechanistic Events Underlying Odorant Binding Protein Chemoreception

Odorant binding proteins (OBPs) are small hydrophilic proteins, belonging to the lipocalin family dedicated to bind and transport small hydrophobic ligands. Despite many works, the mechanism of ligand binding, together with the functional role of these proteins remains a topic of debate and little is known at the atomic level. The present work reports a computational study of odorants capture and release by an OBP, using both constrained and unconstrained simulations, giving a glimpse on the molecular mechanism of chemoreception. The residues at the origin of the regulation of the protein door opening are identified and a tyrosine amino‐acid together with other nearby residues appear to play a crucial role in allowing this event to occur. The simulations reveal that this tyrosine and the proteins L5 loop are implicated in the ligand contact with the protein and act as an anchoring point for the ligand. The protein structural features required for the ligand entry are highly conserved among many transport proteins, suggesting that this mechanism could somewhat be extended to some members of the larger family of lipocalin. Proteins 2007.

Isolation and functional characterization of a τ-cadinol synthase, a new sesquiterpene synthase from Lavandula angustifolia

Abstract In this paper we characterize three sTPSs: a germacrene D (LaGERDS), a (E)-β-caryophyllene (LaCARS) and a τ-cadinol synthase (LaCADS). τ-cadinol synthase is reported here for the first time and its activity was studied in several biological models including transiently or stably transformed tobacco species. Three dimensional structure models of LaCADS and Ocimum basilicum γ-cadinene synthase were built by homology modeling using the template structure of Gossypium arboreum δ-cadinene synthase. The depiction of their active site organization provides evidence of the global influence of the enzymes on the formation of τ-cadinol: instead of a unique amino-acid, the electrostatic properties and solvent accessibility of the whole active site in LaCADS may explain the stabilization of the cadinyl cation intermediate. Quantitative PCR performed from leaves and inflorescences showed two patterns of expression. LaGERDS and LaCARS were mainly expressed during early stages of flower development and, at these stages, transcript levels paralleled the accumulation of the corresponding terpene products (germacrene D and (E)-β-caryophyllene). By contrast, the expression level of LaCADS was constant in leaves and flowers. Phylogenetic analysis provided informative results on potential duplication process leading to sTPS diversification in lavender.

Binding free energy prediction in strongly hydrophobic biomolecular systems.

We present a comparison of various computational approaches aiming at predicting the binding free energy in ligand-protein systems where the ligand is located within a highly hydrophobic cavity. The relative binding free energy between similar ligands is obtained by means of the thermodynamic integration (TI) method and compared to experimental data obtained through isothermal titration calorimetry measurements. The absolute free energy of binding prediction was obtained on a similar system (a pyrazine derivative bound to a lipocalin) by TI, potential of mean force (PMF) and also by means of the MMPBSA protocols. Although the TI protocol performs poorly either with an explicit or an implicit solvation scheme, the PMF calculation using an implicit solvation scheme leads to encouraging results, with a prediction of the binding affinity being 2 kcal mol(-1) lower than the experimental value. The use of an implicit solvation scheme appears to be well suited for the study of such hydrophobic systems, due to the lack of water molecules within the binding site.

Molecular simulations reveal a new entry site in quercetin 2,3-dioxygenase. A pathway for dioxygen?

Molecular dynamics simulations performed on quercetin 2,3‐dioxygenase have shown the existence of a channel linking the bulk solvent and the cavity of the enzyme. Although much is known about the the oxygenolysis reaction catalyzed by this enzyme, the way dioxygen enters the active site has not been firmly established. The size, orientation and hydrophobic character of this channel suggests that it could provide an entrance for molecular dioxygen into the cavity. Free energy calculations show that such a process is likely to occur. Proteins 2006.

ATTRACT and PTools: open source programs for protein-protein docking.

The prediction of the structure of protein-protein complexes based on structures or structural models of isolated partners is of increasing importance for structural biology and bioinformatics. The ATTRACT program can be used to perform systematic docking searches based on docking energy minimization. It is part of the object-oriented PTools library written in Python and C++. The library contains various routines to manipulate protein structures, to prepare and perform docking searches as well as analyzing docking results. It also intended to facilitate further methodological developments in the area of macromolecular docking that can be easily integrated. Here, we describe the application of PTools to perform systematic docking searches and to analyze the results. In addition, the possibility to perform multi-component docking will also be presented.

Molecular simulations bring new insights into flavonoid/quercetinase interaction modes

Molecular dynamics simulations, using the AMBER force field, were performed to study Quercetin 2,3‐Dioxygenase enzyme (Quercetinase or 2,3QD). We have analyzed the structural modifications of the active site and of the linker region between the native enzyme and the enzyme‐substrate complex. New structural informations, such as an allosteric effect in the presence of the substrate, as well as description of the enzyme‐substrate interactions and values of binding free energies were brought out. All these results confirm the idea that the linker encloses the substrate in the active site and also enlighten the recognition role of the substrate B‐ring by the enzyme. Moreover, a specific interaction scheme has been proposed to explain the relative degradation rate of various flavonoid compounds under the oxygenolysis reaction catalyzed by the Quercetin 2,3‐Dioxygenase enzyme. Proteins 2007.