Blind docking of 260 protein-ligand complexes with EADock 2.0.

Molecular docking softwares are one of the important tools of modern drug development pipelines. The promising achievements of the last 10 years emphasize the need for further improvement, as reflected by several recent publications (Leach et al., J Med Chem 2006, 49, 5851; Warren et al., J Med Chem 2006, 49, 5912). Our initial approach, EADock, showed a good performance in reproducing the experimental binding modes for a set of 37 different ligand-protein complexes (Grosdidier et al., Proteins 2007, 67, 1010). This article presents recent improvements regarding the scoring and sampling aspects over the initial implementation, as well as a new seeding procedure based on the detection of cavities, opening the door to blind docking with EADock. These enhancements were validated on 260 complexes taken from the high quality Ligand Protein Database [LPDB, (Roche et al., J Med Chem 2001, 44, 3592)]. Two issues were identified: first, the quality of the initial structures cannot be assumed and a manual inspection and/or a search in the literature are likely to be required to achieve the best performance. Second the description of interactions involving metal ions still has to be improved. Nonetheless, a remarkable success rate of 65% was achieved for a large scale blind docking assay, when considering only the top ranked binding mode and a success threshold of 2 A RMSD to the crystal structure. When looking at the five-top ranked binding modes, the success rate increases up to 76%. In a standard local docking assay, success rates of 75 and 83% were obtained, considering only the top ranked binding mode, or the five top binding modes, respectively.

Computational study of ligand binding in lipid transfer proteins: Structures, interfaces, and free energies of protein-lipid complexes

Plant nonspecific lipid transfer proteins (nsLTPs) bind a wide variety of lipids, which allows them to perform disparate functions. Recent reports on their multifunctionality in plant growth processes have posed new questions on the versatile binding abilities of these proteins. The lack of binding specificity has been customarily explained in qualitative terms on the basis of a supposed structural flexibility and nonspecificity of hydrophobic protein-ligand interactions. We present here a computational study of protein-ligand complexes formed between five nsLTPs and seven lipids bound in two different ways in every receptor protein. After optimizing geometries inmolecular dynamics calculations, we computed Poisson- Boltzmann electrostatic potentials, solvation energies, properties of the protein-ligand interfaces, and estimates of binding free energies of the resulting complexes. Our results provide the first quantitative information on the ligand abilities of nsLTPs, shed new light into protein-lipid interactions, and reveal new features which supplement commonly held assumptions on their lack of binding specificity.

Probabilistic contact preferences in protein-ligand and protein-protein complexes

