Side chain to main chain hydrogen bonds stabilize a polyglutamine helix in a transcription factor.

Polyglutamine (polyQ) tracts are regions of low sequence complexity frequently found in transcription factors. Tract length often correlates with transcriptional activity and expansion beyond specific thresholds in certain human proteins is the cause of polyQ disorders. To study the structural basis of the association between tract length, transcriptional activity and disease, we addressed how the conformation of the polyQ tract of the androgen receptor, associated with spinobulbar muscular atrophy (SBMA), depends on its length. Here we report that this sequence folds into a helical structure stabilized by unconventional hydrogen bonds between glutamine side chains and main chain carbonyl groups, and that its helicity directly correlates with tract length. These unusual hydrogen bonds are bifurcate with the conventional hydrogen bonds stabilizing α-helices. Our findings suggest a plausible rationale for the association between polyQ tract length and androgen receptor transcriptional activity and have implications for establishing the mechanistic basis of SBMA.

Targeting RNA structure in SMN2 reverses spinal muscular atrophy molecular phenotypes.

Modification of SMN2 exon 7 (E7) splicing is a validated therapeutic strategy against spinal muscular atrophy (SMA). However, a target-based approach to identify small-molecule E7 splicing modifiers has not been attempted, which could reveal novel therapies with improved mechanistic insight. Here, we chose as a target the stem-loop RNA structure TSL2, which overlaps with the 5' splicing site of E7. A small-molecule TSL2-binding compound, homocarbonyltopsentin (PK4C9), was identified that increases E7 splicing to therapeutic levels and rescues downstream molecular alterations in SMA cells. High-resolution NMR combined with molecular modelling revealed that PK4C9 binds to pentaloop conformations of TSL2 and promotes a shift to triloop conformations that display enhanced E7 splicing. Collectively, our study validates TSL2 as a target for small-molecule drug discovery in SMA, identifies a novel mechanism of action for an E7 splicing modifier, and sets a precedent for other splicing-mediated diseases where RNA structure could be similarly targeted.

Computational Prediction of HIV-1 Resistance to Protease Inhibitors.

The development of mutations in HIV-1 protease (PR) hinders the activity of antiretroviral drugs, forcing changes in drug prescription. Most resistance assessments used to date rely on expert-based rules on predefined sets of stereotypical mutations; such an information-driven approach cannot capture new polymorphisms or be applied for new drugs. Computational modeling could provide a more general assessment of drug resistance and could be made available to clinicians through the Internet. We have created a protocol involving sequence comparison and all-atom protein-ligand induced fit simulations to predict resistance at the molecular level. We first compared our predictions with the experimentally determined IC50 values of darunavir, amprenavir, ritonavir, and indinavir from reference PR mutants displaying different resistance levels. We then performed analyses on a large set of variants harboring more than 10 mutations. Finally, several sequences from real patients were analyzed for amprenavir and darunavir. Our computational approach detected all of the genotype changes triggering high-level resistance, even those involving a large number of mutations.

Ab initio study of the structural, tautomeric, pairing, and electronic properties of seleno-derivatives of thymine.

The structural, tautomeric, hydrogen-bonding, stacking, and electronic properties of a seleno-derivative of thymine (T), denoted here as 4SeT and created by replacing O4 in T with Se, are investigated by means of ab initio computational techniques. The structural properties of T and 4SeT are very similar, and the geometrical differences are mainly limited to the adjacent environment of the C-Se bond. The canonical "keto" form is the most stable tautomer, in the gas phase and in aqueous solution, for both T and 4SeT. It is argued that the competition between two opposite trends, i.e., a decrease in the base-pairing ability and an increase of the stacking interaction upon incorporation of 4SeT into a duplex, likely explains the similar experimental melting points of a seleno-derivative duplex (Se-DNA) and its native counterpart. Interestingly, the underlying electronic structure shows that replacement of O4 with Se promotes a reduction in the HOMO-LUMO gap and an increase in interplane coupling, which suggests that Se-DNA could be potentially useful for nanodevice applications. This finding is further supported by the fact that transfer integrals between 4SeT...A stacked base pairs are larger than those determined for similarly stacked natural T...A pairs.

COCO: a simple tool to enrich the representation of conformational variability in NMR structures.

