A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding.

We propose a general mean field model of ligand-protein interactions to determine the thermodynamic equilibrium of a system at finite temperature. The method is employed in structural assessments of two human immuno-deficiency virus type 1 protease complexes where the gross effects of protein flexibility are incorporated by utilizing a data base of crystal structures. Analysis of the energy spectra for these complexes has revealed that structural and thermo-dynamic aspects of molecular recognition can be rationalized on the basis of the extent of frustration in the binding energy landscape. In particular, the relationship between receptor-specific binding of these ligands to human immunodeficiency virus type 1 protease and a minimal frustration principle is analyzed.

Design of 1D and 2D molecule-based magnets with the ligand 4,5-dimethyl-1,2-phenylenebis(oxamato)

Three new bimetallic oxamato-based magnets with the proligand 4,5-dimethyl-1,2-phenylenebis-(oxamato) (dmopba) were synthesized using water or dimethylsulfoxide (DMSO) as solvents. Single crystal X-ray diffraction provided structures for two of them: [MnCu(dmopba)(H(2)O)(3)]n center dot 4nH(2)O (1) and [MnCu(dmopba)(DMSO)(3)](n center dot)nDMSO (2). The crystalline structures for both 1 and 2 consist of linearly ordered oxamato-bridged Mn(II)Cu(II) bimetallic chains. The magnetic characterization revealed a typical behaviour of ferrimagnetic chains for 1 and 2. Least-squares fits of the experimental magnetic data performed in the 300-20 K temperature range led to J(MnCu) = -27.9 cm(-1), g(Cu) = 2.09 and g(Mn) = 1.98 for 1 and J(MnCu) = -30.5 cm(-1), g(Cu) = 2.09 and g(Mn) = 2.02 for 2 (H = -J(MnCu)Sigma S(Mn, i)(S(Cu, i) + S(Cu, i-1))). The two-dimensional ferrimagnetic system [Me(4)N](2n){Co(2)[Cu(dmopba)](3)}center dot 4nDMSO center dot nH(2)O (3) was prepared by reaction of Co(II) ions and an excess of [Cu(dmopba)](2-) in DMSO. The study of the temperature dependence of the magnetic susceptibility as well as the temperature and field dependences of the magnetization revealed a cluster glass-like behaviour for 3.

SwissParam: a fast force field generation tool for small organic molecules.

The drug discovery process has been deeply transformed recently by the use of computational ligand-based or structure-based methods, helping the lead compounds identification and optimization, and finally the delivery of new drug candidates more quickly and at lower cost. Structure-based computational methods for drug discovery mainly involve ligand-protein docking and rapid binding free energy estimation, both of which require force field parameterization for many drug candidates. Here, we present a fast force field generation tool, called SwissParam, able to generate, for arbitrary small organic molecule, topologies, and parameters based on the Merck molecular force field, but in a functional form that is compatible with the CHARMM force field. Output files can be used with CHARMM or GROMACS. The topologies and parameters generated by SwissParam are used by the docking software EADock2 and EADock DSS to describe the small molecules to be docked, whereas the protein is described by the CHARMM force field, and allow them to reach success rates ranging from 56 to 78%. We have also developed a rapid binding free energy estimation approach, using SwissParam for ligands and CHARMM22/27 for proteins, which requires only a short minimization to reproduce the experimental binding free energy of 214 ligand-protein complexes involving 62 different proteins, with a standard error of 2.0 kcal mol(-1), and a correlation coefficient of 0.74. Together, these results demonstrate the relevance of using SwissParam topologies and parameters to describe small organic molecules in computer-aided drug design applications, together with a CHARMM22/27 description of the target protein. SwissParam is available free of charge for academic users at www.swissparam.ch.

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.

Recent advances in the epidermal growth factor receptor/ligand system biology on skin homeostasis and keratinocyte stem cell regulation.

The epidermal growth factor (EGF) receptor/ligand system stimulates multiple pathways of signal transduction, and is activated by various extracellular stimuli and inter-receptor crosstalk signaling. Aberrant activation of EGF receptor (EGFR) signaling is found in many tumor cells, and humanized neutralizing antibodies and synthetic small compounds against EGFR are in clinical use today. However, these drugs are known to cause a variety of skin toxicities such as inflammatory rash, skin dryness, and hair abnormalities. These side effects demonstrate the multiple EGFR-dependent homeostatic functions in human skin. The epidermis and hair follicles are self-renewing tissues, and keratinocyte stem cells are crucial for maintaining these homeostasis. A variety of molecules associated with the EGF receptor/ligand system are involved in epidermal homeostasis and hair follicle development, and the modulation of EGFR signaling impacts the behavior of keratinocyte stem cells. Understanding the roles of the EGF receptor/ligand system in skin homeostasis is an emerging issue in dermatology to improve the current therapy for skin disorders, and the EGFR inhibitor-associated skin toxicities. Besides, controlling of keratinocyte stem cells by modulating the EGF receptor/ligand system assures advances in regenerative medicine of the skin. We present an overview of the recent progress in the field of the EGF receptor/ligand system on skin homeostasis and regulation of keratinocyte stem cells.

