Local selection rules that can determine specific pathways of DNA unknotting by type II DNA topoisomerases.

We performed numerical simulations of DNA chains to understand how local geometry of juxtaposed segments in knotted DNA molecules can guide type II DNA topoisomerases to perform very efficient relaxation of DNA knots. We investigated how the various parameters defining the geometry of inter-segmental juxtapositions at sites of inter-segmental passage reactions mediated by type II DNA topoisomerases can affect the topological consequences of these reactions. We confirmed the hypothesis that by recognizing specific geometry of juxtaposed DNA segments in knotted DNA molecules, type II DNA topoisomerases can maintain the steady-state knotting level below the topological equilibrium. In addition, we revealed that a preference for a particular geometry of juxtaposed segments as sites of strand-passage reaction enables type II DNA topoisomerases to select the most efficient pathway of relaxation of complex DNA knots. The analysis of the best selection criteria for efficient relaxation of complex knots revealed that local structures in random configurations of a given knot type statistically behave as analogous local structures in ideal geometric configurations of the corresponding knot type.

The major genetic determinants of HIV-1 control affect HLA class I peptide presentation.

Infectious and inflammatory diseases have repeatedly shown strong genetic associations within the major histocompatibility complex (MHC); however, the basis for these associations remains elusive. To define host genetic effects on the outcome of a chronic viral infection, we performed genome-wide association analysis in a multiethnic cohort of HIV-1 controllers and progressors, and we analyzed the effects of individual amino acids within the classical human leukocyte antigen (HLA) proteins. We identified >300 genome-wide significant single-nucleotide polymorphisms (SNPs) within the MHC and none elsewhere. Specific amino acids in the HLA-B peptide binding groove, as well as an independent HLA-C effect, explain the SNP associations and reconcile both protective and risk HLA alleles. These results implicate the nature of the HLA-viral peptide interaction as the major factor modulating durable control of HIV infection.

Constitutively active G protein-coupled receptors: potential mechanisms of receptor activation.

Mutations of G protein-coupled receptors can increase their constitutive (agonist-independent) activity. Some of these mutations have been artificially introduced by site-directed mutagenesis, others occur spontaneously in human diseases. The analysis of the constitutively active G protein-coupled receptors has provided important informations about the molecular mechanisms underlying receptor activation and drug action.

Stereospecificity of the siderophore pyochelin outer membrane transporters in fluorescent pseudomonads.

Pyochelin (Pch) and enantio-pyochelin (EPch) are enantiomer siderophores that are produced by Pseudomonas aeruginosa and Pseudomonas fluorescens, respectively, under iron limitation. Pch promotes growth of P. aeruginosa when iron is scarce, and EPch carries out the same biological function in P. fluorescens. However, the two siderophores are unable to promote growth in the heterologous species, indicating that siderophore-mediated iron uptake is highly stereospecific. In the present work, using binding and iron uptake assays, we found that FptA, the Fe-Pch outer membrane transporter of P. aeruginosa, recognized (K(d) = 2.5 +/- 1.1 nm) and transported Fe-Pch but did not interact with Fe-EPch. Likewise, FetA, the Fe-EPch receptor of P. fluorescens, was specific for Fe-EPch (K(d) = 3.7 +/- 2.1 nm) but did not bind and transport Fe-Pch. Growth promotion experiments performed under iron-limiting conditions confirmed that FptA and FetA are highly specific for Pch and EPch, respectively. When fptA and fetA along with adjacent transport genes involved in siderophore uptake were swapped between the two bacterial species, P. aeruginosa became able to utilize Fe-EPch as an iron source, and P. fluorescens was able to grow with Fe-Pch. Docking experiments using the FptA structure and binding assays showed that the stereospecificity of Pch recognition by FptA was mostly due to the configuration of the siderophore chiral centers C4'' and C2'' and was only weakly dependent on the configuration of the C4' carbon atom. Together, these findings increase our understanding of the stereospecific interaction between Pch and its outer membrane receptor FptA.

How do thyroid hormone receptors bind to structurally diverse response elements?

