RESUMEN
A major cause of prion infectivity is the early formation of small, fibril-like aggregates consisting of the heptapeptide GNNQQNY. The prion aggregates exhibit a unique stacking mode in which the hydrophobic tyrosine (Y) is exposed outward, forming a bilayer ß-sheet-stacking zipper structure. This stacking mode of the prion peptides, termed "Y-outward" structure for convenience, goes against the common understanding that, for other amyloid-forming peptides, the hydrophobic residues should be hidden within the peptide fibril, referred to as "Y-inward" structure. To explore the extraordinary stacking behaviors of the prion GNNQQNY peptides, two fibril models are constructed in a fashion of "Y-outward" and "Y-inward" stackings and then studied in silico to examine their thermodynamic stabilities and disaggregation pathways. The "Y-inward" structure indeed exhibits stronger thermodynamic stability than the "Y-outward" structure, according to potential energy and stacking energy calculations. To show how the peptide fibrils dissociate, we illustrated two disaggregation pathways. A dihedral-based free energy landscape was then calculated to examine the conformational degrees of freedom of the GNNQQNY chains in the "Y-outward" and "Y-inward" structures. Peptide chains lose more configurational entropy in the "Y-inward" structure than in the "Y-outward" structure, indicating that the prion peptides are prone to aggregate in a fashion of "Y-outward" stacking pattern due to its low conformational constraints. The prion-like aggregation of the GNNQQNY peptides into amyloid fibrils is primarily governed by the configuration entropy.
Asunto(s)
Nanofibras , Priones , Priones/química , Entropía , Estructura Secundaria de Proteína , Péptidos/química , Amiloide/química , Proteínas AmiloidogénicasRESUMEN
New-generation lithium-ion batteries use ionic liquids (ILs) as electrolyte solutions, greatly enhancing the safety and energy storage capacity of the battery. Fundamental molecular insights are useful for understanding the advantages of high conductivity of IL solvent electrolytes over organic solvent ones. In this work, we computationally studied two organic solvents (DMC and DEC) and four IL solvents ([Cnmim][BF4] and [Cnmim][TFSI] (n = 2, 4)) to examine the physicochemical properties of high concentration electrolytes. As expected, the IL solvent electrolytes exhibit higher density and viscosity, and larger self-diffusion coefficients and conductivity than the organic solvent electrolytes. Further, the microstructures of the lithium salt LiTFSI in various solvent electrolytes were investigated to explore the effect of the organic and IL solvents on the ionic association of the ions Li+ and TFSI-. The structural analysis of LiTFSI revealed that the organic solvents restrict the free motion of the ions, reducing the conductivity of the electrolytes. The [BF4]-type IL electrolytes have higher conductivity than the [TFSI]-type IL electrolytes, especially [C4mim][BF4] with the highest conductivity among the IL-based electrolytes. More importantly, it was proved that the dissolution of LiTFSI in the IL solvents is an anion-driven process.
RESUMEN
Targeted therapy for cancer requires thermosensitive components in drug carriers for controlled drug release against viral cells. The conformational transition characteristic of leucine zipper-structured lipopeptides is utilized in our lab to modulate the phase transition temperature of liposomes, thus achieving temperature-responsive control. In this study, we computationally examined the conformational transition behaviors of leucine zipper-structured lipopeptides that were modified at the N-terminus by distinct functional groups. The conformational transition temperatures of these lipopeptides were determined by structural analysis of the implicit-solvent replica exchange molecular dynamics simulation trajectories using the dihedral angle principal component analysis and the dictionary of protein secondary structure method. Our calculations revealed that the computed transition temperatures of the lipopeptides are in good agreement with the experimental measurements. The effect of hydrogen bonds on the conformational stability of the lipopeptide dimers was examined in conventional explicit-solvent molecular dynamics simulations. A quantitative correlation of the degree of structural dissociation of the dimers and their binding strength is well described by an exponential fit of the binding free energies to the conformation transition temperatures of the lipopeptides.
Asunto(s)
Lipopéptidos/química , Liposomas/química , Transición de Fase , Temperatura , Enlace de Hidrógeno , Leucina Zippers , Simulación de Dinámica Molecular , Estructura Secundaria de ProteínaRESUMEN
A computational strategy that integrates our peptide search algorithm with atomistic molecular dynamics simulation was used to design rational peptide drugs that recognize and bind to the anticodon stem and loop domain (ASL(Lys3)) of human tRNAUUULys3 for the purpose of interrupting HIV replication. The score function of the search algorithm was improved by adding a peptide stability term weighted by an adjustable factor λ to the peptide binding free energy. The five best peptide sequences associated with five different values of λ were determined using the search algorithm and then input in atomistic simulations to examine the stability of the peptides' folded conformations and their ability to bind to ASL(Lys3). Simulation results demonstrated that setting an intermediate value of λ achieves a good balance between optimizing the peptide's binding ability and stabilizing its folded conformation during the sequence evolution process, and hence leads to optimal binding to the target ASL(Lys3). Thus, addition of a peptide stability term significantly improves the success rate for our peptide design search.
Asunto(s)
Biología Computacional/métodos , Diseño de Fármacos , Simulación de Dinámica Molecular , Pliegue de Proteína , Estabilidad Proteica , Proteínas de Unión al ARN , Algoritmos , Humanos , Péptidos/química , Péptidos/metabolismo , Unión Proteica , ARN de Transferencia de Lisina/química , ARN de Transferencia de Lisina/metabolismo , Proteínas de Unión al ARN/química , Proteínas de Unión al ARN/metabolismo , TermodinámicaRESUMEN
Our previously developed peptide-design algorithm was improved by adding an energy minimization strategy which allows the amino acid sidechains to move in a broad configuration space during sequence evolution. In this work, the new algorithm was used to generate a library of 21-mer peptides which could substitute for λ N peptide in binding to boxB RNA. Six potential peptides were obtained from the algorithm, all of which exhibited good binding capability with boxB RNA. Atomistic molecular dynamics simulations were then conducted to examine the ability of the λ N peptide and three best evolved peptides, viz. Pept01, Pept26, and Pept28, to bind to boxB RNA. Simulation results demonstrated that our evolved peptides are better at binding to boxB RNA than the λ N peptide. Sequence searches using the old (without energy minimization strategy) and new (with energy minimization strategy) algorithms confirm that the new algorithm is more effective at finding good RNA-binding peptides than the old algorithm. © 2016 Wiley Periodicals, Inc.
Asunto(s)
Algoritmos , Simulación de Dinámica Molecular , Péptidos/química , ARN Viral/química , Proteínas de Unión al ARN/química , Biblioteca de PéptidosRESUMEN
Recent experiments in our lab (Phys. Chem. Chem. Phys., 2016, 18, 10129-10137) suggested using leucine zipper peptides to enhance the thermosensitivity of liposomes. To understand the mechanisms of temperature-responsive control by the leucine zipper peptide in liposomes, we firstly performed quantum mechanics calculations and implicit-solvent replica exchange molecular dynamics simulations to study the thermo-stability of two leucine zipper peptides, CH3(CH2)4-CO-[VAQLEVK-VAQLESK-VSKLESK-VSSLESK] (termed the capped peptide) and A-[VAQLEVK-VAQLESK-VSKLESK-VSSLESK] (termed the ALA peptide). The analysis of dihedral angle principal components and protein secondary structures was conducted to determine the temperature-dependence conformation transition of the two peptides. Simulation results revealed that our computed transition temperature of the capped peptide is 319.1 K that accords with experimental measurement, 321.1 K. Later, explicit-solvent conventional molecular dynamics simulations were carried out to examine the process of folding and unfolding of the ALA and capped peptides complexed with a lipid bilayer and water in the vicinity of their transition temperatures. A further analysis of conformation and energy of the folded peptides showed that the increase of temperature gives rise to a notable decrease in the number of intra-chain hydrogen bonds and a significant increase in the potential energy of the peptides, thereby reducing the folding stability of the two peptides. As compared to the ALA peptide, a lower transition temperature caused by less intra-chain hydrogen bonds was observed in the capped peptide, which is closer to the temperature of tumor cells. This fact suggests that the capped peptide is more suitable to produce highly sensitive liposomes for the delivery of cancer drugs.
RESUMEN
We extend LIME, an intermediate resolution, implicit solvent model for phospholipids previously used in discontinuous molecular dynamics simulations of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayer formation at 325 K, to the description of the geometry and energetics of 1,2-distearoyl-sn-glycero-3-phospho-L-serine (DSPS) and 1,2-dihenarachidoyl-sn-glycero-3-phosphocholine (21PC) and mixtures thereof at both neutral and low pH at 310 K. A multiscale modeling approach is used to calculate the LIME parameters from atomistic simulation data on a mixed DPPC/DSPS system at different pH values. In the model, 17 coarse-grained sites represent DSPS and 18 coarse-grained sites represent 21PC. Each of these coarse-grained sites is classified as 1 of 9 types. LIME/DMD simulations of equimolar bilayers show the following: (1) 21PC/DSPS bilayers with and without surface area restrictions separate faster at low pH than at neutral pH, (2) 21PC/DSPS systems separate at approximately the same rate regardless of whether they are subjected to surface area restrictions, and (3) bilayers with a molar ratio of 9:1 (21PC:DSPS) phase separate to form heterogeneous domains faster at low pH than at neutral pH. Our results are consistent with experimental findings of Sofou and co-workers (Bandekar et al. Mol. Pharmaceutics, 2013, 10, 152-160; Karve et al. Biomaterials, 2010, 31, 4409-4416) that more doxorubicin is released from 21PC/DSPS liposomes at low pH than at neutral pH, presumably because greater phase separation is achieved at low pH than at neutral pH. These are the first molecular-level simulations of the phase separation in mixed lipid bilayers induced by a change in pH.
Asunto(s)
Membrana Dobles de Lípidos/química , Fosfatidilcolinas/química , Fosfoserina/química , Concentración de Iones de Hidrógeno , Cinética , Simulación de Dinámica Molecular , Temperatura , TermodinámicaRESUMEN
Human tRNA(Lys3)UUU is the primer for HIV replication. The HIV-1 nucleocapsid protein, NCp7, facilitates htRNA(Lys3)UUU recruitment from the host cell by binding to and remodeling the tRNA structure. Human tRNA(Lys3)UUU is post-transcriptionally modified, but until recently, the importance of those modifications in tRNA recognition by NCp7 was unknown. Modifications such as the 5-methoxycarbonylmethyl-2-thiouridine at anticodon wobble position-34 and 2-methylthio-N(6)-threonylcarbamoyladenosine, adjacent to the anticodon at position-37, are important to the recognition of htRNA(Lys3)UUU by NCp7. Several short peptides selected from phage display libraries were found to also preferentially recognize these modifications. Evolutionary algorithms (Monte Carlo and self-consistent mean field) and assisted model building with energy refinement were used to optimize the peptide sequence in silico, while fluorescence assays were developed and conducted to verify the in silico results and elucidate a 15-amino acid signature sequence (R-W-Q/N-H-X2-F-Pho-X-G/A-W-R-X2-G, where X can be most amino acids, and Pho is hydrophobic) that recognized the tRNA's fully modified anticodon stem and loop domain, hASL(Lys3)UUU. Peptides of this sequence specifically recognized and bound modified htRNA(Lys3)UUU with an affinity 10-fold higher than that of the starting sequence. Thus, this approach provides an effective means of predicting sequences of RNA binding peptides that have better binding properties. Such peptides can be used in cell and molecular biology as well as biochemistry to explore RNA binding proteins and to inhibit those protein functions.
Asunto(s)
Aminoácidos/metabolismo , ARN de Transferencia de Lisina/química , ARN de Transferencia de Lisina/metabolismo , Productos del Gen gag del Virus de la Inmunodeficiencia Humana/metabolismo , Algoritmos , Secuencia de Aminoácidos , Humanos , Modelos Moleculares , Método de Montecarlo , Especificidad por SustratoRESUMEN
Advancements have been made to dissipative particle dynamics (DPD), a robust coarse-grained (CG) simulation method, to study the folded structures of four miniproteins (1L2Y, 1WN8, 1YRF, and 2I9M) in explicit and implicit solvents. In this endeavor, we aim to establish model parametrization and enhance computational efficiency. Unlike traditional CG models that use empirical force parameters, ex-force parameters (r0(ex), a~, δd, δp) of DPD particles constructed for specific research purposes can be obtained from atomistic molecular dynamics simulations. On the other hand, im-force parameters (r0(im), c, σ) can be derived from ex-DPD simulations, according to the underlying thermodynamic theory. Based on a mapping scheme proposed for the modeling of amino acids, all-atom proteins can be converted into a CG model. Both ex-/im-DPDs are then carried out to investigate the folding pathways of the four mini-proteins. Structural analysis of the RMSDs shows that the im-simulated proteins have greater structural similarity to native proteins than the ex-simulated ones. The constructed CG models achieve a resolution of Angstrom (Å), a level normally associated with atomic models. Additionally, speed tests reveal that im-DPD accelerates the simulation process and significantly improves simulation efficiency.
Asunto(s)
Simulación de Dinámica Molecular , Pliegue de Proteína , Proteínas , Solventes , Termodinámica , Solventes/química , Proteínas/química , Conformación ProteicaRESUMEN
Peptide self-assembly into amyloid fibrils provides numerous applications in drug delivery and biomedical engineering applications. We augment our previously-established computational screening technique along with experimental biophysical characterization to discover 7-mer peptides that self-assemble into "parallel ß-sheets", that is, ß-sheets with N-terminus-to-C-terminus ð½-strand vectors oriented in parallel. To accomplish the desired ß-strand organization, we applied the PepAD amino acid sequence design software to the Class-1 cross-ß spine defined by Sawaya et al. This molecular configuration includes two layers of parallel ß-sheets stacked such that N-terminus-to-C-terminus vectors are oriented antiparallel for molecules on adjacent ß-sheets. The first cohort of PepAD identified peptides were examined for their fibrillation behavior in DMD/PRIME20 simulations, and the top performing sequence was selected as a prototype for a subsequent round of sequence refinement. The two rounds of design resulted in a library of eight 7-mer peptides. In DMD/PRIME20 simulations, five of these peptides spontaneously formed fibril-like structures with a predominantly parallel ð½-sheet arrangement, two formed fibril-like structure with <50% in parallel ð½-sheet arrangement and one remained a random coil. Among the eight candidate peptides produced by PepAD and DMD/PRIME20, five were synthesized and purified. All five assembled into amyloid fibrils composed of parallel ß-sheets based on Fourier transform infrared spectroscopy, circular dichroism, electron microscopy, and thioflavin-T fluorescence spectroscopy measurements.
Asunto(s)
Método de Montecarlo , Conformación Proteica en Lámina beta , Nanofibras/química , Péptidos/química , Secuencia de Aminoácidos , Estructura Secundaria de Proteína , Amiloide/química , Modelos Moleculares , Simulación de Dinámica MolecularRESUMEN
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.
RESUMEN
Clostridioides difficile ( C. diff .) is a bacterium that causes severe diarrhea and inflammation of the colon. The pathogenicity of C. diff . infection is derived from two major toxins, toxins A (TcdA) and B (TcdB). Peptide inhibitors that can be delivered to the gut to inactivate these toxins are an attractive therapeutic strategy. In this work, we present a new approach that combines a pep tide b inding d esign algorithm (PepBD), molecular-level simulations, rapid screening of candidate peptides for toxin binding, a primary human cell-based assay, and surface plasmon resonance (SPR) measurements to develop peptide inhibitors that block the glucosyltransferase activity of TcdA by targeting its glucosyltransferase domain (GTD). Using PepBD and explicit-solvent molecular dynamics simulations, we identified seven candidate peptides, SA1-SA7. These peptides were selected for specific TcdA GTD binding through a custom solid-phase peptide screening system, which eliminated the weaker inhibitors SA5-SA7. The efficacies of SA1-SA4 were then tested using a trans-epithelial electrical resistance (TEER) assay on monolayers of the human gut epithelial culture model. One peptide, SA1, was found to block TcdA toxicity in primary-derived human jejunum (small intestinal) and colon (large intestinal) epithelial cells. SA1 bound TcdA with a K D of 56.1 ± 29.8 nM as measured by surface plasmon resonance (SPR). Significance Statement: Infections by Clostridioides difficile , a bacterium that targets the large intestine (colon), impact a significant number of people worldwide. Bacterial colonization is mediated by two exotoxins: toxins A and B. Short peptides that can inhibit the biocatalytic activity of these toxins represent a promising strategy to prevent and treat C. diff . infection. We describe an approach that combines a Peptide B inding D esign (PepBD) algorithm, molecular-level simulations, a rapid screening assay to evaluate peptide:toxin binding, a primary human cell-based assay, and surface plasmon resonance (SPR) measurements to develop peptide inhibitors that block Toxin A in small intestinal and colon epithelial cells. Importantly, our designed peptide, SA1, bound toxin A with nanomolar affinity and blocked toxicity in colon cells.
RESUMEN
Infections by Clostridioides difficile, a bacterium that targets the large intestine (colon), impact a large number of people worldwide. Bacterial colonization is mediated by two exotoxins: toxins A and B. Short peptides that can be delivered to the gut and inhibit the biocatalytic activity of these toxins represent a promising therapeutic strategy to prevent and treat C. diff. infection. We describe an approach that combines a Peptide Binding Design (PepBD) algorithm, molecular-level simulations, a rapid screening assay to evaluate peptide:toxin binding, a primary human cell-based assay, and surface plasmon resonance (SPR) measurements to develop peptide inhibitors that block Toxin A in colon epithelial cells. One peptide, SA1, is found to block TcdA toxicity in primary-derived human colon (large intestinal) epithelial cells. SA1 binds TcdA with a KD of 56.1 ± 29.8 nM as measured by surface plasmon resonance (SPR).
Asunto(s)
Clostridioides difficile , Humanos , Colon , Algoritmos , Biocatálisis , Péptidos/farmacologíaRESUMEN
Clostridium difficile infection is mediated by two major exotoxins: toxins A (TcdA) and B (TcdB). Inhibiting the biocatalytic activities of these toxins with targeted peptide-based drugs can reduce the risk of C. difficile infection. In this work, we used a computational strategy that integrates a peptide binding design (PepBD) algorithm and explicit-solvent atomistic molecular dynamics simulation to determine promising toxin A-targeting peptides that can recognize and bind to the catalytic site of the TcdA glucosyltransferase domain (GTD). Our simulation results revealed that two out of three in silico discovered peptides, viz. the neutralizing peptides A (NPA) and B (NPB), exhibit lower binding free energies when bound to the TcdA GTD than the phage-display discovered peptide, viz. the reference peptide (RP). These peptides may serve as potential inhibitors against C. difficile infection. The efficacy of the peptides RP, NPA, and NPB to neutralize the cytopathic effects of TcdA was tested in vitro in human jejunum cells. Both phage-display peptide RP and in silico peptide NPA were found to exhibit strong toxin-neutralizing properties, thereby preventing the TcdA toxicity. However, the in silico peptide NPB demonstrates a relatively low efficacy against TcdA.
Asunto(s)
Toxinas Bacterianas/antagonistas & inhibidores , Simulación por Computador , Enterotoxinas/antagonistas & inhibidores , Péptidos/farmacología , Toxinas Bacterianas/química , Diseño de Fármacos , Enterotoxinas/química , Modelos Moleculares , Péptidos/química , Unión Proteica , Conformación Proteica , Dominios Proteicos , Reproducibilidad de los ResultadosRESUMEN
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.
Asunto(s)
Fragmentos Fab de Inmunoglobulinas , Inmunoglobulina G , Humanos , Fragmentos Fab de Inmunoglobulinas/química , Ligandos , Péptidos , Distribución TisularRESUMEN
The COVID-19 pandemic has caused significant social and economic disruption across the globe. Cellular entry of SARS-CoV-2 into the human body is mediated via binding of the Receptor Binding Domain (RBD) on the viral Spike protein (SARS-CoV-2 RBD) to Angiotensin-Converting Enzyme 2 (ACE2) expressed on host cells. Molecules that can disrupt ACE2:RBD interactions are attractive therapeutic candidates to prevent virus entry into human cells. A computational strategy that combines our Peptide Binding Design (PepBD) algorithm with atomistic molecular dynamics simulations was used to design new inhibitory peptide candidates via sequence iteration starting with a 23-mer peptide, referred to as SBP1. SBP1 is derived from a region of the ACE2 Peptidase Domain α1 helix that binds to the SARS-CoV-2 RBD of the initial Wuhan-Hu-1 strain. Three peptides demonstrated a solution-phase RBD-binding dissociation constant in the micromolar range during tryptophan fluorescence quenching experiments, one peptide did not bind, and one was insoluble at micromolar concentrations. However, in competitive ELISA assays, none of these peptides could outcompete ACE2 binding to SARS-CoV-2-RBD up to concentrations of 50 µM, similar to the parent SBP1 peptide which also failed to outcompete ACE2:RBD binding. Molecular dynamics simulations suggest that P4 would have a good binding affinity for the RBD domain of Beta-B.1.351, Gamma-P.1, Kappa-B.1.617.1, Delta-B.1.617.2, and Omicron-B.1.1.529 variants, but not the Alpha variant. Consistent with this, P4 bound Kappa-B.1.617.1 and Delta-B.1.617.2 RBD with micromolar affinity in tryptophan fluorescence quenching experiments. Collectively, these data show that while relatively short unstructured peptides can bind to SARS-CoV-2 RBD with moderate affinity, they are incapable of outcompeting the strong interactions between RBD and ACE2.
Asunto(s)
Enzima Convertidora de Angiotensina 2 , COVID-19 , Humanos , SARS-CoV-2 , Glicoproteína de la Espiga del Coronavirus/química , Pandemias , Triptófano/metabolismo , Unión Proteica , Péptidos/metabolismoRESUMEN
Screening amino acid sequence space via experiments to discover peptides that self-assemble into amyloid fibrils is challenging. We have developed a computational peptide assembly design (PepAD) algorithm that enables the discovery of amyloid-forming peptides. Discontinuous molecular dynamics (DMD) simulation with the PRIME20 force field combined with the FoldAmyloid tool is used to examine the fibrilization kinetics of PepAD-generated peptides. PepAD screening of â¼10,000 7-mer peptides resulted in twelve top-scoring peptides with two distinct hydration properties. Our studies revealed that eight of the twelve in silico discovered peptides spontaneously form amyloid fibrils in the DMD simulations and that all eight have at least five residues that the FoldAmyloid tool classifies as being aggregation-prone. Based on these observations, we re-examined the PepAD-generated peptides in the sequence pool returned by PepAD and extracted five sequence patterns as well as associated sequence signatures for the 7-mer amyloid-forming peptides. Experimental results from Fourier transform infrared spectroscopy (FTIR), thioflavin T (ThT) fluorescence, circular dichroism (CD), and transmission electron microscopy (TEM) indicate that all the peptides predicted to assemble in silico assemble into antiparallel ß-sheet nanofibers in a concentration-dependent manner. This is the first attempt to use a computational approach to search for amyloid-forming peptides based on customized settings. Our efforts facilitate the identification of ß-sheet-based self-assembling peptides, and contribute insights towards answering a fundamental scientific question: "What does it take, sequence-wise, for a peptide to self-assemble?".
RESUMEN
Peptides' hierarchical coassembly into nanostructures enables controllable fabrication of multicomponent biomaterials. In this work, we describe a computational and experimental approach to design pairs of charge-complementary peptides that selectively coassemble into ß-sheet nanofibers when mixed together but remain unassembled when isolated separately. The key advance is a peptide coassembly design (PepCAD) algorithm that searches for pairs of coassembling peptides. Six peptide pairs are identified from a pool of ~106 candidates via the PepCAD algorithm and then subjected to DMD/PRIME20 simulations to examine their co-/self-association kinetics. The five pairs that spontaneously aggregate in kinetic simulations selectively coassemble in biophysical experiments, with four forming ß-sheet nanofibers and one forming a stable nonfibrillar aggregate. Solid-state NMR, which is applied to characterize the coassembling pairs, suggests that the in silico peptides exhibit a higher degree of structural order than the previously reported CATCH(+/−) peptides.
RESUMEN
Peptide coassembly, wherein at least two different peptides interact to form multicomponent nanostructures, is an attractive approach for generating functional biomaterials. Current efforts seek to design pairs of peptides, A and B, that form nanostructures (e.g., ß-sheets with ABABA-type ß-strand patterning) while resisting self-assembly (e.g., AAAAA-type or BBBBB-type ß-sheets). To confer coassembly behavior, most existing designs have been based on highly charged variants of known self-assembling peptides; like-charge repulsion limits self-assembly while opposite-charge attraction promotes coassembly. Recent analyses using solid-state NMR and coarse-grained simulations reveal that preconceived notions of structure and molecular organization are not always correct. This perspective highlights recent advances and key challenges to understanding and controlling peptide coassembly.
Asunto(s)
Materiales Biocompatibles , Nanoestructuras , Espectroscopía de Resonancia Magnética , Péptidos , Conformación Proteica en Lámina betaRESUMEN
Wearable sensors for human health, performance, and state monitoring, which have a linear response to the binding of biomarkers found in sweat, saliva, or urine, are of current interest for many applications. A critical part of any device is a biological recognition element (BRE) that is able to bind a biomarker at the surface of a sensor with a high affinity and selectivity to produce a measurable signal response. In this study, we discover and compare 12-mer peptides that bind to neuropeptide Y (NPY), a stress and human health biomarker, using independent and complimentary experimental and computational approaches. The affinities of the NPY-binding peptides discovered by both methods are equivalent and below the micromolar level, which makes them suitable for application in sensors. The in silico design protocol for peptide-based BREs is low cost, highly efficient, and simple, suggesting its utility for discovering peptide binders to a variety of biomarker targets.