Patchy Particle Model of the Hierarchical Self-Assembly of π-Conjugated Optoelectronic Peptides

Computational Design of Peptide Assemblies

With the ongoing development of peptide self-assembling materials, there is growing interest in exploring novel functional peptide sequences. From short peptides to long polypeptides, as the functionality increases, the sequence space is also expanding exponentially. Consequently, attempting to explore all functional sequences comprehensively through experience and experiments alone has become impractical. By utilizing computational methods, especially artificial intelligence enhanced molecular dynamics (MD) simulation and de novo peptide design, there has been a significant expansion in the exploration of sequence space. Through these methods, a variety of supramolecular functional materials, including fibers, two-dimensional arrays, nanocages, etc., have been designed by meticulously controlling the inter- and intramolecular interactions. In this review, we first provide a brief overview of the current main computational methods and then focus on the computational design methods for various self-assembled peptide materials. Additionally, we introduce some representative protein self-assemblies to offer guidance for the design of self-assembling peptides.

Evolution of π-Peptide Self-Assembly: From Understanding to Prediction and Control

Supramolecular materials derived from the self-assembly of engineered molecules continue to garner tremendous scientific and technological interest. Recent innovations include the realization of nano- and mesoscale particles (0D), rods and fibrils (1D), sheets (2D), and even extended lattices (3D). Our research groups have focused attention over the past 15 years on one particular class of supramolecular materials derived from oligopeptides with embedded π-electron units, where the oligopeptides can be viewed as substituents or side chains to direct the assembly of the central π-electron cores. Upon assembly, the π-systems are driven into close cofacial architectures that facilitate a variety of energy migration processes within the nanomaterial volume, including exciton transport, voltage transmission, and photoinduced electron transfer. Like many practitioners of supramolecular materials science, many of our initial molecular designs were designed with substantial inspiration from biologically occurring self-assembly coupled with input from chemical intuition and molecular modeling and simulation. In this feature article, we summarize our current understanding of the π-peptide self-assembly process as documented through our body of publications in this area. We address fundamental spectroscopic and computational tools used to extract information regarding the internal structures and energetics of the π-peptide assemblies, and we address the current state of the art in terms of recent applications of data science tools in conjunction with high-throughput computational screening and experimental assays to guide the efficient traversal of the π-peptide molecular design space. The abstract image details our integrated program of chemical synthesis, spectroscopic and functional characterization, multiscale simulation, and machine learning which has advanced the understanding and control of the assembly of synthetic π-conjugated peptides into supramolecular nanostructures with energy and biomedical applications.

Coarse-grained modelling to predict the packing of porous organic cages

How molecules pack has vital ramifications for their applications as functional molecular materials. Small changes in a molecule's functionality can lead to large, non-intuitive, changes in their global solid-state packing, resulting in difficulty in targeted design. Predicting the crystal structure of organic molecules from only their molecular structure is a well-known problem plaguing crystal engineering. Although relevant to the properties of many organic molecules, the packing behaviour of modular porous materials, such as porous organic cages (POCs), greatly impacts the properties of the material. We present a novel way of predicting the solid-state phase behaviour of POCs by using a simplistic model containing the dominant degrees of freedom driving crystalline phase formation. We employ coarse-grained simulations to systematically study how chemical functionality of pseudo-octahedral cages can be used to manipulate the solid-state phase formation of POCs. Our results support those of experimentally reported structures, showing that for cages which pack via their windows forming a porous network, only one phase is formed, whereas when cages pack via their windows and arenes, the phase behaviour is more complex. While presenting a lower computational cost route for predicting molecular crystal packing, coarse-grained models also allow for the development of design rules which we start to formulate through our results.

Computer simulations of self-assembly of anisotropic colloids

Computer simulations have played a significant role in understanding the physics of colloidal self-assembly, interpreting experimental observations, and predicting novel mesoscopic and crystalline structures. Recent advances in computer simulations of colloidal self-assembly driven by anisotropic or orientation-dependent inter-particle interactions are highlighted in this review. These interactions are broadly classified into two classes: entropic and enthalpic interactions. They mainly arise due to shape anisotropy, surface heterogeneity, compositional heterogeneity, external field, interfaces, and confinements. Key challenges and opportunities in the field are discussed.

