Core Residues Research Articles

Dynamics of the bypass polymerase, dinB homolog.

The dinB homolog (Dbh) polymerase is a member of the Y-family of translesion DNA polymerases, able to synthesize DNA across damaged DNA template. Since Dbh is derived from the thermophilic archaeon Sulfolobus acidocaldarius, it is capable of functioning at high temperatures where mesophilic proteins would be denatured. Here we use NMR and Circular Dichroism (CD) to compare protein stability at a wide range of temperatures. 1H-15N HSQC spectra show a change of tertiary structures starting at 25 °C and progression at colder temperatures (15, 5 °C). CD show retentions of secondary structure with a slight decrease percent helicity and an increase in percent beta strand as temperature decrease. We have analyzed hydrogen exchange measurement at 35°C and 50°C to investigate the folding dynamics of Dbh. We observe the moderate modulation of dynamics upon heating. Core residues in the palm and little finger domain, especially those in the interior of the β-sheets and interior face of α-helices, are key in maintaining the stability of Dbh; with the palm providing a stable platform for the active site residues. We report on regions of the protein that become more dynamic as the temperature is increased toward the functional temperature.

Histone H3 core domain in chromatin with different DNA linker lengths studied by 1H-Detected solid-state NMR spectroscopy.

Chromatin, a dynamic protein-DNA complex that regulates eukaryotic genome accessibility and essential functions, is composed of nucleosomes connected by linker DNA with each nucleosome consisting of DNA wrapped around an octamer of histones H2A, H2B, H3 and H4. Magic angle spinning solid-state nuclear magnetic resonance (NMR) spectroscopy can yield unique insights into histone structure and dynamics in condensed nucleosomes and nucleosome arrays representative of chromatin at physiological concentrations. Recently we used J-coupling-based solid-state NMR methods to investigate with residue-specific resolution the conformational dynamics of histone H3 N-terminal tails in 16-mer nucleosome arrays containing 15, 30 or 60 bp DNA linkers. Here, we probe the H3 core domain in the 16-mer arrays as a function of DNA linker length via dipolar coupling-based 1H-detected solid-state NMR techniques. Specifically, we established nearly complete assignments of backbone chemical shifts for H3 core residues in arrays with 15-60 bp DNA linkers reconstituted with 2H,13C,15N-labeled H3. Overall, these chemical shifts were similar irrespective of the DNA linker length indicating no major changes in H3 core conformation. Notably, however, multiple residues at the H3-nucleosomal DNA interface in arrays with 15 bp DNA linkers exhibited relatively pronounced differences in chemical shifts and line broadening compared to arrays with 30 and 60 bp linkers. These findings are consistent with increased heterogeneity in nucleosome packing and structural strain within arrays containing short DNA linkers that likely leads to side-chains of these interfacial residues experiencing alternate conformations or shifts in their rotamer populations relative to arrays with the longer DNA linkers.

The residual structure of acid-denatured β2 -microglobulin is relevant to an ordered fibril morphology.

β2 -Microglobulin (β2m) forms amyloid fibrils in vitro under acidic conditions. Under these conditions, the residual structure of acid-denatured β2m is relevant to seeding and fibril extension processes. Disulfide (SS) bond-oxidized β2m has been shown to form rigid, ordered fibrils, whereas SS bond-reduced β2m forms curvy, less-ordered fibrils. These findings suggest that the presence of an SS bond affects the residual structure of the monomer, which subsequently influences the fibril morphology. To clarify this process, we herein performed NMR experiments. The results obtained revealed that oxidized β2m contained a residual structure throughout the molecule, including the N- and C-termini, whereas the residual structure of the reduced form was localized and other regions had a random coil structure. The range of the residual structure in the oxidized form was wider than that of the fibril core. These results indicate that acid-denatured β2m has variable conformations. Most conformations in the ensemble cannot participate in fibril formation because their core residues are hidden by residual structures. However, when hydrophobic residues are exposed, polypeptides competently form an ordered fibril. This conformational selection phase may be needed for the ordered assembly of amyloid fibrils.

Sequence-based Prediction of the Cellular Toxicity Associated with Amyloid Aggregation within Protein Condensates.

