Conservation of polyproline II helices in homologous proteins: implications for structure prediction by model building.

Left-handed polyproline II (PPII) helices commonly occur in globular proteins in segments of 4-8 residues. This paper analyzes the structural conservation of PPII-helices in 3 protein families: serine proteinases, aspartic proteinases, and immunoglobulin constant domains. Calculations of the number of conserved segments based on structural alignment of homologous molecules yielded similar results for the PPII-helices, the alpha-helices, and the beta-strands. The PPII-helices are consistently conserved at the level of 100-80% in the proteins with sequence identity above 20% and RMS deviation of structure alignments below 3.0 A. The most structurally important PPII segments are conserved below this level of sequence identity. These results suggest that the PPII-helices, in addition to the other 2 secondary structure classes, should be identified as part of structurally conserved regions in proteins. This is supported by similar values for the local RMS deviations of the aligned segments for the structural classes of PPII-helices, alpha-helices, and beta-strands. The PPII-helices are shown to participate in supersecondary elements such as PPII-helix/alpha-helix. The conservation of PPII-helices depends on the conservation of a supersecondary element as a whole. PPII-helices also form links, possibly flexible, in the interdomain regions. The role of the PPII-helices in model building by homology is 2-fold; they serve as additional conserved elements in the structure allowing improvement of the accuracy of a model and provide correct chain geometry for modeling of the segments equivalenced to them in a target sequence. The improvement in model building is demonstrated in 2 test studies.     

A “parallel plate” electrostatic model for bimolecular rate constants applied to electron transfer proteins

A “parallel plate” model describing the electrostatic potential energy of protein-protein interactions is presented that provides an analytical representation of the effect of ionic strength on a bimolecular rate constant. The model takes into account the asymmetric distribution of charge on the surface of the protein and localized charges at the site of electron transfer that are modeled as elements of a parallel plate condenser. Both monopolar and dipolar interactions are included. Examples of simple (monophasic) and complex (biphasic) ionic strength dependencies obtained from experiments with several electron transfer protein systems are presented, all of which can be accommodated by the model. The simple cases do not require the use of both monopolar and dipolar terms (i.e., they can be fit well by either alone). The biphasic dependencies can be fit only by using dipolar and monopolar terms of opposite sign, which is physically unreasonable for the molecules considered. Alternatively, the high ionic strength portion of the complex dependencies can be fit using either the monopolar term alone or the complete equation; this assumes a model in which such behavior is a consequence of electron transfer mechanisms involving changes in orientation or site of reaction as the ionic strength is varied. Based on these analyses, we conclude that the principal applications of the model presented here are to provide information about the structural properties of intermediate electron transfer complexes and to quantify comparisons between related proteins or site-specific mutants. We also conclude that the relative contributions of monopolar and dipolar effects to protein electron transfer kinetics cannot be evaluated from experimental data by present approximations.     

Homology models of two isozymes of manganese peroxidase: Prediction of a Mn(II) binding site

The three-dimensional structures of two isozymes of manganese peroxidase (MnP) have been predicted from homology modeling using lignin peroxidase as a template. Although highly homologous, MnP differs from LiP by the requirement of Mn(II) as an intermediate in its oxidation of substrates. The Mn(II) site is absent in LiP and unique to the MnP family of peroxidases. The model structures were used to identify the unique Mn(II) binding sites, to determine to what extent they were conserved in the two isozymes, and to provide insight into why this site is absent in LiP. For each isozyme of MnP, three candidate Mn(II) binding sites were identified. Energy optimizations of the three possible Mn(II) enzyme complexes allowed the selection of the most favorable Mn(II) binding site as one with the most anionic oxygen moieties best configured to act as ligands for the Mn(II). At the preferred site, the Mn(II) is coordinated to the carboxyl oxygens of Glu-35, Glu-39, and Asp-179, and a propionate group of the heme. The predicted Mn(II) binding site is conserved in both isozymes. Comparison between the residues at this site in MnP and the corresponding residues in LiP shows that two of the three anionic residues in MnP are replaced by neutral residues in LiP, explaining why LiP does not bind Mn(II). © 1994 Wiley-Liss, Inc.     

NAD(P)H-ubiquinone oxidoreductases in plant mitochondria

Plant (and fungal) mitochondria contain multiple NAD(P)H dehydrogenases in the inner membrane all of which are connected to the respiratory chain via ubiquinone. On the outer surface, facing the intermembrane space and the cytoplasm, NADH and NADPH are oxidized by what is probably a single low-molecular-weight, nonproton-pumping, unspecific rotenone-insensitive NAD(P)H dehydrogenase. Exogenous NADH oxidation is completely dependent on the presence of free Ca²⁺ with aK _0.5 of about 1 µM. On the inner surface facing the matrix there are two dehydrogenases: (1) the proton-pumping rotenone-sensitive multisubunit Complex I with properties similar to those of Complex I in mammalian and fungal mitochondria. (2) a rotenone-insensitive NAD(P)H dehydrogenase with equal activity with NADH and NADPH and no proton-pumping activity. The NADPH-oxidizing activity of this enzyme is completely dependent on Ca²⁺ with aK _0.5 of 3 µM. The enzyme consists of a single subunit of 26 kDa and has a native size of 76 kDa, which means that it may form a trimer.     

