Molecular dynamics ensemble refinement of the heterogeneous native state of NCBD using chemical shifts and NOEs

Many proteins display complex dynamical properties that are often intimately linked to their biological functions. As the native state of a protein is best described as an ensemble of conformations, it is important to be able to generate models of native state ensembles with high accuracy. Due to limitations in sampling efficiency and force field accuracy it is, however, challenging to obtain accurate ensembles of protein conformations by the use of molecular simulations alone. Here we show that dynamic ensemble refinement, which combines an accurate atomistic force field with commonly available nuclear magnetic resonance (NMR) chemical shifts and NOEs, can provide a detailed and accurate description of the conformational ensemble of the native state of a highly dynamic protein. As both NOEs and chemical shifts are averaged on timescales up to milliseconds, the resulting ensembles reflect the structural heterogeneity that goes beyond that probed, e.g., by NMR relaxation order parameters. We selected the small protein domain NCBD as object of our study since this protein, which has been characterized experimentally in substantial detail, displays a rich and complex dynamical behaviour. In particular, the protein has been described as having a molten-globule like structure, but with a relatively rigid core. Our approach allowed us to describe the conformational dynamics of NCBD in solution, and to probe the structural heterogeneity resulting from both short- and long-timescale dynamics by the calculation of order parameters on different time scales. These results illustrate the usefulness of our approach since they show that NCBD is rather rigid on the nanosecond timescale, but interconverts within a broader ensemble on longer timescales, thus enabling the derivation of a coherent set of conclusions from various NMR experiments on this protein, which could otherwise appear in contradiction with each other

Amyloid-ß and a-Synuclein Decrease the Level of Metal-Catalyzed Reactive Oxygen Species by Radical Scavenging and Redox Silencing

The formation of reactive oxygen species (ROS) is linked to the pathogenesis of neurodegenerative diseases. Here we have investigated the effect of soluble and aggregated amyloid-ß (Aß) and a-synuclein (aS), associated with Alzheimer''s and Parkinson''s diseases, respectively, on the Cu2+-catalyzed formation of ROS in vitro in the presence of a biological reductant. We find that the levels of ROS, and the rate by which ROS is generated, are significantly reduced when Cu2+ is bound to Aß or aS, particularly when they are in their oligomeric or fibrillar forms. This effect is attributed to a combination of radical scavenging and redox silencing mechanisms. Our findings suggest that the increase in ROS associated with the accumulation of aggregated Aß or aS does not result from a particularly ROS-active form of these peptides, but rather from either a local increase of Cu2+ and other ROS-active metal ions in the aggregates or as a downstream consequence of the formation of the pathological amyloid structures

Biophysical and electrochemical studies of protein-nucleic acid interactions

This review is devoted to biophysical and electrochemical methods used for studying protein-nucleic acid (NA) interactions. The importance of NA structure and protein-NA recognition for essential cellular processes, such as replication or transcription, is discussed to provide background for description of a range of biophysical chemistry methods that are applied to study a wide scope of protein-DNA and protein-RNA complexes. These techniques employ different detection principles with specific advantages and limitations and are often combined as mutually complementary approaches to provide a complete description of the interactions. Electrochemical methods have proven to be of great utility in such studies because they provide sensitive measurements and can be combined with other approaches that facilitate the protein-NA interactions. Recent applications of electrochemical methods in studies of protein-NA interactions are discussed in detail

Local Cooperativity in an Amyloidogenic State of Human Lysozyme Observed at Atomic Resolution

The partial unfolding of human lysozyme underlies its conversion from the soluble state into amyloid fibrils observed in a fatal hereditary form of systemic amyloidosis. To understand the molecular origins of the disease, it is critical to characterize the structural and physicochemical properties of the amyloidogenic states of the protein. Here we provide a high-resolution view of the unfolding process at low pH for three different lysozyme variants, the wild-type protein and the mutants I56T and I59T, which show variable stabilities and propensities to aggregate in vitro. Using a range of biophysical techniques that includes differential scanning calorimetry and nuclear magnetic resonance spectroscopy, we demonstrate that thermal unfolding under amyloidogenic solution conditions involves a cooperative loss of native tertiary structure, followed by progressive unfolding of a compact, molten globule-like denatured state ensemble as the temperature is increased. The width of the temperature window over which the denatured ensemble progressively unfolds correlates with the relative amyloidogenicity and stability of these variants, and the region of lysozyme that unfolds first maps to that which forms the core of the amyloid fibrils formed under similar conditions. Together, these results present a coherent picture at atomic resolution of the initial events underlying amyloid formation by a globular protein

AA-Amyloidosis Can Be Transferred by Peripheral Blood Monocytes

