Computational and artificial neural network based study of functional SNPs of human LEPR protein associated with reproductive function

Genetic polymorphisms are mostly associated with inherited diseases, detecting and analyzing the biological significance of functional single-nucleotide polymorphisms (SNPs) using wet laboratory experiments is an arduous task hence the computational analysis of putative SNPs is essential before conducting a study on a large population. SNP in the leptin receptor (LEPR) could result in the retention of intracellular signalling due to the structural and functional instability of the receptor causing abnormal reproductive function in human. In this first comprehensive computational analysis of LEPR gene mutation, we have identified and analyzed the functional consequence and structural significance of the SNPs in LEPR using recently developed several computational algorithms. Thirteen deleterious mutations such as W13C, S93G, I232R, Q307H, Y354C, E497A, Q571H, R612H, K656N, T690A, T699M V741M, and L760R were identified in the LEPR gene coding region. Backpropagation algorithm has been developed to forestall the deleterious nature of SNP and to validate the outcome of the tested computational tools. From ConSurf prediction three SNPs (Q571H, R612H, and T699M) were highly conserved on LEPR protein and the most deleterious variant R612H had one hydrogen bond abolished and severely reduced protein stability. Molecular docking suggested that the mutant (R612H) LEPR had lowest binding energy than native LEPR with the ligand molecule. Thus the energetically destructive changeover of ARG to HIS in R612H could possibly affect the LEPR protein structural stability and functional constancy due to interruption in the amino acid interactions and could result in reproductive disorders in human and increases the complication in obstetric and pregnancy outcome.     

Genetic polymorphisms and haplotype structures of the CYP4A22 gene in a Japanese population

The CYP4A fatty acid monooxygenases oxidize endogenous arachidonic acid to 20-hydroxyeicosatetraenoic acid that acts as a regulator of blood pressure. Among the isoforms of the CYP4A subfamily, the human CYP4A22 was recently identified. In this study, we report the comprehensive investigation of polymorphisms in the CYP4A22 gene. To investigate genetic variation in CYP4A22 in 191 Japanese subjects, we used denaturing HPLC (DHPLC) and direct sequencing. Our investigation has enabled the identification of 13 sequence variations in the CYP4A22 coding region, thereby demonstrating for the first time that this gene is subject to polymorphism. Two of these sequence variations correspond to silent mutations located in exons 8 (His323His) and 9 (Gly390Gly). Nine of these sequence variations correspond to missense mutations located in exons 1 (Arg11Cys), 3 (Arg126Trp), 4 (Gly130Ser and Asn152Tyr), 5 (Val185Phe), 6 (Cys231Arg), 7 (Lys276Thr), 10 (Leu428Pro), and 12 (Leu509Phe). One of these sequence variations corresponds to nonsense mutations located in exon 9 (Gln368stop). The 13th mutation corresponds to a nucleotide deletion (G7067del) that causes a frameshift and consequently results in a stop codon 80 nucleotides downstream. In addition to the wild-type CYP4A22*1 allele, 20 variants, namely CYP4A22*2-15, were characterized by haplotype analysis. Based on these data, we concluded that allelic variants of the human CYP4A22 gene exist and speculated that some of these variants may be functionally relevant.     

Distinct functional consequences of MUTYH variants associated with colorectal cancer: Damaged DNA affinity,glycosylase activity and interaction with PCNA and Hus1

MUTYH is a base excision repair (BER) enzyme that prevents mutations in DNA associated with 8-oxoguanine (OG) by catalyzing the removal of adenine from inappropriately formed OG:A base-pairs. Germline mutations in the MUTYH gene are linked to colorectal polyposis and a high risk of colorectal cancer, a syndrome referred to as MUTYH-associated polyposis (MAP). There are over 300 different MUTYH mutations associated with MAP and a large fraction of these gene changes code for missense MUTYH variants. Herein, the adenine glycosylase activity, mismatch recognition properties, and interaction with relevant protein partners of human MUTYH and five MAP variants (R295C, P281L, Q324H, P502L, and R520Q) were examined. P281L MUTYH was found to be severely compromised both in DNA binding and base excision activity, consistent with the location of this variation in the iron-sulfur cluster (FCL) DNA binding motif of MUTYH. Both R295C and R520Q MUTYH were found to have low fractions of active enzyme, compromised affinity for damaged DNA, and reduced rates for adenine excision. In contrast, both Q324H and P502L MUTYH function relatively similarly to WT MUTYH in both binding and glycosylase assays. However, P502L and R520Q exhibited reduced affinity for PCNA (proliferation cell nuclear antigen), consistent with their location in the PCNA-binding motif of MUTYH. Whereas, only Q324H, and not R295C, was found to have reduced affinity for Hus1 of the Rad9–Hus1–Rad1 complex, despite both being localized to the same region implicated for interaction with Hus1. These results underscore the diversity of functional consequences due to MUTYH variants that may impact the progression of MAP.     

