Molecular Insight into the Steric Shielding Effect of PEG on the Conjugated Staphylokinase: Biochemical Characterization and Molecular Dynamics Simulation

PEGylation is a successful approach to improve potency of a therapeutic protein. The improved therapeutic potency is mainly due to the steric shielding effect of PEG. However, the underlying mechanism of this effect on the protein is not well understood, especially on the protein interaction with its high molecular weight substrate or receptor. Here, experimental study and molecular dynamics simulation were used to provide molecular insight into the interaction between the PEGylated protein and its receptor. Staphylokinase (Sak), a therapeutic protein for coronary thrombolysis, was used as a model protein. Four PEGylated Saks were prepared by site-specific conjugation of 5 kDa/20 kDa PEG to N-terminus and C-terminus of Sak, respectively. Experimental study suggests that the native conformation of Sak is essentially not altered by PEGylation. In contrast, the bioactivity, the hydrodynamic volume and the molecular symmetric shape of the PEGylated Sak are altered and dependent on the PEG chain length and the PEGylation site. Molecular modeling of the PEGylated Saks suggests that the PEG chain remains highly flexible and can form a distinctive hydrated layer, thereby resulting in the steric shielding effect of PEG. Docking analyses indicate that the binding affinity of Sak to its receptor is dependent on the PEG chain length and the PEGylation site. Computational simulation results explain experimental data well. Our present study clarifies molecular details of PEG chain on protein surface and may be essential to the rational design, fabrication and clinical application of PEGylated proteins.     

PASylation Enhances the Stability,Potency, and Plasma Half‐Life of Interferon α‐2a: A Molecular Dynamics Simulation

In this study, the effectiveness of PASylation in enhancing the potency and plasma half‐life of pharmaceutical proteins has been accredited as an alternative technique to the conventional methods such as PEGylation. Proline, alanine, and serine (PAS) chain has shown some advantages including biodegradability improvement and plasma half‐life enhancement while lacking immunogenicity or toxicity. Although some experimental studies have been performed to find the mechanism behind PASylation, the detailed mechanism of PAS effects on the pharmaceutical proteins has remained obscure, especially at the molecular level. In this study, the interaction of interferon α‐2a (IFN) and PAS chain is investigated using molecular dynamics simulation method. Several important parameters including secondary structure, root‐mean‐square distance, and solvent accessible surface area to investigate the stability, bioavailability, and bioactivity of the PASylated protein are studied. The results demonstrate that IFN conformation is not affected critically through PASylation while it results in improvement of the protein stability and bioactivity. Therefore, PASylation can be considered as a proper biological alternative technique to increase the plasma half‐life of the biopharmaceutical proteins through enlarging apparent volume. The proposed simulation represents a computational approach that would provide a basis for the study of PASylated pharmaceutical proteins for different future applications.     

Synthesis and characterization of PEGylated toll like receptor 7 ligands

Toll-like receptor 7 (TLR7) is located in the endosomal compartment of immune cells. Signaling through TLR7, mediated by the adaptor protein MyD88, stimulates the innate immune system and shapes adaptive immune responses. Previously, we characterized TLR7 ligands conjugated to protein, lipid, or poly(ethylene glycol) (PEG). Among the TLR7 ligand conjugates, the addition of PEG chains reduced the agonistic potency. PEGs are safe in humans and widely used for improvement of pharmacokinetics in existing biologics and some low molecular weight compounds. PEGylation could be a feasible method to alter the pharmacokinetics and pharmacodynamics of TLR7 ligands. In this study, we systematically studied the influence of PEG chain length on the in vitro and in vivo properties of potent TLR7 ligands. PEGylation increased solubility of the TLR7 ligands and modulated protein binding. Adding a 6-10 length PEG to the TLR7 ligand reduced its potency toward induction of interleukin (IL)-6 by murine macrophages in vitro and IL-6 and tumor necrosis factor (TNF) in vivo. However, PEGylation with 18 or longer chain restored, and even enhanced, the agonistic activity of the drug. In human peripheral blood mononuclear cells, similar effects of PEGylation were observed for secretion of proinflammatory cytokines, IL-6, IL-12, TNF-α, IL-1β, and type 1 interferon, as well as for B cell proliferation. In summary, these studies demonstrate that conjugation of PEG chains to a synthetic TLR ligand can impact its potency for cytokine induction depending on the size of the PEG moiety. Thus, PEGylation may be a feasible approach to regulate the pharmacological properties of TLR7 ligands.     

