Molecular Simulation of Stapled Peptides.

Constrained peptides represent a relatively new class of biologic therapeutics, which have the potential to overcome several limitations of small-molecule drugs, and of designed antibodies. Because of their modest size, the rational design of such peptides is becoming increasingly amenable to computer simulation; multi-microsecond molecular dynamic (MD) simulations are now routinely possible on consumer-grade graphical processors (GPUs). Here, we describe the procedures for performing and analyzing MD simulations of hydrocarbon-stapled peptides using the CHARMM energy function, in isolation and in complex with a binding partner, to investigate their conformational properties and to compute changes in their binding affinity upon mutation.

Molecular Dynamics-Guided Design of a Functional Protein-ATRP Conjugate That Eliminates Protein-Protein Interactions.

Even the most advanced protein-polymer conjugate therapeutics do not eliminate antibody-protein and receptor-protein recognition. Next-generation bioconjugate drugs will need to replace stochastic selection with rational design to select desirable levels of protein-protein interaction while retaining function. The "Holy Grail" for rational design would be to generate functional enzymes that are fully catalytic with small molecule substrates while eliminating interaction between the protein surface and larger molecules. Using chymotrypsin, an important enzyme that is used to treat pancreatic insufficiency, we have designed a series of molecular chimeras with varied grafting densities and shapes. Guided by molecular dynamic simulations and next-generation molecular chimera characterization with asymmetric flow field-flow fractionation chromatography, we grew linear, branched, and comb-shaped architectures from the surface of the protein by atom-transfer radical polymerization. Comb-shaped polymers, grafted from the surface of chymotrypsin, completely prevented enzyme inhibition with protein inhibitors without sacrificing the ability of the enzyme to catalyze the hydrolysis of a peptide substrate. Asymmetric flow field-flow fractionation coupled with multiangle laser light scattering including dynamic light scattering showed that nanoarmor designed with comb-shaped polymers was particularly compact and spherical. The polymer structure significantly increased protein stability and reduced protein-protein interactions. Atomistic molecular dynamic simulations predicted that a dense nanoarmor with long-armed comb-shaped polymer would act as an almost perfect molecular sieve to filter large ligands from substrates. Surprisingly, a conjugate that was composed of 99% polymer was needed before the elimination of protein-protein interactions.

Atomistic insight towards the impact of polymer architecture and grafting density on structure-dynamics of PEGylated bovine serum albumin and their applications.

Macromolecules such as proteins conjugated to polyethylene glycol (PEG) have been employed in therapeutic drug applications, and recent research has emphasized the potential of varying polymer architectures and conjugation strategies to achieve improved efficacy. In this study, we performed atomistic molecular dynamics simulations of bovine serum albumin (BSA) conjugated to 5 kDa PEG polymers in an array of schemes, including varied numbers of attached chains, grafting density, and nonlinear architectures. Nonlinear architectures included U-shaped PEG, Y-shaped PEG, and poly(oligoethylene glycol methacrylate) (POEGMA). Buried surface area calculations and polymer volume map analyses revealed that volume exclusion behaviors of the high grafting density conjugate promoted additional protein-polymer interactions when compared to simply increasing numbers of conjugated chains uniformly across the protein surface. Investigation of nonlinear polymer architectures showed that stable polymer-lysine loop-like conformations seen in previous conjugate designs were more variable in prevalence, especially in POEGMA, which contained short oligomer PEG chains. The findings of this comprehensive study of alternate PEGylation schemes of BSA provide critical insight into molecular patterns of interaction within bioconjugates and highlight their importance in the future of controlled modification of conjugate system parameters.

Structure-function-dynamics of α-chymotrypsin based conjugates as a function of polymer charge.

The field of protein-polymer conjugates has suffered from a lack of predictive tools and design guidelines to synthesize highly active and stable conjugates. In order to develop this type of information, structure-function-dynamics relationships must be understood. These relationships depend strongly on protein-polymer interactions and how these influence protein dynamics and conformations. Probing nanoscale interactions is experimentally difficult, but computational tools, such as molecular dynamics simulations, can easily obtain atomic resolution. Atomistic molecular dynamics simulations were used to study α-chymotrypsin (CT) densely conjugated with either zwitterionic, positively charged, or negatively charged polymers. Charged polymers interacted with the protein surface to varying degrees and in different regions of the polymer, depending on their flexibilities. Specific interactions of the negatively charged polymer with CT caused structural deformations in CT's substrate binding pocket and active site while no deformations were observed for zwitterionic and positively charged polymers. Attachment of polymers displaced water molecules from CT's surface into the polymer phase and polymer hydration correlated with the Hofmeister series.

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.

