Secondary Structure Packing Research Articles

Sequence-Dependent Backbone Dynamics of Intrinsically Disordered Proteins.

For intrinsically disordered proteins (IDPs), a pressing question is how sequence codes for function. Dynamics serves as a crucial link, reminiscent of the role of structure in sequence-function relations of structured proteins. To define general rules governing sequence-dependent backbone dynamics, we carried out long molecular dynamics simulations of eight IDPs. Blocks of residues exhibiting large amplitudes in slow dynamics are rigidified by local inter-residue interactions or secondary structures. A long region or an entire IDP can be slowed down by long-range contacts or secondary-structure packing. On the other hand, glycines promote fast dynamics and either demarcate rigid blocks or facilitate multiple modes of local and long-range inter-residue interactions. The sequence-dependent backbone dynamics endows IDPs with versatile response to binding partners, with some blocks recalcitrant while others readily adapting to intermolecular interactions.

Deep learning geometrical potential for high-accuracy ab initio protein structure prediction

SummaryAb initio protein structure prediction has been vastly boosted by the modeling of inter-residue contact/distance maps in recent years. We developed a new deep learning model, DeepPotential, which accurately predicts the distribution of a complementary set of geometric descriptors including a novel hydrogen-bonding potential defined by C-alpha atom coordinates. On 154 Free-Modeling/Hard targets from the CASP and CAMEO experiments, DeepPotential demonstrated significant advantage on both geometrical feature prediction and full-length structure construction, with Top-L/5 contact accuracy and TM-score of full-length models 4.1% and 6.7% higher than the best of other deep-learning restraint prediction approaches. Detail analyses showed that the major contributions to the TM-score/contact-map improvements come from the employment of multi-tasking network architecture and metagenome-based MSA collection assisted with confidence-based MSA selection, where hydrogen-bonding and inter-residue orientation predictions help improve hydrogen-bonding network and secondary structure packing. These results demonstrated new progress in the deep-learning restraint-guided ab initio protein structure prediction.

Exploring the structural significance of molecular switch mechanism alias motif phosphorylation in Wnt/β-catenin and their crucial role in triple-negative breast cancer

β-Catenin, a key transcriptional factor involved in the canonical Wnt signaling pathway, is regulated by a cascade of phosphorylations and plays a major role in the progression of triple-negative breast cancer (TNBC). However, the phosphorylation induced conformational changes in a β-Catenin is still poorly understood. Hence, we adopted a conventional molecular dynamics approach to study phosphorylations present in a sequence motif Ser 552 675 and Tyr670 of the β-Catenin domain and analyzed in terms of structural transitions, bond formation, and folding-misfolding conformations. Our results unveil the β-Catenin linear motif 549–555 (RRTSMGG) of armadillo repeats domain prefers order to disorder state. In contrast, helix C associated with 670–678 (YKKRLSVEL) motif prefers disorder to order upon phosphorylation of Ser 552 675 and Tyr670. In addition, the crucial secondary structural transition from α-helix to coil induced by phospho Ser552 and phospho Tyr670 of β-Catenin ARM domain connecting helix C modifies conformational diversity and binding affinities of the complex interaction in functional regulation significantly. Moreover, the post phosphorylation disrupted the hydrogen bond interactions (Ser552-Arg549, Arg550-Asp546 and Ser675-Lys672) and abolished the residual alliance with hydrophobic interactions (Tyr670-Leu674) that easily interrupt in secondary structure packing as well as folding conformations connecting ARM and helix C (R10, 12 & R1C) compared to unphosphorylation. Our integrated computational analysis may help in shedding light on understanding the induced folding and unfolding pattern due to motif phosphorylations. Overall, our results provide an atomistic structural description of the way phosphorylation facilitates conformational and dynamic changes in β-Catenin, a fundamental molecular switch mechanism in triple-negative breast cancer pathogenesis.

Structural features and oligomeric nature of human podocin domain