Computational insights into the binding pattern of mitochondrial calcium uniporter inhibitor through homology modeling,molecular dynamics simulation,binding free energy prediction and density functional theory calculation

Abstract

The mitochondrial calcium uniporter (MCU) is the critical protein of the inner mitochondrial membrane that is the primary mediator for calcium uptake into the mitochondrial matrix. Herein we built the optimal homology model of human MCU which was refined through all-atom molecular dynamics simulation. Then, the binding mode of known inhibitor was predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the docking result and stability of the protein-inhibitor complex. Finally, density functional theory (DFT) calculation enhanced our understanding of the molecular interaction of MCU inhibitor. Our research would provide a deeper insight into the interactions between human MCU and its inhibitor, which boosts to develop novel therapy against MCU related disease.

Communicated by Ramaswamy H. Sarma     

Thermodynamic study of forces involved in bovine serum albumin and ovalbumin partitioning in aqueous two-phase systems

The partitioning of bovine serum albumin and ovalbumin in different two-phase aqueous polymer systems is investigated using a thermodynamic approach. Systems used were polyethylene glycols (PEGs) of molecular weights 1000 to 10,000 Da and Dextran T500 (500,000 Da). Ovalbumin transfer to the top phase is exothermic, which suggests an electrostatic interaction between the hydroxyl groups of PEG and the hydrophilic side chain of the protein, whereas the bovine serum albumin partition is an endothermic process that is entropically driven, which coincides with its high surface hydrophobicity. The effect of PEG molecular weight on enthalpy and heat capacity changes, associated with the partition of both proteins, is examined on the basis of a preferential interaction of low-molecular-weight PEG with the protein surface.     

PEG修饰对PLL-EGF基因载体靶向结合作用的影响

聚阳离子基因载体系统由于安全性好和便于设计等优点,近年来在基因治疗中的应用发展迅速.在进行基因药物的体内靶向输送时,目前国际上主要通过在基因输送系统中修饰聚乙二醇(PEG)和靶向分子来提高体内输送的稳定性和靶向性.PEG的修饰可能会遮蔽靶向分子的功能呈现,因此建立定量分析方法评价PEG修饰对靶向结合作用的影响非常重要.将连接有表皮生长因子(EGF)的聚赖氨酸(PLL)基因载体作为研究模型,建立BIAcore检测方法,比较PLL-EGF,PEG7000修饰的PLL-EGF,PEG20000修饰的PLL-EGF对表皮生长因子受体(EGFR)的结合和解离速率,评价PEG修饰对PLL-EGF靶向功能呈现的影响.结果表明,PEG7000的修饰降低了EGF和EGFR之间的结合速率,提高了解离速率,整体减弱了靶向分子的靶向结合能力.PEG20000的修饰进一步减弱靶向分子功能的呈现.因此在进行靶向型聚阳离子基因输送系统设计时,考察PEG修饰对靶向结合能力的影响程度非常重要.该研究结果也对其他基因载体系统的设计提供必要的参考.     

Role of hydrophilic and hydrophobic contacts in folding of the second beta-hairpin fragment of protein G: molecular dynamics simulation studies of an all-atom model

Predicting the folding mechanism of the second beta-hairpin fragment of the Ig-binding domain B of streptococcal protein G is unexpectedly challenging for simplified reduced models because the models developed so far indicated a different folding mechanism from what was suggested from high-temperature unfolding and equilibrium free-energy surface analysis based on established all-atom empirical force fields in explicit or implicit solvent. This happened despite the use of empirical residue-based interactions, multibody hydrophobic interactions, and inclusions of hydrogen bonding effects in the simplified models. This article employs a recently developed all-atom (except nonpolar hydrogens) model interacting with simple square-well potentials to fold the peptide fragment by molecular dynamics simulation methods. In this study, 193 out of 200 trajectories are folded at two reduced temperatures (3.5 and 3.7) close to the transition temperature T* approximately 4.0. Each simulation takes <7 h of CPU time on a Pentium 800-MHz PC. Folding of the new all-atom model is found to be initiated by collapse before the formation of main-chain hydrogen bonds. This verifies the mechanism proposed from previous all-atom unfolding and equilibrium simulations. The new model further predicts that the collapse is initiated by two nucleation contacts (a hydrophilic contact between D46 and T49 and a hydrophobic contact between Y45 and F52), in agreement with recent NMR measurements. The results suggest that atomic packing and native contact interactions play a dominant role in folding mechanism.     

Shielding effect of a PEG molecule of a mono-PEGylated peptide varies with PEG chain length