Modeller för intermolekulär växelvärkan utnyttjas brett inom biologin. Analys av kontakter mellan proteiner och läkemedelsforskning representerar typiska tillämpningsområden för dylika modeller. En modell som beskriver sådana molekylära växelverkningar kan utformas med hjälp av biofysisk teori, vilket tenderar att resultera i ytterst tung beräkningsbörda även för enkla tillämpningar. Ett alternativt sätt att formulera modeller är att utnyttja stora databaser som innehåller strukturmätningar gjorda med hjälp av till exempel röntgendiffraktion. Då man använder sig av empiriska mätdata direkt, möjliggör en statistisk modell att osäkerheten och inexaktheten i datat tas till hänsyn på ett adekvat sätt, samtidigt som beräkningsbördan håller sig på en rimligare nivå jämfört med kvantmekaniska metoder som i princip borde ge de optimala resultaten. I avhandlingen utvecklades en 3D modell för numerisk undersökning av intermolekulär växelverkan baserad på Bayesiansk statistik. Modellens syfte är att åstadkomma prognoser för det hurdana eller vilka molekylstrukturer prefereras i en given kontext, d.v.s. är mer sannolika inom ramen för interaktion. Modellen testades i essentiella molekyläromgivningar - en liten molekyl vid sin bindningsplats hos ett protein och en gränsyta mellan proteinerna i ett komplex. De erhållna numeriska resultaten motsvarar väl experimentella resultat som tidigare rapporterats i litteraturen, exempelvis kvalitativa bindningsaffiniteter och kemisk kännedom av vissa aminosyrors rumsliga förmågor att utgöra bindningar. I avhandlingen gjordes ytterligare preliminära tester av den statistiska ansatsen för modellering av den centrala molekylära strukturella anpassningsbarheten. I praktiken är den utvecklade modellen ämnad som ett led i en mer omfattande analysmetod, så som en s.k. farmakofor modell. Molekyylivuorovaikutusten mallintamista hyödynnetään laajasti biologisten kysymysten tarkastelussa. Tyypillisiä esimerkkejä sovelluskohteista ovat proteiinien väliset kontaktit ja lääkesuunnittelu. Vuorovaikutuksia kuvaavan mallin lähtökohta voi olla molekyyleihin liittyvä teoria, jolloin soveltamiseen liittyvä laskenta saattaa olla erityisen raskasta, tai suuri havaintojoukko joka on saatu aikaan esimerkiksi mittaamalla rakenteita röntgendiffraktio menetelmällä. Tilastollinen malli mahdollistaa havaintoaineistossa olevan epätarkkuuden ja epävarmuuden huomioimisen, samalla pitäen laskennallisen kuorman pienempänä verrattuna periaatteessa parhaan tuloksen antavaan kvanttimekaaniseen mallinnukseen. Väitöstyössä kehitettiin bayesiläiseen tilastotieteeseen perustuva 3D malli molekyylien välisten vuorovaikutusten laskennalliseen tarkasteluun. Mallin tehtävä on tuottaa ennusteita sen suhteen, minkä tai millaisten molekyylirakenteiden väliset kompleksit ovat etusijalla, toisin sanoen todennäköisempiä, vuorovaikutustilanteessa. Työssä kehitetyn menetelmän toimivuutta testattiin käyttötarkoituksen suhteen olennaisissa molekyyliympäristöissä - pieni molekyyli sitoutumiskohdassaan proteiinissa sekä rajapinta kahden proteiinin välilllä proteiinikompleksissa. Saadut laskennalliset tulokset vastasivat hyvin vertailuun käytettyjä kirjallisuudesta saatuja kokeellisia tuloksia, kuten laadullisia sitoutumisaffiniteetteja, sekä kemiallista tietoa esimerkiksi tiettyjen aminohappojen avaruudellisesta sidoksenmuodostuksesta. Väitöstyössä myös alustavasti testattiin tilastollista lähestymistapaa tärkeän molekyylien rakenteellisen mukautuvuuden mallintamiseen. Käytännössä malli on tarkoitettu osaksi jotakin laajempaa analyysimenetelmää, kuten farmakoforimallia.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Systematic analysis of protein-detergent complexes applying dynamic light scattering to optimize solutions for crystallization trials

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Development and Application of Covalent-Labeling Strategies for the Large-Scale Thermodynamic Analysis of Protein Folding and Ligand Binding

Thermodynamic stability measurements on proteins and protein-ligand complexes can offer insights not only into the fundamental properties of protein folding reactions and protein functions, but also into the development of protein-directed therapeutic agents to combat disease. Conventional calorimetric or spectroscopic approaches for measuring protein stability typically require large amounts of purified protein. This requirement has precluded their use in proteomic applications. Stability of Proteins from Rates of Oxidation (SPROX) is a recently developed mass spectrometry-based approach for proteome-wide thermodynamic stability analysis. Since the proteomic coverage of SPROX is fundamentally limited by the detection of methionine-containing peptides, the use of tryptophan-containing peptides was investigated in this dissertation. A new SPROX-like protocol was developed that measured protein folding free energies using the denaturant dependence of the rate at which globally protected tryptophan and methionine residues are modified with dimethyl (2-hydroxyl-5-nitrobenzyl) sulfonium bromide and hydrogen peroxide, respectively. This so-called Hybrid protocol was applied to proteins in yeast and MCF-7 cell lysates and achieved a ~50% increase in proteomic coverage compared to probing only methionine-containing peptides. Subsequently, the Hybrid protocol was successfully utilized to identify and quantify both known and novel protein-ligand interactions in cell lysates. The ligands under study included the well-known Hsp90 inhibitor geldanamycin and the less well-understood omeprazole sulfide that inhibits liver-stage malaria. In addition to protein-small molecule interactions, protein-protein interactions involving Puf6 were investigated using the SPROX technique in comparative thermodynamic analyses performed on wild-type and Puf6-deletion yeast strains. A total of 39 proteins were detected as Puf6 targets and 36 of these targets were previously unknown to interact with Puf6. Finally, to facilitate the SPROX/Hybrid data analysis process and minimize human errors, a Bayesian algorithm was developed for transition midpoint assignment. In summary, the work in this dissertation expanded the scope of SPROX and evaluated the use of SPROX/Hybrid protocols for characterizing protein-ligand interactions in complex biological mixtures.

