Antibody humanization is a key step in the preclinical phase of the development of therapeutic antibodies, originally developed and tested in non-human models (most typically, in mouse). The standard technique of Complementarity-Determining Regions (CDR) grafting into human Framework Regions of germline sequences has some important drawbacks, in that the resulting sequences often need further back-mutations to ensure functionality and/or stability. Here we propose a new method to characterize the statistical distribution of the sequences of the variable regions of human antibodies, that takes into account phenotypical correlations between pairs of residues, both within and between chains. We define a “humanness score” of a sequence, comparing its performance in distinguishing human from murine sequences, with that of some alternative scores in the literature. We also compare the score with the experimental immunogenicity of clinically used antibodies. Finally, we use the humanness score as an optimization function and perform a search in the sequence space, starting from different murine sequences and keeping the CDR regions unchanged. Our results show that our humanness score outperforms other methods in sequence classification, and the optimization protocol is able to generate humanized sequences that are recognized as human by standard homology modelling tools.

The authors address the problem of downhill protein folding in the framework of a simple statistical mechanical model, which allows an exact solution for the equilibrium and a semianalytical treatment of the kinetics. Focusing on protein 1BBL, a candidate for downhill folding behavior, and comparing it to the WW domain of protein PIN1, a two-state folder of comparable size, the authors show that there are qualitative differences in both the equilibrium and kinetic properties of the two molecules. However, the barrierless scenario which would be expected if 1BBL were a true downhill folder is observed only at low enough temperature.

Previous research has shown a strong correlation of protein folding rates to the native state geometry, yet a complete explanation for this dependence is still lacking. Here we study the rate-geometry relationship with a simple statistical physics model, and focus on two classes of model geometries, representing ideal parallel and antiparallel structures. We find that the logarithm of the rate shows an almost perfect linear correlation with the “absolute contact order”, but the slope depends on the particular class considered. We discuss these findings in the light of experimental results.