NMR structures are typically deposited in databases such as the PDB in the form of an ensemble of structures. Generally, each of the models in such an ensemble satisfies the experimental data and is equally valid. No unique solution can be calculated because the experimental NMR data is insufficient, in part because it reflects the conformational variability and dynamical behavior of the molecule in solution. Even for relatively rigid molecules, the limited number of structures that are typically deposited cannot completely encompass the structural diversity allowed by the observed NMR data, but they can be chosen to try and maximize its representation. We describe here the adaptation and application of techniques more commonly used to examine large ensembles from molecular dynamics simulations, to the analysis of NMR ensembles. The approach, which is based on principal component analysis, we call COCO ("Complementary Coordinates"). The COCO approach analyses the distribution of an NMR ensemble in conformational space, and generates a new ensemble that fills "gaps" in the distribution. The method is very rapid, and analysis of a 25-member ensemble and generation of a new 25 member ensemble typically takes 1-2 min on a conventional workstation. Applied to the 545 structures in the RECOORD database, we find that COCO generates new ensembles that are as structurally diverse-both from each other and from the original ensemble-as are the structures within the original ensemble. The COCO approach does not explicitly take into account the NMR restraint data, yet in tests on selected structures from the RECOORD database, the COCO ensembles are frequently good matches to this data, and certainly are structures that can be rapidly refined against the restraints to yield high-quality, novel solutions. COCO should therefore be a useful aid in NMR structure refinement and in other situations where a richer representation of conformational variability is desired-for example in docking studies. COCO is freely accessible via the website www.ccpb.ac.uk/COCO.

The structure of human 4F2hc ectodomain provides a model for homodimerization and electrostatic interaction with plasma membrane.

4F2hc (CD98hc) is a multifunctional type II membrane glycoprotein involved in amino acid transport and cell fusion, adhesion, and transformation. The structure of the ectodomain of human 4F2hc has been solved using monoclinic (Protein Data Bank code 2DH2) and orthorhombic (Protein Data Bank code 2DH3) crystal forms at 2.1 and 2.8 A, respectively. It is composed of a (betaalpha)(8) barrel and an antiparallel beta(8) sandwich related to bacterial alpha-glycosidases, although lacking key catalytic residues and consequently catalytic activity. 2DH3 is a dimer with Zn(2+) coordination at the interface. Human 4F2hc expressed in several cell types resulted in cell surface and Cys(109) disulfide bridge-linked homodimers with major architectural features of the crystal dimer, as demonstrated by cross-linking experiments. 4F2hc has no significant hydrophobic patches at the surface. Monomer and homodimer have a polarized charged surface. The N terminus of the solved structure, including the position of Cys(109) residue located four residues apart from the transmembrane domain, is adjacent to the positive face of the ectodomain. This location of the N terminus and the Cys(109)-intervening disulfide bridge imposes space restrictions sufficient to support a model for electrostatic interaction of the 4F2hc ectodomain with membrane phospholipids. These results provide the first crystal structure of heteromeric amino acid transporters and suggest a dynamic interaction of the 4F2hc ectodomain with the plasma membrane.

Exploring the binding mode of semicarbazide-sensitive amine oxidase/VAP-1: identification of novel substrates with insulin-like activity.

We previously reported that substrates of semicarbazide-sensitive amine oxidase in combination with low concentrations of vanadate exert potent insulin-like effects. Here we performed homology modeling of the catalytic domain of mouse SSAO/VAP-1 and searched through chemical databases to identify novel SSAO substrates. The modeling of the catalytic domain revealed that aromatic residues Tyr384, Phe389, and Tyr394 define a pocket of stable size that may participate in the binding of apolar substrates. We identified a number of amines as substrates of human, rat, and mouse SSAO. The compounds PD0119035, 2,3-dimethoxy-benzylamine, and C-naphthalen-1-yl-methylamine showed high affinity as substrates of rat SSAO. C-Naphthalen-1-yl-methylamine was the only substrate that showed high affinity for human SSAO. C-Naphthalen-1-yl-methylamine and 4-aminomethyl-benzenesulfonamide showed the highest capacity to stimulate glucose transport in isolated rat adipocytes. The impact of these findings on the development of new treatments for diabetes is discussed.

Unconventional interactions between water and heterocyclic nitrogens in protein structures.

We report an unusual interaction in which a water molecule approaches the heterocyclic nitrogen of tryptophan and histidine along an axis that is roughly perpendicular to the aromatic plane of the side chain. The interaction is distinct from the well-known conventional aromatic hydrogen-bond, and it occurs at roughly the same frequency in protein structures. Calculations indicate that the water-indole interaction is favorable energetically, and we find several cases in which such contacts are conserved among structural orthologs. The indole-water interaction links side chains and peptide backbone in turn regions, connects the side chains in beta-sheets, and bridges secondary elements from different domains. We suggest that the water-indole interaction can be indirectly responsible for the quenching of tryptophan fluorescence that is observed in the folding of homeodomains and, possibly, many other proteins. We also observe a similar interaction between water and the imidazole nitrogens of the histidine side chain. Taken together, these observations suggest that the unconventional water-indole and water-imidazole interactions provide a small but favorable contribution to protein structures.