Search

Leucine-zipper transcription regulator 1 (LZTR1) is a highly mutated tumor suppressor gene, involved in the pathogenesis of several cancer types and developmental disorders. In proteasomal degradation, it acts as an adaptor protein responsible for the recognition and recruitment of substrates to be ubiquitinated in Cullin3-RING ligase E3 (CRL3) machinery. LZTR1 belongs to the BTB-Kelch family, a multi-domain protein where the Kelch propeller plays as the substrate recognition region and for which no experimental structure has been solved. Recently, large effort mutational analyses pointed to the role of disease-associated LZTR1 mutations in the RAS/MAPK signaling pathway and RIT1, a small Ras-related GTPase protein, has been identified by mass spectroscopy to interact with LZTR1. Hence, a better understanding of native structure, molecular mechanism, and substrate specificity would help clarifying the role of LZTR1 in pathological diseases, thus promoting advancement in the development of novel therapeutic strategies. Here, we address the interaction model between adaptor LZTR1 and substrate RIT1 by applying an integrated computational approach, including molecular modeling and docking techniques. We observe that the interaction model LZTR1-RIT1 is stabilized by an electrostatic bond network established between the two protein surfaces, which is reminiscent of homologous ubiquitin ligases complexes. Then, running MD simulations, we characterize differential conformational dynamics of the multi-domain LZTR1, offering interesting implications on the mechanistic role of specific point mutations. We identify G248R and R283Q as damaging mutations involved in the recognition process of the substrate RIT1 and R412C as a possible allosteric mutation from the Kelch to the C-term BTB-domain. Our findings provide important structural insights on targeting CRL3s for drug discovery.

A strict interplay is known to involve copper and zinc in many cellular processes. For this reason, the results of copper’s interaction with zinc binding proteins are of great interest. For instance, copper interferences with the DNA-binding activity of zinc finger proteins are associated with the development of a variety of diseases. The biological impact of copper depends on the chemical properties of its two common oxidation states (Cu(I) and Cu(II)). In this framework, following the attention addressed to unveil the effect of metal ion replacement in zinc fingers and in zinc-containing proteins, we explore the effects of the Zn(II) to Cu(I) or Cu(II) replacement in the prokaryotic zinc finger domain. The prokaryotic zinc finger protein Ros, involved in the horizontal transfer of genes from A. tumefaciens to a host plant infected by it, belongs to a family of proteins, namely Ros/MucR, whose members have been recognized in different bacteria symbionts and pathogens of mammals and plants. Interestingly, the amino acids of the coordination sphere are poorly conserved in most of these proteins, although their sequence identity can be very high. In fact, some members of this family of proteins do not bind zinc or any other metal, but assume a 3D structure similar to that of Ros with the residues replacing the zinc ligands, forming a network of hydrogen bonds and hydrophobic interactions that surrogates the Zn-coordinating role. These peculiar features of the Ros ZF domain prompted us to study the metal ion replacement with ions that have different electronic configuration and ionic radius. The protein was intensely studied as a perfectly suited model of a metal-binding protein to study the effects of the metal ion replacement; it appeared to tolerate the Zn to Cd substitution, but not the replacement of the wildtype metal by Ni(II), Pb(II) and Hg(II). The structural characterization reported here gives a high-resolution description of the interaction of copper with Ros, demonstrating that copper, in both oxidation states, binds the protein, but the replacement does not give rise to a functional domain.

