On the Calculation of SAXS Profiles of Folded and Intrinsically Disordered Proteins from Computer Simulations

Solution techniques such as small-angle X-ray scattering (SAXS) play a central role in structural studies of intrinsically disordered proteins (IDPs); yet, due to low resolution, it is generally necessary to combine SAXS with additional experimental sources of data and to use molecular simulations. Computational methods for the calculation of theoretical SAXS intensity profiles can be separated into two groups, depending on whether the solvent is modeled implicitly as continuous electron density or considered explicitly. The former offers reduced computational cost but requires the definition of a number of free parameters to account for, for example, the excess density of the solvation layer. Overfitting can thus be an issue, particularly when the structural ensemble is unknown. Here, we investigate and show how small variations of the contrast of the hydration shell, δρ, severely affect the outcome, analysis and interpretation of computed SAXS profiles for folded and disordered proteins. For both the folded and disordered proteins studied here, using a default δρ may, in some cases, result in the calculation of non-representative SAXS profiles, leading to an overestimation of their size and a misinterpretation of their structural nature. The solvation layer of the different IDP simulations also impacts their size estimates differently, depending on the protein force field used. The same is not true for the folded protein simulations, suggesting differences in the solvation of the two classes of proteins, and indicating that different force fields optimized for IDPs may cause expansion of the polypeptide chain through different physical mechanisms.     

Deletion of Phenylalanine 508 in the First Nucleotide-binding Domain of the Cystic Fibrosis Transmembrane Conductance Regulator Increases Conformational Exchange and Inhibits Dimerization

Deletion of Phe-508 (F508del) in the first nucleotide-binding domain (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR) results in destabilization of the domain, intramolecular interactions involving the domain, and the entire channel. The destabilization caused by F508del manifests itself in defective channel processing and channel gating defects. Here, we present NMR studies of the effect of F508del and the I539T stabilizing mutation on NBD1 dynamics, with a view to understanding these changes in stability. Qualitatively, F508del NMR spectra exhibit significantly more peak broadening than WT spectra due to the enhanced intermediate time scale (millisecond to microsecond) motions in the mutant. Unexpectedly, studies of fast (nanosecond to picosecond) motions revealed that F508del NBD1 tumbles more rapidly in solution than WT NBD1. Whereas F508del tumbles at a rate nearly consistent with the monomeric state, the WT protein tumbles significantly more slowly. Paramagnetic relaxation enhancement experiments confirm that NBD1 homodimerizes in solution in the expected head-to-tail orientation. NMR spectra of WT NBD1 reveal significant concentration-dependent chemical shift perturbations consistent with NBD1 dimerization. Chemical shift analysis suggests that the more rapid tumbling of F508del is the result of an impaired ability to dimerize. Based on previously published crystal structures and NMR spectra of various NBD1 mutants, we propose that deletion of Phe-508 affects Q-loop conformational sampling in a manner that inhibits dimerization. These results provide a potential mechanism for inhibition of channel opening by F508del and support the dimer interface as a target for cystic fibrosis therapeutics.     

Development of potent anti-infective agents from Silurana tropicalis: Conformational analysis of the amphipathic,alpha-helical antimicrobial peptide XT-7 and its non-haemolytic analogue [G4K]XT-7

Peptide XT-7 (GLLGP⁵LLKIA¹⁰AKVGS¹⁵NLL.NH₂) is a cationic, leucine-rich peptide, first isolated from skin secretions of the frog, Silurana tropicalis (Pipidae). The peptide shows potent, broad-spectrum antimicrobial activity but its therapeutic potential is limited by haemolytic activity (LC₅₀ = 140 µM). The analogue [G4K]XT-7, however, retains potent antimicrobial activity but is non-haemolytic (LC₅₀ > 500 µM). In order to elucidate the molecular basis for this difference in properties, the three dimensional structures of XT-7 and the analogue have been investigated by proton NMR spectroscopy and molecular modelling. In aqueous solution, both peptides lack secondary structure. In a 2,2,2-trifluoroethanol (TFE-d₃)-H₂O mixed solvent system, XT-7 is characterised by a right handed α-helical conformation between residues Leu³ and Leu¹⁷ whereas [G4K]XT-7 adopts a more restricted α-helical conformation between residues Leu⁶ and Leu¹⁷. A similar conformation for XT-7 in 1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC) micellular media was observed with a helical segment between Leu³ and Leu¹⁷. However, differences in side chain orientations restricting the hydrophilic residues to a smaller patch resulted in an increased hydrophobic surface relative to the conformation in TFE-H₂O. Molecular modelling of the structures obtained in our study demonstrates the amphipathic character of the helical segments. It is proposed that the marked decrease in haemolytic activity produced by the substitution Gly⁴ → Lys in XT-7 arises from a decrease in both helicity and hydrophobicity. These studies may facilitate the development of potent but non-toxic anti-infective agents based upon the structure of XT-7.     

