Molecular Dynamics Simulations to Investigate the Influences of Amino Acid Mutations on Protein Three-Dimensional Structures of Cytochrome P450 2D6.1, 2, 10, 14A, 51, and 62

The genomic landscape of CYP2D6 variation in the Indian population

Aim: The CYP2D6 gene is highly polymorphic, causing large interindividual variability in the metabolism of several clinically important drugs. Materials & methods: The authors investigated the diversity and distribution of CYP2D6 alleles in Indians using whole genome sequences (N = 1518). Functional consequences were assessed using pathogenicity scores and molecular dynamics simulations. Results: The analysis revealed population-specific CYP2D6 alleles (*86, *7, *111, *112, *113, *99) and remarkable differences in variant and phenotype frequencies with global populations. The authors observed that one in three Indians could benefit from a dose alteration for psychiatric drugs with accurate CYP2D6 phenotyping. Molecular dynamics simulations revealed large conformational fluctuations, confirming the predicted reduced function of *86 and *113 alleles. Conclusion: The findings emphasize the utility of comprehensive CYP2D6 profiling for aiding precision public health.

[Metabolic Activities Catalyzed by Human Cytochrome P450 (CYP) 2D6 and CYP3A Subfamily Members and Effect of Various Compounds, Including Endogenous Steroid Hormones, on These Activities]

My research focused on the effects of various drugs on (1) dopamine formation from p-tyramine catalyzed by polymorphic cytochrome P450 (CYP or P450) 2D6 variants and (2) endogenous steroid hormone hydroxylation catalyzed by CYP3A subfamily members (CYP3A4, CYP3A5, CYP3A7). The activation (cooperativity) of metabolic reactions catalyzed by P450s was especially emphasized. The effects of various psychotropic agents on dopamine formation from p-tyramine, catalyzed by wild-type CYP2D6.1 and CYP2D6 variants, including CYP2D6.2 (Arg296Cys;Ser486Thr), CYP2D6.10 (Pro34Ser;Ser486Thr), and CYP2D6.39 (Ser486Thr) were compared. Michaelis (Km) and inhibition (Ki) constants of the psychotropic agents in the presence of CYP2D6.10 were higher than those observed in the presence of other CYP2D6 variants. Fluvoxamine, fluoxetine, milnacipran, and haloperidol activated CYP2D6-catalyzed dopamine formation [decreasing the Km and/or increasing the maximal velocity (kcat)], and this activation was CYP2D6 variant-dependent. Regarding the CYP3A subfamily, the effects of various compounds including endogenous steroid hormones on the 6β-hydroxylation of steroid hormones, such as testosterone, progesterone, and cortisol, were determined; it was found that testosterone, dehydroepiandrosterone, and/or α-naphthoflavone activated 6β-hydroxylation of cortisol and/or progesterone, but the effects varied in the presence of different CYP3A subfamily members. Further studies are required to confirm the mechanisms and therapeutic relevance of these activation phenomena.

Role of P34S, G169R, R296C, and S486T Substitutions in Ligand Access and Catalysis for Cytochrome P450 2D6 Allelic Variants CYP2D6*14A and CYP2D6*14B

BACKGROUND

Genetic polymorphism of cytochrome P450 (CYP) contributes to variability in drug metabolism, clearance, and response. This study aimed to investigate the functional and molecular basis for altered ligand binding and catalysis in CYP2D6*14A and CYP2D6*14B, two unique alleles common in the Asian population.

METHODS

CYP proteins expressed in Escherichia coli were studied using the substrate 3-cyano-7- ethoxycoumarin (CEC) and inhibitor probes (quinidine, fluoxetine, paroxetine, terbinafine) in the enzyme assay. Computer modelling was additionally used to create three-dimensional structures of the CYP2D6*14 variants.

RESULTS

Kinetics data indicated significantly reduced intrinsic clearance in CYP2D6*14 variants, suggesting that P34S, G169R, R296C, and S486T substitutions worked cooperatively to alter the conformation of the active site that negatively impacted the deethylase activity of CYP2D6. For the inhibition studies, IC50 values decreased in quinidine, paroxetine, and terbinafine but increased in fluoxetine, suggesting a varied ligand-specific susceptibility to inhibition. Molecular docking further demonstrated the role of P34S and R296C in altering access channel dimensions, thereby affecting ligand access and binding and subsequently resulting in varied inhibition potencies.

