RESUMO
Despite the promise of silk-based devices, the inherent disorder of native silk limits performance. Here, we report highly ordered two-dimensional silk fibroin (SF) films grown epitaxially on van der Waals (vdW) substrates. Using atomic force microscopy, nano-Fourier transform infrared spectroscopy, and molecular dynamics, we show that the films consist of lamellae of SF molecules that exhibit the same secondary structure as the nanocrystallites of native silk. Increasing the SF concentration results in multilayers that grow either by direct assembly of SF molecules into the lamellae or, at high concentrations, along a two-step pathway beginning with a disordered monolayer that then crystallizes. Scanning Kelvin probe measurements show that these films substantially alter the surface potential; thus, they provide a platform for silk-based electronics on vdW solids.
RESUMO
Peptoids are a class of sequence-defined biomimetic polymers with peptide-like backbones and side chains located on backbone nitrogens rather than alpha carbons. These materials demonstrate a strong ability for precise control of single-chain structure, multiunit self-assembly, and macromolecular assembly through careful tuning of sequence due to the diversity of available side chains, although the driving forces behind these assemblies are often not understood. Prior experimental work has shown that linked 15mer peptoids can mimic the protein helical hairpin structure by leveraging the chirality-inducing nature of bulky side chains and hydrophobicity, but there are still gaps in our understanding of the relationship between sequence, stability, and particular secondary or tertiary structure. We present a molecular dynamics (MD) study on the folding behavior of these polymers into hairpins, discussing the differences in structure from sequences with various characteristics in water and acetonitrile, and then compare the handedness preference of common helical motifs between solvents.
Assuntos
Simulação de Dinâmica Molecular , Peptídeos/química , Polímeros/química , Interações Hidrofóbicas e Hidrofílicas , Estrutura Secundária de ProteínaRESUMO
Microbial lipases constitute a class of biocatalysts with the ability to cleave ester linkages of long-chain triglycerides. This property makes them particularly attractive for industrial applications ranging from food processing to pharmaceutical preparation. Among such enzymes, Candida rugosa lipase (CRL) is one of the most frequently used in biotransformation. A notable feature of CRL, among many lipases, is its propensity for interfacial activation: these enzymes exhibit elevated catalytic rates when acting at the interface between aqueous and hydrophobic phases. Notably, this phenomenon can be attributed to the presence of a mobile lid domain, which in its closed state occludes the enzyme active site. To advance our understanding of interfacial activation, we explore the dynamics of CRL rotation at the octane-water interface in this work. To do so, we employ molecular dynamics and umbrella sampling to evaluate the free energy of rotation of the enzyme at the interface. We identify a global minimum in the rotational landscape that coincides with lid opening at the interface. Additionally, we investigate the role of surface residues outside the lid domain as they serve to instigate rotation of the lid toward the aqueous phase. In doing so, we identify a patch of leucine residues which when mutated to glycine impose a barrier to rotation that maintains the enzyme in the inactive (closed lid) state on the order of 1 µs. Importantly, this study presents a novel quantification of the rotational free energy corresponding to CRL lid opening at the octane-water interface. The accompanying mutagenesis study likewise clarifies the role of hydrophobic surface residues in the transition. As such, this work provides valuable insight into the phenomenon of interfacial activation that might open up new avenues for manipulating the microenvironment of industrially relevant lipases, affording enhanced control over the enzyme state.
RESUMO
The development of methods that allow a structural interpretation of linear and nonlinear vibrational spectra is of great importance, both for spectroscopy and for optimizing force field quality. The experimentally measured signals are ensemble averages over all accessible configurations, which complicates spectral calculations. To account for this, we present a recipe for calculating vibrational amide-I spectra of proteins based on metadynamics molecular dynamics simulations. For each frame, a one-exciton Hamiltonian is set up for the backbone amide groups, in which the couplings are estimated with the transition-charge coupling model for nonnearest neighbors, and with a parametrized map of ab initio calculations that give the coupling as a function of the dihedral angles for nearest neighbors. The local-mode frequency variations due to environmental factors such as hydrogen bonds are modeled by exploiting the linear relationship between the amide C-O bond length and the amide-I frequency. The spectra are subsequently calculated while taking into account the equilibrium statistical weights of the frames that are determined using a previously published reweighting procedure. By implementing all these steps in an efficient Fortran code, the spectra can be averaged over very large amounts of structures, thereby extensively covering the phase space of proteins. Using this recipe, the spectral responses of 2.5 million frames of a metadynamics simulation of the miniprotein Trp-cage are averaged to reproduce the experimental temperature-dependent IR spectra very well. The spectral calculations provide new insight into the origin of the various spectral signatures (which are typically challenging to disentangle in the congested amide-I region), and allow for a direct structural interpretation of the experimental spectra and for validation of the molecular dynamics simulations of ensembles.
Assuntos
Simulação de Dinâmica Molecular , Vibração , Peptídeos/química , Amidas/químicaRESUMO
Polyethylene terephthalate (PET) is a type of polymer frequently used in plastic packaging that significantly adds the amount of plastic waste found in landfills. One of the ways to recover valuable raw materials from postconsumer plastic is by depolymerizing PET into its monomeric constituents, which are dimethyl terephthalate (DMT) and ethylene glycol. PET depolymerization is often done in methanolysis with the help of acidic or base catalysts. Tertiary amine is one of the most attractive base catalysts for PET depolymerization in methanolysis since it does not lead to the generation of potentially environmentally harmful waste, unlike metal-based catalysts. However, the mechanism by which tertiary amines catalyze PET depolymerization in methanolysis remains unexplored. Developing a detailed mechanistic understanding of this process is important for improving plastic upcycling since it opens the possibility of employing various cheaper and more environmentally friendly reaction conditions. Using density functional theory and transition state analysis, we show that in the presence of tertiary amine catalysts, methanolysis of PET consists of multiple discrete-step reactions rather than a single concerted step. Furthermore, by comparing our calculations to recent experimental results, we were able to rationalize the DMT yield from the depolymerization process by relating it to charge polarization within tertiary amine catalysts, thus opening a pathway to identify atomic descriptors for future catalyst design.
RESUMO
The 19-residue silaffin-R5 peptide has been widely studied for its ability to precipitate uniform SiO2 particles through mild temperature and pH pathways, in the absence of any organic solvents. There is consensus that post-translational modification (PTM) of side chains has a large impact on the biomineralization process. Thus, it is imperative to understand the precise mechanisms that dictate the formation of SiO2 from R5 peptide, including the effects of PTM on peptide aggregation and peptide-surface adsorption. In this work, we use molecular dynamics (MD) simulations to study the aggregation of R5 dimer with multiple PTMs, with the presence of different ions in solution. Since this system has strong interactions with deep metastable states, we use parallel bias metadynamics with partitioned families to efficiently sample the different states of the system. We find that peptide aggregation is a prerequisite for biomineralization. We observe that the electrostatic interactions are essential in the R5 dimer aggregation; for wild type R5 that only has positively charged residues, phosphate ions HPO4 2- in the solution form a bridge between two peptides and are essential for peptide aggregation.
Assuntos
Simulação de Dinâmica Molecular , Dióxido de Silício/química , Eletricidade Estática , Peptídeos/química , Peptídeos/metabolismo , Processamento de Proteína Pós-Traducional , Fragmentos de Peptídeos , Precursores de ProteínasRESUMO
Inspired by biomineralization, a naturally occurring, protein-facilitated process, solid-binding peptides (SBPs) have gained much attention for their potential to fabricate various shaped nanocrystals and hierarchical nanostructures. The advantage of SBPs over other traditionally used synthetic polymers or short ligands is their tunable interaction with the solid material surface via carefully programmed sequence and being solution-dependent simultaneously. However, designing a sequence with targeted binding affinity or selectivity often involves intensive processes such as phage display, and only a limited number of sequences can be identified. Other computational efforts have also been introduced, but the validation process remains prohibitively expensive once a suitable sequence has been identified. In this paper, we present a new model to rapidly estimate the binding free energy of any given sequence to a solid surface. We show how the overall binding of a polypeptide can be estimated from the free energy contribution of each residue based on the statistics of the thermodynamically stable structure ensemble. We validated our model using five silica-binding peptides of different binding affinities and lengths and showed that the model is accurate and robust across a wider range of chemistries and binding strengths. The computational cost of this method can be as low as 3% of the commonly used enhanced sampling scheme for similar studies and has a great potential to be used in high-throughput algorithms to obtain larger training data sets for machine learning SBP screening.
Assuntos
Nanoestruturas , Peptídeos , Peptídeos/química , Proteínas , Algoritmos , Nanoestruturas/química , Aprendizado de Máquina , Ligação ProteicaRESUMO
Silk fibroin (SF) is a ß-sheet-rich protein that is responsible for the remarkable tensile strength of silk. In addition to its mechanical properties, SF is biocompatible and biodegradable, making it an attractive candidate for use in biotic/abiotic hybrid materials. A pairing of particular interest is the use of SF with graphene-based nanomaterials (GBNs). The properties of this interface drive the formation of well-ordered nanostructures and can improve the electronic properties of the resulting hybrid. It was previously demonstrated that SF can form lamellar nanostructures in the presence of graphite; however, the equilibrium morphology and associated driving interactions are not fully understood. In this study, we characterize these interactions between SF and SF lamellar with graphite using molecular dynamics (MD) simulations and umbrella sampling (US). We find that SF lamellar nanostructures have strong orientational and spatial preferences on graphite that are driven by the hydrophobic effect, destabilizing solvent-protein interactions and stabilizing protein-protein and protein-graphite interactions. Finally, we show how careful consideration of these underlying interactions can be applied to rationally modify the nanostructure morphology.
Assuntos
Fibroínas , Grafite , Nanoestruturas , Fibroínas/química , Grafite/química , Seda/química , Simulação de Dinâmica Molecular , Materiais Biocompatíveis/químicaRESUMO
The misfolding of α-synuclein (αS) into amyloid aggregates is catalyzed by hydrophobic surfaces and associated with severe brain disorders, such as Parkinson's disease. Despite the important role of interfaces, the three-dimensional structure of αS at the interfaces is still not clear. We report interface-specific sum frequency generation (SFG) experiments of monomeric αS binding to the air-water interface, a model system for the important hydrophobic surfaces. We combine the SFG spectra with calculations of theoretical spectra based on molecular dynamics simulations to show that αS, which is an intrinsically disordered protein in solution, folds into a defined, mostly helical secondary structure at the air-water interface. The binding pose resembles an umbrella shape, where the C-terminus protrudes into the water phase, while the N-terminus and the NAC region span the canopy at the interface. In this binding pose, αS is prone to aggregate, which could explain the catalytic effect of hydrophobic interfaces and air bubbles on αS fibrillation.
Assuntos
Doença de Parkinson , alfa-Sinucleína , Humanos , alfa-Sinucleína/química , Água , Doença de Parkinson/metabolismo , Análise Espectral , Simulação de Dinâmica MolecularRESUMO
The degradation of CH3NH3PbI3 (MAPbI3) hybrid organic inorganic perovskite (HOIP) by water has been the major issue hampering its use in commercial perovskites solar cells (PSCs), as MAPbI3 HOIP has been known to easily degrade in the presence of water. Even though there have been numerous studies investigating this phenomenon, there is still no consensus on the mechanisms of the initial stages of dissolution. Here, we attempt to consolidate differing mechanistic interpretations previously reported in the literature through the use of the first-principles constrained ab initio molecular dynamics (AIMD) to study both the energetics and mechanisms that accompany the degradation of MAPbI3 HOIP in liquid water. By comparing the dissolution free energy barrier between surface species of different surficial types, we find that the dominant dissolution mechanisms of surface species vary widely based on the specific surface features. The high sensitivity of the dissolution mechanism to surface features has contributed to the many dissolution mechanisms proposed in the literature. In contrast, the dissolution free energy barriers are mainly determined by the dissolving species rather than the type of surfaces, and the type of surfaces the ions are dissolving from is inconsequential toward the dissolution free energy barrier. However, the presence of surface defects such as vacancy sites is found to significantly lower the dissolution free energy barriers. Based on the estimated dissolution free energy barriers, we propose that the dissolution of MAPbI3 HOIP in liquid water originates from surface defect sites that propagate laterally along the surface layer of the MAPbI3 HOIP crystal.
RESUMO
Thin metal particles with two-dimensional (2D) symmetry are attractive for multiple applications but are difficult to synthesize in a reproducible manner. Although molecules that selectively adsorb to facets have been used to control nanoparticle shape, there is still limited research into the temporal control of growth processes to control these structural outcomes. Moreover, much of the current research into the growth of thin 2D particles lacks mechanistic details. In this work, we study why the substitution of isoleucine for methionine in a gold-binding peptide (Z2, RMRMKMK) results in an increase in gold nanoparticle anisotropy. Nanoplatelet growth in the presence of Z2M246I (RIRIKIK) is characterized using in situ small-angle X-ray scattering (SAXS) and UV-vis spectroscopy. Fitting time-resolved SAXS profiles reveal that 10 nm-thick particles with 2D symmetry are formed within the first few minutes of the reaction. Next, through a combination of electron diffraction and molecular dynamics simulations, we show that substitution of methionine for isoleucine increases the (111) facet selectivity in Z2M246I, and we conclude that this is key to the growth of nanoplatelets. However, the potential application of nanoplatelets formed using Z2M246I is limited due to their uncontrolled lateral growth, aggregation, and rapid sedimentation. Therefore, we use a liquid-handling robot to perform temporally controlled synthesis and dynamic intervention through the addition of Z2 to nanoplatelets grown in the presence of Z2M246I at different times. UV-vis spectroscopy, dynamic light scattering, and electron microscopy show that dynamic intervention results in control over the mean size and stability of plate-like particles. Finally, we use in situ UV-vis spectroscopy to study plate-like particle growth at different times of intervention. Our results demonstrate that both the selectivity and magnitude of binding free energy toward lattices are important for controlling nanoparticle growth pathways.
Assuntos
Nanopartículas Metálicas , Nanoestruturas , Anisotropia , Nanopartículas Metálicas/química , Ouro/química , Espalhamento a Baixo Ângulo , Isoleucina , Difração de Raios X , Nanoestruturas/química , MetioninaRESUMO
This work introduces a three-dimensional (3D) invariant graph-to-string transformer variational autoencoders (VAE) (Vagrant) for generating molecules with accurate density functional theory (DFT)-level properties. Vagrant learns to model the joint probability distribution of a 3D molecular structure and its properties by encoding molecular structures into a 3D-aware latent space. Directed navigation through this latent space implicitly optimizes the 3D structure of a molecule, and the latent embedding can be used to condition a generative transformer to predict the candidate structure as a one-dimensional (1D) sequence. Additionally, we introduce two novel sampling methods that exploit the latent characteristics of a VAE to improve performance. We show that our method outperforms comparable 3D autoregressive and diffusion methods for predicting quantum chemical property values of novel molecules in terms of both sample quality and computational efficiency.
RESUMO
Peptoids are a class of highly customizable biomimetic foldamers that retain properties from both proteins and polymers. It has been shown that peptoids can adopt peptide-like secondary structures through the careful selection of sidechain chemistries, but the underlying conformational landscapes that drive these assemblies at the molecular level remain poorly understood. Given the high flexibility of the peptoid backbone, it is essential that methods applied to study peptoid secondary structure formation possess the requisite sensitivity to discriminate between structurally similar yet energetically distinct microstates. In this work, a generalizable simulation scheme is used to robustly sample the complex folding landscape of various 12mer polypeptoids, resulting in a predictive model that links sidechain chemistry with preferential assembly into one of 12 accessible backbone motifs. Using a variant of the metadynamics sampling method, four peptoid dodecamers are simulated in water: sarcosine, N-(1-phenylmethyl)glycine (Npm), (S)-N-(1-phenylethyl)glycine (Nspe), and (R)-N-(1-phenylethyl)glycine (Nrpe)âto determine the underlying entropic and energetic impacts of hydrophobic and chiral peptoid sidechains on secondary structure formation. Our results indicate that the driving forces to assemble Nrpe and Nspe sequences into polyproline type-I helices in water are found to be enthalpically driven, with small benefits from an entropic gain for isomerization and steric strain due to the presence of the chiral center. The minor entropic gains from bulky chiral sidechains in Nrpe- and Nspe-containing peptoids can be explained through increased configurational entropy in the cis state. However, overall assembly into a helix is found to be overall entropically unfavorable. These results highlight the importance of considering the many various competing interactions in the rational design of peptoid secondary structure building blocks.
Assuntos
Peptoides , Peptoides/química , Glicina/química , Termodinâmica , Estrutura Secundária de Proteína , ÁguaRESUMO
N-substituted glycines (polypeptoids) containing chiral hydrophobic sidechains are known to fold into biomimetic alpha helices. These helix formers often produce conformationally heterogeneous structures and are difficult to characterize at a sub-nanometer resolution. Previously, peptoid N-1-phenylethyl (S)-enantiomer sidechains (Nspe) were inferred from various experiments to form right-handed helices and (R)-enantiomers (Nrpe), left-handed helices. Prior computational work for N(s/r)pe oligomers has struggled to reproduce this trend. Herein, quantum mechanics calculations and molecular dynamics simulations are used to understand the origins of this discrepancy. Results from DFT and molecular mechanics calculations on a variety of Nspe and Nrpe oligomers as a function of chain length are in agreement, showing that Nspe and Nrpe prefer left- and right-handed helices, respectively. Additional metadynamics simulations are used to study Nrpe and Nspe oligomers folding in water. These results show that the free-energy driving forces for assembly into a helical backbone configuration are very small (within â¼kBT). Lastly, we compare DFT calculations for other experimentally characterized peptoid sidechains, N(r/s)sb, N(r/s)tbe, and N(r/s)npe. In this analysis, we show that peptoid sidechains determined to be more robust experimentally (tbe and npe) have helical preferences opposite the trend seen in less robust assemblies formed by N(r/s)pe and N(r/s)sb chemistries. The more robust tbe and nnpe favor the (S)-enantiomer to right-handed and the (R)-enantiomers to left-handed helices.
Assuntos
Peptoides , Peptoides/química , Simulação de Dinâmica Molecular , Estereoisomerismo , Água , Estrutura Secundária de ProteínaRESUMO
Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH2-capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH3(CH2)5CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young's modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. We construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a "Goldilocks" tube radius at which capillary wave fluctuations in the tube wall are minimized.