Atomically precise (catalytic) particles synthesized by a novel cluster deposition instrument.

We report a new high vacuum instrument which is dedicated to the preparation of well-defined clusters supported on model and technologically relevant supports for catalytic and materials investigations. The instrument is based on deposition of size selected metallic cluster ions that are produced by a high flux magnetron cluster source. The throughput of the apparatus is maximized by collecting and focusing ions utilizing a conical octupole ion guide and a linear ion guide. The size selection is achieved by a quadrupole mass filter. The new design of the sample holder provides for the preparation of multiple samples on supports of various sizes and shapes in one session. After cluster deposition onto the support of interest, samples will be taken out of the chamber for a variety of testing and characterization.

Evolutionary origins of the estrogen signaling system: insights from amphioxus.

Classically, the estrogen signaling system has two core components: cytochrome P450 aromatase (CYP19), the enzyme complex that catalyzes the rate limiting step in estrogen biosynthesis; and estrogen receptors (ERs), ligand activated transcription factors that interact with the regulatory region of target genes to mediate the biological effects of estrogen. While the importance of estrogens for regulation of reproduction, development and physiology has been well-documented in gnathostome vertebrates, the evolutionary origins of estrogen as a hormone are still unclear. As invertebrates within the phylum Chordata, cephalochordates (e.g., the amphioxus of the genus Branchiostoma) are among the closest invertebrate relatives of the vertebrates and can provide critical insight into the evolution of vertebrate-specific molecules and pathways. To address this question, this paper briefly reviews relevant earlier studies that help to illuminate the history of the aromatase and ER genes, with a particular emphasis on insights from amphioxus and other invertebrates. We then present new analyses of amphioxus aromatase and ER sequence and function, including an in silico model of the amphioxus aromatase protein, and CYP19 gene analysis. CYP19 shares a conserved gene structure with vertebrates (9 coding exons) and moderate sequence conservation (40% amino acid identity with human CYP19). Modeling of the amphioxus aromatase substrate binding site and simulated docking of androstenedione in comparison to the human aromatase shows that the substrate binding site is conserved and predicts that androstenedione could be a substrate for amphioxus CYP19. The amphioxus ER is structurally similar to vertebrate ERs, but differs in sequence and key residues of the ligand binding domain. Consistent with results from other laboratories, amphioxus ER did not bind radiolabeled estradiol, nor did it modulate gene expression on an estrogen-responsive element (ERE) in the presence of estradiol, 4-hydroxytamoxifen, diethylstilbestrol, bisphenol A or genistein. Interestingly, it has been shown that a related gene, the amphioxus "steroid receptor" (SR), can be activated by estrogens and that amphioxus ER can repress this activation. CYP19, ER and SR are all primarily expressed in gonadal tissue, suggesting an ancient paracrine/autocrine signaling role, but it is not yet known how their expression is regulated and, if estrogen is actually synthesized in amphioxus, whether it has a role in mediating any biological effects. Functional studies are clearly needed to link emerging bioinformatics and in vitro molecular biology results with organismal physiology to develop an understanding of the evolution of estrogen signaling. This article is part of a Special Issue entitled 'Marine organisms'.

Increased silver activity for direct propylene epoxidation via subnanometer size effects.

Production of the industrial chemical propylene oxide is energy-intensive and environmentally unfriendly. Catalysts based on bulk silver surfaces with direct propylene epoxidation by molecular oxygen have not resolved these problems because of substantial formation of carbon dioxide. We found that unpromoted, size-selected Ag3 clusters and approximately 3.5-nanometer Ag nanoparticles on alumina supports can catalyze this reaction with only a negligible amount of carbon dioxide formation and with high activity at low temperatures. Density functional calculations show that, relative to extended silver surfaces, oxidized silver trimers are more active and selective for epoxidation because of the open-shell nature of their electronic structure. The results suggest that new architectures based on ultrasmall silver particles may provide highly efficient catalysts for propylene epoxidation.

