Folding small proteins via annealing stochastic approximation Monte Carlo

Annealing synchronizes the 70S ribosome into a minimum-energy conformation

Researchers commonly anneal metals, alloys, and semiconductors to repair defects and improve microstructures via recrystallization. Theoretical studies indicate that simulated annealing on biological macromolecules helps predict the final structures with minimum free energy. Experimental validation of this homogenizing effect and further exploration of its applications are fascinating scientific questions that remain elusive. Here, we chose the apo-state 70S ribosome from Escherichia coli as a model, wherein the 30S subunit undergoes a thermally driven intersubunit rotation and exhibits substantial structural flexibility as well as distinct free energy. We experimentally demonstrate that annealing at a fast cooling rate enhances the 70S ribosome homogeneity and improves local resolution on the 30S subunit. After annealing, the 70S ribosome is in a nonrotated state with respect to corresponding intermediate structures in unannealed or heated ribosomes. Manifold-based analysis further indicates that the annealed 70S ribosome takes a narrow conformational distribution and exhibits a minimum-energy state in the free-energy landscape. Our experimental results offer a facile yet robust approach to enhance protein stability, which is ideal for high-resolution cryogenic electron microscopy. Beyond structure determination, annealing shows great potential for synchronizing proteins on a single-molecule level and can be extended to study protein folding and explore conformational and energy landscapes.

Distance-dependent hydrophobic-hydrophobic contacts in protein folding simulations

Successful prediction of protein folding from an amino acid sequence is a challenge in computational biology. In order to reveal the geometric constraints that drive protein folding, highlight those constraints kept or missed by distinct lattices and for establishing which class of intra- and inter-secondary structure element interactions is the most relevant for the correct folding of proteins, we have calculated inter-alpha carbon distances in a set of 42 crystal structures consisting of mainly helix, sheet or mixed conformations. The inter-alpha carbon distances were also calculated in several lattice "hydrophobic-polar" models built from the same protein set. We found that helix structures are more prone to form "hydrophobic-hydrophobic" contacts than beta-sheet structures. At a distance lower than or equal to 3.8 Å (very short-range interactions), "hydrophobic-hydrophobic" contacts are almost absent in the native structures, while they are frequent in all the analyzed lattice models. At distances in-between 3.8 and 9.5 Å (short-/medium-range interactions), the best performing lattice for reproducing mainly helix structures is the body-centered-cubic lattice. If protein structures contain sheet portions, lattice performances get worse, with few exceptions observed for double-tetrahedral and body-centered-cubic lattices. Finally, we can observe that ab initio protein folding algorithms, i.e. those based on the employment of lattices and Monte Carlo simulated annealings, can be improved simply and effectively by preventing the generation of "hydrophobic-hydrophobic" contacts shorter than 3.8 Å, by monitoring the "hydrophobic-hydrophobic/polar-polar" contact ratio in short-/medium distance ranges and by using preferentially a body-centered-cubic lattice.