Conformation of the Poly(ethylene Glycol) Chains in DiPEGylated Hemoglobin Specifically Probed by SANS: Correlation with PEG Length and in Vivo Efficiency

Molecular Dynamics Simulations for Rationalizing Polymer Bioconjugation Strategies: Challenges, Recent Developments, and Future Opportunities

The covalent modification of proteins with polymers is a well-established method for improving the pharmacokinetic properties of therapeutically valuable biologics. The conjugated polymer chains of the resulting hybrid represent highly flexible macromolecular structures. As the dynamics of such systems remain rather elusive for established experimental techniques from the field of protein structure elucidation, molecular dynamics simulations have proven as a valuable tool for studying such conjugates at an atomistic level, thereby complementing experimental studies. With a focus on new developments, this review aims to provide researchers from the polymer bioconjugation field with a concise and up to date overview of such approaches. After introducing basic principles of molecular dynamics simulations, as well as methods for and potential pitfalls in modeling bioconjugates, the review illustrates how these computational techniques have contributed to the understanding of bioconjugates and bioconjugation strategies in the recent past and how they may lead to a more rational design of novel bioconjugates in the future.

Research Progress in Oxygen Carrier Design and Application

Ischemia or hypoxia can lead to pathological changes in the metabolism and function of tissues and then lead to various diseases. Timely and effective blood resuscitation or improvement of hypoxia is very important for the treatment of diseases. However, there is a need to develop stable, nontoxic, and immunologically inert oxygen carriers due to limitations such as blood shortages, different blood types, and the risk of transmitting infections. With the development of various technologies, oxygen carriers based on hemoglobin and perfluorocarbon have been widely studied in recent years. This paper reviews the development and application of hemoglobin and perfluorocarbon oxygen carriers. The design of oxygen carriers was analyzed, and their application as blood substitutes or oxygen carriers in various hypoxic diseases was discussed. Finally, the characteristics and future research of ideal oxygen carriers were prospected to provide reference for follow-up research.

Molecular Sieving with PEGylated Dendron-Protein Conjugates

A series of generation 3-5 dendrons based on a bis(2,2-hydroxymethylpropionic acid) (bis-MPA) scaffold bearing three respective lengths of linear poly(ethylene glycol) at their periphery and a dibenzocyclooctyne unit at their core was prepared. These dendrons were appended to the surface of azide-decorated α-chymotrypsin (α-CT) via strain-promoted azide-alkyne cycloaddition to yield a library of dendron-protein conjugates. These conjugates were characterized by FT-IR and NMR spectroscopy and were imaged using cryo-electron microscopy. The activity of the PEGylated α-CT-dendron conjugates was investigated using a small molecule (benzoyl-l-tyrosine p-nitroanilide) as well as different proteins of different sizes and crystallinities (casein and bovine serum albumin) as substrates. It was found that the activity of the conjugates toward the small molecule was largely retained, while the activity toward the proteins was significantly diminished. Furthermore, the results indicate that for most of the conjugates the PEG length had a more pronounced impact on enzyme activity than the dendron generation. Overall, the highest sieving ratios were found for α-CT-dendron conjugates decorated with G3-PEG2000, G4-PEG2000, and G5-PEG1000, with the latter two structures offering the best combination of sieving ratio and small molecule activity.

Atomic-Scale View of Protein-PEG Interactions that Redirect the Thermal Unfolding Pathway of PEGylated Human Galectin-3

PEGylation is a promising approach to address the central challenge of applying biologics, i.e., lack of protein stability in the demanding environment of the human body. Wider application is hindered by lack of atomic level understanding of protein-PEG interactions, preventing design of conjugates with predicted properties. We deployed an integrative structural and biophysical approach to address this critical challenge with the PEGylated carbohydrate recognition domain of human galectin-3 (Gal3C), a lectin essential for cell adhesion and potential biologic. PEGylation dramatically increased Gal3C thermal stability, forming a stable intermediate and redirecting its unfolding pathway. Structural details revealed by NMR pointed to a potential role of PEG localization facilitated by charged residues. Replacing these residues subtly altered the protein-PEG interface and thermal unfolding behavior, providing insight into rationally designing conjugates while preserving PEGylation benefits.

Moving Protein PEGylation from an Art to a Data Science

