Protein structure generation via folding diffusion.

The ability to computationally generate novel yet physically foldable protein structures could lead to new biological discoveries and new treatments targeting yet incurable diseases. Despite recent advances in protein structure prediction, directly generating diverse, novel protein structures from neural networks remains difficult. In this work, we present a diffusion-based generative model that generates protein backbone structures via a procedure inspired by the natural folding process. We describe a protein backbone structure as a sequence of angles capturing the relative orientation of the constituent backbone atoms, and generate structures by denoising from a random, unfolded state towards a stable folded structure. Not only does this mirror how proteins natively twist into energetically favorable conformations, the inherent shift and rotational invariance of this representation crucially alleviates the need for more complex equivariant networks. We train a denoising diffusion probabilistic model with a simple transformer backbone and demonstrate that our resulting model unconditionally generates highly realistic protein structures with complexity and structural patterns akin to those of naturally-occurring proteins. As a useful resource, we release an open-source codebase and trained models for protein structure diffusion.

Thermodynamic Basis for the Stabilization of Helical Peptoids by Chiral Sidechains.

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.

Origins of Conformational Heterogeneity in Peptoid Helices Formed by Chiral N-1-Phenylethyl Sidechains.

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.

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes.

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.

Cognitive theories of autism based on the interactions between brain functional networks.

Cognitive functions are directly related to interactions between the brain's functional networks. This functional organization changes in the autism spectrum disorder (ASD). However, the heterogeneous nature of autism brings inconsistency in the findings, and specific pattern of changes based on the cognitive theories of ASD still requires to be well-understood. In this study, we hypothesized that the theory of mind (ToM), and the weak central coherence theory must follow an alteration pattern in the network level of functional interactions. The main aim is to understand this pattern by evaluating interactions between all the brain functional networks. Moreover, the association between the significantly altered interactions and cognitive dysfunctions in autism is also investigated. We used resting-state fMRI data of 106 subjects (5-14 years, 46 ASD: five female, 60 HC: 18 female) to define the brain functional networks. Functional networks were calculated by applying four parcellation masks and their interactions were estimated using Pearson's correlation between pairs of them. Subsequently, for each mask, a graph was formed based on the connectome of interactions. Then, the local and global parameters of the graph were calculated. Finally, statistical analysis was performed using a two-sample t-test to highlight the significant differences between autistic and healthy control groups. Our corrected results show significant changes in the interaction of default mode, sensorimotor, visuospatial, visual, and language networks with other functional networks that can support the main cognitive theories of autism. We hope this finding sheds light on a better understanding of the neural underpinning of autism.

Peptide Mimic of the Marine Sponge Protein Silicatein Fabricates Ultrathin Nanosheets of Silicon Dioxide and Titanium Dioxide.

Two-dimensional (2D) materials have attracted attention for potential applications in light harvesting, catalysis, and molecular electronics. Mineral proteins involved in hard tissue biogenesis can produce 2D structures with high fidelity by using sustainable production routes. This study shows that a peptide mimic based on the catalytic triad of the marine sponge protein silicatein catalyzes the formation of nanometer thin and stable sheets of silicon dioxide and titanium dioxide.

Hierarchical Self-Assembly Pathways of Peptoid Helices and Sheets.

Peptoids (N-substituted glycines) are a class of tailorable synthetic peptidomic polymers. Amphiphilic diblock peptoids have been engineered to assemble 2D crystalline lattices with applications in catalysis and molecular separations. Assembly is induced in an organic solvent/water mixture by evaporating the organic phase, but the assembly pathways remain uncharacterized. We conduct all-atom molecular dynamics simulations of Nbrpe6Nc6 as a prototypical amphiphilic diblock peptoid comprising an NH2-capped block of six hydrophobic N-((4-bromophenyl)ethyl)glycine residues conjugated to a polar NH3(CH2)5CO tail. We identify a thermodynamically controlled assembly mechanism by which monomers assemble into disordered aggregates that self-order into 1D chiral helical rods then 2D achiral crystalline sheets. We support our computational predictions with experimental observations of 1D rods using small-angle X-ray scattering, circular dichroism, and atomic force microscopy and 2D crystalline sheets using X-ray diffraction and atomic force microscopy. This work establishes a new understanding of hierarchical peptoid assembly and principles for the design of peptoid-based nanomaterials.

Substitution of distal and active site residues reduces product inhibition of E1 from Acidothermus Cellulolyticus.