Various neurological dysfunctions are associated with cytotoxic amyloid-containing aggregates formed through the irreversible maturation of protein condensates generated by phase separation. Here, we investigate the amino acid code for this cytotoxicity using TDP-43 deep-sequencing data. Within the droplet landscape framework, we analyze the impact of mutations in the amyloid core, aggregation hot-spot, and droplet-promoting residues on TDP-43 cytotoxicity. Our analysis suggests that TDP-43 mutations associated with low cytotoxicity moderately decrease the probability of droplet formation while increasing the probability of multimodal binding. These mutations promote both ordered and disordered binding modes, thus facilitating the conversion between the droplet and amyloid states. Based on this understanding, we develop an extension of the FuzDrop method for the sequence-based prediction of the cytotoxicity of aging condensates and test it over 20,000 TDP-43 variants. Our analysis provides insight into the amino acid code that regulates the cytotoxicity associated with the maturation of liquid-like condensates into amyloid-containing aggregates, suggesting that, at least in the case of TDP-43, mutations that promote aggregation tend to decrease cytotoxicity, while those that promote droplet formation tend to increase cytotoxicity.

Physicochemical and antiproliferative characteristics of RNA and DNA sequence-related G-quadruplexes.

In this article the physicochemical and biological properties of sequence-related G-quadruplex forming oligonucleotides in RNA and DNA series are analyzed and compared. The intermolecular G-quadruplexes vary in loop length, number of G-tetrads and homogeneity of the core. Our studies show that even slight variations in sequence initiate certain changes of G-quadruplex properties. DNA G-quadruplexes are less thermally stable than their RNA counterparts, more topologically diversified and are better candidates as inhibitors of cancer cells proliferation. The most efficient antiproliferative activity within the studied group of molecules was observed for two DNA G-quadruplexes with unperturbed core and lower content of thymidine residues within the loops leading to reduction of cells viability up to 65% and 33% for HeLa and MCF-7 cell lines, respectively.

Synthesis and biological evaluation of substituted aurone derivatives as potential tyrosinase inhibitors: in vitro, kinetic, QSAR, docking and drug-likeness studies

Tyrosinase enzyme plays an essential role in melanin biosynthesis and enzymatic browning of fruits and vegetables. To discover potent tyrosinase inhibitors, the present studies were undertaken. In this context, synthetic aurone derivatives 26–50 were designed, synthesized, and structurally elucidated by various spectroscopic techniques including IR, UV, 1H- & 13C-NMR and mass spectrometry. The target compounds 26–50 were screened for their anti-tyrosinase inhibitory potential, and thus kinetic mechanism was analyzed by Lineweaver-Burk plots. All target compounds exhibited good to excellent IC50 values in the range of 7.12 ± 0.32 μM to 66.82 ± 2.44 μM. These synthesized aurone derivatives were found as potent tyrosinase inhibitors relative to the standard kojic acid (IC50 = 16.69 ± 2.81 μM) and the compound 39 inhibited tyrosinase non-competitively (Ki = 11.8 μM) by forming an enzyme-inhibitor complex. The binding modes of these molecules were ascribed through molecular docking studies against tyrosinase protein (PDB ID: 2Y9X). The quantitative structure-activity relationship studies displayed a good correlation between 26–50 structures and their anti-tyrosinase activity (IC50) with a correlation coefficient (R2) of 0.9926. The computational studies were coherent with experimental results and these ligands exhibited good binding values against tyrosinase and interacted with core residues of target protein. Moreover, the drug-likeness analysis also showed that some compounds have a linear correlation with Lipinski’s rule of five, indicating good drug-likeness and bioactivity scores for pharmacological targets. Communicated by Ramaswamy H. Sarma

The mechanisms of melanogenesis inhibition by glabridin: molecular docking, PKA/MITF and MAPK/MITF pathways