CONCLUSION

In summary, the differential selectivity of CYP2D6*14 variants for the ligands (substrate and inhibitor) was governed by the alteration of the active site and access channel architecture induced by the natural mutations found in the alleles.

A Conservative Replacement in the Transmembrane Domain of SARS-CoV-2 ORF7a as a Putative Risk Factor in COVID-19

Simple Summary

The pathogenicity and transmissibility of the COVID-19 pandemic causative agent, the SARS-CoV-2 virus, is related to the functions of the proteins synthesized intracellularly, as guided by viral RNA. It is vitally important to accurately pinpoint novel variants of concern of the SARS-CoV-2 virus, in order to understand the molecular features of novel mutations and manage the on-going battle against the COVID-19 pandemic. We focused on A105V mutation in the ORF7a accessory protein. Sequencing and clinical data showed that this mutation is associated with increased severity and lethality in a group of Romanian patients, despite a lower viral copy number and a lower number of associated comorbidities. This effect is primarily due to increased protein stability through allosteric effects as shown by molecular dynamics analyses. This behavior manifests especially among residues 39–56, and the ones adjacent to 26–30 loop, placed in direct contact with potential interaction partners. Together, the results provide novel insights into the role of ORF7a in the pathogenicity of SARS-CoV-2.

Abstract

The ongoing COVID-19 pandemic follows an unpredictable evolution, driven by both host-related factors such as mobility, vaccination status, and comorbidities and by pathogen-related ones. The pathogenicity of its causative agent, SARS-CoV-2 virus, relates to the functions of the proteins synthesized intracellularly, as guided by viral RNA. These functions are constantly altered through mutations resulting in increased virulence, infectivity, and antibody-evasion abilities. Well-characterized mutations in the spike protein, such as D614G, N439K, Δ69–70, E484K, or N501Y, are currently defining specific variants; however, some less studied mutations outside the spike region, such as p. 3691 in NSP6, p. 9659 in ORF-10, 8782C > T in ORF-1ab, or 28144T > C in ORF-8, have been proposed for altering SARS-CoV-2 virulence and pathogenicity. Therefore, in this study, we focused on A105V mutation of SARS-CoV-2 ORF7a accessory protein, which has been associated with severe COVID-19 clinical manifestation. Molecular dynamics and computational structural analyses revealed that this mutation differentially alters ORF7a dynamics, suggesting a gain-of-function role that may explain its role in the severe form of COVID-19 disease.

Deciphering Structural Alterations Associated with Activity Reductions of Genetic Polymorphisms in Cytochrome P450 2A6 Using Molecular Dynamics Simulations

Cytochrome P450 (CYP) 2A6 is a monooxygenase involved in the metabolism of various endogenous and exogenous chemicals, such as nicotine and therapeutic drugs. The genetic polymorphisms in CYP2A6 are a cause of individual variation in smoking behavior and drug toxicities. The enzymatic activities of the allelic variants of CYP2A6 were analyzed in previous studies. However, the three-dimensional structures of the mutants were not investigated, and the mechanisms underlying activity reduction remain unknown. In this study, to investigate the structural changes involved in the reduction in enzymatic activities, we performed molecular dynamics simulations for ten allelic mutants of CYP2A6. For the calculated wild type structure, no significant structural changes were observed in comparison with the experimental structure. On the other hand, the mutations affected the interaction with heme, substrates, and the redox partner. In CYP2A6.44, a structural change in the substrate access channel was also observed. Those structural effects could explain the alteration of enzymatic activity caused by the mutations. The results of simulations provide useful information regarding the relationship between genotype and phenotype.

Virtual Alanine Scan of the Main Protease Active Site in Severe Acute Respiratory Syndrome Coronavirus 2

