Thermodynamic and structural characterization of Asn and Ala residues in the disallowed II' region of the Ramachandran plot

Combining Ramachandran plot and molecular dynamics simulation for structural-based variant classification: Using TP53 variants as model

The wide application of new DNA sequencing technologies is generating vast quantities of genetic variation data at unprecedented speed. Developing methodologies to decode the pathogenicity of the variants is imperatively demanding. We hypothesized that as deleterious variants may function through disturbing structural stability of their affected proteins, information from structural change caused by genetic variants can be used to identify the variants with deleterious effects. In order to measure the structural change for proteins with large size, we designed a method named RP-MDS composed of Ramachandran plot (RP) and Molecular Dynamics Simulation (MDS). Ramachandran plot captures the variant-caused secondary structural change, whereas MDS provides a quantitative measure for the variant-caused globular structural change. We tested the method using variants in TP53 DNA binding domain of 219 residues as the model. In total, RP-MDS identified 23 of 38 (60.5%) TP53 known Pathogenic variants and 17 of 42 (41%) TP53 VUS that caused significant changes of P53 structure. Our study demonstrates that RP-MDS method provides a powerful protein structure-based tool to screen deleterious genetic variants affecting large-size proteins.

Bioinformatics analysis of four proteins of Leishmania donovani to guide epitopes vaccine design and drug targets selection

Visceral leishmaniasis (VL) is a serious and widespread parasitic disease caused by Leishmania donovani complex. The threat of this fatal disease continues due to the lack of ideal drugs or vaccines. In this study, we selected Amastin, CaNA2, Kmp-11 and PDI proteins of Leishmania donovani for study, which are VL vaccine candidates or possible drug targets. Eleven bioinformatics tools were used to analyze different aspects of these proteins, including amino acid composition, topology, signal peptide, secondary structure, surface properties, phosphorylation sites and kinases, protein binding sites, 3D homology modeling, B cell epitopes, MHC class Ⅰ and Ⅱ epitopes and protein-protein interactions. Finally, the functionally related amino acid sites and dominant epitopes of these proteins were founded. Some possible relationships between protein structure, phosphorylation sites, protein binding sites and epitopes were also discovered. High flexibility and random coils regions of protein have a tendency to be phosphorylated, bind proteins and present epitopes. Since some phosphorylation sites and their kinases are involved in Leishmania invasion and survival in host cells, they may be potential drug targets. Bioinformatics analysis helps us better understand protein function and find dominant epitopes to guide drug design and vaccine development.

Crystallographic studies on protein misfolding: Domain swapping and amyloid formation in the SH3 domain

Oligomerization by 3D domain swapping is found in a variety of proteins of diverse size, fold and function. In the early 1960s this phenomenon was postulated for the oligomers of ribonuclease A, but it was not until the 1990s that X-ray diffraction provided the first experimental evidence of this special manner of oligomerization. Nowadays, structural information has allowed the identification of these swapped oligomers in over one hundred proteins. Although the functional relevance of this phenomenon is not clear, this alternative folding of protomers into intertwined oligomers has been related to amyloid formation. Studies on proteins that develop 3D domain swapping might provide some clues on the early stages of amyloid formation. The SH3 domain is a small modular domain that has been used as a model to study the basis of protein folding. Among SH3 domains, the c-Src-SH3 domain emerges as a helpful model to study 3D domain swapping and amyloid formation.

The First Step of Amyloidogenic Aggregation

The structural and dynamic characterization of the on-pathway intermediates involved in the mechanism of amyloid fibril formation is one of the major remaining biomedical challenges of our time. In addition to mature fibrils, various oligomeric structures are implicated in both the rate-limiting step of the nucleation process and the neuronal toxicity of amyloid deposition. Single-molecule fluorescence spectroscopy (SMFS) is an excellent tool for extracting most of the relevant information on these molecular systems, especially advanced multiparameter approaches, such as pulsed interleaved excitation (PIE). In our investigations of an amyloidogenic SH3 domain of α-spectrin, we have found dynamic oligomerization, even prior to incubation. Our single-molecule PIE experiments revealed that these species are small, mostly dimeric, and exhibit a loose and dynamic molecular organization. Furthermore, these experiments have allowed us to obtain quantitative information regarding the oligomer stability. These pre-amyloidogenic oligomers may potentially serve as the first target for fibrillization-prevention strategies.

Structure of an intermediate state in protein folding and aggregation