Glabridin is the main ingredient of hydrophobic fraction in licorice extract and has been shown to have anti-melanogenesis activity in skins. However, the underlying mechanism(s) remain not completely understood. The aim of this study is thus to elucidate the possible mechanisms related to the melanogenesis suppression by glabridin in cultured B16 murine melanoma cells and in UVA radiation induced hyperpigmentation model of BALB/c mice as well. Molecular docking simulations revealed that between catalytic core residues and the compound. The treatment by glabridin significantly downregulated both transcriptional and/or protein expression of melanogenesis-related factors including melanocyte stimulating hormone receptor (MC1R), microphthalmia-associated transcription factor (MITF), tyrosinase (TYR), TYR-related protein-1 (TRP-1) and TRP-2 in B16 cells. Both PKA/MITF and MAPK/MITF signaling pathways were found to be involved in the suppression of melanogenesis by glabridin in B16 cells. Also in vivo glabridin therapy significantly reduced hyperpigmentation, epidermal thickening, roughness and inflammation induced by frequent UVA exposure in mice skins, thus beneficial for skin healthcare. These data further look insights into the molecular mechanisms of melanogenesis suppression by glabridin, rationalizing the application of the natural compound for skin healthcare.

Molecular basis of the anchoring and stabilization of human islet amyloid polypeptide in lipid hydroperoxidized bilayers

The molecular structure of membrane lipids is formed by mono- or polyunsaturations on their aliphatic tails that make them susceptible to oxidation, facilitating the incorporation of hydroperoxide (R-OOH) functional groups. Such groups promote changes in both composition and complexity of the membrane significantly modifying its physicochemical properties. Human Langerhans islets amyloid polypeptide (hIAPP) is the main component of amyloid deposits found in the pancreas of patients with type-2 diabetes (T2D). hIAPP in the presence of membranes with oxidized lipid species accelerates the formation of amyloid fibrils or the formation of intermediate oligomeric structures. However, the molecular bases at the initial stage of the anchoring and stabilization of the hIAPP in a hydroperoxidized membrane are not yet well understood. To shed some light on this matter, in this contribution, three bilayer models were modeled: neutral (POPC), anionic (POPS), and oxidized (POPCOOH), and full atom Molecular Dynamics (MD) simulations were performed. Our results show that the POPCOOH bilayer increases the helicity in hIAPP when compared to POPC or POPS bilayer. The modification in the secondary structure covers the residues of the so-called amyloidogenic core of the hIAPP. Overall, the hydroperoxidation of the neutral lipids modifies both the anchoring and the stabilization of the peptide hIAPP by reducing the random conformations of the peptide and increasing of hydrogen bond population with the hydroperoxidized lipids.

Stability profile of the neuronal SNARE complex reflects its potency to drive fast membrane fusion

Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) form the SNARE complex to mediate most fusion events of the secretory pathway. The neuronal SNARE complex is featured by its high stability and half-zippered conformation required for driving robust and fast synaptic exocytosis. However, these two features seem to be thermodynamically mutually exclusive. In this study, we have employed temperature-dependent disassociation assays and single-molecule Förster resonance energy transfer (FRET) experiments to analyze the stability and conformation of the neuronal SNARE complex. We reclassified the amino acids of the SNARE motif into four sub-groups (core, core-side I and II, and non-contact). Our data showed that the core residues predominantly contribute to the complex stability to meet a basal requirement for SNARE-mediated membrane fusion, while the core-side residues exert an unbalanced effect on the N- and C-half bundle stability that determines the half-zippered conformation of the neuronal SNARE complex, which would accommodate essential regulations by complexins and synaptotagmins for fast Ca2+-triggered membrane fusion. Furthermore, our data confirmed a strong coupling of folding energy between the N- and C-half assembly of the neuronal SNARE complex, which rationalizes the strong potency of the half-zippered conformation to conduct robust and fast fusion. Overall, these results uncovered that the stability profile of the neuronal SNARE complex reflects its potency to drive fast and robust membrane fusion. Based on these results, we also developed a new parameter, the stability factor (Fs), to characterize the overall stability of the neuronal SNARE complex and resolved a linear correlation between the stability and inter-residue coulombic interactions of the neuronal SNARE complex, which would help rationally design artificial SNARE complexes and remold functional SNARE complexes with desirable stability.

A fine balance of hydrophobic-electrostatic communication pathways in a pH-switching protein

