We present a novel multi-conformation Monte Carlo simulation method that enables the modeling of protein–protein interactions and aggregation in crowded protein solutions. This approach is relevant to a molecular-scale description of realistic biological environments, including the cytoplasm and the extracellular matrix, which are characterized by high concentrations of biomolecular solutes (e.g., 300–400 mg/mL for proteins and nucleic acids in the cytoplasm of Escherichia coli). Simulation of such environments necessitates the inclusion of a large number of protein molecules. Therefore, computationally inexpensive methods, such as rigid-body Brownian dynamics (BD) or Monte Carlo simulations, can be particularly useful. However, as we demonstrate herein, the rigid-body representation typically employed in simulations of many-protein systems gives rise to certain artifacts in protein–protein interactions. Our approach allows us to incorporate molecular flexibility in Monte Carlo simulations at low computational cost, thereby eliminating ambiguities arising from structure selection in rigid-body simulations. We benchmark and validate the methodology using simulations of hen egg white lysozyme in solution, a well-studied system for which extensive experimental data, including osmotic second virial coefficients, small-angle scattering structure factors, and multiple structures determined by X-ray and neutron crystallography and solution NMR, as well as rigid-body BD simulation results, are available for comparison.

Engineering and bioconjugation of proteins is a critically valuable tool that can facilitate a wide range of biophysical and structural studies. The ability to orthogonally tag or label a domain within a multidomain protein may be complicated by undesirable side reactions to noninvolved domains. Furthermore, the advantages of segmental (or domain-specific) isotopic labeling for NMR, or deuteration for neutron scattering or diffraction, can be realized by an efficient ligation procedure. Common methods—expressed protein ligation, protein trans-splicing, and native chemical ligation—each have specific limitations. Here, we evaluated the use of different variants of Staphylococcus aureus sortase A for a range of ligation reactions and demonstrate that conditions can readily be optimized to yield high efficiency (i.e. completeness of ligation), ease of purification, and functionality in detergents. These properties may enable joining of single domains into multidomain proteins, lipidation to mimic posttranslational modifications, and formation of cyclic proteins to aid in the development of nanodisc membrane mimetics. We anticipate that the method for ligating separate domains into a single functional multidomain protein reported here may enable many applications in structural biology.

We demonstrate that the effect of protein crowding is critically dependent on the stability of the protein's hydration shell, which can dramatically vary between different proteins. In the human eye lens, γS-crystallin (γS-WT) forms a densely packed transparent hydrogel with a high refractive index, making it an ideal system for studying the effects of protein crowding. A single point mutation generates the cataract-related variant γS-G18V, dramatically altering the optical properties of the eye lens. This system offers an opportunity to explore fundamental questions regarding the effect of protein crowding, using γS-WT and γS-G18V: (i) how do the diffusion dynamics of hydration water change as a function of protein crowding?; and (ii) upon hydrogel formation of γS-WT, has a dynamic transition occurred generating a single population of hydration water, or do populations of bulk and hydration water coexist? Using localized spin probes, we separately probe the local translational diffusivity of both surface hydration and interstitial water of γS-WT and γS-G18V in solution. Surprisingly, we find that under the influence of hydrogel formation at highly crowded γS-WT concentrations up to 500 mg/mL, the protein hydration shell remains remarkably dynamic, slowing by less than a factor of 2, if at all, compared to that in dilute protein solutions of ∼5 mg/mL. Upon self-crowding, the population of this robust surface hydration water increases, while a significant bulk-like water population coexists even at ∼500 mg/mL protein concentrations. In contrast, surface water of γS-G18V irreversibly dehydrates with moderate concentration increases or subtle alterations to the solution conditions, demonstrating that the effect of protein crowding is highly dependent on the stability of the protein-specific hydration shell. The core function of γS-crystallin in the eye lens may be precisely its capacity to preserve a robust hydration shell, whose stability is abolished by a single G18V mutation.