Secondary Structure Propensities Research Articles

Assessment of transferable forcefields for protein simulations attests improved description of disordered states and secondary structure propensities, and hints at multi-protein systems as the next challenge for optimization

Continuous assessment of transferable forcefields for molecular simulations is essential to identify their weaknesses and direct improvement efforts. The latest efforts focused on better describing disordered proteins while retaining proper description of folded domains, important because forcefields of the previous generations produce overly compact disordered states. Such improvements should additionally alleviate the related problem of over-stabilized protein–protein interactions, which has been largely overlooked. Here we evaluated three state-of-the-art forcefields, current flagships of their respective developers, optimized for ordered and disordered proteins: CHARMM36m with its recommended corrected TIP3P* water, ff19SB with the recommended OPC water, and the 2019 a99SBdisp forcefield by D. E. Shaw Research with its modified TIP4P water; plus ff14SB with TIP3P as an example of the former generation of forcefields. Our evaluation entailed simulations of (i) multiple copies of a protein that is highly soluble yet undergoes weak dimerization, (ii) a disordered peptide with low, well-characterized alpha helical propensity, and (iii) a peptide known to form insoluble β-aggregates. Our results recapitulate ff14SB-TIP3P over-stabilizing aggregates and secondary structures and place a99SBdisp-TIP4PD at the other end i.e. predicting overly weak intermolecular interactions despite reasonably predicting secondary structure propensities. In-between, CHARMM36m-TIP3P* still over-stabilizes aggregates but predicts residue-wise alpha helical propensities in solution slightly better than ff19SB-OPC, while ff19SB-OPC poses the best prediction of weak dimerization of the soluble protein still predicting aggregation of the β-peptides. This independent assessment shows that the claimed forcefield improvements are real, but also that a right balance between noncovalent attraction and repulsion has not yet been reached. We thus propose developers to consider systems like those tested here in their forcefield tuning protocols. Last, the good performance of CHARMM36m-TIP3P* further shows that tuning 3-point water models might still be an alternative to the more costly 4-point models like OPC and TIP4PD.

The MemMoRF database for recognizing disordered protein regions interacting with cellular membranes.

Protein and lipid membrane interactions play fundamental roles in a large number of cellular processes (e.g. signalling, vesicle trafficking, or viral invasion). A growing number of examples indicate that such interactions can also rely on intrinsically disordered protein regions (IDRs), which can form specific reversible interactions not only with proteins but also with lipids. We named IDRs involved in such membrane lipid-induced disorder-to-order transition as MemMoRFs, in an analogy to IDRs exhibiting disorder-to-order transition upon interaction with protein partners termed Molecular Recognition Features (MoRFs). Currently, both the experimental detection and computational characterization of MemMoRFs are challenging, and information about these regions are scattered in the literature. To facilitate the related investigations we generated a comprehensive database of experimentally validated MemMoRFs based on manual curation of literature and structural data. To characterize the dynamics of MemMoRFs, secondary structure propensity and flexibility calculated from nuclear magnetic resonance chemical shifts were incorporated into the database. These data were supplemented by inclusion of sentences from papers, functional data and disease-related information. The MemMoRF database can be accessed via a user-friendly interface at https://memmorf.hegelab.org, potentially providing a central resource for the characterization of disordered regions in transmembrane and membrane-associated proteins.

Solid state NMR assignments of a human \u03bb-III immunoglobulin light chain amyloid fibril

The aggregation of antibody light chains is linked to systemic light chain (AL) amyloidosis, a disease where amyloid deposits frequently affect the heart and the kidney. We here investigate fibrils from the λ-III FOR005 light chain (LC), which is derived from an AL-patient with severe cardiac involvement. In FOR005, five residues are mutated with respect to its closest germline gene segment IGLV3-19 and IGLJ3. All mutations are located close to the complementarity determining regions (CDRs). The sequence segments responsible for the fibril formation are not yet known. We use fibrils extracted from the heart of this particular amyloidosis patient as seeds to prepare fibrils for solid-state NMR. We show that the seeds induce the formation of a specific fibril structure from the biochemically produced protein. We have assigned the fibril core region of the FOR005-derived fibrils and characterized the secondary structure propensity of the observed amino acids. As the primary structure of the aggregated patient protein is different for every AL patient, it is important to study, analyze and report a greater number of light chain sequences associated with AL amyloidosis.

