The protein sequences found in nature represent a tiny fraction of the potential sequences that could be constructed from the 20-amino-acid alphabet. To help define the properties that shaped proteins to stand out from the space of possible alternatives, we conducted a systematic computational and experimental exploration of random (unevolved) sequences in comparison with biological proteins. In our study, combinations of secondary structure, disorder, and aggregation predictions are accompanied by experimental characterization of selected proteins. We found that the overall secondary structure and physicochemical properties of random and biological sequences are very similar. Moreover, random sequences can be well-tolerated by living cells. Contrary to early hypotheses about the toxicity of random and disordered proteins, we found that random sequences with high disorder have low aggregation propensity (unlike random sequences with high structural content) and were particularly well-tolerated. This direct structure content/aggregation propensity dependence differentiates random and biological proteins. Our study indicates that while random sequences can be both structured and disordered, the properties of the latter make them better suited as progenitors (in both in vivo and in vitro settings) for further evolution of complex, soluble, three-dimensional scaffolds that can perform specific biochemical tasks.

The difference spectroscopy of the Raman optical activity (ROA) provides extended information about molecular structure. However, interpretation of the spectra is based on complex and often inaccurate simulations. Previously, the authors attempted to make the calculations more robust by including the solvent and exploring the role of molecular flexibility for alanine and proline zwitterions. In the current study, they analyze the IR, Raman, and ROA spectra of these molecules with the emphasis on the force field modeling. Vibrational harmonic frequencies obtained with 25 ab initio methods are compared to experimental band positions. The role of anharmonic terms in the potential and intensity tensors is also systematically explored using the vibrational self-consistent field, vibrational configuration interaction (VCI), and degeneracy-corrected perturbation calculations. The harmonic approach appeared satisfactory for most of the lower-wavelength (200–1800cm−1) vibrations. Modern generalized gradient approximation and hybrid density functionals, such as the common B3LYP method, provided a very good statistical agreement with the experiment. Although the inclusion of the anharmonic corrections still did not lead to complete agreement between the simulations and the experiment, occasional enhancements were achieved across the entire region of wave numbers. Not only the transitional frequencies of the C–H stretching modes were significantly improved but also Raman and ROA spectral profiles including N–H and C–H lower-frequency bending modes were more realistic after application of the VCI correction. A limited Boltzmann averaging for the lowest-frequency modes that could not be included directly in the anharmonic calculus provided a realistic inhomogeneous band broadening. The anharmonic parts of the intensity tensors (second dipole and polarizability derivatives) were found less important for the entire spectral profiles than the force field anharmonicities (third and fourth energy derivatives), except for a few weak combination bands which were dominated by the anharmonic tensor contributions.