Sequence and structural analysis of two designed proteins with 88% identity adopting different folds

Hexadecanoic acid analogs as potential CviR-mediated quorum sensing inhibitors in Chromobacterium violaceum: an in silico study

Chromobacterium violaceum is a Gram-negative, rod-shaped and opportunistic human pathogen. C. violaceum is resistant to various antibiotics due to the production of quorum sensing (QS)-controlled virulence factor and biofilm formation. Hence, we need to find alternative strategies to overcome the antimicrobial resistance and biofilm formation in Gram-negative bacteria. QS is a mechanism in which bacteria's ability to regulate the virulence factors and biofilm formations leads to disease progression. Previously, hexadecanoic acid was identified as a CviR-mediated quorum-sensing inhibitor. In this study, we aimed to discover potential analogs of hexadecanoic acid as a CviR-mediated quorum-sensing inhibitor against C. violaceum by using ADME/T prediction, density functional theory, molecular docking, molecular dynamics and free energy binding calculations. ADME/T properties predicted for analogs were acceptable for human oral absorption and feasibility. The highest occupied molecular orbitals and lowest unoccupied molecular orbitals gap energies predicted and found oleic acid with -0.3748 energies. Docosatrienoic acid exhibited the highest binding affinity -8.15 Kcal/mol and strong and stable interactions with the amino acid residues on the active site of the CviR protein. These compounds on MD simulations for 100 ns show strong hydrogen-bonding interactions with the protein and remain stable inside the active site. Our results suggest hexadecanoic acid analogs could serve as anti-QS and anti-biofilm molecules for treating C. violaceum infections. However, further validation and investigation of these inhibitors against CviR are needed to claim their candidacy for clinical trials.Communicated by Ramaswamy H. Sarma.

Multiobjective heuristic algorithm for de novo protein design in a quantified continuous sequence space

Protein design usually involves sequence search process and evaluation criteria. Commonly used methods primarily implement the Monte Carlo or simulated annealing algorithm with a single-energy function to obtain ideal solutions, which is often highly time-consuming and limited by the accuracy of the energy function. In this report, we introduce a multiobjective algorithm named Hydra for protein design, which employs two different energy functions to optimize solutions simultaneously and makes use of the latent quantitative relationship between different amino acid types to facilitate the search process. The framework uses two kinds of prior information to transform the original disordered discrete sequence space into a relatively ordered space, and decoy sequences are searched in this ordered space through a multiobjective swarm intelligence algorithm. This algorithm features high accuracy and a high-speed search process. Our method was tested on 40 targets covering different fold classes, which were computationally verified to be well folded, and it experimentally solved the 1UBQ fold by NMR in excellent agreement with the native structure with a backbone RMSD deviation of 1.074 Å. The Hydra software package can be downloaded from: http://www.csbio.sjtu.edu.cn/bioinf/HYDRA/ for academic use.

Sequence and structural analysis of fibronectin-binding protein reveals importance of multiple intrinsic disordered tandem repeats

The location of certain amino acid sequences like repeats along the polypeptide chain is very important in the context of forming the overall shape of the protein molecule which in fact determines its function. In gram-positive bacteria, fibronectin-binding protein (FnBP) is one such repeat containing protein, and it is a cell wall-attached protein responsible for various acute infections in human. Several studies on sequence, structure, and function of fibronectin-binding regions of FnBPs were reported; however, no detailed study was carried out on the full-length protein sequence. In the present study, we have made a thorough sequence and structure analysis on FnBP_A of Staphylococcus aureus and explored the presence of dual ligand-binding ability of fibrinogen (fg)-binding region and its molecular recognition processes. Multiple sequence alignment and protein-protein docking analysis reveal the regions which are likely involved in dual ligand binding. Further analysis of docking of FnBP_A fg-binding region and fn N-terminal modules suggests that if the latter binds to the fg-binding region of FnBP_A, it would inhibit the subsequent binding of fg because of steric hindrance. The sequence analysis further suggests that the abundance of disorder promoting residue glutamic acid and dual personality (both order/disorder promoting) residue threonine in tandem repeats of FnBP_A and B proteins possibly would help the molecule to undergo a conformational change while binding with fn by β-zipper mechanism. The segment-based power spectral analysis was carried out which helps to understand the distribution of hydrophobic residues along the sequence particularly in intrinsic disordered tandem repeats. The results presented here will help to understand the role of internal repeats and intrinsic disorder in the molecular recognition process of a pathogenic cell surface protein.

