New Monte Carlo Based Technique To Study DNA-Ligand Interactions

Efficient DNA-ligand interaction framework using fuzzy C-means clustering based glowworm swarm optimization (FCMGSO) method

Assessment of DNA and ligand interaction is a great challenge to the medical researchers and drug industries since the accurate mapping of DNA and ligand plays an important role in associating drugs for suitable diseases. The primary objective of this research work is to develop an efficient model for predicting the best DNA and Ligand mapping. In this research work, 500 instances of DNA and drugs used for cancer and non-cancer diseases from the National Centre for Biotechnology Information (NCBI) were considered for analysis. Binding energy is one of the important measures to predict and finalize the best DNA and ligand interaction. Existing methods used for the docking process such as Simulated Annealing (SA), Lamarckian Genetic Algorithm (LGA), Genetic Clustering (GC), Fuzzy C-means clustering (FCM), and Genetic Clustering with Multi swarm Optimization (GCMSO) were applied for all 500 instances. These algorithms failed to produce better binding energy due to a lack of optimization in the existing approaches. Optimization methods play a major role in predicting accurate DNA ligand docking. Hence, this research proposes an efficient architecture using Fuzzy C-Means Clustering with Glowworm Swarm (FCMGSO) optimization method for accurate analysis of the DNA-ligand docking process. Results are proving that the proposed FCMGSO algorithm shows less binding energy than other existing methods in all instances of samples considered from the NCBI dataset.Communicated by Ramaswamy H. Sarma.

Enhanced Monte Carlo Methods for Modeling Proteins Including Computation of Absolute Free Energies of Binding

The generation of a complete ensemble of geometrical configurations is required to obtain reliable estimations of absolute binding free energies by alchemical free energy methods. Molecular dynamics (MD) is the most popular sampling method, but the representation of large biomolecular systems may be incomplete owing to energetic barriers that impede efficient sampling of the configurational space. Monte Carlo (MC) methods can possibly overcome this issue by adapting the attempted movement sizes to facilitate transitions between alternative local-energy minima. In this study, we present an MC statistical mechanics algorithm to explore the protein-ligand conformational space with emphasis on the motions of the protein backbone and side chains. The parameters for each MC move type were optimized to better reproduce conformational distributions of 18 dipeptides and the well-studied T4-lysozyme L99A protein. Next, the performance of the improved MC algorithms was evaluated by computing absolute free energies of binding for L99A lysozyme with benzene and seven analogs. Results for benzene with L99A lysozyme from MD and the optimized MC protocol were found to agree within 0.6 kcal/mol, while a mean unsigned error of 1.2 kcal/mol between MC results and experiment was obtained for the seven benzene analogs. Significant advantages in computation speed are also reported with MC over MD for similar extents of configurational sampling.

Distance-Guided Forward and Backward Chain-Growth Monte Carlo Method for Conformational Sampling and Structural Prediction of Antibody CDR-H3 Loops

Antibodies recognize antigens through the complementary determining regions (CDR) formed by six-loop hypervariable regions crucial for the diversity of antigen specificities. Among the six CDR loops, the H3 loop is the most challenging to predict because of its much higher variation in sequence length and identity, resulting in much larger and complex structural space, compared to the other five loops. We developed a novel method based on a chain-growth sequential Monte Carlo method, called distance-guided sequential chain-growth Monte Carlo for H3 loops (DiSGro-H3). The new method samples protein chains in both forward and backward directions. It can efficiently generate low energy, near-native H3 loop structures using the conformation types predicted from the sequences of H3 loops. DiSGro-H3 performs significantly better than another ab initio method, RosettaAntibody, in both sampling and prediction, while taking less computational time. It performs comparably to template-based methods. As an ab initio method, DiSGro-H3 offers satisfactory accuracy while being able to predict any H3 loops without templates.

Challenges of docking in large, flexible and promiscuous binding sites

After decades of work, the correct determination of the binding mode of a small molecule into a target protein is still a challenging problem, whose difficulty depends on: (i) the sizes of the binding site and the ligand; (ii) the flexibility of both interacting partners, and (iii) the differential solvation of bound and unbound partners. We have evaluated the performance of standard rigid(receptor)/flexible(ligand) docking approaches with respect to last-generation fully flexible docking methods to obtain reasonable poses in a very challenging case: soluble Epoxide Hydrolase (sEH), a flexible protein showing different binding sites. We found that full description of the flexibility of both protein and ligand and accurate description of solvation leads to significant improvement in the ability of docking to reproduce well known binding modes, and at the same time capture the intrinsic binding promiscuity of the protein.