Aromatic Cluster Sensor of Protein Folding: Near‐UV Electronic Circular Dichroism Bands Assigned to Fold Compactness

Both far‐ and near‐UV electronic circular dichroism (ECD) spectra have bands sensitive to thermal unfolding of Trp and Tyr residues containing proteins. Beside spectral changes at 222 nm reporting secondary structural variations (far‐UV range), L_b bands (near‐UV range) are applicable as 3D‐fold sensors of protein's core structure. In this study we show that both L_b(Tyr) and L_b(Trp) ECD bands could be used as sensors of fold compactness. ECD is a relative method and thus requires NMR referencing and cross‐validation, also provided here. The ensemble of 204 ECD spectra of Trp‐cage miniproteins is analysed as a training set for “calibrating” Trp?Tyr folded systems of known NMR structure. While in the far‐UV ECD spectra changes are linear as a function of the temperature, near‐UV ECD data indicate a non‐linear and thus, cooperative unfolding mechanism of these proteins. Ensemble of ECD spectra deconvoluted gives both conformational weights and insight to a protein folding?unfolding mechanism. We found that the L_b²⁹³ band is reporting on the 3D‐structure compactness. In addition, the pure near‐UV ECD spectrum of the unfolded state is described here for the first time. Thus, ECD folding information now validated can be applied with confidence in a large thermal window (5≤T≤85 °C) compared to NMR for studying the unfolding of Trp?Tyr residue pairs. In conclusion, folding propensities of important proteins (RNA polymerase II, ubiquitin protein ligase, tryptase‐inhibitor etc.) can now be analysed with higher confidence.     

Probing the conformational features of a phage display polypeptide sequence directed against single-walled carbon nanohorn surfaces

Single-walled carbon nanohorns (SWNHs) are interesting carbon nanostructures that have applications to science and technology. Using M13 phage display technology, polypeptides directed again SWNHs surfaces have been created for a number of nanotechnology and pharmaceutical purposes, yet the molecular mechanism of polypeptide sequence interaction and binding to SWNHs surfaces is not known. Recently, we identified a linear 12-AA M13 phage pIII sequence, NH-12-5-2 (DYFSSPYYEQLF), that binds with high affinity to SWNHs surfaces. To probe the structure of this pIII tail polypeptide further, we investigated the conformation of a model peptide representing the 12 AA NH-12-5-2 sequence. At neutral pH, the NH-12-5-2 model polypeptide is conformationally labile and exhibits two-state conformational exchange involving the D1-S5 N-terminal segment. Simultaneous with this conformational exchange process is the observation that the P6 residue exhibits imido ring conformational variation. In the presence of the structure-stabilizing solvent, TFE, or at pH 2.5, both the exchange process and Pro ring motion phenomena disappear, indicating that the structure of this peptide sequence can be stabilized by extrinsic factors. Interestingly, we observe NMR parameters (ROEs, (3)J coupling constants) for NH-12-5-2 in 90% v/v TFE that are consistent with the presence of a partial helical structure, similar to what was observed at low pH in our earlier CD experiments. We conclude that the NH-12-5-2 model polypeptide sequence possesses an inherent conformational instability that involves the D1-S5 sequence segment and the P6 residue but that this instability can be offset by extrinsic factors (e.g., charge neutralization, imido ring interconversion, and hydrophobic-hydrophobic interactions). These nonbonding interactions may play a role in the recognition and binding of this phage sequence region to SWNHs surfaces.     

蛋白质折叠类型的分类建模与识别