Substitutions at amino acid positions 143, 148, and 155 of HIV-1 integrase define distinct genetic barriers to raltegravir resistance in vivo

Mutations at amino acids 143, 148, and 155 in HIV-1 integrase (IN) define primary resistance pathways in subjects failing raltegravir (RAL)-containing treatments. Although each pathway appears to be genetically distinct, shifts in the predominant resistant virus population have been reported under continued drug pressure. To better understand this dynamic, we characterized the RAL susceptibility of 200 resistant viruses, and we performed sequential clonal analysis for selected cases. Patient viruses containing Y143R, Q148R, or Q148H mutations consistently exhibited larger reductions in RAL susceptibility than patient viruses containing N155H mutations. Sequential analyses of virus populations from three subjects revealed temporal shifts in subpopulations representing N155H, Y143R, or Q148H escape pathways. Evaluation of molecular clones isolated from different time points demonstrated that Y143R and Q148H variants exhibited larger reductions in RAL susceptibility and higher IN-mediated replication capacity (RC) than N155H variants within the same subject. Furthermore, shifts from the N155H pathway to either the Q148R or H pathway or the Y143R pathway were dependent on the amino acid substitution at position 148 and the secondary mutations in Y143R- or Q148R- or H-containing variants and correlated with reductions in RAL susceptibility and restorations in RC. Our observations in patient viruses were confirmed by analyzing site-directed mutations. In summary, viruses that acquire mutations defining the 143 or 148 escape pathways are less susceptible to RAL and exhibit greater RC than viruses containing 155 pathway mutations. These selective pressures result in the displacement of N155H variants by 143 or 148 variants under continued drug exposure.     

Interactions of the vitamin D receptor with the corepressor hairless: analysis of hairless mutants in atrichia with papular lesions

Atrichia with papular lesions (APL) and hereditary vitamin D-resistant rickets have a similar congenital hair loss disorder caused by mutations in hairless (HR) and vitamin D receptor (VDR) genes, respectively. HR is a VDR corepressor, and it has been hypothesized that VDR.HR suppress gene expression during specific phases of the hair cycle. In this study, we examined the corepressor activity of HR mutants (E583V, C622G, N970S, V1056M, D1012N, V1136D, and Q1176X) previously described as the molecular cause of APL as well as HR variants (P69S, C397Y, A576V, E591G, R620Q, T1022A) due to non-synonymous polymorphisms in the HR gene. We found that the corepressor activities of all but one of the pathogenic HR mutants were completely abolished. HR mutant E583V exhibited normal corepressor activity, suggesting that it may not be pathogenic. In co-immunoprecipitation assays, all of the pathogenic HR mutants bound VDR but exhibited reduced binding to histone deacetylase 1 (HDAC1), suggesting that the impaired corepressor activity is due in part to defective interactions with HDACs. The HR variants exhibited two classes of corepressor activity, those with normal activity (C397Y, E591G, R620Q) and those with partially reduced activity (P69S, A576V, T1022A). All of the variants interacted with VDR and HDAC1 with the exception of P69S, which was degraded. When coexpressed with VDR, all of the HR pathogenic mutants and variants increased the level of VDR protein, demonstrating that this function of HR was not impaired by these mutations. This study of HR mutations provides evidence for the molecular basis of APL.     

Vaccine-escape and fast-growing mutations in the United Kingdom,the United States,Singapore, Spain,India, and other COVID-19-devastated countries