Podocytes are crucial cells of the glomerular filtration unit and plays a vital role at the interface of the blood-urine barrier. Podocyte slit-diaphragm is a modified tight junction that facilitates size and charge-dependent permselectivity. Several proteins including podocin, nephrin, CD2AP, and TRPC6 form a macromolecular assembly and constitute the slit-diaphragm. Podocin is an integral membrane protein attached to the inner membrane of the podocyte via a short transmembrane region (101–125). The cytosolic N- and C-terminus help podocin to attain a hook-like structure. Podocin shares 44% homology with stomatin family proteins and similar to the stomatin proteins, podocin was shown to associate into higher-order oligomers at the site of slit-diaphragm. However, the stoichiometry of the homo-oligomers and how it partakes in the macromolecular assemblies with other slit-diaphragm proteins remains elusive. Here we investigated the oligomeric propensity of a truncated podocin construct (residues:126–350). We show that the podocin domain majorly homo-oligomerizes into a 16-mer. Circular dichroism and fluorescence spectroscopy suggest that the 16-mer oligomer has considerable secondary structure and moderate tertiary packing.

Calcium binding to a disordered domain of a type III-secreted protein from a coral pathogen promotes secondary structure formation and catalytic activity

Strains of the Gram-negative bacterium Vibrio coralliilyticus cause the bleaching of corals due to decomposition of symbiotic microalgae. The V. coralliilyticus strain ATCC BAA-450 (Vc450) encodes a type III secretion system (T3SS). The gene cluster also encodes a protein (locus tag VIC_001052) with sequence homology to the T3SS-secreted nodulation proteins NopE1 and NopE2 of Bradyrhizobium japonicum (USDA110). VIC_001052 has been shown to undergo auto-cleavage in the presence of Ca2+ similar to the NopE proteins. We have studied the hitherto unknown secondary structure, Ca2+-binding affinity and stoichiometry of the “metal ion-inducible autocleavage” (MIIA) domain of VIC_001052 which does not possess a classical Ca2+-binding motif. CD and fluorescence spectroscopy revealed that the MIIA domain is largely intrinsically disordered. Binding of Ca2+ and other di- and trivalent cations induced secondary structure and hydrophobic packing after partial neutralization of the highly negatively charged MIIA domain. Mass spectrometry and isothermal titration calorimetry showed two Ca2+-binding sites which promote structure formation with a total binding enthalpy of −110 kJ mol−1 at a low micromolar Kd. Putative binding motifs were identified by sequence similarity to EF-hand domains and their structure analyzed by molecular dynamics simulations. The stoichiometric Ca2+-dependent induction of structure correlated with catalytic activity and may provide a “host-sensing” mechanism that is shared among pathogens that use a T3SS for efficient secretion of disordered proteins.

Foldamer Tertiary Structure through Sequence-Guided Protein Backbone Alteration.

The prospect of recreating the complex structural hierarchy of protein folding in synthetic oligomers with backbones that are artificial in covalent structure ("foldamers") has long fascinated chemists. Foldamers offer complex functions from biostable scaffolds and have found widespread applications in fields from biomedical to materials science. Most precedent has focused on isolated secondary structures or their assemblies. In considering the goal of complex protein-like tertiary folding patterns, a key barrier became apparent. How does one design a backbone with covalent connectivity and a sequence of side-chain functional groups that will support defined intramolecular packing of multiple artificial secondary structures? Two developments were key to overcoming this challenge. First was the recognition of the power of blending α-amino acid residues with monomers differing in backbone connectivity to create "heterogeneous-backbone" foldamers. Second was the finding that replacing some of the natural α-residues in a biological sequence with artificial-backbone variants can result in a mimic that retains both the fold and function of the native sequence and, in some cases, gains advantageous characteristics. Taken together, these precedents lead to a view of a protein as chemical entity having two orthogonal sequences: a sequence of side-chain functional groups and a separate sequence of backbone units displaying those functional groups. In this Account, we describe our lab's work over the last ∼10 years to leverage the above concept of protein sequence duality in order to develop design principles for constructing heterogeneous-backbone foldamers that adopt complex protein-like tertiary folds. Fundamental to the approach is the utilization of a variety of artificial building blocks (e.g., d-α-residues, Cα-Me-α-residues, N-Me-α-residues, β-residues, γ-residues, δ-residues, polymer segments) in concert, replacing a fraction of α-residues in a given prototype sequence. We provide an overview of the state-of-the-art in terms of design principles for choosing substitutions based on consideration of local secondary structure and retention of key side-chain functional groups. We survey high-resolution structures of backbone-modified proteins to illustrate how diverse artificial moieties are accommodated in tertiary fold contexts. We detail efforts to elucidate how backbone alteration impacts folding thermodynamics and describe how such data informs the development of improved design rules. Collectively, a growing body of results by our lab and others spanning multiple protein systems suggests there is a great deal of plasticity with respect to the backbone chemical structures upon which sequence-encoded tertiary folds can manifest. Moreover, these efforts suggest sequence-guided backbone alteration as a broadly applicable strategy for generating foldamers with complex tertiary folding patterns. We conclude by offering some perspective regarding the near future of this field, in terms of unanswered questions, technological needs, and opportunities for new areas of inquiry.

