Springer Science and Business Media LLC

To understand the effect of the long intracellular loop 3 (ICL3) on the intrinsic dynamics of human β2-adrenergic receptor, molecular dynamics (MD) simulations were performed on two different models, both of which were based on the inactive crystal structure in complex with carazolol (after removal of carazolol and T4-lysozyme). In the so-called loop model, the ICL3 region that is missing in available crystal structures was modeled as an unstructured loop of 32-residues length, whereas in the clipped model, the two open ends were covalently bonded to each other. The latter model without ICL3 was taken as a reference, which has also been commonly used in recent computational studies. Each model was embedded into POPC bilayer membrane with explicit water and subjected to a 1 μs molecular dynamics (MD) simulation at 310 K. After around 600 ns, the loop model started a transition to a “very inactive” conformation, which is characterized by a further movement of the intracellular half of transmembrane helix 6 (TM6) towards the receptor core, and a close packing of ICL3 underneath the membrane completely blocking the G-protein’s binding site. Concurrently, the binding site at the extracellular part of the receptor expanded slightly with the Ser207-Asp113 distance increasing to 18 Å from 11 Å, which was further elaborated by docking studies. The essential dynamics analysis indicated a strong coupling between the extracellular and intracellular parts of the intact receptor, implicating a functional relevance for allosteric regulation. In contrast, no such transition to the “very inactive” state, nor any structural correlation, was observed in the clipped model without ICL3. Furthermore, elastic network analysis using different conformers for the loop model indicated a consistent picture on the specific ICL3 conformational change being driven by global modes.

Transgenic proteins expressed by genetically modified food crops are evaluated for their potential allergenic properties prior to marketing, among others by identification of short identical amino acid sequences that occur both in the transgenic protein and allergenic proteins. A strategy is proposed, in which the positive outcomes of the sequence comparison with a minimal length of six amino acids are further screened for the presence of potential linear IgE-epitopes. This double track approach involves the use of literature data on IgE-epitopes and an antigenicity prediction algorithm. Thirty-three transgenic proteins have been screened for identities of at least six contiguous amino acids shared with allergenic proteins. Twenty-two transgenic proteins showed positive results of six- or seven-contiguous amino acids length. Only a limited number of identical stretches shared by transgenic proteins (papaya ringspot virus coat protein, acetolactate synthase GH50, and glyphosate oxidoreductase) and allergenic proteins could be identified as (part of) potential linear epitopes. Many transgenic proteins have identical stretches of six or seven amino acids in common with allergenic proteins. Most identical stretches are likely to be false positives. As shown in this study, identical stretches can be further screened for relevance by comparison with linear IgE-binding epitopes described in literature. In the absence of literature data on epitopes, antigenicity prediction by computer aids to select potential antibody binding sites that will need verification of IgE binding by sera binding tests. Finally, the positive outcomes of this approach warrant further clinical testing for potential allergenicity.

Dengue virus along with the other members of the flaviviridae family has reemerged as deadly human pathogens. Understanding the mechanistic details of these infections can be highly rewarding in developing effective antivirals. During maturation of the virus inside the host cell, the coat proteins E and M undergo conformational changes, altering the morphology of the viral coat. However, due to low resolution nature of the available 3-D structures of viral assemblies, the atomic details of these changes are still elusive. In the present analysis, starting from Cα positions of low resolution cryo electron microscopic structures the residue level details of protein-protein interaction interfaces of dengue virus coat proteins have been predicted. By comparing the preexisting structures of virus in different phases of life cycle, the changes taking place in these predicted protein-protein interaction interfaces were followed as a function of maturation process of the virus. Besides changing the current notion about the presence of only homodimers in the mature viral coat, the present analysis indicated presence of a proline-rich motif at the protein-protein interaction interface of the coat protein. Investigating the conservation status of these seemingly functionally crucial residues across other members of flaviviridae family enabled dissecting common mechanisms used for infections by these viruses. Thus, using computational approach the present analysis has provided better insights into the preexisting low resolution structures of virus assemblies, the findings of which can be made use of in designing effective antivirals against these deadly human pathogens.

