Unique Structure and Dynamics of the EphA5 Ligand Binding Domain Mediate Its Binding Specificity as Revealed by X-ray Crystallography, NMR and MD Simulations

10.1371/journal.pone.0074040PLoS ONE89-POLN

Ethical Issues in Engineering Models: An Operations Researcher’s Reflections

This article starts with an overview of the author’s personal involvement—as an Operations Research consultant—in several engineering case-studies that may raise ethical questions; e.g., case-studies on nuclear waste, water management, sustainable ecology, military tactics, and animal welfare. All these case studies employ computer simulation models. In general, models are meant to solve practical problems, which may have ethical implications for the various stakeholders; namely, the modelers, the clients, and the public at large. The article further presents an overview of codes of ethics in a variety of disciples. It discusses the role of mathematical models, focusing on the validation of these models’ assumptions. Documentation of these model assumptions needs special attention. Some ethical norms and values may be quantified through the model’s multiple performance measures, which might be optimized. The uncertainty about the validity of the model leads to risk or uncertainty analysis and to a search for robust models. Ethical questions may be pressing in military models, including war games. However, computer games and the related experimental economics may also provide a special tool to study ethical issues. Finally, the article briefly discusses whistleblowing. Its many references to publications and websites enable further study of ethical issues in modeling

Accelerated Profile HMM Searches

Profile hidden Markov models (profile HMMs) and probabilistic inference methods have made important contributions to the theory of sequence database homology search. However, practical use of profile HMM methods has been hindered by the computational expense of existing software implementations. Here I describe an acceleration heuristic for profile HMMs, the “multiple segment Viterbi” (MSV) algorithm. The MSV algorithm computes an optimal sum of multiple ungapped local alignment segments using a striped vector-parallel approach previously described for fast Smith/Waterman alignment. MSV scores follow the same statistical distribution as gapped optimal local alignment scores, allowing rapid evaluation of significance of an MSV score and thus facilitating its use as a heuristic filter. I also describe a 20-fold acceleration of the standard profile HMM Forward/Backward algorithms using a method I call “sparse rescaling”. These methods are assembled in a pipeline in which high-scoring MSV hits are passed on for reanalysis with the full HMM Forward/Backward algorithm. This accelerated pipeline is implemented in the freely available HMMER3 software package. Performance benchmarks show that the use of the heuristic MSV filter sacrifices negligible sensitivity compared to unaccelerated profile HMM searches. HMMER3 is substantially more sensitive and 100- to 1000-fold faster than HMMER2. HMMER3 is now about as fast as BLAST for protein searches

Stability mechanisms of a thermophilic laccase probed by molecular dynamics.

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T opt (∼350 K) from 300 K, this change correlated with a beginning "unzipping" of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F(-) (but not Cl(-)) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F(-) intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes

Computational Fragment-Based Binding Site Identification by Ligand Competitive Saturation

Fragment-based drug discovery using NMR and x-ray crystallographic methods has proven utility but also non-trivial time, materials, and labor costs. Current computational fragment-based approaches circumvent these issues but suffer from limited representations of protein flexibility and solvation effects, leading to difficulties with rigorous ranking of fragment affinities. To overcome these limitations we describe an explicit solvent all-atom molecular dynamics methodology (SILCS: Site Identification by Ligand Competitive Saturation) that uses small aliphatic and aromatic molecules plus water molecules to map the affinity pattern of a protein for hydrophobic groups, aromatic groups, hydrogen bond donors, and hydrogen bond acceptors. By simultaneously incorporating ligands representative of all these functionalities, the method is an in silico free energy-based competition assay that generates three-dimensional probability maps of fragment binding (FragMaps) indicating favorable fragment∶protein interactions. Applied to the two-fold symmetric oncoprotein BCL-6, the SILCS method yields two-fold symmetric FragMaps that recapitulate the crystallographic binding modes of the SMRT and BCOR peptides. These FragMaps account both for important sequence and structure differences in the C-terminal halves of the two peptides and also the high mobility of the BCL-6 His116 sidechain in the peptide-binding groove. Such SILCS FragMaps can be used to qualitatively inform the design of small-molecule inhibitors or as scoring grids for high-throughput in silico docking that incorporate both an atomic-level description of solvation and protein flexibility

A practical guide to the simultaneous determination of protein structure and dynamics using metainference

Accurate protein structural ensembles can be determined with metainference, a Bayesian inference method that integrates experimental information with prior knowledge of the system and deals with all sources of uncertainty and errors as well as with system heterogeneity. Furthermore, metainference can be implemented using the metadynamics approach, which enables the computational study of complex biological systems requiring extensive conformational sampling. In this chapter, we provide a step-by-step guide to perform and analyse metadynamic metainference simulations using the ISDB module of the open-source PLUMED library, as well as a series of practical tips to avoid common mistakes. Specifically, we will guide the reader in the process of learning how to model the structural ensemble of a small disordered peptide by combining state-of-the-art molecular mechanics force fields with nuclear magnetic resonance data, including chemical shifts, scalar couplings and residual dipolar couplings.Comment: 49 pages, 9 figure

How Thioredoxin Dissociates Its Mixed Disulfide