Role of intrinsic flexibility in signal transduction mediated by the cell cycle regulator, p27 Kip1

p27^Kip1 (p27), which controls eukaryotic cell division through interactions with cyclin-dependent kinases (Cdks), integrates and transduces promitogenic signals from various nonreceptor tyrosine kinases by orchestrating its own phosphorylation, ubiquitination and degradation. Intrinsic flexibility allows p27 to act as a “conduit” for sequential signaling mediated by tyrosine and threonine phosphorylation and ubiquitination. While the structural features of the Cdk/cyclin-binding domain of p27 are understood, how the C-terminal regulatory domain coordinates multistep signaling leading to p27 degradation is poorly understood. We show that the 100-residue p27 C-terminal domain is extended and flexible when p27 is bound to Cdk2/cyclin A. We propose that the intrinsic flexibility of p27 provides a molecular basis for the sequential signal transduction conduit that regulates p27 degradation and cell division. Other intrinsically unstructured proteins possessing multiple sites of posttranslational modification may participate in similar signaling conduits.     

Structure and Disorder in an Unfolded State under Nondenaturing Conditions from Ensemble Models Consistent with a Large Number of Experimental Restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include ¹H-¹⁵N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O₂-induced ¹³C paramagnetic shifts, hydrogen-exchange protection factors, and ¹⁵N R₂ data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu²⁺-Ni²⁺ binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative α-helical structure and both native-like and nonnative tertiary structure.     

Backbone conformational flexibility of the lipid modified membrane anchor of the human N-Ras protein investigated by solid-state NMR and molecular dynamics simulation

The lipid modified human N-Ras protein, implicated in human cancer development, is of particular interest due to its membrane anchor that determines the activity and subcellular location of the protein. Previous solid-state NMR investigations indicated that this membrane anchor is highly dynamic, which may be indicative of backbone conformational flexibility. This article aims to address if a dynamic exchange between three structural models exist that had been determined previously. We applied a combination of solid-state nuclear magnetic resonance (NMR) methods and replica exchange molecular dynamics (MD) simulations using a Ras peptide that represents the terminal seven amino acids of the human N-Ras protein. Analysis of correlations between the conformations of individual amino acids revealed that Cys 181 and Met 182 undergo collective conformational exchange. Two major structures constituting about 60% of all conformations could be identified. The two conformations found in the simulation are in rapid exchange, which gives rise to low backbone order parameters and nuclear spin relaxation as measured by experimental NMR methods. These parameters were also determined from two 300 ns conventional MD simulations, providing very good agreement with the experimental data.     

Structure determination of UL49.5 transmembrane protein from bovine herpesvirus 1 by NMR spectroscopy and molecular dynamics

The transporter associated with antigen processing (TAP) directly participates in the immune response as a key component of the cytosolic peptide to major histocompatibility complex (MHC) class I protein loading machinery. This makes TAP an important target for viruses avoiding recognition by CD8+ T lymphocytes. Its activity can be suppressed by the UL49.5 protein produced by bovine herpesvirus 1, although the mechanism of this inhibition has not been understood so far.Therefore, the main goal of our study was to investigate the 3D structure of bovine herpesvirus 1 - encoded UL49.5 protein. The final structure of the inhibitor was established using circular dichroism (CD), 2D nuclear magnetic resonance (NMR), and molecular dynamics (MD) in membrane mimetic environments. In NMR studies, UL49.5 was represented by two fragments: the extracellular region (residues 1–35) and the transmembrane-intracellular fragment (residues 36–75), displaying various functions during viral invasion. After the empirical structure determination, a molecular docking procedure was used to predict the complex of UL49.5 with the TAP heterodimer.Our results revealed that UL49.5 adopted a highly flexible membrane-proximal helical structure in the extracellular part. In the transmembrane region, we observed two short α-helices. Furthermore, the cytoplasmic part had an unordered structure. Finally, we propose three different orientations of UL49.5 in the complex with TAP. Our studies provide, for the first time, the experimental structural information on UL49.5 and structure-based insight in its mechanism of action which might be helpful in designing new drugs against viral infections.     