Dihedral angle preferences of amino acid residues forming various non-local interactions in proteins

In theory, a polypeptide chain can adopt a vast number of conformations, each corresponding to a set of backbone rotation angles. Many of these conformations are excluded due to steric overlaps. Ramachandran and coworkers were the first to look into this problem by plotting backbone dihedral angles in a two-dimensional plot. The conformational space in the Ramachandran map is further refined by considering the energetic contributions of various non-bonded interactions. Alternatively, the conformation adopted by a polypeptide chain may also be examined by investigating interactions between the residues. Since the Ramachandran map essentially focuses on local interactions (residues closer in sequence), out of interest, we have analyzed the dihedral angle preferences of residues that make non-local interactions (residues far away in sequence and closer in space) in the folded structures of proteins. The non-local interactions have been grouped into different types such as hydrogen bond, van der Waals interactions between hydrophobic groups, ion pairs (salt bridges), and ππ-stacking interactions. The results show the propensity of amino acid residues in proteins forming local and non-local interactions. Our results point to the vital role of different types of non-local interactions and their effect on dihedral angles in forming secondary and tertiary structural elements to adopt their native fold.

Pairwise contact energy statistical potentials can help to find probability of point mutations

To adopt a particular fold, a protein requires several interactions between its amino acid residues. The energetic contribution of these residue-residue interactions can be approximated by extracting statistical potentials from known high resolution structures. Several methods based on statistical potentials extracted from unrelated proteins are found to make a better prediction of probability of point mutations. We postulate that the statistical potentials extracted from known structures of similar folds with varying sequence identity can be a powerful tool to examine probability of point mutation. By keeping this in mind, we have derived pairwise residue and atomic contact energy potentials for the different functional families that adopt the (α/β)₈ TIM-Barrel fold. We carried out computational point mutations at various conserved residue positions in yeast Triose phosphate isomerase enzyme for which experimental results are already reported. We have also performed molecular dynamics simulations on a subset of point mutants to make a comparative study. The difference in pairwise residue and atomic contact energy of wildtype and various point mutations reveals probability of mutations at a particular position. Interestingly, we found that our computational prediction agrees with the experimental studies of Silverman et al. (Proc Natl Acad Sci 2001;98:3092-3097) and perform better prediction than i_Mutant and Cologne University Protein Stability Analysis Tool. The present work thus suggests deriving pairwise contact energy potentials and molecular dynamics simulations of functionally important folds could help us to predict probability of point mutations which may ultimately reduce the time and cost of mutation experiments. Proteins 2016; 85:54-64. © 2016 Wiley Periodicals, Inc.

Molecular dynamics simulations on the conformational transitions from the GA 98 (GA 88) to GB 98 (GB 88) proteins

We performed conventional and targeted molecular dynamics simulations to address the dynamic transition mechanisms of the conformational transitions from the G_A 98 protein with only 1 mutation of Leu45Tyr to G_B 98 and from the G_A 88 protein with 7 mutations of Gly24Ala, Ile25Thr, Ile30Phe, Ile33Tyr, Leu45Tyr, Ile49Thr, and Leu50Lys to G_B 88. The results show that the conformational transition mechanism from the mutated 3α G_A 98 (G_A 88) state to the α+4β G_B 98 (G_B 88) state via several intermediate conformations involves the bending of loops at the N and C termini firstly, the unfolding of αA and αC, then the traversing of αB, and the formation of the 4β layer with the conversion of the hydrophobic core. The bending of loops at the N and C termini and the formation of the crucial transition conformation with the full unfolded structure are key factors in their transition processes. The communication of the interaction network, the bending directions of loops, and the traversing site of αB in the transition of G_A 98 to G_B 98 are markedly different from those in G_A 88 to G_B 88 because of the different mutated residues. The analysis of the correlations and the calculated mass center distances between some segments further supported their conformational transition mechanisms. These results could help people to better understand the Paracelsus challenge. Copyright © 2016 John Wiley & Sons, Ltd.

