Search

Modern computational chemistry methods provide a powerful tool for use in refining the geometry of proteins determined by X-ray crystallography. Specifically, computational methods can be used to correctly place hydrogen atoms unresolved by this experimental method and improve bond geometry accuracy. Using the semiempirical method PM7, the structure of the nucleotide-sanitizing enzyme MTH1, complete with hydrolyzed substrate 8-oxo-dGMP, was optimized and the resulting geometry compared with the original X-ray structure of MTH1. After determining hydrogen atom placement and the identification of ionized sites, the charge distribution in the binding site was explored. Where comparison was possible, all the theoretical predictions were in good agreement with experimental observations. However, when these were combined with additional predictions for which experimental observations were not available, the result was a new and alternative description of the substrate-binding site interaction. An estimate was made of the strengths and weaknesses of the PM7 method for modeling proteins on varying scales, ranging from overall structure to individual interatomic distances. An attempt to correct a known fault in PM7, the under-estimation of steric repulsion, is also described. This work sheds light on the specificity of the enzyme MTH1 toward the substrate 8-oxo-dGTP; information that would facilitate drug development involving MTH1. Graphical Abstract Overlay of the backbone traces of the two MTH1 protein chains (green and orange respectively) in PDB 3ZR0 and the equivalent PM7 structures (magenta and cyan respectively) each optimized separately.

Substrate permissiveness has long been regarded as the raw materials for the evolution of new enzymatic functions. In land plants, hydroxycinnamoyltransferase (HCT) is an essential enzyme of the phenylpropanoid metabolism. Although essential enzymes are normally associated with high substrate specificity, HCT can utilize a variety of non-native substrates. To examine the structural and dynamic basis of substrate permissiveness in this enzyme, we report the crystal structure of HCT from Selaginella moellendorffii and molecular dynamics (MD) simulations performed on five orthologous HCTs from several major lineages of land plants. Through altogether 17-μs MD simulations, we demonstrate the prevalent swing motion of an arginine handle on a submicrosecond timescale across all five HCTs, which plays a key role in native substrate recognition by these intrinsically promiscuous enzymes. Our simulations further reveal how a non-native substrate of HCT engages a binding site different from that of the native substrate and diffuses to reach the catalytic center and its co-substrate. By numerically solving the Smoluchowski equation, we show that the presence of such an alternative binding site, even when it is distant from the catalytic center, always increases the reaction rate of a given substrate. However, this increase is only significant for enzyme-substrate reactions heavily influenced by diffusion. In these cases, binding non-native substrates ‘off-center’ provides an effective rationale to develop substrate permissiveness while maintaining the native functions of promiscuous enzymes., Author summary Examples abound of enzymes that can process substrates other than their native ones. However, the structural and dynamic basis of this promiscuity remains to be fully understood. In this work, we examine HCT, an intrinsically promiscuous acyltransferase with conserved function in all land plants. We uncover the sub-microsecond swing motion of a key arginine residue facilitating the recognition of both native and non-native substrates of HCT. We also quantify the impact of an off-center binding site on the non-native reaction rate. Although our calculations were inspired by HCT, the results apply in general, i.e., for enzymes heavily influenced by diffusion, binding non-native substrates ‘off-center’, even with rather weak affinity, can accelerate non-native reactions to appreciable levels.

RIP2, one of the RIP kinases, interacts with p75 neurotrophin receptor, regulating the neuron survival, and with NOD1 and NOD2 proteins, causing the innate immune response against gram-negative and gram-positive bacteria via its caspase recruitment domain (CARD). This makes RIP2 a prospective target for novel therapies, aimed to modulate the inflammatory diseases and neurogenesis/neurodegeneration. Several studies report the problems with the stability of human RIP2 CARD and its production in bacterial hosts, which is a prerequisite for the structural investigation with solution NMR spectroscopy. In the present work, we report the high yield production and refolding protocols and resolve the structure of rat RIP2 CARD. The structure reveals the important differences to the previously published conformation of the homologous human protein. Using solution NMR, we characterized the intramolecular mobility and pH-dependent behavior of RIP2 CARD, and found the propensity of the protein to form high-order oligomers at physiological pH while being monomeric under acidic conditions. The oligomerization of protein may be explained, based on the electrostatic properties of its surface. Analysis of the structure and sequences of homologous proteins reveals the residues which are significant for the unusual fold of RIP2 CARD domains from different species. The high-throughput protein production/refolding protocols and proposed explanation for the protein oligomerization, provide an opportunity to design the stabilized variants of RIP2 CARD, which could be used to study the structural details of RIP2/NOD1/NOD2 interaction and perform the rational drug design.

[EN] Comparison of the proteins of thermophilic, mesophilic, and psychrophilic prokaryotes has revealed several features characteristic to proteins adapted to high temperatures, which increase their thermostability. These characteristics include a profusion of disulfide bonds, salt bridges, hydrogen bonds, and hydrophobic interactions, and a depletion in intrinsically disordered regions. It is unclear, however, whether such differences can also be observed in eukaryotic proteins or when comparing proteins that are adapted to temperatures that are more subtly different. When an organism is exposed to high temperatures, a subset of its proteins is overexpressed (heat-induced proteins), whereas others are either repressed (heat-repressed proteins) or remain unaffected. Here, we determine the expression levels of all genes in the eukaryotic model system Arabidopsis thaliana at 22 and 37 degrees C, and compare both the amino acid compositions and levels of intrinsic disorder of heat-induced and heat-repressed proteins. We show that, compared to heat-repressed proteins, heat-induced proteins are enriched in electrostatically charged amino acids and depleted in polar amino acids, mirroring thermophile proteins. However, in contrast with thermophile proteins, heat-induced proteins are enriched in intrinsically disordered regions, and depleted in hydrophobic amino acids. Our results indicate that temperature adaptation at the level of amino acid composition and intrinsic disorder can be observed not only in proteins of thermophilic organisms, but also in eukaryotic heat-induced proteins; the underlying adaptation pathways, however, are similar but not the same., D.A.-P. and F.F. were supported by funds from the University of Nevada, Reno, and by pilot grants from Nevada INBRE (P20GM103440) and the Smooth Muscle Plasticity COBRE from the University of Nevada, Reno (5P30GM110767-04), both funded by the National Institute of General Medical Sciences (National Institutes of Health). M.X.R.-G. and M.A.F. were supported by grants from Science Foundation Ireland (12/IP/1637) and the Spanish Ministerio de Economia y Competitividad, Spain (MINECO-FEDER; BFU201236346 and BFU2015-66073-P) to MAF. MXRG was supported by a JAE DOC fellowship from the MINECO, Spain. F.V.-S. and M.A.P.-A. were supported by grant BIO2014-55946-P from MINECO-FEDER.