pH Influence on the stability of foams with protein-polysaccharide complexes at their interfaces

Food foams such as marshmallow, Chantilly and mousses have behavior and stability directly connected with their microstructure, bubble size distribution and interfacial properties. A high interfacial tension inherent to air/liquid foams interfaces affects its stability, and thus it has a direct impact on processing, storage and product handling. In this work, the interactions of egg albumin with various types of polysaccharides were investigated by drop tensiometry, interfacial rheology and foam stability. The progressive addition of egg albumin and polysaccharide in water induced a drop of the air-water surface tension which was dependent on the pH and polysaccharide type. At pH 4, that is below the isoeletric point of egg albumen (pI = 4.5) the surface tension was decreased from 70 mN/m to 42 mN/m by the presence of the protein, and from 70 mN/m to 43 mN/m, 40 mN/m and 38 mN/m by subsequent addition of xanthan, guar gum and kappa-carrageenan, respectively. At pH 7.5 the surface tension was decreased from 70 mN/m to 43 mN/m by the simultaneous presence of the protein and kappa-carrageenan. However, a higher surface tension of 48 and 50 mN/m was found when xanthan and guar gum were added, respectively, when compared with carrageenan addition. The main role on the stabilization of protein-polysaccharide stabilized interfaces was identified on the elasticity of the interface. Foam stability experiments confirmed that egg-albumin/kappa-carrageenan at pH below the protein isoeletric point are the most efficient systems to stabilize air/water interfaces. These results clearly indicate that protein-polysaccharide coacervation at the air/water interface is an efficient process to increase foam stability. (C) 2009 Elsevier Ltd. All rights reserved.

Protein-ligand docking study: diterpenes from Juniperus brevifolia as anticancer and antimicrobial agentes

REDCAT: Natural Products and related Redox Catalysts: Basic Research and Applications in Medicine and Agriculture, Aveiro, 25-27 Novembro de 2012.

Expression of human estrogen receptor mutants in Xenopus oocytes: correlation between transcriptional activity and ability to form protein-DNA complexes.

The human estrogen receptor (hER) is a trans-acting regulatory protein composed of a series of discrete functional domains. We have microinjected an hER expression vector (HEO) into Xenopus oocyte nuclei and demonstrate, using Western blot assay, that the hER is synthesized. When nuclear extracts from oocytes were prepared and incubated in the presence of a 2.7 kb DNA fragment comprising the 5' end of the vitellogenin gene B2, formation of estrogen-dependent complexes could be visualized by electron microscopy over the estrogen responsive element (ERE). Of crucial importance is the observation that the complex formation is inhibited by the estrogen antagonist tamoxifen, is restored by the addition of the hormone and does not take place with extracts from control oocytes injected with the expression vector lacking the sequences encoding the receptor. The presence of the biologically active hER is confirmed in co-injection experiments, in which HEO is co-introduced with a CAT reporter gene under the control of a vitellogenin promoter containing or lacking the ERE. CAT assays and primer extensions analyses reveal that both the receptor and the ERE are essential for estrogen induced stimulation of transcription. The same approach was used to analyze selective hER mutants. We find that the DNA binding domain (region C) is essential for protein--DNA complex formation at the ERE but is not sufficient by itself to activate transcription from the reporter gene. In addition to region C, both the hormone binding (region E) and amino terminal (region A/B) domains are needed for an efficient transcription activation.(ABSTRACT TRUNCATED AT 250 WORDS)