Interaction of the FAFB-containing subunit with the Photosystem 1 core heterodimer

The structure of the predicted amino acid sequence in the F_X domain of Photosystem 1 was studied by molecular modeling and a working hypothesis was developed for the functional interaction of PsaC with the core heterodimer. We propose that the intervening sequences between homologous cysteines in the F_X cluster form two flexible loops and participate in the binding of PsaC, and that the arginine residues in the two surface-exposed loops may promote the interaction between the P₇₀₀–F_X core and the subunit. The model was tested experimentally; chemical modification of arginine residues in the P₇₀₀–F_X core using phenylglyoxal prevented reconstitution of the core with PsaC and PsaD after insertion of FeS clusters in vitro. Treatment of the P₇₀₀–F_X core with trypsin also prevented reconstitution of terminal electron transfer to F_AF_B, although neither treatments affected the electron transfer to F_X as judged by flash kinetic spectrophotometry. Electron transfer in the P₇₀₀–F_AF_B complex was not impaired by either phenylglyoxal or trypsin treatment indicating that the small subunit(s) protect the arginine residues that become chemically modified or cleaved. The data are consistent with the working model and point to additional experiments designed to identify the specific residues involved in the interaction between the P₇₀₀–F_X core and PsaC.Abbreviations PG- phenylglyoxal - PS 1- Photosystem 1     

Investigation of subpopulation heterogeneity and plasmid stability in recombinant escherichia coli via a simple segregated model

Many microbial and cell cultures exhibit phenomena that can best be described using a segregated modeling approach. Heterogeneties are more marked in recombinant cell cultures because subpopulations, which often exhibit different growth and productivity characteristics, are more easily identified by selective markers. A simple segregated mathematical model that simulates the growth of recombinant Escherichia coli cells is developed. Subpopulations of different growth rate, plasmid replication rate, and plasmid segregation probability are explicitly considered. Results indicate that a third mechanism of plasmid instability, referred to here as a "downward selective pressure," is significant when describing plasmid loss in batch and chemostat cultures. Also, the model agrees well with experimental data from cultures under antibiotic selective pressure. Finally, model simulations of chemostat cultures reveal the importance of initial conditions on culture stability and the possible presence of nonrandom partitioning functions. (c) 1993 John Wiley & Sons, Inc.     

Strategies for selecting mutation sites for methionine enhancement in the bean seed storage protein phaseolin

The complete three-dimensional structure of the bean seed storage protein phaseolin was generated from -carbon coordinates by using molecular mechanic calculations. This structure was used as a template to simulate modifications aimed at increasing the methionine content of phaseolin. A hydrophilic, methionine-rich looping insert sequence was designed. Simulated mutagenesis shows that the insert might be accommodated in turn and loop regions of the protein, but not within an -helix. Methionine content was also increased by the replacement of hydrophobic amino acids with methionine in the central core -barrels of the phaseolin protein. Calculations indicated that methionine can effectively replace conserved or variant leucine, isolecuine, and valine residues. However, alanine residues were much more sensitive to substitution, and demonstrated high variability in the effects of methionine replacement. Introduction of multiple substitutions in the barrel interior demonstrated that the replaced residues could interact favorably to relieve local perturbations caused by individual substitutions. Molecular dynamics simulations were also utilized to study the structural organization of phaseolin. The calculations indicate that there are extensive packing interactions between the major domains of phaseolin, which have important implications for protein folding and stability. Since the proposed mutant proteins can be produced and studied, the results presented here provide an ideal test to determine if there is a correlation between the effects obtained by computer simulation and the effects of the mutations on the protein structure expressedin vivo.     

Human-modified landscapes driving the global primate extinction crisis

The world's primates have been severely impacted in diverse and profound ways by anthropogenic pressures. Here, we evaluate the impact of various infrastructures and human-modified landscapes on spatial patterns of primate species richness, at both global and regional scales. We overlaid the International Union for the Conservation of Nature (IUCN) range maps of 520 primate species and applied a global 100 km² grid. We used structural equation modeling and simultaneous autoregressive models to evaluate direct and indirect effects of six human-altered landscapes variables (i.e., human footprint [HFP], croplands [CROP], road density [ROAD], pasture lands [PAST], protected areas [PAs], and Indigenous Peoples' lands [IPLs]) on global primate species richness, threatened and non-threatened species, as well as on species with decreasing and non-decreasing populations. Two-thirds of all primate species are classified as threatened (i.e., Critically Endangered, Endangered, and Vulnerable), with ~86% experiencing population declines, and ~84% impacted by domestic or international trade. We found that the expansion of PAST, HFP, CROP, and road infrastructure had the most direct negative effects on primate richness. In contrast, forested habitat within IPLs and PAs was positively associated in safeguarding primate species diversity globally, with an even stronger effect at the regional level. Our results show that IPLs and PAs play a critical role in primate species conservation, helping to prevent their extinction; in contrast, HFP growth and expansion has a dramatically negative effect on primate species worldwide. Our findings support predictions that the continued negative impact of anthropogenic pressures on natural habitats may lead to a significant decline in global primate species richness, and likely, species extirpations. We advocate for stronger national and international policy frameworks promoting alternative/sustainable livelihoods and reducing persistent anthropogenic pressures to help mitigate the extinction risk of the world's primate species.