The control of expression in genetic regulation is a fundamental process for cell life. In RNA-mediated silencing, human Argonaute-2 protein (hAgo2) uses sequence information encoded in small RNAs (guide) to identify complementary sites in messenger RNAs (target) for repression. The specificity of this molecular recognition lies at the basis of the mechanisms that control the expression of thousands of genes, which necessarily requires a fine tuning of complex events. Among these, the binding of the first nucleotide of the target RNA (t1) is emerging as an important modulator of hAgo2-mediated machinery. Using atomistic molecular dynamics-derived analyses, we address the mechanism behind t1-dependent regulation and study the impact of different t1 nucleotides (t1A, t1C, t1G, t1U) on the conformational dynamics of both hAgo2 and guide-target RNAs. Only when an adenine is found at this position, t1 directly interacts with a specific hAgo2 binding pocket, favoring the stabilization of target binding. Our findings show that hAgo2 exploits a dynamic recognition mechanism of the t1-target thanks to a modulation of RNA conformations. Here, t1-adenine is the only nucleobase endowed with a dual binding mode: a T-shape and a co-planar conformation, respectively, orthogonal and parallel to the following base-pairs of guide-target duplex. This triggers a composite set of molecular interactions that stabilizes distinctive conformational ensembles. Our comparative analyses show characteristic traits of local and global dynamic interplay between hAgo2 and the RNA molecules and highlight how t1A binding acts as a molecular switch for target recognition and complex stabilization. Implications for future mechanistic studies are discussed.

Ligand-based control of protein functional motions can provide novel opportunities in the study of fundamental biological mechanisms and in the development of novel therapeutics. In this work we addressed the ligand-based modulation of integrin functions. Inhibitors of integrin α v β 3 are interesting anticancer agents but their molecular mechanisms are still unclear: Peptides and peptidomimetics characterized by the Arg-Gly-Asp (RGD) or isoAsp-Gly-Arg (isoDGR) binding motifs have shown controversial agonist/antagonist effects. We have investigated the differential mechanisms of integrin activation/deactivation by three distinct ligands (cyclo-RGDf(NMe)V (Cilengitide), cyclo[DKP3-RGD], cyclo[DKP3-isoDGR]; DKP=diketopiperazine) through a comparative analysis of ligand-controlled protein internal dynamics: Although RGD facilitates the onset of dynamic states leading to activation, isoDGR induces a diffuse rigidification of the complex consistent with antagonist activities. Computational predictions have been experimentally probed by showing that the antibody AP5, which is capable of recognizing the active form of integrin, binds specifically to the RGD complexes and not to the isoDGR complex, which supports opposite functional roles of the two motifs targeting the same binding site.

The crucial role of integrin in pathological processes such as tumor progression and metastasis formation has inspired intense efforts to design novel pharmaceutical agents modulating integrin functions in order to provide new tools for potential therapies. In the past decade, we have investigated the biological proprieties of the chimeric peptide RGDechi, containing a cyclic RGD motif linked to an echistatin C-terminal fragment, able to specifically recognize αvβ3 without cross reacting with αvβ5 and αIIbβ3 integrin. Additionally, we have demonstrated using two RGDechi-derived peptides, called RGDechi1-14 and ψRGDechi, that chemical modifications introduced in the C-terminal part of the peptide alter or abolish the binding to the αvβ3 integrin. Here, to shed light on the structural and dynamical determinants involved in the integrin recognition mechanism, we investigate the effects of the chemical modifications by exploring the conformational space sampled by RGDechi1-14 and ψRGDechi using an integrated natural-abundance NMR/MD approach. Our data demonstrate that the flexibility of the RGD-containing cycle is driven by the echistatin C-terminal region of the RGDechi peptide through a coupling mechanism between the N- and C-terminal regions.

Schistosomiasis is the most significant neglected tropical parasitic disease caused by helminths in terms of morbidity and mortality caused by helminths. In this work, we present the antischistosomal activity against Schistosoma mansoni of a rationally selected small set of thiazinoquinone derivatives, some of which were previously found to be active against Plasmodium falciparum and others synthesized ad hoc. The effects on larvae, juvenile, and adult parasite viability as well as on egg production and development were investigated, resulting in the identification of new multistage antischistosomal hit compounds. The most promising compounds 6, 8, 13, and 14 with a LC50 value on schistosomula from ∼5 to ∼15 μM also induced complete death of juvenile (28 days old) and adult worm pairs (7 weeks old) and a detrimental effect on egg production and development in vitro. Structure-activity relationships (SARs) were analyzed by means of computational studies leading to the hypothesis of a redox-based mechanism of action with a one-electron reduction bioactivation step and the subsequent formation of a toxic semiquinone species, similarly to what was previously observed for the antiplasmodial activity. Our results also evidenced that the selective toxicity against mammalian cells or parasites as well as specific developmental stages of a parasite can be addressed by varying the nature of the introduced substituents.

