Accurate prediction of protein folding rates from sequence and sequence‐derived residue flexibility and solvent accessibility

Protein folding rates vary by several orders of magnitude and they depend on the topology of the fold and the size and composition of the sequence. Although recent works show that the rates can be predicted from the sequence, allowing for high‐throughput annotations, they consider only the sequence and its predicted secondary structure. We propose a novel sequence‐based predictor, PFR‐AF, which utilizes solvent accessibility and residue flexibility predicted from the sequence, to improve predictions and provide insights into the folding process. The predictor includes three linear regressions for proteins with two‐state, multistate, and unknown (mixed‐state) folding kinetics. PFR‐AF on average outperforms current methods when tested on three datasets. The proposed approach provides high‐quality predictions in the absence of similarity between the predicted and the training sequences. The PFR‐AF's predictions are characterized by high (between 0.71 and 0.95, depending on the dataset) correlation and the lowest (between 0.75 and 0.9) mean absolute errors with respect to the experimental rates, as measured using out‐of‐sample tests. Our models reveal that for the two‐state chains inclusion of solvent‐exposed Ala may accelerate the folding, while increased content of Ile may reduce the folding speed. We also demonstrate that increased flexibility of coils facilitates faster folding and that proteins with larger content of solvent‐exposed strands may fold at a slower pace. The increased flexibility of the solvent‐exposed residues is shown to elongate folding, which also holds, with a lower correlation, for buried residues. Two case studies are included to support our findings. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Sources and significance of variation in basal,summit and maximal metabolic rates in birds