Unconventional integrin-ligand interactions

Integrins are cell surface adhesion and signaling receptors. Cells use integrins to attach to the extracellular matrix and to other cells, as well as for sensing their environment. In addition to adhesion and migration, integrins have been shown to be important for many biological processes including apoptosis, cell proliferation, and differentiation into specific tissues. Many important next generation biological drugs inhibit integrin functions. Thus, research into interactions between integrins and their ligands under different physiological and pathological conditions is not only of academic interest, but is also important for the field of drug discovery. In this Ph.D. project, the functions of integrin-ligand interactions were studied under different physiologically interesting conditions including 1) human echovirus 1 binding to integrin α2β1, 2) integrin α2β1 binding to collagen under flow conditions, 3) integrin α2β1 binding to a ligand in the presence of the angiogenesis inhibitor histidine rich glycoprotein (HRG) and 4) integrin binding to posttranslationally citrullinated ligands. As a result of the project, we could show that for each condition the integrin-ligand interaction is somewhat unconventional. 1) Echovirus 1 binds only to non-activated conformations of integrin α2β1. 2) Surprisingly, the non-activated conformation is also the primary conformation of integrin α2β1 when it binds to collagen under flow conditions, like when platelets adhere to subendothelial collagen in vascular injuries. In addition, the pre-activation of integrin α2β1 does not increase adhesion under flow. 3) HRG binds to integrin α2β1 through a low-affinity interaction that inhibits integrin binding to collagen. This shows that low affinity interactions could be biologically relevant and possibly regulate angiogenesis. 4) The citrullination of collagen, a posttranslational modification reported to occur in rheumatoid arthritis, specifically inhibits the binding of integrin α10β1 and α11β1, but does not affect the binding of α1β1 ja α2β1. On the other hand, the citrullination of isoDGR in fibronectin and RGD in pro-TGF- β:n inhibit integrin binding completely. Citrullination seems to be an inflammation related process and integrin ligands become citrullinated frequently in vivo. This Ph.D. thesis suggests that unconventional interaction mechanisms between integrins and their ligands, such as posttranslational modifications, low affinity interactions, and non-activated integrin conformations, can have an important role in pathological processes. The study of these kinds of integrin-ligand interactions is important for understanding biological phenomena more deeply. The research might also be beneficial for the development of integrin based therapies for treating diseases.

Ligand Design for Dual Property Single Molecule Magnets

This thesis describes two different approaches for the preparation of polynuclear clusters with interesting structural, magnetic and optical properties. Firstly, exploiting p-tert-butylcalix[4]arene (TBC4) macrocycles together with selected Ln(III) ions for the assembly of emissive single molecule magnets, and secondly the preparation and coordination of a chiral mpmH ligand with selected 3d transition metal ions, working towards the discovery of chiral polynuclear clusters. In Project 1, the coordination chemistry of the TBC4 macrocycle together with Dy(III) and Tb(III) afforded two Ln6[TBC4]2 complexes that have been structurally, magnetically and optically characterized. X-ray diffraction studies reveal that both complexes contain an octahedral core of Ln6 ions capped by two fully deprotonated TBC4 macrocycles. Although the unit cells of the two complexes are very similar, the coordination geometries of their Ln(III) ions are subtly different. Variable temperature ac magnetic susceptibility studies reveal that both complexes display single molecule magnet (SMM) behaviour in zero dc field and the energy barriers and associated pre-exponential factors for each relaxation process have been determined. Low temperature solid state photoluminescence studies reveal that both complexes are emissive; however, the f-f transitions within the Dy6 complex were masked by broad emissions from the TBC4 ligand. In contrast, the Tb(III) complex displayed green emission with the spectrum comprising four sharp bands corresponding to 5D4 → 7FJ transitions (where J = 3, 4, 5 and 6), highlighting that energy transfer from the TBC4 macrocycle to the Tb(III) ion is more effective than to Dy. Examples of zero field Tb(III) SMMs are scarce in the chemical literature and the Tb6[TBC4]2 complex represents the first example of a Tb(III) dual property SMM assembled from a p-tert-butylcalix[4]arene macrocycle with two magnetically derived energy barriers, Ueff of 79 and 63 K. In Project 2, the coordination of both enantiomers of the chiral ligand, α-methyl-2-pyridinemethanol (mpmH) to Ni(II) and Co(II) afforded three polynuclear clusters that have been structurally and magnetically characterized. The first complex, a Ni4 cluster of stoichiometry [Ni4(O2CCMe3)4(mpm)4]·H2O crystallizes in a distorted cubane topology that is well known in Ni(II) cluster chemistry. The final two Co(II) complexes crystallize as a linear mixed valence trimer with stoichiometry [Co3(mpm)6]·(ClO4)2, and a Co4 mixed valence complex [Co(II)¬2Co(III)2(NO3)2(μ-mpm)4(ONO2)2], whose structural topology resembles that of a defective double cubane. All three complexes crystallize in chiral space groups and circular dichroism experiments further confirm that the chirality of the ligand has been transferred to the respective coordination complex. Magnetic susceptibility studies reveal that for all three complexes, there are competing ferro- and antiferromagnetic exchange interactions. The [Co(II)¬2Co(III)2(NO3)2(μ-mpm)4(ONO2)2] complex represents the first example of a chiral mixed valence Co4 cluster with a defective double cubane topology.