Made available in DSpace on 2018-11-26T17:51:44Z (GMT). No. of bitstreams: 0 Previous issue date: 2018-05-21 Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) Redox Processes in Biomedicine Organic hydroperoxide resistance (Ohr) enzymes are highly efficient Cys-based peroxidases that play central roles in bacterial response to fatty acid hydroperoxides and peroxynitrite, two oxidants that are generated during host-pathogen interactions. In the active site of Ohr proteins, the conserved Arg (Arg19 in Ohr from Xylella fastidiosa) and Glu (Glu51 in Ohr from Xylella fastidiosa) residues, among other factors, are involved in the extremely high reactivity of the peroxidatic Cys (C-p) toward hydroperoxides. In the closed state, the thiolate of C-p is in close proximity to the guanidinium group of Arg19. Ohr enzymes can also assume an open state, where the loop containing the catalytic Arg is far away from C-p and Glu51. Here, we aimed to gain insights into the putative structural switches of the Ohr catalytic cycle. First, we describe the crystal structure of Ohr from Xylella fastidiosa (XfOhr) in the open state that, together with the previously described XfOhr structure in the closed state, may represent two snapshots along the coordinate of the enzyme-catalyzed reaction. These two structures were used for the experimental validation of molecular dynamics (MD) simulations. MD simulations employing distinct protonation states and in silico mutagenesis indicated that the polar interactions of Arg19 with Glu51 and C-p contributed to the stabilization of XfOhr in the closed state. Indeed, C-p oxidation to the disulfide state facilitated the switching of the Arg19 loop from the closed to the open state. In addition to the Arg19 loop, other portions of XfOhr displayed high mobility, such as a loop rich in Gly residues. In summary, we obtained a high correlation between crystallographic data, MD simulations and biochemical/enzymatic assays. The dynamics of the Ohr enzymes are unique among the Cys-based peroxidases, in which the active site Arg undergoes structural switches throughout the catalytic cycle, while C-p remains relatively static. Univ Sao Paulo, Inst Quim, Dept Quim Fundamental, Sao Paulo, SP, Brazil Univ Sao Paulo, Inst Biociencias, Dept Genet & Biol Evolut, Sao Paulo, SP, Brazil Univ Estadual Paulista, Inst Biociencias, Campus Litoral Paulista, Sao Vicente, SP, Brazil Univ Fed Sao Paulo, Dept Bioquim, Sao Paulo, SP, Brazil Univ Estadual Paulista, Inst Biociencias, Campus Litoral Paulista, Sao Vicente, SP, Brazil Redox Processes in Biomedicine: 13/07937-8 FAPESP: 2008/07971-3 FAPESP: 2009/12885-1 FAPESP: 2005/50056-6 FAPESP: 2014/01614-5 FAPESP: 2012/06633-2 FAPESP: 2016/12392-9

The protein kinase PKR activated by viral dsRNA, phosphorylates the eIF2α, which inhibit the mechanism of translation initiation. Viral evolved proteins mimicking the eIF2α block its phosphorylation and help in the viral replication. To decipher the molecular basis for the PKR’s substrate and inhibitor interaction mechanisms, we carried the molecular dynamics studies on the catalytic domain of PKR in complex with substrate eIF2α, and inhibitors TAT and K3L. The studies conducted show the altered domain movements of N lobe, which confers open and close state to the substrate-binding cavity. In addition, PKR exhibits variations in the secondary structural transition of the activation loop residues, and inter molecular contacts with the substrate and the inhibitors. Phosphorylation of the P+1 loop at the Thr-451 increases the affinity of the binding proteins exhibiting its role in the phosphorylation events. The implications of structural mechanisms uncovered will help to understand the basis of the evolution of the host-viral and the viral replication mechanisms.

Generating the immune response requires the discrimination of peptides presented by the human leukocyte antigen complex (HLA) through the T-cell receptor (TCR). However, how a single amino acid substitution in the antigen bonded to HLA affects the response of T cells remains uncertain. Hence, we used molecular dynamics computations to analyze the molecular interactions between peptides, HLA and TCR. We compared immunologically reactive complexes with non-reactive and weakly reactive complexes. MD trajectories were produced to simulate the behavior of isolated components of the various p-HLA-TCR complexes. Analysis of the fluctuations showed that p-HLA binding barely restrains TCR motions, and mainly affects the CDR3 loops. Conversely, inactive p-HLA complexes displayed significant drop in their dynamics when compared with its free versus ternary forms (p-HLA-TCR). In agreement, the free non-reactive p-HLA complexes showed a lower amount of salt bridges than the responsive ones. This resulted in differences between the electrostatic potentials of reactive and inactive p-HLA species and larger vibrational entropies in non-elicitor complexes. Analysis of the ternary p-HLA-TCR complexes also revealed a larger number of salt bridges in the responsive complexes. To summarize, our computations indicate that the affinity of each p-HLA complex towards TCR is intimately linked to both, the dynamics of its free species and its ability to form specific intermolecular salt-bridges in the ternary complexes. Of outstanding interest is the emerging concept of antigen reactivity involving its interplay with the HLA head sidechain dynamics by rearranging its salt-bridges, Financial support was provided by the Spanish Ministry of Economy and Competiveness (Grant No. BFU2012-31670/BMC); The Andalusian Government (Grant PAI, BIO198); The Spanish Fund for Health Research (FIS; code PI11/02366, FI12/00189) and the Ramón Areces Foundation.

