# Analysis of protein missense alterations by combining sequence- and structure-based methods

Aram Gyulkhandanyan, Alireza R. Rezaie, Lubka Roumenina, Nathalie Lagarde, Veronique Fremeaux-bacchi, Maria A. Miteva, Bruno O. Villoutreix

► To cite this version:

Aram Gyulkhandanyan, Alireza R. Rezaie, Lubka Roumenina, Nathalie Lagarde, Veronique Fremeaux-bacchi, et al.. Analysis of protein missense alterations by combining sequence- and structure-based methods. Molecular Genetics & Genomic Medicine, 2020, 8 (4), pp.e1166. 10.1002/mgg3.1166. hal-03246847

HAL Id: hal-03246847

https://hal.science/hal-03246847

Submitted on 21 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Analysis of protein missense alterations by combining sequence- and structure-based methods

Aram Gyulkhandanyan<fnref n="1" />,<fnref n="2" />

Alireza R. Rezaie<fnref n="3" />,<fnref n="4" />

Lubka Roumenina<fnref n="5" />,<fnref n="6" />,<fnref n="7" />

Nathalie Lagarde<fnref n="1" />,<fnref n="8" />

Veronique Fremeaux-Bacchi<fnref n="5" />,<fnref n="6" />,<fnref n="7" />,<fnref n="9" />

Maria A. Miteva<fnref n="1" />,<fnref n="10" />

Bruno O. Villoutreix<fnref n="1" />,<fnref n="11" />

<figure><img src="image_2.png" /></figure>

## Correspondence

Bruno O. Villoutreix, INSERM U973,
Laboratory MTi, University Paris Diderot,
Paris, France.

Email: bruno.villoutreix@inserm.fr

## Funding information

This work was funded by the ANR
COMPC3, ANR ToxME, I-Site Lille and
Heart, Lung, and Blood Institute of the
National Institutes of Health HL062565.

## Abstract

**Background:** Different types of in silico approaches can be used to predict the phe-notypic consequence of missense variants. Such algorithms are often categorized as sequence based or structure based, when they necessitate 3D structural information. In addition, many other in silico tools, not dedicated to the analysis of variants, can be used to gain additional insights about the possible mechanisms at play.

**Methods:** Here we applied different computational approaches to a set of 20 known missense variants present on different proteins (CYP, complement factor B, antithrombin and blood coagulation factor VIII). The tools that were used include fast computational approaches and web servers such as PolyPhen-2, PopMusic, DUET, MaestroWeb, SAAFEC, Missense3D, VarSite, FlexPred, PredyFlexy, Clustal Omega, meta-PPISP, FTMap, ClusPro, pyDock, PPM, RING, Cytoscape, and ChannelsDB.

**Results:** We observe some conflicting results among the methods but, most of the time, the combination of several engines helped to clarify the potential impacts of the amino acid substitutions.

**Conclusion:** Combining different computational approaches including some that were not developed to investigate missense variants help to predict the possible

<footnote n="1">¹INSERM U973, Laboratory MTi, University Paris Diderot, Paris, France</footnote>

<footnote n="2">²Laboratory SABNP, University of Evry, INSERM U1204, Université Paris-Saclay, Evry, France</footnote>

<footnote n="3">³Cardiovascular Biology Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA</footnote>

<footnote n="4">⁴Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK, USA</footnote>

<footnote n="5">⁵INSERM, UMR_S 1138, Centre de Recherche des Cordeliers, Paris, France</footnote>

<footnote n="6">⁶Sorbonne Universités, Paris, France</footnote>

<footnote n="7">⁷Université Paris Descartes, Sorbonne Paris Cité, Paris, France</footnote>

<footnote n="8">⁸Laboratoire GBCM, EA7528, Conservatoire national des arts et métiers, Hesam Université, Paris, France</footnote>

<footnote n="9">⁹Assistance Publique-Hôpitaux de Paris, Service d'Immunologie Biologique, Hôpital Européen Georges Pompidou, Paris, France</footnote>

<footnote n="10">¹⁰Inserm U1268 MCTR, CNRS UMR 8038 CITCoM, Faculté de Pharmacie de Paris, Univ. de Paris, Paris, France</footnote>

<footnote n="11">¹¹INSERM, Institut Pasteur de Lille, U1177-Drugs and Molecules for Living Systems, Université de Lille, Lille, France</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.</footnote>

<footnote>© 2020 The Authors. *Molecular Genetics & Genomic Medicine* published by Wiley Periodicals, Inc.</footnote>
impact of the amino acid substitutions. Yet, when the modified residues are involved in a salt-bridge, the tools tend to fail, even when the analysis is performed in 3D. Thus, interactive structural analysis with molecular graphics packages such as Chimera or PyMol or others are still needed to clarify automatic prediction.

KEYWORDS

Antithrombin, CYP, Factor B, Factor VIII, missense variants, PolyPhen-2, structural analysis, structural bioinformatics

# 1 | INTRODUCTION

Analysis of human genetic variations and its relationship to disease and drug response has gained remarkable attention these recent years. Genome-wide association studies and candidate gene association studies tend to associate single nucleotide alterations with diseases. While these approaches are used at the DNA level, it is important to analyze, when possible, variations at the protein structural level as well. Structural analysis does indeed assist the development of rational hypotheses about possible impacts of substitutions and their possible links to disease states. Furthermore, if the protein under investigation is a therapeutic target, amino acid substitutions can cause drastic changes in drug target phenotypes, thereby resulting in dysfunctional drugs. Such investigations are therefore important in the field of drug discovery.

It is known that amino acid changes can affect both the function and the structure of a protein (i.e., alter catalysis, induce posttranslational modification, folding and stability, perturb ligand binding, favor multimerization, etc). Numerous studies have established the importance of combining experimental and in silico approaches for studying the influence of amino acid substitutions on the structure and function of proteins (Martiny & Miteva, 2013; Ramamoorthy & Skaar, 2011; Singh, Kashyap, Pandey, & Saini, 2011; Takano et al., 2012; Villoutreix, 2002; Witham, Takano, Schwartz, & Alexov, 2011). Nevertheless, it is not feasible to take an experimental approach for studying all the amino acid substitutions identified in patients' proteins. To assist the process and to gain knowledge about the possible impact of amino acid substitutions with the aim of reducing wet laboratory experiments, numerous bioinformatics methods have been developed (Kucukkal, Petukh, Li, & Alexov, 2015; Thusberg & Vihinen, 2009). These tools are usually categorized as sequence based, structure based, or involve the combination of both. The sequence-based tools, in general, provide important information about conserved residues and can give some insights about how amino acid changes impact a protein. On the other hand, the second set of methods uses the 3D structures (experimental or high-quality homology models) of the proteins, and the consensus belief/observation is that analyses performed at the structural levels are more reliable (Thusberg

& Vihinen, 2009). Indeed, when 3D structures of proteins are known, several characteristics can be analyzed including stability and/or the dynamics of the mutant proteins and the impact of amino acid changes in different types of molecular interactions can be analyzed (Kucukkal et al., 2015). Yet, if numerous amino acid changes have to be analyzed, it is important to select fast 3D computational approaches. Overall, it is still difficult at present to decide which tools can reliably predict the pathogenic character of new variants (Grimm et al., 2015).

In this study, we investigated 20 amino acid substitutions present in 6 different proteins involved in the health and disease states (Figure 1): (a) three members of the Cytochrome P450 (CYP) family (drug-metabolizing enzymes responsible for the metabolism of most human drugs); (b) a protein from the complement system (a system involved in immune response and intended to protect the body from foreign agents); and (c) two proteins from the blood coagulation system.

The cytochrome P450 superfamily is responsible for the metabolism of about 90% of the common therapeutic drugs (Fujikura, Ingelman-Sundberg, & Lauschke, 2015; Singh et al., 2011). Among human CYPs, the isoforms CYP1A2, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4 are responsible for the oxidative catalysis of most human drugs and many exhibit genetic polymorphisms (Fujikura et al., 2015; Singh et al., 2011). The genetic variations resulting in missense alterations can lead to loss-of-function variants or gain-of-function variants of CYP and can cause toxicity or lack of efficacy. Several CYP isoforms of the class 1 family are highly polymorphic (Daly, 2015; Martiny & Miteva, 2013; Zhou, Ingelman-Sundberg, & Lauschke, 2017). Here, we have investigated missense alterations on CYP2B6, CYP2C9, and CYP1A2 proteins (Fujikura et al., 2015). The impact of missense alterations on the state of various CYP isoforms and their effect on drug metabolism has already been conducted by several research groups (Elfaki, Mir, Almutairi, & Duhier, 2018; Isvoran et al., 2017; Louet et al., 2018; Rydberg & Olsen, 2012; Simonetti et al., 2018) but in this article, we focused on a total of 10 missense alterations present on the CYP2B6, CYP2C9, and CYP1A2 isoforms. Of particular importance for the present study, CYP proteins have internal channels
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** 3D structures of the investigated proteins. The proteins are shown as cartoon and the substituted amino acids are highlighted in orange. A-CYP 1A2 (Heme in red), B-CYP 2B6 (Heme in red), C-CYP 2C9 (Heme in red), D-Complement Factor B (serine protease [SP] domain in blue and the von Willebrand factor [VWF] domain in green). E-The factor VIII (FVIII) protein structure with the C2 domain colored in blue. F-Antithrombin (AT) with the reactive center loop colored in red</figcaption></figure>

that play critical structural roles in the function of these proteins. These channels or tunnels tend to connect internal spaces of the proteins with the exterior surface, thereby enabling substrate/product transport toward the catalytic site. Some of these channels have been characterized, but there are still others that are not known.

The complement system proteins are key players in the innate immunity with over 20 protein components (Merle, Church, Church, Fremeaux-Bacchi, & Roumenina, 2015; Merle, Noe, Noe, Halbwachs-Mecarelli, Fremeaux-Bacchi, & Roumenina, 2015). Factor B (FB) is a major component of the alternative pathway of complement activation. FB associates with the central complement component C3b to form the C3-cleaving enzyme C3 convertase. The activity of this enzyme leads to the generation of novel derivatives C3, allowing the opsonization of pathogens, the generation of pro-inflammatory anaphylatoxine C3a, and promoting

the terminal membrane lytic complex of the complement pathway. Crystal structures of FB, of its fragments and domains (Ba, Bb, vWF-A domain, and SP-domain) and of the C3 convertase (C3bBb) have been reported, allowing detailed structural analyses of its functional properties (Bhattacharya, Lupher, Staunton, & Liddington, 2004; Forneris et al., 2010; Jing et al., 2000; Milder et al., 2007; Pedersen et al., 2017; Ponnuraj et al., 2004; Rooijakkers et al., 2009). Missense alterations causing changes in the FB function can have severe consequences, and in particular, lead to atypical hemolytic uremic syndrome (aHUS). aHUS is a rare genetic renal thrombotic microangiopathy disease, associated with genetic abnormalities and the over activation of the alternative complement pathway, causing a glomerular endothelial cells damage. Mutations in FB have been described in about 1% of the aHUS cases and extensive research toward understanding the functional
consequences of specific mutations have been reported (Marinozzi et al., 2014). In the present study we investigated four selected amino acid substitutions in FB.

Blood coagulation comprises of more than 30 proteins that interact with each other with a high degree of specificity (Cardenas, Rein-Smith, & Church, 2016; Palta, Saroa, & Palta, 2014; Versteeg, Heemskerk, Levi, & Reitsma, 2013; Villoutreix, 2002). Missense mutations identified in blood coagulation proteins can lead to life-threatening illnesses. For example, amino acid changes in distinct regions of factor VIII (FVIII) can cause bleeding disorders (Hemophilia A) (Liu et al., 2000; Pratt et al., 1999; Spiegel, Jacquemin, Saint-Remy, Stoddard, & Pratt, 2001). Of particular interest for the present study, some of the coagulation proteins including FVIII function properly only when they are anchored in an appropriate membrane surface. Another example is Antithrombin (AT). This plasma inhibitory protein is a key regulator of the coagulation system (Huntington, Olson, Fan, & Gettins, 1996; Izaguirre et al., 2014; Muszbek, Bereczky, Kovács, & Komáromi, 2010; Olson, Richard, Izaguirre, Schedin-Weiss, & Gettins, 2010). This protein acts as an essential inhibitor of activated factor X (Xa) and thrombin (IIa), and a number of other activated coagulation factors (Abildgaard, 2007; Huntington et al., 1996; Izaguirre et al., 2014; Muszbek et al., 2010; Olson et al., 2010). AT becomes an effective inhibitor only upon binding to heparin or heparan sulfate proteoglycans lining the vascular endothelium (Cooper, Coath, Daly, & Makris, 2011; Huntington et al., 1996). Numerous alterations have been identified in the antithrombin gene, which are mainly missense alterations (Bayton & Lane, 22003; Luxembourg et al., 2011; Stenson et al., 2014). In the present study, we investigated two amino acid substitutions in FVIII and four in antithrombin.

The goal of this article is to illustrate how 2D and 3D in silico analysis of amino acid changes can complement each other in developing hypotheses about the possible impacts of the substitutions on the structure and function of the investigated proteins and also highlight their limitations. Furthermore, we suggest that it can be valuable to also use in silico methods not dedicated to the prediction of missense alterations. In addition, we wanted to investigate the benefit of using automatic 3D mapping online services and if 3D interactive analysis of the amino acid substitution using molecular graphics packages could provide additional insights. We have selected 20 illustrative amino acid changes with associated experimental/clinical data. Clearly, there are numerous in silico tools that can be used to explore the possible impact of amino acid changes on the structure and function of a protein (Hu et al., 2019; Villoutreix, Lagorce, Labbé, Sperandio, & Miteva, 2013), and here we have used a few diverse, fast, relatively recent and user-friendly sequence-based and structure-based computational approaches (see Table 1).

## 2 | METHODS

Ethical Compliance. The manuscript is a retrospective case report that does not require ethics committee approval.

Tools and key databases that were used in this study are reported in Table 1.

### 2.1 | Protein structures and variant databases

The following experimental structures were downloaded from the Protein Data Bank (PDB) (Berman et al., 2000) and used in our study: CYP 1A2 (PDB ID: 2HI4) (Sansen et al., 2007), CYP 2B6 (PDB ID: 3IBD) (Shah et al., 2010), CYP 2C9 (PDB ID: 1OG2) (Williams et al., 2003), complement factor B (PDB ID: 2OK5) (Milder et al., 2007), coagulation factor VIII (PDB ID: 2R7E) (Shen et al., 2008), and antithrombin (PDB ID: 2BEH) (Johnson et al., 2006). Heteroatoms were removed except for the heme group of the CYP proteins. Earlier, the nomenclature and properties of allelic variants of human CYP could be found at http://www.cypalleles.ki.se (Sim & Ingelman-Sundberg, 2013). However, the Human Cytochrome P450 Allele Nomenclature Database can now be found at the Pharmacogene Variation (PharmVar) Consortium (Gaedigk et al., 2018). Thus, all investigated amino acid changes for the three types of CYP proteins were taken from the PharmVar updated list. Some missing side chains or engineered residues present in the experimental structures were added or changed back to the wild-type sequence with the Chimera software (Pettersen et al., 2004). The PyMOL Molecular Graphics System (Version 1.8.2.2 Schrödinger, LLC) and Chimera were both used for the structural analysis.

### 2.2 | Computational approaches and databases

#### 2.2.1 | PolyPhen-2