Abstract

‘Shielding’ effect of a conjugated PEG molecule could cause a change in the electrostatic interaction characteristics of a PEGylate. We investigated how PEG chain length (or molecular weight) alters the electrostatic interaction potential of exenatide variants using their mono-PEGylates in a branched and linear form as model PEGylates. First, we performed the experiments to demonstrate the elution time changes of the mono-PEGylates conjugated with various MW PEGs (5, 10, 20, and 40 kD) using cation exchange chromatography (HiTrap^® SP) at various pHs (2.5, 3.0, 3.5, and 4.0). Then, we calculated the net surface charge of each mono-PEGylate to propose the PEG molecule’s shielding range in terms of the number of amino acids adjacent to the conjugation residue, assuming that a PEG molecule in solution sweeps out a spherical space and an exenatide molecule have a secondary structure. The net charge calculation result was well-correlated with the experimental elution time data, where 5, 10, 20, and 40 kD PEG hindered the electrostatic potential of 5, 8, 12, and 17 amino acid residues in maximum, respectively, on each side of the conjugation point.     

Crowding affects structural dynamics and contributes to membrane association of the NS3/4A complex

Using molecular dynamics simulations, we describe how crowded environments affect the internal dynamics and diffusion of the hepatitis C virus proteases NS3/4A. This protease plays a key role in viral replication and is successfully used as a target for antiviral treatment. The NS3 enzyme requires a peptide cofactor, called NS4A, with its central part interacting with the NS3 β-sheet, and flexible, protruding terminal tails that are unstructured in water solution. The simulations describe the enzyme and water molecules at atomistic resolution, whereas crowders are modeled via either all-atom or coarse-grained models to emphasize different aspects of crowding. Crowders reflect the polyethylene glycol (PEG) molecules used in the experiments to mimic the crowded surrounding. A bead-shell model of folded coarse-grained PEG molecules considers mainly the excluded volume effect, whereas all-atom PEG models afford more protein-like crowder interactions. Circular dichroism spectroscopy experiments of the NS4A N-terminal tail show that a helical structure is formed in the presence of PEG crowders. The simulations suggest that crowding may assist in the formation of an NS4A helical fragment, positioned exactly where a transmembrane helix would fold upon the NS4A contact with the membrane. In addition, partially interactive PEGs help the NS4A N-tail to detach from the protease surface, thus enabling the process of helix insertion and potentially helping the virus establish a replication machinery needed to produce new viruses. Results point to an active role of crowding in assisting structural changes in disordered protein fragments that are necessary for their biological function.     

Elucidation of interactions regulating conformational stability and dynamics of SARS-CoV-2 S-protein

The ongoing COVID-19 pandemic caused by the new coronavirus, SARS-CoV-2, calls for urgent developments of vaccines and antiviral drugs. The spike protein of SARS-CoV-2 (S-protein), which consists of trimeric polypeptide chains with glycosylated residues on the surface, triggers the virus entry into a host cell. Extensive structural and functional studies on this protein have rapidly advanced our understanding of the S-protein structure at atomic resolutions, although most of these structural studies overlook the effect of glycans attached to the S-protein on the conformational stability and functional motions between the inactive down and active up forms. Here, we performed all-atom molecular dynamics simulations of both down and up forms of a fully glycosylated S-protein in solution as well as targeted molecular dynamics simulations between them to elucidate key interdomain interactions for stabilizing each form and inducing the large-scale conformational transitions. The residue-level interaction analysis of the simulation trajectories detects distinct amino acid residues and N-glycans as determinants on conformational stability of each form. During the conformational transitions between them, interdomain interactions mediated by glycosylated residues are switched to play key roles on the stabilization of another form. Electrostatic interactions, as well as hydrogen bonds between the three receptor binding domains, work as driving forces to initiate the conformational transitions toward the active form. This study sheds light on the mechanisms underlying conformational stability and functional motions of the S-protein, which are relevant for vaccine and antiviral drug developments.     

How can free energy component analysis explain the difference in protein stability caused by amino acid substitutions? Effect of three hydrophobic mutations at the 56th residue on the stability of human lysozyme