Graphical abstract, Structural investigations of receptor-ligand interactions on living cells surface by high-resolution Nuclear Magnetic Resonance (NMR) are problematic due to their short lifetime, which often prevents the acquisition of experiments longer than few hours. To overcome these limitations, we developed an on-cell NMR-based approach for exploring the molecular determinants driving the receptor-ligand recognition mechanism under native conditions. Our method relies on the combination of high-resolution structural and dynamics NMR data with Molecular Dynamics simulations and Molecular Docking studies. The key point of our strategy is the use of Non Uniform Sampling (NUS) and T1ρ-NMR techniques to collect atomic-resolution structural and dynamics information on the receptor-ligand interactions with living cells, that can be used as conformational constraints in computational studies. In fact, the application of these two NMR methodologies allows to record spectra with high S/N ratio and resolution within the lifetime of cells. In particular, 2D NUS [1H–1H] trNOESY spectra are used to explore the ligand conformational changes induced by receptor binding; whereas T1ρ-based experiments are applied to characterize the ligand binding epitope by defining two parameters: T1ρ Attenuation factor and T1ρ Binding Effect. This approach has been tested to characterize the molecular determinants regulating the recognition mechanism of αvβ5-integrin by a selective cyclic binder peptide named RGDechi15D. Our data demonstrate that the developed strategy represents an alternative in-cell NMR tool for studying, at atomic resolution, receptor-ligand recognition mechanism on living cells surface. Additionally, our application may be extremely useful for screening of the interaction profiling of drugs with their therapeutic targets in their native cellular environment.

Ros/MucR is a widespread family of bacterial zinc-finger (ZF) containing proteins that integrate multiple functions such as virulence, symbiosis and/or cell cycle transcription. NMR solution structure of Ros DNA-binding domain (region 56–142, i.e. Ros87) has been solved by our group and shows that the prokaryotic ZF domain shows interesting structural and functional features that differentiate it from its eukaryotic counterpart as it folds in a significantly larger zinc-binding globular domain. We have recently proposed a novel functional model for this family of proteins suggesting that they may act as H-NS-‘like’ gene silencers. Indeed, the N-terminal region of this family of proteins appears to be responsible for the formation of functional oligomers. No structural characterization of the Ros N-terminal domain (region 1–55) is available to date, mainly because of serious solubility problems of the full-length protein. Here we report the first structural characterization of the N-terminal domain of the prokaryotic ZF family examining by means of MD and NMR the structural preferences of the full-length Ros protein from Agrobacterium tumefaciens.

The present study provides new evidence that cationic porphyrins may be considered as tunable platforms to interfere with the structural &ldquo, key code&rdquo, present on the 20S proteasome &alpha, rings and, by consequence, with its catalytic activity. Here, we describe the functional and conformational effects on the 20S proteasome induced by the cooperative binding of the tri-cationic 5-(phenyl)-10,15,20-(tri N-methyl-4-pyridyl) porphyrin (Tris-T4). Our integrated kinetic, NMR, and in silico analysis allowed us to disclose a complex effect on the 20S catalytic activity depending on substrate/porphyrin concentration. The analysis of the kinetic data shows that Tris-T4 shifts the relative populations of the multiple interconverting 20S proteasome conformations leading to an increase in substrate hydrolysis by an allosteric pathway. Based on our Tris-T4/h20S interaction model, Tris-T4 is able to affect gating dynamics and substrate hydrolysis by binding to an array of negatively charged and hydrophobic residues present on the protein surface involved in the 20S molecular activation by the regulatory proteins (RPs). Accordingly, despite the fact that Tris-T4 also binds to the &alpha, 3&Delta, N mutant, allosteric modulation is not observed since the molecular mechanism connecting gate dynamics with substrate hydrolysis is impaired. We envisage that the dynamic view of the 20S conformational equilibria, activated through cooperative Tris-T4 binding, may work as a simplified model for a better understanding of the intricate network of 20S conformational/functional states that may be mobilized by exogenous ligands, paving the way for the development of a new generation of proteasome allosteric modulators.