Pi*-pi* bonding interactions generated by halogen oxidation of zirconium(IV) redox-active ligand complexes.

The new complex, [Zr(pda)2]n (1, pda2- = N,N'-bis(neo-pentyl)-ortho-phenylenediamide, n = 1 or 2), prepared by the reaction of 2 equiv of pdaLi2 with ZrCl4, reacts rapidly with halogen oxidants to afford the new product ZrX2(disq)2 (3, X = Cl, Br, I; disq- = N,N'-bis(neo-pentyl)-ortho-diiminosemiquinonate) in which each redox-active ligand has been oxidized by one electron. The oxidation products 3a-c have been structurally characterized and display an unusual parallel stacked arrangement of the disq- ligands in the solid state, with a separation of approximately 3 A. Density functional calculations show a bonding-type interaction between the SOMOs of the disq- ligands to form a unique HOMO while the antibonding linear combination forms a unique LUMO. This orbital configuration leads to a closed-shell-singlet ground-state electron configuration (S = 0). Temperature-dependent magnetism measurements indicate a low-lying triplet excited state at approximately 750 cm-1. In solution, 3a-c show strong disq--based absorption bands that are invariant across the halide series. Taken together these spectroscopic measurements provide experimental values for the one- and two-electron energies that characterize the pi-stacked bonding interaction between the two disq- ligands.

Crystal packing modifies ligand binding affinity: the case of aldose reductase.

The relationship between the structures of protein-ligand complexes existing in the crystal and in solution, essential in the case of fragment-based screening by X-ray crystallography (FBS-X), has been often an object of controversy. To address this question, simultaneous co-crystallization and soaking of two inhibitors with different ratios, Fidarestat (FID; K(d) = 6.5 nM) and IDD594 (594; K(d) = 61 nM), which bind to h-aldose reductase (AR), have been performed. The subatomic resolution of the crystal structures allows the differentiation of both inhibitors, even when the structures are almost superposed. We have determined the occupation ratio in solution by mass spectrometry (MS) Occ(FID)/Occ(594) = 2.7 and by X-ray crystallography Occ(FID)/Occ(594) = 0.6. The occupancies in the crystal and in solution differ 4.6 times, implying that ligand binding potency is influenced by crystal contacts. A structural analysis shows that the Loop A (residues 122-130), which is exposed to the solvent, is flexible in solution, and is involved in packing contacts within the crystal. Furthermore, inhibitor 594 contacts the base of Loop A, stabilizing it, while inhibitor FID does not. This is shown by the difference in B-factors of the Loop A between the AR-594 and AR-FID complexes. A stable loop diminishes the entropic energy barrier to binding, favoring 594 versus FID. Therefore, the effect of the crystal environment should be taken into consideration in the X-ray diffraction analysis of ligand binding to proteins. This conclusion highlights the need for additional methodologies in the case of FBS-X to validate this powerful screening technique, which is widely used.

Use of the FACTS solvation model for protein-ligand docking calculations. Application to EADock.

Protein-ligand docking has made important progress during the last decade and has become a powerful tool for drug development, opening the way to virtual high throughput screening and in silico structure-based ligand design. Despite the flattering picture that has been drawn, recent publications have shown that the docking problem is far from being solved, and that more developments are still needed to achieve high successful prediction rates and accuracy. Introducing an accurate description of the solvation effect upon binding is thought to be essential to achieve this goal. In particular, EADock uses the Generalized Born Molecular Volume 2 (GBMV2) solvent model, which has been shown to reproduce accurately the desolvation energies calculated by solving the Poisson equation. Here, the implementation of the Fast Analytical Continuum Treatment of Solvation (FACTS) as an implicit solvation model in small molecules docking calculations has been assessed using the EADock docking program. Our results strongly support the use of FACTS for docking. The success rates of EADock/FACTS and EADock/GBMV2 are similar, i.e. around 75% for local docking and 65% for blind docking. However, these results come at a much lower computational cost: FACTS is 10 times faster than GBMV2 in calculating the total electrostatic energy, and allows a speed up of EADock by a factor of 4. This study also supports the EADock development strategy relying on the CHARMM package for energy calculations, which enables straightforward implementation and testing of the latest developments in the field of Molecular Modeling.

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.