Allostery is the phenomenon of coupling between distal binding sites in a protein. Such coupling is at the crux of protein function and regulation in a myriad of scenarios, yet determining the molecular mechanisms of coupling networks in proteins remains a major challenge. Here, we report mechanisms governing pH-dependent myristoyl switching in monomeric hisactophilin, whereby the myristoyl moves between a sequestered state, i.e., buried within the core of the protein, to an accessible state, in which the myristoyl has increased accessibility for membrane binding. Measurements of the pH and temperature dependence of amide chemical shifts reveal protein local structural stability and conformational heterogeneity that accompany switching. An analysis of these measurements using a thermodynamic cycle framework shows that myristoyl-proton coupling at the single-residue level exists in a fine balance and extends throughout the protein. Strikingly, small changes in the stereochemistry or size of core and surface hydrophobic residues by point mutations readily break, restore, or tune myristoyl switch energetics. Synthesizing the experimental results with those of molecular dynamics simulations illuminates atomistic details of coupling throughout the protein, featuring a large network of hydrophobic interactions that work in concert with key electrostatic interactions. The simulations were critical for discerning which of the many ionizable residues in hisactophilin are important for switching and identifying the contributions of nonnative interactions in switching. The strategy of using temperature-dependent NMR presented here offers a powerful, widely applicable way to elucidate the molecular mechanisms of allostery in proteins at high resolution.

Minimal Increments of Hydrophobic Collapse within the N-Terminus of the Neuropeptide Galanin.

The neuropeptide galanin has a 35-year history as an intriguing target in drug design owing to its implication as a potential anticonvulsant and neuronal trophic factor among many other therapeutically interesting functions including analgesia and mood alteration. In this study, we report the structural characterization of three synthetic fragments of the galanin N-terminus in buffered aqueous solution: hGal(2-12)KK, hGal(1-12)KK, and hGal(1-17)KK. High-field two-dimensional 1H-1H nuclear magnetic resonance (NMR) data were acquired for these fragments and used to derive distance restraints. We further utilized modified hydrogen bonding and dihedral restraints to reflect chemical shift patterns in the data, which revealed the signature of a weakly folded helix. Together, these sets of restraints were used to generate NMR structures of all three fragments, which depict a core of hydrophobic residues that cluster together regardless of the presence of a helical structure, and correspond to residues in the N-terminus of galanin that have been previously shown to be critical for binding its receptors. The helical structure only appears following the inclusion of Gly(1) in the sequence, and at longer sequence lengths, unlike many other peptides, the helix does not propagate. Rather, a few turns of poorly ordered helix appear to be a secondary consequence of clusters of hydrophobic sidechains that are conserved across all of the peptides in this study; the helices themselves appear ordered as a consequence of this clustering, and these clusters compare directly to those observed recently to make contacts between galanin and two of its receptor subtypes. Collapsed hydrophobic residues therefore organize and compose the functional core of human galanin and raise interesting questions about the nature of the conformational order in ligands that bind cell surface receptors.

Markov State Modeling Analysis Captures Changes in the Temperature-Sensitive N-Terminal and β-Turn Regions of the p53 DNA-Binding Domain.

The transcription factor p53 is one of the most widely studied cancer proteins. Its temperature-sensitive nature suggests reduction in functionality at physiological temperatures. Temperature-induced conformational variations and their impact on its functional ability still remain unexplored. A total of 20.8 μs molecular dynamics simulations of wildtype p53 in the apo and the DNA-bound states have been performed at 300 K and 310 K. Further, Markov State Modeling (MSM) analyses were performed, considering Cα-Cα distances as reaction coordinates. Filtering of these distances based on correlation with the time-independent components (tICs) resulted in 16 and 32 distances for apo and DNA-bound systems, respectively. Individual MSM analyses using these filtered distances were performed for both p53 systems. These Cα-Cα residue pairs belonged to the N-terminal, S6/7 β-turn, loop L2, loop L3, and hydrophobic core residues. At physiological temperatures, apo-p53 exhibits exposure of its hydrophobic core, where the temperature-sensitive hotspot residues were also located. This exposure was the result of the S6/7 β-turn and N-terminal moving apart. In the DNA-bound p53 system, loop L1 attains an open conformation at physiological temperatures, which weakens the DNA binding. It is already known that p53 mutants that lack DNA binding also tend to show similar conformational variations. The S6/7 β-turn along with the already known functionally important loop L2 may pose as regions to be targeted to overcome the loss in binding of temperature-sensitive wildtype p53. Rescue strategies directed toward these temperature-sensitive regions may be useful to recuperate its strong binding at physiological temperatures.

