Machine learning and molecular design of self-assembling -conjugated oligopeptides

Tips and Tricks in the Modeling of Supramolecular Peptide Assemblies

Supramolecular peptide assemblies (SPAs) hold promise as materials for nanotechnology and biomedicine. Although their investigation often entails adapting experimental techniques from their protein counterparts, SPAs are fundamentally distinct from proteins, posing unique challenges for their study. Computational methods have emerged as indispensable tools for gaining deeper insights into SPA structures at the molecular level, surpassing the limitations of experimental techniques, and as screening tools to reduce the experimental search space. However, computational studies have grappled with issues stemming from the absence of standardized procedures and relevant crystal structures. Fundamental disparities between SPAs and protein simulations, such as the absence of experimentally validated initial structures and the importance of the simulation size, number of molecules, and concentration, have compounded these challenges. Understanding the roles of various parameters and the capabilities of different models and simulation setups remains an ongoing endeavor. In this review, we aim to provide readers with guidance on the parameters to consider when conducting SPA simulations, elucidating their potential impact on outcomes and validity.

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.

Versatile in silico modelling of microplastics adsorption capacity in aqueous environment based on molecular descriptor and machine learning

To comprehensively evaluate the hazards of microplastics and their coexisting organic pollutants, the sorption capacity of microplastics is a major issue that is quantified through the microplastic-aqueous sorption coefficient (K_d). Almost all quantitative structure-property relationship (QSPR) models that describe K_d apply only to narrow, relatively homogeneous groups of reactants. Herein, non-hybrid QSPR-based models were developed to predict PE-water (K_PE-w), PE-seawater (K_PE-sw), PVC-water (K_PVC-w) and PP-seawater (K_PP-sw) sorption coefficients at different temperatures, with eight machine learning algorithms. Moreover, novel hybrid intelligent models for predicting K_d more accurately were innovatively developed by applying GA, PSO and AdaBoost algorithms to optimize MLP and ELM models. The results indicated that all three optimization algorithms could improve the robustness and predictability of the standalone MLP and ELM models. In all models trained with K_PE-w, K_PE-sw, K_PVC-w and K_PP-sw data sets, GBDT-1 and XGBoost-1 models, MLP-GA-2 and MLP-PSO-2 models, MLR-3 and MLR-4 models performed better in terms of goodness of fit (R_adj²: 0.907-0.999), robustness (Q_BOOT²: 0.900-0.937) and predictability (R_ext²: 0.889-0.970), respectively. Analyzing the descriptors revealed that temperature, lipophilicity, ionization potential and molecular size were correlated closely with the adsorption capacity of microplastics to organic pollutants. The proposed QSPR models may assist in initial environmental exposure assessments without imposing heavy costs in the early experimental phase.

In Silico Antiprotozoal Evaluation of 1,4-Naphthoquinone Derivatives against Chagas and Leishmaniasis Diseases Using QSAR, Molecular Docking, and ADME Approaches

Chagas and leishmaniasis are two neglected diseases considered as public health problems worldwide, for which there is no effective, low-cost, and low-toxicity treatment for the host. Naphthoquinones are ligands with redox properties involved in oxidative biological processes with a wide variety of activities, including antiparasitic. In this work, in silico methods of quantitative structure–activity relationship (QSAR), molecular docking, and calculation of ADME (absorption, distribution, metabolism, and excretion) properties were used to evaluate naphthoquinone derivatives with unknown antiprotozoal activity. QSAR models were developed for predicting antiparasitic activity against Trypanosoma cruzi, Leishmania amazonensis, and Leishmania infatum, as well as the QSAR model for toxicity activity. Most of the evaluated ligands presented high antiparasitic activity. According to the docking results, the family of triazole derivatives presented the best affinity with the different macromolecular targets. The ADME results showed that most of the evaluated compounds present adequate conditions to be administered orally. Naphthoquinone derivatives show good biological activity results, depending on the substituents attached to the quinone ring, and perhaps the potential to be converted into drugs or starting molecules.

Computationally Guided Tuning of Peptide-Conjugated Perylene Diimide Self-Assembly

Peptide-π-conjugated materials are important for biointerfacing charge-transporting applications due to their aqueous compatibility and formation of long-range π-electron networks. Perylene diimides (PDIs), well-established charge-transporting π systems, can self-assemble in aqueous solutions when conjugated with amino acids. In this work, we leveraged computational guidance from our previous work to access two different self-assembled architectures from PDI-amino acid conjugates. Furthermore, we expanded the design rule to other sequences to learn that the closest amino acids to the π core have a significant effect on the photophysical properties of the resulting assemblies. By simply altering glycine to alanine at the closest residue position, we observed significantly different electronic properties as revealed through UV-vis, photoluminescence, and circular dichroism spectroscopies. Accompanying molecular dynamics simulations revealed two distinct types of self-assembled architectures: cofacial structures when the smaller glycine residue is at the closest residue position to the π core versus rotationally shifted structures when glycine is substituted for the larger alanine. This study illustrates the use of tandem computations and experiments to unearth and understand new design rules for supramolecular materials and exposes a modest amino acid substitution as a means to predictably modulate the supramolecular organization and engineer the photophysical properties of π-conjugated peptidic materials.