Recently, the SARS-CoV-2 variants from the United Kingdom (UK), South Africa, and Brazil have received much attention for their increased infectivity, potentially high virulence, and possible threats to existing vaccines and antibody therapies. The question remains if there are other more infectious variants transmitted around the world. We carry out a large-scale study of 506,768 SARS-CoV-2 genome isolates from patients to identify many other rapidly growing mutations on the spike (S) protein receptor-binding domain (RBD). We reveal that essentially all 100 most observed mutations strengthen the binding between the RBD and the host angiotensin-converting enzyme 2 (ACE2), indicating the virus evolves toward more infectious variants. In particular, we discover new fast-growing RBD mutations N439K, S477N, S477R, and N501T that also enhance the RBD and ACE2 binding. We further unveil that mutation N501Y involved in United Kingdom (UK), South Africa, and Brazil variants may moderately weaken the binding between the RBD and many known antibodies, while mutations E484K and K417N found in South Africa and Brazilian variants, L452R and E484Q found in India variants, can potentially disrupt the binding between the RBD and many known antibodies. Among these RBD mutations, L452R is also now known as part of the California variant B.1.427. Finally, we hypothesize that RBD mutations that can simultaneously make SARS-CoV-2 more infectious and disrupt the existing antibodies, called vaccine escape mutations, will pose an imminent threat to the current crop of vaccines. A list of most likely vaccine escape mutations is given, including S494P, Q493L, K417N, F490S, F486L, R403K, E484K, L452R, K417T, F490L, E484Q, and A475S. Mutation T478K appears to make the Mexico variant B.1.1.222 the most infectious one. Our comprehensive genetic analysis and protein-protein binding study show that the genetic evolution of SARS-CoV-2 on the RBD, which may be regulated by host gene editing, viral proofreading, random genetic drift, and natural selection, gives rise to more infectious variants that will potentially compromise existing vaccines and antibody therapies.     

Molecular model of human CYP21 based on mammalian CYP2C5: structural features correlate with clinical severity of mutations causing congenital adrenal hyperplasia

Enhanced understanding of structure-function relationships of human 21-hydroxylase, CYP21, is required to better understand the molecular causes of congenital adrenal hyperplasia. To this end, a structural model of human CYP21 was calculated based on the crystal structure of rabbit CYP2C5. All but two known allelic variants of missense type, a total of 60 disease-causing mutations and six normal variants, were analyzed using this model. A structural explanation for the corresponding phenotype was found for all but two mutants for which available clinical data are also discrepant with in vitro enzyme activity. Calculations of protein stability of modeled mutants were found to correlate inversely with the corresponding clinical severity. Putative structurally important residues were identified to be involved in heme and substrate binding, redox partner interaction, and enzyme catalysis using docking calculations and analysis of structurally determined homologous cytochrome P450s (CYPs). Functional and structural consequences of seven novel mutations, V139E, C147R, R233G, T295N, L308F, R366C, and M473I, detected in Scandinavian patients with suspected congenital adrenal hyperplasia of different severity, were predicted using molecular modeling. Structural features deduced from the models are in good correlation with clinical severity of CYP21 mutants, which shows the applicability of a modeling approach in assessment of new CYP21 mutations.     

Structure-function analysis of Friedreich's ataxia mutants reveals determinants of frataxin binding and activation of the Fe-S assembly complex

Friedreich's ataxia (FRDA) is a progressive neurodegenerative disease associated with the loss of function of the protein frataxin (FXN) that results from low FXN levels due to a GAA triplet repeat expansion or, occasionally, from missense mutations in the FXN gene. Here biochemical and structural properties of FXN variants, including three FRDA missense mutations (N146K, Q148R, and R165C) and three related mutants (N146A, Q148G, and Q153A), were determined in an effort to understand the structural basis for the loss of function. In vitro assays revealed that although the three FRDA missense mutations exhibited similar losses of cysteine desulfurase and Fe-S cluster assembly activities, the causes for these activation defects were distinct. The R165C variant exhibited a k(cat)/K(M) higher than that of native FXN but weak binding to the NFS1, ISD11, and ISCU2 (SDU) complex, whereas the Q148R variant exhibited the lowest k(cat)/K(M) of the six tested FXN variants and only a modest binding deficiency. The order of the FXN binding affinities for the SDU Fe-S assembly complex was as follows: FXN > Q148R > N146A > Q148G > N146K > Q153A > R165C. Four different classes of FXN variants were identified on the basis of their biochemical properties. Together, these structure-function studies reveal determinants for the binding and allosteric activation of the Fe-S assembly complex and provide insight into how FRDA missense mutations are functionally compromised.     

The complete amino acid sequence of coagulogen isolated from Southeast Asian horseshoe crab, Carcinoscorpius rotundicauda