N-Terminus-Mediated Solution Structure of Dimerization Domain of PRC1.

Simple SummaryThe solution structure of the N-terminal domain of Protein Regulator of Cytokinesis 1 (PRC1) was determined, and compared with the previously published crystal structure, significant differences were found. Extensive analyses were carried out to find the true reason for the differences between the solution and crystal structures, we discovered that this might be related to the conformation of residue M1, which is buried in the protein core of the solution structure, while situated outside of the hydrophobic core in the crystal structure. In this study, we have carried out a series of examinations using various methods and confirmed that the N terminal conformation is the key point in better describing the structure of PRC1 dimerization domain under solution conditions.Microtubule-associated proteins (MAPs) are essential for the accurate division of a cell into two daughter cells. These proteins target specific microtubules to be incorporated into the spindle midzone, which comprises a special array of microtubules that initiate cytokinesis during anaphase. A representative member of the MAPs is Protein Regulator of Cytokinesis 1 (PRC1), which self-multimerizes to cross-link microtubules, the malfunction of which might result in cancerous cells. The importance of PRC1 multimerization makes it a popular target for structural studies. The available crystal structure of PRC1 has low resolution (>3 Å) and accuracy, limiting a better understanding of the structure-related functions of PRC1. Therefore, we used NMR spectroscopy to better determine the structure of the dimerization domain of PRC1. The NMR structure shows that the PRC1 N terminus is crucial to the overall structure integrity, but the crystal structure bespeaks otherwise. We systematically addressed the role of the N terminus by generating a series of mutants in which N-terminal residues methionine (Met1) and arginine (Arg2) were either deleted, extended or substituted with other rationally selected amino acids. Each mutant was subsequently analyzed by NMR spectroscopy or fluorescence thermal shift assays for its structural or thermal stability; we found that N-terminal perturbations indeed affected the overall protein structure and that the solution structure better reflects the conformation of PRC1 under solution conditions. These results reveal that the structure of PRC1 is governed by its N terminus through hydrophobic interactions with other core residues, such hitherto unidentified N-terminal conformations might shed light on the structure–function relationships of PRC1 or other proteins. Therefore, our study is of major importance in terms of identifying a novel structural feature and can further the progress of protein folding and protein engineering.

Molecular insights into the binding model and response mechanisms of triclosan with lysozyme

Triclosan is an emerging pollutant, which belongs to Pharmaceutical and Personal Care Products (PPCPs). The high exposure risk of triclosan is triggered by its high persistence and accumulation levels. Accordingly, we investigated the interaction mechanism of triclosan with lysozyme by biophysical methods, including multiple spectroscopic techniques, enzyme activity measurements, isothermal titration calorimetry (ITC) and molecular docking. The results indicated that the binding of triclosan and lysozyme contributed to the fluorescence sensitization and repositioned the fluorophore, such as Trp108, in a weaker polar microenvironment. Meanwhile, the binding of lysozyme and triclosan broke the secondary structure of lysozyme, as well as induced a small amount of α-Helix structure. The local surrounding of amide bond (CO) in lysozyme was altered and the content of α-Helix structure increased, suggesting that the protein scaffold and polypeptide chain of lysozyme were stretched and a lysozyme agglomerated with a larger particle size. In addition, the thermodynamic results indicated that triclosan docked to lysozyme at relatively moderate affinity, relying on hydrophobic forces, which was an endothermic reaction. The LysozymeTriclosan system became more disordered, indicating that the conformation of lysozyme was broken. Molecular simulations demonstrated that the docking pocket of lysozyme and triclosan contained 8 residues, consisting of Ala 107, Arg 112, Trp 108, Asn 59, Ile 98, Gln 57, Val 109 and Asp 52, which was core residues of lysozyme activity. Triclosan at high concentrations could occupy the native binding site of substrate and disturb Asp 52 residue, showing the antagonism with natural substrate, which inhibited the activity of lysozyme. In summary, this study elucidated the structural and functional responses of lysozyme after triclosan exposure simulated physiological conditions in vitro, as well as the interaction mode between triclosan and lysozyme.

Metagenomic mining and structure-function studies of a hyper-thermostable cellobiohydrolase from hot spring sediment

