Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Folding simulations of a three-stranded antiparallel β-sheet peptide

Protein folding is a grand challenge of the postgenomic era. In this paper, 58 folding events sampled during 47 molecular dynamics trajectories for a total simulation time of more than 4 μs provide an atomic detail picture of the folding of a 20-residue synthetic peptide with a stable three-stranded antiparallel β-sheet fold. The simulations successfully reproduce the NMR solution conformation, irrespective of the starting structure. The sampling of the conformational space is sufficient to determine the free energy surface and localize the minima and transition states. The statistically predominant folding pathway involves the formation of contacts between strands 2 and 3, starting with the side chains close to the turn, followed by association of the N-terminal strand onto the preformed 2–3 β-hairpin. The folding mechanism presented here, formation of a β-hairpin followed by consolidation, is in agreement with a computational study of the free energy surface of another synthetic three-stranded antiparallel β-sheet by Bursulaya and Brooks [(1999) J. Am. Chem. Soc. 121, 9947–9951]. Hence, it might hold in general for antiparallel β-sheets with short turns.

Plastid-localized acetyl-CoA carboxylase of bread wheat is encoded by a single gene on each of the three ancestral chromosome sets

5′-End fragments of two genes encoding plastid-localized acetyl-CoA carboxylase (ACCase; EC 6.4.1.2) of wheat (Triticum aestivum) were cloned and sequenced. The sequences of the two genes, Acc-1,1 and Acc-1,2, are 89% identical. Their exon sequences are 98% identical. The amino acid sequence of the biotin carboxylase domain encoded by Acc-1,1 and Acc-1,2 is 93% identical with the maize plastid ACCase but only 80–84% identical with the cytosolic ACCases from other plants and from wheat. Four overlapping fragments of cDNA covering the entire coding region were cloned by PCR and sequenced. The wheat plastid ACCase ORF contains 2,311 amino acids with a predicted molecular mass of 255 kDa. A putative transit peptide is present at the N terminus. Comparison of the genomic and cDNA sequences revealed introns at conserved sites found in the genes of other plant multifunctional ACCases, including two introns absent from the wheat cytosolic ACCase genes. Transcription start sites of the plastid ACCase genes were estimated from the longest cDNA clones obtained by 5′-RACE (rapid amplification of cDNA ends). The untranslated leader sequence encoded by the Acc-1 genes is at least 130–170 nucleotides long and is interrupted by an intron. Southern analysis indicates the presence of only one copy of the gene in each ancestral chromosome set. The gene maps near the telomere on the short arm of chromosomes 2A, 2B, and 2D. Identification of three different cDNAs, two corresponding to genes Acc-1,1 and Acc-1,2, indicates that all three genes are transcriptionally active.

Comparing in vitro, in situ, and in vivo experimental data in a three-dimensional model of mammalian cochlear mechanics

Normal mammalian hearing is refined by amplification of the motion of the cochlear partition. This partition, comprising the organ of Corti sandwiched between the basilar and tectorial membranes, contains the outer hair cells that are thought to drive this amplification process. Force generation by outer hair cells has been studied extensively in vitro and in situ, but, to understand cochlear amplification fully, it is necessary to characterize the role played by each of the components of the cochlear partition in vivo. Observations of cochlear partition motion in vivo are severely restricted by its inaccessibility and sensitivity to surgical trauma, so, for the present study, a computer model has been used to simulate the operation of the cochlea under different experimental conditions. In this model, which uniquely retains much of the three-dimensional complexity of the real cochlea, the motions of the basilar and tectorial membranes are fundamentally different during in situ- and in vivo-like conditions. Furthermore, enhanced outer hair cell force generation in vitro leads paradoxically to a decrease in the gain of the cochlear amplifier during sound stimulation to the model in vivo. These results suggest that it is not possible to extrapolate directly from experimental observations made in vitro and in situ to the normal operation of the intact organ in vivo.

Kinematic geometry of mass-triangles and reduction of Schrödinger’s equation of three-body systems to partial differential equations solely defined on triangular parameters

Schrödinger’s equation of a three-body system is a linear partial differential equation (PDE) defined on the 9-dimensional configuration space, ℝ9, naturally equipped with Jacobi’s kinematic metric and with translational and rotational symmetries. The natural invariance of Schrödinger’s equation with respect to the translational symmetry enables us to reduce the configuration space to that of a 6-dimensional one, while that of the rotational symmetry provides the quantum mechanical version of angular momentum conservation. However, the problem of maximizing the use of rotational invariance so as to enable us to reduce Schrödinger’s equation to corresponding PDEs solely defined on triangular parameters—i.e., at the level of ℝ6/SO(3)—has never been adequately treated. This article describes the results on the orbital geometry and the harmonic analysis of (SO(3),ℝ6) which enable us to obtain such a reduction of Schrödinger’s equation of three-body systems to PDEs solely defined on triangular parameters.