Cellulases are largely afflicted by inhibition from their reaction products, especially at high-substrate loading, which represents a major challenge for biomass processing. This challenge was overcome for endoglucanase 1 (E1) from Acidothermus cellulolyticus by identifying a large conformational change involving distal residues upon binding cellobiose. Having introduced alanine substitutions at each of these residues, we identified several mutations that reduced cellobiose inhibition of E1, including W212A, W213A, Q247A, W249A and F250A. One of the mutations (W212A) resulted in a 47-fold decrease in binding affinity of cellobiose as well as a 5-fold increase in the kcat. The mutation further increased E1 activity on Avicel and dilute-acid treated corn stover and enhanced its productivity at high-substrate loadings. These findings were corroborated by funnel metadynamics, which showed that the W212A substitution led to reduced affinity for cellobiose in the +1 and +2 binding sites due to rearrangement of key cellobiose-binding residues.

Ice-nucleating proteins are activated by low temperatures to control the structure of interfacial water.

Ice-nucleation active (INA) bacteria can promote the growth of ice more effectively than any other known material. Using specialized ice-nucleating proteins (INPs), they obtain nutrients from plants by inducing frost damage and, when airborne in the atmosphere, they drive ice nucleation within clouds, which may affect global precipitation patterns. Despite their evident environmental importance, the molecular mechanisms behind INP-induced freezing have remained largely elusive. We investigate the structural basis for the interactions between water and the ice-nucleating protein InaZ from the INA bacterium Pseudomonas syringae. Using vibrational sum-frequency generation (SFG) and two-dimensional infrared spectroscopy, we demonstrate that the ice-active repeats of InaZ adopt a ß-helical structure in solution and at water surfaces. In this configuration, interaction between INPs and water molecules imposes structural ordering on the adjacent water network. The observed order of water increases as the interface is cooled to temperatures close to the melting point of water. Experimental SFG data combined with molecular-dynamics simulations and spectral calculations show that InaZ reorients at lower temperatures. This reorientation can enhance water interactions, and thereby the effectiveness of ice nucleation.

Orientation and Conformation of Proteins at the Air-Water Interface Determined from Integrative Molecular Dynamics Simulations and Sum Frequency Generation Spectroscopy.

Understanding the assembly of proteins at the air-water interface (AWI) informs the formation of protein films, emulsion properties, and protein aggregation. Determination of protein conformation and orientation at an interface is difficult to resolve with a single experimental or simulation technique alone. To date, the interfacial structure of even one of the most widely studied proteins, lysozyme, at the AWI remains unresolved. In this study, molecular dynamics (MD) simulations are used to determine if the protein adopts a side-on, head-on, or axial orientation at the AWI with two different forcefields, GROMOS-53a6 + SPC/E and a99SB-disp + TIP4P-D. Vibrational sum frequency generation (SFG) spectroscopy experiments and spectral SFG calculations validate consistency between the structure determined from MD and experiments. Overall, we show with strong agreement that lysozyme adopts an axial conformation at pH 7. Further, we provide molecular-level insight as to how pH influences the binding domains of lysozyme resulting in side-on adsorption near the isoelectric point of the lysozyme.

MARTINI-Compatible Coarse-Grained Model for the Mesoscale Simulation of Peptoids.

Peptoids (poly-N-substituted glycines) are a class of synthetic polymers that are regioisomers of peptides (poly-C-substituted glycines), in which the point of side-chain connectivity is shifted from the backbone C to the N atom. Peptoids have found diverse applications as peptidomimetic drugs, protein mimetic polymers, surfactants, and catalysts. Computational modeling is valuable in the understanding and design of peptoid-based nanomaterials. In this work, we report the bottom-up parameterization of coarse-grained peptoid force fields based on the MARTINI peptide force field against all-atom peptoid simulation data. Our parameterization pipeline iteratively refits coarse-grained bonded interactions using iterative Boltzmann inversion and nonbonded interactions by matching the potential of mean force for chain extension. We assure good sampling of the amide bond cis/trans isomerizations in the all-atom simulation data using parallel bias metadynamics. We develop coarse-grained models for two representative peptoids-polysarcosine (poly(N-methyl glycine)) and poly(N-((4-bromophenyl)ethyl)glycine)-and show their structural and thermodynamic properties to be in excellent accord with all-atom calculations but up to 25-fold more efficient and compatible with MARTINI force fields. This work establishes a new rigorously parameterized coarse-grained peptoid force field for the understanding and design of peptoid nanomaterials at length and time scales inaccessible to all-atom calculations.

Sequence-Structure-Binding Relationships Reveal Adhesion Behavior of the Car9 Solid-Binding Peptide: An Integrated Experimental and Simulation Study.

