Daniel R. Ripoll

Transposase-derived transcription factors regulate light signaling in Arabidopsis.

Plants use light to optimize growth and development. The photoreceptor phytochrome A (phyA) mediates various far-red light–induced responses. We show that Arabidopsis FHY3 and FAR1, which encode two proteins related to Mutator-like transposases, act together to modulate phyA signaling by directly activating the transcription of FHY1 and FHL, whose products are essential for light-induced phyA nuclear accumulation and subsequent light responses. FHY3 and FAR1 have separable DNA binding and transcriptional activation domains that are highly conserved in Mutator-like transposases. Further, expression of FHY3 and FAR1 is negatively regulated by phyA signaling. We propose that FHY3 and FAR1 represent transcription factors that have been co-opted from an ancient Mutator-like transposase(s) to modulate phyA-signaling homeostasis in higher plants.

Recent improvements in prediction of protein structure by global optimization of a potential energy function.

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsds are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Functional Dissection of Naturally Occurring Amino Acid Substitutions in eIF4E That Confers Recessive Potyvirus Resistance in Plants

Naturally existing variation in the eukaryotic translation initiation factor 4E (eIF4E) homolog encoded at the pvr1 locus in Capsicum results in recessively inherited resistance against several potyviruses. Previously reported data indicate that the physical interaction between Capsicum-eIF4E and the viral genome-linked protein (VPg) is required for the viral infection in the Capsicum-Tobacco etch virus (TEV) pathosystem. In this study, the potential structural role(s) of natural variation in the eIF4E protein encoded by recessive resistance alleles and their biological consequences have been assessed. Using high-resolution three-dimensional structural models based on the available crystallographic structures of eIF4E, we show that the amino acid substitution G107R, found in many recessive plant virus resistance genes encoding eIF4E, is predicted to result in a substantial modification in the protein binding pocket. The G107R change was shown to not only be responsible for the interruption of VPg binding in planta but also for the loss of cap binding ability in vitro, the principal function of eIF4E in the host. Overexpression of the Capsicum-eIF4E protein containing the G107R amino acid substitution in Solanum lycopersicum indicated that this polymorphism alone is sufficient for the acquisition of resistance against several TEV strains.

Calculation of protein conformation by global optimization of a potential energy function.

A novel hierarchical approach to protein folding has been applied to compute the unknown structures of seven target proteins provided by CASP3. The approach is based exclusively on the global optimization of a potential energy function for a united‐residue model by conformational space annealing, followed by energy refinement using an all‐atom potential. Comparison of the submitted models for five globular proteins with the experimental structures shows that the conformations of large fragments (∼60 aa) were predicted with rmsds of 4.2–6.8 Å for the Cα atoms. Our lowest‐energy models for targets T0056 and T0061 were particularly successful, producing the correct fold of approximately 52% and 80% of the structures, respectively. These results support the thermodynamic hypothesis that protein structure can be computed solely by global optimization of a potential energy function for a given amino acid sequence. Proteins Suppl 1999;3:204–208.

Structural modules for receptor dimerization in the S-locus receptor kinase extracellular domain

The highly polymorphic S-locus receptor kinase (SRK) is the stigma determinant of specificity in the self-incompatibility response of the Brassicaceae. SRK spans the plasma membrane of stigma epidermal cells, and it is activated in an allele-specific manner on binding of its extracellular region (eSRK) to its cognate pollen coat-localized S-locus cysteine-rich (SCR) ligand. SRK, like several other receptor kinases, forms dimers in the absence of ligand. To identify domains in SRK that mediate ligand-independent dimerization, we assayed eSRK for self-interaction in yeast. We show that SRK dimerization is mediated by two regions in eSRK, primarily by a C-terminal region inferred by homology modeling/fold recognition techniques to assume a PAN_APPLE-like structure, and secondarily by a region containing a signature sequence of the S-domain gene family, which might assume an EGF-like structure. We also show that eSRK exhibits a marked preference for homodimerization over heterodimerization with other eSRK variants and that this preference is mediated by a small, highly variable region within the PAN_APPLE domain. Thus, the extensive polymorphism exhibited by the eSRK not only determines differential affinity toward the SCR ligand, as has been assumed thus far, but also underlies a previously unrecognized allelic specificity in SRK dimerization. We propose that preference for SRK homodimerization explains the codominance exhibited by a majority of SRKs in the typically heterozygous stigmas of self-incompatible plants, whereas an increased propensity for heterodimerization combined with reduced affinity of heterodimers for cognate SCRs might underlie the dominant–recessive or mutual weakening relationships exhibited by some SRK allelic pairs.

Atomically detailed folding simulation of the B domain of staphylococcal protein A from random structures

The conformational space of the 10–55 fragment of the B-domain of staphylococcal protein A has been investigated by using the electrostatically driven Monte Carlo (EDMC) method. The ECEPP/3 (empirical conformational energy program for peptides) force-field plus two different continuum solvation models, namely SRFOPT (Solvent Radii Fixed with atomic solvation parameters OPTimized) and OONS (Ooi, Oobatake, Némethy, and Scheraga solvation model), were used to describe the conformational energy of the chain. After an exhaustive search, starting from two different random conformations, three of four runs led to native-like conformations. Boltzmann-averaged root-mean-square deviations (RMSD) for all of the backbone heavy atoms with respect to the native structure of 3.35 Å and 4.54 Å were obtained with SRFOPT and OONS, respectively. These results show that the protein-folding problem can be solved at the atomic detail level by an ab initio procedure, starting from random conformations, with no knowledge except the amino acid sequence. To our knowledge, the results reported here correspond to the largest protein ever folded from a random conformation by an initial-value formulation with a full atomic potential, without resort to knowledge-based information.