Amphipathic helical peptides hamper protein-protein interactions of the intrinsically disordered chromatin nuclear protein 1 (NUPR1)

Background

NUPR1 is a multifunctional intrinsically disordered protein (IDP) involved, among other functions, in chromatin remodelling, and development of pancreatic ductal adenocarcinoma (PDAC). It interacts with several biomolecules through hydrophobic patches around residues Ala33 and Thr68. The drug trifluoperazine (TFP), which hampers PDAC development in xenografted mice, also binds to those regions. Because of the large size of the hot-spot interface of NUPR1, small molecules could not be adequate to modulate its functions.

Structural basis of the signal transduction via transmembrane domain of the human growth hormone receptor

Background

Prior studies of the human growth hormone receptor (GHR) revealed a distinct role of spatial rearrangements of its dimeric transmembrane domain in signal transduction across membrane. Detailed structural information obtained in the present study allowed elucidating the bases of such rearrangement and provided novel insights into receptor functioning.

Improved structural characterizations of the drkN SH3 domain unfolded state suggest a compact ensemble with native-like and non-native structure

Due to their dynamic ensemble nature and a deficiency of experimental restraints, disordered states of proteins are difficult to characterize structurally. Here, we have expanded upon our previous work on the unfolded state of the Drosophila drk N-terminal (drkN) SH3 domain with our program ENSEMBLE, which assigns population weights to pregenerated conformers in order to calculate ensembles of structures whose properties are collectively consistent with experimental measurements. The experimental restraint set has been enlarged with newly measured paramagnetic relaxation enhancements from Cu(2+) bound to an amino terminal Cu(2+)-Ni(2+) binding (ATCUN) motif as well as nuclear Overhauser effect (NOE) and hydrogen exchange data from recent studies. In addition, two new pseudo-energy minimization algorithms have been implemented that have dramatically improved the speed of ENSEMBLE population weight assignment. Finally, we have greatly improved our conformational sampling by utilizing a variety of techniques to generate both random structures and structures that are biased to contain elements of native-like or non-native structure. Although it is not possible to uniquely define a representative structural ensemble, we have been able to assess various properties of the drkN SH3 domain unfolded state by performing ENSEMBLE minimizations of different conformer pools. Specifically, we have found that the experimental restraint set enforces a compact structural distribution that is not consistent with an overall native-like topology but shows preference for local non-native structure in the regions corresponding to the diverging turn and the beta5 strand of the folded state and for local native-like structure in the region corresponding to the beta6 and beta7 strands. We suggest that this approach could be generally useful for the structural characterization of disordered states.     

Characterization of conformational and dynamic properties of natively unfolded human and mouse alpha-synuclein ensembles by NMR: implication for aggregation

Conversion of human α-synuclein (aS) from the free soluble state to the insoluble fibrillar state has been implicated in the etiology of Parkinson's disease. Human aS is highly homologous in amino acid sequence to mouse aS, which contains seven substitutions including the A53T that has been linked to familial Parkinson's disease, and including five substitutions in the C-terminal region. It has been shown that the rate of fibrillation is highly dependent on the exact sequence of the protein, and mouse aS is reported to aggregate more rapidly than human aS in vitro. Nuclear magnetic resonance experiments of mouse and human aS at supercooled temperatures (263 K) are used to understand the effect of sequence on conformational fluctuations in the disordered ensembles and to relate these to differences in propensities to aggregate. We show that both aS are natively unfolded at low temperature with different propensities to secondary structure, backbone dynamics and long-range contacts across the protein. Mouse aS exhibits a higher propensity to helical conformation around the C-terminal substitutions as well as the loss of transient long-range contacts from the C- to the N-terminal end and hydrophobic central regions of the protein relative to human aS. Lack of back-folding from the C-terminal end of mouse aS exposes the N-terminal region, which is shown, by ¹⁵N relaxation experiments, to be very restricted in mobility relative to human aS. We propose that the restricted mobility in the N-terminal region may arise from transient interchain interactions, suggesting that the N-terminal KTK(E/Q)GV repeats may serve as initiation sites for aggregation in mouse aS. These transient interchain interactions coupled with a non-Aβ amyloid component (NAC) region that is both more exposed and has a higher propensity to β structure may accelerate the rate of fibril formation of aS.     

Simulating structural and thermodynamic properties of carcinogen-damaged DNA

