PubMed Health⌕ Search

Biomedical subjects

Xiaoqin Zou

Publications and source records attributed to Xiaoqin Zou.

11 recordsLinked to original sources
PubMed2007

Ensemble docking of multiple protein structures: considering protein structural variations in molecular docking.

One approach to incorporate protein flexibility in molecular docking is the use of an ensemble consisting of multiple protein structures. Sequentially docking each ligand into a large number of protein structures is computationally too expensive to allow large-scale database screening. It is challenging to achieve a good balance between docking accuracy and computational efficiency. In this work, we have developed a fast, novel docking algorithm utilizing multiple protein structures, referred to as ensemble docking, to account for protein structural variations. The algorithm can simultaneously dock a ligand into an ensemble of protein structures and automatically select an optimal protein structure that best fits the ligand by optimizing both ligand coordinates and the conformational variable m, where m represents the m-th structure in the protein ensemble. The docking algorithm was validated on 10 protein ensembles containing 105 crystal structures and 87 ligands in terms of binding mode and energy score predictions. A success rate of 93% was obtained with the criterion of root-mean-square deviation <2.5 A if the top five orientations for each ligand were considered, comparable to that of sequential docking in which scores for individual docking are merged into one list by re-ranking, and significantly better than that of single rigid-receptor docking (75% on average). Similar trends were also observed in binding score predictions and enrichment tests of virtual database screening. The ensemble docking algorithm is computationally efficient, with a computational time comparable to that for docking a ligand into a single protein structure. In contrast, the computational time for the sequential docking method increases linearly with the number of protein structures in the ensemble. The algorithm was further evaluated using a more realistic ensemble in which the corresponding bound protein structures of inhibitors were excluded. The results show that ensemble docking successfully predicts the binding modes of the inhibitors, and discriminates the inhibitors from a set of noninhibitors with similar chemical properties. Although multiple experimental structures were used in the present work, our algorithm can be easily applied to multiple protein conformations generated by computational methods, and helps improve the efficiency of other existing multiple protein structure(MPS)-based methods to accommodate protein flexibility.

Algorithms↗
PubMed2006

An iterative knowledge-based scoring function to predict protein-ligand interactions: I. Derivation of interaction potentials.

Using a novel iterative method, we have developed a knowledge-based scoring function (ITScore) to predict protein-ligand interactions. The pair potentials for ITScore were derived from a training set of 786 protein-ligand complex structures in the Protein Data Bank. Twenty-six atom types were used based on the atom type category of the SYBYL software. The iterative method circumvents the long-standing reference state problem in the derivation of knowledge-based scoring functions. The basic idea is to improve pair potentials by iteration until they correctly discriminate experimentally determined binding modes from decoy ligand poses for the ligand-protein complexes in the training set. The iterative method is efficient and normally converges within 20 iterative steps. The scoring function based on the derived potentials was tested on a diverse set of 140 protein-ligand complexes for affinity prediction, yielding a high correlation coefficient of 0.74. Because ITScore uses SYBYL-defined atom types, this scoring function is easy to use for molecular files prepared by SYBYL or converted by software such as BABEL.

Ligands↗
PubMed2006

An iterative knowledge-based scoring function to predict protein-ligand interactions: II. Validation of the scoring function.

We have developed an iterative knowledge-based scoring function (ITScore) to describe protein-ligand interactions. Here, we assess ITScore through extensive tests on native structure identification, binding affinity prediction, and virtual database screening. Specifically, ITScore was first applied to a test set of 100 protein-ligand complexes constructed by Wang et al. (J Med Chem 2003, 46, 2287), and compared with 14 other scoring functions. The results show that ITScore yielded a high success rate of 82% on identifying native-like binding modes under the criterion of rmsd < or = 2 A for each top-ranked ligand conformation. The success rate increased to 98% if the top five conformations were considered for each ligand. In the case of binding affinity prediction, ITScore also obtained a good correlation for this test set (R = 0.65). Next, ITScore was used to predict binding affinities of a second diverse test set of 77 protein-ligand complexes prepared by Muegge and Martin (J Med Chem 1999, 42, 791), and compared with four other widely used knowledge-based scoring functions. ITScore yielded a high correlation of R2 = 0.65 (or R = 0.81) in the affinity prediction. Finally, enrichment tests were performed with ITScore against four target proteins using the compound databases constructed by Jacobsson et al. (J Med Chem 2003, 46, 5781). The results were compared with those of eight other scoring functions. ITScore yielded high enrichments in all four database screening tests. ITScore can be easily combined with the existing docking programs for the use of structure-based drug design.

Acetylcholinesterase↗
PubMed2006

Efficient molecular docking of NMR structures: application to HIV-1 protease.

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.

Algorithms↗
PubMed2006

Similarities and differences between organic cation inhibition of the Na,K-ATPase and PMCA.