Recently, inhibitors of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro) have been proposed as potential therapeutic agents for COVID-19. Studying effects of amino acid mutations in the conformation of drug targets is necessary for anticipating drug resistance. In this study, with the structure of the SARS-CoV-2 Mpro complexed with a non-covalent inhibitor, we performed molecular dynamics (MD) simulations to determine the conformation of the complex when single amino acid residue in the active site is mutated. As a model of amino acid mutation, we constructed mutant proteins with one residue in the active site mutated to alanine. This method is called virtual alanine scan. The results of the MD simulations showed that the conformation and configuration of the ligand was changed for mutants H163A and E166A, although the structure of the whole protein and of the catalytic dyad did not change significantly, suggesting that mutations in His163 and Glu166 may be linked to drug resistance.

KEAP1 Cancer Mutants: A Large-Scale Molecular Dynamics Study of Protein Stability

We have performed 280 μs of unbiased molecular dynamics (MD) simulations to investigate the effects of 12 different cancer mutations on Kelch-like ECH-associated protein 1 (KEAP1) (G333C, G350S, G364C, G379D, R413L, R415G, A427V, G430C, R470C, R470H, R470S and G476R), one of the frequently mutated proteins in lung cancer. The aim was to provide structural insight into the effects of these mutants, including a new class of ANCHOR (additionally NRF2-complexed hypomorph) mutant variants. Our work provides additional insight into the structural dynamics of mutants that could not be analyzed experimentally, painting a more complete picture of their mutagenic effects. Notably, blade-wise analysis of the Kelch domain points to stability as a possible target of cancer in KEAP1. Interestingly, structural analysis of the R470C ANCHOR mutant, the most prevalent missense mutation in KEAP1, revealed no significant change in structural stability or NRF2 binding site dynamics, possibly indicating an covalent modification as this mutant's mode of action.

Molecular dynamic study on PTEN frameshift mutations in breast cancer provide c2 domain as a potential biomarker

PTEN is a tumour suppressor gene known for regulating apoptosis, cell growth, and many other pathways. It is one of the most frequently mutated genes comprising the phosphatase domain (PD) and C terminal domain (C2). Direct therapeutic methods are not applicable for targeting PTEN because once gets mutated, it needs restoration. For mutant detection and restoration using PTEN mRNA there is a need to explore various mutations taking place in PTEN, identify their particular domains, and study their interactions within the cellular system. Here, we have tried to highlight a few such regions in the mutated PTEN of breast cancer patients. In this study, we have selected the top-most-occurring PTEN mutation in breast cancer and compared them to determine the specific properties of each mutation and its effect on functionality. Molecular dynamic simulation for 50 ns was performed on five structures to compare the structural behaviour of mutated PTEN in the system. Our finding suggests that frameshift mutations are more damaging and affect the c2 domain. Frameshift mutant fs_ACTT is the highest occurring as well as the most damaging mutation in all the compared structures. Docking study shows that substitution mutations D92H and R130Q causes loss of binding ability towards PIP2 in normal PTEN, interfering the dephosphorylation process. Overall, the C2 domain is more frequently mutated, and the amino acid residues in the C2 domain show more fluctuations compared to the other regions. Our study can provide the basis for selecting frequently mutated C2 domain as a potential therapeutic marker.Communicated by Ramaswamy H. Sarma.

Structural dynamics of inosine triphosphate pyrophosphatase (ITPA) protein and two clinically relevant mutants: molecular dynamics simulations

The inosine triphosphate pyrophosphatase (ITPA) protein is responsible for removing noncanonical purine nucleoside triphosphates from intracellular nucleotide pools. Absence of ITPA results in genomic instability and increased levels of inosine in DNA and RNA. The proline to threonine substitution at position 32 (P32T) affects roughly 15% of the global population and can modulate treatment outcomes for cancer, lupus, and hepatitis C patients. The substitution of arginine with cysteine at position 178 (R178C) is extremely uncommon and has only been reported in a small cohort of early infantile encephalopathy patients suggesting that a functional ITPA protein is required for life in humans. Here we present molecular dynamic simulations that describe the structure and dynamics of the wild-type ITPA homodimer and two of its clinically relevant mutants, P32T and R178C. The simulation results indicate that both the P32T and R178C mutations alter the structure and dynamic properties of the protein and provide a possible explanation of the experimentally observed effect of the mutations on ITPA activity. Specifically, the mutations increased the overall flexibility of the protein and changed the dominant collective motions of the top lobe as well as the helix 2 of the lower lobe. Moreover, we have identified key active-site residues that are classified as essential or intermediate for inosine triphosphate (ITP) hydrolyzing activity based on their hydrogen bond occupancy. Here we also present biochemical data indicating that the R178C mutant has very low ITP hydrolyzing activity.Communicated by Ramaswamy H. Sarma.