A pair of stereoisomeric covalent adducts to guanine in double-stranded DNA, derived from the reaction of mutagenic and tumorigenic metabolites of benzo[a]pyrene, have been well characterized structurally and thermodynamically. Both high-resolution NMR solution structures and an array of thermodynamic data are available for these 10S (+)- and 10R (-)-trans-anti -[BP]-N(2)-dG adducts in double-stranded deoxyoligonucleotides. The availability of experimentally well-characterized duplexes containing these two stereoisomeric guanine adducts provides an opportunity for evaluating the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method for computing thermodynamic properties from molecular dynamics ensembles. We have carried out 3-ns molecular dynamics simulations, using NMR solution structures as the starting models for the 10S (+)- and 10R (-)-trans-anti-dG adducts in a DNA duplex 11-mer using AMBER 6.0. We employed the MM-PBSA method to compute the free energies, enthalpies, and entropies of the two adducts. Our complete thermodynamic analysis agrees quite well with the full experimental thermodynamic characterization of these adducts, showing essentially equal stabilities of the two adducts. We also calculated the nuclear Overhauser effect (NOE) distances from the molecular dynamics trajectories, and compared them against the experimental NMR-derived NOE distances. Our results showed that the simulated structures are in good agreement with the NMR experimental NOE data. Furthermore, the molecular dynamics simulations provided new structural and biological insights. Specifically, the puzzling observation that the BP aromatic ring system in the 10S (+)-trans-anti-dG adduct is more exposed to the aqueous solvent than the 10R (-)-trans-anti-dG adduct, is rationalized in terms of the adduct structures. The structural and thermodynamic features of these stereoisomeric adducts are also discussed in relation to their reported low susceptibilities to nucleotide excision repair.     

Three-dimensional solution structure and conformational plasticity of the N-terminal scavenger receptor cysteine-rich domain of human CD5

The lymphocyte receptor CD5 influences cell activation by modifying the strength of the intracellular response initiated by antigen engagement. Regulation through CD5 involves the interaction of one or more of its three scavenger receptor cysteine-rich domains present in the extracellular region. Here, we present the 3D solution structure of a non-glycosylated double mutant of the N-terminal domain of human CD5 expressed in Escherichia coli (eCD5d1m), which has enhanced solubility compared to the non-glycosylated wild-type (eCD5d1). In common with a glycosylated form expressed in Pichia pastoris, the [¹⁵N,¹H]-correlation spectra of both eCD5d1 and eCD5d1m exhibit non-uniform temperature-dependent signal intensities, indicating extensive conformational fluctuations on the micro-millisecond timescale. Although approximately one half of the signals expected for the domain are absent at 298 K, essentially complete resonance assignments and a solution structure could be obtained at 318 K. Because of the sparse nature of the experimental restraint data and the potentially important contribution of conformational exchange to the nuclear Overhauser effect peak intensity, we applied inferential structure determination to calculate the eCD5d1m structure. The inferential structure determination ensemble has similar features to that obtained by traditional simulated annealing methods, but displays superior definition and structural quality. The eCD5d1m structure is similar to other members of the scavenger receptor cysteine-rich superfamily, but the position of the lone α helix differs due to interactions with the unique N-terminal region of the domain. The availability of an experimentally tractable form of CD5d1, together with its 3D structure, provides new tools for further investigation of its function within intact CD5.     

Tryptophan- and arginine-rich antimicrobial peptides: Structures and mechanisms of action

Antimicrobial peptides encompass a number of different classes, including those that are rich in a particular amino acid. An important subset are peptides rich in Arg and Trp residues, such as indolicidin and tritrpticin, that have broad and potent antimicrobial activity. The importance of these two amino acids for antimicrobial activity was highlighted through the screening of a complete combinatorial library of hexapeptides. These residues possess some crucial chemical properties that make them suitable components of antimicrobial peptides. Trp has a distinct preference for the interfacial region of lipid bilayers, while Arg residues endow the peptides with cationic charges and hydrogen bonding properties necessary for interaction with the abundant anionic components of bacterial membranes. In combination, these two residues are capable of participating in cation-π interactions, thereby facilitating enhanced peptide-membrane interactions. Trp sidechains are also implicated in peptide and protein folding in aqueous solution, where they contribute by maintaining native and nonnative hydrophobic contacts. This has been observed for the antimicrobial peptide from human lactoferrin, possibly restraining the peptide structure in a suitable conformation to interact with the bacterial membrane. These unique properties make the Arg- and Trp-rich antimicrobial peptides highly active even at very short peptide lengths. Moreover, they lead to structures for membrane-mimetic bound peptides that go far beyond regular α-helices and β-sheet structures. In this review, the structures of a number of different Trp- and Arg-rich antimicrobial peptides are examined and some of the major mechanistic studies are presented.