Transcription initiation is controlled by three core promoter elements in the hgl5 gene of the protozoan parasite Entamoeba histolytica

Entamoeba histolytica is a single cell eukaryote that is the etiologic agent of amoebic colitis. Core promoter elements of E. histolytica protein encoding genes include a TATA-like sequence (GTATTTAAAG/C) at −30, a novel element designated GAAC (GAACT) that has a variable location between TATA and the site of transcription initiation, and a putative initiator (Inr) element (AAAAATTCA) overlying the site of transcription initiation. The presence of three separate conserved sequences in a eukaryotic core promoter is unprecedented and prompted examination of their roles in regulating transcription initiation. Alterations of all three regions in the hgl5 gene decreased reporter gene activity with the greatest effect seen by mutation of the GAAC element. Positional analysis of the TATA box demonstrated that transcription initiated consistently 30–31 bases downstream of the TATA region. Mutation of either the TATA or GAAC elements resulted in the appearance of new transcription start sites upstream of +1 in the promoter of the hgl5 gene. Mutation of the Inr element resulted in no change in the site of transcription initiation; however, in the presence of a mutated TATA and GAAC regions, the Inr element controlled the site of transcription initiation. We conclude that all three elements play a role in determining the site of transcription initiation. The variable position of the GAAC element relative to the site of transcription initiation, and the multiple transcription initiations that resulted from its mutation, indicate that the GAAC element has an important and apparently novel role in transcriptional control in E. histolytica.

Three novel families of miniature inverted-repeat transposable elements are associated with genes of the yellow fever mosquito, Aedes aegypti

Three novel families of transposable elements, Wukong, Wujin, and Wuneng, are described in the yellow fever mosquito, Aedes aegypti. Their copy numbers range from 2,100 to 3,000 per haploid genome. There are high degrees of sequence similarity within each family, and many structural but not sequence similarities between families. The common structural characteristics include small size, no coding potential, terminal inverted repeats, potential to form a stable secondary structure, A+T richness, and putative 2- to 4-bp A+T-biased specific target sites. Evidence of previous mobility is presented for the Wukong elements. Elements of these three families are associated with 7 of 16 fully or partially sequenced Ae. aegypti genes. Characteristics of these mosquito elements indicate strong similarities to the miniature inverted-repeat transposable elements (MITEs) recently found to be associated with plant genes. MITE-like elements have also been reported in two species of Xenopus and in Homo sapiens. This characterization of multiple families of highly repetitive MITE-like elements in an invertebrate extends the range of these elements in eukaryotic genomes. A hypothesis is presented relating genome size and organization to the presence of highly reiterated MITE families. The association of MITE-like elements with Ae. aegypti genes shows the same bias toward noncoding regions as in plants. This association has potentially important implications for the evolution of gene regulation.

Functional specificity of MutL homologs in yeast: Evidence for three Mlh1-based heterocomplexes with distinct roles during meiosis in recombination and mismatch correction

The yeast genome encodes four proteins (Pms1 and Mlh1–3) homologous to the bacterial mismatch repair component, MutL. Using two hybrid-interaction and coimmunoprecipitation studies, we show that these proteins can form only three types of complexes in vivo. Mlh1 is the common component of all three complexes, interacting with Pms1, Mlh2, and Mlh3, presumptively as heterodimers. The phenotypes of single deletion mutants reveal distinct functions for the three heterodimers during meiosis: in a pms1 mutant, frequent postmeiotic segregation indicates a defect in the correction of heteroduplex DNA, whereas the frequency of crossing-over is normal. Conversely, crossing-over in the mlh3 mutant is reduced to ≈70% of wild-type levels but correction of heteroduplex is normal. In a mlh2 mutant, crossing-over is normal and postmeiotic segregation is not observed but non-Mendelian segregation is elevated and altered with respect to parity. Finally, to a first approximation, the mlh1 mutant represents the combined single mutant phenotypes. Taken together, these data imply modulation of a basic Mlh1 function via combination with the three other MutL homologs and suggest specifically that Mlh1 combines with Mlh3 to promote meiotic crossing-over.

Analysis of three structurally related antiviral compounds in complex with human rhinovirus 16