The effects of three classes of organic cations on the inhibition of the plasma membrane Ca pump (PMCA) were determined and compared to inhibition of the Na pump. Quaternary amines (tetramethylammonium, tetraethylammonium, and tetrapropylammonium, TMA, TEA, and TPA, respectively) did not inhibit PMCA. This is not to imply that PMCA is inherently selective against monovalent cations because guanidine and tetramethylguanidine inhibited PMCA by competing with Ca(2+). The divalent organic cation, ethyl diamine, inhibited PMCA but was not competitive with Ca(2+). In contrast, propyl diamine did compete with Ca(2+) and was about 10-fold more potent than butyl diamine in inhibiting PMCA. For the Na pump, both TEA and TPA inhibited, but TMA did not. TEA, guanidine, and tetramethylguanidine inhibition was competitive with Na(+) for ATPase activation and with K(+) for pNPPase activation, both of which are cytoplasmic substrate cation effects. Thus, these findings are consistent with TEA, guanidine, and tetramethylguanidine inhibiting from the cytoplasmic side of the Na pump; in contrast, we have previously shown that TPA did not inhibit from the cytoplasmic side. The divalent alkane diamines ethyl, propyl, and butyl diamine all inhibited the Na pump and all competed at the intracellular surface. The order of potency was ED > PD > BD consistent with an optimal size for binding; similarly, for the quaternary amines TMA is apparently too small to make appropriate contacts, and TPA is too large. Homology models based upon the high-resolution SERCA structure are included to contextualize the kinetic observations.

Amino Acid Sequence↗
PubMed2006

The two ATP binding sites of cystic fibrosis transmembrane conductance regulator (CFTR) play distinct roles in gating kinetics and energetics.

Cystic fibrosis transmembrane conductance regulator (CFTR), a member of the ABC (ATP binding cassette) transporter family, is a chloride channel whose activity is controlled by protein kinase-dependent phosphorylation. Opening and closing (gating) of the phosphorylated CFTR is coupled to ATP binding and hydrolysis at CFTR's two nucleotide binding domains (NBD1 and NBD2). Recent studies present evidence that the open channel conformation reflects a head-to-tail dimerization of CFTR's two NBDs as seen in the NBDs of other ABC transporters (Vergani et al., 2005). Whether these two ATP binding sites play an equivalent role in the dynamics of NBD dimerization, and thus in gating CFTR channels, remains unsettled. Based on the crystal structures of NBDs, sequence alignment, and homology modeling, we have identified two critical aromatic amino acids (W401 in NBD1 and Y1219 in NBD2) that coordinate the adenine ring of the bound ATP. Conversion of the W401 residue to glycine (W401G) has little effect on the sensitivity of the opening rate to [ATP], but the same mutation at the Y1219 residue dramatically lowers the apparent affinity for ATP by >50-fold, suggesting distinct roles of these two ATP binding sites in channel opening. The W401G mutation, however, shortens the open time constant. Energetic analysis of our data suggests that the free energy of ATP binding at NBD1, but not at NBD2, contributes significantly to the energetics of the open state. This kinetic and energetic asymmetry of CFTR's two NBDs suggests an asymmetric motion of the NBDs during channel gating. Opening of the channel is initiated by ATP binding at the NBD2 site, whereas separation of the NBD dimer at the NBD1 site constitutes the rate-limiting step in channel closing.

Adenosine Triphosphate↗
PubMed2006

Electrostatics of ligand binding: parametrization of the generalized Born model and comparison with the Poisson-Boltzmann approach.

An accurate and fast evaluation of the electrostatics in ligand-protein interactions is crucial for computer-aided drug design. The pairwise generalized Born (GB) model, a fast analytical method originally developed for studying the solvation of organic molecules, has been widely applied to macromolecular systems, including ligand-protein complexes. However, this model involves several empirical scaling parameters, which have been optimized for the solvation of organic molecules, peptides, and nucleic acids but not for energetics of ligand binding. Studies have shown that a good solvation energy does not guarantee a correct model of solvent-mediated interactions. Thus, in this study, we have used the Poisson-Boltzmann (PB) approach as a reference to optimize the GB model for studies of ligand-protein interactions. Specifically, we have employed the pairwise descreening approximation proposed by Hawkins et al.(1) for GB calculations and DelPhi for PB calculations. The AMBER all-atom force field parameters have been used in this work. Seventeen protein-ligand complexes have been used as a training database, and a set of atomic descreening parameters has been selected with which the pairwise GB model and the PB model yield comparable results on atomic Born radii, the electrostatic component of free energies of ligand binding, and desolvation energies of the ligands and proteins. The energetics of the 15 test complexes calculated with the GB model using this set of parameters also agrees well with the energetics calculated with the PB method. This is the first time that the GB model has been parametrized and thoroughly compared with the PB model for the electrostatics of ligand binding.

Ligands↗
PubMed2005

High affinity ATP/ADP analogues as new tools for studying CFTR gating.

