Ligand tunnels in T. brucei and human CYP51: Insights for parasite-specific drug design

Identification of Potent and Selective Inhibitors of Acanthamoeba: Structural Insights into Sterol 14α-Demethylase as a Key Drug Target

Acanthamoeba are free-living pathogenic protozoa that cause blinding keratitis, disseminated infection, and granulomatous amebic encephalitis, which is generally fatal. The development of efficient and safe drugs is a critical unmet need. Acanthamoeba sterol 14α-demethylase (CYP51) is an essential enzyme of the sterol biosynthetic pathway. Repurposing antifungal azoles for amoebic infections has been reported, but their inhibitory effects on Acanthamoeba CYP51 enzymatic activity have not been studied. Here, we report catalytic properties, inhibition, and structural characterization of CYP51 from Acanthamoeba castellanii. The enzyme displays a 100-fold substrate preference for obtusifoliol over lanosterol, supporting the plant-like cycloartenol-based pathway in the pathogen. The strongest inhibition was observed with voriconazole (1 h IC50 0.45 μM), VT1598 (0.25 μM), and VT1161 (0.20 μM). The crystal structures of A. castellanii CYP51 with bound VT1161 (2.24 Å) and without an inhibitor (1.95 Å), presented here, can be used in the development of azole-based scaffolds to achieve optimal amoebicidal effectiveness.

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

Proteomic analysis capsule synthesis and redox mechanisms in the intracellular survival of group B Streptococcus in fish microglia

Group B Streptococcus (GBS) causes meningitis in neonates and Nile tilapia (Oreochromis niloticus). The molecular mechanisms regulating the intracellular survival of this pathogen in the host cell are complex and crucial for the progression of infection. Thus, we propose the use of GBS-infected Nile tilapia microglia as an in vitro model system simulating infection caused by homologous bacteria in humans. We used this model to evaluate the phagocytic activity, as well as the functional aspects of the capsular proteins A, B, C, and D and the major redox enzymes, and the synergistic role of mechanisms/proteins involved in blocking phagocytic process. We observed that in the intracellular phase, GBS showed enhanced synthesis of the polysaccharide capsule and used superoxide dismutase, thioredoxin, NADH oxidase, and alkyl hydroperoxide reductase to scavenge reactive oxygen species and reactive nitrogen species produced by the host cell. Furthermore, although these virulence mechanisms were effective during the initial hours of infection, they were not able to subvert microglial responses, which partially neutralized the infection. Altogether, our findings provided important information regarding the intracellular survival mechanisms of GBS and perspectives for the production of new drugs and vaccines, through the druggability analysis of specific proteins. In conclusion, tilapia microglia serve as a potent in vitro experimental model for the study of meningitis.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Molecular Dynamics Investigations of Binding Mechanism for Triazoles Inhibitors to CYP51

The sterol 14α demethylase enzyme (CYP51) is an important target of fungal infections. However, the molecular mechanism between triazoles inhibitors and CYP51 remains obscure. In this study, we have investigated the binding mechanism and tunnel characteristic upon four triazoles inhibitors with CYP51 based on the molecular docking and molecular dynamics simulations. The results indicate the four inhibitors stabilize in the binding cavity of CYP51 in a similar binding mode. We discover a hydrophobic cavity (F58, Y64, Y118, L121, Y132, L376, S378, S506, S507, and M508) and the hydrophobic interaction is the main driving force for inhibitors binding to CYP51. The long-tailed inhibitors (posaconazole and itraconazole) have stronger binding affinities than short-tailed inhibitors (fluconazole and voriconazole) because long-tailed inhibitors can form more hydrophobic interactions with CYP51. The tunnel 2f is the predominant pathway for inhibitors ingress/egress protein, which is similar to the other works of CYP51. This study could provide the theoretical basis for the development of efficient azoles inhibitors and may lead a better insight into structure-function relationships of CYP51.

The Effect of Force-Field Parameters on Cytochrome P450-Membrane Interactions: Structure and Dynamics

The simulation of membrane proteins requires compatible protein and lipid force fields that reproduce the properties of both the protein and the lipid bilayer. Cytochrome P450 enzymes are bitopic membrane proteins with a transmembrane helical anchor and a large cytosolic globular domain that dips into the membrane. As such, they are representative and challenging examples of membrane proteins for simulations, displaying features of both peripheral and integral membrane proteins. We performed molecular dynamics simulations of three cytochrome P450 isoforms (2C9, 2C19 and 1A1) in a 2-oleoyl-1-palmitoyl-sn-glycerol-3-phosphocholine bilayer using two AMBER force field combinations: GAFF-LIPID with ff99SB for the protein, and LIPID14 with ff14SB for the protein. Comparison of the structural and dynamic properties of the proteins, the lipids and the protein-membrane interactions shows differing sensitivity of the cytochrome P450 isoforms to the choice of force field, with generally better agreement with experiment for the LIPID14 + ff14SB combination.

Advances in preclinical approaches to Chagas disease drug discovery