Rhinoviruses are a frequent cause of the common cold. A series of antirhinoviral compounds have been developed that bind into a hydrophobic pocket in the viral capsid, stabilizing the capsid and interfering with cell attachment. The structures of a variety of such compounds, complexed with rhinovirus serotypes 14, 16, 1A, and 3, previously have been examined. Three chemically similar compounds, closely related to a drug that is undergoing phase III clinical trials, were chosen to determine the structural impact of the heteroatoms in one of the three rings. The compounds were found to have binding modes that depend on their electronic distribution. In the compound with the lowest efficacy, the terminal ring is displaced by 1 Å and rotated by 180° relative to the structure of the other two. The greater polarity of the terminal ring in one of the three compounds leads to a small displacement of its position relative to the other compounds in the hydrophobic end of the antiviral compound binding pocket to a site where it makes fewer interactions. Its lower efficacy is likely to be the result of the reduced number of interactions. A region of conserved residues has been identified near the entrance to the binding pocket where there is a corresponding conservation of the mode of binding of these compounds to different serotypes. Thus, variations in residues lining the more hydrophobic end of the pocket are primarily responsible for the differences in drug efficacies.

The interacting domains of three MutL heterodimers in man: hMLH1 interacts with 36 homologous amino acid residues within hMLH3, hPMS1 and hPMS2

In human cells, hMLH1, hMLH3, hPMS1 and hPMS2 are four recognised and distinctive homologues of MutL, an essential component of the bacterial DNA mismatch repair (MMR) system. The hMLH1 protein forms three different heterodimers with one of the other MutL homologues. As a first step towards functional analysis of these molecules, we determined the interacting domains of each heterodimer and tried to understand their common features. Using a yeast two-hybrid assay, we show that these MutL homologues can form heterodimers by interacting with the same amino acid residues of hMLH1, residues 492–742. In contrast, three hMLH1 partners, hMLH3, hPMS1 and hPMS2 contain the 36 homologous amino acid residues that interact strongly with hMLH1. Contrary to the previous studies, these homologous residues reside at the N-terminal regions of three subdomains conserved in MutL homologues in many species. Interestingly, these residues in hPMS2 and hMLH3 may form coiled-coil structures as predicted by the MULTICOIL program. Furthermore, we show that there is competition for the interacting domain in hMLH1 among the three other MutL homologues. Therefore, the quantitative balance of these three MutL heterodimers may be important in their functions.

Prediction of protein deamidation rates from primary and three-dimensional structure

A method for the quantitative estimation of instability with respect to deamidation of the asparaginyl (Asn) residues in proteins is described. The procedure involves the observation of several simple aspects of the three-dimensional environment of each Asn residue in the protein and a calculation that includes these observations, the primary amino acid residue sequence, and the previously reported complete set of sequence-dependent rates of deamidation for Asn pentapeptides. This method is demonstrated and evaluated for 23 proteins in which 31 unstable and 167 stable Asn residues have been reported and for 7 unstable and 63 stable Asn residues that have been reported in 61 human hemoglobin variants. The relative importance of primary structure and three-dimensional structure in Asn deamidation is estimated.

Differential Expression of Three Members of the 1-Aminocyclopropane-1-Carboxylate Synthase Gene Family in Carnation1

We investigated the expression patterns of three 1-aminocyclopropane-1-carboxylate (ACC) synthase genes in carnation (Dianthus caryophyllus cv White Sim) under conditions previously shown to induce ethylene biosynthesis. These included treatment of flowers with 2,4-dichlorophenoxyacetic acid, ethylene, LiCl, cycloheximide, and natural and pollination-induced flower senescence. Accumulation of ACC synthase transcripts in leaves following mechanical wounding and treatment with 2,4-dichlorophenoxyacetic acid or LiCl was also determined by RNA gel-blot analysis. As in other species, the carnation ACC synthase genes were found to be differentially regulated in a tissue-specific manner. DCACS2 and DCACS3 were preferentially expressed in styles, whereas DCACS1 mRNA was most abundant in petals. Cycloheximide did not induce increased accumulation of ACC synthase transcripts in carnation flowers, whereas the expression of ACC synthase was up-regulated by auxin, ethylene, LiCl, pollination, and senescence in a floral-organ-specific manner. Expression of the three ACC synthases identified in carnation did not correspond to elevated ethylene biosynthesis from wounded or auxin-treated leaves, and there are likely additional members of the carnation ACC synthase gene family responsible for ACC synthase expression in vegetative tissues.

Regulation of Oleoresinosis in Grand Fir (Abies grandis)1 : Differential Transcriptional Control of Monoterpene, Sesquiterpene, and Diterpene Synthase Genes in Response to Wounding

