Relative Binding Affinity Prediction of Charge-Changing Sequence Mutations with FEP in Protein–Protein Interfaces

Building on the substantial progress that has been made in using free energy perturbation (FEP) methods to predict the relative binding affinities of small molecule ligands to proteins, we have previously shown that results of similar quality can be obtained in predicting the effect of mutations on the binding affinity of protein–protein complexes. However, these results were restricted to mutations which did not change the net charge of the side chains due to known difficulties with modeling perturbations involving a change in charge in FEP. Various methods have been proposed to address this problem. Here we apply the co-alchemical water approach to study the efficacy of FEP calculations of charge changing mutations at the protein–protein interface for the antibody–gp120 system investigated previously and three additional complexes. We achieve an overall root mean square error of 1.2 kcal/mol on a set of 106 cases involving a change in net charge selected by a simple suitability filter using side-chain predictions and solvent accessible surface area to be relevant to a biologic optimization project. Reasonable, although less precise, results are also obtained for the 44 more challenging mutations that involve buried residues, which may in some cases require substantial reorganization of the local protein structure, which can extend beyond the scope of a typical FEP simulation. We believe that the proposed prediction protocol will be of sufficient efficiency and accuracy to guide protein engineering projects for which optimization and/or maintenance of a high degree of binding affinity is a key objective.     

Isolation and Characterization of Mutations in the β-Tubulin Gene of SACCHAROMYCES CEREVISIAE

Of 173 mutants of Saccharomyces cerevisiae resistant to the antimitotic drug benomyl (BenR), six also conferred cold-sensitivity for growth and three others conferred temperature-sensitivity for growth in the absence of benomyl. All of the benR mutations tested, including the nine conditional-lethal mutations, were shown to be in the same gene. This gene, TUB2, has previously been molecularly cloned and identified as the yeast structural gene encoding beta-tubulin. Four of the conditional-lethal alleles of TUB2 were mapped to particular restriction fragments within the gene. One of these mutations was cloned and sequenced, revealing a single amino acid change, from arginine to histidine at amino acid position 241, which is responsible for both the BenR and the cold-sensitive lethal phenotypes. The terminal arrest morphology of conditional-lethal alleles of TUB2 at their restrictive temperature showed a characteristic cell-division-cycle defect, suggesting a requirement for tubulin function primarily in mitosis during the vegetative growth cycle. The TUB2 gene was genetically mapped to the distal left arm of chromosome VI, very near the actin gene, ACT1; no CDC (cell-division-cycle) loci have been mapped previously to this location. TUB2 is thus the first cell-division-cycle gene known to encode a cytoskeletal protein that has been identified in S. cerevisiae.     

Protein–Protein and Protein–Membrane Associations in the Lignin Pathway

Supramolecular organization of enzymes is proposed to orchestrate metabolic complexity and help channel intermediates in different pathways. Phenylpropanoid metabolism has to direct up to 30% of the carbon fixed by plants to the biosynthesis of lignin precursors. Effective coupling of the enzymes in the pathway thus seems to be required. Subcellular localization, mobility, protein–protein, and protein–membrane interactions of four consecutive enzymes around the main branch point leading to lignin precursors was investigated in leaf tissues of Nicotiana benthamiana and cells of Arabidopsis thaliana. CYP73A5 and CYP98A3, the two Arabidopsis cytochrome P450s (P450s) catalyzing para- and meta-hydroxylations of the phenolic ring of monolignols were found to colocalize in the endoplasmic reticulum (ER) and to form homo- and heteromers. They moved along with the fast remodeling plant ER, but their lateral diffusion on the ER surface was restricted, likely due to association with other ER proteins. The connecting soluble enzyme hydroxycinnamoyltransferase (HCT), was found partially associated with the ER. Both HCT and the 4-coumaroyl-CoA ligase relocalized closer to the membrane upon P450 expression. Fluorescence lifetime imaging microscopy supports P450 colocalization and interaction with the soluble proteins, enhanced by the expression of the partner proteins. Protein relocalization was further enhanced in tissues undergoing wound repair. CYP98A3 was the most effective in driving protein association.     

Protein–protein interactions of huntingtin in the hippocampus

Molecular Biology - Huntingtin (HTT) occurs in the neuronal cytoplasm and can interact with structural elements of synapses. Huntington’s disease (HD) results from pathological expansion of a...     

Investigation of Gene Insertion–Deletion Polymorphism of the Gene for Angiotensin-Converting Enzyme in Populations of the Volga–Ural Region