Previous kinetic investigations of fungal-peroxidase catalyzed oxidation of N-aryl hydroxamic acids (AHAs) and N-aryl-N- hydroxy urethanes (AHUs) revealed that the rate of reaction was independent of the formal redox potential of substrates. Moreover, the oxidation rate was 3–5 orders of magnitude less than for oxidation of physiological phenol substrates, though the redox potential was similar. To explain the unexpectedly low reactivity of AHAs and AHUs we made ab initio calculations of the molecular structure of the substrates following in silico docking in the active center of the enzyme. AHAs and AHUs were docked at the distal side of heme in the sites formed by hydrophobic amino acid residues that retarded a proton transfer and finally the oxidation rate. The analogous phenol substrates were docked at different sites permitting fast proton transfer in the relay of distal His and water that helped fast substrate oxidation.

Superoxide dismutases (SODs) are important enzymes in defence against oxidative stress. In Plasmodium falciparum, they may be expected to have special significance since part of the parasite life cycle is spent in red blood cells where the formation of reactive oxygen species is likely to be promoted by the products of haemoglobin breakdown. Thus, inhibitors of P. falciparum SODs have potential as anti-malarial compounds. As a step towards their development we have determined the crystal structure of the parasite's cytosolic iron superoxide dismutase. The cytosolic iron superoxide dismutase from P. falciparum (Pf FeSOD) has been overexpressed in E. coli in a catalytically active form. Its crystal structure has been solved by molecular replacement and refined against data extending to 2.5 Å resolution. The structure reveals a two-domain organisation and an iron centre in which the metal is coordinated by three histidines, an aspartate and a solvent molecule. Consistent with ultracentrifugation analysis the enzyme is a dimer in which a hydrogen bonding lattice links the two active centres. The tertiary structure of Pf FeSOD is very similar to those of a number of other iron-and manganese-dependent superoxide dismutases, moreover the active site residues are conserved suggesting a common mechanism of action. Comparison of the dimer interfaces of Pf FeSOD with the human manganese-dependent superoxide dismutase reveals a number of differences, which may underpin the design of parasite-selective superoxide dismutase inhibitors.

Protein deformation has been extensively analysed through global methods based on RMSD, torsion angles and Principal Components Analysis calculations. Here we use a local approach, able to distinguish among the different backbone conformations within loops, α-helices and β-strands, to address the question of secondary structures' shape variation within proteins and deformation at interface upon complexation. Using a structural alphabet, we translated the 3 D structures of large sets of protein-protein complexes into sequences of structural letters. The shape of the secondary structures can be assessed by the structural letters that modeled them in the structural sequences. The distribution analysis of the structural letters in the three protein compartments (surface, core and interface) reveals that secondary structures tend to adopt preferential conformations that differ among the compartments. The local description of secondary structures highlights that curved conformations are preferred on the surface while straight ones are preferred in the core. Interfaces display a mixture of local conformations either preferred in core or surface. The analysis of the structural letters transition occurring between protein-bound and unbound conformations shows that the deformation of secondary structure is tightly linked to the compartment preference of the local conformations. The conformation of secondary structures can be further analysed and detailed thanks to a structural alphabet which allows a better description of protein surface, core and interface in terms of secondary structures' shape and deformation. Induced-fit modification tendencies described here should be valuable information to identify and characterize regions under strong structural constraints for functional reasons.