To elucidate the molecular mechanism of thermal stability, it is essential to determine what are the major free energy components that contribute significantly to the total free energy difference caused by amino acid mutations. In this work, we carried out free energy calculations based on all-atom molecular dynamics simulations to investigate the effect of three hydrophobic mutations at the same position, I56A, I56V and I56F of human lysozyme. The calculated free energy differences are in good agreement with the experimental values in all cases. From free energy component analysis, we found that small changes in stability in the I56A and I56V mutants originate from the short-range Lennard-Jones interactions, whereas the I56F mutant is largely destabilized owing to the changes in the long-range electrostatic interactions. The calculated results are also compared with the free energy components determined by an empirical relationship based on the native-state structure and thermodynamic data. Although this relationship has been shown to be very successful in reproducing the stability changes caused by various amino acid substitutions in several proteins, the changes of stability in I56V and I56F mutants are not reproduced very well. By comparing the free energy components calculated by these two approaches, we showed that the effect of the long-range interaction on the stability changes may be underestimated in the empirical relationships when the structural change caused by mutation is relatively small, as in I56F. It is also suggested that estimation of the change in accessible surface area, deltadeltaASA, may be overestimated if the structure around the mutation site in the denatured state is native-like, which would cause overestimation of the free energy change as in the case of I56V. Our results clearly show that the combined approach of the free energy calculation based on the all-atom molecular dynamics simulation and the empirical relationships is very useful for understanding the detailed mechanism of protein stability.     

Human serum albumin interaction,in silico and anticancer evaluation of Pine-Gold nanoparticles

Unique characteristics displayed by phytoconstituent conjugated nanoparticles and their crucial interactions with proteins serve to develop nanoparticle-bio-interface platform. Gold nanoparticles of 16 nm in size were generated using aqueous extracts of pine bark and further conjugated to human serum albumin. The gold nanoparticles-protein complex was characterized by surface plasmon resonance, UV–vis and emission spectroscopy techniques. Further, it was characterized for surface morphology and elemental composition, crystallographic quality, nanoparticles size, shape, stability, structural determination and the identification of capping agent. Moreover, the interaction of gold nanoparticles with human serum albumin was investigated using conventional spectroscopy techniques. Fluorescence quenching and absorption studies demonstrated an effective binding of human serum albumin with oleamide capped gold nanoparticles. The molecular docking study showed a binding affinity of -6.1 kcal/mol whereas the molecular dynamics simulation indicated that the binding of oleamide to human serum albumin. A biological evaluation of pine bark extract-gold nanoparticles showed cytotoxicity with increased cell mortality in lung cancer cells and minimal toxicity on non-cancerous human embryonic kidney cells, respectively.     

Constructing side chains on near-native main chains for ab initio protein structure prediction

Is there value in constructing side chains while searching protein conformational space during an ab initio simulation? If so, what is the most computationally efficient method for constructing these side chains? To answer these questions, four published approaches were used to construct side chain conformations on a range of near-native main chains generated by ab initio protein structure prediction methods. The accuracy of these approaches was compared with a naive approach that selects the most frequently observed rotamer for a given amino acid to construct side chains. An all-atom conditional probability discriminatory function is useful at selecting conformations with overall low all-atom root mean square deviation (r.m.s.d.) and the discrimination improves on sets that are closer to the native conformation. In addition, the naive approach performs as well as more sophisticated methods in terms of the percentage of chi(1) angles built accurately and the all-atom r. m.s.d., between the native and near-native conformations. The results suggest that the naive method would be extremely useful for fast and efficient side chain construction on vast numbers of conformations for ab initio prediction of protein structure.     

Solid-phase PEGylation of recombinant interferon alpha-2a for site-specific modification: process performance, characterization, and in vitro bioactivity

'Solid-phase' PEGylation, in which a conjugation reaction attaches proteins to a solid matrix, has distinct advantages over the conventional, solution-phase process. We report a case study in which recombinant interferon (rhIFN) alpha-2a was adsorbed to a cation-exchange resin and PEGylated at the N-terminus by 5, 10, and 20 kDa mPEG aldehydes through reductive alkylation. After PEGylation, a salt gradient elution efficiently purified the mono-PEGylate of unwanted species such as unmodified IFN and unreacted PEG. Mono-PEGylation and purification were integrated into a single, chromatographic step. Depending on the molecular weight of the mPEG aldehyde, the mono-PEGylation yield ranged from 50 to 65%. Major problems associated with the solution-phase process such as random or uncontrollable multi-PEGylation and post-PEGylation purification difficulties were overcome. N-terminus sequencing and MALDI-TOF mass spectrophometry confirmed that the PEG molecule was conjugated only to the N-terminus. A cell proliferation study indicated reduced antiviral activity of the mono-PEGylate compared to that of the unmodified IFN. As higher molecular weight PEG was conjugated, in vitro bioactivity and antibody binding activity, as measured by a surface plasmon resonance biosensor, decreased. Nevertheless, trypsin resistance and thermal stability were considerably improved .     

