Biomedical subjects

Niu Huang

Publications and source records attributed to Niu Huang.

12 recordsLinked to original sources

PubMed2006

Benchmarking sets for molecular docking.

Ligand enrichment among top-ranking hits is a key metric of molecular docking. To avoid bias, decoys should resemble ligands physically, so that enrichment is not simply a separation of gross features, yet be chemically distinct from them, so that they are unlikely to be binders. We have assembled a directory of useful decoys (DUD), with 2950 ligands for 40 different targets. Every ligand has 36 decoy molecules that are physically similar but topologically distinct, leading to a database of 98,266 compounds. For most targets, enrichment was at least half a log better with uncorrected databases such as the MDDR than with DUD, evidence of bias in the former. These calculations also allowed 40x40 cross-docking, where the enrichments of each ligand set could be compared for all 40 targets, enabling a specificity metric for the docking screens. DUD is freely available online as a benchmarking set for docking at http://blaster.docking.org/dud/.

Adenosine Deaminase↗

PubMed2006

Molecular mechanics methods for predicting protein-ligand binding.

Ligand binding affinity prediction is one of the most important applications of computational chemistry. However, accurately ranking compounds with respect to their estimated binding affinities to a biomolecular target remains highly challenging. We provide an overview of recent work using molecular mechanics energy functions to address this challenge. We briefly review methods that use molecular dynamics and Monte Carlo simulations to predict absolute and relative ligand binding free energies, as well as our own work in which we have developed a physics-based scoring method that can be applied to hundreds of thousands of compounds by invoking a number of simplifying approximations. In our previous studies, we have demonstrated that our scoring method is a promising approach for improving the discrimination between ligands that are known to bind and those that are presumed not to, in virtual screening of large compound databases. In new results presented here, we explore several improvements to our computational method including modifying the dielectric constant used for the protein and ligand interiors, and empirically scaling energy terms to compensate for deficiencies in the energy model. Future directions for further improving our physics-based scoring method are also discussed.

Algorithms↗

PubMed2006

Conserved patterns in backbone torsional changes allow for single base flipping from duplex DNA with minimal distortion of the double helix.

Base flipping is a structural mechanism common to many DNA processing and repair enzymes. Changes in the local backbone torsions that occur during base flipping and the effect of environment on their behavior are of particular interest in understanding different base flipping mechanisms. In the present study, structures sampled during umbrella sampling molecular dynamics (MD) simulations of base flipping in aqueous and protein-bound environments, carried out with two different MD simulation strategies, are analyzed to find the most significant phosphodiester backbone distortions in the vicinity of the flipping base. Torsional sampling on the 5' side of the flipping base during flipping through the major groove shows similarities to the torsional sampling on the 3' side during flipping through the minor groove and vice versa. In differing environments, this behavior varies only marginally. These compensating torsional changes in the DNA backbone on 5' and 3' sides of the flipping base limit overall distortion of the DNA double helix during single base flipping. Rotameric intermediate states observed during base flipping are identified and postulated to be metastable states implicated in both large-scale structural changes and functional effects of chemical modifications in DNA.

DNA↗

PubMed2005

Specificity in protein-DNA interactions: energetic recognition by the (cytosine-C5)-methyltransferase from HhaI.

Sequence-specific interactions between proteins and DNA are essential for a variety of biological functions. The (cytosine-C5)-methyltransferase from HhaI (M.HhaI) specifically modifies the second base in GCGC sequences, employing a base flipping mechanism to access the target base being chemically modified. The mechanism of sequence-specific recognition of M.HhaI is not evident based on crystallographic structures, leading to the suggestion that recognition is linked to the flipping event itself, a process that may be referred to as energetic recognition. Using computational methods, it is shown that the free energy barriers to flipping are significantly higher in non-cognate versus the cognate sequence, supporting the energetic recognition mechanism. Energetic recognition is imparted by two protein "selectivity filters" that function via a "web" of protein-DNA interactions in short-lived, high energy states present along the base flipping pathway. Other sequence-specific DNA binding proteins whose function involves significant distortion of DNA's conformation may use a similar recognition mechanism.

Crystallography, X-Ray↗

PubMed2004

Caught in the act: visualization of an intermediate in the DNA base-flipping pathway induced by HhaI methyltransferase.

