Force Response of Polypeptide Chains from Water-Explicit MD Simulations

Quantitative Prediction of Protein-Polyelectrolyte Binding Thermodynamics: Adsorption of Heparin-Analog Polysulfates to the SARS-CoV-2 Spike Protein RBD

Interactions of polyelectrolytes (PEs) with proteins play a crucial role in numerous biological processes, such as the internalization of virus particles into host cells. Although docking, machine learning methods, and molecular dynamics (MD) simulations are utilized to estimate binding poses and binding free energies of small-molecule drugs to proteins, quantitative prediction of the binding thermodynamics of PE-based drugs presents a significant obstacle in computer-aided drug design. This is due to the sluggish dynamics of PEs caused by their size and strong charge-charge correlations. In this paper, we introduce advanced sampling methods based on a force-spectroscopy setup and theoretical modeling to overcome this barrier. We exemplify our method with explicit solvent all-atom MD simulations of the interactions between anionic PEs that show antiviral properties, namely heparin and linear polyglycerol sulfate (LPGS), and the SARS-CoV-2 spike protein receptor binding domain (RBD). Our prediction for the binding free-energy of LPGS to the wild-type RBD matches experimentally measured dissociation constants within thermal energy, k B T, and correctly reproduces the experimental PE-length dependence. We find that LPGS binds to the Delta-variant RBD with an additional free-energy gain of 2.4 k B T, compared to the wild-type RBD, due to the additional presence of two mutated cationic residues contributing to the electrostatic energy gain. We show that the LPGS-RBD binding is solvent dominated and enthalpy driven, though with a large entropy-enthalpy compensation. Our method is applicable to general polymer adsorption phenomena and predicts precise binding free energies and reconfigurational friction as needed for drug and drug-delivery design.

Magnetic tweezers characterization of the entropic elasticity of intrinsically disordered proteins and peptoids

Understanding the conformational behavior of biopolymers is essential to unlocking knowledge of their biophysical mechanisms and functional roles. Single-molecule force spectroscopy can provide a unique perspective on this by exploiting entropic elasticity to uncover key biopolymer structural parameters. A particularly powerful approach involves the use of magnetic tweezers, which can easily generate lower stretching forces (0.1-20 pN). For forces at the low end of this range, the elastic response of biopolymers is sensitive to excluded volume effects, and they can be described by Pincus blob elasticity model that allow robust extraction of the Flory polymer scaling exponent. Here, we detail protocols for the use of magnetic tweezers for force-extension measurements of intrinsically disordered proteins and peptoids. We also discuss procedures for fitting low-force elastic curves to the predictions of polymer physics models to extract key conformational parameters.

Pincus blob elasticity in an intrinsically disordered protein

Understanding the dynamic structure of intrinsically disordered proteins (IDPs) is important to deciphering their biological functions. Here, we exploit precision entropic elasticity measurements to infer the conformational behavior of a model IDP construct formed from the disordered tail of the neurofilament low molecular weight protein. The IDP construct notably displays a low-force power-law elastic regime, consistent with the Pincus blob model, which allows direct extraction of the Flory exponent, [Formula: see text], from the force-extension relationship. We find [Formula: see text] increases with added denaturant, transitioning from a nearly ideal chain to a swollen chain in a manner quantitatively consistent with measurements of IDP dimensions from other experimental techniques. We suggest that measurements of entropic elasticity could be broadly useful in the study of IDP structure.

Searching for Aquamelt Behavior among Silklike Biomimetics during Fibrillation under Flow

In this paper, we elucidate a generic mechanism behind strain-induced phase transition in aqueous solutions of silk-inspired biomimetics by atomistic molecular dynamics simulations. We show the results of modeling of homopeptides polyglycine Gly₃₀ and polyalanine Ala₃₀ and a heteropeptide (Gly-Ala-Gly-Ala-Gly-Ser)₅, i.e., the simplest and yet relevant sequences that could mimic the behavior of natural silk under stress conditions. First, we analyze hydrophobicities of the sequences by calculating the Gibbs free energy of hydration and inspecting the interchain hydrogen bonding and hydration by water. Second, the force-extension profiles are scanned and compared with the results for poly(ethylene oxide), the synthetic polymer for which the aquamelt behavior has been proved recently. Additionally, the conformational transitions of oligopeptides from coiled to extended states are characterized by a generalized order parameter and by the dependence of the solvent-accessible surface area of the chains on applied stretching. Fibrillation itself is surveyed using both the two-dimensional interchain pair correlation function and the SAXS/WAXS patterns for the aggregates formed under stress. These are compared with experimental data found in the literature on fibril structure of silk composite materials doped with oligoalanine peptides. Our results show that tensile stress introduced into aqueous oligopeptide solutions facilitates interchain interactions. The oligopeptides display both a greater resistance to extension as compared to poly(ethylene oxide) and a reduced ability for hydrogen bonding of the stretched chains between oligomers and with water. Fiber formation is proved for all simulated objects, but the most structured one is made of a heteropeptide (Gly-Ala-Gly-Ala-Gly-Ser)₅: For this sequence, we obtain the highest degree of the secondary structure motifs in the fiber. We conclude that this is the most promising candidate among considered sequences to find the aquamelt behavior in further experimental studies.