Consensus alignment for reliable framework prediction in homology modeling.

MOTIVATION: Even the best sequence alignment methods frequently fail to correctly identify the framework regions for which backbones can be copied from the template into the target structure. Since the underprediction and, more significantly, the overprediction of these regions reduces the quality of the final model, it is of prime importance to attain as much as possible of the true structural alignment between target and template. RESULTS: We have developed an algorithm called Consensus that consistently provides a high quality alignment for comparative modeling. The method follows from a benchmark analysis of the 3D models generated by ten alignment techniques for a set of 79 homologous protein structure pairs. For 20-to-40% of the targets, these methods yield models with at least 6 A root mean square deviation (RMSD) from the native structure. We have selected the top five performing methods, and developed a consensus algorithm to generate an improved alignment. By building on the individual strength of each method, a set of criteria was implemented to remove the alignment segments that are likely to correspond to structurally dissimilar regions. The automated algorithm was validated on a different set of 48 protein pairs, resulting in 2.2 A average RMSD for the predicted models, and only four cases in which the RMSD exceeded 3 A. The average length of the alignments was about 75% of that found by standard structural superposition methods. The performance of Consensus was consistent from 2 to 32% target-template sequence identity, and hence it can be used for accurate prediction of framework regions in homology modeling.

Semiglobal simplex optimization and its application to determining the preferred solvation sites of proteins.

The classical simplex method is extended into the Semiglobal Simplex (SGS) algorithm. Although SGS does not guarantee finding the global minimum, it affords a much more thorough exploration of the local minima than any traditional minimization method. The basic idea of SGS is to perform a local minimization in each step of the simplex algorithm, and thus, similarly to the Convex Global Underestimator (CGU) method, the search is carried out on a surface spanned by local minima. The SGS and CGU methods are compared by minimizing a set of test functions of increasing complexity, each with a known global minimum and many local minima. Although CGU delivers substantially better success rates in simple problems, the two methods become comparable as the complexity of the problems increases. Because SGS is generally faster than CGU, it is the method of choice for solving optimization problems in which function evaluation is computationally inexpensive and the search region is large. The extreme simplicity of the method is also a factor. The SGS method is applied here to the problem of finding the most preferred (i.e., minimum free energy) solvation sites on a streptavidin monomer. It is shown that the SGS method locates the same lowest free energy positions as an exhaustive multistart Simplex search of the protein surface, with less than one-tenth the number of minizations. The combination of the two methods, i.e.. multistart simplex and SGS, provides a reliable procedure for predicting all potential solvation sites of a protein.

A streptavidin mutant useful for directed immobilization on solid surfaces.

A streptavidin mutant has been designed and produced that allows the specific, covalent immobilization of streptavidin on solid surfaces. This streptavidin mutant was constructed by fusing a six-residue sequence, containing a single cysteine, to the carboxyl terminus of streptavidin. Because this mutant has no other cysteine residues, the reactive sulfhydryl group of the cysteine residue serves as a unique immobilization site for conjugation using sulfhydryl chemistry. This streptavidin mutant was efficiently immobilized on maleimide-coated solid surfaces via its unique immobilization site. Characterization of the immobilized streptavidin mutant for the ability to bind to biotinylated macromolecules and the dissociation rates of bound biotin showed that the biotin-binding properties of this mutant were minimally affected by immobilization on solid surfaces. This streptavidin could be readily incorporated into a wide variety of solid-phase diagnostic tests and biomedical assays. This could enhance the performance of streptavidin-based solid-phase assay systems.

Protein docking along smooth association pathways.