Insertion–deletion polymorphism at the angiotensin I-converting enzyme (ACE) gene in populations of the Volga–Ural region was examined by means of polymerase chain reaction. The populations studied belong to the Finno-Ugric (Komis, Maris, Mordovians, and Udmurts), Turkic (Chuvashes, Tatars, and Bashkirs), and Eastern-Slavic (Russians) ethnic groups. Distribution patterns of allele and genotype frequencies of this polymorphic system in the examined region were characterized. Comparison of the obtained results with the literature data on the ACE gene polymorphism in other Caucasoid and Mongoloid populations revealed some trends in the ACE genotype frequency dynamics depending on the ethnicity of the populations.     

A Mutation in the LDL Receptor–Related Protein 5 Gene Results in the Autosomal Dominant High–Bone-Mass Trait

Osteoporosis is a complex disease that affects >10 million people in the United States and results in 1.5 million fractures annually. In addition, the high prevalence of osteopenia (low bone mass) in the general population places a large number of people at risk for developing the disease. In an effort to identify genetic factors influencing bone density, we characterized a family that includes individuals who possess exceptionally dense bones but are otherwise phenotypically normal. This high–bone-mass trait (HBM) was originally localized by linkage analysis to chromosome 11q12-13. We refined the interval by extending the pedigree and genotyping additional markers. A systematic search for mutations that segregated with the HBM phenotype uncovered an amino acid change, in a predicted β-propeller module of the low-density lipoprotein receptor–related protein 5 (LRP5), that results in the HBM phenotype. During analysis of >1,000 individuals, this mutation was observed only in affected individuals from the HBM kindred. By use of in situ hybridization to rat tibia, expression of LRP5 was detected in areas of bone involved in remodeling. Our findings suggest that the HBM mutation confers a unique osteogenic activity in bone remodeling, and this understanding may facilitate the development of novel therapies for the treatment of osteoporosis.     

Multiple Protein Tyrosine Phosphatases in Sponges and Explosive Gene Duplication in the Early Evolution of Animals Before the Parazoan–Eumetazoan Split

Protein tyrosine phosphatases (PTPs) regulate various physiological events in animal cells. They comprise a diverse family which are classified into two categories, receptor type and nonreceptor type. From the domain organization and phylogenetic tree, we have classified known PTPs into 17 subtypes (9 receptor-type and 8 nonreceptor-type PTPs) which are characterized by different organization of functional domain and independent cluster in tree. The receptor type PTPs are thought to be implicated in cell–cell adhesion by association of cell adhesion molecules. Since sponges are the most primitive multicellular animals and are thought to be lacking cell cohesiveness and coordination typical of eumetazoans, cloning and sequencing of PTP cDNAs of Ephydatia fluviatilis (freshwater sponge) have been conducted by RT-PCR to determine whether or not sponges have PTP genes in their genomes. We have isolated nine PTPs, of which five are possibly receptor type. A phylogenetic tree including the sponge PTPs revealed that most of the gene duplications that gave rise to the 17 subtypes had been completed in the very early evolution of animals before the parazoan–eumetazoan split, the earliest branching among extant animal phyla. The family tree also revealed the rapid evolutionary rate of PTP subtypes in the early stage of animal evolution. Received: 22 October 1998 / Accepted: 27 November 1998     

Residues 240–250 in the C-Terminus of the Pirh2 Protein Complement the Function of the RING Domain in Self-Ubiquitination of the Pirh2 Protein

Pirh2 is a p53 inducible gene that encodes a RING-H2 domain and is proposed to be a main regulator of p53 protein, thus fine tuning the DNA damage response. Pirh2 interacts physically with p53 and promotes its MDM2-independent ubiquitination and subsequent degradation as well as participates in an auto-regulatory feedback loop that controls p53 function. Pirh2 also self-ubiquitinates. Interestingly, Pirh2 is overexpressed in a wide range of human tumors. In this study, we investigated the domains and residues essential for Pirh2 self-ubiquitination. Deletions were made in each of the three major domains of Pirh2: the N-terminal domain (NTD), Ring domain (RING), and C-terminal domain (CTD). The effects of these deletions on Pirh2 self-ubiquitination were then assessed using in vitro ubiquitination assays. Our results demonstrate that the RING domain is essential, but not sufficient, for Pirh2 self-ubiquitination and that residues 240–250 of the C-terminal domain are also essential. Our results demonstrate that Pirh2 mediated p53 polyubiquitination occurs mainly through the K48 residue of ubiquitin in vitro. Our data further our understanding of the mechanism of Pirh2 self-ubiquitination and may help identify valuable therapeutic targets that play roles in reducing the effects of the overexpression of Pirh2, thus maximizing p53''s response to DNA damage.     