The effect of N-terminal diproline segment and charged side chains on the stabilization of helical conformation in alaninebased short peptides are examined using molecular dynamics (MD) simulations. The cationic peptides, Ac-Pro(1)-Pro(2)-Ala(3)-Lys(4)- Ala(5)-Lys(6)-Ala(7)-Lys(8)-Ala(9)-NH2 (Ia) and Ac-(D)Pro(1)-Pro(2)-Ala(3)-Lys(4)- Ala(5)-Lys(6)-Ala(7)-Lys(8)-Ala(9)-NH2 (IIa) are examined for the role of lysine side chains on the inducement of helical conformation in alanine-based short peptides. To examine the influence of lysine and glutamic acid in the i, i + 4 arrangement on the stabilization of helical conformation, cationic peptides, Ia and IIa, are modified as ion-pair peptides, Ac-Pro(1)-Pro(2)-Glu(3)- Glu(4)-Ala(5)-Ala(6)-Lys(7)-Lys(8)-Ala(9)-NH2 (Ib) and Ac-(D)Pro(1)Pro(2)- Glu(3)-Glu(4)-Ala(5)-Ala(6)-Lys(7)-Lys(8)-Ala(9)-NH2 (IIb), respectively. MD simulations manifest enhanced occupancies in the a basin of phi, Psi space for ion-pair peptides as compare to cationic peptides. The radial distribution function (RDF) analysis highlight that large side chain substituents of lysine and glutamic acid assist in helix formation by blocking water molecules from solvating backbone CO and NH groups.

To test the existence of the salt bridge and stability of the HIV-1 p17 matrix protein, an E12A (mutated at helix 1) was established to abolish possible electrostatic interactions. The chemical shift perturbation from the comparison between wild type and E12A suggested the existence of an electrostatic interaction in wild type between E12 and H89 (located in helix 4). Unexpectedly, the studies using urea denaturation indicated that the E12A substitution slightly stabilized the protein. The dynamic structure of E12A was examined under physiological conditions by both amide proton exchange and relaxation studies. The quick exchange method of amide protons revealed that the residues with faster exchange were located at the mutated region, around A12, compared to those of the wild-type protein. In addition, some residues at the region of helix 4, including H89, exhibited faster exchange in the mutant. In contrast, the average values of the kinetic rate constants for amide proton exchange for residues located in all loop regions were slightly lower in E12A than in wild type. Furthermore, the analyses of the order parameter revealed that less flexible structures existed at each loop region in E12A. Interestingly, the structures of the regions including the alpha1-2 loop and helix 5 of E12A exhibited more significant conformational exchanges with the NMR time-scale than those of wild type. Under lower pH conditions, for further destabilization, the helix 1 and alpha2-3 loop in E12A became more fluctuating than at physiological pH. Because the E12A mutant lacks the activities for trimer formation on the basis of the analytical ultra-centrifuge studies on the sedimentation distribution of p17 (Fledderman et al. Biochemistry 49, 9551-9562, 2010), it is possible that the changes in the dynamic structures induced by the absence of the E12-H89 interaction in the p17 matrix protein contributes to a loss of virus assembly.

Protein kinases are key regulatory nodes in cellular networks and their function has been shown to be intimately coupled with their structural flexibility. However, understanding the key structural mechanisms of large conformational transitions remains a difficult task. CDK2 is a crucial regulator of cell cycle. Its activity is finely tuned by Cyclin E/A and the catalytic segment phosphorylation, whereas its deregulation occurs in many types of cancer. ATP competitive inhibitors have failed to be approved for clinical use due to toxicity issues raised by a lack of selectivity. However, in the last few years type III allosteric inhibitors have emerged as an alternative strategy to selectively modulate CDK2 activity. In this study we have investigated the conformational variability of CDK2. A low dimensional conformational landscape of CDK2 was modeled using classical multidimensional scaling on a set of 255 crystal structures. Microsecond-scale plain and accelerated MD simulations were used to populate this landscape by using an out-of-sample extension of multidimensional scaling. CDK2 was simulated in the apo-form and in complex with the allosteric inhibitor 8-anilino-1-napthalenesulfonic acid (ANS). The apo-CDK2 landscape analysis showed a conformational equilibrium between an Src-like inactive conformation and an active-like form. These two states are separated by different metastable states that share hybrid structural features with both forms of the kinase. In contrast, the CDK2/ANS complex landscape is compatible with a conformational selection picture where the binding of ANS in proximity of the αC helix causes a population shift toward the inactive conformation. Interestingly, the new metastable states could enlarge the pool of candidate structures for the development of selective allosteric CDK2 inhibitors. The method here presented should not be limited to the CDK2 case but could be used to systematically unmask similar mechanisms throughout the human kinome.

Made available in DSpace on 2018-12-11T17:03:17Z (GMT). No. of bitstreams: 0 Previous issue date: 2016-06-01 Alpha solenoid proteins play a key role in regulating the classical nuclear import pathway, recognizing a target protein and transporting it into the nucleus. Importin-α (Impα) is the solenoid responsible for cargo protein recognition, and it has been extensively studied by X-ray crystallography to understand the binding specificity. To comprehend the main motions of Impα and to extend the information about the critical interactions during carrier-cargo recognition, we surveyed different conformational states based on molecular dynamics (MD) and normal mode (NM) analyses. Our model of study was a crystallographic structure of Impα complexed with the classical nuclear localization sequence (cNLS) from nucleoplasmin (Npl), which was submitted to multiple 100 ns of MD simulations. Representative conformations were selected for calculating the 87 lowest frequencies NMs of vibration, and a displacement approach was applied along each NM. Based on geometric criteria, using the radius of curvature and inter-repeat angles as the reference metrics, the main motions of Impα were described. Moreover, we determined the salt bridges, hydrogen bonds and hydrophobic interactions in the Impα-NplNLS interface. Our results show the bending and twisting motions participating in the recognition of nuclear proteins, allowing the accommodation and adjustment of a classical bipartite NLS sequence. The essential contacts for the nuclear import were also described and were mostly in agreement with previous studies, suggesting that the residues in the cNLS linker region establish important contacts with Impα adjusting the cNLS backbone. The MD simulations combined with NM analysis can be applied to the Impα-NLS system to help understand interactions between Impα and cNLSs and the analysis of non-classic NLSs. Laborat�rio de Bioinform�tica e Biof�sica Computacional Departamento de F�sica e Biof�sica Instituto de Bioci�ncias de Botucatu UNESP - Universidade Estadual Paulista Instituto de Biotecnologia (IBTEC) UNESP - Universidade Estadual Paulista Laborat�rio de Biologia Computacional e Bioinform�tica Centro de Ci�ncias Naturais e Humanas UFABC - Universidade Federal Do ABC Laborat�rio de Bioinform�tica e Biof�sica Computacional Departamento de F�sica e Biof�sica Instituto de Bioci�ncias de Botucatu UNESP - Universidade Estadual Paulista Instituto de Biotecnologia (IBTEC) UNESP - Universidade Estadual Paulista