The rates at which birds use energy may have profound effects on fitness, thereby influencing physiology, behavior, ecology and evolution. Comparisons of standardized metabolic rates (e.g., lower and upper limits of metabolic power output) present a method for elucidating the effects of ecological and evolutionary factors on the interface between physiology and life history in birds. In this paper we review variation in avian metabolic rates [basal metabolic rate (BMR; minimum normothermic metabolic rate), ...     

Resource distribution influences mating system in the bobuck (<Emphasis Type="Italic">Trichosurus cunninghami</Emphasis>: Marsupialia)

Mammalian mating systems are thought to be shaped by the spatial distribution and abundance of key resources, which in turn influence the spacing behaviour of individuals. In particular, female home range size is predicted to reflect the availability of key resources. We documented the availability and distribution of food and shelter resources for two neighbouring populations of bobucks, or mountain brushtail possums, Trichosurus cunninghami, that were characterised by different mating systems: our “forest population” was socially monogamous, whereas the “roadside population” was polygynous. Both silver wattle, Acacia dealbata, the main food resource for bobucks, and den-trees, which provided shelter, occurred at significantly higher density at the roadside site. The pattern of distribution of these two resources also differed between the sites. Both food and den-trees were scattered evenly throughout the roadside habitat. In contrast, den-trees were located predominantly at one end of the forest site, while silver wattle trees were located at the other. There was no significant difference in the amount of silver wattle, or in the number of den-trees, located within the home ranges of individual females at the two sites. However, forest females had home ranges, on average, almost three times the size of those of roadside females. At the roadside site, the size of female home ranges varied inversely with the density of silver wattle, indicating that these females ranged over as large an area as necessary to gain access to sufficient silver wattle trees. There was no such relationship among forest females. These populations provide a clear example of resource distribution determining female home range size. This influenced the number of female home ranges a male’s home range overlapped with, which in turn determined the mating system. Such clear links between resource availability and mating system have not previously been established in a marsupial.     

Anti-predator defence and the complexity-stability relationship of food webs

The mechanism for maintaining complex food webs has been a central issue in ecology because theory often predicts that complexity (higher the species richness, more the interactions) destabilizes food webs. Although it has been proposed that prey anti-predator defence may affect the stability of prey-predator dynamics, such studies assumed a limited and relatively simpler variation in the food-web structure. Here, using mathematical models, I report that food-web flexibility arising from prey anti-predator defence enhances community-level stability (community persistence and robustness) in more complex systems and even changes the complexity-stability relationship. The model analysis shows that adaptive predator-specific defence enhances community-level stability under a wide range of food-web complexity levels and topologies, while generalized defence does not. Furthermore, while increasing food-web complexity has minor or negative effects on community-level stability in the absence of defence adaptation, or in the presence of generalized defence, in the presence of predator-specific defence, the connectance-stability relationship may become unimodal. Increasing species richness, in contrast, always lowers community-level stability. The emergence of a positive connectance-stability relationship however necessitates food-web compartmentalization, high defence efficiency and low defence cost, suggesting that it only occurs under a restricted condition.     

Genetic population structure, queen supersedure and social polymorphism in a social Hymenoptera

In social insects, the emergence of multiple queening is linked to changes in a suite of traits such as the reproductive life span of queens, mating patterns and population structure. We investigated queen turnover, colony longevity, spatial distribution patterns and genetic differentiation in a population of the socially polymorphic ant Formica fusca. Genetic differentiation between the social forms was absent, and mating patterns were similar in the two forms. The spatial distribution of single- and multi-queen colonies indicated an absence of colony reproduction by budding in both colony types. However, the rate of queen supersedure was high in multi-queen colonies and absent in single-queen ones. The social structure of colonies remained stable across years, but colony mortality did not differ between the two social forms. These results imply that differences between social types may appear and persist also in sympatry, and that these differences may occur in some traits, but not others, despite the presence of homogenizing gene flow.     

Local motions in a benchmark of allosteric proteins

Allosteric proteins have been studied extensively in the last 40 years, but so far, no systematic analysis of conformational changes between allosteric structures has been carried out. Here, we compile a set of 51 pairs of known inactive and active allosteric protein structures from the Protein Data Bank. We calculate local conformational differences between the two structures of each protein using simple metrics, such as backbone and side-chain Cartesian displacement, and torsion angle change and rearrangement in residue-residue contacts. Thresholds for each metric arise from distributions of motions in two control sets of pairs of protein structures in the same biochemical state. Statistical analysis of motions in allosteric proteins quantifies the magnitude of allosteric effects and reveals simple structural principles about allostery. For example, allosteric proteins exhibit substantial conformational changes comprising about 20% of the residues. In addition, motions in allosteric proteins show strong bias toward weakly constrained regions such as loops and the protein surface. Correlation functions show that motions communicate through protein structures over distances averaging 10-20 residues in sequence space and 10-20 A in Cartesian space. Comparison of motions in the allosteric set and a set of 21 nonallosteric ligand-binding proteins shows that nonallosteric proteins also exhibit bias of motion toward weakly constrained regions and local correlation of motion. However, allosteric proteins exhibit twice as much percent motion on average as nonallosteric proteins with ligand-induced motion. These observations may guide efforts to design flexibility and allostery into proteins.     

Solution Structural Studies of Pre-amyloid Oligomer States of the Biofilm Protein Aap

Staphylococcus epidermidis is a commensal bacterium on human skin that is also the leading cause of medical device-related infections. The accumulation-associated protein (Aap) from S. epidermidis is a critical factor for infection via its ability to mediate biofilm formation. The B-repeat superdomain of Aap is composed of 5 to 17 Zn²⁺-binding B-repeats, which undergo rapid, reversible assembly to form dimer and tetramer species. The tetramer can then undergo a conformational change and nucleate highly stable functional amyloid fibrils. In this study, multiple techniques including analytical ultracentrifugation (AUC) and small-angle X-ray scattering (SAXS) are used to probe a panel of B-repeat mutant constructs that assemble to distinct oligomeric states to define the structural characteristics of B-repeat dimer and tetramer species. The B-repeat region from Aap forms an extremely elongated conformation that presents several challenges for standard SAXS analyses. Specialized approaches, such as cross-sectional analyses, allowed for in-depth interpretation of data, while explicit-solvent calculations via WAXSiS allowed for accurate evaluation of atomistic models. The resulting models suggest mechanisms by which Aap functional amyloid fibrils form, illuminating an important contributing factor to recurrent staphylococcal infections.     

Conformational dynamics of nonsynonymous variants at protein interfaces reveals disease association

Recent studies have shown that the protein interface sites between individual monomeric units in biological assemblies are enriched in disease‐associated non‐synonymous single nucleotide variants (nsSNVs). To elucidate the mechanistic underpinning of this observation, we investigated the conformational dynamic properties of protein interface sites through a site‐specific structural dynamic flexibility metric (dfi) for 333 multimeric protein assemblies. dfi measures the dynamic resilience of a single residue to perturbations that occurred in the rest of the protein structure and identifies sites contributing the most to functionally critical dynamics. Analysis of dfi profiles of over a thousand positions harboring variation revealed that amino acid residues at interfaces have lower average dfi (31%) than those present at non‐interfaces (50%), which means that protein interfaces have less dynamic flexibility. Interestingly, interface sites with disease‐associated nsSNVs have significantly lower average dfi (23%) as compared to those of neutral nsSNVs (42%), which directly relates structural dynamics to functional importance. We found that less conserved interface positions show much lower dfi for disease nsSNVs as compared to neutral nsSNVs. In this case, dfi is better as compared to the accessible surface area metric, which is based on the static protein structure. Overall, our proteome‐wide conformational dynamic analysis indicates that certain interface sites play a critical role in functionally related dynamics (i.e., those with low dfi values), therefore mutations at those sites are more likely to be associated with disease. Proteins 2015; 83:428–435. © 2014 Wiley Periodicals, Inc.     

Microsecond simulations of mdm2 and its complex with p53 yield insight into force field accuracy and conformational dynamics

The p53‐MDM2 complex is both a major target for cancer drug development and a valuable model system for computational predictions of protein‐ligand binding. To investigate the accuracy of molecular simulations of MDM2 and its complex with p53, we performed a number of long (200 ns to 1 µs) explicit‐solvent simulations using a range of force fields. We systematically compared nine popular force fields (AMBER ff03, ff12sb, ff14sb, ff99sb, ff99sb‐ildn, ff99sb‐ildn‐nmr, ff99sb‐ildn‐phi, CHARMM22*, and CHARMM36) against experimental chemical shift data, and found similarly accurate results, with microsecond simulations achieving better agreement compared to 200‐ns trajectories. Although the experimentally determined apo structure has a closed binding cleft, simulations in all force fields suggest the apo state of MDM2 is highly flexible, and able to sample holo‐like conformations, consistent with a conformational selection model. Initial structuring of the MDM2 lid region, known to competitively bind the binding cleft, is also observed in long simulations. Taken together, these results show molecular simulations can accurately sample conformations relevant for ligand binding. We expect this study to inform future computational work on folding and binding of MDM2 ligands. Proteins 2015; 83:1665–1676. © 2015 Wiley Periodicals, Inc.     

Vibrational entropy differences between mesophile and thermophile proteins and their use in protein engineering

We recently introduced ENCoM, an elastic network atomic contact model, as the first coarse-grained normal mode analysis method that accounts for the nature of amino acids and can predict the effect of mutations on thermostability based on changes vibrational entropy. In this proof-of-concept article, we use pairs of mesophile and thermophile homolog proteins with identical structures to determine if a measure of vibrational entropy based on normal mode analysis can discriminate thermophile from mesophile proteins. We observe that in around 60% of cases, thermophile proteins are more rigid at equivalent temperatures than their mesophile counterpart and this difference can guide the design of proteins to increase their thermostability through series of mutations. We observe that mutations separating thermophile proteins from their mesophile orthologs contribute independently to a decrease in vibrational entropy and discuss the application and implications of this methodology to protein engineering.