PEGylation is a well-established and clinically proven half-life extension strategy for protein delivery. Protein modification with amine-reactive poly(ethylene glycol) (PEG) generates heterogeneous and complex bioconjugate mixtures, often composed of several PEG positional isomers with varied therapeutic efficacy. Laborious and costly experiments for reaction optimization and purification are needed to generate a therapeutically useful PEG conjugate. Kinetic models which accurately predict the outcome of so-called “random” PEGylation reactions provide an opportunity to bypass extensive wet lab experimentation and streamline the bioconjugation process. In this study, we propose a protein tertiary structure-dependent reactivity model that describes the rate of protein-amine PEGylation and introduces “PEG chain coverage” as a tangible metric to assess the shielding effect of PEG chains. This structure-dependent reactivity model was implemented into three models (linear, structure-based, and machine-learned) to gain insight into how protein-specific molecular descriptors (exposed surface areas, pK_a, and surface charge) impacted amine reactivity at each site. Linear and machine-learned models demonstrated over 75% prediction accuracy with butylcholinesterase. Model validation with Somavert, PEGASYS, and phenylalanine ammonia lyase showed good correlation between predicted and experimentally determined degrees of modification. Our structure-dependent reactivity model was also able to simulate PEGylation progress curves and estimate “PEGmer” distribution with accurate predictions across different proteins, PEG linker chemistry, and PEG molecular weights. Moreover, in-depth analysis of these simulated reaction curves highlighted possible PEG conformational transitions (from dumbbell to brush) on the surface of lysozyme, as a function of PEG molecular weight.

Unraveling the Molar Mass Dependence of Shearing-Induced Aggregation Structure of a High-Mobility Polymer Semiconductor

Aggregation-structure formation of conjugated polymers is a fundamental problem in the field of organic electronics and remains poorly understood. Herein, the molar mass dependence of the aggregation structure of a high-mobility conjugated copolymer (TDPP-Se) comprising thiophene-flanked diketopyrrolopyrrole and selenophene is thoroughly shown. Five batches of TDPP-Se are prepared with the number-average molecular weights (M_n ) varied greatly from 21 to 135 kg mol^-1 . Small-angle neutron scattering and transmission electron microscopy are combined to probe the solution structure of these polymers, consistently using a deuterated solvent. All the polymers adopt 1D rod-like aggregation structures and the radius of the 1D rods is not sensitive to the M_n , while the length increases monotonically with M_n . By utilizing the ordered packing of the aggregated structure in solution, a highly aligned and ordered film is prepared and, thereafter, a reliable hole mobility of 13.8 cm² V^-1 s^-1 is realized in organic thin-film transistors with the moderate M_n batch via bar coating. The hole mobility is among the highest values reported for diketopyrrolopyrrole-based polymers. This work paves the way to visualize the real aggregated structure of polymer semiconductors in solution and sheds light on the microstructure control of high-performance electronic devices.

From hemoglobin allostery to hemoglobin-based oxygen carriers

Hemoglobin (Hb) plays its vital role through structural and functional properties evolutionarily optimized to work within red blood cells, i.e., the tetrameric assembly, well-defined oxygen affinity, positive cooperativity, and heterotropic allosteric regulation by protons, chloride and 2,3-diphosphoglycerate. Outside red blood cells, the Hb tetramer dissociates into dimers, which exhibit high oxygen affinity and neither cooperativity nor allosteric regulation. They are prone to extravasate, thus scavenging endothelial NO and causing hypertension, and cause nephrotoxicity. In addition, they are more prone to autoxidation, generating radicals. The need to overcome the adverse effects associated with cell-free Hb has always been a major hurdle in the development of substitutes of allogeneic blood transfusions for all clinical situations where blood is unavailable or cannot be used due to, for example, religious objections. This class of therapeutics, indicated as hemoglobin-based oxygen carriers (HBOCs), is formed by genetically and/or chemically modified Hbs. Many efforts were devoted to the exploitation of the wealth of biochemical and biophysical information available on Hb structure, function, and dynamics to design safe HBOCs, overcoming the negative effects of free plasma Hb. Unfortunately, so far, no HBOC has been approved by FDA and EMA, except for compassionate use. However, the unmet clinical needs that triggered intensive investigations more than fifty years ago are still awaiting an answer. Recently, HBOCs "repositioning" has led to their successful application in organ perfusion fluids.

Diffusivelike Motions in a Solvent-Free Protein-Polymer Hybrid

The interaction between proteins and hydration water stabilizes protein structure and promotes functional dynamics, with water translational motions enabling protein flexibility. Engineered solvent-free protein-polymer hybrids have been shown to preserve protein structure, function, and dynamics. Here, we used neutron scattering, protein and polymer perdeuteration, and molecular dynamics simulations to explore how a polymer dynamically replaces water. Even though relaxation rates and vibrational properties are strongly modified in polymer coated compared to hydrated proteins, liquidlike polymer dynamics appear to plasticize the conjugated protein in a qualitatively similar way as do hydration-water translational motions.