Grand fir (Abies grandis Lindl.) has been developed as a model system for the study of wound-induced oleoresinosis in conifers as a response to insect attack. Oleoresin is a roughly equal mixture of turpentine (85% monoterpenes [C10] and 15% sesquiterpenes [C15]) and rosin (diterpene [C20] resin acids) that acts to seal wounds and is toxic to both invading insects and their pathogenic fungal symbionts. The dynamic regulation of wound-induced oleoresin formation was studied over 29 d at the enzyme level by in vitro assay of the three classes of synthases directly responsible for the formation of monoterpenes, sesquiterpenes, and diterpenes from the corresponding C10, C15, and C20 prenyl diphosphate precursors, and at the gene level by RNA-blot hybridization using terpene synthase class-directed DNA probes. In overall appearance, the shapes of the time-course curves for all classes of synthase activities are similar, suggesting coordinate formation of all of the terpenoid types. However, closer inspection indicates that the monoterpene synthases arise earlier, as shown by an abbreviated time course over 6 to 48 h. RNA-blot analyses indicated that the genes for all three classes of enzymes are transcriptionally activated in response to wounding, with the monoterpene synthases up-regulated first (transcripts detectable 2 h after wounding), in agreement with the results of cell-free assays of monoterpene synthase activity, followed by the coordinately regulated sesquiterpene synthases and diterpene synthases (transcription beginning on d 3–4). The differential timing in the production of oleoresin components of this defense response is consistent with the immediate formation of monoterpenes to act as insect toxins and their later generation at solvent levels for the mobilization of resin acids responsible for wound sealing.

Identification of three physically and functionally distinct binding sites for C3b in human complement factor H by deletion mutagenesis.

Human complement factor H controls spontaneous activation of complement in plasma and appears to play a role in distinguishing host cells from activators of the alternative pathway of complement. In both mice and humans, the protein is composed of 20 homologous short consensus repeat (SCR) domains. The size of the protein suggests that portions of the structure outside the known C3b binding site (SCR 1-4) possess a significant biological role. We have expressed the full-length cDNA of factor H in the baculovirus system and have shown the recombinant protein to be fully active. Mutants of this full-length protein have now been prepared, purified, and examined for cofactor activity and binding to C3b and heparin. The results demonstrate (i) that factor H has at least three sites that bind C3b, (ii) that one of these sites is located in SCR domains 1-4, as has been shown by others, (iii) that a second site exists in the domain 6-10 region, (iv) that a third site resides in the SCR 16-20 region, and (v) that two heparin binding sites exist in factor H, one near SCR 13 and another in the SCR 6-10 region. Functional assays demonstrated that only the first C3b site located in SCR 1-4 expresses factor I cofactor activity. Mutant proteins lacking any one of the three C3b binding sites exhibited 6- to 8-fold reductions in affinity for C3b on sheep erythrocytes, indicating that all three sites contribute to the control of complement activation on erythrocytes. The identification of multiple functionally distinct sites on factor H clarifies many of the heretofore unexplainable behaviors of this protein, including the heterogeneous binding of factor H to surface-bound C3b, the effects of trypsin cleavage, and the differential control of complement activation on activators and nonactivators of the alternative pathway of complement.

Peptide methionine sulfoxide reductase contributes to the maintenance of adhesins in three major pathogens.

Pathogenic bacteria rely on adhesins to bind to host tissues. Therefore, the maintenance of the functional properties of these extracellular macromolecules is essential for the pathogenicity of these microorganisms. We report that peptide methionine sulfoxide reductase (MsrA), a repair enzyme, contributes to the maintenance of adhesins in Streptococcus pneumoniae, Neisseria gonorrhoeae, and Escherichia coli. A screen of a library of pneumococcal mutants for loss of adherence uncovered a MsrA mutant with 75% reduced binding to GalNAcbeta1-4Gal containing eukaryotic cell receptors that are present on type II lung cells and vascular endothelial cells. Subsequently, it was shown that an E. coli msrA mutant displayed decreased type I fimbriae-mediated, mannose-dependent, agglutination of erythrocytes. Previous work [Taha, M. K., So, M., Seifert, H. S., Billyard, E. & Marchal, C. (1988) EMBO J. 7, 4367-4378] has shown that mutants with defects in the pilA-pilB locus from N. gonorrhoeae were altered in their production of type IV pili. We show that pneumococcal MsrA and gonococcal PilB expressed in E. coli have MsrA activity. Together these data suggest that MsrA is required for the proper expression or maintenance of functional adhesins on the surfaces of these three major pathogenic bacteria.