Enzymatic breakdown is an attractive cellulose utilisation method with a low environmental load. Its high temperature operation could promote saccharification and lower contamination risk. Here we report a hyper-thermostable cellobiohydrolase (CBH), named HmCel6A and its variant HmCel6A-3SNP that were isolated metagenomically from hot spring sediments and expressed in Escherichia coli. They are classified into glycoside hydrolases family 6 (GH6). HmCel6A-3SNP had three amino acid replacements to HmCel6A (P88S/L230F/F414S) and the optimum temperature at 95 °C, while HmCel6A did it at 75 °C. Crystal structure showed conserved features among GH6, a (β/α)8-barrel core and catalytic residues, and resembles TfCel6B, a bacterial CBH II of Thermobifida fusca, that had optimum temperature at 60 °C. From structure-function studies, we discuss unique structural features that allow the enzyme to reach its high thermostability level, such as abundance of hydrophobic and charge-charge interactions, characteristic metal bindings and disulphide bonds. Moreover, structure and surface plasmon resonance analysis with oligosaccharides suggested that the contribution of an additional tryptophan located at the tunnel entrance could aid in substrate recognition and thermostability. These results may help to design efficient enzymes and saccharification methods for cellulose working at high temperatures.

Self-Fluorescent Lone Tryptophan Nanoparticles as Theranostic Agents Against Alzheimer's Disease.

Aggregation of β-amyloid (Aβ42) peptide in the neural extracellular space leads to cellular dysfunction, resulting in Alzheimer's disease (AD). The hydrophobic core of the amyloidogenic Aβ42 peptide contains aromatic residues that play an important role in the self-assembly and subsequent aggregation of the peptide. Hence, targeting these hydrophobic core residues by potent low molecular agents can be a promising therapeutic approach toward AD. In the current work, we have developed self-fluorescent solo tryptophan nanoparticles (TNPs) as nanotheranostic systems against AD. We demonstrated that TNPs could significantly inhibit as well as disrupt the fibrils formed by both Aβ42 peptide and another reductionist approach-based amyloid model dipeptide, phenylalanine-phenylalanine (FF). More importantly, these nanostructures were nontoxic to neural cells and could protect the neurons from Aβ42 peptide and FF aggregate-induced cytotoxicity. In addition, efficacy studies performed in animal model further revealed that the TNPs could rescue spatial and learning memory in intracerebroventricular streptozotocin-administration-induced AD phenotype in rats. Moreover, our pharmacokinetics study further established the BBB permeability and brain delivery potency of TNPs. The inherent excellent fluorescent properties of these nanoparticles could be exploited further to use them as imaging modalities for tagging and detecting FF and Aβ42 peptide fibrils. Overall, our results clearly illustrated that the solo TNPs could serve as promising nanotheranostic agents for AD therapy.

NMR Structure and Localization of the Host Defense Peptide ThanatinM21F in Zwitterionic Dodecylphosphocholine Micelle: Implications in Antimicrobial and Hemolytic Activity.

Non-hemolytic antimicrobial peptides (AMPs) are vital lead molecules for the designing and development of peptide-based antibiotics. Thanatin a 21-amino acid long single disulfide bonded AMP is known to be highly non-hemolytic with a limited toxicity to human cells and model animals. Thanatin demonstrates a potent antibacterial activity against multidrug-resistant Gram-negative pathogens. A single mutated variant of thanatin replaced last residue Met21 to Phe or thanatin M21F has recently been found to be more active compared to the native peptide. In order to gain mechanistic insights toward bacterial cell lysis versus non-hemolysis, here, we report atomic resolution structure and mode insertion of thanatinM21F reconstituted into zwitterionic detergent micelle by use of solution NMR spectroscopy. The 3D structure of thanatinM21F in DPC micelle is defined by an anti-parallel β-sheet between residues I9-F21 with a central cationic loop, residues N12-R14. PRE NMR studies revealed hydrophobic core residues of thanatinM21F are deeply inserted in the DPC micelle, while residues at the extended N-terminal half of the peptide are appeared to be mostly surface localized. Marked structural differences of thanatin and thanatinM21F in negatively charged LPS and DPC micelle could be correlated with non-hemolytic and antibacterial activity.

Can the jigsaw puzzle model of protein folding re-assemble a hydrophobic core?