An Alternatively Packed Dry Molten Globule-like Intermediate in the Native State Ensemble of a Multidomain Protein.

It has been difficult to quantify the degree of side-chain conformational heterogeneity in the native (N) state ensemble of proteins and the relative energetic contributions of the side-chain packing and the hydrophobic effect in protein stability. Here, we show using multiple site-specific spectroscopic probes and tools of thermodynamics that the N state ensemble of a multidomain protein contains an equilibrium intermediate (I) whose interdomain region resembles a dry molten globule. In the I state, a tryptophan residue in the interdomain region is alternatively packed, but its secondary structure and intradomain packing are N-like. The I state also has a larger interdomain distance, but the domain-domain interface is dry and molten. Our results indicate that hydrophobic desolvation and side-chain packing are decoupled during protein folding and that interdomain packing interactions have an important energetic contribution in protein stability. Dynamic interconversion between alternatively packed N-like states could be important for multiple allosteric and ligand binding functions of this protein.

An amino acid code for irregular and mixed protein packing.

To advance our understanding of protein tertiary structure, the development of the knob-socket model is completed in an analysis of the packing in irregular coil and turn secondary structure packing as well as between mixed secondary structure. The knob-socket model simplifies packing based on repeated patterns of two motifs: a three-residue socket for packing within secondary (2°) structure and a four-residue knob-socket for tertiary (3°) packing. For coil and turn secondary structure, knob-sockets allow identification of a correlation between amino acid composition and tertiary arrangements in space. Coil contributes almost as much as α-helices to tertiary packing. In irregular sockets, Gly, Pro, Asp, and Ser are favored, while in irregular knobs, the preference order is Arg, Asp, Pro, Asn, Thr, Leu, and Gly. Cys, His,Met, and Trp are not favored in either. In mixed packing, the knob amino acid preferences are a function of the socket that they are packing into, whereas the amino acid composition of the sockets does not depend on the secondary structure of the knob. A unique motif of a coil knob with an XYZ β-sheet socket may potentially function to inhibit β-sheet extension. In addition, analysis of the preferred crossing angles for strands within a β-sheet and mixed α-helice/β-sheet identifies canonical packing patterns useful in protein design. Lastly, the knob-socket model abstracts the complexity of protein tertiary structure into an intuitive packing surface topology map.

Morphology-Dependent HIV-Enhancing Effect of Semen-Derived Enhancer of Viral Infection

PAP248–286 is a 39-residue fragment (residues 248 to 286) derived from protease cleavage of prostatic acidic phosphatase in semen. The amyloid fibrils formed in vitro by PAP248–286 can dramatically enhance human immunodeficiency virus (HIV) infection. To our knowledge, we present the first report that the HIV-enhancing potency of fibrils formed by PAP248–286 is morphology dependent. We identified pleomorphic fibrils by transmission electron microscopy in two buffer conditions. Our solid-state NMR data showed that these fibrils consist of molecules in distinct conformations. In agreement with NMR, fluorescence measurements confirmed that they are assembled along different pathways, with distinct molecular structures. Furthermore, our cell-based infectivity tests detected distinct HIV-enhancing potencies for fibrils in distinct morphologies. In addition, our transmission electron microscopy and NMR results showed that semen-derived enhancer of viral infection fibrils formed in sodium bicarbonate buffer remain stable over time, but semen-derived enhancer of viral infection fibrils formed in phosphate buffered saline keep evolving after the initial 7 days incubation period. Given time, most of the assemblies in phosphate buffered saline will turn into elongated thin fibrils. They have similar secondary structure but different packing than thin fibrils formed initially after 7 days incubation.

Resonance Raman spectroscopic measurements delineate the structural changes that occur during tau fibril formation.