Versatile binding properties of a di-2-pyridyl ketone nicotinoylhydrazone ligand: Crystal structure of a Cu(II) complex

Six new copper complexes of di-2-pyridyl ketone nicotinoylhydrazone (HDKN) have been synthesized. The complexes have been characterized by a variety of spectroscopic techniques and the structure of [Cu(DKN)2]·H2O has been determined by single crystal X-ray diffraction. The compound [Cu(DKN)2]·H2O crystallized in the monoclinic space group P21 and has a distorted octahedral geometry. The IR spectra revealed the presence of variable modes of chelation for the investigated ligand. The EPR spectra of compounds [Cu2(DKN)2( -N3)2] and [Cu2(DKN)2( -NCS)2] in polycrystalline state suggest a dimeric structure as they exhibited a half field signal, which indicate the presence of a weak interaction between two Cu(II) ions in these complexes

A new cyano-bridged one-dimensional (GdFeIII)-Fe-III coordination polymer with o-phenanthroline as the blocking ligand: synthesis, structure, and magnetic properties

Synthesis, structural characterization, and magnetic properties of a new cyano-bridged one-dimensional iron (III)-gadolinium (III) compound, trans-[Gd(o-phen)(2)(H2O)(2)(mu-CN)(2)Fe(CN)(4)], - 2no-phen (o-phen = 1,10-phenanthroline), have been described. The compound crystallizes in the triclinic P (1) over bar space group with the following unit cell parameters: a = 10.538(14) angstrom, b = 12.004(14) angstrom, c = 20.61(2) angstrom, alpha = 92.41(1)degrees, beta = 92.76(1)degrees, gamma = 11 2.72(1)degrees, and Z = 2. In this complex, each gadolinium (III) is coordinated to two nitrile nitrogens of the CN groups coming from two different ferricyanides, the mutually trans cyanides of each of which links another different Gd-III to create -NC-Fe(CN)(4)-CN-Gd-NC- type 1-D chain structure. The one-dimensional chains are self-assembled in two-dimensions via weak C-H center dot center dot center dot N hydrogen bonds. Both the variable-temperature (2-300 K, 0.01 T and 0.8 T) and variable-field (0-50 000 Gauss, 2 K) magnetic measurements reveal the existence of very weak interaction in this molecule. The temperature dependence of the susceptibilities has been analyzed using a model for a chain of alternating classic (7/2) and quantum (1/2) spins. (c) 2005 Elsevier B.V. All rights reserved.

Use of a reduced Schiff-Base ligand to prepare novel chloro-bridged chains of rare Cu(II) trinuclear complexes with mixed azido/oxo and chloro/oxo bridges