Robust classification-based molecular modelling of diverse chemical entities as potential SARS-CoV-2 3CLpro inhibitors: theoretical justification in light of experimental evidences

COVID-19 is the most unanticipated incidence of 2020 affecting the human population worldwide. Currently, it is utmost important to produce novel small molecule anti-SARS-CoV-2 drugs urgently that can save human lives globally. Based on the earlier SARS-CoV and MERS-CoV infection along with the general characters of coronaviral replication, a number of drug molecules have been proposed. However, one of the major limitations is the lack of experimental observations with different drug molecules. In this article, 70 diverse chemicals having experimental SARS-CoV-2 3CL^proinhibitory activity were accounted for robust classification-based QSAR analysis statistically validated with 4 different methodologies to recognize the crucial structural features responsible for imparting the activity. Results obtained from all these methodologies supported and validated each other. Important observations obtained from these analyses were also justified with the ligand-bound crystal structure of SARS-CoV-2 3CL^pro enzyme. Our results suggest that molecules should contain a 2-oxopyrrolidine scaffold as well as a methylene (hydroxy) sulphonic acid warhead in proper orientation to achieve higher inhibitory potency against SARS-CoV-2 3CL^pro. Outcomes of our study may be able to design and discover highly effective SARS-CoV-2 3CL^pro inhibitors as potential anticoronaviral therapy to crusade against COVID-19.

Membrane Binding of Antimicrobial Peptides Is Modulated by Lipid Charge Modification

Peptide interactions with lipid bilayers play a key role in a range of biological processes and depend on electrostatic interactions between charged amino acids and lipid headgroups. Antimicrobial peptides (AMPs) initiate the killing of bacteria by binding to and destabilizing their membranes. The multiple peptide resistance factor (MprF) provides a defense mechanism for bacteria against a broad range of AMPs. MprF reduces the negative charge of bacterial membranes through enzymatic conversion of the anionic lipid phosphatidyl glycerol (PG) to either zwitterionic alanyl-phosphatidyl glycerol (Ala-PG) or cationic lysyl-phosphatidyl glycerol (Lys-PG). The resulting change in the membrane charge is suggested to reduce the binding of AMPs to membranes, thus impeding downstream AMP activity. Using coarse-grained molecular dynamics to investigate the effects of these modified lipids on AMP binding to model membranes, we show that AMPs have substantially reduced affinity for model membranes containing Ala-PG or Lys-PG. More than 5000 simulations in total are used to define the relationship between lipid bilayer composition, peptide sequence (using five different membrane-active peptides), and peptide binding to membranes. The degree of interaction of a peptide with a membrane correlates with the membrane surface charge density. Free energy profile (potential of mean force) calculations reveal that the lipid modifications due to MprF alter the energy barrier to peptide helix penetration of the bilayer. These results will offer a guide to the design of novel peptides, which addresses the issue of resistance via MprF-mediated membrane modification.

Improving the Accuracy of Nearest-Neighbor Classification Using Principled Construction and Stochastic Sampling of Training-Set Centroids

A conceptually simple way to classify images is to directly compare test-set data and training-set data. The accuracy of this approach is limited by the method of comparison used, and by the extent to which the training-set data cover configuration space. Here we show that this coverage can be substantially increased using coarse-graining (replacing groups of images by their centroids) and stochastic sampling (using distinct sets of centroids in combination). We use the MNIST and Fashion-MNIST data sets to show that a principled coarse-graining algorithm can convert training images into fewer image centroids without loss of accuracy of classification of test-set images by nearest-neighbor classification. Distinct batches of centroids can be used in combination as a means of stochastically sampling configuration space, and can classify test-set data more accurately than can the unaltered training set. On the MNIST and Fashion-MNIST data sets this approach converts nearest-neighbor classification from a mid-ranking- to an upper-ranking member of the set of classical machine-learning techniques.

In search of a novel chassis material for synthetic cells: emergence of synthetic peptide compartment