The aggregation of the microtubule-associated protein, tau, into amyloid fibrils is a hallmark of neurodegenerative diseases such as the tauopathies and Alzheimer's disease. Since monomeric tau is an intrinsically disordered protein and the polymeric fibrils possess an ordered cross-β core, the aggregation process is known to involve substantial conformational conversion besides growth. The aggregation mechanism of tau in the presence of inducers such as heparin, deciphered using probes such as thioflavin T/S fluorescence, light scattering, and electron microscopy assays, has been shown to conform to ligand-induced nucleation-dependent polymerization. These probes do not, however, distinguish between the processes of conformational conversion and fibril growth. In this study, UV resonance Raman spectroscopy is employed to look directly at signatures of changes in secondary structure and side-chain packing during fibril formation by the four repeat functional domain of tau in the presence of the inducer heparin, at pH 7 and at 37 °C. Changes in the positions and intensities of the amide Raman bands are shown to occur in two distinct stages during the fibril formation process. The first stage of UVRR spectral changes corresponds to the transformation of monomer into early fibrillar aggregates. The second stage corresponds to the transformation of these early fibrillar aggregates into the final, ordered, mature fibrils and during this stage; the processes of conformational conversion and the consolidation of the fibril core occur simultaneously. Delineation of these secondary structural changes accompanying the formation of tau fibrils holds significance for the understanding of generic and tau-specific principles of amyloid assembly.

Molecular determinants of expansivity of native globular proteins: a pressure perturbation calorimetry study.

There is a growing interest in understanding how hydrostatic pressure (P) impacts the thermodynamic stability (ΔG) of globular proteins. The pressure dependence of stability is defined by the change in volume upon denaturation, ΔV = (∂ΔG/∂P)T. The temperature dependence of change in volume upon denaturation itself is defined by the changes in thermal expansivity (ΔE), ΔE = (∂ΔV/∂T)P. The pressure perturbation calorimetry (PPC) allows direct experimental measurement of the thermal expansion coefficient, α = E/V, of a protein in the native, αN(T), and unfolded, αU(T), states as a function of temperature. We have shown previously that αU(T) is a nonlinear function of temperature but can be predicted well from the amino acid sequence using α(T) values for individual amino acids (J. Phys. Chem. B 2010, 114, 16166-16170). In this work, we report PPC results on a diverse set of nine proteins and discuss molecular factors that can potentially influence the thermal expansion coefficient, αN(T), and the thermal expansivity, EN(T), of proteins in the native state. Direct experimental measurements by PPC show that αN(T) and EN(T) functions vary significantly for different proteins. Using comparative analysis and site-directed mutagenesis, we have eliminated the role of various structural or thermodynamic properties of these proteins such as the number of amino acid residues, secondary structure content, packing density, electrostriction, dynamics, or thermostability. We have also shown that αN(T) and EN,sp(T) functions for a given protein are rather insensitive to the small changes in the amino acid sequence, suggesting that αN(T) and EN(T) functions might be defined by a topology of a given protein fold. This conclusion is supported by the similarity of αN(T) and EN(T) functions for six resurrected ancestral thioredoxins that vary in sequence but have very similar tertiary structure.

Exploring protein domain organization by recognition of secondary structure packing interfaces.

Protein domains are fundamental units of protein structure, function and evolution; thus, it is critical to gain a deep understanding of protein domain organization. Previous works have attempted to identify key residues involved in organization of domain architecture. Because one of the most important characteristics of domain architecture is the arrangement of secondary structure elements (SSEs), here we present a picture of domain organization through an integrated consideration of SSE arrangements and residue contact networks. In this work, by representing SSEs as main-chain scaffolds and side-chain interfaces and through construction of residue contact networks, we have identified the SSE interfaces well packed within protein domains as SSE packing clusters. In total, 17 334 SSE packing clusters were recognized from 9015 Structural Classification of Proteins domains of <40% sequence identity. The similar SSE packing clusters were observed not only among domains of the same folds, but also among domains of different folds, indicating their roles as common scaffolds for organization of protein domains. Further analysis of 14 small single-domain proteins reveals a high correlation between the SSE packing clusters and the folding nuclei. Consistent with their important roles in domain organization, SSE packing clusters were found to be more conserved than other regions within the same proteins. Supplementary data are available at Bioinformatics online.

An amino acid code for β-sheet packing structure