Transforming protein-polymer conjugate purification by tuning protein solubility.

Almost all commercial proteins are purified using ammonium sulfate precipitation. Protein-polymer conjugates are synthesized from pure starting materials, and the struggle to separate conjugates from polymer, native protein, and from isomers has vexed scientists for decades. We have discovered that covalent polymer attachment has a transformational effect on protein solubility in salt solutions. Here, protein-polymer conjugates with a variety of polymers, grafting densities, and polymer lengths are generated using atom transfer radical polymerization. Charged polymers increase conjugate solubility in ammonium sulfate and completely prevent precipitation even at 100% saturation. Atomistic molecular dynamic simulations show the impact is driven by an anti-polyelectrolyte effect from zwitterionic polymers. Uncharged polymers exhibit polymer length-dependent decreased solubility. The differences in salting-out are then used to simply purify mixtures of conjugates and native proteins into single species. Increasing protein solubility in salt solutions through polymer conjugation could lead to many new applications of protein-polymer conjugates.

Molecular Insight into the Protein?Polymer Interactions in N-Terminal PEGylated Bovine Serum Albumin.

Therapeutic proteins have increasingly been used in modern medical applications, but their effectiveness is limited by factors such as stability and blood circulation time. Recently, there has been significant research into covalently linking polyethylene glycol polymer chains (PEG) to proteins, known as PEGylation, to mitigate these issues. In this work, an atomistic molecular dynamics study of N-terminal conjugated PEG-BSA (bovine serum albumin) was conducted with varying PEG molecular weights (2, 5, 10, and 20 kDa) to probe PEG-BSA interactions and evaluate the effect of polymer length on dynamics. It was found that the affinity of PEG toward the protein surface increased as a function of PEG molecular weight and that a certain weight (around 10 kDa) was required to promote protein?polymer interactions. Additionally, preferential interactions were monitored through formed contacts and hotspots were identified. PEG chains coordinating in looplike conformations were found near lysine residues. Also, it was found that hydrophobic interactions played an important role in promoting PEG-BSA interactions as the PEG molecular weight increased. The results provide insight into underlying mechanisms behind transitions in PEG conformations and will aid in future design of effective PEGylated drug molecules.

Intramolecular Interactions of Conjugated Polymers Mimic Molecular Chaperones to Stabilize Protein-Polymer Conjugates.

The power and elegance of protein-polymer conjugates has solved many vexing problems for society. Rational design of these complex covalent hybrids depends on a deep understanding of how polymer physicochemical properties impact the conjugate structure-function-dynamic relationships. We have generated a large family of chymotrypsin-polymer conjugates which differ in polymer length and charge, using grafting-from atom-transfer radical polymerization, to elucidate how the polymers influenced enzyme structure and function at pHs that would unfold and inactivate the enzyme. We also used molecular dynamics simulations to deepen our understanding of protein-polymer intramolecular interactions. Remarkably, the data revealed that, contrary to current thoughts on how polymers stabilize proteins, appropriately designed polymers actually stabilize partially unfolded intermediates and assist in refolding to an active conformation. Long, hydrophilic polymers minimized interfacial interactions in partially unfolded conjugates leading to increased stabilization. The design of covalently attached intramolecular biomimetic chaperones that drive protein refolding could have far reaching consequences.

Unraveling Binding Interactions between Human RANKL and Its Decoy Receptor Osteoprotegerin.

Recent studies have revealed the importance and the active contribution of the RANKL/OPG/RANK pathway in many bone diseases including different forms of common osteoporosis. In this study, we present an extensive atomistic molecular dynamic study of the OPG/RANKL system. Within the molecular models, we varied the number of OPG molecules bound to the RANKL trimer and carried out a study to determine how the binding affinity of the OPG/RANKL system changes as a function of OPG concentration. The molecular mechanics Poisson-Boltzmann surface area method was used to analyze binding free energies. It is shown that the binding affinity decreases with increasing numbers of OPG molecules. Additionally, conformational changes of RANKL, interactions between the N-terminus outlier module of OPG with RANKL, and residues that play an important role in the binding of OPG to RANKL trimer were investigated. A probable cause for unfavorable binding for a third OPG molecule was found. Along with the currently available experimental studies, this computational study will be valuable for the comprehensive understanding of OPG/RANKL at the atomistic level.