Jiaguo Lü

Three-Dimensional Model of Lanosterol 14α-Demethylase from Cryptococcus neoformans: Active-Site Characterization and Insights into Azole Binding

ABSTRACT Cryptococcus neoformans is one of the most important causes of life-threatening fungal infections in immunocompromised patients. Lanosterol 14α-demethylase (CYP51) is the target of azole antifungal agents. This study describes, for the first time, the 3-dimensional model of CYP51 from Cryptococcus neoformans (CnCYP51). The model was further refined by energy minimization and molecular-dynamics simulations. The active site of CnCYP51 was well characterized by multiple-copy simultaneous-search calculations, and four functional regions important for rational drug design were identified. The mode of binding of the natural substrate and azole antifungal agents with CnCYP51 was identified by flexible molecular docking. A G484S substitution mechanism for azole resistance in CnCYP51, which might be important for the conformation of the heme environment, is suggested.

Homology Modeling of Lanosterol 14α-Demethylase of Candida albicans and Aspergillus fumigatus and Insights into the Enzyme-Substrate Interactions

Abstract The crystal structure of 14α-sterol demethylase from Mycobacterium tuberculosis(MT_14DM) provides a good template for modeling the three dimensional structure of lanosterol 14α-demethylase, which is the target of azole antifungal agents. Homologous 3D models of lanosterol 14α-demethylase from Candida albicans (CA_14DM) and Aspergillus fumigatus (AF_14DM) were built on the basis of the crystal coordinates of MT_14DM in complex with 4-phenylimidazole and fluconazole. The reliability of the two models was assessed by Ramachandran plots, Profile-3D analysis, and by analyzing the consistency of the two models with the experimental data on the P45014DM. The overall structures of the resulting CA_14DM model and AF_14DM model are similar to those of the template structures. The two models remain the core structure characteristic for cytochrome P450s and most of the insertions and deletions expose the molecular surface. The structurally and functionally important residues such as the heme binding residues, the residues lining the substrate access channel, and residues in active site were identified from the model. To explore the binding mode of the substrate with the two models, 24(28)- methylene-24,25-dihydrolanosterol was docked into the active site of the two models and hydrophobic interaction and hydrogen-bonding were found to play an important role in substrate recognition and orientation. These results provided a basis for experiments to probe structure-function relationships in the P45014DM. Although CA_14DM and AF_14DM shared similar core structural character, the active site of the two models were quite different, thus allowing the rational design of specific inhibitors to the target enzyme and the discovery of novel antifungal agents with broad spectrum.

Evolutionary trace analysis of CYP51 family: implication for site-directed mutagenesis and novel antifungal drug design

Lanosterol 14α-demethylase (CYP51) is an essential enzyme in the fungal life cycle and also an important target for the antifungal drug development. Based on the multiple sequence alignments of CYP51 family, an evolutionary tree of the CYP51 family was constructed by the evolutionary trace (ET) method. The identified trace residues could provide a reliable and rational guide to the design of CYP51 mutations and give more information about the detailed mechanism of substrate (drug) recognition and binding. The reliability of ET analysis to identify residues of functional importance was validated by the reported site-directed mutagenesis studies of CYP51s. Several residues in the active site were also validated by our mutagenesis studies. Mapping the identified trace residues onto the active site of the modeled structure of Candida albicans CYP51 (CACYP51) may provide useful information for the design of novel antifungal agents.

Homology modeling and molecular dynamics simulation of N-myristoyltransferase from protozoan parasites: active site characterization and insights into rational inhibitor design

Myristoyl-CoA:protein N-myristoyltransferase (NMT) is a cytosolic monomeric enzyme that catalyzes the transfer of the myristoyl group from myristoyl-CoA to the N-terminal glycine of a number of eukaryotic cellular and viral proteins. Recent experimental data suggest NMT from parasites could be a promising new target for the design of novel antiparasitic agents with new mode of action. However, the active site topology and inhibitor specificity of these enzymes remain unclear. In this study, three-dimensional models of NMT from Plasmodium falciparum (PfNMT), Leishmania major (LmNMT) and Trypanosoma brucei (TbNMT) were constructed on the basis of the crystal structures of fungal NMTs using homology modeling method. The models were further refined by energy minimization and molecular dynamics simulations. The active sites of PfNMT, LmNMT and TbNMT were characterized by multiple copy simultaneous search (MCSS). MCSS functional maps reveal that PfNMT, LmNMT and TbNMT share a similar active site topology, which is defined by two hydrophobic pockets, a hydrogen-bonding (HB) pocket, a negatively-charged HB pocket and a positively-charged HB pocket. Flexible docking approaches were then employed to dock known inhibitors into the active site of PfNMT. The binding mode, structure–activity relationships and selectivity of inhibitors were investigated in detail. From the results of molecular modeling, the active site architecture and certain key residues responsible for inhibitor binding were identified, which provided insights for the design of novel inhibitors of parasitic NMTs.

Evolutionary trace analysis of eukaryotic DNA topoisomerase I superfamily: Identification of novel antitumor drug binding site

The studies of novel inhibitors of DNA topoisomerase I (Topo I) have already become very promising in cancer chemotherapy. Identifying the new drug-binding residues is playing an important role in the design and optimization of Topo I inhibitors. The designed compounds may have novel scaffolds, thus will be helpful to overcome the toxicities of current camptothecin (CPT) drugs and may provide a solution to cross resistance with these drugs. Multiple sequence alignments were performed on eukaryotic DNA topoisomerase I superfamily and thus the evolutionary tree was constructed. The Evolutionary Trace method was applied to identify functionally important residues of human Topo I. It has been demonstrated that class-specific hydrophobic residues Ala351, Met428, Pro431 are located around the 7,9-position of CPT, indicating suitable substitution of hydrophobic group on CPT will increase antitumor activity. The conservative residue Lys436 in the superfamily is of particular interest and new CPT derivatives designed based on this residue may greatly increase water solubility of such drugs. It has also been demonstrated that the residues Asn352 and Arg364 were conservative in the superfamily, whose mutation will render CPT resistance. As our molecular docking studies demonstrated they did not make any direct interaction with CPT, they are important drug-binding site residues for future design of novel non-camptothecin lead compounds. This work provided a strong basis for the design and synthesis of novel highly potent CPT derivatives and virtual screening for novel lead compounds.

A method of active conformation search based on active and inactive analogues and its application to allylamine antimycotics

A new program ACSBAIA (Active Conformation Search Based on Active and Inactive Analogues) for determination of the active conformations was developed based on the rationales that specific functional groups of active analogues could reach and interact with the active site of target receptor by means of the change of conformations, but that of inactive analogues could not interact with the active site owing to conformational restriction. The program consisted of 4 sub-programs: conformation sampling system, active conformation constraint system, inactive conformation exclusion system, and activity prediction system. Pharmacophoric conformation of allylamine antimycotics was studied by this method. Activities of 2 analogues were predicted and tested. The results suggested that the method was scientific and practical. The application of this method was not restricted by the three-dimensional structural knowledge of target receptor. In the absence of structural information about the receptor, the method was particularly applicable.