Molecular simulation of zwitterionic polypeptides on protecting glucagon-like peptide-1 (GLP-1)

Owing to their anti-fouling properties, zwitterionic polypeptides demonstrate great advantage on protecting protein drugs. When conjugated to glucagon-like peptide-1 (GLP-1), a drug for type-II diabetes, zwitterionic polypeptides confer better pharmacokinetics than uncharged counterparts. However, its microscopic mechanism is still unclear due to the complicated conformational space. To address this challenge, this work explored the interaction modes of GLP-1 with the unconnected repeat units, instead of the full-length polypeptides. The three repeat units are two zwitterionic pentapeptides VPKEG and VPREG, and one uncharged control VPGAG. Our molecular simulations revealed that the helical conformation of GLP-1 was stabilized by adding 40 polypeptides. Both VPGAG and VPREG formed dense packing shells around GLP-1, but the driving forces were hydrophobic and electrostatic interactions, respectively. In contrast, the packing shell composed of VPKEG was most loose, while could still stabilize GLP-1. The moderate electrostatic interactions endowed VPKEG an anti-fouling property, thereby avoiding non-specific interaction with other amino acids. The strong electrostatic interactions exerted by arginine promoted atomic contacts between VPREG and other residues, making it as "hydrophobic" as VPGAG. In summary, the combination of hydrophobic and moderate electrostatic interactions in VPKEG brings about a subtle balance between stabilizing GLP-1 and avoiding non-specific interaction.

Protocol for Simulations of PEGylated Proteins with Martini 3

Enhancement of proteins by PEGylation is an active area of research. However, the interactions between polymer and protein are far from fully understood. To gain a better insight into these interactions or even make predictions, molecular dynamics (MD) simulations can be applied to study specific protein-polymer systems at molecular level detail. Here we present instructions on how to simulate PEGylated proteins using the latest iteration of the Martini coarse-grained (CG) force-field. CG MD simulations offer near atomistic information and at the same time allow to study complex biological systems over longer time and length scales than fully atomistic-level simulations.

Insight into Protein-Polymer Conjugate Relaxation Dynamics: The Importance of Polymer Grafting

The bio and chemical physics of protein-polymer conjugates are related to parameters that characterize each component. With this work, it is intended to feature the dynamical properties of the protein-polymer conjugate myoglobin (Mb)-poly(ethyl ethylene phosphate), in the ps and ns time scales, in order to understand the respective roles of the protein and of the polymer size in the dynamics of the conjugate. Elastic and quasi-elastic neutron scattering is performed on completely hydrogenated samples with variable number of polymer chains covalently attached to the protein. The role of the polymer length in the protein solvation and internal dynamics is investigated using two conjugates formed by polymers of different molecular weight. It is confirmed that the flexibility of the complex increases with the number of grafted polymer chains and that a sharp dynamical transition appears when either grafting density or polymer molecular weight are high. It is shown that protein size is crucial for the polymer structural organization and interaction on the protein surface and it is established that the glass properties of the polymer change upon conjugation. The results give a better insight of the equivalence of the polymer coating and the role of water on the surface of proteins.

PEGylation within a confined hydrophobic cavity of a protein

The conjugation of polyethylene glycol (PEG) to proteins, known as PEGylation, has increasingly been employed to expand the efficacy of therapeutic drugs. Recently, research has emphasized the effect of the conjugation site on protein-polymer interactions. In this study, we performed atomistic molecular dynamics (MD) simulations of lysine 116 PEGylated bovine serum albumin (BSA) to illustrate how conjugation near a hydrophobic pocket affects the conjugate's dynamics and observed altered low mode vibrations in the protein. MD simulations were performed for a total of 1.5 μs for each PEG chain molecular mass from 2 to 20 kDa. Analysis of preferential PEG-BSA interactions showed that polymer behavior was also affected as proximity to the attractive protein surface patches promoted interactions in small (2 kDa) PEG chains, while the confined environment of the conjugation site reduced the expected BSA surface coverage when the polymer molecular mass increased to 10 kDa. This thorough analysis of PEG-BSA interactions and polymer dynamics increases the molecular understanding of site-specific PEGylation and enhances the use of protein-polymer conjugates as therapeutics.

Metal-Organic Framework Encapsulating Hemoglobin as a High-Stable and Long-Circulating Oxygen Carriers to Treat Hemorrhagic Shock