According to the "jigsaw puzzle" model of protein folding, the isomorphism between sequence and structure is substantially determined by the specific geometry of side-chain interactions, within the protein interior. In this work, we have attempted to predict the hydrophobic core of cyclophilin (LdCyp) from Leishmania donovani, utilizing a surface complementarity function, which selects for high goodness of fit between hydrophobic side-chain surfaces, rather in the manner of assembling a three-dimensional jigsaw puzzle. The computational core prediction method implemented here has been tried on two distinct scenarios, on the LdCyp polypeptide chain with native non-core residues and all core residues initially set to alanine, on a poly-glycine polypeptide chain. Molecular dynamics simulations appeared to indicate partial destabilization of the two designed sequences. However, experimental characterization of the designed sequences by circular dichroism (CD) spectroscopy and denaturant (GdmCl) induced unfolding, demonstrated disordered proteins. Stepwise reconstruction of the designed cores by cumulative sequential mutations identified the specific mutation (M122L) as primarily responsible for fold collapse and all design objectives were achieved upon rectifying this mutation. In summary, the study demonstrates regions of the core to contain highly specific (jigsaw puzzle-like) interactions sensitive to any perturbations and a predictive algorithm to identify such regions. A mutation within the core has been identified which exercises an inordinate influence on the global fold, reminiscent of metamorphic proteins. In addition, the computational procedure could predict substantial regions of the core (given main-chain coordinates) without any reference to non-core residues.

A sequence space search engine for computational protein design to modulate molecular functionality

De-novo protein design explores the untapped sequence space that is otherwise less discovered during the evolutionary process. This necessitates an efficient sequence space search engine for effective convergence in computational protein design. We propose a greedy simulated annealing-based Monte-Carlo parallel search algorithm for better sequence-structure compatibility probing in protein design. The guidance provided by the evolutionary profile, the greedy approach, and the cooling schedule adopted in the Monte Carlo simulation ensures sufficient exploration and exploitation of the search space leading to faster convergence. On evaluating the proposed algorithm, we find that a dataset of 76 target scaffolds report an average root-mean-square-deviation (RMSD) of 1.07 Å and an average TM-Score of 0.93 with the modeled designed protein sequences. High sequence recapitulation of 48.7% (59.4%) observed in the design sequences for all (hydrophobic) solvent-inaccessible residues again establish the goodness of the proposed algorithm. A high (93.4%) intra-group recapitulation of hydrophobic residues in the solvent-inaccessible region indicates that the proposed protein design algorithm preserves the core residues in the protein and provides alternative residue combinations in the solvent-accessible regions of the target protein. Furthermore, a COFACTOR-based protein functional analysis shows that the design sequences exhibit altered molecular functionality and introduce new molecular functions compared to the target scaffolds. Communicated by Ramaswamy H. Sarma

Influence of the Dynamically Disordered N-Terminal Tail Domain on the Amyloid Core Structure of Human Y145Stop Prion Protein Fibrils.

The Y145Stop mutant of human prion protein (huPrP23-144) is associated with a familial prionopathy and provides a convenient in vitro model for investigating amyloid strains and cross-seeding barriers. huPrP23-144 fibrils feature a compact and relatively rigid parallel in-register β-sheet amyloid core spanning ∼30 C-terminal amino acid residues (∼112–141) and a large ∼90-residue dynamically disordered N-terminal tail domain. Here, we systematically evaluate the influence of this dynamic domain on the structure adopted by the huPrP23-144 amyloid core region, by investigating using magic-angle spinning solid-state nuclear magnetic resonance (NMR) spectroscopy a series of fibril samples formed by huPrP23-144 variants corresponding to deletions of large segments of the N-terminal tail. We find that deletion of the bulk of the N-terminal tail, up to residue 98, yields amyloid fibrils with native-like huPrP23-144 core structure. Interestingly, deletion of additional flexible residues in the stretch 99–106 located outside of the amyloid core yields shorter heterogenous fibrils with fingerprint NMR spectra that are clearly distinct from those for full-length huPrP23-144, suggestive of the onset of perturbations to the native structure and degree of molecular ordering for the core residues. For the deletion variant missing residues 99–106 we show that native huPrP23-144 core structure can be “restored” by seeding the fibril growth with preformed full-length huPrP23-144 fibrils.