Protein-folding intermediates have been implicated in amyloid fibril formation involved in neurodegenerative disorders. However, the structural mechanisms by which intermediates initiate fibrillar aggregation have remained largely elusive. To gain insight, we used relaxation dispersion nuclear magnetic resonance spectroscopy to determine the structure of a low-populated, on-pathway folding intermediate of the A39V/N53P/V55L (A, Ala; V, Val; N, Asn; P, Pro; L, Leu) Fyn SH3 domain. The carboxyl terminus remains disordered in this intermediate, thereby exposing the aggregation-prone amino-terminal β strand. Accordingly, mutants lacking the carboxyl terminus and thus mimicking the intermediate fail to safeguard the folding route and spontaneously form fibrillar aggregates. The structure provides a detailed characterization of the non-native interactions stabilizing an aggregation-prone intermediate under native conditions and insight into how such an intermediate can derail folding and initiate fibrillation.

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain-domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain-domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.

Conformational and functional analysis of molecular dynamics trajectories by self-organising maps

BACKGROUND

Molecular dynamics (MD) simulations are powerful tools to investigate the conformational dynamics of proteins that is often a critical element of their function. Identification of functionally relevant conformations is generally done clustering the large ensemble of structures that are generated. Recently, Self-Organising Maps (SOMs) were reported performing more accurately and providing more consistent results than traditional clustering algorithms in various data mining problems. We present a novel strategy to analyse and compare conformational ensembles of protein domains using a two-level approach that combines SOMs and hierarchical clustering.

RESULTS

The conformational dynamics of the α-spectrin SH3 protein domain and six single mutants were analysed by MD simulations. The Cα's Cartesian coordinates of conformations sampled in the essential space were used as input data vectors for SOM training, then complete linkage clustering was performed on the SOM prototype vectors. A specific protocol to optimize a SOM for structural ensembles was proposed: the optimal SOM was selected by means of a Taguchi experimental design plan applied to different data sets, and the optimal sampling rate of the MD trajectory was selected. The proposed two-level approach was applied to single trajectories of the SH3 domain independently as well as to groups of them at the same time. The results demonstrated the potential of this approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space in a simple 2D visualisation, as well as of effectively highlighting differences in the conformational dynamics directly related to biological functions.

CONCLUSIONS

The use of a two-level approach combining SOMs and hierarchical clustering for conformational analysis of structural ensembles of proteins was proposed. It can easily be extended to other study cases and to conformational ensembles from other sources.

High-resolution structure of an alpha-spectrin SH3-domain mutant with a redesigned hydrophobic core

The alpha-spectrin SH3 domain (Spc-SH3) is a small modular domain which has been broadly used as a model protein in folding studies and these studies have sometimes been supported by structural information obtained from the coordinates of Spc-SH3 mutants. The structure of B5/D48G, a multiple mutant designed to improve the hydrophobic core and as a consequence the protein stability, has been solved at 1 A resolution. The crystals belonged to the orthorhombic space group P2(1)2(1)2(1), with unit-cell parameters a=24.79, b=37.23, c=62.95 A. This mutant also bears a D48G substitution in the distal loop and this mutation has also been reported to increase the stability of the protein by itself. The structure of the B5/D48G mutant shows a highly packed hydrophobic core and a more ordered distal loop compared with previous Spc-SH3 structures.

Comparison of solid-state dipolar couplings and solution relaxation data provides insight into protein backbone dynamics

Analyses of solution (15)N relaxation data and solid-state (1)H(N)-(15)N dipolar couplings from a small globular protein, alpha-spectrin SH3 domain, produce a surprisingly similar pattern of order parameters. This result suggests that there is little or no ns-mus dynamics throughout most of the sequence and, in particular, in the structured portion of the backbone. At the same time, evidence of ns-mus motions is found in the flexible loops and termini. These findings, corroborated by the MD simulations of alpha-spectrin SH3 in a hydrated crystalline environment and in solution, are consistent with the picture of protein dynamics that has recently emerged from the solution studies employing residual dipolar couplings.

The in vivo and in vitro aggregation properties of globular proteins correlate with their conformational stability: the SH3 case

