南宁市HIV-1/AIDS人群治疗前pol区遗传特性及蛋白结构分析

为探讨南宁市某县艾滋病病毒1型(HIV-1)感染人群中治疗前pol区遗传特性及蛋白结构变化情况,本研究通过RT-PCR扩增pol区部分序列并进行测序,将序列同源比对构建系统进化树;分型确定毒株亚型和斯坦福大学HIV耐药性数据库比对,分析耐药相关位点;SWISS-MODEL蛋白质同源数据库进行建模分析氨基酸的突变对蛋白质结构和功能的影响。本研究在90份HIV-1标本中获得46个pol区有效序列,共发现4种亚型,其中CRF01_AE占76.08%(35/46)、CRF08_BC占15.22%(7/46)、CRF07_BC占(3/46)6.52%、CRF59_01B1占2.17%(1/46);46个序列中有4例(8.69%)出现耐药突变位点,没有针对核苷酸反转录酶抑制剂(NRTI)的耐药突变;针对蛋白酶类抑制剂(PIs)1例,PR蛋白酶的柔性部位I47V位点发生突变,β折叠结构的I84V位点发生突变,都是异亮氨酸突变为缬氨酸;针对非核苷酸反转录酶抑制剂(NNRTI)有3例,2例位于活性中心的Y181C位点由酪氨酸突变为半胱氨酸,1例位于转角处的E138G位点由谷氨酸突变为甘氨酸。研究表明,南宁市某县HIV-1病毒CRF01_AE重组亚型比例最大,未经抗病毒治疗HIV1感染者中已经出现pol区耐药突变株,突变位点主要位于活性中心及柔性部位,传播水平已经处于中等流行状态。深入分析蛋白质与抑制剂相互作用机制,有助于为艾滋病抗病毒及耐药性监测方案提供科学依据。     

Systematic variation of amino acid substitutions for stringent assessment of pairwise covariation

During protein evolution, amino acids change due to a combination of functional constraints and genetic drift. Proteins frequently contain pairs of amino acids that appear to change together (covariation). Analysis of covariation from naturally occurring sets of orthologs cannot distinguish between residue pairs retained by functional requirements of the protein and those pairs existing due to changes along a common evolutionary path. Here, we have separated the two types of covariation by independently recombining every naturally occurring amino acid variant within a set of 15 subtilisin orthologs. Our analysis shows that in this family of subtilisin orthologs, almost all possible pairwise combinations of amino acids can coexist. This suggests that amino acid covariation found in the subtilisin orthologs is almost entirely due to common ancestral origin of the changes rather than functional constraints. We conclude that naturally occurring sequence diversity can be used to identify positions that can vary independently without destroying protein function.     

Adaptive evolution of HIV at HLA epitopes is associated with ethnicity in Canada