Three classes of thyroid hormone response elements have been described. They are composed of two half-sites arranged either as a palindromic, a direct repeat or as an inverted palindromic array. Receptor homodimers as well as heterodimers can bind to all three types of response element. While the ligand binding domain of the receptors provides the major dimerization surface, asymmetric contacts between the DNA binding domains are necessary for binding to a direct repeat. Moreover, some recent findings suggest that in TR, compared to RXR, the ligand binding domain has a 180 degrees rotation with respect to the DNA binding domain. This feature could explain the preferential binding of the RXR-TR heterodimer to the direct repeat response element, in which RXR exclusively binds the 5' half-site, and of the TR homodimer to the inverted palindrome response element.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

A template-assembled synthetic protein surface mimetic of the von Willebrand factor A1 domain inhibits botrocetin-induced platelet aggregation.

Platelet adhesion, the initial step of platelet activation, is mediated by the interaction of von Willebrand factor (VWF) with its platelet receptor, the GPIb-IX complex. The binding of VWF to GPIb-IX is induced either by increased shear stress or by exogenous modulators, such as botrocetin. At a molecular level, this interaction takes place between the A1 domain of VWF and the GPIb alpha chain of the GPIb-IX complex. We report here the design and functional characteristics of a VWF template-assembled synthetic protein (TASP), a chimeric four-helix-bundle TASP scaffold mimicking the surface of the A1 domain. Twelve residues located on helices alpha 3 and alpha 4 in the native A1 domain were grafted onto a surface formed by two neighboring helices of the TASP. VWF TASP was found to inhibit specifically botrocetin-induced platelet aggregation and to bind both botrocetin and GPIb alpha. However, in contrast to the native A1 domain, VWF TASP did not bind simultaneously to both ligands. Modeling studies revealed that the relative orientation of the alpha helices in VWF TASP led to a clash of bound botrocetin and GPIb alpha. These results demonstrate that a chimeric four-helix-bundle TASP as a scaffold offers a suitable surface for presenting crucial residues of the VWF A1 domain; the potential of the TASP approach for de novo protein design and mimicry is thereby illustrated.

BAFF, APRIL and their receptors: structure, function and signaling.

BAFF, APRIL and their receptors play important immunological roles, especially in the B cell arm of the immune system. A number of splice isoforms have been described for both ligands and receptors in this subfamily, some of which are conserved between mouse and human, while others are species-specific. Structural and mutational analyses have revealed key determinants of receptor-ligand specificity. BAFF-R has a strong selectivity for BAFF; BCMA has a higher affinity for APRIL than for BAFF, while TACI binds both ligands equally well. The molecular signaling events downstream of BAFF-R, BCMA and TACI are still incompletely characterized. Survival appears to be mediated by upregulation of Bcl-2 family members through NF-kappaB activation, degradation of the pro-apototic Bim protein, and control of subcellular localization of PCKdelta. Very little is known about other signaling events associated with receptor engagement by BAFF and APRIL that lead for example to B cell activation or to CD40L-independent Ig switch.

DNA-DNA interactions in bacteriophage capsids are responsible for the observed DNA knotting.

Recent experiments showed that the linear double-stranded DNA in bacteriophage capsids is both highly knotted and neatly structured. What is the physical basis of this organization? Here we show evidence from stochastic simulation techniques that suggests that a key element is the tendency of contacting DNA strands to order, as in cholesteric liquid crystals. This interaction favors their preferential juxtaposition at a small twist angle, thus promoting an approximately nematic (and apolar) local order. The ordering effect dramatically impacts the geometry and topology of DNA inside phages. Accounting for this local potential allows us to reproduce the main experimental data on DNA organization in phages, including the cryo-EM observations and detailed features of the spectrum of DNA knots formed inside viral capsids. The DNA knots we observe are strongly delocalized and, intriguingly, this is shown not to interfere with genome ejection out of the phage.

Crystal structure of yeast peroxisomal multifunctional enzyme: structural basis for substrate specificity of (3R)-hydroxyacyl-CoA dehydrogenase units.