Protein misfolding and deposition underlie an increasing number of debilitating human disorders and constitute a problem of major concern in biotechnology. In the last years, in vitro studies have provided valuable insights into the physicochemical principles underlying protein aggregation. Nevertheless, information about the determinants of protein deposition within the cell is scarce and only a few systematic studies comparing in vitro and in vivo data have been reported. Here, we have used the SH3 domain of alpha-spectrin as a model globular protein in an attempt to understand the relationship between protein aggregation in the test-tube and in the more complex cellular environment. The investigation of the aggregation in Escherichia coli of this domain and a large set of mutants, together with the analysis of their sequential and conformational properties allowed us to evaluate the contribution of different polypeptidic factors to the cellular deposition of globular proteins. The data presented here suggest that the rules that govern in vitro protein aggregation are also valid in in vivo contexts. They also provide relevant insights into intracellular protein deposition in both conformational diseases and recombinant protein production.

An error analysis for two-state protein-folding kinetic parameters and phi-values: progress toward precision by exploring pH dependencies on Leffler plots

The interpretation of phi-values has led to an understanding of the folding transition state ensemble of a variety of proteins. Although the main guidelines and equations for calculating phi are well established, there remains some controversy about the quality of the numerical values obtained. By analyzing a complete set of results from kinetic experiments with the SH3 domain of alphaspectrin (Spc-SH3) and applying classical error methods and error-propagation formulas, we evaluated the uncertainties involved in two-state-folding kinetic experimental parameters and the corresponding calculated phi-values. We show that kinetic constants in water and m values can be properly estimated from a judicious weighting of fitting errors and describe some procedures to calculate the errors in Gibbs energies and phi-values from a traditional two-point Leffler analysis. Furthermore, on the basis of general assumptions made with the protein engineering method, we show how to generate multipoint Leffler plots via the analysis of pH dependencies of kinetic parameters. We calculated the definitive phi-values for a collection of single mutations previously designed to characterize the folding transition state of the alphaspectrin SH3 domain. The effectiveness of the pH-scanning procedure is also discussed in the context of error analysis. Judging from the magnitudes of the error bars obtained from two-point and multipoint Leffler plots, we conclude that the precision obtained for phi-values should be approximately 25%, a reasonable limit that takes into account the propagation of experimental errors.

Increasing protein conformational stability by optimizing beta-turn sequence

Protein conformational stability is an important concern in many fields. Here, we describe a strategy for significantly increasing conformational stability by optimizing beta-turn sequence. Proline and glycine residues are statistically preferred at several beta-turn positions, presumably because their unique side-chains contribute favorably to conformational stability in certain beta-turn positions. However, beta-turn sequences often deviate from preferred proline or preferred glycine. Therefore, our strategy involves replacing non-proline and non-glycine beta-turn residues with preferred proline or preferred glycine residues. Here, we develop guidelines for selecting appropriate mutations, and present results for five mutations (S31P, S42G, S48P, T76P, and Q77G) that significantly increase the conformational stability of RNase Sa. The increases in stability ranged from 0.7 kcal/mol to 1.3 kcal/mol. The strategy was successful in overlapping or isolated beta-turns, at buried (up to 50%) or completely exposed sites, and at relatively flexible or inflexible sites. Considering the significant number of beta-turn residues in every globular protein and the frequent deviation of beta-turn sequences from preferred proline and preferred glycine residues, this simple, efficient strategy will be useful for increasing the conformational stability of proteins.

Cooperative propagation of local stability changes from low-stability and high-stability regions in a SH3 domain

Site-directed mutagenesis has been used to produce local stability changes at two regions of the binding site surface of the alpha-spectrin SH3 domain (Spc-SH3) differing in their intrinsic stability. Mutations were made at residue 56, located at the solvent-exposed side of the short 3(10) helix, and at residue 21 in the tip of the flexible RT-loop. NMR chemical-shift analysis and X-ray crystallography indicated negligible changes produced by the mutations in the native structure limited to subtle rearrangements near the mutated residue and at flexible loops. Additionally, mutations do not alter importantly the SH3 binding site structure, although produce significant changes in its affinity for a proline-rich decapeptide. The changes in global stability measured by differential scanning calorimetry are consistent the local energy changes predicted by theoretical models, with the most significant effects observed for the Ala-Gly mutations. Propagation of the local stability changes throughout the domain structure has been studied at a per-residue level of resolution by NMR-detected amide hydrogen-deuterium exchange (HX). Stability propagation is remarkably efficient in this small domain, apparently due to its intrinsically low stability. Nevertheless, the HX-core of the domain is not fully cooperative, indicating the existence of co-operative subunits within the core, which is markedly polarized. An equilibrium phi-analysis of the changes in the apparent Gibbs energies of HX per residue produced by the mutations has allowed us to characterize structurally the conformational states leading to HX. Some of these states resemble notably the folding transition state of the Spc-SH3 domain, suggesting a great potential of this approach to explore the folding energy landscape of proteins. An energy perturbation propagates more effectively from a flexible region to the core than in the opposite direction, because the former affects a broader region of the energy landscape than the latter. This might be of importance in understanding the special thermodynamic signature of the SH3-peptide interaction and the relevance of the dual character of SH3 binding sites.