蛋白质的氨基酸序列如何决定空间结构是当今生命科学研究中的核心问题之一. 折叠类型反映了蛋白质核心结构的拓扑模式, 折叠识别是蛋白质序列-结构研究的重要内容. 我们以占Astral 1.65序列数据库中α, β和α/β三类蛋白质总量41.8%的36个无法独立建模的折叠类型为研究对象, 选取其中序列一致性小于25%的样本作为训练集, 以均方根偏差(RMSD)为指标分别进行系统聚类, 生成若干折叠子类, 并对各子类建立基于多结构比对算法(MUSTANG)结构比对的概形隐马尔科夫模型(profile-HMM). 将Astral 1.65中序列一致性小于95%的9505个样本作为检验集, 36个折叠类型的平均识别敏感性为90%, 特异性为99%, 马修斯相关系数(MCC)为0.95. 结果表明: 对于成员较多, 无法建立统一模型的折叠类型, 基于RMSD的系统分类建模均可实现较高准确率的识别, 为蛋白质折叠识别拓展了新的方法和思路, 为进一步研究奠定了基础.     

OB-fold: growing bigger with functional consistency

It was predicted that the folding space for various protein sequences is restricted and a maximum of 1000 protein folds could be expected. Although, there were about 648 folds identified, general functional features of individual folds is not thoroughly studied. We selected OB-fold, which is supposed to be an oligonucleotide and oligosaccharide binding fold to study the general functional features. OB-fold is a small beta-barrel fold formed from 5 strands connected by modulating loops. We observed consistently 2 or 3 loops on the same face of barrel acting as clamps to bind to their ligands. Depending on the ligand, which could be a single or double stranded DNA/RNA or an oligosaccharide, and their conformational properties the loops change in length and sequence to accommodate various ligands. Different classes of OB-folded proteins were analyzed and found that the functional features are retained in spite of negligible sequence homology among various proteins studied.     

Quantitative approach to conformations of proteins

A quantitative conformational theory of proteins is developed that enables one to predict the native structure of a protein from its amino acid sequence. The theory is based on the following principles: (1) the spatial structure and conformational properties of a protein are predetermined by its amino acid sequence; (2) the native conformation of a protein corresponds to the free energy minimum; (3) all interactions within a protein molecule are specified as short-, mediumy-, and long-range types, interactions of different types being consistent with each other. The role of the short-, medium-, and long-range interactions in the spatial organization of a protein globule is discussed, and a step-by-step analysis of amino acid sequences with gradually increasing lengths is presented. The proposed theory is based on a semiempirical computational method that involves quantitative evaluation of all pairwise atomic interactions within a protein molecule in an aqueous medium. Examples illustrating the suggested approach are presented.     

Structural and Functional Analysis of the Pyridoxal Phosphate Homeostasis Protein YggS from Fusobacterium nucleatum

Pyridoxal 5′-phosphate (PLP) is the active form of vitamin B6, but it is highly reactive and poisonous in its free form. YggS is a PLP-binding protein found in bacteria and humans that mediates PLP homeostasis by delivering PLP to target enzymes or by performing a protective function. Several biochemical and structural studies of YggS have been reported, but the mechanism by which YggS recognizes PLP has not been fully elucidated. Here, we report a functional and structural analysis of YggS from Fusobacterium nucleatum (FnYggS). The PLP molecule could bind to native FnYggS, but no PLP binding was observed for selenomethionine (SeMet)-derivatized FnYggS. The crystal structure of FnYggS showed a type III TIM barrel fold, exhibiting structural homology with several other PLP-dependent enzymes. Although FnYggS exhibited low (<35%) amino acid sequence similarity with previously studied YggS proteins, its overall structure and PLP-binding site were highly conserved. In the PLP-binding site of FnYggS, the sulfate ion was coordinated by the conserved residues Ser201, Gly218, and Thr219, which were positioned to provide the binding moiety for the phosphate group of PLP. The mutagenesis study showed that the conserved Ser201 residue in FnYggS was the key residue for PLP binding. These results will expand the knowledge of the molecular properties and function of the YggS family.     

Structure and photoreaction of photoactive yellow protein, a structural prototype of the PAS domain superfamily