We selected PolyPhen-2 (Polymorphism Phenotyping) (Adzhubei et al., 2010) to investigate the amino acid changes because this method is fast, user-friendly and used in many laboratories. The tool is mainly sequence based but it also incorporates some structural information when available to produce qualitative predictions. In the current 2.2 version, the mutations are labeled as "benign," "possibly damaging," and as "probably damaging."

#### 2.2.2 | Protein stability prediction in 3D

For the assessment of the protein stability, we used four different programs: PopMusic, DUET, MaestroWeb, and SAAFEC.
<table><caption>**TABLE 1** In silico tools and databases used in the present study</caption><thead><tr><td>Tool</td><td>URL</td><td>Input</td></tr></thead><tbody><tr><td>PolyPhen-2 (investigation of variants)</td><td>http://genetics.bwh.harvard.edu/pph2</td><td>Sequence</td></tr><tr><td>PopMusic (ΔΔG)</td><td>http://dezyme.com</td><td>Uploaded 3D structure</td></tr><tr><td>DUET (ΔΔG)</td><td>http://biosig.unimelb.edu.au/duet</td><td>Uploaded 3D structure</td></tr><tr><td>MaestroWeb (ΔΔG)</td><td>https://biwww.che.sbg.ac.at/maestro/web</td><td>Uploaded 3D structure</td></tr><tr><td>SAAFEC (ΔΔG)</td><td>http://compbio.clemson.edu/SAAFEC</td><td>Uploaded 3D structure</td></tr><tr><td>Missense3D</td><td>http://www.sbg.bio.ic.ac.uk/~missense3d</td><td>Protein sequence or 3D structure/ PDB ID</td></tr><tr><td>VarSite</td><td>https://www.ebi.ac.uk/thornton-srv/databases/ VarSite</td><td>UniProt ID, or search terms or disease</td></tr><tr><td>FlexPred (flexibility prediction)</td><td>http://kiharalab.org/flexPred</td><td>Uploaded 3D structure</td></tr><tr><td>PredyFlexy (flexibility prediction)</td><td>http://www.dsimb.inserm.fr/dsimb_tools/predyflexy</td><td>Sequence</td></tr><tr><td>Clustal Omega (multiple sequence alignment)</td><td>https://www.ebi.ac.uk/Tools/msa</td><td>Sequence</td></tr><tr><td>meta-PPISP (prediction of PPI sites)</td><td>http://pipe.scs.fsu.edu/meta-ppisp</td><td>Uploaded 3D structure</td></tr><tr><td>FTMap (binding pocket prediction)</td><td>http://ftmap.bu.edu/home.php</td><td>Uploaded 3D structure</td></tr><tr><td>ClusPro (macromolecular docking)</td><td>https://cluspro.org</td><td>Uploaded 3D structure</td></tr><tr><td>PPM (protein membrane interaction prediction)</td><td>http://opm.phar.umich.edu</td><td>Uploaded 3D structure</td></tr><tr><td>pyDockWeb (macromolecular docking)</td><td>https://life.bsc.es/pid/pydock</td><td>Uploaded 3D structure</td></tr><tr><td>RING-2.0 (residue interaction network)</td><td>http://protein.bio.unipd.it/ring</td><td>Uploaded 3D structure</td></tr><tr><td>Cytoscape (data visualization)</td><td>http://www.cytoscape.org</td><td>RING output XLM file</td></tr><tr><td>ChannelsDB (investigate channels)</td><td>https://webchemdev.ncbr.muni.cz/ChannelsDB</td><td>PDB ID</td></tr><tr><td>PharmVar (variation database)</td><td>https://www.pharmvar.org/htdocs/archive/index_ original.htm</td><td>—</td></tr><tr><td>gnomAD (Genome Aggregation Database)</td><td>https://gnomad.broadinstitute.org/</td><td>—</td></tr><tr><td>ClinVar (genomic variation database)</td><td>https://www.ncbi.nlm.nih.gov/clinvar/</td><td>—</td></tr><tr><td>Chimera (molecular graphics)</td><td>http://www.cgl.ucsf.edu/chimera</td><td>—</td></tr><tr><td>PyMol (molecular graphics)</td><td>http://www.pymol.org</td><td>—</td></tr><tr><td>PDB (protein structure database)</td><td>https://www.rcsb.org</td><td>—</td></tr><tr><td>UniProt (annotated protein database)</td><td>http://www.uniprot.org</td><td>—</td></tr></tbody></table>

These web servers compute ΔΔG values between the wild-type and the variant proteins. *PoPMuSiC* (Gonnelli, Rooman, & Dehouck, 2012) is a software which evaluates the changes in the stability of a given protein upon amino acid changes. The server predicts the thermodynamic stability changes caused by a single site substitution using a linear combination of statistical potentials whose coefficients depend on the solvent accessibility of the modified residue (Dehouck, Kwasigroch, Gilis, & Rooman, 2011). *DUET* is a software that evaluates protein stability with an optimized predictor that makes use of support vector machine approaches (Pires, Ascher, & Blundell, 2014aa). *DUET* consolidates two complementary approaches SDM (Pandurangan, Ochoa-Montaño, Ascher, & Blundell, 2017) and mCSM (Pires, Ascher, & Blundell, 2014bb) in a consensus prediction. This is achieved by combining the results of the two separate methods with an optimized predictor. *MAESTRO* is a method for predicting changes in stability. It is a structure-based approach that provides predicted free energy change (ΔΔG) values as well as corresponding confidence estimation values for the predictions

while at the same time allowing for high-throughput scanning of multi-point amino acid changes (Laimer, Hiebl-Flach, Lengauer, & Lackner, 2016; Laimer, Hofer, Fritz, Wegenkittl, & Lackner, 2015). The *SAAFEC* (Single Amino Acid Folding Free Energy Changes) method (Getov, Petukh, & Alexov, 2016) is designed for calculating the folding free energy changes caused by missense alterations. Based on the MM-PBSA method (Homeyer & Gohlke, 2012) with weight coefficients, this approach was optimized using experimental data from the ProTherm database (Bava, Gromiha, Uedaira, Kitajima, & Sarai, 2004). With this approach, the protein structures undergo an energy minimization using the NAMD software (Phillips et al., 2005).

### 2.2.3 | Automatic 3D online structural mapping of missense variants

Two web servers were used: Missense3D (Ittisoponpisan et al., 2019) and VarSite (Laskowski, Stephenson, Sillitoe,
Orengo, & Thornton, 2020). The *Missense3D* pipeline maps and analyze amino acid changes on experimental and homology model protein 3D structures. *VarSite* maps known disease-associated variants from UniProt (UniProt Consortium, 2019), ClinVar (Landrum et al., 2014), and gnomAD (genome aggregation database) (Karczewski et al., 2019) or data provided by the users onto protein experimental 3D structures. A disease propensity score is also reported, the value quantifies how much more often a variant is observed in diseases than in the natural variant data obtained from gnomAD. The value ranges from very low (propensity = 0.25) to very high (propensity = 3.27). On both servers, users obtain a report card with information about the amino acid substitution.

## 2.2.4 | Multiple sequence alignment

To investigate sequence conservation for the abovementioned proteins, multiple sequence alignments (MSA) were performed with the EMBL-EBI *Clustal Omega* server (Sievers et al., 2011) using as input sequences from different species downloaded from the UniProtKB database (UniProt Consortium, 2019).

## 2.2.5 | Protein flexibility

Some regions of proteins can be moderately to highly flexible. Flexibility can be inferred in some cases from X-ray experiments, obtained from NMR studies or explored using long molecular simulation approaches. Yet, some very fast approaches have been reported to provide relatively accurate predictions without a need for CPU/GPU intensive calculations. We here employed the predicted B-factor (relative vibrational motion) and RMSFs (root-mean-square fluctuations) obtained from the prediction program *PredyFlexy* (fast computations carried out only over the protein sequences) (de Brevern, Bornot, Craveur, Etchebest, & Gelly, 2012). Three types of flexibility are proposed by this approach. PredyFlexy classifies amino acid residues into rigid, intermediate or flexible sites. We also used the program *FlexPred*, a fast method that uses the protein 3D structure as input and predicting fluctuations using easily computed static structural features (Kuznetsov & McDuffie, 2008). Overall, the tool determines which amino acid residues are located in flexible sites or in more rigid regions (Jamroz, Kolinski, & Kihara, 2012; Kuznetsov & McDuffie, 2008).

## 2.2.6 | Residue interaction network and visualization

The Residue Interaction Network Generator (RING) software was used to gain additional insights into the structures of the

selected proteins through visualization of nonbonded interactions (Martin et al., 2011). The *RING-2.0* server was used in our study (Piovesan, Minervini, & Tosatto, 2016). The generated RING network files were analyzed with *Cytoscape*, a tool that provides a basic set of features for data integration, analysis, and visualization (Shannon et al., 2003).

## 2.2.7 | Predictions of ligand-binding pockets

To investigate binding pockets and ligand-binding hot spots we used the *FTMap* fragment mapping server (Brenke et al., 2009; Kozakov et al., 2015, 2011; Ngan et al., 2012). FTMap positions 16 small organic probe molecules of varying size, shape, and polarity over the protein surface to predict binding pockets and key residues involved. This approach can also provide insights into protein-protein interface regions.

## 2.2.8 | Protein-Protein interaction site prediction

To predict protein-protein interaction sites, we used the *meta-PPISP* web server (Qin & Zhou, 2007). A number of complementary methods have been developed for predicting protein-protein interaction sites. By combining results from different predictors, the meta-PPISP method was found to improve prediction robustness and accuracy. Its operating system is built on three different web servers: cons-PPISP (Chen & Zhou, 2005; Zhou & Shan, 2001), PINUP (Liang, Zhang, Liu, & Zhou, 2006), and Promate (Neuvirth, Raz, & Schreiber, 2004).

## 2.2.9 | Channel-related analyses for the CYP proteins

For these investigations we used the *ChannelsDB* database (Pravda et al., 2018). Channels, tunnels, and pores are very important structural features within biomacromolecules. Tunnels connect internal spaces with exterior regions. Enzyme active sites can be connected to the exterior environment by one or more channels passing through the protein (Louet et al., 2018; Pravda et al., 2014). Here, CYP proteins were investigated with such tools.

## 2.2.10 | Antibody-protein and heparin-protein docking

Protein-protein and heparin-protein docking computations were carried out with the *ClusPro* server. ClusPro offers a number of options for docking experiments in particular modes tuned for antibody and heparin molecules (Kozakov
et al., 2017; Mottarella et al., 2014). To study the interaction of an antibody with the FVIII C2 domain the ClusPro server was used (Brenke et al., 2012). Additional docking experiments between a Fab fragment and FVIII C2 domain were performed with the *pyDockWeb* server (Jiménez-García, Pons, & Fernandez-Recio, 2013). Docking of heparin on antithrombin was performed using the advanced options of the ClusPro server (Kozakov et al., 2017; Mottarella et al., 2014).

## 2.2.11 | Protein-membrane interaction prediction

To position FVIII at the surface of a membrane, we used the *PPM* web server (Lomize, Pogozheva, Joo, Mosberg, & Lomize, 2012). The location and orientation of the protein are obtained by various rotations and translations with the goal of optimizing the protein transfer energy from water to a virtual lipid bilayer.

# 3 | RESULTS AND DISCUSSION

We investigated several CYP450, complement, and blood coagulation proteins. To analyze the potential impacts of the selected amino acid changes, we first used PolyPhen-2 and then additional sequence- and structure-based approaches associated with other structural bioinformatics tools not developed specifically to investigate point mutations and with interactive 3D structural analysis. Also, two automatic 3D mapping tools (Missense3D, VarSite) were used to gather additional information about the possible impact of the amino acid substitutions and to compare the outputs with the 3D interactive investigations carried out with standalone molecular graphics engine PyMol and Chimera. Our analysis focuses on highlighting the differences and complementarity of the in silico methods as well as limitations. Warnings are reported based on our observations together with some recommendations.

The mainly sequence-based approach PolyPhen-2 was used with other structure-based approaches (DUET, PopMusic, SAAFEC, MAESTROweb). Different computations not fully dedicated to the study of amino acid substitutions were also used. Residue interaction network computation is used to gain knowledge about protein structures and functions. In this context, amino acid residues are referred to as nodes while edges represent noncovalent interactions. Structural features can be reported for each node such as degree, here the number of noncovalent interactions with surrounding amino acids or cofactors. Intuitively, a substitution that involves a highly connected residue is likely to perturb the structure and/or the function of the protein. We used the RING server (Martin et al., 2011; Piovesan et al., 2016) to carry out such investigation followed by visualization performed with Cytoscape (Shannon et

al., 2003). It is also important to estimate whether the amino acid change is located in a very rigid area or in a flexible region since this can change the dynamics of the system. When dealing with several substitutions, obviously the use of fast computational approaches is recommended. For these reasons we decided to use the structure-based FlexPred approach (Kuznetsov & McDuffie, 2008) as well as the sequence-based PredyFlexy tool (de Brevern et al., 2012). Damaging amino acid changes are often located in the core interior of proteins, but it is also known that when they are solvent exposed, damaging amino acid substitutions tend to be located in protein-protein interaction sites (Yates & Sternberg, 2013). We used the meta-PPISP web server (Qin & Zhou, 2007) to predict such regions. Other types of information can be investigated in silico, depending on the molecular functions that need to be addressed, via for instance, protein-protein docking or protein-heparin docking experiments, prediction of putative ligand-binding sites or of channels, or predictions of zones potentially important for membrane interactions. Of importance, we note that many laboratories working with amino acid substitutions worldwide make little use of 3D structural data, we investigated the benefit of using automatic 3D mapping tools such as VarSite and Missense3D.

## 3.1 | CYP1A2

For CYP1A2 (PDB ID: 2HI4, UniProt: P05177), we investigated the possible molecular effects of two clinically important alterations (Asp348Asn and Ile386Phe) (Table 2). Human cytochrome P450 1A2 catalyzes important reactions during the metabolism of xenobiotics, including N-hydroxylation of carcinogenic aromatic amines. Variants Asp348Asn and Ile386Phe were expressed at levels less than half of the wild-type proteins and these differences (vs. wild-type) were statistically significant ($p < 0.05$) (Zhou, Josephy, Kim, & Guengerich, 2004). Table 2 shows results of the in silico analyses for CYP1A2 (PDB ID: 2HI4). CYP1A2 is not processed yet on the ChannelsDB server. Yet, from our previous work, it seems possible that Ile386 is located in one of the substrate-binding channel (see Figure 6, channel 2f in Louet et al., 2018). Potential ligand-binding pockets were also investigated with FTMap (Brenke et al., 2009; Kozakov et al., 2015, 2011; Ngan et al., 2012) and 14 binding pockets were suggested (Table S1). Only the amino acid Ile386 is predicted to be located in a ligand-binding pocket (Table 2).

### 3.1.1 | Asp348Asn (allele CYP1A2*3)

