Aflatoxin B1 (AFB1)—the most potent natural carcinogen known to men—is metabolized by cytochrome P450 3A4 (CYP3A4), either to the genotoxic AFB1 exo-8,9-epoxide or to the detoxified 3α-hydroxy AFB1. The activation of the procarcinogen proceeds in a highly cooperative fashion, which differs from common allosteric regulation in the sense that it can be attributed to simultaneous occupancy of a single large and malleable active site by multiple ligand molecules. Unfortunately, unlike in the case of ketoconazole, there is currently no experimental structure available for the doubly ligated CYP3A4-AFB1 complex. Therefore, we employed a sequential molecular docking protocol to create various possible doubly ligated complexes and subsequently performed molecular dynamics simulations and free-energy calculations to check for their consistency with the available experimental data on regio- and stereoselectivity of both AFB1 oxidations as well as with available kinetic data. Only the system in which the first AFB1 molecule was bound in a face-on C8–C9 epoxidation mode and the second AFB1 molecule was bound in a side-on 3α-hydroxylation mode—a result of an unconstrained molecular docking protocol—has successfully fulfilled all the imposed criteria and is therefore proposed as the most likely structure of the doubly ligated complex of CYP3A4 with AFB1. The empirical Linear Interaction Energy method revealed that shape complementarity through nonpolar dispersion interactions between the two bound AFB1 molecules is the main source of the experimentally observed positive homotropic cooperativity. The reported study represents a nice example of how state-of-the-art molecular modeling techniques can be used to study complicated macromolecular complexes, whose structures have not yet been experimentally determined, and to validate these against the available experimental data. The proposed structure will facilitate future studies on the rational design of successful AFB1 modulators or on human subpopulations characterized by specific CYP3A4 polymorphisms that are especially sensitive to AFB1.
COBISS.SI-ID: 5609242
The outbreak of avian influenza A (H5N1) virus has raised a global concern for both the animal as well as human health. Besides vaccination, that may not achieve full protection in certain groups of patients, inhibiting neuraminidase or the transmembrane protein M2 represents the main measure of controlling the disease. Due to alarming emergence of influenza virus strains resistant to the currently available drugs, development of new neuraminidase N1 inhibitors is of utmost importance. The present paper provides an overview of the recent advances in the design of new antiviral drugs against avian influenza. It also reports findings in binding free energy calculations for nine neuraminidase N1 inhibitors (oseltamivir, zanamivir, and peramivir -carboxylate, -phosphonate, and -sulfonate) using the Linear Interaction Energy method. Molecular dynamics simulations of these inhibitors were performed in a free and two bound states – the so called ˝open˝ and ˝closed˝ conformations of neuraminidase N1. Obtained results successfully reproduce the experimental binding affinities of the already known neuraminidase N1 inhibitors, i.e. peramivir being a stronger binder than zanamivir that is in turn stronger binder than oseltamivir, or phosphonate inhibitors being stronger binders than their carboxylate analogues. In addition, the newly proposed sulfonate inhibitors are predicted to be the strongest binders – a fact to be confirmed by their chemical synthesis and a subsequent test of their biological activity. Finally, contributions of individual inhibitor moieties to the overall binding affinity are explicitly evaluated to assist further drug development towards inhibition of the H5N1 avian influenza A virus.
COBISS.SI-ID: 5396250
Sterol 14alpha-demethylase (CYP51) is the main drug target for the treatment of fungal infections. The discovery of new efficient fungal CYP51 inhibitors requires an understanding of the structural requirements for selectivity for the fungal over the human ortholog. In this study, a binding mode of the pyridylethanol(phenylethyl)amine type CYP51 inhibitor to the human ortholog was determined at the atomic level. We isolated and purified a full-length human CYP51. The inhibitor-specific binding and its conformational and dynamic properties were evaluated using UV-visible and NMR spectroscopy. Considering the experimental data in docking calculations and molecular dynamics simulations, the location of the inhibitor moieties and their interactions with the enzyme active site were determined. The inhibitor binds to the enzyme in two diastereomeric forms, which have common location of aromatic ring moieties, while the less bulky propyl chain can adapt to various hydrophobic regions of the enzyme active site. The halogenated phenyl ring binds in the substrate access channel making numerous contacts with the hydrophobic side chains, and its interactions with the unconserved residues are especially informative. The results reveal the unique binding properties of the investigated inhibitor in comparison to the azoles and provide novel directions for the design of selective fungal inhibitors.
COBISS.SI-ID: 31710937