We propose a docking method that mimics the way proteins bind. The method accounts for the dominant driving forces at the different length scales of the protein binding process, allowing for an efficient selection of a downhill path on the evolving receptor-ligand-free energy landscape. Starting from encounter complexes with as much as 10 A rms deviation from the native conformation, the method locally samples the six dimensional space of rigid-body receptor-ligand structures subject to a van der Waals constraint. The sampling is initially biased only by the desolvation and electrostatic components of the free energy, which capture the partial affinity of unbound structures that are more than 4 A away from the native state. Below this threshold, improved discrimination is attained by adding an increasing fraction of the van der Waals energy to the force field. The method, with no free parameters, was tested in eight different sets of independently crystallized receptor-ligand structures consistently predicting bound conformations with the lowest free energies and appropriate stability gap around 2 A from the native complex. This multistage approach is consistent with the underlying kinetics and internal structure of the free energy funnel to the bound state. Implications for the nature of the protein binding pathways are also discussed.

Dynamical view of the positions of key side chains in protein-protein recognition.

When a complex is constructed from the separately determined rigid structures of a receptor and its ligand, some key side chains are usually in wrong positions. These distortions of the interface yield an apparent loss in affinity and would unfavorably affect the kinetics of association. It is generally assumed that the interacting proteins should drive the appropriate conformational changes, leading to their complementarity, but this hypothesis does not explain their fast association rates. However, nanosecond explicit solvent molecular dynamics simulations of misfolded surface side chains from the independently solved structures of barstar, bovine pancreatic trypsin inhibitor, and lysozyme show that even before any receptor-ligand interaction, key side chains frequently visit the rotamer conformations seen in the complex. We show that these simple structural motifs can reconcile most of the binding affinity required for a rapid and highly specific association process. Side chains amenable to induced fit are also identified. These results corroborate that solvent-side chain interactions play a critical role in the recognition process. Our findings are also supported by crystallographic data.

Discrimination of near-native protein structures from misfolded models by empirical free energy functions.

Free energy potentials, combining molecular mechanics with empirical solvation and entropic terms, are used to discriminate native and near-native protein conformations from slightly misfolded decoys. Since the functional forms of these potentials vary within the field, it is of interest to determine the contributions of individual free energy terms and their combinations to the discriminative power of the potential. This is achieved in terms of quantitative measures of discrimination that include the correlation coefficient between RMSD and free energy, and a new measure labeled the minimum discriminatory slope (MDS). In terms of these criteria, the internal energy is shown to be a good discriminator on its own, which implies that even well-constructed decoys are substantially more strained than the native protein structure. The discrimination improves if, in addition to the internal energy, the free energy expression includes the electrostatic energy, calculated by assuming non-ionized side chains, and an empirical solvation term, with the classical atomic solvation parameter model providing slightly better discrimination than a structure-based atomic contact potential. Finally, the inclusion of a term representing the side chain entropy change, and calculated by an established empirical scale, is so inaccurate that it makes the discrimination worse. It is shown that both the correlation coefficient and the MDS value (or its dimensionless form) are needed for an objective assessment of a potential, and that together they provide much more information on the origins of discrimination than simple inspection of the RMSD-free energy plots.

Scoring docked conformations generated by rigid-body protein-protein docking.

Rigid-body methods, particularly Fourier correlation techniques, are very efficient for docking bound (co-crystallized) protein conformations using measures of surface complementarity as the target function. However, when docking unbound (separately crystallized) conformations, the method generally yields hundreds of false positive structures with good scores but high root mean square deviations (RMSDs). This paper describes a two-step scoring algorithm that can discriminate near-native conformations (with less than 5 A RMSD) from other structures. The first step includes two rigid-body filters that use the desolvation free energy and the electrostatic energy to select a manageable number of conformations for further processing, but are unable to eliminate all false positives. Complete discrimination is achieved in the second step that minimizes the molecular mechanics energy of the retained structures, and re-ranks them with a combined free-energy function which includes electrostatic, solvation, and van der Waals energy terms. After minimization, the improved fit in near-native complex conformations provides the free-energy gap required for discrimination. The algorithm has been developed and tested using docking decoys, i.e., docked conformations generated by Fourier correlation techniques. The decoy sets are available on the web for testing other discrimination procedures. Proteins 2000;40:525-537.