Schistosomiasis is the most significant neglected tropical parasitic disease caused by helminths in terms of morbidity and mortality caused by helminths. In this work, we present the antischistosomal activity against

Ligand-based control of protein functional motions can provide novel opportunities in the study of fundamental biological mechanisms and in the development of novel therapeutics. In this work we addressed the ligand-based modulation of integrin functions. Inhibitors of integrin α

[Image: see text] Covalent inhibitors of PfGAPDH characterized by a 3-bromoisoxazoline warhead were developed, and their mode of interaction with the target enzyme was interpreted by means of molecular modeling studies: some of them displayed a submicromolar antiplasmodial activity against both chloroquine sensitive and resistant strains of Plasmodium falciparum, with good selectivity indices.

A small library of antiplasmodial methoxy-thiazinoquinones, rationally designed on the model of the previously identified hit 1, has been prepared by a simple and inexpensive procedure. The synthetic derivatives have been subjected to in vitro pharmacological screening, including antiplasmodial and toxicity assays. These studies afforded a new lead candidate, compound 9, endowed with higher antiplasmodial potency compared to 1, a good selectivity index when tested against a panel of mammalian cells, no toxicity against RBCs, a synergistic antiplasmodial action in combination with dihydroartemisinin, and a promising inhibitory activity on stage V gametocyte growth. Computational studies provided useful insights into the structural requirements needed for the antiplasmodial activity of thiazinoquinone compounds and on their putative mechanism of action.

Lactobacillus plantarum is one of the most predominant species in the human gut microbiota of healthy individuals. We have previously characterized some probiotic features of L. plantarum LM3, as the high resistance to different stress, the binding ability toward some extracellular matrix proteins and plasminogen and the immunomodulatory role of the surface expressed adhesin EnoA1. We have also identified the flmA, flmB and flmC genes, coding for putative proteins named FlmA, FlmB and FlmC, whose null mutations partially impaired biofilm development, the L. plantarum LM3&ndash, 6 strain, carrying a deletion in flmC, showed a high rate of autolysis, supporting the hypothesis that FlmC might be involved in cell wall integrity. Here, we report the in-silico characterization of &Delta, TM-FlmC, a portion of the FlmC protein. The protein has been also expressed, purified and characterized by means of CD spectroscopy, ICP-mass and UHPLC-HRMS. The obtained experimental data validated the predicted model unveiling also the presence of a bound lipid molecule and of a Mg(II) ion. Overall, we provide strong evidences that &Delta, TM-FlmC belongs to the LytR-CpsA-Psr (LCP) family of domains and is involved in cell envelope biogenesis.

We present a novel comparative analysis of representative protein kinases to characterize the main dynamic and energetic determinants of functional regulation shared among different families. The relationships between stability and plasticity are also used to rationalize kinase tendencies to interact with the molecular chaperone Hsp90. These questions are tackled through newly developed molecular-dynamics-based methods of analysis of internal energy and dynamics applied to a total of 37 different systems, which represent wild-type and mutated proteins, including active and inactive states. Energetic decomposition analysis is coupled to multiple structural alignments and dynamic decomposition methods and identifies, across different families, common elements that underlie fold stabilization and conformational regulation. This analysis also exposes which substructures play a key role in determining chaperone dependence. Overall, the results highlight common interaction networks that underpin kinase stabilization, are modulated by mutations (even if located at a distance), and underlie their tendencies to act as clients or nonclients of Hsp90.

