Cyclic peptides (CPs) are a promising class of drugs because of their high biological activity and specificity. However, the design of CP remains challenging due to their conformational flexibility and difficulties in designing stable binding conformation. Herein, we present a high-throughput MD screening (HTMDS) process for the iterative design of stable CP binders with a combinatorial CP library composed of canonical and non-canonical amino acids. As a proof of concept, we apply our methods to design CP inhibitors for the bromodomain (BrD) of ATAD2B. 698,800 CP candidates with a total of 25,570 ns MD simulations were performed to study the protein-ligand binding interactions. The binding free energies (ΔGbind) estimated by MM/PBSA approach for eight lead CP designs were found to be low. CP-1st.43 was the best CP candidate with an estimated ΔGbind of -28.48 kcal/mol when compared to the standard inhibitor C-38 which has been experimentally validated and shown to exhibit ΔGbind of -17.11 kcal/mol. The major contribution of binding sites for BrD of ATAD2B involved the hydrogen-bonding anchor within the Aly-binding pocket, salt bridging, and hydrogen-bonding mediated stabilization of the ZA loop and BC loop, and the complementary Van der Waals attraction. Our methods demonstrate encouraging results by yielding conformationally stable and high-potential CP binders that should have potential applicability in future CP drug development.Communicated by Ramaswamy H. Sarma.

Erratum in J Biomol Struct Dyn. 2023 May 27;:1.

In order to know the insights of a unique naturally existing trimodular licheninase from GH16 family, sub-family 21 (RfGH16_21) from Ruminococcus flavefaciens, its structure was modeled to understand its functional relations to reveal information regarding modifying the enzyme for improved properties with enhanced catalytic efficiency. Homology modeling revealed three tandem repeats of β-jelly roll like folds linked by natural linkers. Catalytic pockets and the catalytically important amino acids in each tandem repeat of RfGH16_21 determined by multiple sequence alignment and structure superposition with its homologues indicated that two Glu residues are involved in a retaining-type of catalytic mechanism. Sequential molecular docking revealed maximum binding energy with mixed linked cellotriose showing that cellotriose is the lowest oligomeric hydrolysed product formed by the catalytic action of endo-β-1,3-1,4-glucanase. Molecular dynamic (MD) simulation of RfGH16_21-cellotriose complex confirmed the structural specificity of catalytic residues and increased stability of enzyme in presence of ligand as compared to simulated RfGH16_21 alone. The binding affinity of cellotriose towards the three tandem repeats of RfGH16_21 was also confirmed by calculating total binding Gibbs free energy, i.e. -100.8 ± 2.6 KJ/mol, by using g_mmpbsa tool. The stability of the protein was determined by protein melting analysis that showed Ca2+ and Mg2+ ions imparted structural stability to RfGH16_21. Dynamic light scattering analysis of RfGH16_21 showed monodispersity and hydrodynamic radius of 4.0 nm at 2.0 mg/mL protein concentration, which was comparable with the radius of gyration of 3.2 nm determined by MD simulation showing the protein to be in monomeric form.Communicated by Ramaswamy H. Sarma.