Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the “euagarics.” Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.

(4-Aminomethyl)phenylguanidine derivatives as nonpeptidic highly selective inhibitors of human urokinase

Increased expression of the serine protease urokinase-type plasminogen activator (uPA) in tumor tissues is highly correlated with tumor cell migration, invasion, proliferation, progression, and metastasis. Thus inhibition of uPA activity represents a promising target for antimetastatic therapy. So far, only the x-ray crystal structure of uPA inactivated by H-Glu-Gly-Arg-chloromethylketone has been reported, thus limited data are available for a rational structure-based design of uPA inhibitors. Taking into account the trypsin-like arginine specificity of uPA, (4-aminomethyl)phenylguanidine was selected as a potential P1 residue and iterative derivatization of its amino group with various hydrophobic residues, and structure–activity relationship-based optimization of the spacer in terms of hydrogen bond acceptor/donor properties led to N-(1-adamantyl)-N′-(4-guanidinobenzyl)urea as a highly selective nonpeptidic uPA inhibitor. The x-ray crystal structure of the uPA B-chain complexed with this inhibitor revealed a surprising binding mode consisting of the expected insertion of the phenylguanidine moiety into the S1 pocket, but with the adamantyl residue protruding toward the hydrophobic S1′ enzyme subsite, thus exposing the ureido group to hydrogen-bonding interactions. Although in this enzyme-bound state the inhibitor is crossing the active site, interactions with the catalytic residues Ser-195 and His-57 are not observed, but their side chains are spatially displaced for steric reasons. Compared with other trypsin-like serine proteases, the S2 and S3/S4 pockets of uPA are reduced in size because of the 99-insertion loop. Therefore, the peculiar binding mode of the new type of uPA inhibitors offers the possibility of exploiting optimized interactions at the S1′/S2′ subsites to further enhance selectivity and potency. Because crystals of the uPA/benzamidine complex allow inhibitor exchange by soaking procedures, the structure-based design of new generations of uPA inhibitors can rely on the assistance of x-ray analysis.

Dynamic weighting in Monte Carlo and optimization

Dynamic importance weighting is proposed as a Monte Carlo method that has the capability to sample relevant parts of the configuration space even in the presence of many steep energy minima. The method relies on an additional dynamic variable (the importance weight) to help the system overcome steep barriers. A non-Metropolis theory is developed for the construction of such weighted samplers. Algorithms based on this method are designed for simulation and global optimization tasks arising from multimodal sampling, neural network training, and the traveling salesman problem. Numerical tests on these problems confirm the effectiveness of the method.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Structure-assisted design of mechanism-based irreversible inhibitors of human rhinovirus 3C protease with potent antiviral activity against multiple rhinovirus serotypes

Human rhinoviruses, the most important etiologic agents of the common cold, are messenger-active single-stranded monocistronic RNA viruses that have evolved a highly complex cascade of proteolytic processing events to control viral gene expression and replication. Most maturation cleavages within the precursor polyprotein are mediated by rhinovirus 3C protease (or its immediate precursor, 3CD), a cysteine protease with a trypsin-like polypeptide fold. High-resolution crystal structures of the enzyme from three viral serotypes have been used for the design and elaboration of 3C protease inhibitors representing different structural and chemical classes. Inhibitors having α,β-unsaturated carbonyl groups combined with peptidyl-binding elements specific for 3C protease undergo a Michael reaction mediated by nucleophilic addition of the enzyme’s catalytic Cys-147, resulting in covalent-bond formation and irreversible inactivation of the viral protease. Direct inhibition of 3C proteolytic activity in virally infected cells treated with these compounds can be inferred from dose-dependent accumulations of viral precursor polyproteins as determined by SDS/PAGE analysis of radiolabeled proteins. Cocrystal-structure-assisted optimization of 3C-protease-directed Michael acceptors has yielded molecules having extremely rapid in vitro inactivation of the viral protease, potent antiviral activity against multiple rhinovirus serotypes and low cellular toxicity. Recently, one compound in this series, AG7088, has entered clinical trials.

How optimization of potential functions affects protein folding.

The relationship between the optimization of the potential function and the foldability of theoretical protein models is studied based on investigations of a 27-mer cubic-lattice protein model and a more realistic lattice model for the protein crambin. In both the simple and the more complicated systems, optimization of the energy parameters achieves significant improvements in the statistical-mechanical characteristics of the systems and leads to foldable protein models in simulation experiments. The foldability of the protein models is characterized by their statistical-mechanical properties--e.g., by the density of states and by Monte Carlo folding simulations of the models. With optimized energy parameters, a high level of consistency exists among different interactions in the native structures of the protein models, as revealed by a correlation function between the optimized energy parameters and the native structure of the model proteins. The results of this work are relevant to the design of a general potential function for folding proteins by theoretical simulations.

Chaperonin-facilitated protein folding: optimization of rate and yield by an iterative annealing mechanism.

We develop a heuristic model for chaperonin-facilitated protein folding, the iterative annealing mechanism, based on theoretical descriptions of "rugged" conformational free energy landscapes for protein folding, and on experimental evidence that (i) folding proceeds by a nucleation mechanism whereby correct and incorrect nucleation lead to fast and slow folding kinetics, respectively, and (ii) chaperonins optimize the rate and yield of protein folding by an active ATP-dependent process. The chaperonins GroEL and GroES catalyze the folding of ribulose bisphosphate carboxylase at a rate proportional to the GroEL concentration. Kinetically trapped folding-incompetent conformers of ribulose bisphosphate carboxylase are converted to the native state in a reaction involving multiple rounds of quantized ATP hydrolysis by GroEL. We propose that chaperonins optimize protein folding by an iterative annealing mechanism; they repeatedly bind kinetically trapped conformers, randomly disrupt their structure, and release them in less folded states, allowing substrate proteins multiple opportunities to find pathways leading to the most thermodynamically stable state. By this mechanism, chaperonins greatly expand the range of environmental conditions in which folding to the native state is possible. We suggest that the development of this device for optimizing protein folding was an early and significant evolutionary event.

Protein structure-based design of potent orally bioavailable, nonpeptide inhibitors of human immunodeficiency virus protease.

A class of potent nonpeptidic inhibitors of human immunodeficiency virus protease has been designed by using the three-dimensional structure of the enzyme as a guide. By employing iterative protein cocrystal structure analysis, design, and synthesis the binding affinity of the lead compound was incrementally improved by over four orders of magnitude. An inversion in inhibitor binding mode was observed crystallographically, providing information critical for subsequent design and highlighting the utility of structural feedback in inhibitor optimization. These inhibitors are selective for the viral protease enzyme, possess good antiviral activity, and are orally available in three species.