PreFRP: Prediction and visualization of fluctuation residues in proteins

Aim:

The PreFRP web server extracts sequence and basic information of a protein structure and groups amino acid residues in a protein into three important types such as high, moderate, and weak fluctuating residues.

Materials and Methods:

The server takes a protein data bank file or an amino acid sequence as input and prints the probability of amino acid residues to fluctuate. The server also provides a link to Jmol, a molecular visualization program to visualize the high, moderate, and weak fluctuating residues in three different colors.

Results:

Prediction and visualization of fluctuating amino acid residues in proteins may help to understand the complex three-dimensional structure of proteins and may further help in docking and mutation experiments.

Availability:

The web server is freely accessible through the web page of the author's institution http://www.mpi.edu.in/prefrp/link.html.

Detailed computational analysis revealed mutation V210I on PrP induced conformational conversion on β2–α2 loop and α2–α3

The development of fatal transmissible spongiform encephalopathies (TSE) is associated with the conformational conversion of the normal cellular prion protein, PrP^C, into its pathogenic isoform, PrP^Sc.

Analysis of secondary structural and physicochemical changes in protein-protein complexes

Conformation switching in protein-protein complexes is considered important for the molecular recognition process. Overall analysis of 123 protein-protein complexes in a benchmark data-set showed that 6.8% of residues switched over their secondary structure conformation upon complex formation. Amino acid residue-wise preference for conformation change has been analyzed in binding and non-binding site residues separately. In this analysis, residues such as Ser, Leu, Glu, and Lys had higher frequency of secondary structural conformation change. The change of helix to coil and sheet to coil conformation and vice versa has been observed frequently, whereas the conformation change of helix to extended sheet occurred rarely in the studied complexes. Influence of conformation change toward the N and C terminal on either side of the binding site residues has been analyzed. Further, analysis on φ and ψ angle variation, conservation, stability, and solvent accessibility have been performed on binding site residues. Knowledge obtained from the present study could be effectively employed in the protein-protein modeling and docking studies.

Better theoretical models and protein design experiments can help to understand protein folding

In our study, we have concluded that two proteins with 88% homology choose different energetically favorable pathways in the very early stage of the folding process to attain their native folds. Subsequent reports from other investigators by performing folding and unfolding kinetics experiments concur with our findings. We herewith discuss the key papers revealing computational and experimental analysis of two designed proteins with similar sequence distant folds. Further we suggest that the theoretical/computational analysis of protein sequences and structures along with the relevant experiments provide a better understanding of the relationship between protein sequence, folding, and structure.

An evolution-based approach to De Novo protein design and case study on Mycobacterium tuberculosis

Computational protein design is a reverse procedure of protein folding and structure prediction, where constructing structures from evolutionarily related proteins has been demonstrated to be the most reliable method for protein 3-dimensional structure prediction. Following this spirit, we developed a novel method to design new protein sequences based on evolutionarily related protein families. For a given target structure, a set of proteins having similar fold are identified from the PDB library by structural alignments. A structural profile is then constructed from the protein templates and used to guide the conformational search of amino acid sequence space, where physicochemical packing is accommodated by single-sequence based solvation, torsion angle, and secondary structure predictions. The method was tested on a computational folding experiment based on a large set of 87 protein structures covering different fold classes, which showed that the evolution-based design significantly enhances the foldability and biological functionality of the designed sequences compared to the traditional physics-based force field methods. Without using homologous proteins, the designed sequences can be folded with an average root-mean-square-deviation of 2.1 Å to the target. As a case study, the method is extended to redesign all 243 structurally resolved proteins in the pathogenic bacteria Mycobacterium tuberculosis, which is the second leading cause of death from infectious disease. On a smaller scale, five sequences were randomly selected from the design pool and subjected to experimental validation. The results showed that all the designed proteins are soluble with distinct secondary structure and three have well ordered tertiary structure, as demonstrated by circular dichroism and NMR spectroscopy. Together, these results demonstrate a new avenue in computational protein design that uses knowledge of evolutionary conservation from protein structural families to engineer new protein molecules of improved fold stability and biological functionality.