Proteolytically stabilizing fibronectin without compromising cell and gelatin binding activity

Excessive proteolytic degradation of fibronectin (FN) has been implicated in impaired tissue repair in chronic wounds. We previously reported two strategies for stabilizing FN against proteolytic degradation; the first conjugated polyethylene glycol (PEG) through cysteine residues and the second conjugated PEG chains of varying molecular weight on lysine residues. PEGylation of FN via lysine residues resulted in increased resistance to proteolysis with increasing PEG size, but an overall decrease in biological activity, as characterized by cell and gelatin binding. Our latest method to stabilize FN against proteolysis masks functional regions in the protein during lysine PEGylation. FN is PEGylated while it is bound to gelatin Sepharose beads with 2, 5, and 10 kDa PEG precursors. This results in partially PEGylated FN that is more stable than native FN and whose proteolytic stability increases with PEG molecular weight. Unlike completely PEGylated FN, partially PEGylated FN has cell adhesion, gelatin binding, and matrix assembly responses that are comparable to native FN. This is new evidence of how PEGylation variables can be used to stabilize FN while retaining its activity. The conjugates developed herein can be used to dissect molecular mechanisms mediated by FN stability and functionality, and address the problem of FN degradation in chronic wounds. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 31:277–288, 2015     

The Influence of 150-Cavity Binders on the Dynamics of Influenza A Neuraminidases as Revealed by Molecular Dynamics Simulations and Combined Clustering

Neuraminidase inhibitors are the main pharmaceutical agents employed for treatments of influenza infections. The neuraminidase structures typically exhibit a 150-cavity, an exposed pocket that is adjacent to the catalytic site. This site offers promising additional contact points for improving potency of existing pharmaceuticals, as well as generating entirely new candidate inhibitors. Several inhibitors based on known compounds and designed to interact with 150-cavity residues have been reported. However, the dynamics of any of these inhibitors remains unstudied and their viability remains unknown. This work reports the outcome of long-term, all-atom molecular dynamics simulations of four such inhibitors, along with three standard inhibitors for comparison. Each is studied in complex with four representative neuraminidase structures, which are also simulated in the absence of ligands for comparison, resulting in a total simulation time of 9.6µs. Our results demonstrate that standard inhibitors characteristically reduce the mobility of these dynamic proteins, while the 150-binders do not, instead giving rise to many unique conformations. We further describe an improved RMSD-based clustering technique that isolates these conformations – the structures of which are provided to facilitate future molecular docking studies – and reveals their interdependence. We find that this approach confers many advantages over previously described techniques, and the implications for rational drug design are discussed.     

Chemical camouflage of nanospheres with a poorly reactive surface: towards development of stealth and target-specific nanocarriers

A two-step approach is described to chemically camouflage the inert surface of model polystyrene nanospheres of 60 nm in diameter against recognition by the body's defenses. The first step was based on the strong protein adsorbing potency of polystyrene, and the second step utilized the chemical reactivity of the adsorbed proteins for conjugation with cyanuric chloride-activated methoxypoly(ethyleneglycol)5000, mPEG5000. Bovine serum albumin (BSA) and rat IgG were used as models of non-immune and immune proteins, respectively. The maximum adsorbance values for both proteins were near expectation for a close-packed monolayer. Adsorption isotherms studies and analysis of the hydrodynamic thickness of the adsorbed protein layer confirmed the close-packed side-on mode of adsorption for BSA and the end-on mode of adsorption for IgG, respectively. Nucleophiles on the adsorbed proteins were then reacted with cyanuric chloride activated mPEG5000. The average poly(ethyleneglycol) (PEG) content for a 60-nm nanospheres was found to be 13.7+/-0.4 micromol PEG/micromol BSA and 3.6+/-0.3 micromol PEG/micromol IgG. The interaction of both PEG-bearing nanospheres with the hydrophobic column material octyl-agarose indicated surface heterogeneity among the nanospheres. Only nanospheres with the most hydrophilic phenotype (approximately 70% of the total population) exhibited stealth properties after intravenous injection to rats. In contrast to the described approach, incubation of uncoated nanospheres with preformed BSA-mPEG5000 conjugates failed to produce long circulating entities. For design of splenotropic particles cyanuric chloride-activated mPEG5000 was conjugated to BSA-coated polystyrene beads of 225 nm in diameter. Despite their stealth property to hepatic Kupffer cell recognition, these nanospheres were cleared by the splenic red pulp macrophages.     