The dissociation mechanism of the thioredoxin (Trx) mixed disulfide complexes is unknown and has been debated for more than twenty years. Specifically, opposing arguments for the activation of the nucleophilic cysteine as a thiolate during the dissociation of the complex have been put forward. As a key model, the complex between Trx and its endogenous substrate, arsenate reductase (ArsC), was used. In this structure, a Cys29Trx-Cys89ArsC intermediate disulfide is formed by the nucleophilic attack of Cys29Trx on the exposed Cys82ArsC-Cys89ArsC in oxidized ArsC. With theoretical reactivity analysis, molecular dynamics simulations, and biochemical complex formation experiments with Cys-mutants, Trx mixed disulfide dissociation was studied. We observed that the conformational changes around the intermediate disulfide bring Cys32Trx in contact with Cys29Trx. Cys32Trx is activated for its nucleophilic attack by hydrogen bonds, and Cys32Trx is found to be more reactive than Cys82ArsC. Additionally, Cys32Trx directs its nucleophilic attack on the more susceptible Cys29Trx and not on Cys89ArsC. This multidisciplinary approach provides fresh insights into a universal thiol/disulfide exchange reaction mechanism that results in reduced substrate and oxidized Trx

Global profiling of co- and post-translationally N-myristoylated proteomes in human cells

Protein N-myristoylation is a ubiquitous co- and post-translational modification that has been implicated in the development and progression of a range of human diseases. Here, we report the global N-myristoylated proteome in human cells determined using quantitative chemical proteomics combined with potent and specific human N-myristoyltransferase (NMT) inhibition. Global quantification of N-myristoylation during normal growth or apoptosis allowed the identification of >100 N-myristoylated proteins, >95% of which are identified for the first time at endogenous levels. Furthermore, quantitative dose response for inhibition of N-myristoylation is determined for >70 substrates simultaneously across the proteome. Small-molecule inhibition through a conserved substrate-binding pocket is also demonstrated by solving the crystal structures of inhibitor-bound NMT1 and NMT2. The presented data substantially expand the known repertoire of co- and post-translational N-myristoylation in addition to validating tools for the pharmacological inhibition of NMT in living cells

Molecular Evolution of the Primate Antiviral Restriction Factor Tetherin

Background: Tetherin is a recently identified antiviral restriction factor that restricts HIV-1 particle release in the absence of the HIV-1 viral protein U (Vpu). It is reminiscent of APOBEC3G and TRIM5a that also antagonize HIV. APOBEC3G and TRIM5a have been demonstrated to evolve under pervasive positive selection throughout primate evolution, supporting the redqueen hypothesis. Therefore, one naturally presumes that Tetherin also evolves under pervasive positive selection throughout primate evolution and supports the red-queen hypothesis. Here, we performed a detailed evolutionary analysis to address this presumption. Methodology/Principal Findings: Results of non-synonymous and synonymous substitution rates reveal that Tetherin as a whole experiences neutral evolution rather than pervasive positive selection throughout primate evolution, as well as in non-primate mammal evolution. Sliding-window analyses show that the regions of the primate Tetherin that interact with viral proteins are under positive selection or relaxed purifying selection. In particular, the sites identified under positive selection generally focus on these regions, indicating that the main selective pressure acting on the primate Tetherin comes from virus infection. The branch-site model detected positive selection acting on the ancestral branch of the New World Monkey lineage, suggesting an episodic adaptive evolution. The positive selection was also found in duplicated Tetherins in ruminants. Moreover, there is no bias in the alterations of amino acids in the evolution of the primate Tetherin, implyin

Structural, Stability, Dynamic and Binding Properties of the ALS-Causing T46I Mutant of the hVAPB MSP Domain as Revealed by NMR and MD Simulations

T46I is the second mutation on the hVAPB MSP domain which was recently identified from non-Brazilian kindred to cause a familial amyotrophic lateral sclerosis (ALS). Here using CD, NMR and molecular dynamics (MD) simulations, we characterized the structure, stability, dynamics and binding capacity of the T46I-MSP domain. The results reveal: 1) unlike P56S which we previously showed to completely eliminate the native MSP structure, T46I leads to no significant disruption of the native secondary and tertiary structures, as evidenced from its far-UV CD spectrum, as well as Cα and Cβ NMR chemical shifts. 2) Nevertheless, T46I does result in a reduced thermodynamic stability and loss of the cooperative urea-unfolding transition. As such, the T46I-MSP domain is more prone to aggregation than WT at high protein concentrations and temperatures in vitro, which may become more severe in the crowded cellular environments. 3) T46I only causes a 3-fold affinity reduction to the Nir2 peptide, but a significant elimination of its binding to EphA4. 4) EphA4 and Nir2 peptide appear to have overlapped binding interfaces on the MSP domain, which strongly implies that two signaling networks may have a functional interplay in vivo. 5) As explored by both H/D exchange and MD simulations, the MSP domain is very dynamic, with most loop residues and many residues on secondary structures highly fluctuated or/and exposed to bulk solvent. Although T46I does not alter overall dynamics, it does trigger increased dynamics of several local regions of the MSP domain which are implicated in binding to EphA4 and Nir2 peptide. Our study provides the structural and dynamic understanding of the T46I-causing ALS; and strongly highlights the possibility that the interplay of two signaling networks mediated by the FFAT-containing proteins and Eph receptors may play a key role in ALS pathogenesis