Flexibility and mobility in mesophilic and thermophilic homologous proteins from molecular dynamics and FoldUnfold method

In silico analysis of alpha-synuclein protein variants and posttranslational modifications related to Parkinson's disease

Parkinson's disease (PD) is among the most prevalent neurodegenerative disorders, affecting over 10 million people worldwide. The protein encoded by the SNCA gene, alpha-synuclein (ASYN), is the major component of Lewy body (LB) aggregates, a histopathological hallmark of PD. Mutations and posttranslational modifications (PTMs) in ASYN are known to influence protein aggregation and LB formation, possibly playing a crucial role in PD pathogenesis. In this work, we applied computational methods to characterize the effects of missense mutations and PTMs on the structure and function of ASYN. Missense mutations in ASYN were compiled from the literature/databases and underwent a comprehensive predictive analysis. Phosphorylation and SUMOylation sites of ASYN were retrieved from databases and predicted by algorithms. ConSurf was used to estimate the evolutionary conservation of ASYN amino acids. Molecular dynamics (MD) simulations of ASYN wild-type and variants A30G, A30P, A53T, and G51D were performed using the GROMACS package. Seventy-seven missense mutations in ASYN were compiled. Although most mutations were not predicted to affect ASYN stability, aggregation propensity, amyloid formation, and chaperone binding, the analyzed mutations received relatively high rates of deleterious predictions and predominantly occurred at evolutionarily conserved sites within the protein. Moreover, our predictive analyses suggested that the following mutations may be possibly harmful to ASYN and, consequently, potential targets for future investigation: K6N, T22I, K34E, G36R, G36S, V37F, L38P, G41D, and K102E. The MD analyses pointed to remarkable flexibility and essential dynamics alterations at nearly all domains of the studied variants, which could lead to impaired contact between NAC and the C-terminal domain triggering protein aggregation. These alterations may have functional implications for ASYN and provide important insight into the molecular mechanism of PD, supporting the design of future biomedical research and improvements in existing therapies for the disease.

Recent Advances in β-Glucosidase Sequence and Structure Engineering: A Brief Review

β-glucosidases (BGLs) play a crucial role in the degradation of lignocellulosic biomass as well as in industrial applications such as pharmaceuticals, foods, and flavors. However, the application of BGLs has been largely hindered by issues such as low enzyme activity, product inhibition, low stability, etc. Many approaches have been developed to engineer BGLs to improve these enzymatic characteristics to facilitate industrial production. In this article, we review the recent advances in BGL engineering in the field, including the efforts from our laboratory. We summarize and discuss the BGL engineering studies according to the targeted functions as well as the specific strategies used for BGL engineering.

Protein stability and dynamics influenced by ligands in extremophilic complexes - a molecular dynamics investigation

In this study, we explore the structural and dynamic adaptations of the Tryptophan synthase α-subunit in a ligand bound state in psychrophilic, mesophilic and hyperthermophilic organisms at different temperatures by MD simulations. We quantify the global and local fluctuations in the 40 ns time scale by analyzing the root mean square deviation/fluctuations. The distinct behavior of the active site and loop 6 is observed with the elevation of temperature. Protein stability relies more on electrostatic interactions, and these interactions might be responsible for the stability of varying temperature evolved proteins. The paper also focuses on the effect of temperature on protein dynamics and stability governed by the distinct behavior of the ligand associated with its retention, binding and dissociation over the course of time. The integration of principle component analysis and a free energy landscape was useful in identifying the conformational space accessible to ligand bound homologues and how the presence of the ligand alters the conformational and dynamic properties of the protein.

Protein flexibility in the light of structural alphabets

Protein structures are valuable tools to understand protein function. Nonetheless, proteins are often considered as rigid macromolecules while their structures exhibit specific flexibility, which is essential to complete their functions. Analyses of protein structures and dynamics are often performed with a simplified three-state description, i.e., the classical secondary structures. More precise and complete description of protein backbone conformation can be obtained using libraries of small protein fragments that are able to approximate every part of protein structures. These libraries, called structural alphabets (SAs), have been widely used in structure analysis field, from definition of ligand binding sites to superimposition of protein structures. SAs are also well suited to analyze the dynamics of protein structures. Here, we review innovative approaches that investigate protein flexibility based on SAs description. Coupled to various sources of experimental data (e.g., B-factor) and computational methodology (e.g., Molecular Dynamic simulation), SAs turn out to be powerful tools to analyze protein dynamics, e.g., to examine allosteric mechanisms in large set of structures in complexes, to identify order/disorder transition. SAs were also shown to be quite efficient to predict protein flexibility from amino-acid sequence. Finally, in this review, we exemplify the interest of SAs for studying flexibility with different cases of proteins implicated in pathologies and diseases.

Vibrational entropy differences between mesophile and thermophile proteins and their use in protein engineering

We recently introduced ENCoM, an elastic network atomic contact model, as the first coarse-grained normal mode analysis method that accounts for the nature of amino acids and can predict the effect of mutations on thermostability based on changes vibrational entropy. In this proof-of-concept article, we use pairs of mesophile and thermophile homolog proteins with identical structures to determine if a measure of vibrational entropy based on normal mode analysis can discriminate thermophile from mesophile proteins. We observe that in around 60% of cases, thermophile proteins are more rigid at equivalent temperatures than their mesophile counterpart and this difference can guide the design of proteins to increase their thermostability through series of mutations. We observe that mutations separating thermophile proteins from their mesophile orthologs contribute independently to a decrease in vibrational entropy and discuss the application and implications of this methodology to protein engineering.