The rational design of novel biomolecules with desired functional properties is one of the most fascinating challenges in science, with implications at the fundamental and practical levels. From the fundamental point of view, the design of proteins able to support nonnatural reactivities represents the decisive test on our understanding of the molecular mechanisms through which biomolecules operate. From the practical point of view, new designs may open the way to applications in a wide variety of fields, ranging from health to life science, and from catalysis to material sciences. During the past decades, we have witnessed an amazing transition in the application of computational methods to protein and enzyme design, from simple fold prediction to ab initio structural design. Herein, we review key areas and fundamental aspects of research in the design of protein structures, interactions, and reactivities. We also provide our perspective on the exciting range of developments that are made possible by the integration of innovations in hardware, software, and theory, while keeping an eye on the applications, challenges, and opportunities that can open up in many different domains of science. For further resources related to this article, please visit the WIREs website.

Understanding how binding events modulate functional motions of multidomain proteins is a major issue in chemical biology. We address several aspects of this problem by analyzing the differential dynamics of αvβ3 integrin bound to wild type (wtFN10, agonist) or high affinity (hFN10, antagonist) mutants of fibronectin. We compare the dynamics of complexes from large-scale domain motions to inter-residue coordinated fluctuations to characterize the distinctive traits of conformational evolution and shed light on the determinants of differential αvβ3 activation induced by different FN sequences. We propose an allosteric model for ligand-based integrin modulation: the conserved integrin binding pocket anchors the ligand, while different residues on the two FN10’s act as the drivers that reorganize relevant interaction networks, guiding the shift towards inactive (hFN10-bound) or active states (wtFN10-bound). We discuss the implications of results for the design of integrin inhibitors., Author Summary Interactions between proteins are at the basis of all biological processes in the cell. In this context, the study of conformational responses of protein receptors to the binding of endogenous ligands may be a source of inspiration for the design of small molecule modulators that permit to control such biological processes. To progress towards this goal, we need to unravel the molecular determinants that underlie the correlations between sequences, structures and the onset of functional motions in the receptor, eventually illuminating the links between the fine atomic-scale protein-ligand interactions and the large-scale protein motions. Herein, we have concentrated on the multidomain receptor αvβ3 integrin bound to two different sequences of the endogenous ligand fibronectin: the wild type one, wtFN10, which acts as an agonist activating the receptor, and a high affinity mutant, hFN10, which acts as a true antagonist inhibiting the receptor. Through the comparative analysis of several dynamic descriptors at different levels of resolution, from the residue to domain level, we shed light on the salient conformational dynamics differences determined by fibronectin sequence mutations: we show that it is possible to identify interaction hotspots in the integrin binding site that specifically respond to the fibronectin sequence variations, and allosterically drive conformational changes towards integrin activation (in the case of wtFN10 binding) or inhibition (hFN10 binding). Finally we propose an allosteric model of integrin regulation that can be used in the design of small molecule integrin inhibitors or modulators.