Giant lipid vesicles have been used extensively as a synthetic cell model to recapitulate various life-like processes, including in vitro protein synthesis, DNA replication, and cytoskeleton organization. Cell-sized lipid vesicles are mechanically fragile in nature and prone to rupture due to osmotic stress, which limits their usability. Recently, peptide vesicles have been introduced as a synthetic cell model that would potentially overcome the aforementioned limitations. Peptide vesicles are robust, reasonably more stable than lipid vesicles and can withstand harsh conditions including pH, thermal, and osmotic variations. This mini-review summarizes the current state-of-the-art in the design, engineering, and realization of peptide-based chassis materials, including both experimental and computational work. We present an outlook for simulation-aided and data-driven design and experimental realization of engineered and multifunctional synthetic cells.

Theoretical design of novel antimalarial agents against P. falciparum strain, Dd2 through the QSAR modeling of synthesized 2'-substituted triclosan derivatives

In an attempt to design compounds with higher antimalarial activities, quantitative structure-activity relationship (QSAR) technique was utilized in the development of a molecular model for some synthesized 2′-substituted triclosan derivatives through a hybrid of the GA-MLR method. The model was found to have excellent statistical parameters (R² = 0.8919, R²_Adj = 0.8728, LOF = 0.2563). The descriptors mean effect (MF) revealed BCUTw-1l, which increases with an increase in molecular weight, to be the most contributive to the antimalarial activity. Consequently, compound 3, with the highest activities (pEC₅₀ = 6.9586) was deployed as the design template. The molecular weight of the template was increasing through substitutions of its atoms at several positions with heavier atoms/groups to increases the descriptor (BCUTw-1l) value. Twelves (12) theoretical derivatives of the template were designed where six of the designed derivatives have better activity than the design template. The most active designed compound, 3L was found to have the highest antimalarial activity (pEC₅₀ = 7.930) than that of the design template.

Evolutionary reinforcement learning of dynamical large deviations

We show how to bound and calculate the likelihood of dynamical large deviations using evolutionary reinforcement learning. An agent, a stochastic model, propagates a continuous-time Monte Carlo trajectory and receives a reward conditioned upon the values of certain path-extensive quantities. Evolution produces progressively fitter agents, potentially allowing the calculation of a piece of a large-deviation rate function for a particular model and path-extensive quantity. For models with small state spaces, the evolutionary process acts directly on rates, and for models with large state spaces, the process acts on the weights of a neural network that parameterizes the model's rates. This approach shows how path-extensive physics problems can be considered within a framework widely used in machine learning.

Learning to grow: Control of material self-assembly using evolutionary reinforcement learning

We show that neural networks trained by evolutionary reinforcement learning can enact efficient molecular self-assembly protocols. Presented with molecular simulation trajectories, networks learn to change temperature and chemical potential in order to promote the assembly of desired structures or choose between competing polymorphs. In the first case, networks reproduce in a qualitative sense the results of previously known protocols, but faster and with higher fidelity; in the second case they identify strategies previously unknown, from which we can extract physical insight. Networks that take as input the elapsed time of the simulation or microscopic information from the system are both effective, the latter more so. The evolutionary scheme we have used is simple to implement and can be applied to a broad range of examples of experimental self-assembly, whether or not one can monitor the experiment as it proceeds. Our results have been achieved with no human input beyond the specification of which order parameter to promote, pointing the way to the design of synthesis protocols by artificial intelligence.

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.

Revealing the Sequence-Structure-Electronic Property Relation of Self-Assembling π-Conjugated Oligopeptides by Molecular and Quantum Mechanical Modeling

Self-assembled nanoaggregates of π-conjugated synthetic peptides present a biocompatible and highly tunable alternative to silicon-based optical and electronic materials. Understanding the relationship between structural morphology and electronic properties of these assemblies is critical for understanding and controlling their mechanical, optical, and electronic responses. In this work, we combine all-atom classical molecular simulations with quantum mechanical electronic structure calculations to ascertain the sequence-structure-electronic property relationship within a family of Asp-X-X-quaterthiophene-X-X-Asp (DXX-OT4-XXD) oligopeptides in which X is one of the five amino acids {Ala, Phe, Gly, Ile, Val} ({A, F, G, I, V}). Molecular dynamics simulations reveal that smaller amino acid substituents (A, G) favor linear stacking within a peptide dimer, whereas larger groups (F, I, V) induce larger twist angles between the peptides. Density functional theory calculations on the dimer show the absorption spectrum to be dominated by transitions between carbon and sulfur p orbitals. Although the absorption spectrum is largely insensitive to the relative twist angle, the highest occupied molecular orbital strongly localizes onto one molecule within the dimer at large twist angles, impeding the efficiency of transport between molecules. Our results provide a fundamental understanding of the relation between peptide orientation and electronic structure and offer design precepts for rational engineering of these systems.