Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and prime/MM-GBSA calculation

Performance of the MM/GBSA scoring using a binding site hydrogen bond network-based frame selection: the protein kinase case

Conformational clustering using hydrogen bond network analysis improved the MM/GBSA scoring for some protein-kinase–ligand systems used as case studies.

Deciphering the glycan preference of bacterial lectins by glycan array and molecular docking with validation by microcalorimetry and crystallography

Recent advances in glycobiology revealed the essential role of lectins for deciphering the glycocode by specific recognition of carbohydrates. Integrated multiscale approaches are needed for characterizing lectin specificity: combining on one hand high-throughput analysis by glycan array experiments and systematic molecular docking of oligosaccharide libraries and on the other hand detailed analysis of the lectin/oligosaccharide interaction by x-ray crystallography, microcalorimetry and free energy calculations. The lectins LecB from Pseudomonas aeruginosa and BambL from Burkholderia ambifaria are part of the virulence factors used by the pathogenic bacteria to invade the targeted host. These two lectins are not related but both recognize fucosylated oligosaccharides such as the histo-blood group oligosaccharides of the ABH(O) and Lewis epitopes. The specificities were characterized using semi-quantitative data from glycan array and analyzed by molecular docking with the Glide software. Reliable prediction of protein/oligosaccharide structures could be obtained as validated by existing crystal structures of complexes. Additionally, the crystal structure of BambL/Lewis x was determined at 1.6 Å resolution, which confirms that Lewis x has to adopt a high-energy conformation so as to bind to this lectin. Free energies of binding were calculated using a procedure combining the Glide docking protocol followed by free energy rescoring with the Prime/Molecular Mechanics Generalized Born Surface Area (MM-GBSA) method. The calculated data were in reasonable agreement with experimental free energies of binding obtained by titration microcalorimetry. The established predictive protocol is proposed to rationalize large sets of data such as glycan arrays and to help in lead discovery projects based on such high throughput technology.

In silico design: extended molecular dynamic simulations of a new series of dually acting inhibitors against EGFR and HER2

Based on the hit structures that have been identified in our previous studies against EGFR and HER2, new potential inhibitors that share the same scaffold of the hit structures are designed and screened in silico. Insights into understanding the potential inhibitory effect of the new inhibitors against both EGFR and HER2 receptors is obtained using extended molecular dynamics (MD) simulations and different scoring techniques. The binding mechanisms and dynamics are detailed with respect to two approved inhibitors against EGFR (lapatinib) and HER2 (SYR127063). The best scoring inhibitor (T9) is chosen for additional in silico investigation against both the wild-type and T790M mutant strain of EGFR and the wild-type HER2. The results reveal that certain substitution patterns increase the stability and assure stronger binding and higher H-bond occupancy of the conserved water molecule that is commonly observed with kinase crystal structures. Furthermore, the new inhibitor (T9) forms stable interactions with the mutant strain as a direct consequence of the enhanced ability to form additional hydrogen bonding interactions with binding site residues.

Genome-based approaches to develop epitope-driven subunit vaccines against pathogens of infective endocarditis

Infective endocarditis (IE) has emerged as a public health problem due to changes in the etiologic spectrum and due to involvement of resistant bacterial strains with increased virulence. Developing potent vaccine is an important strategy to tackle IE. Complete genome sequences of eight selected pathogens of IE paved the way to design common T-cell driven subunit vaccines. Comparative genomics and subtractive genomic analysis were applied to identify adinosine tri phosphate (ATP)-binding cassette (ABC) transporter ATP-binding protein from Streptococcus mitis (reference organism) as common vaccine target. Reverse vaccinology technique was implemented using computational tools such as ProPred, SYFPEITHI, and Immune epitope database. Twenty-one T-cell epitopes were predicted from ABC transporter ATP-binding protein. Multiple sequence alignment of ABC transporter ATP-binding protein from eight selected IE pathogens was performed to identify six conserved T-cell epitopes. The six selected T-cell epitopes were further evaluated at structure level for HLA-DRB binding through homology modeling and molecular docking analysis using Maestro v9.2. The proposed six T-cell epitopes showed better binding affinity with the selected HLA-DRB alleles. Subsequently, the docking complexes of T-cell epitope and HLA-DRBs were ranked based on XP Gscore. The T-cell epitope (208-LNYITPDVV-216)-HLA-DRB1(∗)0101 (1T5 W) complex having the best XP Gscore (-13.25 kcal/mol) was assessed for conformational stability and interaction stability through molecular dynamic simulation for 10 ns using Desmond v3.2. The simulation results revealed that the HLA-DRB-epitope complex was stable throughout the simulation time. Thus, the epitope would be ideal candidate for T-cell driven subunit vaccine design against infective endocarditis.

Quantum Mechanics Approaches to Drug Research in the Era of Structural Chemogenomics

The rapid growth of the available crystallographic information about proteins and binding pockets creates remarkable opportunities for enriching the drug research pipelines with computational prediction of novel protein-ligand interactions. While ab initio quantum mechanical approaches are known to provide unprecedented accuracy in structure-based binding energy calculations, they are limited to only small systems of dozens of atoms. In the structural chemogenomics era, it is critical that new approaches are developed that enable application of QM methodologies to non-covalent interactions in systems as large as protein-ligand complexes and conformational ensembles. This perspective highlights recent advances towards bridging the gap between high accuracy and high volume computations in drug research.

Discovery of new HER2/EGFR dual kinase inhibitors based on the anilinoquinazoline scaffold as potential anti-cancer agents

Herein, we designed and synthesized certain anilinoquinazoline derivatives bearing bulky arylpyridinyl, arylpropenoyl and arylpyrazolyl moieties at the 4' position of the anilinoquinazoline, as potential dual HER2/EGFR kinase inhibitors. A detailed molecular modeling study was performed by docking the synthesized compounds in the active site of the epidermal growth factor receptor (EGFR). The synthesized compounds were further tested for their inhibitory activity on EGFR and HER2 tyrosine kinases. The aryl 2-imino-1,2-dihydropyridine derivatives 5d and 5e displayed the most potent inhibitory activity on EGFR with IC50 equal to 2.09 and 1.94 μM, respectively, and with IC50 equal to 3.98 and 1.04 μM on HER2, respectively. Furthermore, the anti-proliferative activity of these most active compounds on MDA-MB-231 breast cancer cell lines, known to overexpress EGFR, showed an IC50 range of 2.4 and 2.5 μM, respectively.

Elucidation of allosteric inhibition mechanism of 2-Cys human peroxiredoxin by molecular modeling

We used molecular dynamics (MD) simulations and protein docking to elucidate the mechanism of allosteric inhibition of the human form of peroxiredoxin (Prx), 2-Cys proliferation associated gene (PAG). Beginning by using the rat form of Prx, 2-Cys heme-binding protein as a template, we used homology modeling to find the structure of human 2-Cys PAG, which is in dimeric form. Molecular dynamics simulations showed that the structure of the reduced form of the 2-Cys PAG dimer fluctuates as the two monomers drift away and approach each other. We then used SiteMap to search for binding sites on the surface of this dimer. A binding site between the two monomers was found, and virtual screening with docking was performed to identify a ligand binding to this site. Subsequent MD simulation revealed that with this ligand in the binding site, the dimer structure of 2-Cys PAG becomes stabilized such that two cysteine residues from two monomers, which are partners of a disulfide bond of the oxidized form, remain separated. This mechanism can be used as an allosteric inhibition of Prx as a hydrogen peroxide reducer, the role of which has been studied as an anticancer drug target.