Introduction: Chagas disease affects 8-10 million people worldwide, mainly in Latin America. The current therapy for Chagas disease is limited to nifurtimox and benznidazole, which are effective in treating only the acute phase of the disease but with severe side effects. Therefore, there is an unmet need for new drugs and for the exploration of innovative approaches which may lead to the discovery of new effective and safe drugs for its treatment. Areas covered: The authors report and discuss recent approaches including structure-based design that have led to the discovery of new promising small molecule candidates for Chagas disease which affect prime targets that intervene in the sterol pathway of T. cruzi. Other trypanosome targets, phenotypic screening, the use of artificial intelligence and the challenges with Chagas disease drug discovery are also discussed. Expert opinion: The application of recent scientific innovations to the field of Chagas disease have led to the discovery of new promising drug candidates for Chagas disease. Phenotypic screening brought new hits and opportunities for drug discovery. Artificial intelligence also has the potential to accelerate drug discovery in Chagas disease and further research into this is warranted.

Differing Membrane Interactions of Two Highly Similar Drug-Metabolizing Cytochrome P450 Isoforms: CYP 2C9 and CYP 2C19

The human cytochrome P450 (CYP) 2C9 and 2C19 enzymes are two highly similar isoforms with key roles in drug metabolism. They are anchored to the endoplasmic reticulum membrane by their N-terminal transmembrane helix and interactions of their cytoplasmic globular domain with the membrane. However, their crystal structures were determined after N-terminal truncation and mutating residues in the globular domain that contact the membrane. Therefore, the CYP-membrane interactions are not structurally well-characterized and their dynamics and the influence of membrane interactions on CYP function are not well understood. We describe herein the modeling and simulation of CYP 2C9 and CYP 2C19 in a phospholipid bilayer. The simulations revealed that, despite high sequence conservation, the small sequence and structural differences between the two isoforms altered the interactions and orientations of the CYPs in the membrane bilayer. We identified residues (including K72, P73, and I99 in CYP 2C9 and E72, R73, and H99 in CYP 2C19) at the protein-membrane interface that contribute not only to the differing orientations adopted by the two isoforms in the membrane, but also to their differing substrate specificities by affecting the substrate access tunnels. Our findings provide a mechanistic interpretation of experimentally observed effects of mutagenesis on substrate selectivity.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways

The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.

CYP51 as drug targets for fungi and protozoan parasites: past, present and future

The efficiency of treatment of human infections with the unicellular eukaryotic pathogens such as fungi and protozoa remains deeply unsatisfactory. For example, the mortality rates from nosocomial fungemia in critically ill, immunosuppressed or post-cancer patients often exceed 50%. A set of six systemic clinical azoles [sterol 14α-demethylase (CYP51) inhibitors] represents the first-line antifungal treatment. All these drugs were discovered empirically, by monitoring their effects on fungal cell growth, though it had been proven that they kill fungal cells by blocking the biosynthesis of ergosterol in fungi at the stage of 14α-demethylation of the sterol nucleus. This review briefs the history of antifungal azoles, outlines the situation with the current clinical azole-based drugs, describes the attempts of their repurposing for treatment of human infections with the protozoan parasites that, similar to fungi, also produce endogenous sterols, and discusses the most recently acquired knowledge on the CYP51 structure/function and inhibition. It is our belief that this information should be helpful in shifting from the traditional phenotypic screening to the actual target-driven drug discovery paradigm, which will rationalize and substantially accelerate the development of new, more efficient and pathogen-oriented CYP51 inhibitors.

Crystal Structures of Full-Length Lanosterol 14α-Demethylases of Prominent Fungal Pathogens Candida albicans and Candida glabrata Provide Tools for Antifungal Discovery

Targeting lanosterol 14α-demethylase (LDM) with azole drugs provides prophylaxis and treatments for superficial and disseminated fungal infections, but cure rates are not optimal for immunocompromised patients and individuals with comorbidities. The efficacy of azole drugs has also been reduced due to the emergence of drug-resistant fungal pathogens. We have addressed the need to improve the potency, spectrum, and specificity for azoles by expressing in Saccharomyces cerevisiae functional, recombinant, hexahistidine-tagged, full-length Candida albicans LDM (CaLDM6×His) and Candida glabrata LDM (CgLDM6×His) and determining their X-ray crystal structures. The crystal structures of CaLDM6×His, CgLDM6×His, and ScLDM6×His have the same fold and bind itraconazole in nearly identical conformations. The catalytic domains of the full-length LDMs have the same fold as the CaLDM6×His catalytic domain in complex with posaconazole, with minor structural differences within the ligand binding pocket. Our structures give insight into the LDM reaction mechanism and phenotypes of single-site CaLDM mutations. This study provides a practical basis for the structure-directed discovery of novel antifungals that target LDMs of fungal pathogens.