Improved retention of idarubicin after intravenous injection obtained for cholesterol-free liposomes

To date there has been a focus on the application of sterically stabilized liposomes, composed of saturated diacylphospholipid, polyethylene glycol (PEG) conjugated lipids (5-10 mole%) and cholesterol (CH) (>30 mole%), for the systemic delivery of drugs. However, we are now exploring the utility of liposome formulations composed of diacylphospholipid conjugated PEG mixtures prepared in the absence of added cholesterol, with the primary objective of developing formulations that retain encapsulated drug better than comparable formulations prepared with cholesterol. In this report the stability of cholesterol-free distearoylphosphatidylcholine (DSPC):distearoylphosphatidylethanolamine (DSPE)-PEG(2000) (95:5 mol/mol) liposomes was characterized in comparison to cholesterol-containing formulations DSPC:CH (55:45 mol/mol) and DSPC:CH:DSPE-PEG(2000) (50:45:5 mol/mol/mol), in vivo. Circulation longevity of these formulations was determined in consideration of variables that included varying phospholipid acyl chain length, PEG content and molecular weight. The application of cholesterol-free liposomes as carriers for the hydrophobic anthracycline antibiotic, idarubicin (IDA), was assessed. IDA was encapsulated using a transmembrane pH gradient driven process. To determine stability in vivo, pharmacokinetic studies were performed using 'empty' and drug-loaded [(3)H]cholesteryl hexadecyl ether radiolabeled liposomes administered intravenously to Balb/c mice. Inclusion of 5 mole% of DSPE-PEG(2000) or 45 mole% cholesterol to DSPC liposomes increased the mean plasma area under the curve (AUC(0-24h)) 19-fold and 10-fold, respectively. Cryo-transmission electron micrographs of IDA loaded liposomes indicated that the drug formed a precipitate within liposomes. The mean AUC(0-4h) for free IDA was 0.030 micromole h/ml as compared to 1.38 micromole h/ml determined for the DSPC:DSPE-PEG(2000) formulation, a 45-fold increase, demonstrating that IDA was retained better in cholesterol-free compared to cholesterol-containing liposomes.     

Coarse-grained modeling study of nonpeptide RGD ligand density and PEG molecular weight on the conformation of poly(γ-glutamyl-glutamate) paclitaxel conjugates

Molecular shape, flexibility, and surface hydrophilicity are thought to influence the ability of nanoparticles to cross biological barriers during drug delivery. In this study, coarse-grained (CG) molecular dynamics (MD) simulations were used to study these properties of a polymer-drug construct in potential clinical development: poly(γ-glutamyl-glutamate)-paclitaxel-poly(ethylene glycol) nonpeptide RGD (PGG-PTX-PEG-npRGD), a linear glutamyl-glutamate polymer with paclitaxel and poly(ethylene glycol)-nonpeptide RGD side groups. It was hypothesized that the PEG molecular weight (MW) (500 Da; 1,000 Da; and 2,000 Da) and nonpeptide RGD ligand density (4, 8, 12, and 16 per molecule), respectively, may have advantageous effects on the shape, flexibility, and surface hydrophilicity of PGG-PTX-PEG-npRGD. Circular dichroism spectroscopy was used to suggest initial structures for the all-atom (AA) models of PGG-PTX-PEG-npRGD, which were further converted to CG models using a commercially available mapping algorithm. Due to its semi-flexibility, PGG-PTX-PEG-npRGD is not limited to one specific conformation. Thus, CG MD simulations were run until statistical equilibrium, at which PGG-PTX-PEG-npRGD is represented as an ensemble of statistically similar conformations. The size of a PGG-PTX-PEG-npRGD molecule is not affected by the PEG MW or the nonpeptide RGD density, but higher PEG MW results in increased surface density of a PGG-PTX-PEG-npRGD molecule. Most PGG-PTX-PEG-npRGD shapes are globular, although filamentous shapes were also observed in the PEG500 and PEG1000 molecules. PEG500 and PEG1000 molecules are more flexible than PEG2000 systems. A higher presence of npRGD ligands results in decrease surface hydrophilicity of PGG-PTX-PEG-npRGD. These results indicate that the PGG-PTX-PEG1000-npRGD₄ and PGG-PTX-PEG1000-npRGD₈ molecules are the most efficacious candidates and are further recommended for experimental preclinical studies.