Channeling by Proximity: The Catalytic Advantages of Active Site Colocalization Using Brownian Dynamics

Nature often colocalizes successive steps in a metabolic pathway. Such organization is predicted to increase the effective concentration of pathway intermediates near their recipient active sites and to enhance catalytic efficiency. Here, the pathway of a two-step reaction is modeled using a simple spherical approximation for the enzymes and substrate particles. Brownian dynamics are used to simulate the trajectory of a substrate particle as it diffuses between the active site zones of two different enzyme spheres. The results approximate distances for the most effective reaction pathways, indicating that the most effective reaction pathway is one in which the active sites are closely aligned. However, when the active sites are too close, the ability of the substrate to react with the first enzyme was hindered, suggesting that even the most efficient orientations can be improved for a system that is allowed to rotate or change orientation to optimize the likelihood of reaction at both sites

The Native 3D Organization of Bacterial Polysomes

SummaryRecent advances have led to insights into the structure of the bacterial ribosome, but little is known about the 3D organization of ribosomes in the context of translating polysomes. We employed cryoelectron tomography and a template-matching approach to map 70S ribosomes in vitrified bacterial translation extracts and in lysates of active E. coli spheroplasts. In these preparations, polysomal arrangements were observed in which neighboring ribosomes are densely packed and exhibit preferred orientations. Analysis of characteristic examples of polysomes reveals a staggered or pseudohelical organization of ribosomes along the mRNA trace, with the transcript being sequestered on the inside, the tRNA entrance sites being accessible, and the polypeptide exit sites facing the cytosol. Modeling of elongating nascent polypeptide chains suggests that this arrangement maximizes the distance between nascent chains on adjacent ribosomes, thereby reducing the probability of intermolecular interactions that would give rise to aggregation and limit productive folding

Capturing the essence of folding and functions of biomolecules using Coarse-Grained Models

The distances over which biological molecules and their complexes can function range from a few nanometres, in the case of folded structures, to millimetres, for example during chromosome organization. Describing phenomena that cover such diverse length, and also time scales, requires models that capture the underlying physics for the particular length scale of interest. Theoretical ideas, in particular, concepts from polymer physics, have guided the development of coarse-grained models to study folding of DNA, RNA, and proteins. More recently, such models and their variants have been applied to the functions of biological nanomachines. Simulations using coarse-grained models are now poised to address a wide range of problems in biology.Comment: 37 pages, 8 figure

Influence of Nanoparticle Size and Shape on Oligomer Formation of an Amyloidogenic Peptide

Understanding the influence of macromolecular crowding and nanoparticles on the formation of in-register $\beta$-sheets, the primary structural component of amyloid fibrils, is a first step towards describing \emph{in vivo} protein aggregation and interactions between synthetic materials and proteins. Using all atom molecular simulations in implicit solvent we illustrate the effects of nanoparticle size, shape, and volume fraction on oligomer formation of an amyloidogenic peptide from the transthyretin protein. Surprisingly, we find that inert spherical crowding particles destabilize in-register $\beta$-sheets formed by dimers while stabilizing $\beta$-sheets comprised of trimers and tetramers. As the radius of the nanoparticle increases crowding effects decrease, implying smaller crowding particles have the largest influence on the earliest amyloid species. We explain these results using a theory based on the depletion effect. Finally, we show that spherocylindrical crowders destabilize the ordered $\beta$-sheet dimer to a greater extent than spherical crowders, which underscores the influence of nanoparticle shape on protein aggregation

Diffusion, Crowding & Protein Stability in a Dynamic Molecular Model of the Bacterial Cytoplasm

A longstanding question in molecular biology is the extent to which the behavior of macromolecules observed in vitro accurately reflects their behavior in vivo. A number of sophisticated experimental techniques now allow the behavior of individual types of macromolecule to be studied directly in vivo; none, however, allow a wide range of molecule types to be observed simultaneously. In order to tackle this issue we have adopted a computational perspective, and, having selected the model prokaryote Escherichia coli as a test system, have assembled an atomically detailed model of its cytoplasmic environment that includes 50 of the most abundant types of macromolecules at experimentally measured concentrations. Brownian dynamics (BD) simulations of the cytoplasm model have been calibrated to reproduce the translational diffusion coefficients of Green Fluorescent Protein (GFP) observed in vivo, and “snapshots” of the simulation trajectories have been used to compute the cytoplasm's effects on the thermodynamics of protein folding, association and aggregation events. The simulation model successfully describes the relative thermodynamic stabilities of proteins measured in E. coli, and shows that effects additional to the commonly cited “crowding” effect must be included in attempts to understand macromolecular behavior in vivo