A survey of RNA secondary structural propensity encoded within human herpesvirus genomes: global comparisons and local motifs.

There are nine herpesviruses known to infect humans, of which Epstein–Barr virus (EBV) is the most widely distributed (>90% of adults infected). This ubiquitous virus is implicated in a variety of cancers and autoimmune diseases. Previous analyses of the EBV genome revealed numerous regions with evidence of generating unusually stable and conserved RNA secondary structures and led to the discovery of a novel class of EBV non-coding (nc)RNAs: the stable intronic sequence (sis)RNAs. To gain a better understanding of the roles of RNA structure in EBV biology and pathogenicity, we revisit EBV using recently developed tools for genome-wide motif discovery and RNA structural characterization. This corroborated previous results and revealed novel motifs with potential functionality; one of which has been experimentally validated. Additionally, since many herpesviruses increasingly rival the seroprevalence of EBV (VZV, HHV-6 and HHV-7 being the most notable), analyses were expanded to include all sequenced human Herpesvirus RefSeq genomes, allowing for genomic comparisons. In total 10 genomes were analyzed, for EBV (types 1 and 2), HCMV, HHV-6A, HHV-6B, HHV-7, HSV-1, HSV-2, KSHV, and VZV. All resulting data were archived in the RNAStructuromeDB (https://structurome.bb.iastate.edu/herpesvirus) to make them available to a wide array of researchers.

Predicting Secondary Structure Propensities in IDPs Using Simple Statistics from Three-Residue Fragments

Intrinsically disordered proteins (IDPs) play key functional roles facilitated by their inherent plasticity. In most of the cases, IDPs recognize their partners through partially structured elements inserted in fully disordered chains. The identification and characterization of these elements is fundamental to understand the functional mechanisms of IDPs. Although several computational methods have been developed to identify order within disordered chains, most of the current secondary structure predictors are focused on globular proteins and are not necessarily appropriate for IDPs. Here, we present a comprehensible method, called Local Structural Propensity Predictor (LS2P), to predict secondary structure elements from IDP sequences. LS2P performs statistical analyses from a database of three-residue fragments extracted from coil regions of high-resolution protein structures. In addition to identifying scarcely populated helical and extended regions, the method pinpoints short stretches triggering β-turn formation or promoting α-helices. The simplicity of the method enables a direct connection between experimental observations and structural features encoded in IDP sequences.

New facets of larger Nest motifs in proteins.

The Nest is a concave-shaped structural motif in proteins formed by consecutive enantiomeric left-handed (L) and right-handed (R) helical conformation of the backbone. This important motif subsumes many turn and helix capping structures and binds electron-rich ligands. Simple Nests are either RL or LR. Larger Nests (>2 residues long) may be RLR, LRL, RLRL, and so forth, being considered as composed of overlapping simple Nests. The larger Nests remain under-explored despite their widely known contributions to protein function. In our study, we address whether the recurrence of enantiomeric geometry in the larger Nests constrains the peptide backbone such that distinct compositional and conformational preferences are seen compared to simple Nests. Our analysis reveals the critical role of the L helical torsion angle in the formation of larger Nests. This can be observed through the higher propensity of residue or secondary structure combinations in LR and LRL backbone conformation in comparison to RL or RLR, although LR/LRL is considerably lower by occurrence. We also find that the most abundant doublets and triplets in Nests have a propensity for particular secondary structures, suggesting a strong sequence-structure relationship in the larger Nest. Overall, our analysis corroborates distinct features of simple and the larger Nests. Such insights would be helpful towards in-vitro design of peptides and peptidomimetic studies.

Role of Tyr-39 for the Structural Features of α-Synuclein and for the Interaction with a Strong Modulator of Its Amyloid Assembly.

Recent studies suggest that Tyr-39 might play a critical role for both the normal function and the pathological dysfunction of α-synuclein (αS), an intrinsically disordered protein involved in Parkinson’s disease. We perform here a comparative analysis between the structural features of human αS and its Y39A, Y39F, and Y39L variants. By the combined application of site-directed mutagenesis, biophysical techniques, and enhanced sampling molecular simulations, we show that removing aromatic functionality at position 39 of monomeric αS leads to protein variants populating more compact conformations, conserving its disordered nature and secondary structure propensities. Contrasting with the subtle changes induced by mutations on the protein structure, removing aromaticity at position 39 impacts strongly on the interaction of αS with the potent amyloid inhibitor phthalocyanine tetrasulfonate (PcTS). Our findings further support the role of Tyr-39 in forming essential inter and intramolecular contacts that might have important repercussions for the function and the dysfunction of αS.

RNA genome conservation and secondary structure in SARS-CoV-2 and SARS-related viruses: a first look.

As the COVID-19 outbreak spreads, there is a growing need for a compilation of conserved RNA genome regions in the SARS-CoV-2 virus along with their structural propensities to guide development of antivirals and diagnostics. Here we present a first look at RNA sequence conservation and structural propensities in the SARS-CoV-2 genome. Using sequence alignments spanning a range of betacoronaviruses, we rank genomic regions by RNA sequence conservation, identifying 79 regions of length at least 15 nt as exactly conserved over SARS-related complete genome sequences available near the beginning of the COVID-19 outbreak. We then confirm the conservation of the majority of these genome regions across 739 SARS-CoV-2 sequences subsequently reported from the COVID-19 outbreak, and we present a curated list of 30 “SARS-related-conserved” regions. We find that known RNA structured elements curated as Rfam families and in prior literature are enriched in these conserved genome regions, and we predict additional conserved, stable secondary structures across the viral genome. We provide 106 “SARS-CoV-2-conserved-structured” regions as potential targets for antivirals that bind to structured RNA. We further provide detailed secondary structure models for the extended 5′ UTR, frameshifting stimulation element, and 3′ UTR. Lastly, we predict regions of the SARS-CoV-2 viral genome that have low propensity for RNA secondary structure and are conserved within SARS-CoV-2 strains. These 59 “SARS-CoV-2-conserved-unstructured” genomic regions may be most easily accessible by hybridization in primer-based diagnostic strategies.

PH-dependent secondary structure propensity of the influenza A virus M2 cytoplasmic tail

The cytoplasmic C-terminal tail of the matrix protein 2 (M2) from influenza A virus has a well conserved sequence and is involved in interactions with several host proteins as well as the influenza matrix protein 1 (M1). Whereas the transmembrane domain of M2 has been well characterised structurally and functionally, high resolution information about the distal cytoplasmic tail is lacking. Here we report the chemical shifts of the cytoplasmic tail of M2 and the chemical shift perturbations at low pH and in the presence of membrane mimetics. The cytoplasmic tail residues are mostly disordered but an extended backbone conformation is adopted by the LC3 binding motif and the putative M1 interaction site has partial helical content with a small pH-dependence. The chemical shift assignments provide a basis for further investigations into interactions of the M2 cytoplasmic tail with viral and host cell factors.

The Dynamic Search Mode of a Disordered Transcription Factor

The bZIPs are the largest family of transcription factors to operate exclusively in eukaryotes. DNA binding domains are typically thought of as highly structured, yet the bZIP domain is intrinsically disordered in isolation. We have performed a combined structural and kinetic analysis of the process by which CREB (cAMP response-element binding protein), a model homodimeric bZIP, finds its DNA target site amongst vast excesses of non-specific DNA. Using NMR we have characterised for both monomer and dimer the secondary structure propensities and dynamics with both target and non-specific DNA.

Atypical structural tendencies among low-complexity domains in the Protein Data Bank proteome.

A variety of studies have suggested that low-complexity domains (LCDs) tend to be intrinsically disordered and are relatively rare within structured proteins in the Protein Data Bank (PDB). Although LCDs are often treated as a single class, we previously found that LCDs enriched in different amino acids can exhibit substantial differences in protein metabolism and function. Therefore, we wondered whether the structural conformations of LCDs are likewise dependent on which specific amino acids are enriched within each LCD. Here, we directly examined relationships between enrichment of individual amino acids and secondary structure tendencies across the entire PDB proteome. Secondary structure tendencies varied as a function of the identity of the amino acid enriched and its degree of enrichment. Furthermore, divergence in secondary structure profiles often occurred for LCDs enriched in physicochemically similar amino acids (e.g. valine vs. leucine), indicating that LCDs composed of related amino acids can have distinct secondary structure tendencies. Comparison of LCD secondary structure tendencies with numerous pre-existing secondary structure propensity scales resulted in relatively poor correlations for certain types of LCDs, indicating that these scales may not capture secondary structure tendencies as sequence complexity decreases. Collectively, these observations provide a highly resolved view of structural tendencies among LCDs parsed by the nature and magnitude of single amino acid enrichment.

Prediction of polyproline II secondary structure propensity in proteins.

Background: The polyproline II helix (PPIIH) is an extended protein left-handed secondary structure that usually but not necessarily involves prolines. Short PPIIHs are frequently, but not exclusively, found in disordered protein regions, where they may interact with peptide-binding domains. However, no readily usable software is available to predict this state. Results: We developed PPIIPRED to predict polyproline II helix secondary structure from protein sequences, using bidirectional recurrent neural networks trained on known three-dimensional structures with dihedral angle filtering. The performance of the method was evaluated in an external validation set. In addition to proline, PPIIPRED favours amino acids whose side chains extend from the backbone (Leu, Met, Lys, Arg, Glu, Gln), as well as Ala and Val. Utility for individual residue predictions is restricted by the rarity of the PPIIH feature compared to structurally common features. Conclusion: The software, available at http://bioware.ucd.ie/PPIIPRED, is useful in large-scale studies, such as evolutionary analyses of PPIIH, or computationally reducing large datasets of candidate binding peptides for further experimental validation.

Hamiltonian Reweighing To Refine Protein Backbone Dihedral Angle Parameters in the GROMOS Force Field.

Molecular dynamics simulations of proteins depend critically on the underlying force field, which may be parameterized against experimental data or high-quality quantum calculations. Here, we develop search algorithms based on Monte Carlo and steepest descent calculations to optimize the backbone dihedral angle parameters from a single reference simulation. We apply these tools to improve the agreement between simulations of single, capped amino acids and experimentally determined J values and secondary structure propensities of these molecules. The parameters are further refined based on simulations of a set of seven proteins and finally validated in simulations on a large set of 52 protein structures. Improvements in the dihedral angle distributions are observed, and structural propensities of the proteins are reproduced very well. Overall, the GROMOS 54A8_bb parameter set forms an improvement to previous parameter sets, both for small molecules and for protein simulations.

Structure of an Unfolding Intermediate of an RRM Domain of ETR-3 Reveals Its Native-like Fold

Prevalence of one or more partially folded intermediates during protein unfolding with different secondary and ternary conformations has been identified as an integral character of protein unfolding. These transition-state species need to be characterized structurally for elucidation of their folding pathways. We have determined the three-dimensional structure of an intermediate state with increased conformational space sampling under urea-denaturing condition. The protein unfolds completely at 10 M urea but retains residual secondary structural propensities with restricted motion. Here, we describe the native state, observable intermediate state, and unfolded state for ETR-3 RRM-3, which has canonical RRM fold. These observations can shed more light on unfolding events for RRM-containing proteins.

Quasi-native transition and self-diffusion of proteins in water-glycerol mixture

An orderly investigation of the levels of secondary and tertiary structures, kinetics of tertiary structural changes, and self diffusion coefficient of lysozyme and cytochrome c in the 0―70% (weight/volume) range of glycerol is reported. While secondary structural propensity of both proteins is larger in glycerol, results for tertiary structure and translational diffusion coefficient with increasing glycerol provide two contrasting depictions ― lysozyme becomes increasingly compact, plausibly due to disulfide bridge constraints, but cytochrome c expands and loses the tertiary structure. The chain expansion and contraction corresponding to loss and reformation of tertiary structure of cytochrome c are ultrafast that occur in the submillisecond bin. Changes in protein conformation appear in as little as 2% glycerol, and the results suggest that glycerol does not unfold the protein but reversibly destabilizes to quasi-native state(s). These observations make one ponder whether results of studies on protein dynamics, relaxation, and conformational substates reported in literature can be associated with native-state properties.

Do sequence neighbours of intrinsically disordered regions promote structural flexibility in intrinsically disordered proteins?

Intrinsically disordered proteins (IDPs) are crucial players in various cellular activities. Several experimental and computational analyses have been conducted to study structural pliability and functional potential of IDPs. In spite of active research in past few decades, what induces structural disorder in IDPs and how is still elusive. Many studies testify that sequential and spatial neighbours often play important roles in determining structural and functional behaviour of proteins. Considering this fact, we assessed sequence neighbours of intrinsically disordered regions (IDRs) to understand if they have any role to play in inducing structural flexibility in IDPs. Our analysis includes 97% eukaryotic IDPs and 3% from bacteria and viruses. Physicochemical and structural parameters including amino acid propensity, hydrophobicity, secondary structure propensity, relative solvent accessibility, B-factor and atomic packing density are used to characterise the neighbouring residues of IDRs (NRIs). We show that NRIs exhibit a unique nature, which makes them stand out from both ordered and disordered residues. They show correlative occurrences of residue pairs like Ser-Thr and Gln-Asn, indicating their tendency to avoid strong biases of order or disorder promoting amino acids. We also find differential preferences of amino acids between N- and C-terminal neighbours, which might indicate a plausible directional effect on the dynamics of adjacent IDRs. We designed an efficient prediction tool using Random Forest to distinguish the NRIs from the ordered residues. Our findings will contribute to understand the behaviour of IDPs, and may provide potential lead in deciphering the role of IDRs in protein folding and assembly.

Molecular mechanics and quantum chemical calculations unveil the combating effect of baicalein on human islet amyloid polypeptide aggregates

ABSTRACTFormation of insoluble toxic aggregates by hIAPP polypeptide is found to be a core component for disease pathogenicity of patients suffering from type II diabetes. Naturally occurring polyphenols that possess anti-aggregation property are being majorly studied in the current scenario for treating various conformational diseases. Herein, we study the disaggregation mechanism of naturally occurring polyphenol baicalein, obtained from the roots of S. baicalensis and Indian trumpet flower on hIAPP dimer using quantum chemical calculation and discrete molecular dynamics. Our study reported that a drastic loss in the secondary structural propensity of hIAPP was seen upon binding of baicalein. Notably, the hydrophobic core and the phenolic groups present in the ends of baicalein molecule play a key role in inhibiting the aggregates formed upon binding to the amyloidogenic core region of hIAPP. Thus, our study provides a comprehensive understanding over the disaggregation effect of baicalein on hIAPP dimer from a computational point of view and thereby bridging the gap for future therapeutic strategy in designing the anti-aggregation compounds that inhibit hIAPP amyloids.

Structural and Aggregation Features of a Human κ-Casein Fragment with Antitumor and Cell-Penetrating Properties.

Intrinsically disordered proteins play a central role in dynamic regulatory and assembly processes in the cell. Recently, a human κ-casein proteolytic fragment called lactaptin (8.6 kDa) was found to induce apoptosis of human breast adenocarcinoma MCF-7 and MDA-MB-231 cells with no cytotoxic activity toward normal cells. Earlier, we had designed some recombinant analogs of lactaptin and compared their biological activity. Among these analogs, RL2 has the highest antitumor activity, but the amino acid residues and secondary structures that are responsible for RL2′s activity remain unclear. To elucidate the structure–activity relations of RL2, we studied the structural and aggregation features of this fairly large intrinsically disordered fragment of human milk κ-casein by a combination of physicochemical methods: NMR, paramagnetic relaxation enhancement (PRE), Electron Paramagnetic Resonance (EPR), circular dichroism, dynamic light scattering, atomic force microscopy, and a cytotoxic activity assay. It was found that in solution, RL2 exists as stand-alone monomeric particles and large aggregates. Whereas the disulfide-bonded homodimer turned out to be more prone to assembly into large aggregates, the monomer predominantly forms single particles. NMR relaxation analysis of spin-labeled RL2 showed that the RL2 N-terminal region, which is essential not only for multimerization of the peptide but also for its proapoptotic action on cancer cells, is more ordered than its C-terminal counterpart and contains a site with a propensity for α-helical secondary structure.

Mapping Hidden Residual Structure within the Myc bHLH-LZ Domain Using Chemical Denaturant Titration

Intrinsically disordered proteins (IDPs) underpin biological regulation and hence are highly desirable drug-development targets. NMR is normally the tool of choice for studying the conformational preferences of IDPs, but the association of regions with residual structure into partially collapsed states can lead to poor spectral quality. The bHLH-LZ domain of the oncoprotein Myc is an archetypal example of such behavior. To circumvent spectral limitations, we apply chemical denaturant titration (CDT)-NMR, which exploits the predictable manner in which chemical denaturants disrupt residual structure and the rapid exchange between conformers in IDP ensembles. The secondary structure propensities and tertiary interactions of Myc are determined for all bHLH-LZ residues, including those with poor NMR properties under native conditions. This reveals conformations that are not predictable using existing crystal structures. The CDT-NMR method also maps sites perturbed by the prototype Myc inhibitor, 10058-F4, to areas of residual structure.

Functional interplay between DEAD-box RNA helicases Ded1 and Dbp1 in preinitiation complex attachment and scanning on structured mRNAs in vivo.

RNA structures that impede ribosome binding or subsequent scanning of the 5'-untranslated region (5'-UTR) for the AUG initiation codon reduce translation efficiency. Yeast DEAD-box RNA helicase Ded1 appears to promote translation by resolving 5'-UTR structures, but whether its paralog, Dbp1, performs similar functions is unknown. Furthermore, direct in vivo evidence was lacking that Ded1 or Dbp1 resolves 5'-UTR structures that impede attachment of the 43S preinitiation complex (PIC) or scanning. Here, profiling of translating 80S ribosomes reveals that the translational efficiencies of many more mRNAs are reduced in a ded1-ts dbp1Δ double mutant versus either single mutant, becoming highly dependent on Dbp1 or Ded1 only when the other helicase is impaired. Such 'conditionally hyperdependent' mRNAs contain unusually long 5'-UTRs with heightened propensity for secondary structure and longer transcript lengths. Consistently, overexpressing Dbp1 in ded1 cells improves the translation of many such Ded1-hyperdependent mRNAs. Importantly, Dbp1 mimics Ded1 in conferring greater acceleration of 48S PIC assembly in a purified system on mRNAs harboring structured 5'-UTRs. Profiling 40S initiation complexes in ded1 and dbp1 mutants provides direct evidence that Ded1 and Dbp1 cooperate to stimulate both PIC attachment and scanning on many Ded1/Dbp1-hyperdependent mRNAs in vivo.