Investigating the role of non-covalent interactions in conformation and assembly of triazine-based sequence-defined polymers

Mechanistic Insights into the Pressure-Induced Polymerization of Aryl/Perfluoroaryl Co-Crystals

Recently discovered diamond nanothreads offer a stiff, sp3-hybridized backbone unachievable in conventional polymer synthesis that is formed through the solid-state pressure-induced polymerization of simple aromatics. This method enables monomeric A-B...

Gaussian-Accelerated Molecular Dynamics with the Weighted Ensemble Method: A Hybrid Method Improves Thermodynamic and Kinetic Sampling

Gaussian-accelerated molecular dynamics (GaMD) is a well-established enhanced sampling method for molecular dynamics simulations that effectively samples the potential energy landscape of the system by adding a boost potential, which smoothens the surface and lowers the energy barriers between states. GaMD is unable to give time-dependent properties such as kinetics directly. On the other hand, the weighted ensemble (WE) method can efficiently sample transitions between states with its many weighted trajectories, which directly yield rates and pathways. However, convergence to equilibrium conditions remains a challenge for the WE method. Hence, we have developed a hybrid method that combines the two methods, wherein GaMD is first used to sample the potential energy landscape of the system and WE is subsequently used to further sample the potential energy landscape and kinetic properties of interest. We show that the hybrid method can sample both thermodynamic and kinetic properties more accurately and quickly compared to using either method alone.

Hierarchical Nanomaterials Assembled from Peptoids and Other Sequence-Defined Synthetic Polymers

In nature, the self-assembly of sequence-specific biopolymers into hierarchical structures plays an essential role in the construction of functional biomaterials. To develop synthetic materials that can mimic and surpass the function of these natural counterparts, various sequence-defined bio- and biomimetic polymers have been developed and exploited as building blocks for hierarchical self-assembly. This review summarizes the recent advances in the molecular self-assembly of hierarchical nanomaterials based on peptoids (or poly-N-substituted glycines) and other sequence-defined synthetic polymers. Modern techniques to monitor the assembly mechanisms and characterize the physicochemical properties of these self-assembly systems are highlighted. In addition, discussions about their potential applications in biomedical sciences and renewable energy are also included. This review aims to highlight essential features of sequence-defined synthetic polymers (e.g., high stability and protein-like high-information content) and how these unique features enable the construction of robust biomimetic functional materials with high programmability and predictability, with an emphasis on peptoids and their self-assembled nanomaterials.

Ranking of Ligand Binding Kinetics Using a Weighted Ensemble Approach and Comparison with a Multiscale Milestoning Approach

To improve lead optimization efforts in finding the right ligand, pharmaceutical industries need to know the ligand's binding kinetics, such as binding and unbinding rate constants, which often correlate with the ligand's efficacy in vivo. To predict binding kinetics efficiently, enhanced sampling methods, such as milestoning and the weighted ensemble (WE) method, have been used in molecular dynamics (MD) simulations of these systems. However, a comparison of these enhanced sampling methods in ranking ligands has not been done. Hence, a WE approach called the concurrent adaptive sampling (CAS) algorithm that uses MD simulations was used to rank seven ligands for β-cyclodextrin, a system in which a multiscale milestoning approach called simulation enabled estimation of kinetic rates (SEEKR) was also used, which uses both MD and Brownian dynamics simulations. Overall, the CAS algorithm can successfully rank ligands using the unbinding rate constant k_off values and binding free energy ΔG values, as SEEKR did, with reduced computational cost that is about the same as SEEKR. We compare the CAS algorithm simulations with different parameters and discuss the impact of parameters in ranking ligands and obtaining rate constant and binding free energy estimates. We also discuss similarities and differences and advantages and disadvantages of SEEKR and the CAS algorithm for future use.