Molecular Structure of a Novel Membrane Protease Specific for a Stomatin Homolog from the Hyperthermophilic Archaeon Pyrococcus horikoshii

Membrane-bound proteases are involved in various regulatory functions. Our previous study indicated that the N-terminal region of an open reading frame, PH1510 (residues 16-236, designated as 1510-N) from the hyperthermophilic archaeon Pyrococcus horikoshii, is a serine protease with a catalytic Ser-Lys dyad that specifically cleaves the C-terminal hydrophobic residues of a membrane protein, the stomatin-homolog PH1511. In humans, an absence of stomatin is associated with a form of hemolytic anemia known as hereditary stomatocytosis, but the function of stomatin is not fully understood. Here, we report the crystal structure of 1510-N in dimeric form. Each active site of 1510-N is rich in hydrophobic residues, which accounts for the substrate-specificity. The monomer of 1510-N shows structural similarity to one monomer of Escherichia coli ClpP, an ATP-dependent tetradecameric protease. But, their oligomeric forms are different. Major contributors to dimeric interaction in 1510-N are the alpha7 helix and beta9 strand, both of which are missing from ClpP. While the long handle region of ClpP contributes to the stacking of two heptameric rings, the corresponding L2 loop of 1510-N is disordered because the region has little interaction with other residues of the same molecule. The catalytic Ser97 of 1510-N is in almost the same location as the catalytic Ser97 of E.coli ClpP, whereas another residue, Lys138, presumably forming the catalytic dyad, is located in the disordered L2 region of 1510-N. These findings suggest that the binding of the substrate to the catalytic site of 1510-N induces conformational changes in a region that includes loop L2 so that Lys138 approaches the catalytic Ser97.

A single mutation induces amyloid aggregation in the alpha-spectrin SH3 domain: analysis of the early stages of fibril formation

The Src-homology region 3 domain of chicken alpha-spectrin (Spc-SH3) is a small two-state folding protein, which has never been described to form amyloid fibrils under any condition investigated so far. We show here that the mutation of asparagine 47 to alanine at the distal loop, which destabilises similarly the native and folding transition states of the domain, induces the formation of amyloid fibrils under mild acid conditions. Amyloid aggregation of the mutant is enhanced by the increase in temperature, protein concentration and NaCl concentration. The early stages of amyloid formation have been monitored as a function of time and temperature using a variety of biophysical methods. Differential scanning calorimetry experiments under conditions of amyloid formation have allowed the identification of different thermal transitions corresponding to conformational and aggregation processes as well as to the high-temperature disaggregation and unfolding of the amyloid fibrils. Aggregation is preceded by a rapid conformational change in the monomeric domain involving about 40% of the global unfolding enthalpy, considerable change in secondary structure, large loss of tertiary structure and exposure of hydrophobic patches to the solvent. The conformational change is followed by formation of a majority of oligomeric species with apparent hydrodynamic radius between 2.5 nm and 10nm, depending on temperature, together with the appearance and progressive growth of protofibrillar aggregates. After these early aggregation stages, long and curved fibrils of up to several micrometers start to develop by elongation of the protofibrils. The calorimetric data indicate that the specific enthalpy of fibril disaggregation and unfolding is relatively low, suggesting a low density of interactions within the fibril structure as compared to the native protein and a main entropy contribution to the stability of the amyloid fibrils.

Detection of dynamic water molecules in a microcrystalline sample of the SH3 domain of alpha-spectrin by MAS solid-state NMR