Solid-binding peptides (SBPs) recognizing inorganic and synthetic interfaces have enabled a broad range of materials science applications and hold promise as adhesive or morphogenetic control units that can be genetically encoded within desirable or designed protein frameworks. To date, the underlying relationships governing both SBP-surface and SBP-SBP interactions and how they give rise to different adsorption mechanisms remain unclear. Here, we combine protein engineering, surface plasmon resonance characterization, and molecular dynamics (MD) simulations initiated from Rosetta predictions to gain insights on the interplay of amino acid composition, structure, self-association, and adhesion modality in a panel of variants of the Car9 silica-binding peptide (DSARGFKKPGKR) fused to the C-terminus of superfolder green fluorescent protein (sfGFP). Analysis of kinetics, energetics, and MD-predicted structures shows that the high-affinity binding of Car9 to the silanol-rich surface of silica is dominated by electrostatic contributions and a spectrum of several persistent interactions that, along with a high surface population of bound molecules, promote cooperative interactions between neighboring SBPs and higher order structure formation. Transition from cooperative to Langmuir adhesion in sfGFP-Car9 variants occurs in concert with a reduction of stable surface interactions and self-association, as confirmed by atomic force microscopy imaging of proteins exhibiting the two different binding behaviors. We discuss the implications of these results for the de novo design of SBP-surface binding systems.

Impact of Glutamate Carboxylation in the Adsorption of the α-1 Domain of Osteocalcin to Hydroxyapatite and Titania.

One proposed mechanism of implant fouling is attributed to the nonspecific adsorption of non-collagenous bone matrix proteins (NCPs) onto a newly implanted interface. With the goal of capturing the fundamental mechanistic and thermodynamic forces that govern changes in these NCP recognition domains as a function of γ-carboxyglutamic acid (Gla) post-translational modification and surface chemistry, we probe the adsorption process of the most commonly occurring NCP, osteocalcin, onto a mineral and metal oxide surface. Here, we apply two enhanced sampling methods to independently probe the effects of post-translational modification and peptide structure on adsorption. First, well-tempered metadynamics was used to capture the binding of acetyl and N-methylamide capped glutamic acid and Gla single amino acids onto crystalline hydroxyapatite and titania model surfaces at physiological pH. Following this, parallel tempering metadynamics in the well-tempered ensemble (PTMetaD-WTE) was used to study adsorption of the α-1 domain of osteocalcin onto hydroxyapatite and titania. Simulations were performed for the α-1 domain of osteocalcin in both its fully decarboxylated (dOC) and fully carboxylated (OC) form. Our simulations find that increased charge density due to carboxylation results in increased interactions at the interface, and stronger adsorption of the single amino acids to both surfaces. Interestingly, the role of Gla in promoting compact and helical structure in the α-1 domain resulted in disparate binding modes at the two surfaces, which is attributed to differences in interfacial water behavior. Overall, this work provides a benchmark for understanding the mechanisms that drive adsorption of Gla-containing mineralizing proteins onto different surface chemistries.

Ice-binding site of surface-bound type III antifreeze protein partially decoupled from water.

Type III antifreeze proteins (AFP III) have been widely recognized as one class of ice-binding proteins produced by several biological organisms to withstand freezing conditions. Besides their ability to restrict ice growth through their ice-binding site (IBS), AFP III have also been shown to possess a great propensity for hydrophobic surfaces such as the air-water interface. Yet, it is not known whether AFP III adsorb with a specific orientation and how hydrophobic interactions affect the IBS. Molecular insights on the accessibility of the IBS and its interactions with water are important for understanding AFP III action in vivo but also for their application as ice-inhibiting agents for deicing, frozen food storage, as well as for long-term blood and organ cryo-preservation. Here, the orientation of fish AFP III adsorbed at the air-water interface has been studied using a combination of molecular dynamics (MD) simulations and vibrational sum-frequency generation (SFG) spectroscopy together with spectral calculations. The SFG/MD analysis indicated that when AFP III adsorbs at the air-water interface, it mostly retains its native state and orients with a tilt angle of 120° with respect to the surface normal. We found that the IBS is only partially solvated, leaving the pyramidal ice plane binding domain exposed to the vapor phase. These findings suggest that interactions with hydrophobic interfaces (e.g., cell membranes, polymers) could lead to the partial decoupling of the IBS from water and, to some extent, to a loss of AFP III antifreezing activity.

Spin-orbit engineering in transition metal dichalcogenide alloy monolayers.

Binary transition metal dichalcogenide monolayers share common properties such as a direct optical bandgap, spin-orbit splittings of hundreds of meV, light-matter interaction dominated by robust excitons and coupled spin-valley states. Here we demonstrate spin-orbit-engineering in Mo(1-x)WxSe2 alloy monolayers for optoelectronics and applications based on spin- and valley-control. We probe the impact of the tuning of the conduction band spin-orbit spin-splitting on the bright versus dark exciton population. For MoSe2 monolayers, the photoluminescence intensity decreases as a function of temperature by an order of magnitude (4-300 K), whereas for WSe2 we measure surprisingly an order of magnitude increase. The ternary material shows a trend between these two extreme behaviours. We also show a non-linear increase of the valley polarization as a function of tungsten concentration, where 40% tungsten incorporation is sufficient to achieve valley polarization as high as in binary WSe2.