Analyse von Protein-Ligand Interaktionen zur Charakterisierung zyklonukleotidbindender Proteine

Hyperpolarisations-aktivierte zyklonukleotid-gesteuerte (HCN) Kanäle übernehmen wichtige Funktionen in der Regulation der Herz- und Neuronalaktivität und können über einen dualen Mechanismus aus Membranhyperpolarisation und der Bindung von zyklischen Nukleotiden aktiviert werden. Ein großes Ziel der aktuellen Forschung ist die Entwicklung neuartiger Inhibitoren, die einer Fehlregulation der Kanäle entgegenwirken. In der vorliegenden Arbeit wurde die Regulation von HCN Kanälen durch zyklische Nukleotide im Detail analysiert, indem erstmals ein umfassender Screen mit 48 unterschiedlichen Zyklonukleotid-Analoga am C-terminalen Bereich (bestehend aus C-Linker und Zyklonukleotid-Bindedomäne) der drei Isoformen HCN1, HCN2 und HCN4 durchgeführt wurde. Mit Hilfe eines Fluoreszenzpolarisations-Assays wurde der Einfluss von Modifikationen in der Base, der Ribose und dem zyklischen Phosphat auf die Bindungsaffinitäten innerhalb der Zyklonukleotid-Bindedomäne untersucht. Zyklonukleotid-Analoga mit Modifikationen an der Position 7 und 8 der Base verschoben die apparenten Affinitäten im Vergleich zu den beiden natürlich vorkommenden Zyklonukleotiden cAMP und cGMP vom mikromolaren in den nanomolaren Bereich. Selektiv für die HCN4 Isoform erwiesen sich Zyklonukleotid-Analoga mit Modifikationen an der Position 6 der Base, während Modifikationen an der Position 8 der Base zu einer höheren Affinität für die HCN2 Isoform führten. Im Gegensatz zu HCN2 und HCN4 zeigte die HCN1 Isoform besonders hohe Affinitäten für Zyklonukleotid-Analoga mit Modifikationen an der Position 8 von cGMP. Eine Substitution der 2’-Hydroxylgruppe erlaubte keine Bindung an die HCN Kanäle. Mit 7-CH-cAMP konnte ein hochaffines Bindemolekül für HCN Kanäle identifiziert werden, denn der Austausch eines Stickstoffs gegen eine CH-Einheit an Position 7 der Base führte zu einer 100-fachen Steigerung der Affinität im Vergleich zu cAMP. In Übereinstimmung mit der hochaffinen Bindung konnte in kinetischen Analysen eine langsamere Dissoziationsrate für 7-CH-cAMP gemessen werden. Anhand thermodynamischer Messungen konnte ein entropisch favorisierter Bindungsmodus für 7-CH-cAMP im Vergleich zu cAMP identifiziert werden. Basierend auf einer Kristallstruktur des HCN4 CNBD:7-CH-cAMP Komplexes (2,5 Å) lässt sich erklären, dass 7-CH-cAMP durch seine höhere Lipohilie im Vergleich zu cAMP eine stärkere Präferenz für das hydrophobe Netzwerk zwischen Protein und Base besitzt. In detaillierten, vergleichenden Analysen mit den zyklonukleotidbindenden Proteinen PKA Typ I und II, hPKGIβ und Epac 1 und 2 konnte gezeigt werden, dass 7-CH-cAMP die höchsten Affinitäten für die drei Isoformen der HCN Kanäle aufwies. Somit könnte sich 7-CH-cAMP als vielversprechender Kandidat für die selektive Regulation von HCN Kanälen in vitro und in lebenden Zellen eignen und möglicherweise einen wichtigen Beitrag als krankheitsrelevanter Effektor leisten.