The goal of computational protein design is to create new protein sequences of desirable structure and biological function. Most protein design methods are developed to search for sequences with the lowest free-energy based on physics-based force fields following Anfinsen's thermodynamic hypothesis. A major obstacle of such approaches is the inaccuracy of the force-field design, which cannot accurately describe atomic interactions or correctly recognize protein folds. We propose a novel method which uses evolutionary information, in the form of sequence profiles from structure families, to guide the sequence design. Since sequence profiles are generally more accurate than physics-based potentials in protein fold recognition, a unique advantage lies on that it targets the design procedure to a family of protein sequence profiles to enhance the robustness of designed sequences. The method was tested on 87 proteins and the designed sequences can be folded by I-TASSER to models with an average RMSD 2.1 Å. As a case study of large-scale application, the method is extended to redesign all structurally resolved proteins in the human pathogenic bacteria, Mycobacterium tuberculosis. Five sequences varying in fold and sizes were characterized by circular dichroism and NMR spectroscopy experiments and three were shown to have ordered tertiary structure.

Search and analysis of identical reverse octapeptides in unrelated proteins

For the past few decades, intensive studies have been carried out in an attempt to understand how the amino acid sequences of proteins encode their three dimensional structures to perform their specific functions. In order to understand the sequence-structure relationship of proteins, several sub-sequence search studies in non-redundant sequence-structure databases have been undertaken which have given some fruitful clues. In our earlier work, we analyzed a set of 3124 non-redundant protein sequences from the Protein Data Bank (PDB) and retrieved 30 identical octapeptides having different secondary structures. These octapeptides were characterized by using different computational procedures. This prompted us to explore the presence of octapeptides with reverse sequences and to analyze whether these octapeptides would adopt similar structures as that of their parent octapeptides. Our identical reverse octapeptide search resulted in the finding of eight octapeptide pairs (octapeptide and reverse octapeptide) with similar secondary structure and 23 octapeptide pairs with different secondary structures. In the present work, the geometrical and biophysical characteristics of identical reverse octapeptides were explored and compared with unrelated octapeptide pairs by using various computational tools. We thus conclude that proteins containing identical reverse octapeptides are not very abundant and residues in the octapeptide pairs do not contribute to the stability of the protein. Furthermore, compared to unrelated octapeptides, identical reverse octapeptides do not show certain biophysical and geometrical properties.

Search for identical octapeptides in unrelated proteins: Structural plasticity revisited

Since proteins are dynamic in nature, they can alter their local structure in response to changes in their environment factors such as temperature, pH, phosphorylation, and binding of other small molecules. These conformational changes are extremely important for the correct folding and functioning of proteins. There are also a number of diseases associated with protein conformational change such as amyloid diseases. To stimulate research into the above factors which specify one conformation over another, different theoretical models have been proposed and tested against sequence similar distant structure protein fragments. In order to simplify the computational complexity of identifying conformational changes in proteins, various local sequence search algorithms were employed and the structural plasticity in unrelated proteins was examined by various research groups. In the present work, we revisit the mechanism of structural plasticity in unrelated proteins with increased number of structures in Protein Data Bank by comparing identical octapeptides in unrelated proteins with dictionary of protein secondary structure extracted from existing experimental data. Our goal is to bring out the influence of hydrophobic residues, hydrophilic residues, flanking residues, difference in secondary structural propensities of surrounding residues, difference in phi-psi angles and local and nonlocal interactions in identical octapeptides adopting different conformations. Also we have used surrounding hydrophobicity, environment dependent interaction energy, atomic mean force potential, structural unit contacts and difference profiles models to explore the factors which cause structural plasticity. The results discussed here may provide insights into protein folding, design and function.