Analysis of the Deleterious Single-Nucleotide Polymorphisms Associated With Antidepressant Efficacy in Major Depressive Disorder

Major depressive disorder (MDD) is a serious mental disease with negative effects on both mental and physical health of the patient. Currently, antidepressants are among the major ways to ease or treat MDD. However, the existing antidepressants have limited efficacy in treating MDD, with a large fraction of patients either responding inadequately or differently to antidepressants during the treatment. Pharmacogenetics studies have found that the genetic features of some genes are associated with the antidepressant efficacy. In order to obtain a better understanding on the relationship between the genetic factors and antidepressant treatment response, we compiled a list of 233 single-nucleotide polymorphisms (SNPs) significantly associated with the antidepressant efficacy in treating MDD. Of the 13 non-synonymous SNPs in the list, three (rs1065852, rs3810651, and rs117986340) may influence the structures and function of the corresponding proteins. Besides, the influence of rs1065852 on the structure of CYP2D6 was further investigated via molecular dynamics simulations. Our results showed that compared to the native CYP2D6 the flexibility of the F-G loop was reduced in the mutant. As a portion of the substrate access channel, the lower flexibility of F-G loop may reduce the ability of the substrates to enter the channel, which may be the reason for the lower enzyme activity of mutant. This study may help us to understand the impact of genetic variation on antidepressant efficacy and provide clues for developing new antidepressants.

Elucidating the Role of Val-Asn 95 and Arg-Gly 52 Mutations on Structure and Stability of Fibroblast Growth Factor Homologous Factor 2

Background:

Fibroblast growth Factor Homologous Factors (FHFs) belong to a subclass of Fibroblast Growth Factor (FGF) family owing to their high sequence and structural similarities with FGFs. However, despite these similarities, there are properties which set them apart from FGFs. FHFs lack the secretion signal sequence unlike other FGF members, except FGF1 and 2. Unlike FGFs, FHFs are not able to bind to FGF Receptors (FGFRs) and instead have been implicated in binding to Voltage-Gated Sodium Channels (VGSCs), neuronal MAP kinase scaffold protein and islet-brain-2 (IB2). The two amino acids Arg-52 and Val95 are conserved in all FHFs and mutation of these residues lead to its inability to bind with VGSC/IB2. However, it is not clear whether the loss of binding is due to destabilization of the protein on mutation or due to involvement of Arg52 and Val95 in conferring functionality to FHFs.

Objective:

In the present study, we have mutated these two conserved residues of FHF2 with its corresponding FGF counterpart amino acids and studied the effects of the mutations on the structure and stability of the protein.

Methods:

Several biophysical methods like isothermal equilibrium denaturation study, ANS fluorescence, intrinsic fluorescence, acrylamide quenching, circular dichroism studies as well as using computational approaches were employed.

Results:

The single mutations were found to affect the overall stability, conformation and functionality of the protein.

Conclusion:

Thus, the studies throw light on the role of specific amino acids in deciding the stability, structure and functionality of proteins and will be useful for development of therapeutically engineered proteins.

Dynamic Effects of CYP2D6 Genetic Variants in a Set of Poor Metaboliser Patients with Infiltrating Ductal Cancer Under Treatment with Tamoxifen

Breast cancer is a group of multigenic diseases. It is the most common cancer diagnosed among women worldwide and is often treated with tamoxifen. Tamoxifen is catalysed by cytochrome P450 2D6 (CYP2D6), and inter-individual variations in the enzyme due to single nucleotide polymorphisms (SNPs) could alter enzyme activity. We evaluated SNPs in patients from Colombia in South America who were receiving tamoxifen treatment for breast cancer. Allelic diversity in the CYP2D6 gene was found in the studied population, with two patients displaying the poor-metaboliser phenotype. Molecular dynamics and trajectory analyses were performed for CYP2D6 from these two patients, comparing it with the common allelic form (CYP2D6*1). Although we found no significant structural change in the protein, its dynamics differ significantly from those of CYP2D6*1, the effect of such differential dynamics resulting in an inefficient enzyme with serious implications for tamoxifen-treated patients, increasing the risk of disease relapse and ineffective treatment.