To understand the relationship between protein sequence and structure, this work extends the knob-socket model in an investigation of β-sheet packing. Over a comprehensive set of β-sheet folds, the contacts between residues were used to identify packing cliques: sets of residues that all contact each other. These packing cliques were then classified based on size and contact order. From this analysis, the two types of four-residue packing cliques necessary to describe β-sheet packing were characterized. Both occur between two adjacent hydrogen bonded β-strands. First, defining the secondary structure packing within β-sheets, the combined socket or XY:HG pocket consists of four residues i, i+2 on one strand and j, j+2 on the other. Second, characterizing the tertiary packing between β-sheets, the knob-socket XY:H+B consists of a three-residue XY:H socket (i, i+2 on one strand and j on the other) packed against a knob B residue (residue k distant in sequence). Depending on the packing depth of the knob B residue, two types of knob-sockets are found: side-chain and main-chain sockets. The amino acid composition of the pockets and knob-sockets reveal the sequence specificity of β-sheet packing. For β-sheet formation, the XY:HG pocket clearly shows sequence specificity of amino acids. For tertiary packing, the XY:H+B side-chain and main-chain sockets exhibit distinct amino acid preferences at each position. These relationships define an amino acid code for β-sheet structure and provide an intuitive topological mapping of β-sheet packing.

Coevolution signals capture the specific packing of secondary structures in protein architecture

Coevolution signals capture the specific packing of secondary structures in protein architecture

Conformational conversion, fibril growth and the ordering of the fibril core occur simultaneously during tau aggregation into amyloid fibrils

Amyloid fibril formation by the microtubule‐associated protein tau is a hallmark of neurodegenerative diseases such as Alzheimer's disease. The aggregation mechanism is known to involve substantial conformational conversion and growth since the intrinsically disordered monomer converts into ordered, β‐sheet rich fibrils. The focus of this work is to determine if these processes occur simultaneously or sequentially. The polymerization mechanism of tau in the presence of ligands, deciphered using probes such as thioflavin T fluorescence and light scattering, has been described in the past as ligand‐induced nucleation‐dependent polymerization. These probes cannot, however, distinguish the processes of conformational conversion and fibril growth. In this study, we employ UV Resonance Raman spectroscopy to look directly at signatures of changes in secondary structure and side chain packing, during fibril formation by the functional domain of tau in the presence of the inducer heparin, at pH 7 and at 37 °C. By examining changes in the positions and intensities of the amide and aromatic Raman bands, we find direct evidence for the conclusion that the processes of conformational conversion, fibril growth and the ordering of the fibril core occur simultaneously.This work was funded by the Tata Institute of Fundamental Research and the Department of Biotechnology, Government of India.

MICAN : a protein structure alignment algorithm that can handle Multiple-chains, Inverse alignments, C α only models, Alternative alignments, and Non-sequential alignments

BackgroundProtein pairs that have the same secondary structure packing arrangement but have different topologies have attracted much attention in terms of both evolution and physical chemistry of protein structures. Further investigation of such protein relationships would give us a hint as to how proteins can change their fold in the course of evolution, as well as a insight into physico-chemical properties of secondary structure packing. For this purpose, highly accurate sequence order independent structure comparison methods are needed.ResultsWe have developed a novel protein structure alignment algorithm, MICAN (a structure alignment algorithm that can handle Multiple-chain complexes, Inverse direction of secondary structures, Cα only models, Alternative alignments, and Non-sequential alignments). The algorithm was designed so as to identify the best structural alignment between protein pairs by disregarding the connectivity between secondary structure elements (SSE). One of the key feature of the algorithm is utilizing the multiple vector representation for each SSE, which enables us to correctly treat bent or twisted nature of long SSE. We compared MICAN with other 9 publicly available structure alignment programs, using both reference-dependent and reference-independent evaluation methods on a variety of benchmark test sets which include both sequential and non-sequential alignments. We show that MICAN outperforms the other existing methods for reproducing reference alignments of non-sequential test sets. Further, although MICAN does not specialize in sequential structure alignment, it showed the top level performance on the sequential test sets. We also show that MICAN program is the fastest non-sequential structure alignment program among all the programs we examined here.ConclusionsMICAN is the fastest and the most accurate program among non-sequential alignment programs we examined here. These results suggest that MICAN is a highly effective tool for automatically detecting non-trivial structural relationships of proteins, such as circular permutations and segment-swapping, many of which have been identified manually by human experts so far. The source code of MICAN is freely download-able at http://www.tbp.cse.nagoya-u.ac.jp/MICAN.

Folded functional lipid-poor apolipoprotein A-I obtained by heating of high-density lipoproteins: relevance to high-density lipoprotein biogenesis