Spongiform encephalopathies have been reported to be transmitted by blood transfusion even prior to the clinical onset. Experimental AA-amyloidosis shows similarities with prion disease and amyloid-containing organ-extracts can prime a recipient for the disease. In this systemic form of amyloidosis N-terminal fragments of the acute-phase reactant apolipoprotein serum amyloid A are the main amyloid protein. Initial amyloid deposits appear in the perifollicular region of the spleen, followed by deposits in the liver. We used the established murine model and induced AA-amyloidosis in NMRI mice by intravenous injections of purified amyloid fibrils (‘amyloid enhancing factor’) combined with inflammatory challenge (silver nitrate subcutaneously). Blood plasma and peripheral blood monocytes were isolated, sonicated and re-injected into new recipients followed by an inflammatory challenge during a three week period. When the animals were sacrificed presence of amyloid was analyzed in spleen sections after Congo red staining. Our result shows that some of the peripheral blood monocytes, isolated from animals with detectable amyloid, contained amyloid-seed that primed for AA-amyloid. The seeding material seems to have been phagocytosed by the cells since the AA-precursor (SAA1) was found not be expressed by the monocytes. Plasma recovered from mice with AA amyloidosis lacked seeding capacity. Amyloid enhancing activity can reside in monocytes recovered from mice with AA-amyloidosis and in a prion-like way trigger amyloid formation in conjunction with an inflammatory disorder. Human AA-amyloidosis resembles the murine form and every individual is expected to be exposed to conditions that initiate production of the acute-phase reactant. The monocyte-transfer mechanism should be eligible for the human disease and we point out blood transfusion as a putative route for transfer of amyloidosis

Crystallographic and Molecular Dynamics Analysis of Loop Motions Unmasking the Peptidoglycan-Binding Site in Stator Protein MotB of Flagellar Motor

Background: The C-terminal domain of MotB (MotB-C) shows high sequence similarity to outer membrane protein A and related peptidoglycan (PG)-binding proteins. It is believed to anchor the power-generating MotA/MotB stator unit of the bacterial flagellar motor to the peptidoglycan layer of the cell wall. We previously reported the first crystal structure of this domain and made a puzzling observation that all conserved residues that are thought to be essential for PG recognition are buried and inaccessible in the crystal structure. In this study, we tested a hypothesis that peptidoglycan binding is preceded by, or accompanied by, some structural reorganization that exposes the key conserved residues. Methodology/Principal Findings: We determined the structure of a new crystalline form (Form B) of Helicobacter pylori MotB-C. Comparisons with the existing Form A revealed conformational variations in the petal-like loops around the carbohydrate binding site near one end of the b-sheet. These variations are thought to reflect natural flexibility at this site required for insertion into the peptidoglycan mesh. In order to understand the nature of this flexibility we have performed molecular dynamics simulations of the MotB-C dimer. The results are consistent with the crystallographic data and provide evidence that the three loops move in a concerted fashion, exposing conserved MotB residues that have previously been implicated in binding of the peptide moiety of peptidoglycan. Conclusion/Significance: Our structural analysis provides a new insight into the mechanism by which MotB inserts into th

A Common Model for Cytokine Receptor Activation: Combined Scissor-Like Rotation and Self-Rotation of Receptor Dimer Induced by Class I Cytokine

The precise mechanism by which the binding of a class I cytokine to the extracellular domain of its corresponding receptor transmits a signal through the cell membrane remains unclear. Receptor activation involves a cytokine-receptor complex with a 1∶2 stoichiometry. Previously we used our transient-complex theory to calculate the rate constant of the initial cytokine-receptor binding to form a 1∶1 complex. Here we computed the binding pathway leading to the 1∶2 activation complex. Three cytokine systems (growth hormone, erythropoietin, and prolactin) were studied, and the focus was on the binding of the extracellular domain of the second receptor molecule after forming the 1∶1 complex. According to the transient-complex theory, translational and rotation diffusion of the binding entities bring them together to form a transient complex, which has near-native relative separation and orientation but not the short-range specific native interactions. Subsequently conformational rearrangement leads to the formation of the native complex. We found that the changes in relative orientations between the two receptor molecules from the transient complex to the 1∶2 native complex are similar for the three cytokine-receptor systems. We thus propose a common model for receptor activation by class I cytokines, involving combined scissor-like rotation and self-rotation of the two receptor molecules. Both types of rotations seem essential: the scissor-like rotation separates the intracellular domains of the two receptor molecules to make room for the associated Janus kinase molecules, while the self-rotation allows them to orient properly for transphosphorylation. This activation model explains a host of experimental observations. The transient-complex based approach presented here may provide a strategy for designing antagonists and prove useful for elucidating activation mechanisms of other receptors

Farseer-NMR: automatic treatment, analysis and plotting of large, multi-variable NMR data