Photoactive yellow protein (PYP) is a water-soluble photosensor protein found in purple photosynthetic bacteria. Unlike bacterial rhodopsins, photosensor proteins composed of seven transmembrane helices and a retinal chromophore in halophilic archaebacteria, PYP is a highly soluble globular protein. The alpha/beta fold structure of PYP is a structural prototype of the PAS domain superfamily, many members of which function as sensors for various kinds of stimuli. To absorb a photon in the visible region, PYP has a p-coumaric acid chromophore binding to the cysteine residue via a thioester bond. It exists in a deprotonated trans form in the dark. The primary photochemical event is photo-isomerization of the chromophore from trans to cis form. The twisted cis chromophore in early intermediates is relaxed and finally protonated. Consequently, the chromophore becomes electrostatically neutral and rearrangement of the hydrogen-bonding network triggers overall structural change of the protein moiety, in which local conformational change around the chromophore is propagated to the N-terminal region. Thus, it is an ideal model for protein conformational changes that result in functional change, responding to stimuli and expressing physiological activity. In this paper, recent progress in investigation of the photoresponse of PYP is reviewed.     

Protein model refinement using structural fragment tessellation

A method is described for the refinement of rough protein models based on finding a selection of structural fragments that match the model. Unlike most fragment-based methods, these are not necessarily contiguous in the sequence and form a tiling (tessellation) that covers most of the structure. The residue positions of the fragments are then used as a target for the model atoms to generate a revised model which is used as the basis of a subsequent pattern definition and search. The method was shown to improve the recognition of the native fold in a series of decoys largely as a result of improved secondary structure representation.     

Ribonuclease T1: Structure,Function, and Stability

Proteins carry out the most important and difficult tasks in all living organisms. To do so, they must often interact specifically with other small and large molecules. This requires that they fold to a globular conformation with a unique active site that is used for the specific interaction. Consequently, protein folding can be regarded as the “secret of life”. Biochemists and chemists have a great interest in elucidating the mechanism by which proteins fold and in predicting the folded conformation and its stability given just the amino acid sequence. This challenge is sometimes called the “protein folding problem”. The ability to construct proteins differing in sequence by one or more amino acids and to analyze their three-dimensional structures by X-ray crystallography and NMR spectroscopy is a powerful tool for investigating the conformational stability and folding of proteins. Several proteins are now under intensive study by this approach. One of these is ribonuclease T1.     

Conformational flexibility in the peripheral site of Torpedo californica acetylcholinesterase revealed by the complex structure with a bifunctional inhibitor

The X-ray crystallographic structure of Torpedo californica acetylcholinesterase (TcAChE) in complex with the bifunctional inhibitor NF595, a potentially new anti-Alzheimer drug, has been solved. For the first time in TcAChE, a major conformational change in the peripheral-site tryptophan residue is observed upon complexation. The observed conformational flexibility highlights the dynamic nature of protein structures and is of importance for structure-based drug design.     

Protein secondary structure preferencs

This paper addresses the question to what extent steric properties of sequence neighbors effect the preferences of an amino acid residue to assume the-helical or some other secondary structure conformation. We find that an amino acid has increased tendency to be in-helical conformation when its sequence neighbors are bulky. This result is an outcome of our automated method for finding conformational preferences as functions of physical parameters important for protein folding. The steric environment for a given residue in a protein is defined as an average of water-accessible surface areas of its primary structure neighbors in extended conformation for model tripeptides. For all amino acids, including non-helix formers like glycine and arginine, the preference for the helical structure increases if their primary structure neighbors form a larger steric environment.     

Selecting folded proteins from a library of secondary structural elements