The complete amino acid sequence of coagulogen purified from the hemocytes of the horseshoe crab Carcinoscorpius rotundicauda was determined by characterization of the NH2-terminal sequence and the peptides generated after digestion of the protein with lysyl endopeptidase, Staphylococcal aureus protease V8 and trypsin. Upon sequencing the peptides by the automated Edman method, the following sequence was obtained: A D T N A P L C L C D E P G I L G R N Q L V T P E V K E K I E K A V E A V A E E S G V S G R G F S L F S H H P V F R E C G K Y E C R T V R P E H T R C Y N F P P F V H F T S E C P V S T R D C E P V F G Y T V A G E F R V I V Q A P R A G F R Q C V W Q H K C R Y G S N N C G F S G R C T Q Q R S V V R L V T Y N L E K D G F L C E S F R T C C G C P C R N Y Carcinoscorpius coagulogen consists of a single polypeptide chain with a total of 175 amino acid residues and a calculated molecular weight of 19,675. The secondary structure calculated by the method of Chou and Fasman reveals the presence of an alpha-helix region in the peptide C segment (residue Nos. 19 to 46), which is released during the proteolytic conversion of coagulogen to coagulin gel. The beta-sheet structure and the 16 half-cystines found in the molecule appear to yield a compact protein stable to acid and heat. The amino acid sequences of coagulogen of four species of limulus have been compared and the interspecies evolutionary differences are discussed.     

Selective modulation of protein C affinity for EPCR and phospholipids by Gla domain mutation

Uniquely amongst vitamin K-dependent coagulation proteins, protein C interacts via its Gla domain both with a receptor, the endothelial cell protein C receptor (EPCR), and with phospholipids. We have studied naturally occurring and recombinant protein C Gla domain variants for soluble (s)EPCR binding, cell surface activation to activated protein C (APC) by the thrombin-thrombomodulin complex, and phospholipid dependent factor Va (FVa) inactivation by APC, to establish if these functions are concordant. Wild-type protein C binding to sEPCR was characterized with surface plasmon resonance to have an association rate constant of 5.23 x 10(5) m(-1).s(-1), a dissociation rate constant of 7.61 x 10(-2) s(-1) and equilibrium binding constant (K(D)) of 147 nm. It was activated by thrombin over endothelial cells with a K(m) of 213 nm and once activated to APC, rapidly inactivated FVa. Each of these interactions was dramatically reduced for variants causing gross Gla domain misfolding (R-1L, R-1C, E16D and E26K). Recombinant variants Q32A, V34A and D35A had essentially normal functions. However, R9H and H10Q/S11G/S12N/D23S/Q32E/N33D/H44Y (QGNSEDY) variants had slightly reduced (< twofold) binding to sEPCR, arising from an increased rate of dissociation, and increased K(m) (358 nm for QGNSEDY) for endothelial cell surface activation by thrombin. Interestingly, these variants had greatly reduced (R9H) or greatly enhanced (QGNSEDY) ability to inactivate FVa. Therefore, protein C binding to sEPCR and phospholipids is broadly dependent on correct Gla domain folding, but can be selectively influenced by judicious mutation.     

Identification and characterization of a novel variant of the human P2X<Subscript>7</Subscript> receptor resulting in gain of function

The P2X₇ receptor exhibits significant allelic polymorphism in humans, with both loss and gain of function variants potentially impacting on a variety of infectious and inflammatory disorders. At least five loss-of-function polymorphisms (G150R, R307Q, T357S, E496A, and I568N) and two gain-of-function polymorphisms (H155Y and Q460R) have been identified and characterized to date. In this study, we used RT-PCR cloning to isolate and characterize P2X₇ cDNA clones from human PBMCs and THP-1 cells. A previously unreported variant with substitutions of V80M and A166G was identified. When expressed in HEK293 cells, this variant exhibited heightened sensitivity to the P2X₇ agonist (BzATP) relative to the most frequent allele, as shown by pore formation measured by fluorescent dye uptake into cells. Mutational analyses showed that A166G alteration was critical for the gain-of-function change, while V80M was not. Full-length variants with multiple previously identified nonsynonymous SNPs (H155Y, H270R, A348T, and E496A) were also identified. Distinct functional phenotypes of the P2X₇ variants or mutants constructed with multiple polymorphisms were observed. Gain-of-function variations (A166G or H155Y) could not rescue the loss-of-function E496A polymorphism. Synergistic effects of the gain-of-function variations were also observed. We also identified the A348T alteration as a weak gain-of-function variant. Thus, these results identify the new gain-of-function variant A166G and demonstrate that multiple-gene polymorphisms contribute to functional phenotypes of the human P2X₇ receptor. Furthermore, the results demonstrate that the C-terminal of the cysteine-rich domain 1 of P2X₇ is critical for regulation of P2X₇-mediated pore formation.     

Molecular cloning and characterization of the circumsporozoite protein gene of Plasmodium inui isolated from Formosan macaques (Macaca cyclopis) in Taiwan