Host immune selection pressure influences the development of mutations that allow for HIV escape. Mutation patterns induced in HIV by the human leukocyte antigen (HLA) are HLA-allele specific. As ethnic groups have distinct and characteristic HLA allele frequencies, we can expect divergent viral evolution within ethnicities. Here, we have sequenced and analyzed the HIV pol gene from 1248 subtype B infected, treatment-na?ve individuals in Canada. Phylogenetic analysis showed no separation between pol sequences from five self-identified ethnic groups, yet fixation index (F(ST)) values showed significant divergence between ethnicities. A total of 17 amino acid sites showed an ethnic-specific fixation pattern (0.015相似文献   

A novel method for detecting intramolecular coevolution: adding a further dimension to selective constraints analyses

Protein evolution depends on intramolecular coevolutionary networks whose complexity is proportional to the underlying functional and structural interactions among sites. Here we present a novel approach that vastly improves the sensitivity of previous methods for detecting coevolution through a weighted comparison of divergence between amino acid sites. The analysis of the HIV-1 Gag protein detected convergent adaptive coevolutionary events responsible for the selective variability emerging between subtypes. Coevolution analysis and functional data for heat-shock proteins, Hsp90 and GroEL, highlight that almost all detected coevolving sites are functionally or structurally important. The results support previous suggestions pinpointing the complex interdomain functional interactions within these proteins and we propose new amino acid sites as important for interdomain functional communication. Three-dimensional information sheds light on the functional and structural constraints governing the coevolution between sites. Our covariation analyses propose two types of coevolving sites in agreement with previous reports: pairs of sites spatially proximal, where compensatory mutations could maintain the local structure stability, and clusters of distant sites located in functional domains, suggesting a functional dependency between them. All sites detected under adaptive evolution in these proteins belong to coevolution groups, further underlining the importance of testing for coevolution in selective constraints analyses.     

Detection of Pairwise Residue Proximity by Covariation Analysis for 3D-Structure Prediction of G-Protein-Coupled Receptors

G-protein-coupled receptor (GPCR) is one of the most important targets for medicines. Homology modeling based on the crystal structure of bovine rhodopsin is currently the most frequently used method for GPCR targeted drug design. Information about residue-residue contacts and the structural specificity in the subfamily is essential for constructing more precise 3D structures, to distinguish the structural differences between the template and targets. In this study, we adopted the covariation analysis to extract information about residue-residue interactions from the amino acid sequence. In the opsin family, a large number of adjacent covarying residue pairs were detected. The detected residue pairs have a strong tendency to gather in some regions important for the structure and function. These results suggest that the covariation analysis is practically utilized to detect adjacent residue pairs and also to apply for predicting functional sites. Analyses of other GPCR subfamilies, olfactory receptor and chemokine receptor families, demonstrated that some adjacent covarying residue pairs were common. Thus, the covariation analysis has possibilities in the substantial improvement of the 3D-structure modeling of GPCRs and in the detection of functional sites such as the ligand-binding sites.     

Isolation and characterization of new proliferating cell nuclear antigen (POL30) mutator mutants that are defective in DNA mismatch repair

A number of studies have suggested a role for proliferating cell nuclear antigen (PCNA) in DNA mismatch repair (MMR). However, the majority of mutations in the POL30 gene encoding PCNA that cause MMR defects also cause replication and other repair defects that contribute to the increased mutation rate caused by these mutations. Here, 20 new pol30 mutants were identified and screened for MMR and other defects, resulting in the identification of two mutations, pol30-201 and pol30-204, that appear to cause MMR defects but little if any other defects. The pol30-204 mutation altered an amino acid (C81R) in the monomer-monomer interface region and resulted in a partial general MMR defect and a defect in MSH2-MSH6 binding in vitro. The pol30-201 mutation altered an amino acid (C22Y) located on the surface of the PCNA trimer that slides over the DNA but did not cause a defect in MSH2-MSH6 binding in vitro. The pol30-201 mutation caused an intermediate mutator phenotype. However, the pol30-201 mutation caused almost a complete defect in the repair of AC and GT mispairs and only a small defect in the repair of a "+T" insertion, an effect similar to that caused by an msh6Delta mutation, indicating that pol30-201 primarily effects MSH6-dependent MMR. The chromosomal double mutant msh3-FF>AA msh6-FF>AA eliminating the conserved FF residues of the PCNA interacting motif of these proteins caused a small (<10%) defect in MMR but showed synergistic interactions with mutations in POL30, indicating that the FF>AA substitution may not eliminate PCNA interactions in vivo. These results indicate that the interaction between PCNA and MMR proteins is more complex than was previously appreciated.     

Structure-function defects of human mitochondrial DNA polymerase in autosomal dominant progressive external ophthalmoplegia

Progressive external ophthalmoplegia (PEO) is a mitochondrial disorder associated with mutations in the POLG gene encoding the mitochondrial DNA polymerase (pol gamma). Four autosomal dominant mutations that cause PEO encode the amino acid substitutions G923D, R943H, Y955C and A957S in the polymerase domain of pol gamma. A homology model of the pol gamma catalytic domain in complex with DNA was developed to investigate the effects of these mutations. Two mutations causing the most severe disease phenotype, Y955C and R943H, change residues that directly interact with the incoming dNTP. Polymerase mutants exhibit 0.03-30% wild-type polymerase activity and a 2- to 35-fold decrease in nucleotide selectivity in vitro. The reduced selectivity and catalytic efficiency of the autosomal dominant PEO mutants predict in vivo dysfunction, and the extent of biochemical defects correlates with the clinical severity of the disease.     

Identification of HCV protease inhibitor resistance mutations by selection pressure-based method

A major challenge to successful antiviral therapy is the emergence of drug-resistant viruses. Recent studies have developed several automated analyses of HIV sequence polymorphism based on calculations of selection pressure (K_a/K_s) to predict drug resistance mutations. Similar resistance analysis programs for HCV inhibitors are not currently available. Taking advantage of the recently available sequence data of patient HCV samples from a Phase II clinical study of protease inhibitor boceprevir, we calculated the selection pressure for all codons in the HCV protease region (amino acid 1–181) to identify potential resistance mutations. The correlation between mutations was also calculated to evaluate linkage between any two mutations. Using this approach, we identified previously known major resistant mutations, including a recently reported mutation V55A. In addition, a novel mutation V158I was identified, and we further confirmed its resistance to boceprevir in protease enzyme and replicon assay. We also extended the approach to analyze potential interactions between individual mutations and identified three pairs of correlated changes. Our data suggests that selection pressure-based analysis and correlation mapping could provide useful tools to analyze large amount of sequencing data from clinical samples and to identify new drug resistance mutations as well as their linkage and correlations.     

Codon-substitution models for heterogeneous selection pressure at amino acid sites

Comparison of relative fixation rates of synonymous (silent) and nonsynonymous (amino acid-altering) mutations provides a means for understanding the mechanisms of molecular sequence evolution. The nonsynonymous/synonymous rate ratio (omega = d(N)d(S)) is an important indicator of selective pressure at the protein level, with omega = 1 meaning neutral mutations, omega < 1 purifying selection, and omega > 1 diversifying positive selection. Amino acid sites in a protein are expected to be under different selective pressures and have different underlying omega ratios. We develop models that account for heterogeneous omega ratios among amino acid sites and apply them to phylogenetic analyses of protein-coding DNA sequences. These models are useful for testing for adaptive molecular evolution and identifying amino acid sites under diversifying selection. Ten data sets of genes from nuclear, mitochondrial, and viral genomes are analyzed to estimate the distributions of omega among sites. In all data sets analyzed, the selective pressure indicated by the omega ratio is found to be highly heterogeneous among sites. Previously unsuspected Darwinian selection is detected in several genes in which the average omega ratio across sites is <1, but in which some sites are clearly under diversifying selection with omega > 1. Genes undergoing positive selection include the beta-globin gene from vertebrates, mitochondrial protein-coding genes from hominoids, the hemagglutinin (HA) gene from human influenza virus A, and HIV-1 env, vif, and pol genes. Tests for the presence of positively selected sites and their subsequent identification appear quite robust to the specific distributional form assumed for omega and can be achieved using any of several models we implement. However, we encountered difficulties in estimating the precise distribution of omega among sites from real data sets.     

Deep Sequencing of Protease Inhibitor Resistant HIV Patient Isolates Reveals Patterns of Correlated Mutations in Gag and Protease

While the role of drug resistance mutations in HIV protease has been studied comprehensively, mutations in its substrate, Gag, have not been extensively cataloged. Using deep sequencing, we analyzed a unique collection of longitudinal viral samples from 93 patients who have been treated with therapies containing protease inhibitors (PIs). Due to the high sequence coverage within each sample, the frequencies of mutations at individual positions were calculated with high precision. We used this information to characterize the variability in the Gag polyprotein and its effects on PI-therapy outcomes. To examine covariation of mutations between two different sites using deep sequencing data, we developed an approach to estimate the tight bounds on the two-site bivariate probabilities in each viral sample, and the mutual information between pairs of positions based on all the bounds. Utilizing the new methodology we found that mutations in the matrix and p6 proteins contribute to continued therapy failure and have a major role in the network of strongly correlated mutations in the Gag polyprotein, as well as between Gag and protease. Although covariation is not direct evidence of structural propensities, we found the strongest correlations between residues on capsid and matrix of the same Gag protein were often due to structural proximity. This suggests that some of the strongest inter-protein Gag correlations are the result of structural proximity. Moreover, the strong covariation between residues in matrix and capsid at the N-terminus with p1 and p6 at the C-terminus is consistent with residue-residue contacts between these proteins at some point in the viral life cycle.     

Persistence of transmitted drug resistance among subjects with primary human immunodeficiency virus infection

Following interruption of antiretroviral therapy among individuals with acquired drug resistance, preexisting drug-sensitive virus emerges relatively rapidly. In contrast, wild-type virus is not archived in individuals infected with drug-resistant human immunodeficiency virus (HIV) and thus cannot emerge rapidly in the absence of selective drug pressure. Fourteen recently HIV-infected patients with transmitted drug-resistant virus were followed for a median of 2.1 years after the estimated date of infection (EDI) without receiving antiretroviral therapy. HIV drug resistance and pol replication capacity (RC) in longitudinal plasma samples were assayed. Resistance mutations were characterized as pure populations or mixtures. The mean time to first detection of a mixture of wild-type and drug-resistant viruses was 96 weeks (1.8 years) (95% confidence interval, 48 to 192 weeks) after the EDI. The median time to loss of detectable drug resistance using population-based assays ranged from 4.1 years (conservative estimate) to longer than the lifetime of the individual (less conservative estimate). The transmission of drug-resistant virus was not associated with virus with reduced RC. Sexual transmission of HIV selects for highly fit drug-resistant variants that persist for years. The prolonged persistence of transmitted drug resistance strongly supports the routine use of HIV resistance genotyping for all newly diagnosed individuals.     

Mutation patterns and structural correlates in human immunodeficiency virus type 1 protease following different protease inhibitor treatments

Although many human immunodeficiency virus type 1 (HIV-1)-infected persons are treated with multiple protease inhibitors in combination or in succession, mutation patterns of protease isolates from these persons have not been characterized. We collected and analyzed 2,244 subtype B HIV-1 isolates from 1,919 persons with different protease inhibitor experiences: 1,004 isolates from untreated persons, 637 isolates from persons who received one protease inhibitor, and 603 isolates from persons receiving two or more protease inhibitors. The median number of protease mutations per isolate increased from 4 in untreated persons to 12 in persons who had received four or more protease inhibitors. Mutations at 45 of the 99 amino acid positions in the protease-including 22 not previously associated with drug resistance-were significantly associated with protease inhibitor treatment. Mutations at 17 of the remaining 99 positions were polymorphic but not associated with drug treatment. Pairs and clusters of correlated (covarying) mutations were significantly more likely to occur in treated than in untreated persons: 115 versus 23 pairs and 30 versus 2 clusters, respectively. Of the 115 statistically significant pairs of covarying residues in the treated isolates, 59 were within 8 A of each other-many more than would be expected by chance. In summary, nearly one-half of HIV-1 protease positions are under selective drug pressure, including many residues not previously associated with drug resistance. Structural factors appear to be responsible for the high frequency of covariation among many of the protease residues. The presence of mutational clusters provides insight into the complex mutational patterns required for HIV-1 protease inhibitor resistance.     

Structural constraints identified with covariation analysis in ribosomal RNA

Covariation analysis is used to identify those positions with similar patterns of sequence variation in an alignment of RNA sequences. These constraints on the evolution of two positions are usually associated with a base pair in a helix. While mutual information (MI) has been used to accurately predict an RNA secondary structure and a few of its tertiary interactions, early studies revealed that phylogenetic event counting methods are more sensitive and provide extra confidence in the prediction of base pairs. We developed a novel and powerful phylogenetic events counting method (PEC) for quantifying positional covariation with the Gutell lab's new RNA Comparative Analysis Database (rCAD). The PEC and MI-based methods each identify unique base pairs, and jointly identify many other base pairs. In total, both methods in combination with an N-best and helix-extension strategy identify the maximal number of base pairs. While covariation methods have effectively and accurately predicted RNAs secondary structure, only a few tertiary structure base pairs have been identified. Analysis presented herein and at the Gutell lab's Comparative RNA Web (CRW) Site reveal that the majority of these latter base pairs do not covary with one another. However, covariation analysis does reveal a weaker although significant covariation between sets of nucleotides that are in proximity in the three-dimensional RNA structure. This reveals that covariation analysis identifies other types of structural constraints beyond the two nucleotides that form a base pair.     

Contributions of residue pairing to beta-sheet formation: conservation and covariation of amino acid residue pairs on antiparallel beta-strands

In an effort to better understand beta-sheet assembly, we have investigated the evolutionary behavior of neighboring residues on adjacent antiparallel beta-strands. Residue pairs were classified according to solvent exposure as well as by whether their backbone NH and C==O groups are hydrogen bonded. The conservation and covariation of 19,241 pairs in 219 sequence alignments was analyzed. Buried pairs were found to be the most conserved, while stronger covariation was detected in the solvent-exposed pairs. However, residues on neighboring strands showed a degree of conservation and covariation similar to that of well-separated residues on the same strand, suggesting that evolutionary pressure to maintain complementarity between pairs on neighboring strands is weak. Moreover, in spite of the preference of certain amino acid pairs to occupy neighboring positions on adjacent strands, such favored pairs are neither more strongly mutually conserved nor covary more strongly than pairs of the same type in non-interacting positions. Although the beta-sheet pairs did not show outstanding evolutionary coupling, in many protein families significant conservation and covariation patterns were detected for some of the residue pairs. Overall, the weak evolutionary conservation and covariation of the beta-sheet pairs indicates that sheet structure is unlikely to be dictated by specific side-chain interactions.     

Adaptive covariation between the coat and movement proteins of prunus necrotic ringspot virus

The relative functional and/or structural importance of different amino acid sites in a protein can be assessed by evaluating the selective constraints to which they have been subjected during the course of evolution. Here we explore such constraints at the linear and three-dimensional levels for the movement protein (MP) and coat protein (CP) encoded by RNA 3 of prunus necrotic ringspot ilarvirus (PNRSV). By a maximum-parsimony approach, the nucleotide sequences from 46 isolates of PNRSV varying in symptomatology, host tree, and geographic origin have been analyzed and sites under different selective pressures have been identified in both proteins. We have also performed covariation analyses to explore whether changes in certain amino acid sites condition subsequent variation in other sites of the same protein or the other protein. These covariation analyses shed light on which particular amino acids should be involved in the physical and functional interaction between MP and CP. Finally, we discuss these findings in the light of what is already known about the implication of certain sites and domains in structure and protein-protein and RNA-protein interactions.     

Unique tertiary and neighbor interactions determine conservation patterns of Cis Watson-Crick A/G base-pairs

X-ray, phylogenetic and quantum chemical analysis of molecular interactions and conservation patterns of cis Watson-Crick (W.C.) A/G base-pairs in 16S rRNA, 23S rRNA and other molecules was carried out. In these base-pairs, the A and G nucleotides interact with their W.C. edges with glycosidic bonds oriented cis relative to each other. The base-pair is stabilised by two hydrogen bonds, the C1'-C1' distance is enlarged and the G(N2) amino group is left unpaired. Quantum chemical calculations show that, in the absence of other interactions, the unpaired amino group is substantially non-planar due to its partial sp(3) pyramidalization, while the whole base-pair is internally propeller twisted and very flexible. The unique molecular properties of the cis W.C. A/G base-pairs make them distinct from other base-pairs. They occur mostly at the ends of canonical helices, where they serve as interfaces between the helix and other motifs. The cis W.C. A/G base-pairs play crucial roles in natural RNA structures with salient sequence conservation patterns. The key contribution to conservation is provided by the unpaired G(N2) amino group that is involved in a wide range of tertiary and neighbor contacts in the crystal structures. Many of them are oriented out of the plane of the guanine base and utilize the partial sp(3) pyramidalization of the G(N2). There is a lack of A/G to G/A covariation, which, except for the G(N2) position, would be entirely isosteric. On the contrary, there is a rather frequent occurrence of G/A to G/U covariation, as the G/U wobble base-pair has an unpaired amino group in the same position as the cis W.C. G/A base-pair. The cis W.C. A/G base-pairs are not conserved when there is no tertiary or neighbor interaction. Obtaining the proper picture of the interactions and phylogenetic patterns of the cis W.C. A/G base-pairs requires a detailed analysis of the relation between the molecular structures and the energetics of interactions at a level of single H-bonds and contacts.