The chemical structure of RNA and DNA is very similar; however, the three-dimensional conformation of these two nucleic acids is very different. Whereas the DNA adopts a repetitive structure of a double-stranded helix, RNA is primarily single stranded with a complex three-dimensional structure in which the hairpin is the most common secondary structure. Apart from the difference between uracil and thymine, the difference in the chemical structure between RNA and DNA is the presence of a hydroxyl group at position 2' of the sugar (ribose) instead of a hydrogen (deoxyribose). In this paper, we present molecular dynamics simulations addressing the contribution of 2'-hydroxyls to the stability of a GCUAA pentaloop motif. The results indicate that the 2'-hydroxyls stabilize the hairpin conformation of the GCUAA pentaloop relative to an analogous oligonucleotide in which the ribose sugars in the loop region were substituted with deoxyriboses. The magnitude of the stabilization was found to be 23.8 ± 4.1 kJ/mol using an alchemical mutations free energy method and 4.2 ± 6.5 kJ/mol using potential of mean force calculations. The latter indicates that in addition to its larger thermodynamic stability the RNA hairpin is also kinetically more stable. We find that the excess stability is a result of intrahairpin hydrogen bonds in the loop region between the 2'-hydroxyls and sugars, bases, and phosphates. The hydrogen bonds with the sugars and phosphates involve predominantly interactions with adjacent nucleotides. However, the hydrogen bonds with the bases involve also interactions between groups on opposite sides of the loop or with the middle base of the loop and are therefore likely to contribute significantly to the stability of the loop. Of these hydrogen bonds, the most frequent is observed between the 2'-hydroxyl at the first position of the pentaloop with N6/N7 of adenine at the forth position, as well as between the 2'-hydroxyl at position -1 with N6 of adenine at the fifth position. Our results contribute to the notion that one of the important roles of the ribose sugars in RNA is to facilitate hairpin formation.

Proteins are not static objects. To carry out their functions in the cells and participate in biochemical interaction networks, proteins have to explore different conformational substates, which favor the adaptation to different partners and ultimately allow them to respond to changes in the environment. In this paper we discuss the implications of including the atomistic description of protein dynamics and flexibility in the context of drug discovery and design. The underlying idea is that a better understanding of the atomistic details of molecular recognition phenomena and conformational cross-talk between a ligand and a receptor can in fact translate in unexplored opportunities for the discovery of new drug like molecules. We will illustrate and discuss dynamics-based pharmacophores, the discovery of cryptic binding sites, the characterization and exploitation of allosteric regulation mechanisms and the definition of potential protein-protein interaction sites as potential sources of new bases for the rational design of small molecules endowed with specific biological functions. Overall, the inclusion of protein flexibility in the drug discovery process is starting to attract attention not only in the academic but also in the industrial community. This is supported by experimental tests that prove the actual feasibility of considering the explicit dynamics of drug-protein interactions at all relevant levels of resolution and the use of multiple receptor conformations in drug discovery, as affordable complements (if not an alternative) to classical High Throughput Screening (HTS) efforts based on static structures.

Racemic and enantiopure benzofuranmethanamines 5a – c have been reacted with N -Boc-3-(4-cyanophenyl)oxaziridine to give N -Boc-hydrazines 7a – c , which have in turn been transformed by deprotection and cyclisation into triazoles 4a – c , potent antiaromatase agents, in good overall yield and with high enantiomeric excess.

RNA oligonucleotides exhibit a large tendency to bend and form a loop conformation which is a major motif contributing to their complex three-dimensional structure. This is in contrast to DNA molecules that predominantly form the double-helix structure. In this paper we investigate by molecular dynamics simulation, as well as, by its combination with the replica-exchange method, the propensity of RNA chains containing the GCUAA pentaloop to form spontaneously a hairpin conformation. The results were then compared with those of analogous hybrid oligonucleotides in which the ribose groups in the loop-region were substituted by deoxyriboses. We find that the RNA oligomers exhibit a marginal excess stability to form loop structures. The equilibrium constant for opening the loop to an extended conformation is twice as large in the hybrid than it is in the RNA chain. Analyses of the hydrogen bonds indicate that the excess stability for forming a hairpin is a result of hydrogen bonds the 2′-hydroxyls in the loop region form with other groups in the loop. Of these hydrogen bonds, the most important is the hydrogen bond donated from the 2′-OH at the first position of the loop to N7 of adenine at the forth position. RNA and DNA backbones are characterized by different backbone dihedral angles and sugar puckering that can potentially facilitate or hamper the hydrogen bonds involving the 2′-OH. Nevertheless, the sugar puckerings of all the pentaloop nucleotides were not significantly different between the two chains displaying the C3′-endo conformation characteristic to the A-form double helix. All of the other backbone dihedrals also did not show any considerable difference in the loop-region except of the δ -dihedral. In this case, the RNA loop exhibited bimodal distributions corresponding to, both, the RNA and DNA backbones, whereas the loop of the hybrid chain behaved mostly as that of a DNA backbone. Thus, it is possible that the behavior of the δ -dihedrals in the loop-region of the RNA adopts conformations that facilitate the intra-nucleotide hydrogen bondings of the 2′-hydroxyls, and consequently renders loop structures in RNA more stable.