Maturity-onset diabetes of the young type 3 (MODY3) is a non-ketotic form of diabetes associated with poor insulin secretion. Over the past years, several studies have reported the association of missense mutations in the Hepatocyte Nuclear Factor 1 Alpha (HNF1A) with MODY3. Missense mutations in the POU homeodomain (POUH) of HNF1A hinder binding to the DNA, thereby leading to a dysfunctional protein. Missense mutations of the HNF1A were retrieved from public databases and subjected to a three-step computational mutational analysis to identify the underlying mechanism. First, the pathogenicity and stability of the mutations were analyzed to determine whether they alter protein structure and function. Second, the sequence conservation and DNA-binding sites of the mutant positions were assessed; as HNF1A protein is a transcription factor. Finally, the biochemical properties of the biological system were validated using molecular dynamic simulations in Gromacs 4.6.3 package. Two arginine residues (131 and 203) in the HNF1A protein are highly conserved residues and contribute to the function of the protein. Furthermore, the R131W, R131Q, and R203C mutations were predicted to be highly deleterious by in silico tools and showed lower binding affinity with DNA when compared to the native protein using the molecular docking analysis. Triplicate runs of molecular dynamic (MD) simulations (50ns) revealed smaller changes in patterns of deviation, fluctuation, and compactness, in complexes containing the R131Q and R131W mutations, compared to complexes containing the R203C mutant complex. We observed reduction in the number of intermolecular hydrogen bonds, compactness, and electrostatic potential, as well as the loss of salt bridges, in the R203C mutant complex. Substitution of arginine with cysteine at position 203 decreases the affinity of the protein for DNA, thereby destabilizing the protein. Based on our current findings, the MD approach is an important tool for elucidating the impact and affinity of mutations in DNA-protein interactions and understanding their function. Qatar University grant QUUG-CAS-DHS-14/15-3.

The heat-tolerance mechanisms of (hyper)thermophilic proteins provide a unique opportunity to investigate the unsolved protein folding problem. In an attempt to determine whether the interval between residues in sequence might play a role in determining thermostability, we constructed a sequence interval-dependent value function to calculate the residue pair frequency. Additionally, we identified a new sequence arrangement pattern, where like-charged residues tend to be adjacently assembled, while unlike-charged residues are distributed over longer intervals, using statistical analysis of a large sequence database. This finding indicated that increasing the intervals between unlike-charged residues can increase protein thermostability, with the arrangement patterns of these charged residues serving as thermodynamically favorable nucleation points for protein folding. Additionally, we identified that the residue pairs K-E, R-E, L-V and V-V involving long sequence intervals play important roles involving increased protein thermostability. This work demonstrated a novel approach for considering sequence intervals as keys to understanding protein folding. Our findings of novel relationships between residue arrangement and protein thermostability can be used in industry and academia to aid the design of thermostable proteins.

Protein design is a useful method to create novel artificial proteins. A rational approach to design a heterodimeric protein using domain swapping for horse myoglobin (Mb) was developed. As confirmed by X-ray crystallographic analysis, a heterodimeric Mb with two different active sites was produced efficiently from two surface mutants of Mb, in which the charges of two amino acids involved in the dimer salt bridges were reversed in each mutant individually, with the active site of one mutant modified. This study shows that the method of constructing heterodimeric Mb with domain swapping is useful for designing artificial multiheme proteins.

The 20 residue Trp-cage mini-protein is one of smallest proteins that adopt a stable folded structure containing also well-defined secondary structure elements. The hydrophobic core is arranged around a single central Trp residue. Despite several experimental and simulation studies the detailed folding mechanism of the Trp-cage protein is still not completely understood. Starting from fully extended as well as from partially folded Trp-cage structures a series of molecular dynamics simulations in explicit solvent and using four different force fields was performed. All simulations resulted in rapid collapse of the protein to on average relatively compact states. The simulations indicate a significant dependence of the speed of folding to near-native states on the side chain rotamer state of the central Trp residue. Whereas the majority of intermediate start structures with the central Trp side chain in a near-native rotameric state folded successfully within less than 100 ns only a fraction of start structures reached near-native folded states with an initially non-native Trp side chain rotamer state. Weak restraining of the Trp side chain dihedral angles to the state in the folded protein resulted in significant acceleration of the folding both starting from fully extended or intermediate conformations. The results indicate that the side chain conformation of the central Trp residue can create a significant barrier for controlling transitions to a near native folded structure. Similar mechanisms might be of importance for the folding of other protein structures.