The PolyPhen-2 result suggests that this change is benign but three structure-based tools indicate that the substitution is moderately destabilizing (specially SAAFEC) while
<table><caption>TABLE 2 In silico analysis of CYP1A2 (PDB ID: 2H14, UniProt: P05177; GenBank: NM_000761.5)</caption><thead><tr><th colspan="2">CYP1A2</th><th colspan="2">PolyPhen-2</th><th colspan="2">Involved in a predicted ligand-binding pocket</th><th colspan="2">Involved in known or predicted channel</th><th colspan="2">Involved in predicted PPI sites (meta-PPISP)</th><th colspan="2">Predicted fluctuation value (RING-2.0)</th><th colspan="2">Flexibility class<sup>d</sup> (PredyFlexy)</th><th colspan="2">Decreased metabolic activity or protein expression</th><th rowspan="2">References</th></tr><tr><th>Amino acid/ mutation/ alleles</th><th>DUET</th><th>PopMusic<sup>a</sup></th><th>SAAFEC</th><th>ΔΔG, kcal/mol</th><th>ΔΔG, kcal/mol</th><th>MSA<sup>c</sup></th><th>aa</th><th>MSA<sup>c</sup></th><th>aa</th><th>Node degree (RING-2.0)</th><th>FlexPred</th><th>0</th><th>D348N was expressed at levels less than half of the wild type (in vitro)</th><th>Zhou et al. (2004)</th></tr><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td><b>Asp348</b></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><b>D348N</b></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><b>CYP1A2*3</b></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><b>Ile386</b></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><b>I386F</b></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><b>CYP1A2*4</b></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

MaestroWeb proposes a small stabilization effect. Asp348 is strictly conserved in our multiple sequence alignment (Figure S1) and partially solvent exposed (SA = 48.03%), located on a helical structure (Figure 1) and predicted to be in a relatively rigid area of the protein. It forms a salt-bridge with Arg353 that will be lost upon its replacement with an asparagine (data not shown). Possibly, this substitution destabilizes this region of the protein. This residue makes several additional noncovalent interactions with its surrounding (node degree = 8), it is not thought to be part of a ligand-binding pocket or a channel (Louet et al., 2018), and it is not predicted to be involved in a protein-protein interaction site (Table 2). Overall, three of the four 3D web servers suggest that this substitution could be problematic, in overall agreement with our 3D interactive structural analysis suggesting a possible destabilizing effect for this substitution. This in turn is in agreement with the in vitro experimental data since it was found that the expression level of the mutant was lower than that of the wild type, while the activity level of the mutant was also decreased (Zhou et al., 2004). Missense3D, reports that no structural damage is detected for this substitution (i.e., the salt-bridge, about 3.6 Å between the charged groups, is not found) while VarSite underlines that the residue is highly conserved (in 110 aligned sequences), that the protein function could be perturbed and that there are not natural variants recorded in the gnomAD database at this position. The disease propensity score to the Asp to Asn change is in the scale of the method considered as low (value = 0.99). Taken together, for CYP1A2 Asp348Asn, the PolyPhen-2 prediction is not supported by the ΔΔG stability studies nor with the interactive analyses performed with the PyMol/ Chimera molecular graphics systems. The 3D analyses tend to support the experimental data and provide some hypotheses about the possible mechanisms resulting from the amino acid change. The automatic mapping Missense3D method does not flag this substitution while VarSite cautions about a possible impact on the structure/function of the modified protein.

### 3.1.2 | Ile386Phe (allele CYP1A2*4)

The PolyPhen-2 data suggest that this substitution is most likely damaging, in agreement with three ΔΔG structure-based approaches (except for MaestroWeb that suggests more a neutral effect) (Table 2). This residue is strictly conserved in our multiple sequence alignment (Figure S1), it is located in a loop and is mainly buried (SA = 13.52%) in a tightly packed hydrophobic and aromatic environment with several noncovalent interactions with its surrounding (node degree = 7) while being close to the Heme group (about 4 Å, Figure 1). It is predicted to be in a relatively rigid area. The residue is located in the

<footnote><sup>a</sup>For the program PopMusic solvent accessibility (SA) values are shown (in percent).</footnote>

<footnote><sup>b</sup>For the program MaestroWeb, the confidence estimation Cpred is shown (0.0-not reliable and 1.0-highly reliable).</footnote>

<footnote><sup>c</sup>MSA-Multiple sequence alignment.</footnote>

<footnote><sup>d</sup>Flexibility class was determined by the program PredyFlexy (rigid-0, intermediate-1, flexible-2).</footnote>

<footnote>Wild-type residue (bold) and amino acid substitution (Underlined).</footnote>
<table><caption>TABLE 3 In silico analysis of CYP2B6 (PDB ID: 3IBD, UniProt: P20813; GenBank: NM_000767.5)</caption><thead><tr><th>CYP2B6<br>Amino acid/<br>mutation/<br>alleles</th><th>DUET<br>ΔΔG, kcal/mol</th><th>PopMusic<sup>a</sup><br>ΔΔG, kcal/mol</th><th>SAAFEC<br>ΔΔG, kcal/mol</th><th>MAESTROWeb<sup>b</sup><br>ΔΔG, kcal/mol</th><th>PolyPhen-2<br>Score/<br>mutation<br>prediction</th><th>MSA<sup>c</sup><br>aa<br>conservation</th><th>Involved in<br>predicted<br>ligand-<br>binding<br>pockets/<br>pocket No.<br>(FTMap)</th><th>Involved in<br>predicted<br>predicted<br>or known<br>channels</th><th>Involved<br>in<br>predicted<br>PPI sites<br>(meta-<br>PPISP)</th><th>Node<br>degree<br>(RING-<br>2.0)</th><th>Predicted<br>fluctuation<br>value<br>(FlexPred)</th><th>Flexi-bility<br>class<sup>d</sup><br>(PredyFlexy)</th><th>Decreased<br>metabolic<br>activity<br>or protein<br>expression</th><th>References</th></tr></thead><tbody><tr><td>Gly99<br>G99E<br>CYP2B6*12</td><td>-2.245<br>Destabilizing</td><td>1.49<br>Destabilizing<br>SA = 3.3%</td><td>-0.491<br>Destabilizing</td><td>0.025<br>Destabilizing<br>Cpred = 0.924</td><td>1.000<br>Probably<br>damaging</td><td>High</td><td>No</td><td>Yes</td><td>No</td><td>3</td><td>1.052</td><td>1</td><td>Amino acid<br>changes<br>resulted<br>in almost<br>undetectable<br>enzyme<br>activity (in<br>vitro)</td><td>Lang et al.<br>(2004)</td></tr><tr><td>Arg140<br>R140Q<br>CYP2B6*14</td><td>-0.562<br>Destabilizing</td><td>0.32<br>Destabilizing<br>SA = 40.9%</td><td>-3.798<br>Destabilizing</td><td>0.471<br>Destabilizing<br>Cpred = 0.786</td><td>0.020<br>Benign</td><td>High</td><td>Yes/pocket<br>9</td><td>No</td><td>Yes</td><td>7</td><td>1.747</td><td>1</td><td>Reduced<br>expression<br>and/or<br>function<br>activity of<br>protein (in<br>vitro)</td><td>Lang et al.<br>(2004)</td></tr><tr><td>Met198<br>M198T<br>CYP2B6*27</td><td>-1.374<br>Destabilizing</td><td>1.12<br>Destabilizing<br>SA = 26.2%</td><td>-0.580<br>Destabilizing</td><td>-0.278<br>Stabilizing<br>Cpred = 0.832</td><td>0.000<br>Benign</td><td>High</td><td>Yes/pocket<br>1</td><td>No</td><td>Yes</td><td>10</td><td>1.223</td><td>0</td><td>Loss-of<br>function<br>allele *27<br>results in<br>85% decrease<br>in enzyme<br>activity</td><td>Roiger et al.<br>(2007)</td></tr><tr><td>Ile391<br>I39N<br>CYP2B6*15</td><td>-3.149<br>Destabilizing</td><td>2.71<br>Destabilizing<br>SA = 0.0%</td><td>-1.012<br>Destabilizing</td><td>0.992<br>Destabilizing<br>Cpred = 0.890</td><td>1.000<br>Probably<br>damaging</td><td>High</td><td>No</td><td>No</td><td>No</td><td>11</td><td>0.897</td><td>1</td><td>Amino acid<br>changes<br>resulted<br>in almost<br>undetectable<br>enzyme<br>activity (in<br>vitro)</td><td>Lang et al.<br>(2004)</td></tr></tbody></table>

<footnote><sup>a</sup>For the program PopMusic solvent accessibility (SA) values are shown (in percent).</footnote>

<footnote><sup>b</sup>For the program MaestroWeb the confidence estimation Cpred is shown (0.0-not reliable and 1.0-highly reliable).</footnote>

<footnote><sup>c</sup>MSA-Multiple sequence alignment.</footnote>

<footnote><sup>d</sup>Flexibility class was determined by the program PredyFlexy (rigid-0, intermediate-1, flexible-2).</footnote>

<footnote>Wild-type residue (bold).</footnote>
catalytic pocket and in the vicinity of the substrate access channel. Interactive analysis suggests that the Ile to Phe substitution should generate steric clashes, possibly locally destabilize the protein and impair the catalytic activity. Missense3D results do not highlight structural damage while VarSite underlines an aliphatic to aromatic substitution, that the residue is highly conserved (in 111 aligned sequences) and that the protein function could be perturbed as residue 386 has contact with a ligand in several CYP 3D structures. The disease propensity value (=1.13) is in this case labeled as high. For CYP1A2 Ile386Phe, the most likely damaging effect predicted by PolyPhen-2 is in good agreement with most 3D ΔΔG investigations and with our interactive 3D studies, further highlighted by the VarSite report card but not by the Missense3D results. Overall, most in silico investigations are in good agreement with the in vitro experimental observations, all supporting the basis for a decreased level of expression of the mutant protein (Chevalier et al., 2001; Zhou et al., 2004).

## 3.2 | CYP2B6

For CYP2B6 (PDB ID: 3IBD, UniProt: P20813), we investigated the possible molecular effect of four clinically important alterations (Table 3). The drugs metabolized by CYP2B6 include the prodrug cyclophosphamide, the antimalarial artemisinin, the anesthetics ketamine and propofol, and the HIV-1 reverse transcriptase inhibitors nevirapine and efavirenz (Lang et al., 2004). In vitro studies revealed that these CYP2B6 protein variants (Arg140Gln, Gly99Glu, and Ile391Asn) have a reduced expression level and/or activity. Two amino acid changes (Gly99Glu and Ile391Asn) resulted

in almost undetectable enzyme activity, despite the presence of residual protein level (Lang et al., 2004). Resequencing defined loss-of-function allele *27 (Met198Thr), resulted in 85% decrease in enzyme activity (Rotger et al., 2007).

The results of our computations are shown in Table 3. Four channels were found for the CYP2B6 protein with ChannelsDB (Pravda et al., 2018). From the four amino acids investigated here, only Gly99 seems to be involved in forming a channel. We predicted 19 ligand-binding pockets with the FTMap server (Brenke et al., 2009; Kozakov et al., 2015, 2011; Ngan et al., 2012) and Arg140 and Met198 belong to the predicted binding pocket numbers 9 and 1, respectively (Table S2).

### 3.2.1 | Gly99Glu (allele CYP2B6*12)

PolyPhen-2 suggests that this amino acid change is most likely damaging. This conclusion matches the results obtained with three structure-based programs that predict a destabilizing effect (MAESTROweb suggests more a neutral effect). This residue is buried (SA = 3.29%) and relatively close to the Heme group (Figure 1). It is fully conserved in our MSA (Figure S2) and part of the channel 2e (Figure 2) but it is not predicted to be part of a putative ligand-binding pocket nor is present in a putative PPI site. There is no room to accommodate a glutamate in this region of the protein and the amino acid change would in addition possibly put the substituted residue next to a negatively charged residue, Glu387. Most in silico results suggest that this substitution could be destabilizing in agreement with the in vitro experiments in which the amino acid change resulted in an almost undetectable enzyme activity (Lang et al., 2004). The automatic 3D

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** Overview of the CYP2B6 protein (PDB ID: 3IBD) with the predicted channels 2e (magenta, solid surface), 2b and 2f (yellow, solid surface). A solvent channel is shown in brown (solid surface)</figcaption></figure>
mapping tool, Missense3D, reports that the substitution re-
places a buried uncharged and small residue by a large and
charged one while VarSite underlines a large difference in
term of amino acid substitution and a likely change in the
protein's function (the residue is found highly conserved in
about 175 protein sequences). The disease propensity score is
high with a value of 1.67. For CYP2B6 Gly99Glu, essentially
all the silico tools tend to be in good agreement and assist the
analysis of the experimental studies.

**3.2.2** | **Arg140Gln (allele CYP2B6*14)**

PolyPhen-2 suggests that this mutation is benign while all the
structural approaches propose more a destabilizing effect. This
residue is not in the catalytic pocket and is in part solvent exposed
(SA = 41.28%) (Figure 1). It forms a salt-bridge with Glu148
(Figure 3). This residue is fully conserved in the MSA (Figure S2)
underlining its likely importance. Arg140 is located in a slightly
flexible region (Table 3; Figure S3). The residue is predicted to
be involved in a PPI interaction site (Figure S4) and indeed it
may be involved in electrostatic interaction with the P450 oxi-
doreductase protein (POR) (Lee et al., 2014). It is not predicted to
be involved in a channel (Table 3). The residue R140 has several
noncovalent interactions with its surrounding (node degree = 7)
(Figure 2). Thus, this substitution would likely damage the pro-
tein structure and function as predicted by most of the in silico
tools used with the exception of PolyPhen-2. The predictions of
the structure-based approaches are in agreement with the experi-
mental data regarding the reduced expression and/or activity of
the mutant protein (in vitro) (Lang et al., 2004). Missense3D, re-
ports no structural damage (the salt-bridge is not found but clearly
visible in Figure 3, the distance between the charged groups is
about 3.2 Å) while VarSite results underline that the substitu-
tion is somewhat conservative but suggest a likely change in the
protein's function (the residue is found highly conserved in about
190 protein sequences). Two natural variants at this position are
known (Trp and Pro). The disease propensity score is high with
a value of 1.19. Residue Arg140 is also reported by VarSite as
involved in PPI. For the CYP2B6 Arg140Gln, most 3D in silico
outputs tend to be in good agreement with the experimental stud-
ies and help to propose possible molecular mechanisms associ-
ated with the substitution. These 3D results are not in agreement
with the PolyPhen-2 prediction.

**3.2.3** | **Met198Thr (allele CYP2B6*27)**

PolyPhen-2 suggests that this alteration is benign while three
structure-based approaches propose a destabilizing effect
(MaestroWeb indicates a possible stabilizing effect) (Table 3).
Residue Met198 is relatively large and is essentially hydro-
phobic while Thr is small and polar. This residue is mainly

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** Part of CYP2B6 (PDB ID: 3IBD, cartoon representation) is presented with a zoom on residue Arg140 (blue). The protein is in gray and the heme (part of the catalytic site) in red-orange. Arg140 is partially solvent exposed, forming a salt-bridge (~3.2 Å) with Glu148 (red)</figcaption></figure>

buried (SA = 26.20%), located far from the catalytic site and
tightly packed in a very hydrophobic and aromatic environ-
ment (Figure 1). Although Met198 is not fully conserved in
our MSA, a hydrophobic residue is always present at this po-
sition. It belongs to a possible ligand-binding pocket as pre-
dicted by FTMap (pocket No. 1, Table S2) but it is not present
in the predicted channels. The residue makes several nonco-
valent interactions with its surrounding (node degree = 10;
Figure 4). Interactive 3D analysis indicates that the small Thr
residue is likely to create a destabilizing hole that could alter
proper packing of this region. Overall, the stability investiga-
tions and the structural analysis seem in agreement with the
experimental data that reports 87% decrease in protein activity
level for this substitution (Rotger et al., 2007). Missense3D
results suggest that no structural damage is expected due to
the substitution while VarSite underlines that the residue is
highly conserved in about 190 protein sequences. Yet, the dis-
ease propensity score is low with a value of 0.69. For CYP2B6
Met198Thr, most 3D stability investigations and the interac-
tive analysis tend to be in good agreement with the experimen-
tal studies. These investigations are not in agreement with the
PolyPhen-2 prediction.
### 3.2.4 | Ile391Asn (allele CYP2B6*15)