Spectrum of Mutations in the ATP7B Gene in Russian Patients with Wilson’s Disease

Russian Journal of Genetics - Wilson’s disease (WD) is an autosomal recessive disease caused by an excessive accumulation of copper. The molecular genetic etiology of the disease is due to...     

Crystal Structure of the Bacteriophage Qβ Coat Protein in Complex with the RNA Operator of the Replicase Gene

The coat proteins of single-stranded RNA bacteriophages specifically recognize and bind to a hairpin structure in their genome at the beginning of the replicase gene. The interaction serves to repress the synthesis of the replicase enzyme late in infection and contributes to the specific encapsidation of phage RNA. While this mechanism is conserved throughout the Leviviridae family, the coat protein and operator sequences from different phages show remarkable variation, serving as prime examples for the co-evolution of protein and RNA structure. To better understand the protein–RNA interactions in this virus family, we have determined the three-dimensional structure of the coat protein from bacteriophage Qβ bound to its cognate translational operator. The RNA binding mode of Qβ coat protein shares several features with that of the widely studied phage MS2, but only one nucleotide base in the hairpin loop makes sequence-specific contacts with the protein. Unlike in other RNA phages, the Qβ coat protein does not utilize an adenine-recognition pocket for binding a bulged adenine base in the hairpin stem but instead uses a stacking interaction with a tyrosine side chain to accommodate the base. The extended loop between β strands E and F of Qβ coat protein makes contacts with the lower part of the RNA stem, explaining the greater length dependence of the RNA helix for optimal binding to the protein. Consequently, the complex structure allows the proposal of a mechanism by which the Qβ coat protein recognizes and discriminates in favor of its cognate RNA.     

Polymorphism of Human Glutatione S-Transferase Gene in the Populations of the Volga–Ural Region

The frequency of the GstM1 gene deletion homozygotes in eight populations of the Volga–Ural region belonging according to linguistic classification to Turkic (Bashkirs, Tatars, and Chuvashs), Finno–Ugric (Maris, Komis, Mordovians, and Udmurts), and Eastern–Slavic (Russians) ethnic groups, was examined by means of PCR technique. The frequency of the deletion homozygotes varied from 41.4% in Bashkirs to 61.3% in Mordovians. The mean deletion frequency comprised 50.1%, which was consistent with the data for European populations (²= 0.009).     

Spectroscopic Characterization of the Excitation Energy Transfer in the Fucoxanthin–Chlorophyll Protein of Diatoms

We characterized the energy transfer pathways in the fucoxanthin–chlorophyll protein (FCP) complex of the diatom Cyclotella meneghiniana by conducting ultrafast transient absorption measurements. This light harvesting antenna has a distinct pigment composition and binds chlorophyll a (Chl-a), fucoxanthin and chlorophyll c (Chl-c) molecules in a 4:4:1 ratio. We find that upon excitation of fucoxanthin to its S₂ state, a significant amount of excitation energy is transferred rapidly to Chl-a. The ensuing dynamics illustrate the presence of a complex energy transfer network that also involves energy transfer from the unrelaxed or ‘hot’ intermediates. Chl-c to Chl-a energy transfer occurs on a timescale of a 100 fs. We observe no significant spectral evolution in the Chl-a region of the spectrum. We have applied global and target analysis to model the measured excited state dynamics and estimate the spectra of the states involved; the energy transfer network is discussed in relation to the pigment organization of the FCP complex.     

Mutations in the 5’ NTR and the Non-Structural Protein 3A of the Coxsackievirus B3 Selectively Attenuate Myocarditogenicity

The 5’ non-translated region (NTR) is an important molecular determinant that controls replication and virulence of coxsackievirus B (CVB)3. Previous studies have reported many nucleotide (nt) sequence differences in the Nancy strain of the virus, including changes in the 5’ NTR with varying degrees of disease severity. In our studies of CVB3-induced myocarditis, we sought to generate an infectious clone of the virus for routine in vivo experimentation. By determining the viral nt sequence, we identified three new nt substitutions in the clone that differed from the parental virus strain: C97U in the 5’ NTR; a silent mutation, A4327G, in non-structural protein 2C; and C5088U (resulting in P1449L amino acid change) in non-structural protein 3A of the virus leading us to evaluate the role of these changes in the virulence properties of the virus. We noted that the disease-inducing ability of the infectious clone-derived virus in three mouse strains was restricted to pancreatitis alone, and the incidence and severity of myocarditis were significantly reduced. We then reversed the mutations by creating three new clones, representing 1) U97C; 2) G4327A and U5088C; and 3) their combination together in the third clone. The viral titers obtained from all the clones were comparable, but the virions derived from the third clone induced myocarditis comparable to that induced by wild type virus; however, the pancreatitis-inducing ability remained unaltered, suggesting that the mutations described above selectively influence myocarditogenicity. Because the accumulation of mutations during passages is a continuous process in RNA viruses, it is possible that CVB3 viruses containing such altered nts may evolve naturally, thus favoring their survival in the environment.     

