From isotropic to anisotropic side chain representations: comparison of three models for residue contact estimation

Long-range Regulation of Partially Folded Amyloidogenic Peptides

Neurodegeneration involves abnormal aggregation of intrinsically disordered amyloidogenic peptides (IDPs), usually mediated by hydrophobic protein-protein interactions. There is mounting evidence that formation of α-helical intermediates is an early event during self-assembly of amyloid-β42 (Aβ42) and α-synuclein (αS) IDPs in Alzheimer’s and Parkinson’s disease pathogenesis, respectively. However, the driving force behind on-pathway molecular assembly of partially folded helical monomers into helical oligomers assembly remains unknown. Here, we employ extensive molecular dynamics simulations to sample the helical conformational sub-spaces of monomeric peptides of both Aβ42 and αS. Our computed free energies, population shifts, and dynamic cross-correlation network analyses reveal a common feature of long-range intra-peptide modulation of partial helical folds of the amyloidogenic central hydrophobic domains via concerted coupling with their charged terminal tails (N-terminus of Aβ42 and C-terminus of αS). The absence of such inter-domain fluctuations in both fully helical and completely unfolded (disordered) states suggests that long-range coupling regulates the dynamicity of partially folded helices, in both Aβ42 and αS peptides. The inter-domain coupling suggests a form of intra-molecular allosteric regulation of the aggregation trigger in partially folded helical monomers. This approach could be applied to study the broad range of amyloidogenic peptides, which could provide a new path to curbing pathogenic aggregation of partially folded conformers into oligomers, by inhibition of sites far from the hydrophobic core.

Characteristics Analyses and Comparisons of the Protein Structure Networks Constructed by Different Methods

Protein structure networks (PSNs) were widely used in analyses of protein structure and function. In this work, we analyzed and compared the characters of PSNs by different methods. The degrees of the different types of the nodes were found to be associated with the amino acid characters, including SAS, secondary structure, hydropathy and the volume of amino acids. It showed that PSNs by the methods of CA10, SC10 and AT5 inherited more amino acid characters and had higher correlations with the original protein structures. And PSNs by these three methods would be powerful tools in understanding the characters of protein structures.

Residue interaction network analysis of Dronpa and a DNA clamp

Topology is an essential aspect of protein structure. The network paradigm is increasingly used to describe the topology and dynamics of proteins. In this paper, the effect of topology on residue interaction network was investigated for two different proteins: Dronpa and a DNA clamp, which have cylindrical and toroidal topologies, respectively. Network metrics including characteristic path lengths, clustering coefficients, and diameters were calculated to investigate their global topology parameters such as small-world properties and packing density. Measures of centrality including betweenness, closeness, and residue centrality were computed to predict residues critical to function. Additionally, the detailed topology of the hydrophobic pocket in Dronpa, and communication pathways across the interface in the DNA clamp, were investigated using the network. The results are presented and discussed with regard to existing residue interaction network properties of globular proteins and elastic network models on Dronpa and the DNA clamp. The topological principle underlying residue interaction networks provided insight into the architectural organization of proteins.

Design of a rotamer library for coarse-grained models in protein-folding simulations

Rotamer libraries usually contain geometric information to trace an amino acid side chain, atom by atom, onto a protein backbone. These libraries have been widely used in protein design, structure refinement and prediction, homology modeling, and X-ray and NMR structure validation. However, they usually present too much information and are not always fully compatible with the coarse-grained models of the protein geometry that are frequently used to tackle the protein-folding problem through molecular simulation. In this work, we introduce a new backbone-dependent rotamer library for side chains compatible with low-resolution models in polypeptide chains. We have dispensed with an atomic description of proteins, representing each amino acid side chain by its geometric center (or centroid). The resulting rotamers have been estimated from a statistical analysis of a large structural database consisting of high-resolution X-ray protein structures. As additional information, each rotamer includes the frequency with which it has been found during the statistical analysis. More importantly, the library has been designed with a careful control to ensure that the vast majority of side chains in protein structures (at least 95% of residues) are properly represented. We have tested our library using an independent set of proteins, and our results support a good correlation between the reconstructed centroids from our rotamer library and those in the experimental structures. This new library can serve to improve the definition of side chain centroids in coarse-grained models, avoiding at the same time an excessive additional complexity in a geometric model for the polypeptide chain.

A beta-complex statistical four body contact potential combined with a hydrogen bond statistical potential recognizes the correct native structure from protein decoy sets

We present a new four-body knowledge-based potential for recognizing the native state of proteins from their misfolded states. This potential was extracted from a large set of protein structures determined by X-ray crystallography using BetaMol, a software based on the recent theory of the beta-complex (β-complex) and quasi-triangulation of the Voronoi diagram of spheres. This geometric construct reflects the size difference among atoms in their full Euclidean metric; property not accounted for in a typical 3D Delaunay triangulation. The ability of this potential to identify the native conformation over a large set of decoys was evaluated. Experiments show that this potential outperforms a potential constructed with a classical Delaunay triangulation in decoy discrimination tests. The addition of a statistical hydrogen bond potential to our four-body potential allows a significant improvement in the decoy discrimination, in such a way that we are able to predict successfully the native structure in 90% of cases.

Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review

Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.

The effect of edge definition of complex networks on protein structure identification

The main objective of this study is to explore the contribution of complex network together with its different definitions of vertexes and edges to describe the structure of proteins. Protein folds into a specific conformation for its function depending on interactions between residues. Consequently, in many studies, a protein structure was treated as a complex system comprised of individual components residues, and edges were interactions between residues. What is the proper time for representing a protein structure as a network? To confirm the effect of different definitions of vertexes and edges in constructing the amino acid interaction networks, protein domains and the structural unit of proteins were described using this method. The identification performance of 2847 proteins with domain/domains proved that the structure of proteins was described well when R_Cα was around 5.0–7.5 Å, and the optimal cutoff value for constructing the protein structure networks was 5.0 Å (C_α-C_α distances) while the ideal community division method was community structure detection based on edge betweenness in this study.

A modified amino acid network model contains similar and dissimilar weight

For a more detailed description of the interaction between residues, this paper proposes an amino acid network model, which contains two types of weight—similar weight and dissimilar weight. The weight of the link is based on a self-consistent statistical contact potential between different types of amino acids. In this model, we can get a more reasonable representation of the distance between residues. Furthermore, with the network parameter, average shortest path length, we can get a more accurate reflection of the molecular size. This amino acid network is a “small-world” network, and the network parameter is sensitive to the conformation change of protein. For some disease-related proteins, the highly central residues of the amino acid network are highly correlated with the hot spots. In the compound with the related drug, these residues either interacted directly with the drug or with the residue which is in contact with the drug.

Three-dimensional structures of membrane proteins from genomic sequencing

We show that amino acid covariation in proteins, extracted from the evolutionary sequence record, can be used to fold transmembrane proteins. We use this technique to predict previously unknown 3D structures for 11 transmembrane proteins (with up to 14 helices) from their sequences alone. The prediction method (EVfold_membrane) applies a maximum entropy approach to infer evolutionary covariation in pairs of sequence positions within a protein family and then generates all-atom models with the derived pairwise distance constraints. We benchmark the approach with blinded de novo computation of known transmembrane protein structures from 23 families, demonstrating unprecedented accuracy of the method for large transmembrane proteins. We show how the method can predict oligomerization, functional sites, and conformational changes in transmembrane proteins. With the rapid rise in large-scale sequencing, more accurate and more comprehensive information on evolutionary constraints can be decoded from genetic variation, greatly expanding the repertoire of transmembrane proteins amenable to modeling by this method.