As an oxygen-transporting protein, free hemoglobin (Hb) often suffers from the disadvantage of undesirable stability and short blood circulation, which severely impairs the potential clinical applications as the blood substitute. In this work, Hb was facilely encapsulated into a kind of metal-organic frameworks (MOFs) (ZIF-8) inspired by the natural biomineralization process. The obtained ZIF-8 encapsulating Hb (ZIF-8@Hb) showed the small hydrodynamic size of 180.8 nm and neutral zeta potential of -2.1 mV by adjusting the ratio of Hb in ZIF-8 frameworks. Intriguingly, Hb encapsulated by ZIF-8 exhibited significantly enhanced stability in alkaline, oxidation, high temperature, or enzymatic environment compared with free Hb because of the excellent protective MOF coatings. More importantly, the negative charge of Hb neutralized the original positive charge of ZIF-8, which led to the better biocompatibility, lower protein adsorption, and macrophage uptake of ZIF-8@Hb than bare ZIF-8 nanoparticles. Furthermore, ZIF-8@Hb displayed extended blood circulation with the elimination half-life of 13.9 h as well as reduced nonspecific distribution in normal organs compared with free Hb or ZIF-8 nanoparticles. With the above advantages, ZIF-8@Hb showed significantly extended survival time of mice in a disease model of hemorrhagic shock compared with free Hb or bare ZIF-8 nanoparticles. Overall, this work offers a high-stable and long-circulating oxygen carrier platform, which may find wide applications as a blood substitute to treat various oxygen-relevant diseases.

Human Dystrophin Structural Changes upon Binding to Anionic Membrane Lipids

Scaffolding proteins play important roles in supporting the plasma membrane (sarcolemma) of muscle cells. Among them, dystrophin strengthens the sarcolemma through protein-lipid interactions, and its absence due to gene mutations leads to the severe Duchenne muscular dystrophy. Most of the dystrophin protein consists of a central domain made of 24 spectrin-like coiled-coil repeats (R). Using small angle neutron scattering (SANS) and the contrast variation technique, we specifically probed the structure of the three first consecutive repeats 1-3 (R1-3), a part of dystrophin known to physiologically interact with membrane lipids. R1-3 free in solution was compared to its structure adopted in the presence of phospholipid-based bicelles. SANS data for the protein/lipid complexes were obtained with contrast-matched bicelles under various phospholipid compositions to probe the role of electrostatic interactions. When bound to anionic bicelles, large modifications of the protein three-dimensional structure were detected, as revealed by a significant increase of the protein gyration radius from 42 ± 1 to 60 ± 4 Å. R1-3/anionic bicelle complexes were further analyzed by coarse-grained molecular dynamics simulations. From these studies, we report an all-atom model of R1-3 that highlights the opening of the R1 coiled-coil repeat when bound to the membrane lipids. This model is totally in agreement with SANS and click chemistry/mass spectrometry data. We conclude that the sarcolemma membrane anchoring that occurs during the contraction/elongation process of muscles could be ensured by this coiled-coil opening. Therefore, understanding these structural changes may help in the design of rationalized shortened dystrophins for gene therapy. Finally, our strategy opens up new possibilities for structure determination of peripheral and integral membrane proteins not compatible with different high-resolution structural methods.

Unraveling the Solution-State Supramolecular Structures of Donor-Acceptor Polymers and their Influence on Solid-State Morphology and Charge-Transport Properties

Polymer self-assembly in solution prior to film fabrication makes solution-state structures critical for their solid-state packing and optoelectronic properties. However, unraveling the solution-state supramolecular structures is challenging, not to mention establishing a clear relationship between the solution-state structure and the charge-transport properties in field-effect transistors. Here, for the first time, it is revealed that the thin-film morphology of a conjugated polymer inherits the features of its solution-state supramolecular structures. A "solution-state supramolecular structure control" strategy is proposed to increase the electron mobility of a benzodifurandione-based oligo(p-phenylene vinylene) (BDOPV)-based polymer. It is shown that the solution-state structures of the BDOPV-based conjugated polymer can be tuned such that it forms a 1D rod-like structure in good solvent and a 2D lamellar structure in poor solvent. By tuning the solution-state structure, films with high crystallinity and good interdomain connectivity are obtained. The electron mobility significantly increases from the original value of 1.8 to 3.2 cm² V^-1 s^-1 . This work demonstrates that "solution-state supramolecular structure" control is critical for understanding and optimization of the thin-film morphology and charge-transport properties of conjugated polymers.