Water molecules are a major determinant of protein stability and are important for understanding protein-protein interactions. We present two experiments which allow to measure first the effective T(2) decay rate of individual amide proton, and second the magnetization build-up rates for a selective transfer from H(2)O to H(N) using spin diffusion as a mixing element. The experiments are demonstrated for a uniformly (2)H, (15)N labeled sample of a microcrystalline SH3 domain in which exchangeable deuterons were back-substituted with protons. In order to evaluate the NMR experimental data, as X-ray structure of the protein was determined using the same crystallization protocol as for the solid-state NMR sample. The NMR experimental data are correlated with the dipolar couplings calculated from H(2)O-H(N) distances which were extracted from the X-ray structure of the protein. We find that the H(N) T(2) decay rates and H(2)O-H(N) build-up rates are sensitive to distance and dynamics of the detected water molecules with respect to the protein. We show that qualitative information about localization and dynamics of internal water molecules can be obtained in the solid-state by interpretation of the spin dynamics of a reporter amide proton.

Folding specificity induced by loop stiffness

To test the importance of loop stiffness in restricting the heterogeneity of transition state ensemble, we relaxed the distal loop of 10 unstable redesigned hydrophobic core mutants of alpha-spectrin SH3 domain. This was achieved by replacing Asp48 by Gly at the tip of the distal hairpin. Although the change was local in nature, the effect on stabilization was not uniform across the core mutants tested. There is an inverse rough correlation between the stabilization and the increase in buried hydrophobic volume, with respect to the wild type. Interestingly enough, proteins that although unstable are properly folded become molten globule-like after relaxation of the distal loop. These results highlight the importance of stiffness in restricting the conformational heterogeneity of a protein during the folding reaction. An interplay between unspecific hydrophobic interactions and constraint induced by polar interactions, or in this case local stiffness, is essential to achieve a well-ordered folded structure.

Structural cooperativity in the SH3 domain studied by site-directed mutagenesis and amide hydrogen exchange

We have studied the effects produced by site-directed mutagenesis upon energetic and structural cooperativity in the Src homology region 3 domain of alpha-spectrin. The mutation of Asn47 to Gly or Ala in the distal loop brings about significant changes to the global stability of the domain in spite of not affecting its structure to any great extent. The binding affinity for a proline-rich peptide is also largely diminished in both mutant domains. We have compared the apparent Gibbs energies of the amide hydrogen-deuterium exchange (HX) between the wild-type and the Gly47 mutant. The observed changes in the Gibbs energy of HX indicate a remarkable energetic cooperativity in this small domain. Regions of the domain's core have a high cooperativity with the position of the mutation, indicating that their HX occurs mainly in states in which the distal loop is unstructured. More flexible regions, which undergo HX mainly by local motions, show a lower but still considerable cooperativity with the distal loop. We conclude that there is an important correlation between regional stability and cooperativity in this small domain.

Amino acid residue mutations uncouple cooperative effects in Escherichia coli D-3-phosphoglycerate dehydrogenase

d-3-Phosphoglycerate dehydrogenase from Escherichia coli contains two Gly-Gly sequences that occur at junctions between domains. A previous study (Grant, G. A., Xu, X. L., and Hu, Z. (2000) Biochemistry 39, 7316-7319) determined that the Gly-Gly sequence at the junction between the regulatory and substrate binding domain functions as a hinge between the domains. Mutations in this area significantly decrease the ability of serine to inhibit activity but have little effect on the K(m) and k(cat). Conversely, the present study shows that mutations to the Gly-Gly sequence at the junction of the substrate and nucleotide binding domains, which form the active site cleft, have a significant effect on the k(cat) of the enzyme without substantially altering the enzyme's sensitivity to serine. In addition, mutation of Gly-294, but not Gly-295, has a profound effect on the cooperativity of serine inhibition. Interestingly, even though cooperativity of inhibition can be reduced significantly, there is little apparent effect on the cooperativity of serine binding itself. An additional mutant, G336V,G337V, also reduces the cooperativity of inhibition, but in this case serine binding also is reduced to the point at which it cannot be measured by equilibrium dialysis. The double mutant G294V,G336V demonstrates that strain imposed by mutation at one hinge can be relieved partially by mutation at the other hinge, demonstrating linkage between the two hinge regions. These data show that the two cooperative processes, serine binding and catalytic inhibition, can be uncoupled. Consideration of the allowable torsional angles for the side chains introduced by the mutations yields a range of values for these angles that the glycine residues likely occupy in the native enzyme. A comparison of these values with the torsional angles found for the inhibited enzyme from crystal coordinates provides potential beginning and ending orientations for the transition from active to inhibited enzyme, which will allow modeling of the dynamics of domain movement.