Microsecond MD simulations of human CYP2D6 wild-type and five allelic variants reveal mechanistic insights on the function

Characterization of cytochrome P450 2D6 (CYP2D6) and the impact of the major identified allelic variants on the activity of one of the most dominating drug-metabolising enzymes is essential to increase drug safety and avoid adverse reactions. Microsecond molecular dynamics simulations have been performed to capture the dynamic signatures of this complex enzyme and five allelic variants with diverse enzymatic activity. In addition to the apo simulations, three substrates (bufuralol, veliparib and tamoxifen) and two inhibitors (prinomastat and quinidine) were included to explore their influence on the structure and dynamical features of the enzyme. Our results indicate that the altered enzyme activity can be attributed to changes in the hydrogen bonding network within the active site, and local structural differences in flexibility, position and shape of the binding pocket. In particular, the increased (CYP2D6*53) or the decreased (CYP2D6*17) activity seems to be related to a change in dynamics of mainly the BC loop due to a modified hydrogen bonding network around this region. In addition, the smallest active site volume was found for CYP2D6*4 (no activity). CYP2D6*2 (normal activity) showed no major differences in dynamic behaviour compared to the wild-type.

Molecular Modeling and Structural Stability of Wild-Type and Mutant CYP51 from Leishmania major: In Vitro and In Silico Analysis of a Laboratory Strain

Cutaneous leishmaniasis is a neglected tropical disease and a major public health in the most countries. Leishmania major is the most common cause of cutaneous leishmaniasis. In the Leishmania parasites, sterol 14α-demethylase (CYP51), which is involved in the biosynthesis of sterols, has been identified as an attractive target for development of new therapeutic agents. In this study, the sequence and structure of CYP51 in a laboratory strain (MRHO/IR/75/ER) of L. major were determined and compared to the wild-type strain. The results showed 19 mutations including seven non-synonymous and 12 synonymous ones in the CYP51 sequence of strain MRHO/IR/75/ER. Importantly, an arginine to lysine substitution at position of 474 resulted in destabilization of CYP51 (ΔΔG = 1.17 kcal/mol) in the laboratory strain; however, when the overall effects of all substitutions were evaluated by 100 ns molecular dynamics simulation, the final structure did not show any significant changes (p-value < 0.05) in stability parameter of the strain MRHO/IR/75/ER compared to the wild-type protein. The energy level for the CYP51 of wild-type and MRHO/IR/75/ER strain were −40,027.1 and −39,706.48 Kcal/mol respectively. The overall Root-mean-square deviation (RMSD) deviation between two proteins was less than 1 Å throughout the simulation and Root-mean-square fluctuation (RMSF) plot also showed no substantial differences between amino acids fluctuation of the both protein. The results also showed that, these mutations were located on the protein periphery that neither interferes with protein folding nor with substrate/inhibitor binding. Therefore, L. major strain MRHO/IR/75/ER is suggested as a suitable laboratory model for studying biological role of CYP51 and inhibitory effects of sterol 14α-demethylase inhibitors.

Evolution of methicillin-resistant Staphylococcus aureus: Evidence of positive selection in a penicillin-binding protein (PBP) 2a coding gene mecA

