Elena Papaleo

Free-energy landscape, principal component analysis, and structural clustering to identify representative conformations from molecular dynamics simulations: The myoglobin case

Several molecular dynamics (MD) simulations were used to sample conformations in the neighborhood of the native structure of holo-myoglobin (holo-Mb), collecting trajectories spanning 0.22 micros at 300 K. Principal component (PCA) and free-energy landscape (FEL) analyses, integrated by cluster analysis, which was performed considering the position and structures of the individual helices of the globin fold, were carried out. The coherence between the different structural clusters and the basins of the FEL, together with the convergence of parameters derived by PCA indicates that an accurate description of the Mb conformational space around the native state was achieved by multiple MD trajectories spanning at least 0.14 micros. The integration of FEL, PCA, and structural clustering was shown to be a very useful approach to gain an overall view of the conformational landscape accessible to a protein and to identify representative protein substates. This method could be also used to investigate the conformational and dynamical properties of Mb apo-, mutant, or delete versions, in which greater conformational variability is expected and, therefore identification of representative substates from the simulations is relevant to disclose structure-function relationship.

ETHE1 mutations are specific to ethylmalonic encephalopathy

Mutations in ETHE1, a gene located at chromosome 19q13, have recently been identified in patients affected by ethylmalonic encephalopathy (EE). EE is a devastating infantile metabolic disorder, characterised by widespread lesions in the brain, hyperlactic acidaemia, petechiae, orthostatic acrocyanosis, and high levels of ethylmalonic acid in body fluids. To investigate to what extent ETHE1 is responsible for EE, we analysed this gene in 29 patients with typical EE and in 11 patients presenting with early onset progressive encephalopathy with ethylmalonic aciduria (non-EE EMA). Frameshift, stop, splice site, and missense mutations of ETHE1 were detected in all the typical EE patients analysed. Western blot analysis of the ETHE1 protein indicated that some of the missense mutations are associated with the presence of the protein, suggesting that the corresponding wild type amino acid residues have a catalytic function. No ETHE1 mutations were identified in non-EE EMA patients. Experiments based on two dimensional blue native electrophoresis indicated that ETHE1 protein works as a supramolecular, presumably homodimeric, complex, and a three dimensional model of the protein suggests that it is likely to be a mitochondrial matrix thioesterase acting on a still unknown substrate. Finally, the 625G→A single nucleotide polymorphism in the gene encoding the short chain acyl-coenzyme A dehydrogenase (SCAD) was previously proposed as a co-factor in the aetiology of EE and other EMA syndromes. SNP analysis in our patients ruled out a pathogenic role of SCAD variants in EE, but did show a highly significant prevalence of the 625A alleles in non-EE EMA patients.

Protein flexibility in psychrophilic and mesophilic trypsins. Evidence of evolutionary conservation of protein dynamics in trypsin-like serine-proteases

Psychrophilic trypsins present fewer interdomain interactions and enhanced localized flexibility in regions close to the catalytic site. Notably, these regions fit well with the pattern of protein flexibility previously reported for psychrophilic elastases. Our results indicate that specific sites within the serine‐protease fold can be considered hot spots of cold‐adaptation and that psychrophilic trypsins and elastases have independently discovered similar molecular strategies to optimize flexibility at low temperatures.

Compaction Properties of an Intrinsically Disordered Protein: Sic1 and Its Kinase-Inhibitor Domain

IDPs in their unbound state can transiently acquire secondary and tertiary structure. Describing such intrinsic structure is important to understand the transition between free and bound state, leading to supramolecular complexes with physiological interactors. IDP structure is highly dynamic and, therefore, difficult to study by conventional techniques. This work focuses on conformational analysis of the KID fragment of the Sic1 protein, an IDP with a key regulatory role in the cell-cycle of Saccharomyces cerevisiae. FT-IR spectroscopy, ESI-MS, and IM measurements are used to capture dynamic and short-lived conformational states, probing both secondary and tertiary protein structure. The results indicate that the isolated Sic1 KID retains dynamic helical structure and populates collapsed states of different compactness. A metastable, highly compact species is detected. Comparison between the fragment and the full-length protein suggests that chain length is crucial to the stabilization of compact states of this IDP. The two proteins are compared by a length-independent compaction index.

Dynamic properties of extremophilic subtilisin-like serine-proteases

The investigation of the structural determinants of enzymatic temperature adaptation is a crucial pre-requisite both in terms of fundamental research and industrial applications to develop new biocatalysts active at different temperature ranges. In several cases, the differences related to cold- or warm-adaptation are related to subtle structural and aminoacidic differences at the molecular level, often hard to detect. In this context, we present a comparative study of psychrophilic, mesophilic and thermophilic subtilisin-like serine proteases by all-atom molecular dynamics (MD) simulations in explicit solvent using a multiple-replica approach. Our results strongly enforce the current view on localized flexibility in crucial functional regions for cold-adapted serine proteases and point out a different optimization and usage of salt-bridge interactions and networks in cold- and warm-adapted enzymes. The analyses allow to identify a subset of structural and dynamic features strictly associated to cold adaptation and which change from cold- to heat-active subtilisins. In particular, the thermophilic subtilisin presents a high affinity calcium binding site which is not structurally conserved in the mesophilic and psychrophilic counterparts, which, as it turns out from the MD analyses, at the same position show a stable salt bridge network and no stabilizing intra-molecular interactions, respectively. These aspects, along with differential flexibility in regions close to the active site or substrate binding pocket, can be an indication of evolution at this protein site toward a lower stability moving from high to low temperature conditions.

An acidic loop and cognate phosphorylation sites define a molecular switch that modulates ubiquitin charging activity in Cdc34-like enzymes.