Sirtuin family, in humans as well as in all mammalia, is composed by seven different homologous proteins with NAD-dependent deacetylase/ADP-ribosyltransferase activity. Numerous studies have determined their cellular location and their biological functions. In particular, Sirt-1 is defined as a nuclear protein involved in the molecular mechanisms of inflammation and neurodegeneration through the de-acetylation of many different substrates (PGC-α, FOXOs, NFκB). However experimental data in mouse suggest both its cytoplasmatic presence and nucleo-cytoplasmic shuttling upon oxidative stress. Recently we have modeled the three-dimensional structure of human Sirt-1, and highlighted that it is composed by four different regions: N-terminal region, allosteric site, catalytic core and C-terminal region and underlined that the two terminal regions have high intrinsic disorder propensity. Since Sirt-1 is implicated in various diseases and cancers, many different papers report experimental studies related to its functional activators. The aim of this article is i) to present interactomic studies on human SIRT-1 to understand its most important functional relationships in the light of the gene-protein interactions that control major metabolic pathways and ii) to show how this protein binds some activator molecules in order to evidence structural determinants, physico-chemical features and those residues involved in the formation of complexes. It is believed that these data will be useful to synthesize new effectors through drug design approaches.

Lukasz Kowalczyk et al., Transporters of the amino acid, polyamine and organocation (APC) superfamily play essential roles in cell redox balance, cancer, and aminoacidurias. The bacterial L-arginine/agmatine antiporter, AdiC, is the main APC structural paradigm and shares the >5 + 5 inverted repeat> fold found in other families like the Na+-coupled neurotransmitter transporters. The available AdiC crystal structures capture two states of its transport cycle: the open-to-out apo and the outward-facing Arg +-bound occluded. However, the role of Arg+ during the transition between these two states remains unknown. Here, we report the crystal structure at 3.0 Å resolution of an Arg+-bound AdiC mutant (N101A) in the open-to-out conformation, completing the picture of the major conformational states during the transport cycle of the 5 + 5 inverted repeat fold-transporters. The N101A structure is an intermediate state between the previous known AdiC conformations. The Arg+-guanidinium group in the current structure presents high mobility and delocalization, hampering substrate occlusion and resulting in a low translocation rate. Further analysis supports that proper coordination of this group with residues Asn101 and Trp293 is required to transit to the occluded state, providing the first clues on the molecular mechanism of substrate-induced fit in a 5 + 5 inverted repeat fold-transporter. The pseudosymmetry found between repeats in AdiC, and in all fold-related transporters, restraints the conformational changes, in particular the transmembrane helices rearrangements, which occur during the transport cycle. In AdiC these movements take place away from the dimer interface, explaining the independent functioning of each subunit., Este trabajo ha sido financiado parcialmente por el Ministerio de España de Ciencia e Innovación (SAF2009-12606-C02-01 (M.P.), BFU2008-04637 (J.L.V.-I.), y BFU2009-09268 (I.F., X.C.)), por el Proyecto Consolider E-Science y la Fundacion Marcelino Botin (M.O.). Financiación adicional de la Generalitat de Catalunya (SGR2009-1355 (M.P.) y BP-B2008-00239 (X.C.) )