Methicillin-resistant Staphylococcus aureus (S. aureus) (MRSA) represents more and more S. aureus infections. MecA, the novel coding gene of penicillin-binding protein (PBP) 2a of MRSA, is the key resistance factor of β-lactam, but little is known about the evolution of this gene. Given the crucial role of mecA in S. aureus physiology and β-lactam resistance, the selective forces may contribute to adaptation of the bacteria to the special environments such as its host or antibiotics. To understand the evolution of this gene, we screened GenBank database and analyzed mecA of 249 S. aureus strains. Twenty-nine unique alleles with 26 unique amino acid sequences were identified. Phylogenetic analysis showed three main groups of mecA in the global S. aureus strains. Analysis of these alleles using codon-substitution models (M8, M3, and M2a) and likelihood ratio tests (LRTs) of the codeML package and a random-effects likelihood (REL) method of HyPhy package for the site-specific ratio of nonsynonymous to synonymous substitution rates suggested that fourteen sites in the allosteric domain of PBP2a have been subjected to strong positive selection pressure. Mutations of two positive selection sites (N146K and E239K) were reported to be essential for ceftaroline- or L-695, 256-resistant. Further study indicated that the positive selection pressure might be more likely related to the host's inflammatory or immune response during S. aureus infection. Our studies provide the first evidence of positive Darwinian selection in the mecA of S. aureus, contributing to a better understanding of the adaptive mechanism of this bacterium.

Implementation of Pharmacogenomics in Everyday Clinical Settings

Currently, germline pharmacogenomics (PGx) is successfully implemented within certain specialties in clinical care. With the integration of PGx in pharmacotherapy multiple stakeholders are involved, which are identified in this chapter. Clinically relevant pharmacogenes with their related PGx test are discussed, along with diagnostic test criteria to guide clinicians and policy makers in PGx test selection. The chapter further reviews the similarities and the differences between the guidelines of the Dutch Pharmacogenetics Working Group and the Clinical Pharmacogenetics Implementation Consortium which both support healthcare professionals in understanding PGx test results and help guiding pharmacotherapy by providing evidence-based dosing recommendations. Finally, clinical studies which provide scientific evidence and information on cost-effectiveness supporting clinical implementation of PGx in clinical care are discussed along with the remaining barriers for adoption of PGx testing by healthcare professionals.

Investigation of substrate recognition for cytochrome P450 1A2 mediated by water molecules using docking and molecular dynamics simulations

The role of water molecules in the active site of cytochrome P450 1A2 (CYP1A2) was investigated using an explicit water model to simulate biological environments. Moreover, differences in ligand recognition between the inhibitor α-naphthoflavone (ANF) and the substrate 7-ethoxyresorufin (7ER) in the CYP1A2 complex were examined. More than 200-ns molecular dynamics (MD) simulations were performed for each complex structure of CYP1A2. In the complex structure with 7ER obtained after MD simulation, some water molecules existed in the active site and formed hydrogen bonds between 7ER and some residues. However, in the complex structure with ANF, the hydrogen bond network differed. These results suggest that CYP1A2 requires water molecules in its active site for substrate recognition. The observed differences in the hydrogen bond network in the complex with ANF or 7ER may be due to the fact that ANF is an inhibitor.

Drug Metabolism in Preclinical Drug Development: A Survey of the Discovery Process, Toxicology, and Computational Tools

BACKGROUND

While establishing efficacy in translational models and humans through clinically-relevant endpoints for disease is of great interest, assessing the potential toxicity of a putative therapeutic drug is critical. Toxicological assessments in the pre-clinical discovery phase help to avoid future failure in the clinical phases of drug development. Many in vitro assays exist to aid in modular toxicological assessment, such as hepatotoxicity and genotoxicity. While these methods have provided tremendous insight into human toxicity by investigational new drugs, they are expensive, require substantial resources, and do not account for pharmacogenomics as well as critical ADME properties. Computational tools can fill this niche in toxicology if in silico models are accurate in relating drug molecular properties to toxicological endpoints as well as reliable in predicting important drug-target interactions that mediate known adverse events or adverse outcome pathways (AOPs).

METHODS

We undertook an unstructured search of multiple bibliographic databases for peer-reviewed literature regarding computational methods in predictive toxicology for in silico drug discovery. As this review paper is meant to serve as a survey of available methods for the interested reader, no focused criteria were applied. Literature chosen was based on the writers' expertise and intent in communicating important aspects of in silico toxicology to the interested reader.

CONCLUSION

This review provides a purview of computational methods of pre-clinical toxicologic assessments for novel small molecule drugs that may be of use for novice and experienced investigators as well as academic and commercial drug discovery entities.