Spatially Resolved Nucleated Polymerization: A Free-Boundary Model of Protein Aggregation in Concentrated Solutions
In the authors' words
Kinetic models of protein aggregation describe populations by size, not by spatial organization or morphology. We extend Lumry-Eyring nucleated polymerization to a model in which the monomer is a density field and each aggregate is a region bounded by a level set. Growth is a flux condition on the available sites of a surface. Condensation is a reaction between the bonding sites of two surfaces in contact, at a rate set by the bond rate and the contact geometry. The availability of those sites is a field on the interface, and its equilibrium value follows from Wertheim's perturbation theory. The collision efficiency and the Fuchs stability ratio are therefore computed, not fitted. In a well-mixed limit the model's spatial averages satisfy the rate equations term by term; the monomer fraction agrees to eight parts in ten thousand, a difference that arises from equating aggregate size with volume. The condensation kernel's exponent is against the fitted to a monoclonal antibody. The computed stability ratio reproduces thirteen of fourteen published conditions at twelve , but only with the bond rate at the top of its range. In a many-body box, aggregates merge at to times the two-body rate.
Appeared: Thursday, September 24. arXiv. Preprint, not yet peer-reviewed.
Authors' comment: 36 pages, 9 figures. Supporting Information (189 pages) included as an ancillary file. Submitted to The Journal of Physical Chemistry B