PolyPhen-2 suggests that this alteration is damaging. This is in agreement with the four structure-based computational packages that indicate a relatively strong destabilizing effect (Table 3). This residue is fully buried (SA = 0.00%) and located in a hydrophobic environment. The residue is essentially conserved or replaced with similar hydrophobic residues in our MSA (Figure S2). This residue is not predicted to be in PPI-binding surface nor it is part of a predicted ligand-binding pocket (Table 3). It is situated close to the channel 2f (Figure 2). It has several nonbonded interactions with its surrounding (node degree = 11). Taken together, the hypothesis here is that this amino acid replacement most likely destabilize the protein, in agreement with the in vitro experimental data that shows that this substitution resulted in an almost undetectable enzymatic activity (Lang et al., 2004). Missense3D reports that the substitution will bury a hydrophilic residue while VarSite notes that the residue is highly conserved in 181 protein sequences and that the substitution could impede the protein's function. The disease propensity score is high with a value of 1.99. For CYP2B6 Ile391Asn, most in silico predictions are in good agreement with the experimental studies.

### 3.3 | CYP2C9

CYP2C9 metabolizes approximately 15% of the clinically used drugs (Hirota, Eguchi, & Ieiri, 2013; Yiannakopoulou, 2013), including hypoglycemic agents, anticonvulsants,

anticoagulants (e.g., warfarin), nonsteroidal anti-inflammatory (e.g., diclofenac), antihypertensive (e.g., losartan), and diuretic drugs (Isvoran et al., 2017). We selected four clinically important alterations (Arg150His, Ile359Leu, Ile434Phe, and Arg132Gln) associated with altered drug metabolism by >50%. Table 4 shows the results of the in silico analysis for CYP2C9 (PDB ID: 1OG2, UniProt: P11712).

Six channels were found for the CYP2C9 protein (PDB ID: 1OG2) with ChannelsDB (Pravda et al., 2018). None of the studied amino acids is involved in these channels. With regard to the predicted ligand-binding pockets, 12 pockets were found for CYP2C9 using the FTMap approach (Brenke et al., 2009; Kozakov et al., 2015, 2011; Ngan et al., 2012). From the four investigated mutations, only one amino acid (Ile434) is a part of a predicted binding pocket, this pocket is indeed involved in the binding Heme (predicted pockets 1 and 5, Table S3).

### 3.3.1 | Arg132Gln (allele CYP2C9*33)

PolyPhen-2 predicts this alteration has most likely a damaging effect (Table 4). This is essentially in agreement with the four structure-based programs. Arg132 tends to be solvent exposed (SA = 54.88%). The residue is strictly conserved in our MSA (Figure S5) and predicted to be in a relatively rigid region, with several noncovalent interactions with its surrounding (node degree = 8). This residue is not predicted

<figure><img src="image_6.png" /><figcaption>**FIGURE 4** CYP2B6 residues interaction networks. This is a degree sorted (few interactions, small red dots, toward many interactions, large red dot, following the black arrow) circular layout representation of the computed RING network visualized in Cytoscape. The noncovalent interactions are shown as gray lines and each dot on the circle represents a residue (red). Here, we focus only on Met198 (yellow dot) to illustrate the use of such visualization. This residue (node) makes several noncovalent interactions with its surrounding (node degree = 10). The substitution by a smaller residue likely creates a destabilizing hole in the structure</figcaption></figure>
to be in a ligand-binding pocket nor in a channel but could be in a protein-protein interaction region (Table 4). The residue makes several noncovalent interactions with its surrounding (Table 4). Arg132 is located at the C-terminus region of an alpha helix (Figure 1). A hydrogen bond is detected between the Arg132 side chain (epsilon N) and the backbone carbonyl of Met129 (a residue located in the helix) that likely stabilizes this region of the protein. Furthermore, alpha helices have a large macroscopic dipole moment (–0.5 charge unit at the C-terminal). The arginine residue possibly contributes to the stabilization of the helix dipole, a role that the glutamine cannot play. Overall, it seems possible that this amino acid change locally destabilizes the protein and impedes its function. This would be in agreement with the decreased activity observed in vitro (Yin et al., 2008). In addition, Arg132 is predicted to be part of a PPI-binding region and indeed, it plays an important role in electrostatic interaction with the P450 oxidoreductase (POR) (Lee et al., 2014). The Missense3D report card indicates that the substitution should not damage the protein structure while on the VarSite report it is noted that a Arg to Gln is not dramatic change in term of properties, yet the residue is highly conserved in 181 protein sequences. The disease propensity score is labeled as high with a value of 1.19 and VarSite indicates also that residues equivalent to Arg132 are found at protein-protein interfaces in this family of proteins. For CYP2C9 Arg132Gln, most in silico predictions are in good agreement with the experimental studies.

### 3.3.2 | Arg150His (allele CYP2C9*8)

PolyPhen-2 labels this change as benign but the structural approaches (more strongly for DUET and SAAFEC) suggest that this amino acid substitution has destabilizing effects (Table 4). Arg150 is essentially solvent exposed (SA = 60.43%), and is involved in several noncovalent interactions with its surrounding (node degree = 8). It is located in the middle of an alpha helix (Figure 1). This residue is strictly conserved in our MSA (Figure S5) and it is not predicted to be in a very flexible region. It is involved in a salt-bridge network (with Asp143 and Glu147, both also interacting with Arg139, data not shown). Arg150 is not predicted to be located in a binding pocket, in a channel or in a PPI-binding region. Most likely the shorter histidine residue cannot form the above noted ionic and electrostatic interactions. The substitution could thus destabilize the structure locally. The results obtained with the 3D approaches are in agreement with the experimental data that shows a decreased activity for the variant protein in both in vitro and in vivo systems (Allabi, Gala, & Horsmans, 2005). Missense3D reports that the substitution will not damage the structure (the salt-bridge network is not seen) while the VarSite report card indicates that a Arg to His may not be dramatic change in term of physicochemical

properties, yet the residue is highly conserved in 181 protein sequences. Two natural variants have been identified (Leu and Cys) at position 150. The disease propensity score is high with a value of 1.45. For CYP2C9 Arg150His, most in silico stability prediction methods as well as the interactive structural analysis are in good agreement with the experimental studies. The results of PolyPhen-2 are not in agreement with these 3D methods nor they support the experimental data.

### 3.3.3 | Ile359Leu (allele CYP2C9*3)

PolyPhen-2 suggests that this alteration is benign but two structure-based approaches flag this substitution as possibly destabilizing while PopMusic and MaestroWeb propose a stabilizing effect (Table 4). Ile359 is buried (SA = 1.4%) and located on the C-terminus of a helical segment (Figure 1). This residue is strictly conserved in our MSA (Figure S5), it is not predicted to be in a highly flexible region. The residue has several noncovalent interactions with its surrounding (node degree = 9). It is not far from the Heme group, yet it is not predicted to be in a ligand-binding pocket, channels or in a PPI-binding region (Table 4). Interestingly, interactive structural analysis shows that the Ile359Leu substitution leads to steric clashes in the Heme region although the substitution is very conservative (data not shown). This would seem in agreement with experimental data which indicates a decreased enzymatic activity level for the variant protein (affecting the warfarin dose) (in vitro and in vivo) (King, Khan, Aithal, Kamali, & Daly, 2004; Shintani et al., 2001). Missense3D reports that the substitution will not damage the structure while VarSite mentions that the substitution should be tolerated, although it is also noted that the residue is highly conserved in 182 protein sequences. Two natural variants have been identified (Val and Thr) at position 359. The disease propensity score is low with a value of 0.39 but VarSite still highlights the proximity to the Heme region. For CYP2C9 Ile359Leu, stability prediction methods as well as the interactive structural analysis underline a possible impact on the protein catalytic activity in good agreement with the experimental studies. The results of PolyPhen-2 are not in agreement with the stability prediction methods, the interactive analysis nor with the experimental data.

### 3.3.4 | Ile434Phe (allele CYP2C9*59)

PolyPhen-2 suggests that this alteration is most likely damaging with a high score, in agreement with the four structure-based approaches (more strongly with DUET and PopMusic predictions; Table 4). This residue tends to be buried (SA = 22.05%) and is located near the Heme group. Ile434 makes some noncovalent interactions with its surrounding
(node degree = 6; Table 4). The residue is strictly conserved in our MSA (Figure S5) and predicted to be in a binding pocket and in a relatively rigid region. The substitution by a larger aromatic amino acid could damage the folding locally. It should also perturb the binding of the Heme group since it most likely changes the orientation of Cys435, a residue involved in the interaction with ring-shaped Heme molecule. Moreover, it is also predicted to be next to a PPI-binding region (not shown) and indeed Ile434 is close to the P450 oxidoreductase (POR)-binding area (Lee et al., 2014). All these in silico and structural analyses are in good agreement with the greatly decreased enzymatic activity of the variant protein observed in in vitro and in vivo models (Dai et al., 2015). Missense3D reports that the substitution will not damage the structure while the VarSite report card mentions that the substitution should be tolerated, but that the residue is fully conserved in 146 protein sequences. Two natural variants have been identified (Val and Thr) at position 359. The disease propensity score is high with a value of 1.13. VarSite highlights the proximity with the ligand-binding site and Heme and points to known protein-protein interactions in this protein family. For CYP2C9 Ile434Phe, basically all in silico predictions and the interactive analysis are all in good agreement and help to rationalize the experimental data.

## 3.4 | Complement Factor B

For Complement FB (PDB ID: 2OK5, UniProt ID: P00751), we investigated the possible molecular effect of four clinically important alterations involved in strengthening and facilitating the formation of an overactive C3 convertase (Table 5). The mutations were identified in patients with aHUS (Marinozzi et al., 2014; Roumenina et al., 2009). For this protein, we predicted putative ligand-binding pockets with the FTMap. Eleven pockets were identified (Table S4). None of the four amino acids investigated here were found to be in a predicted ligand-binding pockets.

### 3.4.1 | Arg113Trp

PolyPhen-2 suggests that this alteration is benign but three structure-based tools indicate that the substitution is destabilizing except for MaestroWeb (Table 5). This residue is highly conserved in the FB sequences from different species (Figure S6) and is not predicted to be in a highly flexible region. It is located at the C-terminus of a short beta strand (Figure 1), tends to be buried (SA = 36.5%; Table 5) and makes a salt-bridge with Asp134 (Figure 5). The large W side chain is expected to have local clashes with the surrounding, destroying a salt-bridge and possibly locally destabilizing the structure of the mutant protein. Arg113 does

make several noncovalent interactions with its surrounding (node degree 10, Figure 5), it does not seem to be located in a ligand-binding pocket but could be part of a protein-protein interaction site. It is difficult to predict if these structural changes will have significant impacts on the 3D structure, because the residue is located nearby the Cys133–Cys106 disulfide bond. The substitution could interfere with the formation of this covalent bond. Alternatively, this covalent bond may compensate the likely destabilizing effect of the substitution. Experimentally, the recombinant CFB was expressed normally without exhibiting a severe phenotype (Marinozzi et al., 2014). The lack of detectable functional consequences is likely related to the localization of the mutation in the Ba part of the molecule, which is cleaved and released in the process of the formation of the active C3 convertase C3bBb. This could also explain the benign phenotype, despite the fact that most predictive tools suggest a local destabilization of the protein domain. Missense3D reports that the substitution will not damage the structure while VarSite mentions that the substitution is important and that the residue is highly conserved in 45 protein sequences. VarSite also points to known protein-protein interactions in this region in this protein family. The disease propensity is high with a value of 1.93. For FB Arg113Trp, basically all 3D in silico predictions and the interactive analysis are all in good agreement with a likely effect on the protein structure. As this region of the protein is cleaved, the phenotype is benign. The 3D tools essentially give a different view as compared to PolyPhen-2. In this case the PolyPhen-2 prediction is in good agreement with the experimental data but this also seems to be a lucky prediction. This substitution also highlights the complexity of the biological system. Biophysical studies would be required to fully clarify the impact of the Arg113Trp substitution on the 3D structure.

### 3.4.2 | Arg178Gln

PolyPhen-2 suggests that this amino acid change is most likely damaging in agreement with all the structure-based programs (Table 5). This substitution is located in the middle of a short beta strand (Figure 1), is fully conserved in the MSA (Figure S6) and is almost fully buried (SA = 2.88%; Table 5). The residue is present on a complement control protein domain (CCP), in a region that is predicted to be rigid. The residue has numerous noncovalent interactions with its surrounding (node degree = 16, Figure 5) and forms a buried salt-bridge (predicted to be energetically strong) with Asp422. Arg178 does not seem to be in a ligand-binding pocket nor involved in a protein-protein interaction site (Table 5). The amino acid change to Gln should alter a strong salt-bridge and may also perturb the nearby Cys180–Cys140 disulfide bond. The in silico analysis and structural
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** Factor B. The salt-bridge involving Arg113 and Asp134 (~2.5 Å) is shown. Substitution with a Trp was performed with PyMol and some steric clashes are noticed. A nearby disulfide bond involving Cys106 and Cys133 is also highlighted</figcaption></figure>

analysis suggest a destabilizing effect, which is in agreement with some experimental data showing that Arg178Gln mutant was expressed partially cleaved and nonfunctional (Marinozzi et al., 2014). Missense3D reports that the substitution will replace a buried charge by an uncharged residue while VarSite mentions that the substitution may or may not be tolerated. The residue is very highly conserved in 45 protein sequences. VarSite gives a relatively high disease propensity score for this substitution with a value of 1.19. For FB Arg178Gln, basically all in silico predictions and the interactive analysis are all in good agreement with a likely effect of the substitution on the protein structure and with the experimental data.

### 3.4.3 | Asp254Gly

PolyPhen-2 suggests that this change is benign while three structure-based tools indicate that the substitution is modestly destabilizing except for SAAFEC that predicts a stabilization effect (Table 5). The residue is located in a loop on the von Willebrand domain of CFB (Figure 1), next to a Mg++ binding site, and in a relatively rigid segment. The residue is strictly conserved in the MSA (Figure S6) and partially solvent exposed (SA = 38.3%; Table 5). Its node degree was found to be 6 and this residue makes hydrogen bonds with Gln28 of the nearby complement control protein domain. The substitution could alter these interactions and could increase the flexibility of this loop. The residue is not predicted to be in a ligand-binding pocket but is expected to be part of a protein-protein binding site. Indeed, experimentally, the Asp254Gly is a gain-of-function mutation

leading to the formation of an overactive C3 convertase with enhanced binding to C3b (Hourcade, Mitchell, & Oglesby, 1999; Roumenina et al., 2009). This residue is indeed found to be in the binding region for C3b (Marinozzi et al., 2014; Rooijakkers et al., 2009; Roumenina et al., 2009). Interactive analysis suggests a stability change but the main role of this substitution is to affect a PPI-binding site. The functional consequence, without access to the structure of the C3 convertase (Rooijakkers et al., 2009) could not be predicted by the in silico tools used here, but some hypotheses could be made since the residue is predicted to be in a PPI-binding region. Interestingly, the type A domain of factor B is similar in structure to the type A domain of the complement receptor and integrin, CR3, in which the residue homologous to Asp254 is also a G. Therefore, the Asp254Gly substitution was studied long before its discovery in aHUS patients and served to delineate the C3b-binding site (Hourcade et al., 1999) and to create a more stable system, suitable for crystallization (Forneris et al., 2010; Pedersen et al., 2017). Missense3D report card indicates that the substitution should not damage the 3D structure while VarSite mentions that the substitution is significant and that the residue is very highly conserved in 46 protein sequences. VarSite gives a low disease propensity score for this substitution with a value of 1.00 and highlights that the substitution is part of a protein-protein interaction region in this family. For FB Asp254Gly, the in silico predictions and the interactive analysis have difficulties in the prediction of the possible impact of this change. Without experimental data, the key hypothesis that could be made is that the region is involved in protein-protein interactions.
### 3.4.4 | Phe261Leu