Kinetics of desolvation-mediated protein-protein binding.

The role of desolvation in protein binding kinetics is investigated using Brownian dynamics simulations in complexes in which the electrostatic interactions are relatively weak. We find that partial desolvation, modeled by a short-range atomic contact potential, is not only a major contributor to the binding free energy but also substantially increases the diffusion-limited rate for complexes in which long-range electrostatics is weak. This rate enhancement is mostly due to weakly specific pathways leading to a low free-energy attractor, i.e., a precursor state before docking. For alpha-chymotrypsin and human leukocyte elastase, both interacting with turkey ovomucoid third domain, we find that the forward rate constant associated with a collision within a solid angle phi around their corresponding attractor approaches 10(7) and 10(6) M(-1)s(-1), respectively, in the limit phi approximately 2 degrees. Because these estimates agree well with experiments, we conclude that the final bound conformation must be preceded by a small set of well-defined diffusion-accessible precursor states. The inclusion of the otherwise repulsive desolvation interaction also explains the lack of aggregation in proteins by restricting nonspecific association times to approximately 4 ns. Under the same reaction conditions but without short range forces, the association rate would be only approximately 10(3) M(-1)s(-1). Although desolvation increases these rates by three orders of magnitude, desolvation-mediated association is still at least 100-fold slower than the electrostatically assisted binding in complexes such as barnase and barstar.

Continuum electrostatic analysis of preferred solvation sites around proteins in solution.

To understand water-protein interactions in solution, the electrostatic field is calculated by solving the Poisson-Boltzmann equation, and the free energy surface of water is mapped by translating and rotating an explicit water molecule around the protein. The calculation is applied to T4 lysozyme with data available on the conservation of solvent binding sites in 18 crystallographically independent molecules. The free energy maps around the ordered water sites provide information on the relationship between water positions in crystal structure and in solution. Results show that almost all conserved sites and the majority of nonconserved sites are within 1.3 A of local free energy minima. This finding is in sharp contrast to the behavior of randomly placed water molecules in the boundary layer, which, on the average, must travel more than 3 A to the nearest free energy minimum. Thus, the solvation sites are at least partially determined by protein-water interactions rather than by crystal packing alone. The characteristic water residence times, obtained from the free energies at the local minima, are in good agreement with nuclear magnetic resonance experiments. Only about half of the potential sites show up as ordered water in the 1.7 A resolution X-ray structure. Crystal packing interactions can stabilize weak or mobile potential sites (in fact, some ordered water positions are not close to free energy minima) or can prevent water from occupying certain sites. Apart from a few buried water molecules that are strong binders, the free energies are not very different for conserved and nonconserved sites. We show that conservation of a water site between two crystals occurs if the positions of protein atoms, primarily contributing to the free energy at the local minimum, do not substantially change from one structure to the other. This requirement can be correlated with the nature of the side chain contacting the water molecule in the site.

Human transaldolase and cross-reactive viral epitopes identified by autoantibodies of multiple sclerosis patients.

Multiple sclerosis is mediated by an autoimmune process causing selective destruction of oligodendrocytes. Transaldolase, which is expressed in the brain selectively in oligodendrocytes, is a target of high affinity autoantibodies in serum and cerebrospinal fluid of multiple sclerosis patients. A three-dimensional model of human transaldolase was developed based on the crystal structure of the enzyme from Escherichia coli. To identify immunodominant epitopes, 33 peptides overlapping human transaldolase by 5 amino acids were synthesized. Ab 12484, raised against enzymatically active human transaldolase, recognized antigenic determinants corresponding to linear epitopes (residues 27-31 and 265-290) and alpha helices (residues 75-98 and 302-329). Four immunodominant peptides harboring charged amino acid residues with topographically exposed side chains were identified by sera from 13 multiple sclerosis patients with predetermined autoreactivity to transaldolase. Autoantibodies binding to the most prominent human transaldolase epitope, between residues 271 and 285, showed cross-reactivity with Epstein-Barr and herpes simplex virus type 1 capsid-derived peptides. Molecular mimicry between immunodominant autoepitopes and viral Ags may be a decisive factor in directing autoimmunity to transaldolase in multiple sclerosis patients.