Sterol 14α-Demethylase Structure-Based Design of VNI (( R)- N-(1-(2,4-Dichlorophenyl)-2-(1 H-imidazol-1-yl)ethyl)-4-(5-phenyl-1,3,4-oxadiazol-2-yl)benzamide)) Derivatives To Target Fungal Infections: Synthesis, Biological Evaluation, and Crystallographic Analysis

Because of the increase in the number of immunocompromised patients, the incidence of invasive fungal infections is growing, but the treatment efficiency remains unacceptably low. The most potent clinical systemic antifungals (azoles) are the derivatives of two scaffolds: ketoconazole and fluconazole. Being the safest antifungal drugs, they still have shortcomings, mainly because of pharmacokinetics and resistance. Here, we report the successful use of the target fungal enzyme, sterol 14α-demethylase (CYP51), for structure-based design of novel antifungal drug candidates by minor modifications of VNI [( R)- N-(1-(2,4-dichlorophenyl)-2-(1 H-imidazol-1-yl)ethyl)-4-(5-phenyl-1,3,4-oxadiazol-2-yl)benzamide)], an inhibitor of protozoan CYP51 that cures Chagas disease. The synthesis of fungi-oriented VNI derivatives, analysis of their potencies to inhibit CYP51s from two major fungal pathogens ( Aspergillus fumigatus and Candida albicans), microsomal stability, effects in fungal cells, and structural characterization of A. fumigatus CYP51 in complexes with the most potent compound are described, offering a new antifungal drug scaffold and outlining directions for its further optimization.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

Trypanosoma brucei CYP51: Essentiality and Targeting Therapy in an Experimental Model

Trypanosoma brucei gambiense is the main causative agent of Human African Trypanosomiasis (HAT), also known as sleeping sickness. Because of limited alternatives and treatment toxicities, new therapeutic options are urgently needed for patients with HAT. Sterol 14alpha-demethylase (CYP51) is a potential drug target but its essentiality has not been determined in T. brucei. We used a tetracycline-inducible RNAi system to assess the essentiality of CYP51 in T. brucei bloodstream form (BSF) cells and we evaluated the effect of posaconazole, a well-tolerated triazole drug, within a panel of virulent strains in vitro and in a murine model. Expression of CYP51 in several T. brucei cell lines was demonstrated by western blot and its essentiality was demonstrated by RNA interference (CYP51RNAi) in vitro. Following reduction of TbCYP51 expression by RNAi, cell growth was reduced and eventually stopped compared to WT or non-induced cells, showing the requirement of CYP51 in T. brucei. These phenotypes were rescued by addition of ergosterol. Additionally, CYP51RNAi induction caused morphological defects with multiflagellated cells (p<0.05), suggesting cytokinesis dysfunction. The survival of CYP51RNAi Doxycycline-treated mice (p = 0.053) and of CYP51RNAi 5-day pre-induced Doxycycline-treated mice (p = 0.008) were improved compared to WT showing a CYP51 RNAi effect on trypanosomal virulence in mice. The posaconazole concentrations that inhibited parasite growth by 50% (IC50) were 8.5, 2.7, 1.6 and 0.12 μM for T. b. brucei 427 90-13, T. b. brucei Antat 1.1, T. b. gambiense Feo (Feo/ITMAP/1893) and T. b. gambiense Biyamina (MHOM/SD/82), respectively. During infection with these last three virulent strains, posaconazole-eflornithine and nifurtimox-eflornithine combinations showed similar improvement in mice survival (p≤0.001). Our results provide support for a CYP51 targeting based treatment in HAT. Thus posaconazole used in combination may represent a therapeutic alternative for trypanosomiasis.

Human sterol 14α-demethylase as a target for anticancer chemotherapy: towards structure-aided drug design

Rapidly multiplying cancer cells synthesize greater amounts of cholesterol to build their membranes. Cholesterol-lowering drugs (statins) are currently in clinical trials for anticancer chemotherapy. However, given at higher doses, statins cause serious side effects by inhibiting the formation of other biologically important molecules derived from mevalonate. Sterol 14α-demethylase (CYP51), which acts 10 steps downstream, is potentially a more specific drug target because this portion of the pathway is fully committed to cholesterol production. However, screening a variety of commercial and experimental inhibitors of microbial CYP51 orthologs revealed that most of them (including all clinical antifungals) weakly inhibit human CYP51 activity, even if they display high apparent spectral binding affinity. Only one relatively potent compound, (R)-N-(1-(3,4'-difluorobiphenyl-4-yl)-2-(1H-imidazol-1-yl)ethyl)-4-(5-phenyl-1,3,4-oxadiazol-2-yl)benzamide (VFV), was identified. VFV has been further tested in cellular experiments and found to decrease proliferation of different cancer cell types. The crystal structures of human CYP51-VFV complexes (2.0 and 2.5 Å) both display a 2:1 inhibitor/enzyme stoichiometry, provide molecular insights regarding a broader substrate profile, faster catalysis, and weaker susceptibility of human CYP51 to inhibition, and outline directions for the development of more potent inhibitors.