Miraculin differs from other sweet-tasting proteins because it is a taste-modifier having the unusual property of modifying sourness into sweetness. Its dimer is covalently linked by an inter-chain disulphide bond, and shows its taste-modifying activity at acidic pH, with maximum at pH 3.0, while it is flat at neutral pH. Previous studies suggested the importance of two histidine residues for the taste-modifying activity of miraculin. In this work, we have conducted molecular dynamics simulations on wild type miraculin and on three mutated dimers (H29A, H59A and H29A/H59A) both at neutral and acidic pH to investigate the structural and functional role of these two His residues. Our results suggested that at acidic pH the presence of two charged His at the interface induced a structural rearrangement of the two monomers, thus leading to their relative opening and the following adaptation of their conformation to the receptor surface. On the other hand the simulations on three mutants showed that the mutated dimers had a closed form, and highlighted the important role of H29 in stabilizing/destabilizing the dimer arrangement and also a cooperative effect of the two histidines.

Neurotrophins (NTs) represent a family of proteins that play an important role in the survival, development, and function of neurons. Extensive efforts are currently being made to develop small molecules endowed with agonist or antagonist NT activity. The structurally versatile N-termini of these proteins are considered regions of interest for the design of new molecules. By combining experimental and computational approaches, we analyzed the intrinsic conformational preferences of the N-termini of two of the most important NTs: NGF (NGF-Nter) and NT4 (NT4-Nter). Circular dichroism spectra clearly indicate that both peptides show a preference for random coil states. Because this finding does not preclude the possibility that structured forms may occur in solution as minor conformational states, we performed molecular-dynamics simulations to gain insights into the structural features of populated species. In line with the circular dichroism analysis, the simulations show a preference for unstructured states for both peptides. However, the simulations also show that for NT4-Nter, and to a lesser extent for NGF-Nter, helical conformations, which are required for binding to the Trk receptor, are present in the repertoire of structures that are intrinsically accessible to these peptides. Accordingly, molecular recognition of NTs by the Trk receptor is accomplished by the general mechanism known as population shift. These findings provide a structural rationale for the observed activity of synthetic peptides based on these NT regions. They also suggest strategies for the development of biologically active peptide-based compounds.

Nine genetically inherited neurodegenerative diseases are linked to abnormal expansions of a polyglutamine (polyQ) encoding region. Over the years, several structural models for polyQ regions have been proposed and confuted. The cross-β-spine steric zipper motif, identified recently for the GNNQQNY peptide, represents an attractive model for amyloid fibers formed by polyQ fragments. Here we report a detailed molecular dynamics investigation of polyQ models assembled by cross-β-spine steric zipper motifs. Our simulations indicate clearly that these assemblies are very stable. Glutamine side chains contribute strongly to the overall stability of the models by fitting perfectly within the zipper. In contrast to GNNQQNY zipper motifs, hydrogen bonding interactions provide a significant contribution to the overall stability of polyQ models. Molecular dynamics simulations carried out on monomeric polyQ forms (composed by 40–60 residues) show clearly that they can also assume structures stabilized by steric zipper motifs. Based on these findings, we build monomeric polyQ models that can explain recent data on the toxicity exerted by these species. In a more general context, our data suggests that polyQ models with interdigitated side chains can provide a structural rationale to several literature experiments on polyQ formation, stability, and toxicity.

The center of mass of a protein is an artificial point useful for detecting important and simple features of proteins structure, shape and association. CALCOM is a software which calculates the center of mass of a protein, starting from PDB protein structure files. In the case of protein complexes and of protein-small ligand complexes, the position of protein residues or of ligand atoms respect to each protein subunit can be evaluated, as well as the distance among the center of mass of the protein subunits, in order to compare different conformations and evaluate the relative motion of subunits. Availability The service is available at the URL: http://bioinformatica.isa.cnr.it/CALCOM/.