The Martini Model in Materials Science

The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here, the main applications to date of the Martini model in materials science are highlighted, and a perspective for the future developments in this field is given, particularly in light of recent developments such as the new version of the model, Martini 3.

Molecular Modeling and Simulations of Peptide–Polymer Conjugates

Peptide–polymer conjugates are a class of soft materials composed of covalently linked blocks of protein/polypeptides and synthetic/natural polymers. These materials are practically useful in biological applications, such as drug delivery, DNA/gene delivery, and antimicrobial coatings, as well as nonbiological applications, such as electronics, separations, optics, and sensing. Given their broad applicability, there is motivation to understand the molecular and macroscale structure, dynamics, and thermodynamic behavior exhibited by such materials. We focus on the past and ongoing molecular simulation studies aimed at obtaining such fundamental understanding and predicting molecular design rules for the target function. We describe briefly the experimental work in this field that validates or motivates these computational studies. We also describe the various models (e.g., atomistic, coarse-grained, or hybrid) and simulation methods (e.g., stochastic versus deterministic, enhanced sampling) that have been used and the types of questions that have been answered using these computational approaches.

100th Anniversary of Macromolecular Science Viewpoint: Modeling and Simulation of Macromolecules with Hydrogen Bonds: Challenges, Successes, and Opportunities

Macromolecular materials with directional interactions such as hydrogen bonds exhibit numerous attractive features in terms of structure, thermodynamics, and dynamics. Besides enabling precise tuning of desirable geometries in the assembled state (e.g., programmable coordination numbers depending on the valency of the directional interaction), mixing in a blend/composite through stabilization via hydrogen bonds between the various components, hydrogen bonds can also impart responsiveness to external stimuli (e.g., temperature, pH). In biomacromolecules (e.g., proteins, DNA, polysaccharides), hydrogen bonds play a key role in stabilizing secondary and tertiary structures, which in turn define the function of these macromolecules. In this Viewpoint, I present the challenges, successes, and opportunities for molecular modeling and simulations to conduct fundamental and application-focused research on macromolecular materials with hydrogen bonding interactions. The past successes and limitations of atomistic simulations are discussed first, followed by highlights from recent developments in coarse-grained modeling and their use in studies of (synthetic and biologically relevant) macromolecular materials. Model development focused on polynucleotides (e.g., DNA, RNA, etc.), polypeptides, polysaccharides, and synthetic polymers at experimentally relevant conditions are highlighted. This viewpoint ends with potential future directions for macromolecular modeling and simulations with other types of directional interactions beyond hydrogen bonding.

Discovery of Self-Assembling π-Conjugated Peptides by Active Learning-Directed Coarse-Grained Molecular Simulation

Electronically active organic molecules have demonstrated great promise as novel soft materials for energy harvesting and transport. Self-assembled nanoaggregates formed from π-conjugated oligopeptides composed of an aromatic core flanked by oligopeptide wings offer emergent optoelectronic properties within a water-soluble and biocompatible substrate. Nanoaggregate properties can be controlled by tuning core chemistry and peptide composition, but the sequence-structure-function relations remain poorly characterized. In this work, we employ coarse-grained molecular dynamics simulations within an active learning protocol employing deep representational learning and Bayesian optimization to efficiently identify molecules capable of assembling pseudo-1D nanoaggregates with good stacking of the electronically active π-cores. We consider the DXXX-OPV3-XXXD oligopeptide family, where D is an Asp residue and OPV3 is an oligophenylenevinylene oligomer (1,4-distyrylbenzene), to identify the top performing XXX tripeptides within all 20³ = 8000 possible sequences. By direct simulation of only 2.3% of this space, we identify molecules predicted to exhibit superior assembly relative to those reported in prior work. Spectral clustering of the top candidates reveals new design rules governing assembly. This work establishes new understanding of DXXX-OPV3-XXXD assembly, identifies promising new candidates for experimental testing, and presents a computational design platform that can be generically extended to other peptide-based and peptide-like systems.