A better understanding of the binding preference of amino acids for gold nanoparticles of different diameters could aid in the design of peptides that bind specifically to nanoparticles of a given ...diameter. Here we identify the binding preference of 19 natural amino acids for three gold nanoparticles with diameters of 1.0, 2.0, and 4.0 nm, and investigate the mechanisms that govern these preferences. We calculate potentials of mean force between 36 entities (19 amino acids and 17 side chains) and the three gold nanoparticles in explicit water using well-tempered metadynamics simulations. Comparing these potentials of mean force determines the amino acids’ nanoparticle binding preferences and if these preferences are controlled by the backbone, the side chain, or both. Twelve amino acids prefer to bind to the 4.0 nm gold nanoparticle, and seven prefer to bind to the 2.0 nm one. We also use atomistic molecular dynamics simulations to investigate how water molecules near the nanoparticle influence the binding of the amino acids. The solvation shells of the larger nanoparticles have higher water densities than those of the smaller nanoparticles while the orientation distributions of the water molecules in the shells of all three nanoparticles are similar. The nanoparticle preferences of the amino acids depend on whether their binding free energy is determined mainly by their ability to replace or to reorient water molecules in the nanoparticle solvation shell. The amino acids whose binding free energy depends mainly on the replacement of water molecules are likely to prefer to bind to the largest nanoparticle and tend to have relatively simple side chain structures. Those whose binding free energy depends mainly on their ability to reorient water molecules prefer a smaller nanoparticle and tend to have more complex side chain structures.
The ability of nanoparticles to alter protein structure and dynamics plays an important role in their medical and biological applications. We investigate allosteric effects of gold nanoparticles on ...human serum albumin protein using molecular simulations. The extent to which bound nanoparticles influence the structure and dynamics of residues distant from the binding site is analyzed. The root mean square deviation, root mean square fluctuation and variation in the secondary structure of individual residues on a human serum albumin protein are calculated for four protein-gold nanoparticle binding complexes. The complexes are identified in a brute-force search process using an implicit-solvent coarse-grained model for proteins and nanoparticles. They are then converted to atomic resolution and their structural and dynamic properties are investigated using explicit-solvent atomistic molecular dynamics simulations. The results show that even though the albumin protein remains in a folded structure, the presence of a gold nanoparticle can cause more than 50% of the residues to decrease their flexibility significantly, and approximately 10% of the residues to change their secondary structure. These affected residues are distributed on the whole protein, even on regions that are distant from the nanoparticle. We analyze the changes in structure and flexibility of amino acid residues on a variety of binding sites on albumin and confirm that nanoparticles could allosterically affect the ability of albumin to bind fatty acids, thyroxin and metals. Our simulations suggest that allosteric effects must be considered when designing and deploying nanoparticles in medical and biological applications that depend on protein-nanoparticle interactions.
Directing the assembly of colloidal particles through the use of external electric or magnetic fields shows promise for the creation of reconfigurable materials. Self-propelled particles can also be ...used to dynamically drive colloidal systems to nonequilibrium steady-states. We investigate colloidal systems that combine both of these methods of directed assembly, simulating mixtures of passive dipolar colloids and active soft spheres in an external magnetic field using Brownian dynamics in two dimensions. In these systems, the dipolar particles align in the direction of the external field but the active particles are unaffected by the field. The phase behaviors exhibited included percolated dipolar-network, dipolar string-fluid, isotropic fluid, and a phase-separated state. We find that the external field allows the dipolar particles to form a percolated network more easily compared to when no external field is present. Additionally, the mixture phase separates at lower active particle velocity in an external field than when no field is present. Our results suggest that combining multiple methods of directing colloidal assembly could lead to new pathways to fabricate reconfigurable materials.
The adsorption of proteins on nanoparticles results in the formation of the protein corona, the composition of which determines how nanoparticles influence their biological surroundings. We seek to ...better understand corona formation by developing models that describe protein adsorption on nanoparticles using computer simulation results as data. Using a coarse-grained protein model, discontinuous molecular dynamics simulations are conducted to investigate the adsorption of two small proteins (Trp-cage and WW domain) on a model nanoparticle of diameter 10.0 nm at protein concentrations ranging from 0.5 to 5 mM. The resulting adsorption isotherms are well described by the Langmuir, Freundlich, Temkin and Kiselev models, but not by the Elovich, Fowler-Guggenheim and Hill-de Boer models. We also try to develop a generalized model that can describe protein adsorption equilibrium on nanoparticles of different diameters in terms of dimensionless size parameters. The simulation results for three proteins (Trp-cage, WW domain, and GB3) on four nanoparticles (diameter = 5.0, 10.0, 15.0, and 20.0 nm) illustrate both the promise and the challenge associated with developing generalized models of protein adsorption on nanoparticles.
Elastin-like polypeptides (ELPs) with the repeat sequence of VPGVG are widely used as a model system for investigation of lower critical solution temperature (LCST) transition behavior. In this ...paper, the effect of temperature on the structure, dynamics and association of (VPGVG)18 in aqueous solution is investigated using atomistic molecular dynamics simulations. Our simulations show that as the temperature increases the ELP backbones undergo gradual conformational changes, which are attributed to the formation of more ordered secondary structures such as β-strands. In addition, increasing temperature changes the hydrophobicity of the ELP by exposure of hydrophobic valine-side chains to the solvent and hiding of proline residues. Based on our simulations, we conclude that the transition behavior of (VPGVG)18 can be attributed to a combination of thermal disruption of the water network that surrounds the polypeptide, reduction of solvent accessible surface area of the polypeptide, and increase in its hydrophobicity. Simulations of the association of two (VPGVG)18 molecules demonstrated that the observed gradual changes in the structural properties of the single polypeptide chain are enough to cause the aggregation of polypeptides above the LCST. These results lead us to propose that the LCST phase behavior of poly(VPGVG) is a collective phenomenon that originates from the correlated gradual changes in single polypeptide structure and the abrupt change in properties of hydration water around the peptide and is a result of a competition between peptide–peptide and peptide–water interactions. This is a computational study of an important intrinsically disordered peptide system that provides an atomic-level description of structural features and interactions that are relevant in the LCST phase behavior.
•Computational protocol to discover peptide ligands that capture Fab-κ and λ.•Five peptides identified via our in-house PepBD algorithm and atomistic MD simulations.•Binding capacities up to 32 mg ...Fab per mL resin and mild Fab binding affinity (KD ∼ 10−5 M).•Dynamic studies returned product yield up to ∼90%.•Lays groundwork for future development of biomanufacturing ligands.
Antibody fragments and their engineered variants show true potential as next-generation therapeutics as they combine excellent targeting with superior biodistribution and blood clearance. Unlike full antibodies, however, antibody fragments do not yet have a standard platform purification process for large-scale production. Short peptide ligands are viable alternatives to protein ligands in affinity chromatography. In this work, an integrated computational and experimental scheme is described to de novo design 9-mer peptides that bind to Fab fragments. The first cohort of designed sequences was tested experimentally using human polyclonal Fab, and the top performing sequence was selected as a prototype for a subsequent round of ligand refinement in silico. The resulting peptides were conjugated to chromatographic resins and evaluated via equilibrium and dynamic binding studies using human Fab-κ and Fab-λ. The equilibrium studies returned values of binding capacities up to 32 mg of Fab per mL of resin with mild affinity (KD ∼ 10−5 M) that are conducive to high product capture and recovery. Dynamic studies returned values of product yield up to ∼90%. Preliminary purification studies provided purities of 83–93% and yields of 11–89%. These results lay the groundwork for future development of these ligands towards biomanufacturing translation.
Prion diseases are infectious neurodegenerative diseases that are capable of cross‐species transmission, thus arousing public health concerns. Seed‐templating propagation of prion protein is believed ...to underlie prion cross‐species transmission pathology. Understanding the molecular fundamentals of prion propagation is key to unravelling the pathology of prion diseases. In this study, we use coarse‐grained molecular dynamics to investigate the seeding and cross‐seeding aggregation of three prion protein fragments PrP(120–144) originating from human (Hu), bank vole (BV), and Syrian hamster (SHa). We find that the seed accelerates the aggregation of the monomer peptides by eliminating the lag phase. The monomer aggregation kinetics are mainly determined by the structure of the seed. The stronger the hydrophobic residues on the seed associate with each other, the higher the probability that the seed recruits monomer peptides to its surface/interface. For cross‐seeding aggregation, we show that Hu has a strong tendency to adopt the conformation of the BV seed and vice versa; the Hu and BV monomers have a weak tendency to adopt the conformation of the SHa seed. These two findings are consistent with Apostol et al.'s experimental findings on PrP(138–143) and partially consistent with Jones et al.'s finding on PrP(23–144). We also identify several conformational mismatches when SHa cross‐seeds BV and Hu peptides, indicating the existence of a cross‐seeding barrier between SHa and the other two sequences. This study sheds light on the molecular mechanism of seed‐templating aggregation of prion protein fragments underlying the sequence‐dependent transmission barrier in prion diseases.
The physical origin of the lower critical solution temperature (LCST) behavior of a variety of fluids, including elastin-like polypeptides (ELPs), has been studied for the past few decades. As is the ...case for polymer solutions, LCST behavior of ELPs is invariably reported for large systems of molecules and is considered evidence for collective behavior. In contrast, we find evidence for properties changes associated with LCST behavior in a single molecule by performing long atomic-level molecular dynamics simulation on the ELP sequences (Val-Pro-Gly-Val-Gly)n for four different length peptides over a wide range of temperatures. We observe a sharp transition in the number of hydrogen bonds between peptide and water and in the number of water molecules within the first hydration shell as temperature rises; this is used to locate the transition temperature. The dependence of the transition temperatures of ELPs on their lengths agrees well with experiments in that both have the same power law exponents. Our simulations reveal that the tendency for pentamers (VPGVG) in ELPs of all lengths to lose H-bonds with water or to gain H-bonds with themselves as temperature rises is independent of the length of the chain in which they are embedded. Thus, the transition temperature of ELPs in pure water is determined by two factors: the hydrogen bonding tendency of the pentamers and the number of pentamers per ELP. Moreover, the hydrogen bonding tendency of pentamers depends only on their sequences, not on the ELP chain length.
Peptides that bind to inorganic materials can be used to functionalize surfaces, control crystallization, or assist in interfacial self-assembly. In the past, inorganic-binding peptides have been ...found predominantly through peptide library screening. While this method has successfully identified peptides that bind to a variety of materials, an alternative design approach that can intelligently search for peptides and provide physical insight for peptide affinity would be desirable. In this work, we develop a computational, physics-based approach to design inorganic-binding peptides, focusing on peptides that bind to the common plastics polyethylene, polypropylene, polystyrene, and poly(ethylene terephthalate). The PepBD algorithm, a Monte Carlo method that samples peptide sequence and conformational space, was modified to include simulated annealing, relax hydration constraints, and an ensemble of conformations to initiate design. These modifications led to the discovery of peptides with significantly better scores compared to those obtained using the original PepBD. PepBD scores were found to improve with increasing van der Waals interactions, although strengthening the intermolecular van der Waals interactions comes at the cost of introducing unfavorable electrostatic interactions. The best designs are enriched in amino acids with bulky side chains and possess hydrophobic and hydrophilic patches whose location depends on the adsorbed conformation. Future work will evaluate the top peptide designs in molecular dynamics simulations and experiment, enabling their application in microplastic pollution remediation and plastic-based biosensors.