E2 ubiquitin-conjugating enzymes are crucial mediators of protein ubiquitination, which strongly influence the ultimate fate of the target substrates. Recently, it has been shown that the activity of several enzymes of the ubiquitination pathway is finely tuned by phosphorylation, an ubiquitous mechanism for cellular regulation, which modulates protein conformation. In this contribution, we provide the first rationale, at the molecular level, of the regulatory mechanism mediated by casein kinase 2 (CK2) phosphorylation of E2 Cdc34-like enzymes. In particular, we identify two co-evolving signature elements in one of the larger families of E2 enzymes: an acidic insertion in β4α2 loop in the proximity of the catalytic cysteine and two conserved key serine residues within the catalytic domain, which are phosphorylated by CK2. Our investigations, using yeast Cdc34 as a model, through 2.5 µs molecular dynamics simulations and biochemical assays, define these two elements as an important phosphorylation-controlled switch that modulates opening and closing of the catalytic cleft. The mechanism relies on electrostatic repulsions between a conserved serine phosphorylated by CK2 and the acidic residues of the β4α2 loop, promoting E2 ubiquitin charging activity. Our investigation identifies a new and unexpected pivotal role for the acidic loop, providing the first evidence that this loop is crucial not only for downstream events related to ubiquitin chain assembly, but is also mandatory for the modulation of an upstream crucial step of the ubiquitin pathway: the ubiquitin charging in the E2 catalytic cleft.

xPyder: A PyMOL Plugin To Analyze Coupled Residues and Their Networks in Protein Structures.

A versatile method to directly identify and analyze short- or long-range coupled or communicating residues in a protein conformational ensemble is of extreme relevance to achieve a complete understanding of protein dynamics and structural communication routes. Here, we present xPyder, an interface between one of the most employed molecular graphics systems, PyMOL, and the analysis of dynamical cross-correlation matrices (DCCM). The approach can also be extended, in principle, to matrices including other indexes of communication propensity or intensity between protein residues, as well as the persistence of intra- or intermolecular interactions, such as those underlying protein dynamics. The xPyder plugin for PyMOL 1.4 and 1.5 is offered as Open Source software via the GPL v2 license, and it can be found, along with the installation package, the user guide, and examples, at http://linux.btbs.unimib.it/xpyder/.

Molecular Dynamics of Mesophilic-Like Mutants of a Cold-Adapted Enzyme: Insights into Distal Effects Induced by the Mutations

Networks and clusters of intramolecular interactions, as well as their “communication” across the three-dimensional architecture have a prominent role in determining protein stability and function. Special attention has been dedicated to their role in thermal adaptation. In the present contribution, seven previously experimentally characterized mutants of a cold-adapted α-amylase, featuring mesophilic-like behavior, have been investigated by multiple molecular dynamics simulations, essential dynamics and analyses of correlated motions and electrostatic interactions. Our data elucidate the molecular mechanisms underlying the ability of single and multiple mutations to globally modulate dynamic properties of the cold-adapted α-amylase, including both local and complex unpredictable distal effects. Our investigation also shows, in agreement with the experimental data, that the conversion of the cold-adapted enzyme in a warm-adapted variant cannot be completely achieved by the introduction of few mutations, also providing the rationale behind these effects. Moreover, pivotal residues, which are likely to mediate the effects induced by the mutations, have been identified from our analyses, as well as a group of suitable candidates for protein engineering. In fact, a subset of residues here identified (as an isoleucine, or networks of mesophilic-like salt bridges in the proximity of the catalytic site) should be considered, in experimental studies, to get a more efficient modification of the features of the cold-adapted enzyme.

Molecular Determinants of Enzyme Cold Adaptation: Comparative Structural and Computational Studies of Cold- and Warm-Adapted Enzymes

The identification of molecular mechanisms underlying enzyme cold adaptation is a hot-topic both for fundamental research and industrial applications. In the present contribution, we review the last decades of structural computational investigations on cold-adapted enzymes in comparison to their warm-adapted counterparts. Comparative sequence and structural studies allow the definition of a multitude of adaptation strategies. Different enzymes carried out diverse mechanisms to adapt to low temperatures, so that a general theory for enzyme cold adaptation cannot be formulated. However, some common features can be traced in dynamic and flexibility properties of these enzymes, as well as in their intra- and inter-molecular interaction networks. Interestingly, the current data suggest that a family-centered point of view is necessary in the comparative analyses of cold- and warm-adapted enzymes. In fact, enzymes belonging to the same family or superfamily, thus sharing at least the three-dimensional fold and common features of the functional sites, have evolved similar structural and dynamic patterns to overcome the detrimental effects of low temperatures.

Validation of protein models by a neural network approach

BackgroundThe development and improvement of reliable computational methods designed to evaluate the quality of protein models is relevant in the context of protein structure refinement, which has been recently identified as one of the bottlenecks limiting the quality and usefulness of protein structure prediction.ResultsIn this contribution, we present a computational method (Artificial Intelligence Decoys Evaluator: AIDE) which is able to consistently discriminate between correct and incorrect protein models. In particular, the method is based on neural networks that use as input 15 structural parameters, which include energy, solvent accessible surface, hydrophobic contacts and secondary structure content. The results obtained with AIDE on a set of decoy structures were evaluated using statistical indicators such as Pearson correlation coefficients, Znat, fraction enrichment, as well as ROC plots. It turned out that AIDE performances are comparable and often complementary to available state-of-the-art learning-based methods.ConclusionIn light of the results obtained with AIDE, as well as its comparison with available learning-based methods, it can be concluded that AIDE can be successfully used to evaluate the quality of protein structures. The use of AIDE in combination with other evaluation tools is expected to further enhance protein refinement efforts.