Two mixed bridged one-dimensional (1D) polynuclear complexes, [Cu3L2(mu(1,1)-N-3)(2)(mu-Cl)Cl](n) (1) and {[Cu3L2(mu-Cl)(3)Cl]center dot 0.46CH(3)OH}(n), (2), have been synthesized using the tridentate reduced Schiff-base ligand HL (2-[(2-dimethylamino-ethylamino)-methyl]-phenol). The complexes have been characterized by X-ray structural analyses and variable-temperature magnetic susceptibility measurements. In both complexes the basic trinuclear angular units are joined together by weak chloro bridges to form a 1D chain. The trinuclear structure of 1 is composed of two terminal square planar [Cu(L)(mu(1,1)-N-3)] units connected by a central Cu(II) atom through bridging nitrogen atoms of end-on azido ligands and the phenoxo oxygen atom of the tridentate ligand. These four coordinating atoms along with a chloride ion form a distorted trigonal bipyramidal geometry around the central Cu(II). The structure of 2 is similar; the only difference being a Cl bridge replacing the mu(1,1)-N-3 bridge in the trinuclear unit. The magnetic properties of both trinuclear complexes can be very well reproduced with a simple linear symmetrical trimer model (H = JS(i)S(i+1)) with only one intracluster exchange coupling (J) including a weak intertrimer interaction (.j) reproduced with the molecular field approximation. This model provides very satisfactory fits for both complexes in the whole temperature range with the following parameters: g = 2.136(3), J = 93.9(3) cm(-1) and zj= -0.90(3) cm(-1) (z = 2) for 1 and g = 2.073(7), J = -44.9(4) cm(-1) and zJ = -1.26(6) cm(-1) (z = 2) for 2.

Highly efficient separation of actinides from lanthanides by a phenanthroline-derived bis-triazine ligand

The synthesis, lanthanide complexation, and solvent ex- traction of actinide(III) and lanthanide(III) radiotracers from nitric acid solutions by a phenanthroline-derived quadridentate bis-triazine ligand are described. The ligand separates Am(III) and Cm(III) from the lanthanides with remarkably high efficiency, high selectivity, and fast extraction kinetics compared to its 2,2'-bipyridine counterpart. Structures of the 1:2 bis-complexes of the ligand with Eu(III) and Yb(III) were elucidated by X-ray crystallography and force field calculations, respec-tively. The Eu(III) bis-complex is the first 1:2 bis-complex of a quadridentate bis-triazine ligand to be characterized by crystallography. The faster rates of extraction were verified by kinetics measurements using the rotating membrane cell technique in several diluents. The improved kinetics of metal ion extraction are related to the higher surface activity of the ligand at the phase interface. The improvement in the ligand's properties on replacing the bipyridine unit with a phenanthroline unit far exceeds what was anticipated based on ligand design alone.

Ligand substituents effect on 10Dq

Reaction of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline (L) with Ni(ClO4)(2)center dot 6H(2)O in methanol in 3:1 M proportion at room temperature yields [NiL3](ClO4)(2)center dot 2H(2)O. The X-ray crystal structure of the cation Nil(3)(2+) has been determined. Aminolysis of the three epoxide rings in NiL32+ by 4-substituted anilines in boiling water without any Lewis acid catalyst gives a family of Ni(II) complexes with octahedral NiL62+ core. In these complexes, crystal field splitting 10Dq varies from 11601 to 15798 cm(-1) in acetonitrile. The variation in 10Dq is found to be satisfactorily linear (r(2) = 0.951) with the Hammett sigma(R) parameter of the substituent on the anilino fragment. 10Dq increases with the increase in the electron donation ability of the substituent.

Structure and magnetic properties of a tetranuclear Cu(II) complex containing the 2-(pyridine-2-yliminomethyl)-phenol ligand

A tetranuclear Cu(II) complex [Cu4L4(H2O)4](ClO4)4 has been synthesized using the terdentate Schiff base 2-(pyridine-2-yliminomethyl)-phenol (HL) (the condensation product of salicylaldehyde and 2-aminopyridine) and copper perchlorate. Chemical characterizations such as IR and UV/Vis of the complex have been carried out. A single-crystal diffraction study shows that the complex contains a nearly planar tetranuclear core containing four copper atoms, which occupy four equivalent five-coordinate sites with a square pyramidal environment. Magnetic measurements have been carried out over the temperature range 2–300K and with 100Oe field strengths. Analysis of magnetic susceptibility data indicates a strong antiferromagnetic (J1=−638cm−1) exchange interaction between diphenoxo-bridged Cu(II) centers and a moderate antiferromagnetic (J2=−34cm−1) interaction between N–C–N bridged Cu(II) centers. Magnetic exchange interactions (J’s) are also discussed on the basis of a computational study using DFT methodology. The spin density distribution (singlet ground state) is calculated to visualize the effect of delocalization of spin density through bridging groups.

In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.