A protein evolution strategy is described by which double-stranded DNA fragments encoding defined Escherichia coli protein secondary structural elements (alpha-helices, beta-strands, and loops) are assembled semirandomly into sequences comprised of as many as 800 amino acid residues. A library of novel polypeptides generated from this system was inserted into an enhanced green fluorescent protein (EGFP) fusion vector. Library members were screened by fluorescence activated cell sorting (FACS) to identify those polypeptides that fold into soluble, stable structures in vivo that comprised a subset of shorter sequences ( approximately 60 to 100 residues) from the semirandom sequence library. Approximately 108 clones were screened by FACS, a set of 1149 high fluorescence colonies were characterized by dPCR, and four soluble clones with varying amounts of secondary structure were identified. One of these is highly homologous to a domain of aspartate racemase from a marine bacterium (Polaromonas sp.) but is not homologous to any E. coli protein sequence. Several other selected polypeptides have no global sequence homology to any known protein but show significant alpha-helical content, limited dispersion in 1D nuclear magnetic resonance spectra, pH sensitive ANS binding and reversible folding into soluble structures. These results demonstrate that this strategy can generate novel polypeptide sequences containing secondary structure.     

Transform and relax sampling for highly anisotropic systems: application to protein domain motion and folding

Transform and relax sampling (TRS) is proposed as a conformational sampling method to enhance "soft" fluctuation in highly anisotropic systems using molecular dynamics simulation. This method consists of three stages; transform, relax, and sampling. In the transform stage, molecular dynamics simulation is performed with randomly assigned force bias to enhance the fluctuations along relatively soft collective movements, as expected from the linear response theory. After relaxing the heated system to equilibrium without force bias in the relax stage, Monte Carlo-type determination is made as to whether the generated state is accepted or not. The sampling stage is then conducted for conformational sampling by conventional molecular dynamics simulation. TRS is first applied for the idealized multidimensional double-well C(α) model to mimic protein open-close transition. Subsequently, it is applied to three different all-atom protein systems in an explicit solvent model; T4 lysozyme, glutamine binding protein, and a mini-protein chignolin. Investigation of structural variations in the hinge angle of T4 lysozyme in crystals is demonstrated by TRS. The liganded close structure of the glutamine binding protein is sampled starting from the unliganded open form. Chignolin is shown to fold into a native structure multiple times starting from highly extended structures within 100 ns. It is concluded that TRS sampled a reasonable conformational space within a relatively short simulation time in these cases. Possible future extensions of TRS are also discussed.     

Subtle pH differences trigger single residue motions for moderating conformations of calmodulin

This study reveals the essence of ligand recognition mechanisms by which calmodulin (CaM) controls a variety of Ca(2+) signaling processes. We study eight forms of calcium-loaded CaM each with distinct conformational states. Reducing the structure to two degrees of freedom conveniently describes main features of the conformational changes of CaM via simultaneous twist-bend motions of the two lobes. We utilize perturbation-response scanning (PRS) technique, coupled with molecular dynamics simulations. PRS is based on linear response theory, comprising sequential application of directed forces on selected residues followed by recording the resulting protein coordinates. We analyze directional preferences of the perturbations and resulting conformational changes. Manipulation of a single residue reproduces the structural change more effectively than that of single/pairs/triplets of collective modes of motion. Our findings also give information on how the flexible linker acts as a transducer of binding information to distant parts of the protein. Furthermore, by perturbing residue E31 located in one of the EF hand motifs in a specific direction, it is possible to induce conformational change relevant to five target structures. Independently, using four different pK(a) calculation strategies, we find this particular residue to be the charged residue (out of a total of 52), whose ionization state is most sensitive to subtle pH variations in the physiological range. It is plausible that at relatively low pH, CaM structure is less flexible. By gaining charged states at specific sites at a pH value around 7, such as E31 found in the present study, local conformational changes in the protein will lead to shifts in the energy landscape, paving the way to other conformational states. These findings are in accordance with Fluorescence Resonance Energy Transfer (FRET) measured shifts in conformational distributions towards more compact forms with decreased pH. They also corroborate mutational studies and proteolysis results which point to the significant role of E31 in CaM dynamics.     

A protein fold classifier formed by fusing different modes of pseudo amino acid composition via PSSM