HDL (high-density lipoproteins) remove cell cholesterol and protect from atherosclerosis. The major HDL protein is apoA-I (apolipoprotein A-I). Most plasma apoA-I circulates in lipoproteins, yet ~5% forms monomeric lipid-poor/free species. This metabolically active species is a primary cholesterol acceptor and is central to HDL biogenesis. Structural properties of lipid-poor apoA-I are unclear due to difficulties in isolating this transient species. We used thermal denaturation of human HDL to produce lipid-poor apoA-I. Analysis of the isolated lipid-poor fraction showed a protein/lipid weight ratio of 3:1, with apoA-I, PC (phosphatidylcholine) and CE (cholesterol ester) at approximate molar ratios of 1:8:1. Compared with lipid-free apoA-I, lipid-poor apoA-I showed slightly altered secondary structure and aromatic packing, reduced thermodynamic stability, lower self-associating propensity, increased adsorption to phospholipid surface and comparable ability to remodel phospholipids and form reconstituted HDL. Lipid-poor apoA-I can be formed by heating of either plasma or reconstituted HDL. We propose the first structural model of lipid-poor apoA-I which corroborates its distinct biophysical properties and postulates the lipid-induced ordering of the labile C-terminal region. In summary, HDL heating produces folded functional monomolecular lipid-poor apoA-I that is distinct from lipid-free apoA-I. Increased adsorption to phospholipid surface and reduced C-terminal disorder may help direct lipid-poor apoA-I towards HDL biogenesis.

FlexSnap: Flexible Non-sequential Protein Structure Alignment

BackgroundProteins have evolved subject to energetic selection pressure for stability and flexibility. Structural similarity between proteins that have gone through conformational changes can be captured effectively if flexibility is considered. Topologically unrelated proteins that preserve secondary structure packing interactions can be detected if both flexibility and Sequential permutations are considered. We propose the FlexSnap algorithm for flexible non-topological protein structural alignment.ResultsThe effectiveness of FlexSnap is demonstrated by measuring the agreement of its alignments with manually curated non-sequential structural alignments. FlexSnap showed competitive results against state-of-the-art algorithms, like DALI, SARF2, MultiProt, FlexProt, and FATCAT. Moreover on the DynDom dataset, FlexSnap reported longer alignments with smaller rmsd.ConclusionsWe have introduced FlexSnap, a greedy chaining algorithm that reports both sequential and non-sequential alignments and allows twists (hinges). We assessed the quality of the FlexSnap alignments by measuring its agreements with manually curated non-sequential alignments. On the FlexProt dataset, FlexSnap was competitive to state-of-the-art flexible alignment methods. Moreover, we demonstrated the benefits of introducing hinges by showing significant improvements in the alignments reported by FlexSnap for the structure pairs for which rigid alignment methods reported alignments with either low coverage or large rmsd.AvailabilityAn implementation of the FlexSnap algorithm will be made available online at http://www.cs.rpi.edu/~zaki/software/flexsnap.

Kinetically trapped metastable intermediate of a disulfide‐deficient mutant of the starch‐binding domain of glucoamylase

Refolding of a thermally unfolded disulfide-deficient mutant of the starch-binding domain of glucoamylase was investigated using differential scanning calorimetry, isothermal titration calorimetry, CD, and (1)H NMR. When the protein solution was rapidly cooled from a higher temperature, a kinetic intermediate was formed during refolding. The intermediate was unexpectedly stable compared with typical folding intermediates that have short half-lives. It was shown that this intermediate contained substantial secondary structure and tertiary packing and had the same binding ability with beta-cyclodextrin as the native state, suggesting that the intermediate is highly-ordered and native-like on the whole. These characteristics differ from those of partially folded intermediates such as molten globule states. Far-UV CD spectra showed that the secondary structure was once disrupted during the transition from the intermediate to the native state. These results suggest that the intermediate could be an off-pathway type, possibly a misfolded state, that has to undergo unfolding on its way to the native state.

Protein structure prediction aided by geometrical and probabilistic constraints

Database-assisted ab initio protein structure prediction methods have exhibited considerable promise in the recent past, with several implementations being successful in community-wide experiments (CASP). We have employed combinatorial optimization techniques toward solving the protein structure prediction problem. A Monte Carlo minimization algorithm has been employed on a constrained search space to identify minimum energy configurations. The search space is constrained by using radius of gyration cutoffs, the loop backbone dihedral probability distributions, and various secondary structure packing conformations. Simulations have been carried out on several sequences and 1000 conformations have been initially generated. Of these, 50 best candidates have then been selected as probable conformations. The search for the optimum has been simplified by incorporating various geometrical constraints on secondary structural elements using distance restraint potential functions. The advantages of the reported methodology are its simplicity, and modifiability to include other geometric and probabilistic restraints.