Ensemble docking can be a successful virtual screening technique that addresses the innate conformational heterogeneity of macromolecular drug targets. Yet, lacking a method to identify a subset of conformational states that effectively segregates active and inactive small molecules, ensemble docking may result in the recommendation of a large number of false positives. Here, three knowledge-based methods that construct structural ensembles for virtual screening are presented. Each method selects ensembles by optimizing an objective function calculated using the receiver operating characteristic (ROC) curve: either the area under the ROC curve (AUC) or a ROC enrichment factor (EF). As the number of receptor conformations, N, becomes large, the methods differ in their asymptotic scaling. Given a set of small molecules with known activities and a collection of target conformations, the most resource intense method is guaranteed to find the optimal ensemble but scales as O(2(N)). A recursive approximation to the optimal solution scales as O(N(2)), and a more severe approximation leads to a faster method that scales linearly, O(N). The techniques are generally applicable to any system, and we demonstrate their effectiveness on the androgen nuclear hormone receptor (AR), cyclin-dependent kinase 2 (CDK2), and the peroxisome proliferator-activated receptor δ (PPAR-δ) drug targets. Conformations that consisted of a crystal structure and molecular dynamics simulation cluster centroids were used to form AR and CDK2 ensembles. Multiple available crystal structures were used to form PPAR-δ ensembles. For each target, we show that the three methods perform similarly to one another on both the training and test sets.
The allosteric modulation of the homodimeric H10-03-6 protein to glycan ligands L1 and L2, and the STAB19 protein to glycan ligands L3 and L4, respectively, has been studied by molecular dynamics simulations and free energy calculations. The results revealed that the STAB19 protein has a significantly higher affinity for L3 (-11.38 ± 2.32 kcal/mol) than that for L4 (-5.51 ± 1.92 kcal/mol). However, the combination of the H10-03-6 protein with glycan L2 (1.23 ± 6.19 kcal/mol) is energetically unfavorable compared with that of L1 (-13.96 ± 0.35 kcal/mol). Further, the binding of glycan ligands L3 and L4 to STAB19 would result in the significant closure of the two CH2 domains of the STAB19 conformation with the decrease of the centroid distances between the two CH2 domains compared with the H10-03-6/L1/L2 complex. The CH2 domain closure of STAB19 relates directly to the formation of new hydrogen bonds and hydrophobic interactions between the residues Ser239, Val240, Asp265, Glu293, Asn297, Thr299, Ser337, Asp376, Thr393, Pro395, and Pro396 in STAB19 and glycan ligands L3 and L4, which suggests that these key residues would contribute to the specific regulation of STAB19 to L3 and L4. In addition, the distance analysis revealed that the EF loop in the H10-03-6/L1/L2 model presents a high flexibility and partial disorder compared with the stabilized STAB19/L3/L4 complex. These results will be helpful in understanding the specific regulation through the asymmetric structural characteristics in the CH2 and CH3 domains of the H10-03-6 and STAB19 proteins.
The continuing rise in tuberculosis incidence and the problem of drug resistance strains have prompted the research on new drug candidates and the mechanism of drug resistance. Molecular docking and molecular dynamics simulation (MD) were performed to study the binding of isoniazid onto the active site of Mycobacterium tuberculosis enoyl-acyl carrier protein reductase (InhA) in an attempt to address the mycobacterial resistance against isoniazid. Results show that isonicotinic acyl-NADH (INADH) has an extremely high binding affinity toward the wild type InhA by forming stronger interactions compared to the parent drug (isoniazid) (INH). Due to the increase of hydrophobicity and reduction in the side chain's volume of A94 of mutant type InhA, both INADH and the mutated protein become more mobile. Due to this reason, the molecular interactions of INADH with mutant type are weaker than that observed with the wild type. However, the reduced interaction caused by the fluctuation of INADH and the mutant protein only inflected minor resistance in the mutant strain as inferred from free energy calculation. MD results also showed there exists a water-mediated hydrogen bond between INADH and InhA. However, the bridged water molecule is only present in the INADH-wild type complex, reflecting the putative role of the water molecule in the binding of INADH to the wild type protein. The results support the assumption that the conversion of prodrug isoniazid into its active form INADH is mediated by KatG as a necessary step prior to target binding on InhA. Our findings also contribute to a better understanding of INH resistance in mutant type.
T1 lipase is isolated from the palm Geobacillus zalihae strain T1 in Malaysia, functioning as a secreted protein responsible for the catalyzing hydrolysis of long-chain triglycerides into fatty acids and glycerol at high temperatures. In the current study, using 30 ns molecular dynamics simulations at different temperatures, an aqueous activation was detected for T1 lipase. This aqueous activation in T1 lipase was mainly caused by a double-flap movement mechanism. The double flaps were constituted by the hydrophobic helices 6 and 9. Helix 6 employed two major components with the hydrophilic part at the surface and the hydrophobic part inside. In the aqueous solution, the hydrophobic part could provide enough power for helix 6 to move away, driving the protein into an open configuration and exposing the catalytic triad. Our findings could provide structural evidence to support the double-flap movement, revealing the catalytic mechanism for T1 lipase.
Erythromycin A and roxithromycin are clinically important macrolide antibiotics that selectively act on the bacterial 50S large ribosomal subunit to inhibit bacteria's protein elongation process by blocking the exit tunnel for the nascent peptide away from ribosome. The detailed molecular mechanism of macrolide binding is yet to be elucidated as it is currently known to the most general idea only. In this study, molecular dynamics (MD) simulation was employed to study their interaction at the molecular level, and the binding free energies for both systems were calculated using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The calculated binding free energies for both systems were slightly overestimated compared to the experimental values, but individual energy terms enabled better understanding in the binding for both systems. Decomposition of results into residue basis was able to show the contribution of each residue at the binding pocket toward the binding affinity of macrolides and hence identified several key interacting residues that were in agreement with previous experimental and computational data. Results also indicated the contributions from van der Waals are more important and significant than electrostatic contribution in the binding of macrolides to the binding pocket. The findings from this study are expected to contribute to the understanding of a detailed mechanism of action in a quantitative matter and thus assisting in the development of a safer macrolide antibiotic.
Increased reports of oseltamivir (OTV)-resistant strains of the influenza virus, such as the H274Y mutation on its neuraminidase (NA), have created some cause for concern. Many studies have been conducted in the attempt to uncover the mechanism of OTV resistance in H274Y NA. However, most of the reported studies on H274Y focused only on the drug-bound system, so the direct effects of the mutation on NA itself prior to drug binding still remain unclear. Therefore, molecular dynamics simulations of NA in apo form, followed by principal component analysis and interaction energy calculations, were performed to investigate the structural changes of the NA binding site as a result of the H274Y mutation. It was observed that the disruption of the NA binding site due to the H274Y mutation was initiated by the repulsive effect of Y274 on the 250-loop, which in turn altered the hydrogen-bonding network around residue 274. The rotated W295 side chain caused the upward movement of the 340-loop. Consequently, sliding box docking results suggested that the binding pathway of OTV was compromised because of the disruption of this binding site. This study also highlighted the importance of the functional group at C6 of the sialic acid mimicry. It is hoped that these results will improve the understanding of OTV resistance and shed some light on the design of a novel anti-influenza drug.
ProBiS is a new method to identify the binding site of protein through local structural alignment against the nonredundant Protein Data Bank (PDB), which may result in unique findings compared to the energy-based, geometry-based, and sequence-based predictors. In this work, binding sites of Hemagglutinin (HA), which is an important target for drugs and vaccines in influenza treatment, have been revisited by ProBiS. For the first time, the identification of conserved binding sites by local structural alignment across all subtypes and strains of HA available in PDB is presented. ProBiS finds three distinctive conserved sites on HA's structure (named Site 1, Site 2, and Site 3). Compared to other predictors, ProBiS is the only one that accurately defines the receptor binding site (Site 1). Apart from that, Site 2, which is located slightly above the TBHQ binding site, is proposed as a potential novel conserved target for membrane fusion inhibitor. Lastly, Site 3, located around Helix A at the stem domain and recently targeted by cross-reactive antibodies, is predicted to be conserved in the latest H7N9 China 2013 strain as well. The further exploration of these three sites provides valuable insight in optimizing the influenza drug and vaccine development.