Recently a new class of methods for fast protein structure comparison has emerged. We call the methods in this class projection methods as they rely on a mapping of protein structure into a high-dimensional vector space. Once the mapping is done, the structure comparison is reduced to distance computation between corresponding vectors. As structural similarity is approximated by distance between projections, the success of any projection method depends on how well its mapping function is able to capture the salient features of protein structure. There is no agreement on what constitutes a good projection technique and the three currently known projection methods utilize very different approaches to the mapping construction, both in terms of what structural elements are included and how this information is integrated to produce a vector representation. In this paper we propose a novel projection method that uses secondary structure information to produce the mapping. First, a diverse set of spatial arrangements of triplets of secondary structure elements, a set of structural models, is automatically selected. Then, each protein structure is mapped into a high-dimensional vector of "counts" or footprint, where each count corresponds to the number of times a given structural model is observed in the structure, weighted by the precision with which the model is reproduced. We perform the first comprehensive evaluation of our method together with all other currently known projection methods. The results of our evaluation suggest that the type of structural information used by a projection method affects the ability of the method to detect structural similarity. In particular, our method that uses the spatial conformations of triplets of secondary structure elements outperforms other methods in most of the tests.

Structural switches upon binding of phosphorylated moieties underpin many signalling networks. The ligand activation is a form of allosteric modulation of the protein, where the binding site is remote from the structural change in the protein. Recently this structural switch has been elegantly demonstrated with the crystallisation of the activated form of 3-phosphoinositide-dependent protein kinase-1 (PDK1). The purpose of the present work is to determine whether the allosteric coupling in PDK1 emerges at the level of a simple coarse grained model of protein dynamics. It is shown here that the allosteric effects of the agonist binding to the small lobe upon the activation loop in the large lobe of PDK1 are explainable within a simple 'ball and spring' elastic network model (ENM) of protein dynamics. In particular, the model shows that the bound phospho peptide mimetic fluctuations have a high degree of correlation with the activation loop of PDK1. The ENM approach to small molecule activation of proteins may offer a first pass predictive methodology where affinity is encoded in residues remote from the active site, and aid in the design of specific protein agonists that enhance the allosteric coupling and antagonist that repress it.

Differences in sweet taste perception among species depend on structural variations of the sweet taste receptor. The commercially used isovanillyl sweetener neohesperidin dihydrochalcone activates the human but not the rat sweet receptor TAS1R2+TAS1R3. Analysis of interspecies combinations and chimeras of rat and human TAS1R2+TAS1R3 suggested that the heptahelical domain of human TAS1R3 is crucial for the activation of the sweet receptor by neohesperidin dihydrochalcone. By mutational analysis combined with functional studies and molecular modeling we identified a set of different amino acid residues within the heptahelical domain of human TAS1R3 that forms the neohesperidin dihydrochalcone binding pocket. Sixteen amino acid residues in the transmembrane domains 2 to 7 and one in the extracellular loop 2 of hTAS1R3 influenced the receptor's response to neohesperidin dihydrochalcone. Some of these seventeen residues are also part of the binding sites for the sweetener cyclamate or the sweet taste inhibitor lactisole. In line with this observation, lactisole inhibited activation of the sweet receptor by neohesperidin dihydrochalcone and cyclamate competitively, whereas receptor activation by aspartame, a sweetener known to bind to the N-terminal domain of TAS1R2, was allosterically inhibited. Seven of the amino acid positions crucial for activation of hTAS1R2+hTAS1R3 by neohesperidin dihydrochalcone are thought to play a role in the binding of allosteric modulators of other class C GPCRs, further supporting our model of the neohesperidin dihydrochalcone pharmacophore. From our data we conclude that we identified the neohesperidin dihydrochalcone binding site at the human sweet taste receptor, which overlaps with those for the sweetener cyclamate and the sweet taste inhibitor lactisole. This readily delivers a molecular explanation of our finding that lactisole is a competitive inhibitor of the receptor activation by neohesperidin dihydrochalcone and cyclamate. Some of the amino acid positions crucial for activation of hTAS1R2+hTAS1R3 by neohesperidin dihydrochalcone are involved in the binding of allosteric modulators in other class C GPCRs, suggesting a general role of these amino acid positions in allosterism and pointing to a common architecture of the heptahelical domains of class C GPCRs.