PolyPhen-2 suggests that this substitution is most likely damaging and three structure-based programs also indicate that the substitution is destabilizing except for SAAFEC that predicts a small stabilizing effect (Table 5). This residue is fully buried (SA = 0.34%) in a tightly packed hydrophobic (node degree = 14) and rigid environment (Table 5). The residue is present in the von Willebrand domain of CFB (Figure 1) and is fully conserved in the MSA (Figure S6). The amino acid change appears to induce steric clashes when investigated interactively and it could indeed be destabilizing. It is not predicted to be in a ligand-binding pocket nor seems to be directly involved in a protein-protein interaction site (Table 5). It is, however, located nearby a predicted PPI-binding site. Overall it would seem that this alteration changes the stability in this region of the protein and could affect protein-protein interactions. However, using only the in silico methods, it was not possible to directly predict an increased C3 binding. This mutant was found in a HUS patients and resulted in the formation of a distinct, rapidly cycling C3 convertase, characterized by a faster association but also a faster dissociation rate (Goicoechea de Jorge et al., 2007). Missense3D, reports that the substitution will not damage the 3D structure of the protein while VarSite mentions that the substitution is relatively conservative and that the residue is very highly conserved in 46 protein sequences. VarSite gives a low disease propensity score for this substitution with a value of 0.85. For FB Phe261Leu, the in silico predictions and the interactive analysis have difficulties in the prediction of the possible impact of the mutation. The signal provided by PolyPhen-2 and more explicitly by the approaches that predict a stability change is that the dynamic of this region of the protein could be changed. As the residue is buried but close to the surface, and close to a predicted PPI site, one could propose a possible impact on protein interactions, but here again, this is a very complex case for the in silico tools.

### 3.5 | Factor VIII

The C2 domain (or discoidin domain) of FVIII is known to interact with negatively charged membranes and with other proteins such as von Willebrand factor (VWF) (Liu et al., 2000; Pratt et al., 1999). We investigated the potential molecular effect of two substitutions, Met2238Val (moderate) and Ala2201Pro (mild), located on this domain of FVIII, both associated with hemophilia A (Liu et al., 2000; Pratt et al., 1999; Spiegel et al., 2001; Villoutreix & Miteva, 2016). Table 6 shows the results of the in silico analysis (PDB ID: 2R7E, UniProt ID: P00451). Twelve ligand-binding pockets

for the FVIII C2 domain were predicted using the FtMap server (Brenke et al., 2009; Kozakov et al., 2015, 2011; Ngan et al., 2012) (Table S5) and Met2238 and Ala2201 are not directly located in a predicted ligand-binding pockets. Yet, Met2238 seems next to a predicted PPI binding site while Ala2201 is predicted to be in a PPI-binding region.

### 3.5.1 | Ala2201Pro

PolyPhen-2 labels this change as most likely damaging (high score) in agreement with three structure-based approaches except DUET (Table 6). Ala2201 is well conserved in the sequence (Figure S7; Table 6) with limited noncovalent interactions with its surrounding (node degree = 4) in a region predicted to be relatively rigid. This residue is partially solvent exposed (36%; Table 6) and located on the N-terminal side of a short beta strand (Figure 1). It is predicted to be both, in a PPI-binding site (Table 6) and located in a membrane-binding site as analyzed by the PPM server prediction (Lomize et al., 2012; Figure 6). Interactive structural analysis suggests that its replacement by a Pro should induce a structural change in the nearby loop region of the domain while some steric clashes are also noticed upon introduction of the mutation in the 3D structure. The in silico and interactive structural analysis thus propose a moderate stability change substitution. Experimentally, Spiegel, Murphy, and Stoddard (2004) have shown in vitro that while the substitution does not destabilize significantly the domain, it nevertheless impedes membrane binding, although the interaction with VWF is not affected and the amino acid replacement. It is interesting to note that the in silico analysis predicts this residue to be near a PPI-binding site. Indeed, it is known that some patients with hemophilia A develop inhibitory antibodies against FVIII. The FVIII C2 domain has been cocrystallized with an antibody prepared from a cell line derived from the memory B-cell repertoire of a patient with hemophilia A (Spiegel et al., 2001). This Fab fragment interacts with the membrane-binding loops of the FVIII C2 domain where A2201 is located, thereby impeding the interaction of the cofactor with the membrane. This event most likely occurs when FVIII dissociates from VWF. It is also interesting to note that upon extraction of the Fab fragment from the experimental structure reported by Spiegel et al. (2001) (PDB ID: 1IQD), it was possible to dock the antibody in a manner that is almost identical to the X-ray structure (best docking score) with the pyDockWeb server (Jiménez-García et al., 2013) (data not shown). Further in vivo experimental work has indicated that infusion of FVIII to patients carrying the Ala2201Pro substitution induces a strong immunological reaction and the development of inhibitors against FVIII (Ettinger, James, Kwok,
<table><caption>TABLE 6 In silico analysis of blood Factor VIII (FVIII) (PDB ID: 2R7E, UniProt ID: P00451; GenBank: NM_000132.3)</caption><thead><tr><th>FVIII</th><th>DUET</th><th>PopMusic<sup>a</sup></th><th>SAAFEC</th><th>MAESTROweb<sup>b</sup></th><th>PolyPhen-2</th><th>MSA<sup>c</sup></th><th>Involved in predicted ligand-binding</th><th>Involved in predicted protein-protein interaction</th><th>Node degree (RING-2.0)</th><th>Predicted fluctuation value (FlexPred)</th><th>Flexibility class<sup>e</sup></th><th>Experimental data</th><th>References</th></tr><tr><th>Amino acid/ mutation</th><th>ΔΔG, kcal/mol</th><th>ΔΔG, kcal/mol</th><th>ΔΔG, kcal/mol</th><th>ΔΔG, kcal/mol</th><th>Score/mutation prediction</th><th>aa</th><th>pockets/ pocket No. (FTMap)</th><th>Surface exposure (%)<sup>d</sup></th><th>sites (meta-PPISP)</th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>Ala2201</td><td>0.378</td><td>0.23</td><td>-0.985608</td><td>0.196</td><td>0.988</td><td>High</td><td>No</td><td>36</td><td>Yes</td><td>4</td><td>0.075</td><td>1</td><td>Mutation associated with mild hemophilia A</td></tr><tr><td>A2201P</td><td>Stabilizing</td><td>Stabilizing</td><td>Stabilizing</td><td>Cpred = 0.887</td><td>Probably damaging</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>The mutant protein domain in vitro seems to have molecular functions similar to the wild type (investigated properties: stability, and binding to the vWF) but membrane binding is damaged</td></tr><tr><td>Met2238</td><td>-0.503</td><td>0.27</td><td>0.553310</td><td>0.248</td><td>0.152</td><td>High</td><td>No</td><td>11</td><td>No but very close to a proposed site</td><td>8</td><td>1.487</td><td>1</td><td>Mutation associated with moderate hemophilia A</td></tr><tr><td>M2238V</td><td>Stabilizing</td><td>Stabilizing</td><td>Stabilizing</td><td>Cpred = 0.874</td><td>Benign</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>The mutant protein domain in vitro seems to have molecular functions similar to the wild type (investigated properties: stability, and binding to the vWF) (in vitro)</td></tr></tbody></table>

<footnote><sup>a</sup>For the program PopMusic solvent accessibility (SA) values are shown (in percent).</footnote>

<footnote><sup>b</sup>For the program MaestroWeb the confidence estimation Cpred is shown (0.0-not reliable and 1.0-highly reliable).</footnote>

<footnote><sup>c</sup>MSA-Multiple sequence alignment.</footnote>

<footnote><sup>d</sup>Liu et al. (2000).</footnote>

<footnote><sup>e</sup>Flexibility class was determined by the program PredyFlexy (rigid-0, intermediate-1, flexible-2).</footnote>

<footnote>Wild-type residue (bold) and amino acid substitution (Underlined).</footnote>
<figure><img src="image_8.png" /><figcaption>**FIGURE 6** Factor VIII point mutations. The predicted position of the Factor VIII C2 domain on the membrane (in blue), as well as the point mutations investigated in this study are shown</figcaption></figure>

Thompson, & Pratt, 2010) highlighting the complex anti-
genic changes associated with the 3D structure of the loop
harboring amino acid change. Missense3D reports that the
substitution will not damage the 3D structure of the pro-
tein while VarSite mentions that the residue is very highly
conserved in 28 protein sequences. VarSite gives a high
disease propensity score for this substitution with a value
of 1.58. For FVIII Ala2201Pro, the in silico predictions and
the interactive analysis have difficulties in explaining the
possible impact of the amino acid substitution on the struc-
ture. Some insights come from some 3D prediction tools
highlighting a possible PPI-binding site and protein–mem-
brane interaction site. Clearly, the in silico investigations
can help in generating some hypotheses and in designing
experiments but cannot fully explain such complex mo-
lecular events.

**3.5.2** | Met2238Val

PolyPhen-2 suggests that this change is benign and three
structure-based approaches propose a small destabilizing ef-
fect with the exception of SAAFEC (Table 6). This residue
is essentially buried (SA = 10.0%) with some noncovalent
hydrophobic interactions with its surrounding (node de-
gree = 8). M2238 is located in a loop structure, next to the

interface with the A domains of FVIII and it is in a relatively
rigid area (Figure 1). This residue is well conserved in the
sequences or may be replaced, indeed with a valine (Figure
S7; Table 6). Interactive structural analysis indicates that
the valine side chain can be accommodated in this region of
the domain suggesting that if this substitution is involved in
bleeding, it may act via unclear mechanism such as the splic-
ing of FVIII or other modifications, as suggested by Spiegel
et al. (Spiegel et al., 2004). Missense3D reports that the sub-
stitution will not damage the 3D structure of the protein while
VarSite mentions that the residue is very highly conserved
in 28 protein sequences. VarSite gives a low disease pro-
pensity score for this substitution with a value of 0.58. For
FVIII Met2238Val, the in silico predictions and the interac-
tive analysis have difficulties in explaining the impact of this
amino acid substitution. Taken together, this analysis sug-
gests that more experiments are required to understand if the
mutation is responsible for the bleeding phenotype.

**3.6** | Antithrombin