The relative promiscuity of hub proteins such as postsynaptic density protein-95 (PSD-95) can be achieved by alternative splicing, allosteric regulation, and post-translational modifications, the latter of which is the most efficient method of accelerating cellular responses to environmental changes in vivo. Here, a mutational approach was used to determine the impact of phosphorylation and succinimidation post-translational modifications on the binding affinity of the postsynaptic density protein-95/discs large/zonula occludens-1 (PDZ3) domain of PSD-95. Molecular dynamics simulations revealed that the binding affinity of this domain is influenced by an interplay between salt-bridges linking the α3 helix, the β2–β3 loop and the positively charged Lys residues in its high-affinity hexapeptide ligand KKETAV. The α3 helix is an extra structural element that is not present in other PDZ domains, which links PDZ3 with the following SH3 domain in the PSD-95 protein. This regulatory mechanism was confirmed experimentally via thermodynamic and NMR chemical shift perturbation analyses, discarding intra-domain long-range effects. Taken together, the results presented here reveal the molecular basis of the regulatory role of the α3 extra-element and the effects of post-translational modifications of PDZ3 on its binding affinity, both energetically and dynamically., This research was supported by grants CVI-05915, from the Andalusian Regional Government (http://www.juntadeandalucia.es), BIO2009-13261-C02 and BIO2012-39922-C02, from the Spanish Ministry of Science and Innovation (http://www.idi.mineco.gob.es/portal/site​/MICINN/) and FEDER. JMC received a postdoctoral contract from the Spanish Ministry of Science and Innovation. CCV was a recipient of a Formación de Personal Investigador fellowship from the Spanish Ministry of Science and Innovation.

Site-directed mutagenesis was used to assess the contribution of individual residues and a bound calcium in the 55-69 region of the thermolysin-like protease of Bacillus stearothermophilus (TLP-ste) to thermal stability. The importance of the 55-69 region was reflected by finding that almost all mutations had drastic effects on stability. These effects (both stabilizing and destabilizing) were obtained by mutations affecting main chain flexibility, as well as by mutations affecting the interaction between the 55-69 region and the rest of the protease molecule, The calcium-dependency of stability could be largely abolished by mutating one of its ligands (Asp57 or Asp59). In the case of the Asp57-->Ser mutation the accompanying loss in stability was modest compared with the effects of other destabilizing mutations or the effects of (combinations of) stabilizing mutations, The detailed knowledge of the stability-determining region of TLP-ste permits effective rational design of stabilizing mutations, which, presumably, are also useful for related TLP such as thermolysin, This is demonstrated by the successful design of a stabilizing salt bridge involving residues 65 and 11.

Site-directed mutagenesis was used to assess the contribution of individual residues and a bound calcium in the 55-69 region of the thermolysin-like protease of Bacillus stearothermophilus (TLP-ste) to thermal stability. The importance of the 55-69 region was reflected by finding that almost all mutations had drastic effects on stability. These effects (both stabilizing and destabilizing) were obtained by mutations affecting main chain flexibility, as well as by mutations affecting the interaction between the 55-69 region and the rest of the protease molecule. The calcium-dependency of stability could be largely abolished by mutating one of its ligands (Asp57 or Asp59). In the case of the Asp57-->Ser mutation, the accompanying loss in stability was modest compared with the effects of other destabilizing mutations or the effects of (combinations of) stabilizing mutations. The detailed knowledge of the stability-determining region of TLP-ste permits effective rational design of stabilizing mutations, which, presumably, are also useful for related TLP such as thermolysin. This is demonstrated by the successful design of a stabilizing salt bridge involving residues 65 and 11.

By applying REMD simulations we have performed comparative analysis of the conformational ensembles of amino-truncated Aβ10-40 peptide produced with five force fields, which combine four protein parameterizations (CHARMM36, CHARMM22*, CHARMM22/cmap, and OPLS-AA) and two water models (standard and modified TIP3P). Aβ10-40 conformations were analyzed by computing secondary structure, backbone fluctuations, tertiary interactions, and radius of gyration. We have also calculated Aβ10-40 3JHNHα-coupling and RDC constants and compared them with their experimental counterparts obtained for the full-length Aβ1-40 peptide. Our study led us to several conclusions. First, all force fields predict that Aβ adopts unfolded structure dominated by turn and random coil conformations. Second, specific TIP3P water model does not dramatically affect secondary or tertiary Aβ10-40 structure, albeit standard TIP3P model favors slightly more compact states. Third, although the secondary structures observed in CHARMM36 and CHARMM22/cmap simulations are qualitatively similar, their tertiary interactions show little consistency. Fourth, two force fields, OPLS-AA and CHARMM22* have unique features setting them apart from CHARMM36 or CHARMM22/cmap. OPLS-AA reveals moderate β-structure propensity coupled with extensive, but weak long-range tertiary interactions leading to Aβ collapsed conformations. CHARMM22* exhibits moderate helix propensity and generates multiple exceptionally stable long- and short-range interactions. Our investigation suggests that among all force fields CHARMM22* differs the most from CHARMM36. Fifth, the analysis of 3JHNHα-coupling and RDC constants based on CHARMM36 force field with standard TIP3P model led us to an unexpected finding that in silico Aβ10-40 and experimental Aβ1-40 constants are generally in better agreement than these quantities computed and measured for identical peptides, such as Aβ1-40 or Aβ1-42. This observation suggests that the differences in the conformational ensembles of Aβ10-40 and Aβ1-40 are small and the former can be used as proxy of the full-length peptide. Based on this argument, we concluded that CHARMM36 force field with standard TIP3P model produces the most accurate representation of Aβ10-40 conformational ensemble., Author Summary Dependence of protein conformational ensembles on force field parameterizations limits the predictive power of molecular dynamics simulations. To address this problem, we evaluated five all-atom force fields for their consistency in reproducing the conformational ensemble of Alzheimer’s Aβ10-40 peptide. To generate conformational ensembles, we have used replica exchange molecular dynamics and computed Aβ10-40 secondary and tertiary structures. We found that, although all force fields predict Aβ10-40 unfolded structure, they strongly disagree on helix and β propensities and tertiary structure distributions. We have also calculated Aβ10-40 J-coupling and residual dipolar coupling constants and compared them with the experimental data for the full-length Aβ1-40 peptide. Unexpectedly, we determined that in silico Aβ10-40 and experimental Aβ1-40 constants are in better agreement than these quantities computed and measured previously for identical peptides, such as Aβ1-40 or Aβ1-42. We then concluded that the conformational ensembles of Aβ10-40 and Aβ1-40 are similar and on this basis argue that CHARMM36 force field with standard TIP3P water model provides the most accurate description of Aβ10-40. Although our objective was not to evaluate the biomolecular force fields in general, our study is expected to facilitate their proper selection for the simulations of Alzheimer’s peptides.

The flexible hydrophobic ligand binding pocket (LBP) of estrogen receptor α (ERα) allows the binding of a wide variety of endocrine disruptors. Upon ligand binding, the LBP reshapes around the contours of the ligand and stabilizes the complex by complementary hydrophobic interactions and specific hydrogen bonds with the ligand. Here we present a framework for quantitative analysis of the steric and electronic features of the human ERα-ligand complex using three dimensional (3D) protein-ligand interaction description combined with 3D-QSAR approach. An empirical hydrophobicity density field is applied to account for hydrophobic contacts of ligand within the LBP. The obtained 3D-QSAR model revealed that hydrophobic contacts primarily determine binding affinity and govern binding mode with hydrogen bonds. Several residues of the LBP appear to be quite flexible and adopt a spectrum of conformations in various ERα-ligand complexes, in particular His524. The 3D-QSAR was combined with molecular docking based on three receptor conformations to accommodate receptor flexibility. The model indicates that the dynamic character of the LBP allows accommodation and stable binding of structurally diverse ligands, and proper representation of the protein flexibility is critical for reasonable description of binding of the ligands. Our results provide a quantitative and mechanistic understanding of binding affinity and mode of ERα agonists and antagonists that may be applicable to other nuclear receptors.

International audience; Innate immune receptors NOD1 and NOD2 are activated by bacterial peptidoglycans leading to recruitment of adaptor kinase RIP2, which, upon phosphorylation and ubiquitination, becomes a scaffold for downstream effectors. The kinase domain (RIP2K) is a pharmaceutical target for inflammatory diseases caused by aberrant NOD2-RIP2 signalling. Although structures of active RIP2K in complex with inhibitors have been reported, the mechanism of RIP2K activation remains to be elucidated. Here we analyse RIP2K activation by combining crystal structures of the active and inactive states with mass spectrometric characterization of their phosphorylation profiles. The active state has Helix αC inwardly displaced and the phosphorylated Activation Segment (AS) disordered, whilst in the inactive state Helix αC is outwardly displaced and packed against the helical, non-phosphorylated AS. Biophysical measurements show that the active state is a stable dimer whilst the inactive kinase is in a monomer-dimer equilibrium, consistent with the observed structural differences at the dimer interface. We conclude that RIP2 kinase auto-phosphorylation is intimately coupled to dimerization, similar to the case of BRAF. Our results will help drug design efforts targeting RIP2 as a potential treatment for NOD2-RIP2 related inflammatory diseases.

Eukaryotic protein kinases regulate most cellular functions by phosphorylating targeted protein substrates through a highly conserved catalytic core. In the active state, the catalytic core oscillates between open, intermediate, and closed conformations. Currently, the intramolecular interactions that regulate the active state mechanics are not well understood. Here, using cAMP-dependent protein kinase as a representative model coupled with biochemical, biophysical, and computational techniques, we define a set of highly conserved electrostatic and hydrophobic interactions working harmoniously to regulate these mechanics. These include the previously identified salt bridge between a lysine from the β3-strand and a glutamate from the αC-helix as well as an electrostatic interaction between the phosphorylated activation loop and αC-helix and an ensemble of hydrophobic residues of the Regulatory spine and Shell. Moreover, for over three decades it was thought that the highly conserved β3-lysine was essential for phosphoryl transfer, but our findings show that the β3-lysine is not required for phosphoryl transfer but is essential for the active state mechanics., Author Summary Eukaryotic protein kinases (EPKs) regulate over a third of the human proteome by transferring the γ-phosphate from adenosine triphosphate (ATP) to a protein substrate in a process known as protein phosphorylation. Biochemical and biophysical studies have shown that EPKs undergo multiconformational rearrangements in which the catalytic core is oscillating between open, intermediate, and closed conformations when active. Presently, the intramolecular interactions that regulate this dynamic process are not well understood. In this paper, we show how a set of conserved electrostatic and hydrophobic interactions harmoniously regulate the active state mechanics. The electrostatic interactions involve the highly conserved salt bridge between the lysine from subdomain-II and glutamate from subdomain-III as well as an interaction between the activation loop and αC-helix. The hydrophobic interactions include the nonlinear motifs known as the Regulatory spine and Shell that traverse both lobes of the catalytic core. Furthermore, our findings show that the highly conserved “catalytic lysine” is not directly required for phosphoryl transfer but rather serves as a hub that aligns and positions the dynamic core elements required for catalysis.

Jak-family tyrosine kinases mediate signaling from diverse cytokine receptors. Binding of Jaks to their cognate receptors is mediated by their N-terminal region, which contains FERM and SH2 domains. Here we describe the crystal structure of the FERM-SH2 region of Jak2 at 3.0Å resolution. The structure reveals that these domains and their flanking linker segments interact intimately to form an integrated structural module. The Jak2 FERM-SH2 structure closely resembles that recently described for Tyk2, another member of the Jak family. While the overall architecture and interdomain orientations are preserved between Jak2 and Tyk2, we identify residues in the putative receptor-binding groove that differ between the two and may contribute to the specificity of receptor recognition. Analysis of Jak mutations that are reported to disrupt receptor binding reveals that they lie in the hydrophobic core of the FERM domain, and are thus expected to compromise the structural integrity of the FERM-SH2 unit. Similarly, analysis of mutations in Jak3 that are associated with severe combined immunodeficiency suggests that they compromise Jak3 function by destabilizing the FERM-SH2 structure.

The high-resolution crystal structure of the leucine transporter (LeuT) is frequently used as a template for homology models of the dopamine transporter (DAT). Although similar in structure, DAT differs considerably from LeuT in a number of ways: (i) when compared to LeuT, DAT has very long intracellular amino and carboxyl termini; (ii) LeuT and DAT share a rather low overall sequence identity (22%) and (iii) the extracellular loop 2 (EL2) of DAT is substantially longer than that of LeuT. Extracellular zinc binds to DAT and restricts the transporter‚s movement through the conformational cycle, thereby resulting in a decrease in substrate uptake. Residue H293 in EL2 praticipates in zinc binding and must be modelled correctly to allow for a full understanding of its effects. We exploited the high-affinity zinc binding site endogenously present in DAT to create a model of the complete transmemberane domain of DAT. The zinc binding site provided a DAT-specific molecular ruler for calibration of the model. Our DAT model places EL2 at the transporter lipid interface in the vicinity of the zinc binding site. Based on the model, D206 was predicted to represent a fourth co-ordinating residue, in addition to the three previously described zinc binding residues H193, H375 and E396. This prediction was confirmed by mutagenesis: substitution of D206 by lysine and cysteine affected the inhibitory potency of zinc and the maximum inhibition exerted by zinc, respectively. Conversely, the structural changes observed in the model allowed for rationalizing the zinc-dependent regulation of DAT: upon binding, zinc stabilizes the outward-facing state, because its first coordination shell can only be completed in this conformation. Thus, the model provides a validated solution to the long extracellular loop and may be useful to address other aspects of the transport cycle., Author Summary The dopamine transporter (DAT) regulates dopaminergic neurotransmission in the brain and is implicated in numerous human disease states. DAT is unique among the monoamine neurotransmitter transporter family because its substrate transport is inhibited by extracellular zinc. DAT homology models rely upon the crystal structure of LeuT solved in 2005. LeuT and DAT share a relatively low overall sequence identity of 22%. In addition, the length of the second extracellular loop of DAT exceeds that of LeuT by 21 residues. The zinc binding site cannot be directly modeled from the LeuT template alone because of these differences. Current available homology models of DAT focused on substrate or inhibitor binding rather than on the second extracellular loop. We exploited the specificity of the zinc binding site to build and calibrate a DAT homology model of the complete transmembrane domain. Our model predicted that the zinc binding site in DAT consists of four zinc co-ordinating residues rather than three that had been previously identified. We verified this hypothesis by site-directed mutagenesis and uptake inhibition studies.

Acidianus two-tailed virus (ATV) infects crenarchaea of the genus Acidianus living in terrestrial thermal springs at extremely high temperatures and low pH. ATV is a member of the Bicaudaviridae virus family and undergoes extra-cellular development of two tails, a process that is unique in the viral world. To understand this intriguing phenomenon, we have undertaken structural studies of ATV virion proteins and here we present the crystal structure of one of these proteins, ATV[Formula: see text]. ATV[Formula: see text] forms tetramers in solution and a molecular envelope is provided for the tetramer, computed from small-angle X-ray scattering (SAXS) data. The crystal structure has properties typical of hyperthermostable proteins, including a relatively high number of salt bridges. However, the protein also exhibits flexible loops and surface pockets. Remarkably, ATV[Formula: see text] displays a new [Formula: see text] protein fold, consistent with the absence of homologues of this protein in public sequence databases.

Here we report a comprehensive analysis through alanine-scanning mutagenesis of the contribution of surface ion pairs to the thermal stability of Alicyclobacillus acidocaldarius esterase 2 (EST2). We produced 16 single mutants, 4 double mutants corresponding to selected ion pairs R31/E118, E43/K102, R58/D130, D145/R148, 2 double mutants (R63A/R98A and E50A/D94A) involving residues of a large ion network on the protein surface and the double-mutant R98A/R148A meant to disrupt the R98 interactions within the said network and, contextually, the interaction between R148 and D145. The double-mutant E43A/E273K was obtained by chance. All selected residues were replaced with alanine except E91, which was mutated to a glycine and K102, which was changed to a glutamine. All 24 proteins were over-expressed in Escherichia coli, purified and characterized with respect to the main features. Structural stability data were compared with an in silico prediction of ΔΔG values. Our study of the individual factors involved in thermostability and their structural interpretation reveals that the great stability of this thermophilic protein can be explained by the contribution of a few residues at the protein surface.

The distinguishing property of Sm protein associations is their high stability. In order to understand this property, we analyzed the interface non-covalent interactions and compared the properties of the Sm protein interfaces with those of a test set, Binding Interface Database (BID). The comparison revealed that the main differences between interfaces of Sm proteins and those of the BID set are the content of charged residues, hydrogen bonds, salt bridges, and conservation scores of interface residues. In Sm proteins, the interfaces have more hydrophobic and fewer charged residues than the surface, which is also the case for the BID test set and other proteins. However, in the interfaces, the content of charged residues in Sm proteins (26%) is substantially larger than that in the BID set (22%). Both interfaces of Sm proteins and of test set have a similar number of hydrophobic interactions per 100 angstrom(2). The interfaces of Sm proteins have substantially more hydrogen bonds than the interfaces in test set. The results show clearly that the interfaces of Sm proteins form more salt bridges compared with test set. On average, there are about 16 salt bridges per interface. The high conservation score of amino acids that are involved in non-covalent interactions in protein interfaces is an additional strong argument for their importance. The overriding conclusion from this study is that the non-covalent interactions in Sm protein interfaces considerably contribute to stability of higher order structures. (C) 2010 Elsevier Ltd. All rights reserved.

Sso7d is a small protein especially attractive as a model for characterizing by NMR, the structural determinants responsible for thermal stability. With this aim, structural parameters of the Sso7d protein have been compared with those of single point mutants, in particular with K12L, more stable (T-m around 100 degrees C) and F31A, less stable (T-m 61 degrees C) than the wild-type. In addition, a mutant lacking eight residues of the C-terminal helix (T-m 50 degrees C) has been studied to investigate the role of the helix in structure stabilization. Melting temperatures, Delta Delta G's of unfolding, cavities, electrostatic interactions, solvent accessible areas, hydrophobic cores, and aromatic cluster geometries have been analyzed in order to gain insights into the key structural factors determining the thermostability differences. Our data suggest that absence of cavities and a larger salt bridge network represent the major determinants contributing to the enhanced stability of K12L relative to the Sso7d protein and its mutant F31A.

Experimental evidence suggests that the folding and aggregation of the amyloid $\beta$-protein (A$\beta$) into oligomers is a key pathogenetic event in Alzheimer's disease (AD). Inhibiting the pathologic folding and oligomerization of A$\beta$ could be effective in the prevention and treatment of AD. Here, using all-atom molecular dynamics simulations in explicit solvent, we probe the initial stages of folding of a decapeptide segment of A$\beta$, A$\beta_{21-30}$, shown experimentally to nucleate the folding process. In addition, we examine the folding of a homologous decapeptide containing an amino acid substitution linked to hereditary cerebral hemorrhage with amyloidosis--Dutch type, [Gln22]A$\beta_{21-30}$. We find that: (i) when the decapeptide is in water, hydrophobic interactions and transient salt bridges between Lys28 and either Glu22 or Asp23 are important in the formation of a loop in the Val24--Lys28 region of the wild type decapeptide; (ii) in the presence of salt ions, salt bridges play a more prominent role in the stabilization of the loop; (iii) in water with a reduced density, the decapeptide forms a helix, indicating the sensitivity of folding to different aqueous environments; (iv) the ``Dutch'' peptide in water, in contrast to the wild type peptide, fails to form a long-lived Val24--Lys28 loop, suggesting that loop stability is a critical factor in determining whether A$\beta$ folds into pathologic structures. Our results are relevant to understand the mechanism of A$\beta$ peptide folding in different environments, such as intra- and extracellular milieus or cell membranes, and how amino acid substitutions linked to familial forms of amyloidosis cause disease., Comment: 3 figures, 2 tables

Das Hüllglykoprotein Hämagglutinin (HA) von Influenzavirus ist verantwortlich sowohl für die Bindung als auch für die nachfolgende Fusion der viralen Hülle mit der endosomalen Membran. Eine Analyse der 3D Struktur der HA-Ektodomaine zeigt, dass die Stabilität des Proteins sowohl durch kovalente als auch durch nicht-kovalente Wechselwirkungen bedingt ist. Die Konformationsänderung von HA bei saurem pH-Wert weißt auf eine mögliche Rolle von Protonierungseffekten auf ionisierbare Aminosäuren hin. Untersuchungen zur Bedeutung geladener Aminosäuren und Salzbrücken für die Struktur des HA wurden auf der Grundlage von ‚site directed mutagenesis’ durchgeführt. Der Einfluss der Mutationen auf die Konformationsänderung und die Fusionsaktivität von HA wurden durch einen Proteinase K-Assay bzw. Fluoreszenzmikroskopie erfasst. Die Ergebnisse beider Methoden wurden miteinander korreliert. Abgesehen von der Mutante R109E zeigten Wildtyp-HA und alle anderen Mutanten eine vergleichbare Oberflächenexpression. Die beobachteten Unterschiede in der pH-Abhängigkeit der Konformationumwandlung zwischen Wildtyp-HA und HA-Mutanten zeigen, daß eine Zerstörung von Salzbrücken und ggf. eine Erhöhung der elektrostatischen Abstoßung an den betrachteten Kontakstellen sehr wahrscheinlich eine Herabsetzung der energetischen Barriere der Konformationsumwandlung verursacht. Dieser Ergebnisse erklären die molekularen Grundlagen des erhöhten pH-Schwellwertes der HA-Konformationsumwandlung von Amantadin-resistenten Influenzaviren. Im Gegensatz wurde für Mutanten, die die Stabilität von HA erhöhten, keine Konformationsumwandlung bei einem pH-Wert beobachtet, der typisch für die Konformationumwandlung von Wildtyp-HA war. Aminosäuren, die denen dieser stabiliserenden Mutationen entsprachen, wurden in einer natürlichen Influenzavirusvariante – A/JPN/305/57 – gefunden. Die Bedeutung von Ladungen für die Stabilität der HA-Ektodomaine wird dadurch unterstrichen, dass eine Konservierung einer positiven Ladung und insbesondere eines Argininrestes in der Position 109 (Nummerierung auf der Basis von HA X31) für alle Influenzaviren A und B gefunden wurde. Die Ergebnisse der Arbeit zeigen, dass sehr wahrscheinlich eine komplexe Salzbrücke an der Kontaktfläche zwischen HA1 und HA2 für alle Influenzaviren A evolutionär konserviert ist. Hemagglutinin (HA), a major envelope glycoprotein is responsible for fusing viral and endosomal membranes during influenza virus entry. The analysis of 3D crystal structure of the HA ectodomain shows that the stability of protein is maintained by both non-covalent and covalent interactions. The conformational change of HA at low pH indicates a role for protonation effects of the ionisable amino acids. Structural investigations were done using “site directed mutagenesis” in order to conceive the importance of charged amino acids and more emphatically the involvement of salt bridges. The effect of mutations on the conformational change and fusion activity was probed by proteinase K assay and fluorescence microscopy respectively. It was observed that HA-wt and all the mutants except R109E showed comparable surface expression. The difference in pH threshold between the HA-wt and the mutants showed that breakage of salt bridge and further incorporation of repulsion at the considered interfaces would lower the energy barrier requirements for the conformational change. The results explain the molecular basis of the higher pH threshold for naturally occurring amantadine resistant mutants. On the other hand, mutants designed to stabilise the HA were resistant to conformational changes at those pH values which typically trigger the conformational change of HA-wt. Coincidentally these mutations were found to be existing in the natural variant of H2 Japan subtype (A/JPN/305/57). Interestingly, the study shows that a positive charge and, more specifically, an Arg residue at position 109 (numbering based on X-31 strain) is conserved in all of the influenza A and B viruses underlining the relevance of electrostatic interactions for the HA stability. Aptly a complex salt bridge at the interface of HA1 and HA2 is probably conserved evolutionarily in all the members of influenza A virus.

International audience

SCOPUS: ar.j, FLWIN, info:eu-repo/semantics/published

UNLABELLED: Salt bridges can play important roles in protein structure and function and have stabilizing and destabilizing effects in protein folding. ESBRI is a software available as web tool which analyses the salt bridges in a protein structure, starting from the atomic coordinates. In the case of protein complexes, the salt bridges between protein chains can be evaluated, as well as those among specific charged amino acids and the different protein subunits, in order to obtain useful information regard the protein-protein interaction. AVAILABILITY: The service is available at the URL: http://bioinformatica.isa.cnr.it/ESBRI/