Free energy landscapes of encounter complexes in protein-protein association.

We report the computer generation of a high-density map of the thermodynamic properties of the diffusion-accessible encounter conformations of four receptor-ligand protein pairs, and use it to study the electrostatic and desolvation components of the free energy of association. Encounter complex conformations are generated by sampling the translational/rotational space of the ligand around the receptor, both at 5-A and zero surface-to-surface separations. We find that partial desolvation is always an important effect, and it becomes dominant for complexes in which one of the reactants is neutral or weakly charged. The interaction provides a slowly varying attractive force over a small but significant region of the molecular surface. In complexes with no strong charge complementarity this region surrounds the binding site, and the orientation of the ligand in the encounter conformation with the lowest desolvation free energy is similar to the one observed in the fully formed complex. Complexes with strong opposite charges exhibit two types of behavior. In the first group, represented by barnase/barstar, electrostatics exerts strong orientational steering toward the binding site, and desolvation provides some added adhesion within the local region of low electrostatic energy. In the second group, represented by the complex of kallikrein and pancreatic trypsin inhibitor, the overall stability results from the rather nonspecific electrostatic attraction, whereas the affinity toward the binding region is determined by desolvation interactions.

A streptavidin mutant with altered ligand-binding specificity.

The biotin-binding site of streptavidin was modified to alter its ligand-binding specificity. In natural streptavidin, the side chains of N23 and S27 make two of the three hydrogen bonds with the ureido oxygen of biotin. These two residues were mutated to severely weaken biotin binding while attempting to maintain the affinity for two biotin analogs, 2-iminobiotin and diaminobiotin. Redesigning of the biotin-binding site used the difference in local electrostatic charge distribution between biotin and these biotin analogs. Free energy calculations predicted that the introduction of a negative charge at the position of S27 plus the mutation N23A should disrupt two of the three hydrogen bonds between natural streptavidin and the ureido oxygen of biotin. In contrast, the imino hydrogen of 2-iminobiotin should form a hydrogen bond with the side chain of an acidic amino acid at position 27. This should reduce the biotin-binding affinity by approximately eight orders of magnitude, while leaving the affinities for these biotin analogs virtually unaffected. In good agreement with these predictions, a streptavidin mutant with the N23A and S27D substitutions binds 2-iminobiotin with an affinity (Ka) of 1 x 10(6) M-1, two orders of magnitude higher than that for biotin (1 x 10(4) M-1). In contrast, the binding affinity of this streptavidin mutant for diaminobiotin (2.7 x 10(4) M-1) was lower than predicted (2.9 x 10(5) M-1), suggesting the position of the diaminobiotin in the biotin-binding site was not accurately determined by modeling.

Genetic engineering of streptavidin, a versatile affinity tag.

Streptavidin, a tetrameric protein produced by Streptomyces avidinii, has been used as a useful, versatile affinity tag in a variety of biological applications. The efficacy of streptavidin is derived from its extremely high binding affinity for the vitamin biotin. For the last several years, we have used genetic engineering as a primary means to enhance the properties of streptavidin and to expand the application of streptavidin as an affinity tag. In this review, we describe several genetically engineered streptavidin variants, which include a streptavidin with a reduced biotin-binding affinity, a dimeric streptavidin, and a fusion protein between streptavidin and protein A, along with their potential applications in biological science.

Selecting near-native conformations in homology modeling: the role of molecular mechanics and solvation terms.