We characterized the complete nucleic and amino acid sequences of the Plasmodium inui circumsporozoite protein (Pincsp) gene and analyzed nucleotide diversity across the entire Pincsp gene by using 7 field isolates and strains Taiwan I and II obtained from Formosan macaques (Macaca cyclopis) in Taiwan. The length of the circumsporozoite protein ( CSP ) gene ranged from 1077 to 1125 bp. Size polymorphisms were due to variations in the number of tandem repeat units. The non-repetitive (NR) region exhibited high homology (99.1 ～ 100 and 98.7 ～ 100% at the nucleotide and amino acid levels, respectively) and was conserved among the variants (nucleotide diversities, π, of the 5'NR and 3'NR regions were 0.00364 and 0.00392, respectively). In the central repetitive (CR) region, we decomposed the sequences into 2 kinds of repeating amino acid motifs, i.e., a repeat unit R1, PA(P/A)(P/A)A(E)GG (n = 11-13), and a following repeat unit R2: P(A/G)(A/P/G)(P/Q)AQ(N/K) (n = 9-10). Analyzing these repeat sequences showed evidence of 3 genetic mechanisms for generating variations in the repeats of the Pincsp gene, i.e., point mutation, insertion, and recombination. These findings suggest that polymorphisms in the Pincsp gene are essentially limited to the CR region, which showed much greater variability in terms of length, number of repeats, and sequence.     

Structural basis of thermostability. Analysis of stabilizing mutations in subtilisin BPN'

The crystal structures of two thermally stabilized subtilisin BPN' variants, S63 and S88, are reported here at 1.8 and 1.9 A resolution, respectively. The micromolar affinity calcium binding site (site A) has been deleted (Delta75-83) in these variants, enabling the activity and thermostability measurements in chelating conditions. Each of the variants includes mutations known previously to increase the thermostability of calcium-independent subtilisin in addition to new stabilizing mutations. S63 has eight amino acid replacements: D41A, M50F, A73L, Q206W, Y217K, N218S, S221C, and Q271E. S63 has 75-fold greater stability than wild type subtilisin in chelating conditions (10 mm EDTA). The other variant, S88, has ten site-specific changes: Q2K, S3C, P5S, K43N, M50F, A73L, Q206C, Y217K, N218S, and Q271E. The two new cysteines form a disulfide bond, and S88 has 1000 times greater stability than wild type subtilisin in chelating conditions. Comparisons of the two new crystal structures (S63 in space group P2(1) with A cell constants 41.2, 78.1, 36.7, and beta = 114.6 degrees and S88 in space group P2(1)2(1)2(1) with cell constants 54.2, 60.4, and 82.7) with previous structures of subtilisin BPN' reveal that the principal changes are in the N-terminal region. The structural bases of the stabilization effects of the new mutations Q2K, S3C, P5S, D41A, Q206C, and Q206W are generally apparent. The effects are attributed to the new disulfide cross-link and to improved hydrophobic packing, new hydrogen bonds, and other rearrangements in the N-terminal region.     

Acetylcholine receptor gating at extracellular transmembrane domain interface: the "pre-M1" linker

Charged residues in the beta10-M1 linker region ("pre-M1") are important in the expression and function of neuromuscular acetylcholine receptors (AChRs). The perturbation of a salt bridge between pre-M1 residue R209 and loop 2 residue E45 has been proposed as being a principle event in the AChR gating conformational "wave." We examined the effects of mutations to all five residues in pre-M1 (positions M207-P211) plus E45 in loop 2 in the mouse alpha(1)-subunit. M207, Q208, and P211 mutants caused small (approximately threefold) changes in the gating equilibrium constant (K(eq)), but the changes for R209, L210, and E45 were larger. Of 19 different side chain substitutions at R209 on the wild-type background, only Q, K, and H generated functional channels, with the largest change in K(eq) (67-fold) from R209Q. Various R209 mutants were functional on different E45 backgrounds: H, Q, and K (E45A), H, A, N, and Q (E45R), and K, A, and N (E45L). Phi values for R209 (on the E45A background), L210, and E45 were 0.74, 0.35, and 0.80, respectively. Phi values for R209 on the wt and three other backgrounds could not be estimated because of scatter. The average coupling energy between 209/45 side chains (six different pairs) was only -0.33 kcal/mol (for both alpha subunits, combined). Pre-M1 residues are important for expression of functional channels and participate in gating, but the relatively modest changes in closed- vs. open-state energy caused mutations, the weak coupling energy between these residues and the functional activity of several unmatched-charge pairs are not consistent with the perturbation of a salt bridge between R209 and E45 playing the principle role in gating.     

