A flexible approach for understanding protein stability

A distance constraint model (DCM) is presented that identifies flexible regions within protein structure consistent with specified thermodynamic condition. The DCM is based on a rigorous free energy decomposition scheme representing structure as fluctuating constraint topologies. Entropy non-additivity is problematic for naive decompositions, limiting the success of heat capacity predictions. The DCM resolves non-additivity by summing over independent entropic components determined by an efficient network-rigidity algorithm. A minimal 3-parameter DCM is demonstrated to accurately reproduce experimental heat capacity curves. Free energy landscapes and quantitative stability-flexibility relationships are obtained in terms of global flexibility. Several connections to experiment are made.     

Modulation of S6 fibrillation by unfolding rates and gatekeeper residues

We present a protein engineering analysis of the fibrillation of a protein from a thermophilic organism, the 101 residue S6 from Thermus thermophilus. When agitated, S6 fibrillates at pH 2.0 in 0.4 M NaCl. Under these solvent conditions, S6 has native-like secondary structure and also unfolds and refolds cooperatively. However, its tertiary structure appears to be more plastic than at neutral pH, and some regions of the protein may be partially unstructured. At 42 degrees C, there is a lag phase of several days after which fibrillation takes place over several hours. Data from the fibrillation behaviour of a comprehensive series of single and double mutants of S6 suggests that several factors control the onset of fibrillation. Firstly, there appears to be a contiguous region of "gatekeeper" residues that inhibit fibrillation, since their truncation significantly reduces the duration of the lag phase. This region overlaps extensively with the partially unstructured region of the protein, suggesting that residues with enhanced flexibility and solvent-accessibility are important for the initiation of fibrillation. Secondly, longer lag phases correlate with faster rates of unfolding. We interpret this to mean that kinetic stability also controls fibrillation but in the sense that the quasi-native state, rather than the denatured state, is the species that participates in nucleation. This implies that fibrillation can also occur from a quasi-native state as opposed to an ensemble of highly fluctuating structures, and highlights the delicate balance between flexibility and structure required to form organized assemblies of polypeptide chains.     

Local propensities and statistical potentials of backbone dihedral angles in proteins

The following three issues concerning the backbone dihedral angles of protein structures are presented. (1) How do the dihedral angles of the 20 amino acids depend on the identity and conformation of their nearest residues? (2) To what extent are the native dihedral angles determined by local (dihedral) potentials? (3) How to build a knowledge-based potential for a residue's dihedral angles, considering the identity and conformation of its nearest residues? We find that the dihedral angle distribution for a residue can significantly depend on the identity and conformation of its adjacent residues. These correlations are in sharp contrast to the Flory isolated-pair hypothesis. Statistical potentials are built for all combinations of residue triplets and depend on the dihedral angles between consecutive residues. First, a low-resolution potential is obtained, which only differentiates between the main populated basins in the dihedral angle density plots. Minimization of the dihedral potential for 125 test proteins reveals that most native alpha-helical residues (89%) and a large fraction of native beta-sheet residues (47%) adopt conformations close to their native one. For native loop residues, the percentage is 48%. It is also found that this fraction is higher for residues away from the ends of alpha or beta secondary structure elements. In addition, a higher resolution potential is built as a function of dihedral angles by a smoothing procedure and continuous functions interpolations. Monte Carlo energy minimization with this potential results in a lower fraction for native beta-sheet residues. Nevertheless, because of the higher flexibility and entropy of beta structures, they could be preferred under the influence of non-local interactions. In general, most alpha-helices and many beta-sheets are strongly determined by the local potential, while the conformations in loops and near the end of beta-sheets are more influenced by non-local interactions.     

A structural similarity analysis of double-helical DNA

A database of the structural properties of all 32,896 unique DNA octamer sequences has been calculated, including information on stability, the minimum energy conformation and flexibility. The contents of the database have been analysed using a variety of Euclidean distance similarity measures. A global comparison of sequence similarity with structural similarity shows that the structural properties of DNA are much less diverse than the sequences, and that DNA sequence space is larger and more diverse than DNA structure space. Thus, there are many very different sequences that have very similar structural properties, and this may be useful for identifying DNA motifs that have similar functional properties that are not apparent from the sequences. On the other hand, there are also small numbers of almost identical sequences that have very different structural properties, and these could give rise to false-positives in methods used to identify function based on sequence alignment. A simple validation test demonstrates that structural similarity can differentiate between promoter and non-promoter DNA. Combining structural and sequence similarity improves promoter recall beyond that possible using either similarity measure alone, demonstrating that there is indeed information available in the structure of double-helical DNA that is not readily apparent from the sequence.     

Histone-like protein HU from Deinococcus radiodurans binds preferentially to four-way DNA junctions

The histone-like protein HU from Escherichia coli is involved in DNA compaction and in processes such as DNA repair and recombination. Its participation in these events is reflected in its ability to bend DNA and in its preferred binding to DNA junctions and DNA with single-strand breaks. Deinococcus radiodurans is unique in its ability to reconstitute its genome from double strand breaks incurred after exposure to ionizing radiation. Using electrophoretic mobility shift assays (EMSA), we show that D.radiodurans HU (DrHU) binds preferentially only to DNA junctions, with half-maximal saturation of 18 nM. In distinct contrast to E.coli HU, DrHU does not exhibit a marked preference for DNA with nicks or gaps compared to perfect duplex DNA, nor is it able to mediate circularization of linear duplex DNA. These unexpected properties identify DrHU as the first member of the HU protein family not to serve an architectural role and point to its potential participation in DNA recombination events. Our data also point to a mechanism whereby differential target site selection by HU proteins is achieved and suggest that the substrate specificity of HU proteins should be expected to vary as a consequence of their individual capacity for inducing the required DNA bend.     