For Antithrombin (AT) (PDB ID: 2BEH, UniProt ID: P01008) we investigated the possible molecular effects of four clinically important amino acid substitutions. AT is a major blood coagulation protease inhibitor and becomes an effective inhibitor upon binding to heparin, a negatively charged and highly sulfated polysaccharide that enhances the anticoagulant activity of AT toward its main targets, thrombin, factor IXa, and factor Xa (Huntington, 2011). The crystal structure of antithrombin–heparin and thrombin has been reported several years ago (Li, Johnson, Esmon, & Huntington, 2004). Mutations in the gene encoding AT can cause AT deficiency that predisposes affected individuals to venous and arterial thromboembolism (Ding et al., 2013). According to the phenotype, AT deficiency has been classified into two main categories, type I, characterized by a parallel reduction of the AT antigen and activity levels and type II, defined by a reduction of the AT activity and normal or almost normal antigen levels. The type II AT deficiency has been further subdivided into type IIb with a heparin-binding defect (type IIHBS) and type IIa with a reactive center loop (RCL) defect and type IIc with pleiotropic defects (Patnaik & Moll, 2008). We initially investigated all the amino acid changes reported in (Luxembourg, D'Souza, Körber, & Seifried, 2015) and then selected two variants with type I deficiency and two variants with type II deficiency. Table 7 shows the results of the prediction for these AT mutants. Using FTMap web server (Brenke et al., 2009; Kozakov et al., 2015, 2011; Ngan et al., 2012) we identified 12 essentially unknown putative ligand-binding cavities for AT. None of the selected substituted amino acids seem to be located in these predicted pockets (Table S6). Of importance, the binding site for heparin
<table><caption>TABLE 7 In silico analysis of blood Antithrombin (AT) (PDB ID: 2BEH, UniProt ID: P01008; GenBank: NM_000488.3)</caption><thead><tr><th>AT</th><th>Amino acid/ mutation</th><th>DUET</th><th>PopMusic<sup>a</sup></th><th>SAAFEC</th><th>MAESTROweb<sup>b</sup></th><th>PolyPhen-2 Score/ mutation prediction</th><th>MSA<sup>c</sup> Conservativity level</th><th>Involved in predicted ligand binding pockets/ pocket No. (FTMap)</th><th>Possible interaction with heparin<sup>d</sup> (ClusPro)</th><th>Involved in predicted PPI sites (meta-PPISP)</th><th>Node degree (RING-2.0)</th><th>Predicted fluctuation value (FlexPreds)</th><th>Flexibility class<sup>e</sup> (PredyFlexy)</th><th>Experimental data</th><th>References</th></tr></thead><tbody><tr><td colspan="16" style="text-align: left; font-weight: bold;">AT type I</td></tr><tr><td><b>Ser82</b></td><td>-0.439</td><td>1.64</td><td>-0.509178</td><td>0.025</td><td>1.000</td><td>High</td><td>No</td><td>No</td><td>No</td><td>12</td><td>1.478</td><td>2</td><td>AT type I deficiency</td><td>Ding et al. (2013)</td></tr><tr><td><b>S82R</b></td><td>Destabilizing</td><td>Destabilizing</td><td>Destabilizing</td><td>Destabilizing</td><td>Cpred = 0.885</td><td>Probably damaging</td><td>No</td><td>No</td><td>Yes</td><td>9</td><td>1.016</td><td>2</td><td>AT type I deficiency</td><td>Ding et al. (2013)</td></tr><tr><td><b>Cys95</b></td><td>-0.384</td><td>1.67</td><td>-2.904263</td><td>1.550</td><td>1.000</td><td>High</td><td>No</td><td>No</td><td>Yes</td><td>5</td><td>2.429</td><td>2</td><td>AT type I deficiency</td><td>Olson et al. (2010). Luxembourg et al. (2015)</td></tr><tr><td><b>C95R</b></td><td>Destabilizing</td><td>Destabilizing</td><td>Destabilizing</td><td>Destabilizing</td><td>Cpred = 0.785</td><td>Probably damaging</td><td>No</td><td>No</td><td>Yes</td><td>9</td><td>1.243</td><td>2</td><td>AT type IIHBS deficiency</td><td>Olson et al. (2010). Luxembourg et al. (2015)</td></tr><tr><td colspan="16" style="text-align: left; font-weight: bold;">AT type II HBS</td></tr><tr><td><b>Arg13</b></td><td>-0.358</td><td>0.22</td><td>0.678068</td><td>0.124</td><td>1.000</td><td>High</td><td>No</td><td>No</td><td>Yes</td><td>5</td><td>2.429</td><td>2</td><td>AT type I deficiency</td><td>Olson et al. (2010). Luxembourg et al. (2015)</td></tr><tr><td><b>R13W</b></td><td>Destabilizing</td><td>Destabilizing</td><td>Increase</td><td>Destabilizing</td><td>Cpred = 0.752</td><td>Probably damaging</td><td>No</td><td>No</td><td>Yes</td><td>9</td><td>1.243</td><td>2</td><td>AT type IIHBS deficiency</td><td>Olson et al. (2010). Luxembourg et al. (2015)</td></tr><tr><td><b>Ser116</b></td><td>-0.022</td><td>0.15</td><td>-0.588801</td><td>-0.117</td><td>0.290</td><td>High</td><td>No</td><td>No</td><td>No</td><td>9</td><td>1.243</td><td>2</td><td>AT type IIHBS deficiency</td><td>Olson et al. (2010). Luxembourg et al. (2015)</td></tr><tr><td><b>S116P</b></td><td>Destabilizing</td><td>Destabilizing</td><td>Destabilizing</td><td>Stabilizing</td><td>Cpred = 0.923</td><td>Benign</td><td>No</td><td>No</td><td>No</td><td>9</td><td>1.243</td><td>2</td><td>AT type IIHBS deficiency</td><td>Olson et al. (2010). Luxembourg et al. (2015)</td></tr></tbody></table>

<footnote><sup>a</sup>For the program PopMusic solvent accessibility (SA) values are shown (in percent).</footnote>

<footnote><sup>b</sup>For the program MaestroWeb the confidence estimation Cpred is shown (0.0-not reliable and 1.0-highly reliable).</footnote>

<footnote><sup>c</sup>MSA-Multiple sequence alignment.</footnote>

<footnote><sup>d</sup>Heparin-AT docking were carried out with the ClusPro server (see the Methods section).</footnote>

<footnote><sup>e</sup>Flexibility class was determined by the program PredyFlexy (rigid-0, intermediate-1, flexible-2).</footnote>

<footnote>Wild-type residue (bold) and amino acid substitution (Underlined).</footnote>
differs from what is known regarding pockets involved in
binding small drug-like molecules. As such, we do not expect
to really identify the heparin-binding site with binding pocket
predictors. This is the reason why we used a docking engine
to predict the binding of heparin to AT.

**3.6.1 | Ser82Arg**

PolyPhen-2 suggests that this substitution is most likely dam-
aging with a high score, in agreement with the four structure-
based approaches where the computed destabilization is
relatively strong (Table 7). Ser82 is fully buried (SA = 0.00%)
and is located in a hydrophobic and aromatic environment
on a helical structure (Figure 1). This residue is fully con-
served in our MSA (Figure S8), is not included in a predicted
ligand-binding pockets (Table S6), and cannot be a part of a
PPI-binding site (Table 7). It has several nonbonded interac-
tions with its surrounding (node degree = 12) and packed in
a region where some flexibility is required (i.e., located in a
region where the RCL inserts when it is cleaved upon interac-
tion with the target proteases). The Ser82Arg substitution gen-
erates numerous steric clashes, suggesting it can cause folding
problems. The structural analysis is thus in agreement with
the experimentally observed type I deficiency. Missense3D
reports that the substitution will bury a charged residue while
VarSite mentions that the residue is very highly conserved in
108 protein sequences and could impact the function of the
protein. Yet VarSite gives a low disease propensity score for
this substitution with a value of 0.73. Taken together, all the in
silico prediction tools suggest that the AT Ser82Arg substitu-
tion could impact the protein structure, in agreement with the
experimental data (Ding et al., 2013).

**3.6.2 | Cys95Arg**

PolyPhen-2 labels this change as most likely damaging with
a high score, in agreement with the four structure-based ap-
proaches where the computed destabilization is strong (Table
7). Cys95 is buried (SA = 3.79%) and is located on a moder-
ately flexible loop structure (Figure 1). This residue is fully
conserved in the MSA (Figure S8) underlining its impor-
tance. Cys95 residue is not predicted to be part of putative
ligand-binding pockets (Table S6) but it could be part of a
PPI-binding site. The residue has several nonbonded interac-
tions with its surrounding (node degree = 10) and forms a
disulfide bond with Cys21. This change is predicted to be
destabilizing, which is in agreement with the type I pheno-
type. Missense3D reports that the substitution will disrupt
a disulfide bond while VarSite mentions that the residue is
very highly conserved in 109 protein sequences and that such
substitution is highly unfavorable and will destroy a disulfide

bond. VarSite gives a very high disease propensity score for
this substitution with a value of 3.27. All in silico predictions
suggests that the AT Cys95Arg substitution could impact the
protein, in agreement with the experimental data.

**3.6.3 | Arg13Trp**

PolyPhen-2 labels this change as most likely damaging with a
high score, in agreement with three structure-based approaches
while SAAFEC indicates a possible increased stability (Table
7). The predicted destabilization is, however, relatively mild.
This residue is fully conserved in the AT sequences (Figure
S9), it is not involved in a PPI site and not predicted to be
in a ligand-binding pocket (Table S6). Indeed, although it is
known experimentally that this zone binds heparin (Johnson
et al., 2006), the region does not resemble a typical binding
pocket cavity as this surface is essentially flat. Arg13 residue
is mainly solvent exposed (SA = 74.69%) and is predicted to
be in a relatively flexible region (Figure 1). Docking of hepa-
rin to AT is well predicted via the ClusPro server (Kozakov
et al., 2017; Mottarella et al., 2014) and the best pose found
in silico is very similar to the experimental structure (Li et al.,
2004). The docked heparin indicates that Arg13 (as well as
Arg47 and Arg129 residues) should be part of the binding site
(Figure 7) in agreement with what is known experimentally.
Thus, even if the crystal structure of the AT-heparin complex
was not available, the in silico approach could have reliably
provided insight into the possible molecular impact of the
substitution on the structure, function, and heparin binding
of the variant protein. The replacement of the Arg by a Trp
could indeed alter the binding of heparin to this region of AT
without affecting its heparin-independent inhibitory function,
in agreement with a type II HBS phenotype for the mutation.
Missense3D proposes that the substitution does not damage
the 3D structure while VarSite mentions that the residue is
very highly conserved in 92 protein sequences and that such
substitution is highly unfavorable. The residue is also found to
interact with ligands (i.e., heparin-like molecules) in several
related protein structures. VarSite gives a high disease pro-
pensity score for this substitution with a value of 1.93. Taken
together, most in silico prediction tools directly or indirectly
suggests that the AT Arg13Trp substitution could impact the
function protein, in agreement with the experimental data.

**3.6.4 | Ser116Pro**

PolyPhen-2 suggests that this type IIHBS AT change is be-
nign while three structure-based approaches propose a mod-
erately destabilizing effect while MaestroWeb indicates a
possible stabilizing effect (Table 7). Ser116 residue is fully
conserved in the AT sequences (Figure S9), it is not involved
<figure><img src="image_9.png" /><figcaption>**FIGURE 7** Docking of heparin to AT (PDB ID: 2BEH) using the ClusPro server. The reactive center loop is in red and the docked heparin is in green</figcaption></figure>

in a predicted PPI site and not part of a predicted ligand-binding pockets (Table S6). The residue is predicted to be in a relatively flexible region and does make several noncovalent interactions with its surrounding (node degree = 9). Ser116 is buried (SA = 3.64%) and located in the N-terminal side of an alpha helix involved in heparin binding (Figures 1 and 7). Destabilizing effect and local misfolding is possible since the amino acid proline can lead to a change in the orientation of the protein backbone chain. Here again the docking of heparin on the AT surface would have shed light on the potential impact of the amino acid substitution in the absence of an experimental structure of the complex. Taken together, the in silico analyses in 3D provide valuable information on the effect of this amino acid change. This is in agreement with the observed type IIHBS phenotype although the use of PolyPhen-2 alone here is not sufficient to grasp the possible role of this residue. Missense3D indicates that the substitution will break some hydrogen bonds while the VarSite report card mentions that the residue is very highly conserved in 103 protein sequences and that this residue is also found

to interact with a ligand (i.e., heparin-like molecules) in one experimental protein structure. VarSite gives a high disease propensity score for this substitution with a value of 1.32. Taken together, the 3D in silico prediction tools suggest that the AT Ser116Pro substitution could impact the protein and the function, in agreement with the experimental data. Here PolyPhen-2 suggested a benign effect.

4 | CONCLUSION

Numerous in silico prediction packages have been developed to help prioritize potentially deleterious variants. In addition, other in silico tools that have not been developed for the investigation of variants can also be used in an attempt to provide additional insights into the mechanism by which an amino acid substitution may affect the protein structure and/or function.

The present study investigates the possible impact of different types of amino acid changes on the protein structure and/or function using different computational approaches. Some substitutions are easier to analyze while others involve highly complex mechanisms. Among the 20 substitutions investigated here, the experimental data/clinical data are unclear for the FVIII Met2238Val substitution while the possible structural damage is not known for the FB Arg113Trp change. For the other amino acid substitutions, we observed that the PolyPhen-2 predictions tend to agree with experimental data about 60% of the time and in most cases it is not possible to make suggestion about the possible underlying molecular mechanisms. In silico approaches evaluating stability change (stabilizing or destabilizing), tend to agree most of the time with experimental observations, even though, the molecular mechanisms involved cannot always be explained (e.g., this is for example the case when a substitution favors macromolecular interactions). Conflicting stability predictions among the different tools are known, this can be due to many reasons including scoring function parameters or definition of a threshold to label the substitution as stabilizing or destabilizing. Indeed, it has been shown that for such computations, only the absolute value of the change matters, not the sign (Peng & Alexov, 2016; Petukh, Kucukkal, & Alexov, 2015). Interactive structural analysis can here be used to investigate further the predicted ΔΔG values. When this step is carried out, it is most often possible to clarify the calculation output. When PolyPhen-2 and 3D stability predictions differ, the residue that is substituted is most often involved in a salt-bridge in the wild-type protein structure. The stabilizing or destabilizing nature of a substitution is difficult to estimate with fast computational approaches, but still a stability change can point toward some modification in the dynamics of the system that, if associated with interactive 3D analysis, can provide some insights about the possible impact of the amino acid change. Of importance, when a substitution introduces
a titratable residue in the core interior of a protein, the residue may not be charged and pKa computations can provide some insights. Yet, in such circumstance, there is usually an energetic price to pay. In general, in most situations, interactive structural analysis associated with prediction of stability changes and prediction of binding pockets, channels, protein-protein interaction sites or protein-heparin interactions (i.e., prediction of hotspots or prediction making use of docking) provide the users with some rational hypotheses about molecular mechanisms that 2D approaches cannot offer at this time. Definitively, structure-based programs combined with interactive structural analysis can provide detailed information about the effect of the amino acid change. In this case, the generated hypotheses can in general be rationally translated into assays. 3D approaches complement and even challenge sequence-based approaches. The combined use of 2D and 3D computational tools is therefore highly recommended, even though, in most clinical genomic centers, the 3D approaches are for the time being largely ignored. Online 3D mapping tools such as VarSite and Missense3D should definitively help biologists with no training in structural biology or structural bioinformatics to rationalize data in 3D. Yet we observe in our dataset that for several substitutions and essentially the ones involving salt-bridges or salt-bridge networks, the provided report cards do not report the ionic interactions. For these types of charged substitution, interactive investigations are definitively needed.

## ACKNOWLEDGMENTS

Supports from the Inserm Institute, University of Lille, Lille I-Site and Lille Region are greatly appreciated.

## CONFLICT OF INTEREST

The authors have no conflict of interest to declare.

## DATA AVAILABILITY STATEMENT

Data sharing is not applicable to this article as no new missense mutation data were created in this study.

## ORCID

Bruno O. Villoutreix D https://orcid.org/0000-0002-6456-7730

## REFERENCES

Abildgaard, U. (2007). Antithrombin—Early prophecies and present challenges. *Thrombosis and Haemostasis*, *98*, 97–104. https://doi.org/10.1160/TH07-04-0277

Adzhubei, I. A., Schmidt, S., Peshkin, L., Ramensky, V. E., Gerasimova, A., Bork, P., ... Sunyaev, S. R. (2010). A method and server for predicting damaging missense mutations. *Nature Methods*, *7*, 248–249. https://doi.org/10.1038/nmeth0410-248

Allabi, A. C., Gala, J. L., & Horsmans, Y. (2005). CYP2C9, CYP2C19, ABCB1 (MDR1) genetic polymorphisms and phenytoin metabolism in a Black Beninese population. *Pharmacogenetics and*

Genomics, 15(11), 779–786. https://doi.org/10.1097/01.fpc.0000174787.92861.91

Bava, K. A., Gromiha, M. M., Uedaira, H., Kitajima, K., & Sarai, A. (2004). ProTherm, version 4.0: Thermodynamic database for proteins and mutants. *Nucleic Acids Research*, *32*, D120–D121. https://doi.org/10.1093/nar/gkh082

Bayton, T., & Lane, D. (2003). *Antithrombin mutation data base*. *Departement de Haematology*. London, UK: Imperial College Faculty of Medicine, Charing Cross Hospital Campus. Available at http://www1.imperial.ac.uk/departmentofmedicine/divisions/

Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., ... Bourne, P. E. (2000). The protein data bank. *Nucleic Acids Research*, *28*, 235–242. https://doi.org/10.1093/nar/28.1.235

Bhattacharya, A. A., Lupher, M. L. Jr, Staunton, D. E., & Liddington, R. C. (2004). Crystal structure of the A domain from complement factor B reveals an integrin-like open conformation. *Structure*, *12*(3), 371–378. https://doi.org/10.1016/j.str.2004.02.012

Brenke, R., Hall, D. R., Chuang, G.-Y., Comeau, S. R., Bohnuud, T., Beglov, D., ... Kozakov, D. (2012). Application of asymmetric statistical potentials to antibody-protein docking. *Bioinformatics*, *28*(20), 2608–2614. https://doi.org/10.1093/bioinformatics/bts493

Brenke, R., Kozakov, D., Chuang, G. Y., Beglov, D., Hall, D., Landon, M. R., ... Vajda, S. (2009). Fragment-based identification of drugable 'hot spots' of proteins using Fourier domain correlation techniques. *Bioinformatics*, *25*(5), 621–627. https://doi.org/10.1093/bioinformatics/btp036

Cardenas, J. C., Rein-Smith, C. M., & Church, F. C. (2016). Overview of blood coagulation and the pathophysiology of blood coagulation disorders. *Encyclopedia of Cell Biology*, *1*, 714–722.

Chen, H.-L., & Zhou, H.-X. (2005). Prediction of interface residues in protein-protein complexes by a consensus neural network method: Test against NMR data. *Proteins*, *61*, 21–35. https://doi.org/10.1002/prot.20514

Chevalier, D., Cauffiez, C., Allorge, D., Lo-Guidice, J. M., Lhermitte, M., Lafitte, J. J., & Broly, F. (2001). Five novel natural allelic variants-951A>C, 1042G>A (D348N), 1156A>T (I386F), 1217G>A (C406Y) and 1291C>T (C431Y)-of the human CYP1A2 gene in a French Caucasian population. *Human Mutation*, *17*(4), 355–356. https://doi.org/10.1002/humu.49

Cooper, P. C., Coath, F., Daly, M. E., & Makris, M. (2011). The phenotypic and genetic assessment of antithrombin deficiency. *International Journal of Laboratory Hematology*, *33*, 227–237. https://doi.org/10.1111/j.1751-553X.2011.01307.x

Dai, D. P., Wang, S. H., Li, C. B., Geng, P. W., Cai, J., Wang, H., ... Cai, J. P. (2015). Identification and functional assessment of a new CYP2C9 allelic variant CYP2C9*59. *Drug Metabolism and Disposition*, *43*(8), 1246–1249.

Daly, A. K. (2015). Polymorphic variants of Cytochrome P450: Relevance to cancer and other diseases. *Advances in Pharmacology*, *74*, 85–111. https://doi.org/10.1016/bs.apha.2015.03.001

de Brevern, A. G., Bornot, A., Craveur, P., Etchebest, C., & Gelly, J.-C. (2012). PredyFlexy: Flexibility and local structure prediction from sequence. *Nucleic Acids Research*, *40*, W317–W322. https://doi.org/10.1093/nar/gks482

Dehouck, Y., Kwasigroch, J. M., Gilis, D., & Rooman, M. (2011). PoPMuSiC 2.1: A web server for the estimation of protein stability changes upon mutation and sequence optimality. *BMCBioinformatics*, *12*, 151–162. https://doi.org/10.1186/1471-2105-12-151

Ding, Q., Wang, M., Xu, G., Ye, X., Xi, X., Yu, T., ... Wang, H. (2013). Molecular basis and thrombotic manifestations of Antithrombin
deficiency in 15 unrelated Chinese patients. *Thrombosis Research*, **132**, 367-373. https://doi.org/10.1016/j.thromres.2013.07.013

Elfaki, I., Mir, M., Almutairi, F. M., & Duhier, F. M. A. (2018). Cytochrome P450: Polymorphisms and roles in cancer, diabetes and atherosclerosis. *Asian Pacific Journal of Cancer Prevention*, **19**(8), 2057-2070. https://doi.org/10.22034/APJCP.2018.19.8.2057

Ettinger, R. A., James, E. A., Kwok, W. W., Thompson, A. R., & Pratt, K. P. (2010). HLA-DR-restricted T-cell responses to factor VIII epitopes in a mild haemophilia A family with missense substitution A2201P. *Haemophilia*, **16**(102), 44-55. https://doi.org/10.1111/j.1365-2516.2008.01905.x

Forneris, F., Ricklin, D., Wu, J., Tzekou, A., Wallace, R. S., Lambris, J. D., & Gros, P. (2010). Structures of C3b in complex with factors B and D give insight into complement convertase formation. *Science*, **330**(6012), 1816-1820. https://doi.org/10.1126/science.1195821

Fujikura, K., Ingelman-Sundberg, M., & Lauschke, V. M. (2015). Genetic variation in the human cytochrome P450 supergene family. *Pharmacogenetics and Genomics*, **25**, 584-594. https://doi.org/10.1097/FPC.0000000000000172

Gaedigk, A., Ingelman-Sundberg, M., Miller, N. A., Leeder, J. S., Whirl-Carrillo, M., & Klein, T. E.; PharmVar Steering Committee. (2018). The Pharmacogene Variation (PharmVar) Consortium: Incorporation of the Human Cytochrome P450 (CYP) Allele Nomenclature Database. *Clinical Pharmacology & Therapeutics*, **103**(3), 399-401. https://doi.org/10.1002/cpt.910

Getov, I., Petukh, M., & Alexov, E. (2016). SAAFEC: Predicting the effect of single point mutations on protein folding free energy using a knowledge-modified MM/PBSA approach. *International Journal of Molecular Sciences*, **17**, 512-525. https://doi.org/10.3390/ijms17040512

Goicoechea de Jorge, E., Harris, C. L., Esparza-Gordillo, J., Carreras, L., Arranz, E. A., Garrido, C. A., ... Rodríguez de Córdoba, S. (2007). Gain-of-function mutations in complement factor B are associated with atypical hemolytic uremic syndrome. *Proceedings of the National Academy of Sciences of the United States of America*, **104**, 240-245. https://doi.org/10.1073/pnas.0603420103

Gonnelli, G., Rooman, M., & Dehouck, Y. (2012). Structure-based mutant stability predictions on proteins of unknown structure. *Journal of Biotechnology*, **161**(3), 287-293. https://doi.org/10.1016/j.jbiotec.2012.06.020

Grimm, D. G., Azencott, C.-A., Aicheler, F., Gieraths, U., MacArthur, D. G., Samocha, K. E., ... Borgwardt, K. M. (2015). The evaluation of tools used to predict the impact of missense variants is hindered by two types of circularity. *Human Mutation*, **36**(5), 513-523. https://doi.org/10.1002/humu.22768

Hirota, T., Eguchi, S., & Ieiri, I. (2013). Impact of genetic polymorphisms in CYP2C9 and CYP2C19 on the pharmacokinetics of clinically used drugs. *Drug Metabolism and Pharmacokinetics*, **28**, 28-37. https://doi.org/10.2133/dmpk.DMPK-12-RV-085

Homeyer, N., & Gohlke, H. (2012). Free energy calculations by the molecular mechanics Poisson-Boltzmann surface area method. *Molecular Informatics*, **31**(2), 114-122. https://doi.org/10.1002/minf.201100135

Hourcade, D. E., Mitchell, L. M., & Oglesby, T. J. (1999). Mutations of the type A domain of complement factor B that promote high-affinity C3b-binding. *J Immunology*, **162**, 2906-2911.

Hu, Z., Yu, C., Furutsuki, M., Andreoletti, G., Ly, M., Hoskins, R., ... Brenner, S. E. (2019). VIPdb, a genetic variant impact predictor database. *Human Mutation*, **40**, 1202-1214. https://doi.org/10.1002/humu.23858

Huntington, J. A. (2011). Serpin structure, function and dysfunction. *Journal of Thrombosis and Haemostasis*, **9**(Suppl. 1), 26-34. https://doi.org/10.1111/j.1538-7836.2011.04360.x

Huntington, J. A., Olson, S. T., Fan, B., & Gettins, P. G. (1996). Mechanism of heparin activation of Antithrombin. Evidence for reactive center loop pre-insertion with expulsion upon heparin binding. *Biochemistry*, **35**, 8495-8503.

Isvoran, A., Louet, M., Vladoiu, D. L., Craciun, D., Loriot, M.-A., Villoutreix, B. O., & Miteva, M. A. (2017). Pharmacogenomics of the cytochrome P450 2C family: Impacts of amino acid variations on drug metabolism. *Drug Discovery Today*, **22**(2), 366-376. https://doi.org/10.1016/j.drudis.2016.09.015

Ittisoponpisan, S., Islam, S. A., Khanna, T., Alhuzimi, E., David, A., & Sternberg, M. J. E. (2019). Can Predict Protein 3D Structures Provide Reliable Insights into whether Missense Variants Are Disease Associated? *Journal of Molecular Biology*, **431**, 2197-2212.

Izaguirre, G., Aguila, S., Qi, L., Swanson, R., Roth, R., Rezaie, A. R., ... Olson, S. T. (2014). Conformational activation of antithrombin by heparin involves an altered exosite interaction with protease. *Journal of Biological Chemistry*, **289**, 34049-34064. https://doi.org/10.1074/jbc.M114.611707

Jamroz, M., Kolinski, A., & Kihara, D. (2012). Structural features that predict real-value fluctuations of globular proteins. *Proteins: Structure, Function, and Bioinformatics*, **80**(5), 1425-1435. https://doi.org/10.1002/prot.24040

Jiménez-García, B., Pons, C., & Fernandez-Recio, J. (2013). pyDockWEB: A web server for rigid-body protein-protein docking using electrostatics and desolvation scoring. *Bioinformatics*, **29**(13), 1698-1699. https://doi.org/10.1093/bioinformatics/btt262

Jing, H., Xu, Y., Carson, M., Moore, D., Macon, K. J., Volanakis, J. E., & Narayana, S. V. (2000). New structural motifs on the chymotrypsin fold and their potential roles in complement factor B. *EMBO Journal*, **19**(2), 164-173. https://doi.org/10.1093/emboj/19.2.164

Johnson, D. J., Langdown, J., Li, W., Luis, S. A., Baglin, T. P., & Huntington, J. A. (2006). Crystal structure of monomeric native antithrombin reveals a novel reactive center loop conformation. *Journal of Biological Chemistry*, **281**, 35478-35486. https://doi.org/10.1074/jbc.M607204200

Karczewski, K. J., Francioli, L. C., Tiao, G., Cummings, B. B., Alföldi, J., Wang, Q., ... MacArthur, D. G. (2019). Variation across 141,456 Human Exomes and Genomes Reveals the Spectrum of loss-of-function Intolerance across Human protein-coding Genes. *bioRxiv*. 531210.

King, B. P., Khan, T. I., Aithal, G. P., Kamali, F., & Daly, A. K. (2004). Upstream and coding region CYP2C9 polymorphisms: Correlation with warfarin dose and metabolism. *Pharmacogenomics*, **14**(12), 813-822. https://doi.org/10.1097/00008571-200412000-00004

Kozakov, D., Grove, L. E., Hall, D. R., Bohnuud, T., Mottarella, S. E., Luo, L., ... Vajda, S. (2015). The FTMap family of web servers for determining and characterizing ligand-binding hot spots of proteins. *Nature Protocols*, **10**(5), 733-755. https://doi.org/10.1038/nprot.2015.043

Kozakov, D., Hall, D. R., Chuang, G.-Y., Cencic, R., Brenke, R., Grove, L. E., ... Vajda, S. (2011). Structural conservation of drugable hot spots in protein-protein interfaces. *Proceedings of the National Academy of Sciences*, **108**(33), 13528-13533. https://doi.org/10.1073/pnas.1101835108

Kozakov, D., Hall, D. R., Xia, B., Porter, K. A., Padhorny, D., Yueh, C., ... Vajda, S. (2017). The ClusPro web server for protein-protein
docking. *Nature Protocols*, **12**(2), 255–278. https://doi.org/10.1038/s41540-016-169

Kucukkal, T. G., Petukh, M., Li, L., & Alexov, E. (2015). Structural and physico-chemical effects of disease and non-disease nsSNPs on proteins. *Current Opinion in Structural Biology*, **32**, 18–24. https://doi.org/10.1016/j.sbi.2015.01.003

Kuznetsov, I. B., & McDuffie, M. (2008). FlexPred: A web-server for predicting residue positions involved in conformational switches in proteins. *Bioinformation*, **3**(3), 134–136. https://doi.org/10.6026/97320630003134

Laimer, J., Hiebl-Flach, J., Lengauer, D., & Lackner, P. (2016). MAESTROweb: A web server for structure-based protein stability prediction. *Bioinformatics*, **32**(9), 1414–1416. https://doi.org/10.1093/bioinformatics/btv769

Laimer, J., Hofer, H., Fritz, M., Wegenkittl, S., & Lackner, P. (2015). MAESTRO-multi agent stability prediction upon point mutations. *BMC Bioinformatics*, **16**, 116–128. https://doi.org/10.1186/s12859-015-0548-6

Landrum, M. J., Lee, J. M., Riley, G. R., Jang, W., Rubinstein, W. S., Church, D. M., & Maglott, D. R. (2014). ClinVar: Public archive of relationships among sequence variation and human phenotype. *Nucleic Acids Research*, **42**, D980–D985. https://doi.org/10.1093/nar/gkt1113

Lang, T., Klein, K., Richter, T., Zibat, A., Kerb, R., Eichelbaum, M., ... Zanger, U. M. (2004). Multiple novel nonsynonymous CYP2B6 gene polymorphisms in Caucasians: Demonstration of phenotypic null alleles. *Journal of Pharmacology and Experimental Therapeutics*, **311**(1), 34–43.

Laskowski, R. A., Stephenson, J. D., Sillitoe, I., Orengo, C. A., & Thornton, J. M. (2020). VarSite: Disease variants and protein structure. *Protein Science*, **29**(1), 111–119. https://doi.org/10.1002/pro.3746

Lee, M. Y., Borgiani, P., Johansson, I., Oteri, F., Mkrtchian, S., Falconi, M., & Ingelman-Sundberg, M. (2014). High warfarin sensitivity in carriers of CYP2C9*35 is determined by the impaired interaction with P450 oxidoreductase. *The Pharmacogenomics Journal*, **14**, 343–349.

Li, W., Johnson, D. J., Esmon, C. T., & Huntington, J. A. (2004). Structure of the antithrombin-thrombin-heparin ternary complex reveals the antithrombotic mechanism of heparin. *Nature Structural & Molecular Biology*, **11**, 857–862. https://doi.org/10.1038/nsmb811

Liang, S., Zhang, C., Liu, S., & Zhou, Y. (2006). Protein binding site prediction using an empirical scoring function. *Nucleic Acids Research*, **34**(13), 3698–3707. https://doi.org/10.1093/nar/gkl454

Liu, M.-L., Shen, B. W., Nakaya, S., Pratt, K. P., Fujikawa, K., Davie, E. W., ... Thompson, A. R. (2000). Hemophilic factor VIII C1- and C2-domain missense mutations and their modeling to the 1.5-angstrom human C2-domain crystal structure. *Blood*, **96**, 979–987. https://doi.org/10.1182/blood.V96.3.979

Lomize, M. A., Pogozheva, I. D., Joo, H., Mosberg, H. I., & Lomize, A. L. (2012). OPM database and PPM web server: Resources for positioning of proteins in membranes. *Nucleic Acids Research*, **40**, D370–D376. https://doi.org/10.1093/nar/gkr703

Louet, M., Labbe, C. M., Fagnen, C., Aono, C. M., Homem-de-Mello, P., Villoutreix, B. O., & Miteva, M. A. (2018). Insights into molecular mechanisms of drug metabolism dysfunction of human CYP2C9*30. *PLoS ONE*, **13**(5), e0197249. https://doi.org/10.1371/journal.pone.0197249

Luxembourg, B., Delev, D., Geisen, C., Spannagl, M., Krause, M., Miesbach, W., ... Pavlova, A. (2011). Molecular basis of antithrombin deficiency. *Thrombosis and Haemostasis*, **105**, 635–646.

Luxembourg, B., D'Souza, M., Körber, S., & Seifried, E. (2015). Prediction of the pathogenicity of antithrombin sequence variations by in silico methods. *Thrombosis Research*, **135**(2), 404–409. https://doi.org/10.1016/j.thromres.2014.11.022

Marinozzi, M. C., Vergoz, L., Rybkine, T., Ngo, S., Bettoni, S., Pashov, A., ... Roumenina, L. T. (2014). Complement factor b mutations in atypical hemolytic uremic syndrome-disease-relevant or benign? *Journal of the American Society of Nephrology*, **25**, 2053–2065.

Martin, A. J. M., Vidotto, M., Boscariol, F., Di Domenico, T., Walsh, I., & Tosatto, S. C. E. (2011). RING: Networking interacting residues, evolutionary information and energetics in protein structures. *Bioinformatics*, **27**, 2003–2005. https://doi.org/10.1093/bioinformatics/btr191

Martiny, V. Y., & Miteva, M. A. (2013). Advances in molecular modeling of human cytochrome P450 polymorphism. *Journal of Molecular Biology*, **425**(21), 3978–3992. https://doi.org/10.1016/j.jmb.2013.07.010

Merle, N. S., Church, S. E., Fremeaux-Bacchi, V., & Roumenina, L. T. (2015). Complement System Part I—Molecular mechanisms of activation and regulation. *Frontiers in Immunology*, **6**, Article 262. https://doi.org/10.3389/fimmu.2015.00262

Merle, N. S., Noe, R., Halbwachs-Mecarelli, L., Fremeaux-Bacchi, V., & Roumenina, L. T. (2015). Complement System Part II: Role in immunity. *Frontiers in Immunology*, **6**, Article 257. https://doi.org/10.3389/fimmu.2015.00257

Milder, F. J., Gomes, L., Schouten, A., Janssen, B. J., Huizinga, E. G., Romijn, R. A., ... Gros, P. (2007). Factor B structure provides insights into activation of the central protease of the complement system. *Nature Structural & Molecular Biology*, **14**, 224–228. https://doi.org/10.1038/nsmb1210

Mottarella, S. E., Beglov, D., Beglova, N., Nugent, M. A., Kozakov, D., & Vajda, S. (2014). Docking server for the identification of heparin binding sites on proteins. *Journal of Chemical Information and Modeling*, **54**(7), 2068–2078. https://doi.org/10.1021/ci500115j

Muszbek, L., Bereczky, Z., Kovács, B., & Komáromi, I. (2010). Antithrombin deficiency and its laboratory diagnosis. *Clinical Chemistry and Laboratory Medicine*, **48**(Suppl 1), S67–S78. https://doi.org/10.1515/CCLM.2010.368

Neuvirth, H., Raz, R., & Schreiber, G. (2004). ProMate: A structure based prediction program to identify the location of protein–protein binding sites. *Journal of Molecular Biology*, **338**(1), 181–199. https://doi.org/10.1016/j.jmb.2004.02.040

Ngan, C. H., Bohnuud, T., Mottarella, S. E., Beglov, D., Villar, E. A., Hall, D. R., ... Vajda, S. (2012). FTMAP: Extended protein mapping with user-selected probe molecules. *Nucleic Acids Research*, **40**(W1), W271–W275. https://doi.org/10.1093/nar/gks441

Olson, S. T., Richard, B., Izaguirre, G., Schedin-Weiss, S., & Gettins, P. G. (2010). Molecular mechanisms of antithrombin-heparin regulation of blood clotting proteinases. A paradigm for understanding proteinase regulation by serpin family protein proteinase inhibitors. *Biochimie*, **92**, 1587–1596. https://doi.org/10.1016/j.biochi.2010.05.011

Palta, S., Saroa, R., & Palta, A. (2014). Overview of the coagulation system. *Indian Journal of Anaesthesia*, **58**(5), 515–523. https://doi.org/10.4103/0019-5049.144643

Pandurangan, A. P., Ochoa-Montaño, B., Ascher, D. B., & Blundell, T. L. (2017). SDM: A server for predicting effects of mutations on protein stability. *Nucleic Acids Research*, **45**(W1), W229–W235. https://doi.org/10.1093/nar/gkx439
Patnaik, M. M., & Moll, S. (2008). Inherited antithrombin deficiency. *Haemophilia, **14***, 1229-1239.

Pedersen, D. V., Roumenina, L., Jensen, R. K., Gadeberg, T. A., Marinozzi, C., Picard, C., ... Andersen, G. R. (2017). Functional and structural insight into properdin control of complement alternative pathway amplification. *EMBO Journal, **36***, 1084-1099. https://doi.org/10.15252/embj.201696173

Peng, Y., & Alexov, E. (2016). Investigating the linkage between disease-causing amino acid variants and their effect on protein stability and binding. *Proteins, **84***, 232-239. https://doi.org/10.1002/prot.24968

Pettersen, E. F., Goddard, T. D., Huang, C. C., Couch, G. S., Greenblatt, D. M., Meng, E. C., & Ferrin, T. E. (2004). UCSF Chimera—A visualization system for exploratory research and analysis. *Journal of Computational Chemistry, **25***, 1605-1612. https://doi.org/10.1002/jcc.20084

Petukh, M., Kucukkal, T. G., & Alexov, E. (2015). On human disease-causing amino acid variants: Statistical study of sequence and structural patterns. *Human Mutation, **36***, 524-534. https://doi.org/10.1002/humu.22770

Phillips, J. C., Braun, R., Wang, W., Gumbart, J., Tajkhorshid, E., Villa, E., ... Schulten, K. (2005). Scalable molecular dynamics with NAMD. *Journal of Computational Chemistry, **26***, 1781-1802. https://doi.org/10.1002/jcc.20289

Piovesan, D., Minervini, G., & Tosatto, S. C. E. (2016). The RING 2.0 web server for high quality residue interaction networks. *Nucleic Acids Research, **44*** (W1), W367-W374.

Pires, D. E., Ascher, D. B., & Blundell, T. L. (2014a). DUET: A server for predicting effects of mutations on protein stability via an integrated computational approach. *Nucleic Acids Research, **42*** (W1), W314-W319.

Pires, D. E., Ascher, D. B., & Blundell, T. L. (2014b). mCSM: Predicting the effects of mutations in proteins using graph-based signatures. *Bioinformatics, **30*** (3), 335-342. https://doi.org/10.1093/bioinformatics/btt691

Ponnuraj, K., Xu, Y., Macon, K., Moore, D., Volanakis, J. E., & Narayana, S. V. (2004). Structural analysis of engineered Bb fragment of complement factor B: Insights into the activation mechanism of the alternative pathway C3-convertase. *Molecular Cell, **14*** (1), 17-28. https://doi.org/10.1016/S1097-2765(04)00160-1

Pratt, K. P., Shen, B. W., Takeshima, K., Davie, E. W., Fujikawa, K., & Stoddard, B. L. (1999). Structure of the C2 domain of human factor VIII at 1.5 Å resolution. *Nature, **402***, 439-442.

Pravda, L., Berka, K., Svobodová Vařeková, R., Sehnal, D., Banáš, P., Laskowski, R. A., ... Otyepka, M. (2014). Anatomy of enzyme channels. *BMC Bioinformatics, **15***, 379-386. https://doi.org/10.1186/s12859-014-0379-x

Pravda, L., Sehnal, D., Svobodová Vařeková, R., Navratilova, V., Tousek, D., Berka, K., ... Koča, J. (2018). ChannelsDB: Database of biomacromolecular tunnels and pores. *Nucleic Acids Research, **46*** (D1), D399-D405. https://doi.org/10.1093/nar/gkx868

Qin, S. B., & Zhou, H.-X. (2007). meta-PPISP: A meta web server for protein-protein interaction site prediction. *Bioinformatics, **23***, 3386-3387. https://doi.org/10.1093/bioinformatics/btm434

Ramamoorthy, A., & Skaar, T. C. (2011). In silico identification of microRNAs predicted to regulate the drug metabolizing cytochrome P450 genes. *Drug Metabolism Letters, **5***, 126-131.

Rooijakkers, S. H., Wu, J., Ruyken, M., van Domselaar, R., Planken, K. L., Tzekou, A., ... Gros, P. (2009). Structural and functional implications of the alternative complement pathway C3 convertase

stabilized by a staphylococcal inhibitor. *Nature Immunology, **10*** (7), 721-727. https://doi.org/10.1038/ni.1756

Rotger, M., Tegude, H., Colombo, S., Cavassini, M., Furrer, H., Decosterd, L., ... Zanger, U. M. (2007). Predictive value of known and novel alleles of CYP2B6 for efavirenz plasma concentrations in HIV-infected individuals. *Clinical Pharmacology and Therapeutics, **81*** (4), 557-566. https://doi.org/10.1038/sj.clpt.6100072

Roumenina, L. T., Jablonski, M., Hue, C., Blouin, J., Dimitrov, J. D., Dragon-Durey, M. A., ... Fremeaux-Bacchi, V. (2009). Hyperfunctional C3 convertase leads to complement deposition on endothelial cells and contributes to atypical hemolytic uremic syndrome. *Blood, **114***, 2837-2845. https://doi.org/10.1182/blood-2009-01-197640

Rydberg, P., & Olsen, L. (2012). Predicting drug metabolism by cytochrome P450 2C9: Comparison with the 2D6 and 3A4 isoforms. *ChemMedChem, **7*** (7), 1202-1209. https://doi.org/10.1002/chemmed.20084

Sansen, S., Yano, J. K., Reynald, R. L., Schoch, G. A., Griffin, K. J., Stout, C. D., & Johnson, E. F. (2007). Adaptations for the oxidation of polycyclic aromatic hydrocarbons exhibited by the structure of human P450 1A2. *Journal of Biological Chemistry, **282***, 14348-14355. https://doi.org/10.1074/jbc.M611692200

Shah, M. B., Talakad, J. C., Maekawa, K., Roberts, A. G., Wilderman, P. R., Sun, L., ... Halpert, J. R. (2010). Crystal structure of a cytochrome P450 2B6 genetic variant in complex with the inhibitor 4-(4-chlorophenyl)imidazole at 2.0-A resolution. *Molecular Pharmacology, **77***, 529-538.

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., ... Ideker, T. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. *Genome Research, **13*** (11), 2498-2504. https://doi.org/10.1101/gr.1239303

Shen, B. W., Spiegel, P. C., Chang, C. H., Huh, J. W., Lee, J. S., Kim, J., ... Stoddard, B. L. (2008). The tertiary structure and domain organization of coagulation factor VIII. *Blood, **111***, 1240-1247. https://doi.org/10.1182/blood-2007-08-109918

Shintani, M., Ieiri, I., Inoue, K., Mamiya, K., Ninomiya, H., Tashiro, N., ... Otsubo, K. (2001). Genetic polymorphisms and functional characterization of the 5'-flanking region of the human CYP2C9 gene: In vitro and in vivo studies. *Clinical Pharmacology and Therapeutics, **70*** (2), 175-182. https://doi.org/10.1067/mcp.2001.117367

Sievers, F., Wilm, A., Dineen, D., Gibson, T. J., Karplus, K., Li, W., ... Higgins, D. G. (2011). Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. *Molecular Systems Biology, **7***, 539-544. https://doi.org/10.1038/msb.2011.75

Sim, S. C., & Ingelman-Sundberg, M. (2013). Update on allele nomenclature for human cytochromes P450 and the Human Cytochrome P450 Allele (CYP-Allele) Nomenclature Database. *Methods in Molecular Biology, **987***, 251-259.

Simonetti, L., Bruque, C. D., Fernández, C. S., Benavides-Mori, B., Delea, M., Kolomenski, J. E., ... Dain, L. (2018). CYP21A2 mutation update: Comprehensive analysis of databases and published genetic variants. *Human Mutation, **39***, 61-68.

Singh, D., Kashyap, A., Pandey, R. V., & Saini, K. S. (2011). Novel advances in cytochrome P450 research. *Drug Discov Today, **16***, 793-799. https://doi.org/10.1016/j.drudis.2011.08.003

Spiegel, P. C. Jr, Jacquemin, M., Saint-Remy, J.-M.-R., Stoddard, B. L., & Pratt, K. P. (2001). Structure of a factor VIII C2 domain-immunoglobulin G4k Fab complex: Identification of an inhibitory antibody epitope on the surface of factor VIII. *Blood, **98*** (1), 13-19.
Spiegel, P. C., Murphy, P., & Stoddard, B. L. (2004). Surface-exposed hemophilic mutations across the factor VIII C2 domain have variable effects on stability and binding activities. *Journal of Biological Chemistry, 279*(51), 53691–53698. https://doi.org/10.1074/jbc. M409389200

Stenson, P. D., Mort, M., Ball, E. V., Shaw, K., Phillips, A. D., & Cooper, D. N. (2014). The human gene mutation database: Building a comprehensive mutation repository for clinical and molecular genetics, diagnostic testing and personalized genomic medicine. *Human Genetics, 133*, 1–9. https://doi.org/10.1007/s00439-013-1358-4

Takano, K., Liu, D., Tarpey, P., Gallant, E., Lam, A., Witham, S., ... Dulhunty, A. F. (2012). An X-linked channelopathy with cardio-megaly due to a CLIC2 mutation enhancing ryanodine receptor channel activity. *Human Molecular Genetics, 21*, 4497–4507. https://doi.org/10.1093/hmg/dds292

Thusberg, J., & Vihinen, M. (2009). Pathogenic or not? And if so, then how? Studying of the effect of missense mutations using bioinformatics methods. *Human Mutation, 30*, 703–714. https://doi.org/10.1002/humu.20938

UniProt Consortium. (2019). Uniprot: A worldwide hub of protein knowledge. *Nucleic Acids Research, 47*, D506–D515.

Versteeg, H. H., Heemskerk, J. W. M., Levi, M., & Reitsma, P. H. (2013). New fundamentals in hemostasis. *Physiological Reviews, 93*(1), 327–358. https://doi.org/10.1152/physrev.00016.2011

Villoutreix, B. O. (2002). Structural bioinformatics: Methods, concepts and applications to blood coagulation proteins. *Current Protein and Peptide Science, 3*, 341–364. https://doi.org/10.2174/1389203023 380657

Villoutreix, B. O., Lagorce, D., Labbé, C. M., Sperandio, O., & Miteva, M. A. (2013). One hundred thousand mouse clicks down the road: Selected online resources supporting drug discovery collected over a decade. *Drug Discov Today, 18*, 1081–1089. https://doi.org/10.1016/j.drudis.2013.06.013

Villoutreix, B. O., & Miteva, M. A. (2016). Discoidin domains as emerging therapeutic targets. *Trends in Pharmacological Sciences, 37*(8), 641–659. https://doi.org/10.1016/j.tips.2016.06.003

Williams, P. A., Cosme, J., Ward, A., Angove, H. C., Matak Vinkovic, D., & Jhoti, H. (2003). Crystal structure of human cytochrome P450 2C9 with bound warfarin. *Nature, 424*, 464–468. https://doi.org/10.1038/nature01862

Witham, S., Takano, K., Schwartz, C., & Alexov, E. (2011). A missense mutation in CLIC2 associated with intellectual disability is predicted by in silico modeling to affect protein stability and dynamics. *Proteins, 79*, 2444–2454. https://doi.org/10.1002/prot.23065

Yates, C. M., & Sternberg, M. J. (2013). The effects of non-synonymous single nucleotide polymorphisms (nsSNPs) on protein–protein interactions. *Journal of Molecular Biology, 425*, 3949–3963. https://doi.org/10.1016/j.jmb.2013.07.012

Yiannakopoulou, E. C. (2013). Pharmacogenomics for individualized therapies. In D. Barh, V. Zambare & V. Azevedo (Eds.), *OMICS-Applications in biomedical, agricultural and enviromental sciences* (pp. 95–120). Boca Raton, FL: CRC Press.

Yin, T., Maekawa, K., Kamide, K., Saito, Y., Hanada, H., Miyashita, K., ... Miyata, T. (2008). Genetic variations of CYP2C9 in 724 Japanese individuals and their impact on the antihypertensive effects of losartan. *Hypertension Research, 31*, 1549–1557. https://doi.org/10.1291/hypres.31.1549

Zhou, H., Josephy, P. D., Kim, D., & Guengerich, F. P. (2004). Functional characterization of four allelic variants of human cytochrome P450 1A2. *Archives of Biochemistry and Biophysics, 422*(1), 23–30. https://doi.org/10.1016/j.abb.2003.11.019

Zhou, H.-X., & Shan, Y. (2001). Prediction of protein interaction sites from sequence profiles and residue neighbor list. *Proteins, 44*, 336–343.

Zhou, Y., Ingelman-Sundberg, M., & Lauschke, V. M. (2017). Worldwide Distribution of Cytochrome P450 Alleles: A meta-analysis of population-scale sequencing projects. *Cliniacal Pharmacology and Therapeutics, 102*(4), 688–700. https://doi.org/10.1002/cpt.690

## SUPPORTING INFORMATION

Additional supporting information may be found online in the Supporting Information section.

**How to cite this article:** Gyulkhandanyan A, Rezaie AR, Roumenina L, et al. Analysis of protein missense alterations by combining sequence- and structure-based methods. *Mol Genet Genomic Med.* 2020;8:e1166. https://doi.org/10.1002/mgg3.1166