Previous studies using non-hydrolysable ATP analogues and hydrolysis-deficient cystic fibrosis transmembrane conductance regulator (CFTR) mutants have indicated that ATP hydrolysis precedes channel closing. Our recent data suggest that ATP binding is also important in modulating the closing rate. This latter hypothesis predicts that ATP analogues with higher binding affinities should stabilize the open state more than ATP. Here we explore the possibility of using N6-modified ATP/ADP analogues as high-affinity ligands for CFTR gating, since these analogues have been shown to be more potent than native ATP/ADP in other ATP-binding proteins. Among the three N6-modified ATP analogues tested, N6-(2-phenylethyl)-ATP (P-ATP) was the most potent, with a K(1/2) of 1.6 +/- 0.4 microm (>50-fold more potent than ATP). The maximal open probability (P(o)) in the presence of P-ATP was approximately 30% higher than that of ATP, indicating that P-ATP also has a higher efficacy than ATP. Single-channel kinetic analysis showed that as [P-ATP] was increased, the opening rate increased, whereas the closing rate decreased. The fact that these two kinetic parameters have different sensitivities to changes of [P-ATP] suggests an involvement of two different ATP-binding sites, a high-affinity site modulating channel closing and a low affinity site controlling channel opening. The effect of P-ATP on the stability of open states was more evident when ATP hydrolysis was abolished, either by mutating the nucleotide-binding domain 2 (NBD2) Walker B glutamate (i.e. E1371) or by using the non-hydrolysable ATP analogue AMP-PNP. Similar strategies to develop nucleotide analogues with a modified adenine ring could be valuable for future studies of CFTR gating.

ATP-Binding Cassette Transporters↗
PubMed2005

CFTR gating II: Effects of nucleotide binding on the stability of open states.

Previously, we demonstrated that ADP inhibits cystic fibrosis transmembrane conductance regulator (CFTR) opening by competing with ATP for a binding site presumably in the COOH-terminal nucleotide binding domain (NBD2). We also found that the open time of the channel is shortened in the presence of ADP. To further study this effect of ADP on the open state, we have used two CFTR mutants (D1370N and E1371S); both have longer open times because of impaired ATP hydrolysis at NBD2. Single-channel kinetic analysis of DeltaR/D1370N-CFTR shows unequivocally that the open time of this mutant channel is decreased by ADP. DeltaR/E1371S-CFTR channels can be locked open by millimolar ATP with a time constant of approximately 100 s, estimated from current relaxation upon nucleotide removal. ADP induces a shorter locked-open state, suggesting that binding of ADP at a second site decreases the locked-open time. To test the functional consequence of the occupancy of this second nucleotide binding site, we changed the [ATP] and performed similar relaxation analysis for E1371S-CFTR channels. Two locked-open time constants can be discerned and the relative distribution of each component is altered by changing [ATP] so that increasing [ATP] shifts the relative distribution to the longer locked-open state. Single-channel kinetic analysis for DeltaR/E1371S-CFTR confirms an [ATP]-dependent shift of the distribution of two locked-open time constants. These results support the idea that occupancy of a second ATP binding site stabilizes the locked-open state. This binding site likely resides in the NH2-terminal nucleotide binding domain (NBD1) because introducing the K464A mutation, which decreases ATP binding affinity at NBD1, into E1371S-CFTR shortens the relaxation time constant. These results suggest that the binding energy of nucleotide at NBD1 contributes to the overall energetics of the open channel conformation.

Adenosine Diphosphate↗
PubMed2004

Detailed kinetic studies of an aggregating inhibitor; inhibition of phosphomannomutase/phosphoglucomutase by disperse blue 56.

Phosphomannomutase/phosphoglucomutase occupies a central position in the pathways by which several virulence factors are synthesized in Pseudomonas aeruginosa. Virtual screening was used to identify potential inhibitors of phosphomannomutase/ phosphoglucomutase, and one compound, the anthraquinone-based dye Disperse Blue 56, showed potent inhibition in vitro. The kinetics of inhibition was complex; the time courses for reactions in the presence of the inhibitor were biphasic, suggestive of slow-binding inhibition. Quantitative analysis of the progress curves and preincubation experiments demonstrated that slow-binding inhibition was not occurring, however. Initial velocity kinetic studies indicated that Disperse Blue 56 was a parabolic, noncompetitve inhibitor. Progress curves for reactions in the presence of Disperse Blue 56 could be fitted very well by a model in which 2 equiv of the inhibitor bound to free enzyme or the enzyme-substrate complex. The inhibition was largely relieved by the inclusion of 0.01% Triton X-100 in the assay solutions, which has been suggested to be the hallmark for inhibition by compounds that exert their effect through aggregates [McGovern, S. L., Caselli, E., Grigorieff, N., and Shiochet, B. K. (2002) J. Med. Chem. 45, 1712-1722]. Our kinetic data appear to be consistent with either inhibition by a dimer of Disperse Blue 56 or inhibition by a Disperse Blue 56 aggregate, but the latter appears much more likely. We present a detailed analysis of the system to provide further information that may help in the recognition of inhibition through aggregation.

Anthraquinones↗