Protein function is related to its chemical reaction to the surrounding environment including other proteins. On the other hand, this depends on the spatial shape and tertiary structure of protein and folding of its constituent components in space. The correct identification of protein domain fold solely using extracted information from protein sequence is a complicated and controversial task in the current computational biology. In this article a combined classifier based on the information content of extracted features from the primary structure of protein has been introduced to face this challenging problem. In the first stage of our proposed two-tier architecture, there are several classifiers each of which is trained with a different sequence based feature vector. Apart from the application of the predicted secondary structure, hydrophobicity, van der Waals volume, polarity, polarizability, and different dimensions of pseudo-amino acid composition vectors in similar studies, the position specific scoring matrix (PSSM) has also been used to improve the correct classification rate (CCR) in this study. Using K-fold cross validation on training dataset related to 27 famous folds of SCOP, the 28 dimensional probability output vector from each evidence theoretic K-NN classifier is used to determine the information content or expertness of corresponding feature for discrimination in each fold class. In the second stage, the outputs of classifiers for test dataset are fused using Sugeno fuzzy integral operator to make better decision for target fold class. The expertness factor of each classifier in each fold class has been used to calculate the fuzzy integral operator weights. Results make it possible to provide deeper interpretation about the effectiveness of each feature for discrimination in target classes for query proteins.     

Single molecule conformational dynamics of adenylate kinase: energy landscape, structural correlations, and transition state ensembles

We developed a coarse grained two-well model to study the single molecule protein conformational dynamics in microscopic detail at the residue level, overcoming the often encountered computational bottleneck. In particular, we explored the underlying conformational energy landscape of adenylate kinase, a crucial protein for signal transduction in the cell, and identified two major kinetic pathways for the conformational switch between open and closed states through either the intermediate state or the transient state. Based on the parameters fitted to the room-temperature experimental data, we predicted open and closed kinetic rates at the whole temperature ranges from 10 to 50 degrees C, which agree well with the experimental turnover numbers. After uncovering the underlying mechanism for conformational dynamics and exploring the structural correlations, we found the crucial dynamical interplay between the nucleoside monophosphate binding domain (NMP) and the ATP-binding domain (LID) in controlling the conformational switch. The key residues and contacts responsible for the conformational transitions are identified by following the time evolution of the two-dimensional spatial contact maps and characterizing the transition state as well as intermediate structure ensembles through phi value analysis. Our model provides a general framework to study the conformational dynamics of biomolecules and can be applied to many other systems.     

The Role of Aromatic Residues in Stabilizing the Secondary and Tertiary Structure of Avian Pancreatic Polypeptide

Avian Pancreatic Polypeptide is a 36 residue protein that exhibits a tertiary fold. Results of previous experimental and computational studies indicate that the structure of aPP is stabilized more by non-bonded interactions than by the hydrophobic effect. Aromatic residues are known to participate in a variety of long range non-bonded interactions, with both backbone atoms and the atoms of other side-chains, which could be responsible, in part, for the stability of both the local secondary structure and the tertiary fold. The effect of these aromatic interactions on the stability of aPP was calculated using BHandHLYP/cc-pVTZ. Aromatic residues were shown to participate in multiple hydrogen bonded and weakly polar interactions in the secondary structure. The energies of the weakly polar interactions are comparable with those of hydrogen bonds. Aromatic residues were also shown to participate in multiple weakly polar interactions across the tertiary fold, again with energies similar to those of hydrogen bonds.     

Affinity-Based Screening Technology and HCV Drug Discovery

NS5A is one of the non-structural gene products encoded by Hepatitis C virus (HCV) and related viruses that are essential for viral replication. The amino acid sequence of NS5A is conserved between different HCV genotypes and the primary amino acid sequence of NS5A is unique to HCV and closely related viruses. Importantly, NS5A is unrelated to any human protein. This indicates that drugs designed to block the actions of NS5A could inhibit the replication of HCV without showing toxic side effects in human host cells, thus making NS5A inhibitors ideal anti-viral drugs. However, there are presently no functional assays for this essential viral protein. Therefore, conventional high throughput screening (HTS) approaches can not be used to discover antiviral drugs against NS5A.