Structure of the Glucanase Inhibitor Protein (GIP) Family from Phytophthora Species Suggests Coevolution with Plant Endo-β-1,3-Glucanases

During invasion of their plant hosts, species of the oomycete genus Phytophthora secrete glucanase inhibitor proteins (GIPs) into the plant apoplast, which bind and inhibit the activity of plant extracellular endo-beta-1,3-glucanases (EGases). GIPs show structural homology to the chymotrypsin class of serine proteases (SP) but lack proteolytic activity due to the absence of an intact catalytic triad and, thus, belong to a broader class of proteins called serine protease homologs (SPH). To study the evolutionary relationship between GIPs and functional SP, database searches were used to identify 48 GIP homologs in the P. sojae, P. ramorum, and P. infestans genomes, composing GIPs, SPH, and potentially functional SP. Analyses of P. infestans-inoculated tomato leaves showed that P. infestans GIPs and tomato EGases are present in the apoplast and form stable complexes in planta. Studies of the temporal expression of a four-membered GIP family from P. infestans (PiGIP1 to PiGIP4) further revealed that the genes show distinctly different patterns during an infection timecourse. Codon evolution analyses of GIP homologs identified several positively selected peptide sites and structural modeling revealed them to be in close proximity to rapidly evolving EGase residues, suggesting that the interaction between GIPs and EGases has the hallmarks of a coevolving molecular arms race.

The multiple-minima problem in the conformational analysis of polypeptides. III. An electrostatically driven Monte Carlo method: tests on enkephalin.

The three-dimensional conformation of Met-enkephalin, corresponding to the lowest minimum of the empirical potential energy function ECEPP/2 (empirical conformational energy program for peptides), has been determined using a new algorithm, viz. the Electrostatically Driven Monte Carlo Method. This methodology assumes that a polypeptide or protein molecule is driven toward the native structure by the combined action of electrostatic interactions and stochastic conformational changes associated with thermal movements. These features are included in the algorithm that produces a Monte Carlo search in the conformational hyperspace of the polypeptide, using electrostatic predictions and a random sampling technique to locate low-energy conformations. In addition, we have incorporated an alternative mechanism that allows the structure to escape from some conformational regions representing metastable local energy minima and even from regions of the conformational space with great stability. In 33 test calculations on Met-enkephalin, starting from arbitrary or completely random conformations, the structure corresponding to the global energy minimum was found inall the cases analyzed, with a relatively small search of the conformational space. Some of these starting conformations wereright orleft-handed α-helices, characterized by good electrostatic interactions involving their backbone peptide dipoles; nevertheless, the procedure was able to convert such locally stable structures to the global-minimum conformation.

Tracing the Origin of the Fungal α1 Domain Places Its Ancestor in the HMG-Box Superfamily: Implication for Fungal Mating-Type Evolution

Background Fungal mating types in self-incompatible Pezizomycotina are specified by one of two alternate sequences occupying the same locus on corresponding chromosomes. One sequence is characterized by a gene encoding an HMG protein, while the hallmark of the other is a gene encoding a protein with an α1 domain showing similarity to the Matα1p protein of Saccharomyces cerevisiae. DNA-binding HMG proteins are ubiquitous and well characterized. In contrast, α1 domain proteins have limited distribution and their evolutionary origin is obscure, precluding a complete understanding of mating-type evolution in Ascomycota. Although much work has focused on the role of the S. cerevisiae Matα1p protein as a transcription factor, it has not yet been placed in any of the large families of sequence-specific DNA-binding proteins. Methodology/Principal Findings We present sequence comparisons, phylogenetic analyses, and in silico predictions of secondary and tertiary structures, which support our hypothesis that the α1 domain is related to the HMG domain. We have also characterized a new conserved motif in α1 proteins of Pezizomycotina. This motif is immediately adjacent to and downstream of the α1 domain and consists of a core sequence Y-[LMIF]-x(3)-G-[WL] embedded in a larger conserved motif. Conclusions/Significance Our data suggest that extant α1-box genes originated from an ancestral HMG gene, which confirms the current model of mating-type evolution within the fungal kingdom. We propose to incorporate α1 proteins in a new subclass of HMG proteins termed MATα_HMG.

A Tomato Endo-β-1,4-glucanase, SlCel9C1, Represents a Distinct Subclass with a New Family of Carbohydrate Binding Modules (CBM49)

A critical structural feature of many microbial endo-β-1,4-glucanases (EGases, or cellulases) is a carbohydrate binding module (CBM), which is required for effective crystalline cellulose degradation. However, CBMs are absent from plant EGases that have been biochemically characterized to date, and accordingly, plant EGases are not generally thought to have the capacity to degrade crystalline cellulose. We report the biochemical characterization of a tomato EGase, Solanum lycopersicum Cel8 (SlCel9C1), with a distinct C-terminal noncatalytic module that represents a previously uncharacterized family of CBMs. In vitro binding studies demonstrated that this module indeed binds to crystalline cellulose and can similarly bind as part of a recombinant chimeric fusion protein containing an EGase catalytic domain from the bacterium Thermobifida fusca. Site-directed mutagenesis studies show that tryptophans 559 and 573 play a role in crystalline cellulose binding. The SlCel9C1 CBM, which represents a new CBM family (CBM49), is a defining feature of a new structural subclass (Class C) of plant EGases, with members present throughout the plant kingdom. In addition, the SlCel9C1 catalytic domain was shown to hydrolyze artificial cellulosic polymers, cellulose oligosaccharides, and a variety of plant cell wall polysaccharides.