We present Farseer-NMR (https://git.io/vAueU), a software package to treat, evaluate and combine NMR spectroscopic data from sets of protein-derived peaklists covering a range of experimental conditions. The combined advances in NMR and molecular biology enable the study of complex biomolecular systems such as flexible proteins or large multibody complexes, which display a strong and functionally relevant response to their environmental conditions, e.g. the presence of ligands, site-directed mutations, post translational modifications, molecular crowders or the chemical composition of the solution. These advances have created a growing need to analyse those systems’ responses to multiple variables. The combined analysis of NMR peaklists from large and multivariable datasets has become a new bottleneck in the NMR analysis pipeline, whereby information-rich NMR-derived parameters have to be manually generated, which can be tedious, repetitive and prone to human error, or even unfeasible for very large datasets. There is a persistent gap in the development and distribution of software focused on peaklist treatment, analysis and representation, and specifically able to handle large multivariable datasets, which are becoming more commonplace. In this regard, Farseer-NMR aims to close this longstanding gap in the automated NMR user pipeline and, altogether, reduce the time burden of analysis of large sets of peaklists from days/weeks to seconds/minutes. We have implemented some of the most common, as well as new, routines for calculation of NMR parameters and several publication-quality plotting templates to improve NMR data representation. Farseer-NMR has been written entirely in Python and its modular code base enables facile extension

Non-native hydrophobic interactions detected in unfolded apoflavodoxin by paramagnetic relaxation enhancement

Transient structures in unfolded proteins are important in elucidating the molecular details of initiation of protein folding. Recently, native and non-native secondary structure have been discovered in unfolded A. vinelandii flavodoxin. These structured elements transiently interact and subsequently form the ordered core of an off-pathway folding intermediate, which is extensively formed during folding of this α–β parallel protein. Here, site-directed spin-labelling and paramagnetic relaxation enhancement are used to investigate long-range interactions in unfolded apoflavodoxin. For this purpose, glutamine-48, which resides in a non-native α-helix of unfolded apoflavodoxin, is replaced by cysteine. This replacement enables covalent attachment of nitroxide spin-labels MTSL and CMTSL. Substitution of Gln-48 by Cys-48 destabilises native apoflavodoxin and reduces flexibility of the ordered regions in unfolded apoflavodoxin in 3.4 M GuHCl, because of increased hydrophobic interactions in the unfolded protein. Here, we report that in the study of the conformational and dynamic properties of unfolded proteins interpretation of spin-label data can be complicated. The covalently attached spin-label to Cys-48 (or Cys-69 of wild-type apoflavodoxin) perturbs the unfolded protein, because hydrophobic interactions occur between the label and hydrophobic patches of unfolded apoflavodoxin. Concomitant hydrophobic free energy changes of the unfolded protein (and possibly of the off-pathway intermediate) reduce the stability of native spin-labelled protein against unfolding. In addition, attachment of MTSL or CMTSL to Cys-48 induces the presence of distinct states in unfolded apoflavodoxin. Despite these difficulties, the spin-label data obtained here show that non-native contacts exist between transiently ordered structured elements in unfolded apoflavodoxin

The Carbohydrate-Binding Site in Galectin-3 Is Preorganized To Recognize a Sugarlike Framework of Oxygens: Ultra-High-Resolution Structures and Water Dynamics

The recognition of carbohydrates by proteins is a fundamental aspect of communication within and between living cells. Understanding the molecular basis of carbohydrate-protein interactions is a prerequisite for the rational design of synthetic ligands. Here we report the high- to ultrahigh-resolution crystal structures of the carbohydrate recognition domain of galectin-3 (Gal3C) in the ligand-free state (1.08 angstrom at 100 K, 1.25 angstrom at 298 K) and in complex with lactose (0.86 angstrom) or glycerol (0.9 angstrom). These structures reveal striking similarities in the positions of water and carbohydrate oxygen atoms in all three states, indicating that the binding site of Gal3C is preorganized to coordinate oxygen atoms in an arrangement that is nearly optimal for the recognition of beta-galactosides. Deuterium nuclear magnetic resonance (NMR) relaxation dispersion experiments and molecular dynamics simulations demonstrate that all water molecules in the lactose-binding site exchange with bulk water on a time scale of nanoseconds or shorter. Nevertheless, molecular dynamics simulations identify transient water binding at sites that agree well with those observed by crystallography, indicating that the energy landscape of the binding site is maintained in solution. All heavy atoms of glycerol are positioned like the corresponding atoms of lactose in the Gal3C complexes. However, binding of glycerol to Gal3C is insignificant in solution at room temperature, as monitored by NMR spectroscopy or isothermal titration calorimetry under conditions where lactose binding is readily detected. These observations make a case for protein cryo-crystallography as a valuable screening method in fragment-based drug discovery and further suggest that identification of water sites might inform inhibitor design