A free energy function, combining molecular mechanics energy with empirical solvation and entropic terms, is used for ranking near-native conformations that occur in the conformational search steps of homology modeling, i.e., side-chain search and loop closure calculations. Correlations between the free energy and RMS deviation from the X-ray structure are established. It is shown that generally both molecular mechanics and solvation/entropic terms should be included in the potential. The identification of near-native backbone conformations is accomplished primarily by the molecular mechanics term that becomes the dominant contribution to the free energy if the backbone is even slightly strained, as frequently occurs in loop closure calculations. Both terms become equally important if a sufficiently accurate backbone conformation is found. Finally, the selection of the best side-chain positions for a fixed backbone is almost completely governed by the solvation term. The discriminatory power of the combined potential is demonstrated by evaluating the free energies of protein models submitted to the first meeting on Critical Assessment of techniques for protein Structure Prediction (CASP1), and comparing them to the free energies of the native conformations.

Empirical free energy calculation: comparison to calorimetric data.

An effective free energy potential, developed originally for binding free energy calculation, is compared to calorimetric data on protein unfolding, described by a linear combination of changes in polar and nonpolar surface areas. The potential consists of a molecular mechanics energy term calculated for a reference medium (vapor or nonpolar liquid), and empirical terms representing solvation and entropic effects. It is shown that, under suitable conditions, the free energy function agrees well with the calorimetric expression. An additional result of the comparison is an independent estimate of the side-chain entropy loss, which is shown to agree with a structure-based entropy scale. These findings confirm that simple functions can be used to estimate the free energy change in complex systems, and that a binding free energy evaluation model can describe the thermodynamics of protein unfolding correctly. Furthermore, it is shown that folding and binding leave the sum of solute-solute and solute-solvent van der Waals interactions nearly invariant and, due to this invariance, it may be advantageous to use a nonpolar liquid rather than vacuum as the reference medium.

Engineering subunit association of multisubunit proteins: a dimeric streptavidin.

A dimeric streptavidin has been designed by molecular modeling using effective binding free energy calculations that decompose the binding free energy into electrostatic, desolvation, and side chain entropy loss terms. A histidine-127 --> aspartic acid (H127D) mutation was sufficient to introduce electrostatic repulsion between subunits that prevents the formation of the natural tetramer. However, the high hydrophobicity of the dimer-dimer interface, which would be exposed to solvent in a dimeric streptavidin, suggests that the resulting molecule would have very low solubility in aqueous media. In agreement with the calculations, a streptavidin containing the H127D mutation formed insoluble aggregates. Thus, the major design goal was to reduce the hydrophobicity of the dimer-dimer interface while maintaining the fundamental structure. Free energy calculations suggested that the hydrophobicity of the dimer-dimer interface could be reduced significantly by deleting a loop from G113 through W120 that should have no apparent contact with biotin in a dimeric molecule. The resulting protein, containing both the H127D mutation and the loop deletion, formed a soluble dimeric streptavidin in the presence of biotin.

Beta-turn new classification and its some features in proteins.

An inspection of the phi-psi angle distribution strongly suggests that protein folding is highly constrained. A number of researchers have even suggested that a relatively small set of discrete phi-psi regions might be sufficient to describe most protein conformation. The total of 541 tight turns from 101 non-identical proteins were extracted form Brookhaven DataBank. The dihedral values of tight turns were scattered into the seven regions on the Ramachandran plot. These seven regions were called A1, A2, B1, B2, B22, T1 and T2. A1 and A2 are the traditional alpha-helix regions, B1, B2 and B22 the beta-strand regions, T1 and T2 the beta-turn regions. The A2 and T2 regions were not defined as "discrete" or single points but rather as one dimensional extended states. Based on the geometry of the two central residues of the tight turns, the new classification of beta-turn was defined. This classification of the majority of beta-turns fell into only six of the possible forty nine region combinations and were identifiable with the traditional nomenclature of Venkatachalam(1), but much simpler. The function of beta-turn in the conformation of proteins was studied. The hydrophobicity for different type turns was discussed. It shows that beta-turns have very strong hydrophilic property, so they are usually situated at the folding protein surface. The features of beta-turn and its amino acid distribution in this 541 beta-turn group and different type beta-turn were given.