Rotation of a DNA or RNA nucleotide out of the double helix and into a protein pocket ('base flipping') is a mechanistic feature common to some DNA/RNA-binding proteins. Here, we report the structure of HhaI methyltransferase in complex with DNA containing a south-constrained abasic carbocyclic sugar at the target site in the presence of the methyl donor byproduct AdoHcy. Unexpectedly, the locked south pseudosugar appears to be trapped in the middle of the flipping pathway via the DNA major groove, held in place primarily through Van der Waals contacts with a set of invariant amino acids. Molecular dynamics simulations indicate that the structural stabilization observed with the south-constrained pseudosugar will not occur with a north-constrained pseudosugar, which explains its lowered binding affinity. Moreover, comparison of structural transitions of the sugar and phosphodiester backbone observed during computational studies of base flipping in the M.HhaI-DNA-AdoHcy ternary complex indicate that the south-constrained pseudosugar induces a conformation on the phosphodiester backbone that corresponds to that of a discrete intermediate of the base-flipping pathway. As previous crystal structures of M.HhaI ternary complex with DNA displayed the flipped sugar moiety in the antipodal north conformation, we suggest that conversion of the sugar pucker from south to north beyond the middle of the pathway is an essential part of the mechanism through which flipping must proceed to reach its final destination. We also discuss the possibility of the south-constrained pseudosugar mimicking a transition state in the phosphodiester and sugar moieties that occurs during DNA base flipping in the presence of M.HhaI.

Binding Sites↗

PubMed2004

Atomistic view of base flipping in DNA.

Base flipping is essential for the enzyme-catalysed methylation of DNA. In our previous studies, the flipping of bases out of duplex DNA in DNA alone and when bound to the (cytosine-C5)-methyltransferase from HhaI (M.HhaI) were investigated via potential of mean force calculations. Insights into various experimental observations were obtained. In the present paper we present an overview of previous computational studies of base flipping along with new detailed structural and energetic analysis on atomic events that contribute to the free energy surfaces. The contributions from different intrinsic and environmental effects to the base-flipping process are explored, and experimental data derived from a variety of methods are reconciled. A detailed protein-facilitated base-flipping mechanism is proposed. Ground-state destabilization is achieved via disruption of the target base Watson-Crick interactions by substitution with favourable DNA-protein interactions. In addition, specific DNA-protein interactions and favourable solvation effects further promote target base flipping along the major groove through the protein matrix, and maximal interactions occur between the DNA and the protein upon reaching the fully flipped state. Other DNA binding proteins that involve distortion of DNA's conformation may use a similar mechanism to that by which M.HhaI facilitates base flipping.

Base Pairing↗

PubMed2004

Identification of non-phosphate-containing small molecular weight inhibitors of the tyrosine kinase p56 Lck SH2 domain via in silico screening against the pY + 3 binding site.

The protein p56 lymphoid T cell tyrosine kinase (Lck) is predominantly expressed in T lymphocytes where it plays a critical role in T-cell-mediated immune response. Lck participates in phosphotyrosine-dependent protein-protein interactions through its modular binding unit, the Src homology-2 (SH2) domain. Accordingly, virtual screening methods combined with experimental assays were used to identify small molecular weight nonpeptidic compounds that block Lck SH2 domain-dependent interactions. Virtual screening included scoring normalization procedures and postdocking structural clustering that is shown to facilitate the selection of active compounds. By targeting the well-defined hydrophobic binding pocket known to impart specificity on Lck-protein interactions (i.e., pY + 3 site), inhibitors of the Lck SH2 domain were discovered that omit the phosphotyrosine (pY) or related moieties. The 34 out of 196 computationally selected compounds were shown to inhibit Lck SH2 domain association with phosphorylated immunoreceptor tyrosine based activation motifs peptide. Twenty-four of the active compounds were further tested for their ability to modulate biological function. Thirteen of these compounds showed inhibitory activity in mixed lymphocyte culture assay. Fluorescence titration experiments on four of these active compounds further verified their binding to the SH2 domain. Because of their simple chemical structures, these small organic compounds have the potential to act as lead compounds for the development of novel immunosuppressant drugs.

Animals↗

PubMed2004

[Molecular design and biological activity of BCR-ABL tyrosine kinase inhibitors].

OBJECTIVE: To discover BCR-ABL tyrosine kinase inhibitors through structure based virtual screening. METHODS: Docking screening against the distinctive inactive conformation of the catalytic domain of BCR-ABL tyrosine kinase was performed on 3D database. The MTT assay was performed to assess the viability of the tumor cells treated with selected compounds. The amount and kinase activity of BCR-ABL protein were detected in the presence of compounds by Western blot analysis and immunoprecipitation. RESULTS: From the top 1,000 compounds with the best DOCK energy score, 15 compounds were selected for biological assay. Eight out of 15 compounds showed notable inhibitory activity against Ph+ human K562 cells with IC50 values ranging from 10 to 200 micromol/L. In cell-based assays of ABL tyrosine phosphorylation, the ability of two kinds of novel, structurally diverse, lead compounds to inhibit ABL kinase activity was observed. However, no significant differences in the amount of BCR-ABL protein were noted on the ABL immunoblot in the presence of these lead compounds. CONCLUSIONS: Two promising lead compounds were discovered to inhibit BCR-ABL tyrosine kinase activity. Virtual screening technique has been proven to narrow down the size of screening compound libraries to the most prospective drug candidates with high success rates.