Relative flexibility of DNA and RNA: a molecular dynamics study

State of the art molecular dynamics simulations are used to study the structure, dynamics, molecular interaction properties and flexibility of DNA and RNA duplexes in aqueous solution. Special attention is paid to the deformability of both types of structures, revisiting concepts on the relative flexibility of DNA and RNA duplexes. Our simulations strongly suggest that the concepts of flexibility, rigidity and deformability are much more complex than usually believed, and that it is not always true that DNA is more flexible than RNA.     

X‐ray scattering reveals disordered linkers and dynamic interfaces in complexes and mechanisms for DNA double‐strand break repair impacting cell and cancer biology

Evolutionary selection ensures specificity and efficiency in dynamic metastable macromolecular machines that repair DNA damage without releasing toxic and mutagenic intermediates. Here we examine non‐homologous end joining (NHEJ) as the primary conserved DNA double‐strand break (DSB) repair process in human cells. NHEJ has exemplary key roles in networks determining the development, outcome of cancer treatments by DSB‐inducing agents, generation of antibody and T‐cell receptor diversity, and innate immune response for RNA viruses. We determine mechanistic insights into NHEJ structural biochemistry focusing upon advanced small angle X‐ray scattering (SAXS) results combined with X‐ray crystallography (MX) and cryo‐electron microscopy (cryo‐EM). SAXS coupled to atomic structures enables integrated structural biology for objective quantitative assessment of conformational ensembles and assemblies in solution, intra‐molecular distances, structural similarity, functional disorder, conformational switching, and flexibility. Importantly, NHEJ complexes in solution undergo larger allosteric transitions than seen in their cryo‐EM or MX structures. In the long‐range synaptic complex, X‐ray repair cross‐complementing 4 (XRCC4) plus XRCC4‐like‐factor (XLF) form a flexible bridge and linchpin for DNA ends bound to KU heterodimer (Ku70/80) and DNA‐PKcs (DNA‐dependent protein kinase catalytic subunit). Upon binding two DNA ends, auto‐phosphorylation opens DNA‐PKcs dimer licensing NHEJ via concerted conformational transformations of XLF‐XRCC4, XLF–Ku80, and LigIV^BRCT–Ku70 interfaces. Integrated structures reveal multifunctional roles for disordered linkers and modular dynamic interfaces promoting DSB end processing and alignment into the short‐range complex for ligation by LigIV. Integrated findings define dynamic assemblies fundamental to designing separation‐of‐function mutants and allosteric inhibitors targeting conformational transitions in multifunctional complexes.     

Large‐brained mammals live longer

Many mammals have brains substantially larger than expected for their body size, but the reasons for this remain ambiguous. Enlarged brains are metabolically expensive and require elongated developmental periods, and so natural selection should have favoured their evolution only if they provide counterbalancing advantages. One possible advantage is facilitating the construction of behavioural responses to unusual, novel or complex socio‐ecological challenges. This buffer effect should increase survival rates and favour a longer reproductive life, thereby compensating for the costs of delayed reproduction. Here, using a global database of 493 species, we provide evidence showing that mammals with enlarged brains (relative to their body size) live longer and have a longer reproductive lifespan. Our analysis supports and extends previous findings, accounting for the possible confounding effects of other life history traits, ecological and dietary factors, and phylogenetic autocorrelation. Thus, these findings provide support for the hypothesis that mammals counterbalance the costs of affording large brains with a longer reproductive life.     

AFM visualization of clathrin triskelia under fluid and in air

Atomic force microscopy (AFM) is used to characterize the structure and interactions of clathrin triskelia. Time sequence images of individual, wet triskelia resting on mica surfaces clearly demonstrate conformational fluctuations of the triskelia. AFM of dried samples yields images having nanometric resolution comparable to that obtainable by electron microscopy of shadowed samples. Increased numbers of triskelion dimers and assembly intermediates, as well as structures having dimensions similar to those of clathrin cages, are observed when the triskelia were immersed in a low salt, low pH buffer. These entities have been quantified by AFM protein volume computation.

Structured summary

MINT-7299119, MINT-7299136:Clathrin (uniprotkb:P49951) and Clathrin (uniprotkb:P49951) bind (MI:0407) by atomic force microscopy (MI:0872)     

Role of disulfide bonds in modulating internal motions of proteins to tune their function: molecular dynamics simulation of scorpion toxin Lqh III

A series of 1-ns MD simulations were performed on the scorpion toxin Lqh III in native and disulfide bond broken states. The removal of disulfide bonds has caused hydrogen bond network alteration in the five-residue turn, the long loop, the alpha-helix, the loop connecting strands II and III, and the C-terminal region. In addition and more importantly, it has influenced the amplitude of the fluctuations of five-residue turn, loops, and C-terminal region with a minor effect on the fluctuations of the cysteines in the broken bond sites. These findings suggest that disulfide bonds are not the most important factors in rigidifying their own locations, while they have more important effects at a global scale. Furthermore, our results reveal that disulfide bonds have considerable influence on the functionally important essential modes of motions and the correlations between the motions of the binding site residues. Therefore, we can conclude that disulfide bonds have a crucial role in modulating the function via adjusting the dynamics of scorpion toxin molecules. Although this conclusion cannot be generalized to all peptides and proteins, it demonstrates the importance of more investigations on this aspect of disulfide bond efficacy.