(3R)-hydroxyacyl-CoA dehydrogenase is part of multifunctional enzyme type 2 (MFE-2) of peroxisomal fatty acid beta-oxidation. The MFE-2 protein from yeasts contains in the same polypeptide chain two dehydrogenases (A and B), which possess difference in substrate specificity. The crystal structure of Candida tropicalis (3R)-hydroxyacyl-CoA dehydrogenase AB heterodimer, consisting of dehydrogenase A and B, determined at the resolution of 2.2A, shows overall similarity with the prototypic counterpart from rat, but also important differences that explain the substrate specificity differences observed. Docking studies suggest that dehydrogenase A binds the hydrophobic fatty acyl chain of a medium-chain-length ((3R)-OH-C10) substrate as bent into the binding pocket, whereas the short-chain substrates are dislocated by two mechanisms: (i) a short-chain-length 3-hydroxyacyl group ((3R)-OH-C4) does not reach the hydrophobic contacts needed for anchoring the substrate into the active site; and (ii) Leu44 in the loop above the NAD(+) cofactor attracts short-chain-length substrates away from the active site. Dehydrogenase B, which can use a (3R)-OH-C4 substrate, has a more shallow binding pocket and the substrate is correctly placed for catalysis. Based on the current structure, and together with the structure of the 2-enoyl-CoA hydratase 2 unit of yeast MFE-2 it becomes obvious that in yeast and mammalian MFE-2s, despite basically identical functional domains, the assembly of these domains into a mature, dimeric multifunctional enzyme is very different.

The implications of alternative splicing in the ENCODE protein complement.

Alternative premessenger RNA splicing enables genes to generate more than one gene product. Splicing events that occur within protein coding regions have the potential to alter the biological function of the expressed protein and even to create new protein functions. Alternative splicing has been suggested as one explanation for the discrepancy between the number of human genes and functional complexity. Here, we carry out a detailed study of the alternatively spliced gene products annotated in the ENCODE pilot project. We find that alternative splicing in human genes is more frequent than has commonly been suggested, and we demonstrate that many of the potential alternative gene products will have markedly different structure and function from their constitutively spliced counterparts. For the vast majority of these alternative isoforms, little evidence exists to suggest they have a role as functional proteins, and it seems unlikely that the spectrum of conventional enzymatic or structural functions can be substantially extended through alternative splicing.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling.

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

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.

Rational design of indoleamine 2,3-dioxygenase inhibitors.

Indoleamine 2,3-dioxygenase (IDO) is an important therapeutic target for the treatment of diseases such as cancer that involve pathological immune escape. We have used the evolutionary docking algorithm EADock to design new inhibitors of this enzyme. First, we investigated the modes of binding of all known IDO inhibitors. On the basis of the observed docked conformations, we developed a pharmacophore model, which was then used to devise new compounds to be tested for IDO inhibition. We also used a fragment-based approach to design and to optimize small organic molecule inhibitors. Both approaches yielded several new low-molecular weight inhibitor scaffolds, the most active being of nanomolar potency in an enzymatic assay. Cellular assays confirmed the potential biological relevance of four different scaffolds.

Simulations of action of DNA topoisomerases to investigate boundaries and shapes of spaces of knots.

The configuration space available to randomly cyclized polymers is divided into subspaces accessible to individual knot types. A phantom chain utilized in numerical simulations of polymers can explore all subspaces, whereas a real closed chain forming a figure-of-eight knot, for example, is confined to a subspace corresponding to this knot type only. One can conceptually compare the assembly of configuration spaces of various knot types to a complex foam where individual cells delimit the configuration space available to a given knot type. Neighboring cells in the foam harbor knots that can be converted into each other by just one intersegmental passage. Such a segment-segment passage occurring at the level of knotted configurations corresponds to a passage through the interface between neighboring cells in the foamy knot space. Using a DNA topoisomerase-inspired simulation approach we characterize here the effective interface area between neighboring knot spaces as well as the surface-to-volume ratio of individual knot spaces. These results provide a reference system required for better understanding mechanisms of action of various DNA topoisomerases.