Computer Simulation↗

PubMed2003

Identification of novel inhibitors of BCR-ABL tyrosine kinase via virtual screening.

Inhibition of BCR-ABL tyrosine kinase activity has shown to be essential for the treatment of chronic myelogenous leukemia (CML). However, drug resistance has quickly arisen in recent clinical trials for STI571 (Gleevec), which is the first approved drug of CML by inhibiting ABL tyrosine kinase. It is desirable to develop new types of ABL tyrosine kinase inhibitors that may overcome this drug resistance problem. Here we present the discovery of novel inhibitors targeted at the catalytic domain of ABL tyrosine kinase by using three-dimensional database searching techniques. From a database containing 200,000 commercially available compounds, the top 1000 compounds with the best DOCK energy score were selected and subjected to structural diversity and drug likeness analysis, 15 compounds were submitted for biological assay. Eight out of the 15 showed inhibitory activity against K562 cells with IC(50) value ranging from 10 to 200 microM. Two promising compounds showed inhibition in further ABL tyrosine phosphorylation assay. It is anticipated that those two compounds can serve as lead compounds for further drug design and optimization.

Binding Sites↗

PubMed2002

Protein-facilitated base flipping in DNA by cytosine-5-methyltransferase.

DNA methylation, various DNA repair mechanisms, and possibly early events in the opening of DNA as required for transcription and replication are initiated by flipping of a DNA base out of the DNA double helix. The energetics and structural mechanism of base flipping in the presence of the DNA-processing enzyme, cytosine 5-methyltransferase from HhaI (M.HhaI), were obtained through molecular dynamics based upon free-energy calculations. Free-energy profiles for base flipping show that, when in the closed conformation, M.HhaI lowers the free-energy barrier to flipping by 17 kcalmol and stabilizes the fully flipped state. Flipping is shown to occur via the major groove of the DNA. Structural analysis indicates that flipping is facilitated by destabilization of the DNA double-helical structure and substitution of DNA base-pairing and base-stacking interactions with DNA-protein interactions. The fully flipped state is stabilized by DNA-protein interactions that are enhanced upon binding of coenzyme. This study represents an atomic detail description of the mechanism by which a protein facilitates specific structural distortion in DNA.

Base Composition↗

PubMedDate not supplied

Consideration of molecular weight during compound selection in virtual target-based database screening.

Virtual database screening allows for millions of chemical compounds to be computationally selected based on structural complimentary to known inhibitors or to a target binding site on a biological macromolecule. Compound selection in virtual database screening when targeting a biological macromolecule is typically based on the interaction energy between the chemical compound and the target macromolecule. In the present study it is shown that this approach is biased toward the selection of high molecular weight compounds due to the contribution of the compound size to the energy score. To account for molecular weight during energy based screening, we propose normalization strategies based on the total number of heavy atoms in the chemical compounds being screened. This approach is computationally efficient and produces molecular weight distributions of selected compounds that can be selected to be (1) lower than that of the original database used in the virtual screening, which may be desirable for selection of leadlike compounds or (2) similar to that of the original database, which may be desirable for the selection of drug-like compounds. By eliminating the bias in target-based database screening toward higher molecular weight compounds it is anticipated that the proposed procedure will enhance the success rate of computer-aided drug design.

Computer-Aided Design↗

PubMedDate not supplied

Physics-based scoring of protein-ligand complexes: enrichment of known inhibitors in large-scale virtual screening.

We demonstrate that using an all-atom molecular mechanics force field combined with an implicit solvent model for scoring protein-ligand complexes is a promising approach for improving inhibitor enrichment in the virtual screening of large compound databases. The rescoring method is evaluated by the extent to which known binders for nine diverse, therapeutically relevant enzymes are enriched against a background of approximately 100,000 drug-like decoys. The improvement in enrichment is most robust and dramatic within the top 1% of the ranked database, that is, the first thousand compounds; below the first few percent of the ranked database, there is little overall improvement. The improved early enrichment is likely due to the more realistic treatment of ligand and receptor desolvation in the rescoring procedure. We also present anecdotal but encouraging results assessing the ability of the rescoring method to predict specificity of inhibitors for structurally related proteins.

Biophysical Phenomena↗