Structural basis for type I and type II deficiencies of antithrombotic plasma protein C: Patterns revealed by three-dimensional molecular modelling of mutations of the protease domain

Familial deficiency of protein C is associated with inherited thrombophilia. To explore how specific missense mutations might cause observed clinical phenotypes, known protein C missense mutations were mapped onto three-dimensional homology models of the protein C protease domain, and the implications for domain folding and structure were evaluated. Most Type I missense mutations either replaced internal hydrophobic residues (I201T, L223F, A259V, A267T, A346T, A346V, G376D) or nearby interacting residues (I403M, T298M, Q184H), thus disrupting the packing of internal hydrophobic side chains, or changed hydrophilic residues, thus disrupting ion pairs (N256D, R178W). Mutations (P168L, R169W) at the activation site destabilized the region containing the activation peptide structure. Most Type II mutations involved solvent-exposed residues and were clustered either in a positively charged region (R147W, R157Q, R229Q, R352W) or were located in or near the active site region (S252N, D359N, G381S, G391S, H211Q). The cluster of arginines 147, 157, 229, and 352 may identify a functionally important exosite. Identification of the spatial relationships of natural mutations in the protein C model is helpful for understanding manifestations of protein C deficiency and for identification of novel, functionally important molecular features and exosites. © 1994 John Wiley & Sons, Inc.     

Analysis of <Emphasis Type="Italic">CYP21A2</Emphasis> gene mutations in patients from Ukraine with congenital adrenal hyperplasia

In the article, the data on the distribution of CYP21A2 gene mutations (gene deletion/conversion, c.290-13C>A/G, E110Vfs, I172N, cluster of mutations I236N, V237E, M239K, V281L, Q318X, and R356W) among Ukrainian patients with CAH (congenital adrenal hyperplasia) of different clinical phenotypes are presented. The most common mutation in the studied group (n = 27) is the CYP21A2 gene deletion/conversion. Possible patterns of the studied mutations distribution in different populations of the world and the patients’ genotype–phenotype association are discussed.     

Computational Analysis of Functional Single Nucleotide Polymorphisms Associated with the CYP11B2 Gene

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variations in humans and play a major role in the genetics of human phenotype variation and the genetic basis of human complex diseases. Recently, there is considerable interest in understanding the possible role of the CYP11B2 gene with corticosterone methyl oxidase deficiency, primary aldosteronism, and cardio-cerebro-vascular diseases. Hence, the elucidation of the function and molecular dynamic behavior of CYP11B2 mutations is crucial in current genomics. In this study, we investigated the pathogenic effect of 51 nsSNPs and 26 UTR SNPs in the CYP11B2 gene through computational platforms. Using a combination of SIFT, PolyPhen, I-Mutant Suite, and ConSurf server, four nsSNPs (F487V, V129M, T498A, and V403E) were identified to potentially affect the structure, function, and activity of the CYP11B2 protein. Furthermore, molecular dynamics simulation and structure analyses also confirmed the impact of these nsSNPs on the stability and secondary properties of the CYP11B2 protein. Additionally, utilizing the UTRscan, MirSNP, PolymiRTS and miRNASNP, three SNPs in the 3′UTR region were predicted to exhibit a pattern change in the upstream open reading frames (uORF), and eight microRNA binding sites were found to be highly affected due to 3′UTR SNPs. This cataloguing of deleterious SNPs is essential for narrowing down the number of CYP11B2 mutations to be screened in genetic association studies and for a better understanding of the functional and structural aspects of the CYP11B2 protein.     

Neurovirulence in mice of H5N1 influenza virus genotypes isolated from Hong Kong poultry in 2001

We studied the pathogenicity of five different genotypes (A to E) of highly pathogenic avian H5N1 viruses, which contained HA genes similar to those of the H5N1 virus A/goose/Guangdong/1/96 and five different combinations of "internal" genes, in a mouse model. Highly pathogenic, neurotropic variants of genotypes A, C, D, and E were isolated from the brain after a single intranasal passage in mice. Genotype B virus was isolated from lungs only. The mouse brain variants had amino acid changes in all gene products except PB1, NP, and NS1 proteins but no common sets of mutations. We conclude that the original H5N1/01 isolates of genotypes A, C, D, and E were heterogeneous and that highly pathogenic neurotropic variants can be rapidly selected in mice.