Predicting active site residue annotations in the Pfam database

<p>Abstract</p> <p>Background</p> <p>Approximately 5% of Pfam families are enzymatic, but only a small fraction of the sequences within these families (<0.5%) have had the residues responsible for catalysis determined. To increase the active site annotations in the Pfam database, we have developed a strict set of rules, chosen to reduce the rate of false positives, which enable the transfer of experimentally determined active site residue data to other sequences within the same Pfam family.</p> <p>Description</p> <p>We have created a large database of predicted active site residues. On comparing our active site predictions to those found in UniProtKB, Catalytic Site Atlas, PROSITE and <it>MEROPS </it>we find that we make many novel predictions. On investigating the small subset of predictions made by these databases that are not predicted by us, we found these sequences did not meet our strict criteria for prediction. We assessed the sensitivity and specificity of our methodology and estimate that only 3% of our predicted sequences are false positives.</p> <p>Conclusion</p> <p>We have predicted 606110 active site residues, of which 94% are not found in UniProtKB, and have increased the active site annotations in Pfam by more than 200 fold. Although implemented for Pfam, the tool we have developed for transferring the data can be applied to any alignment with associated experimental active site data and is available for download. Our active site predictions are re-calculated at each Pfam release to ensure they are comprehensive and up to date. They provide one of the largest available databases of active site annotation.</p

Prediction of catalytic residues using Support Vector Machine with selected protein sequence and structural properties

BACKGROUND: The number of protein sequences deriving from genome sequencing projects is outpacing our knowledge about the function of these proteins. With the gap between experimentally characterized and uncharacterized proteins continuing to widen, it is necessary to develop new computational methods and tools for functional prediction. Knowledge of catalytic sites provides a valuable insight into protein function. Although many computational methods have been developed to predict catalytic residues and active sites, their accuracy remains low, with a significant number of false positives. In this paper, we present a novel method for the prediction of catalytic sites, using a carefully selected, supervised machine learning algorithm coupled with an optimal discriminative set of protein sequence conservation and structural properties. RESULTS: To determine the best machine learning algorithm, 26 classifiers in the WEKA software package were compared using a benchmarking dataset of 79 enzymes with 254 catalytic residues in a 10-fold cross-validation analysis. Each residue of the dataset was represented by a set of 24 residue properties previously shown to be of functional relevance, as well as a label {+1/-1} to indicate catalytic/non-catalytic residue. The best-performing algorithm was the Sequential Minimal Optimization (SMO) algorithm, which is a Support Vector Machine (SVM). The Wrapper Subset Selection algorithm further selected seven of the 24 attributes as an optimal subset of residue properties, with sequence conservation, catalytic propensities of amino acids, and relative position on protein surface being the most important features. CONCLUSION: The SMO algorithm with 7 selected attributes correctly predicted 228 of the 254 catalytic residues, with an overall predictive accuracy of more than 86%. Missing only 10.2% of the catalytic residues, the method captures the fundamental features of catalytic residues and can be used as a "catalytic residue filter" to facilitate experimental identification of catalytic residues for proteins with known structure but unknown function

Thermal Adaptation of Dihydrofolate Reductase from the Moderate ThermophileGeobacillus stearothermophilus

The thermal melting temperature of dihydrofolate reductase from Geobacillus stearothermophilus (BsDHFR) is 30 °C higher than that of its homologue from the psychrophile Moritella profunda. Additional proline residues in the loop regions of BsDHFR have been proposed to enhance the thermostability of BsDHFR, but site-directed mutagenesis studies reveal that these proline residues contribute only minimally. Instead, the high thermal stability of BsDHFR is partly due to removal of water-accessible thermolabile residues such as glutamine and methionine, which are prone to hydrolysis or oxidation at high temperatures. The extra thermostability of BsDHFR can be obtained by ligand binding, or in the presence of salts or cosolvents such as glycerol and sucrose. The sum of all these incremental factors allows BsDHFR to function efficiently in the natural habitat of G. stearothermophilus, which is characterized by temperatures that can reach 75 °C