Stability and rigidity/flexibility-two sides of the same coin?

Protein molecules require both flexibility and rigidity for functioning. The fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. We have determined flexible regions for four homologous pairs from thermophilic and mesophilic organisms by two methods: the fast FoldUnfold which uses amino acid sequence and the time consuming MDFirst which uses three-dimensional structures. We demonstrate that both methods allow determining flexible regions in protein structure. For three of the four thermophile-mesophile pairs of proteins, FoldUnfold predicts practically the same flexible regions which have been found by the MD/First method. As expected, molecular dynamics simulations show that thermophilic proteins are more rigid in comparison to their mesophilic homologues. Analysis of rigid clusters and their decomposition provides new insights into protein stability. It has been found that the local networks of salt bridges and hydrogen bonds in thermophiles render their structure more stable with respect to fluctuations of individual contacts. Such network includes salt bridge triads Agr-Glu-Lys and Arg-Glu-Arg, or salt bridges (such as Arg-Glu) connected with hydrogen bonds. This ionic network connects alpha helices and rigidifies the structure. Mesophiles can be characterized by stand alone salt bridges and hydrogen bonds or small ionic clusters. Such difference in the network of salt bridges results in different flexibility of homologous proteins. Combining both approaches allows characterizing structural features in atomic detail that determine the rigidity/flexibility of a protein structure. This article is a part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

Regulation of protein-protein binding by coupling between phosphorylation and intrinsic disorder: analysis of human protein complexes

Phosphorylation offers a dynamic way to regulate protein activity, subcellular localization, and stability. The majority of signaling pathways involve an extensive set of protein-protein interactions, and phosphorylation is widely used to regulate protein-protein binding by affecting the stability, kinetics and specificity of interactions. Previously it was found that phosphorylation sites tend to be located on protein-protein binding interfaces and may orthosterically modulate the strength of interactions. Here we studied the effect of phosphorylation on protein binding in relation to intrinsic disorder for different types of human protein complexes with known structure of the binding interface. Our results suggest that the processes of phosphorylation, binding and disorder-order transitions are coupled to each other, with about one quarter of all disordered interface Ser/Thr/Tyr sites being phosphorylated. Namely, residue site disorder and interfacial states significantly affect the phosphorylation of serine and to a lesser extent of threonine. Tyrosine phosphorylation might not be directly associated with binding through disorder, and is often observed in ordered interface regions which are not predicted to be disordered in the unbound state. We analyze possible mechanisms of how phosphorylation might regulate protein-protein binding via intrinsic disorder, and specifically focus on how phosphorylation could prevent disorder-order transitions upon binding.

Structure fluctuations and conformational changes in protein binding

Structure fluctuations and conformational changes accompany all biological processes involving macromolecules. The paper presents a classification of protein residues based on the normalized equilibrium fluctuations of the residue centers of mass in proteins and a statistical analysis of conformation changes in the side-chains upon binding. Normal mode analysis and an elastic network model were applied to a set of protein complexes to calculate the residue fluctuations and develop the residue classification. Comparison with a classification based on normalized B-factors suggests that the B-factors may underestimate protein flexibility in solvent. Our classification shows that protein loops and disordered fragments are enriched with highly fluctuating residues and depleted with weakly fluctuating residues. Strategies for engineering thermostable proteins are discussed. To calculate the dihedral angles distribution functions, the configuration space was divided into cells by a cubic grid. The effect of protein association on the distribution functions depends on the amino acid type and a grid step in the dihedral angles space. The changes in the dihedral angles increase from the near-backbone dihedral angle to the most distant one, for most residues. On average, one fifth of the interface residues change the rotamer state upon binding, whereas the rest of the interface residues undergo local readjustments within the same rotamer.

PredyFlexy: flexibility and local structure prediction from sequence

Protein structures are necessary for understanding protein function at a molecular level. Dynamics and flexibility of protein structures are also key elements of protein function. So, we have proposed to look at protein flexibility using novel methods: (i) using a structural alphabet and (ii) combining classical X-ray B-factor data and molecular dynamics simulations. First, we established a library composed of structural prototypes (LSPs) to describe protein structure by a limited set of recurring local structures. We developed a prediction method that proposes structural candidates in terms of LSPs and predict protein flexibility along a given sequence. Second, we examine flexibility according to two different descriptors: X-ray B-factors considered as good indicators of flexibility and the root mean square fluctuations, based on molecular dynamics simulations. We then define three flexibility classes and propose a method based on the LSP prediction method for predicting flexibility along the sequence. This method does not resort to sophisticate learning of flexibility but predicts flexibility from average flexibility of predicted local structures. The method is implemented in PredyFlexy web server. Results are similar to those obtained with the most recent, cutting-edge methods based on direct learning of flexibility data conducted with sophisticated algorithms. PredyFlexy can be accessed at http://www.dsimb.inserm.fr/dsimb_tools/predyflexy/.