Analysis of Polymorphisms of the Huntington Disease Gene in Ethnic Populations of the Volga–Ural Region

Eleven populations of the Volga–Ural region were analyzed with respect to three intragenic polymorphisms of the Huntington disease gene (IT15), including highly polymorphic (CAG)_n and moderately polymorphic (CCG)_n of exon 1 and neutral del2642 of exon 58. In the case of (CAG)_n, 101 genotypes were observed, with genotype number varying from 15 in Southeastern Bashkirs to 34 in Mari. Allele diversity R_S ranged from 9.70 in Southeastern Bashkirs to 18.00 in Chuvash, averaging 13.79 ± 2.12. The (CAG)_n allele frequency distribution was unimodal and had a maximum at (CAG)₁₇. In the case of (CCG)_n, six alleles with 6–12 repeats were observed. R_S was 4.13 ± 0.44, ranging from 3.73 in Udmurts to 4.99 in Tatars. In the case of del2642, allele del– was detected at a frequency 0.830 in Mari to 0.932 in Udmurts. Of all Volga–Ural ethnic populations, Finno-Ugric ones proved to be most heterogeneous with respect to the three polymorphisms, whereas Turkic populations and, in particular, Bashkirs were homogeneous. Microdifferentiation of the Volga–Ural populations corresponded to the European type.     

Characterization of Protein–Protein Interfaces

We analyze the characteristics of protein–protein interfaces using the largest datasets available from the Protein Data Bank (PDB). We start with a comparison of interfaces with protein cores and non-interface surfaces. The results show that interfaces differ from protein cores and non-interface surfaces in residue composition, sequence entropy, and secondary structure. Since interfaces, protein cores, and non-interface surfaces have different solvent accessibilities, it is important to investigate whether the observed differences are due to the differences in solvent accessibility or differences in functionality. We separate out the effect of solvent accessibility by comparing interfaces with a set of residues having the same solvent accessibility as the interfaces. This strategy reveals residue distribution propensities that are not observable by comparing interfaces with protein cores and non-interface surfaces. Our conclusions are that there are larger numbers of hydrophobic residues, particularly aromatic residues, in interfaces, and the interactions apparently favored in interfaces include the opposite charge pairs and hydrophobic pairs. Surprisingly, Pro-Trp pairs are over represented in interfaces, presumably because of favorable geometries. The analysis is repeated using three datasets having different constraints on sequence similarity and structure quality. Consistent results are obtained across these datasets. We have also investigated separately the characteristics of heteromeric interfaces and homomeric interfaces.     

Utilization of Methyl Proton Resonances in Cross-Saturation Measurement for Determining the Interfaces of Large Protein–Protein Complexes

Cross-saturation experiments allow the identification of the contact residues of large protein complexes (MW>50 K) more rigorously than conventional NMR approaches which involve chemical shift perturbations and hydrogen-deuterium exchange experiments [Takahashi et al. (2000) Nat. Struct. Biol., 7, 220–223]. In the amide proton-based cross-saturation experiment, the combined use of high deuteration levels for non-exchangeable protons of the ligand protein and a solvent with a low concentration of ¹H₂O greatly enhanced the selectivity of the intermolecular cross-saturation phenomenon. Unfortunately, experimental limitations caused losses in sensitivity. Furthermore, since main chain amide protons are not generally exposed to solvent, the efficiency of the saturation transfer directed to the main chain amide protons is not very high. Here we propose an alternative cross-saturation experiment which utilizes the methyl protons of the side chains of the ligand protein. Owing to the fast internal rotation along the methyl axis, we theoretically and experimentally demonstrated the enhanced efficiency of this approach. The methyl-utilizing cross-saturation experiment has clear advantages in sensitivity and saturation transfer efficiency over the amide proton-based approach. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Analysis of Apolipoprotein E Gene Polymorphism in Populations of the Volga–Ural Region

Polymorphism at the apolipoprotein E gene (apoE) in populations of the Volga–Ural region was studied by means of polymerase chain reaction. In the region examined the population-specific patterns of the apoE alleles and genotypes frequency distribution were established. The results obtained were compared with the literature data on the apoE polymorphism in other world populations. Substantial heterogeneity of different ethnic populations in respect to the apoE genotypes distribution and frequency was revealed.