PubMed Health⌕ Search

Biomedical subjects

Baldomero Oliva

Publications and source records attributed to Baldomero Oliva.

At least 19 recordsLinked to original sources

A supersecondary structure library and search algorithm for modeling loops in protein structures.

We present a fragment-search based method for predicting loop conformations in protein models. A hierarchical and multidimensional database has been set up that currently classifies 105,950 loop fragments and loop flanking secondary structures. Besides the length of the loops and types of bracing secondary structures the database is organized along four internal coordinates, a distance and three types of angles characterizing the geometry of stem regions. Candidate fragments are selected from this library by matching the length, the types of bracing secondary structures of the query and satisfying the geometrical restraints of the stems and subsequently inserted in the query protein framework where their fit is assessed by the root mean square deviation (r.m.s.d.) of stem regions and by the number of rigid body clashes with the environment. In the final step remaining candidate loops are ranked by a Z-score that combines information on sequence similarity and fit of predicted and observed phi/psi main chain dihedral angle propensities. Confidence Z-score cut-offs were determined for each loop length that identify those predicted fragments that outperform a competitive ab initio method. A web server implements the method, regularly updates the fragment library and performs prediction. Predicted segments are returned, or optionally, these can be completed with side chain reconstruction and subsequently annealed in the environment of the query protein by conjugate gradient minimization. The prediction method was tested on artificially prepared search datasets where all trivial sequence similarities on the SCOP superfamily level were removed. Under these conditions it is possible to predict loops of length 4, 8 and 12 with coverage of 98, 78 and 28% with at least of 0.22, 1.38 and 2.47 A of r.m.s.d. accuracy, respectively. In a head-to-head comparison on loops extracted from freshly deposited new protein folds the current method outperformed in a approximately 5:1 ratio an earlier developed database search method.

Algorithms↗

Study of the influence of temperature on the dynamics of the catalytic cleft in 1,3-1,4-beta-glucanase by molecular dynamics simulations.

The dependence of some molecular motions in the enzyme 1,3-1,4-beta-glucanase from Bacillus licheniformis on temperature changes and the role of the calcium ion in them were explored. For this purpose, two molecular dynamics simulated trajectories along 4 ns at low (300 K) and high (325 K) temperatures were generated by the GROMOS96 package. Several structural and thermodynamic parameters were calculated, including entropy values, solvation energies, and essential dynamics (ED). In addition, thermoinactivation experiments to study the influence of the calcium ion and some residues on the activity were conducted. The results showed the release of the calcium ion, which, in turn, significantly affected the movements of loops 1, 2, and 3, as shown by essential dynamics. These movements differ at low and high temperatures and affect dramatically the activity of the enzyme, as observed by thermoinactivation studies.

Bacillus↗

Red cell glucose phosphate isomerase (GPI): a molecular study of three novel mutations associated with hereditary nonspherocytic hemolytic anemia.

Molecular characteristics of red blood cell (RBC) glucose phosphate isomerase (GPI) deficiency are described in two Spanish patients with chronic nonspherocytic hemolytic anemia. One patient, with residual GPI activity in RBCs of around 7% (GPI-Catalonia), is homozygous for the missense mutation c.1648A>G (p.Lys550Glu) in exon 18. The other patient, with residual activity in RBCs of around 20% (GPI-Barcelona), was found to be a compound heterozygote for two different missense mutations: c.341A>T (p.Asp113Val) in exon 4 and c.663T>G (p.Asn220Lys) in exon 7. Molecular modeling using the human crystal structure of GPI as a model was performed to determine how these mutations could affect enzyme structure and function.

Adolescent↗

Prediction of the conformation and geometry of loops in globular proteins: testing ArchDB, a structural classification of loops.

In protein structure prediction, a central problem is defining the structure of a loop connecting 2 secondary structures. This problem frequently occurs in homology modeling, fold recognition, and in several strategies in ab initio structure prediction. In our previous work, we developed a classification database of structural motifs, ArchDB. The database contains 12,665 clustered loops in 451 structural classes with information about phi-psi angles in the loops and 1492 structural subclasses with the relative locations of the bracing secondary structures. Here we evaluate the extent to which sequence information in the loop database can be used to predict loop structure. Two sequence profiles were used, a HMM profile and a PSSM derived from PSI-BLAST. A jack-knife test was made removing homologous loops using SCOP superfamily definition and predicting afterwards against recalculated profiles that only take into account the sequence information. Two scenarios were considered: (1) prediction of structural class with application in comparative modeling and (2) prediction of structural subclass with application in fold recognition and ab initio. For the first scenario, structural class prediction was made directly over loops with X-ray secondary structure assignment, and if we consider the top 20 classes out of 451 possible classes, the best accuracy of prediction is 78.5%. In the second scenario, structural subclass prediction was made over loops using PSI-PRED (Jones, J Mol Biol 1999;292:195-202) secondary structure prediction to define loop boundaries, and if we take into account the top 20 subclasses out of 1492, the best accuracy is 46.7%. Accuracy of loop prediction was also evaluated by means of RMSD calculations.

Models, Molecular↗

Comparative modelling of protein structure and its impact on microbial cell factories.

Comparative modeling is becoming an increasingly helpful technique in microbial cell factories as the knowledge of the three-dimensional structure of a protein would be an invaluable aid to solve problems on protein production. For this reason, an introduction to comparative modeling is presented, with special emphasis on the basic concepts, opportunities and challenges of protein structure prediction. This review is intended to serve as a guide for the biologist who has no special expertise and who is not involved in the determination of protein structure. Selected applications of comparative modeling in microbial cell factories are outlined, and the role of microbial cell factories in the structural genomics initiative is discussed.

Journal Article↗

Detecting remotely related proteins by their interactions and sequence similarity.

The function of an uncharacterized protein is usually inferred either from its homology to, or its interactions with, characterized proteins. Here, we use both sequence similarity and protein interactions to identify relationships between remotely related protein sequences. We rely on the fact that homologous sequences share similar interactions, and, therefore, the set of interacting partners of the partners of a given protein is enriched by its homologs. The approach was bench-marked by assigning the fold and functional family to test sequences of known structure. Specifically, we relied on 1,434 proteins with known folds, as defined in the Structural Classification of Proteins (SCOP) database, and with known interacting partners, as defined in the Database of Interacting Proteins (DIP). For this subset, the specificity of fold assignment was increased from 54% for position-specific iterative BLAST to 75% for our approach, with a concomitant increase in sensitivity for a few percentage points. Similarly, the specificity of family assignment at the e-value threshold of 10(-8) was increased from 70% to 87%. The proposed method would be a useful tool for large-scale automated discovery of remote relationships between protein sequences, given its unique reliance on sequence similarity and protein-protein interactions.

Computational Biology↗

Characterization of the first described mutation of human red blood cell phosphoglycerate mutase.

In a patient with clinical diagnosis of Hereditary Spherocytosis and partial deficiency (50%) of red blood cell phosphoglycerate mutase (PGAM) activity, we have recently reported [A. Repiso, P. Pérez de la Ossa, X. Avilés, B. Oliva, J. Juncá, R. Oliva, E. Garcia, J.L.L. Vives-Corrons, J. Carreras, F. Climent, Red blood cell phosphoglycerate mutase. Description of the first human BB isoenzyme mutation, Haematologica 88 (2003) (03) ECR07] the first described mutation of type B PGAM subunit that as a dimer constitutes the PGAM (EC 5.4.2.1) isoenzyme present in red blood cells. The mutation is the substitution c.690G>A (p.Met230Ile). In this report, we show that the mutated PGAM possesses an abnormal behaviour on ion-exchange chromatography and is more thermo-labile that the native enzyme. We also confirm that, similar to the PGAM isoenzymes from other sources, the BB-PGAM from human erythrocytes has a ping pong or phosphoenzyme mechanism, and that the mutation does not significantly change the K(m) and K(i) values, and the optimum pH of the enzyme. The increased instability of the mutated enzyme can account for the decreased PGAM activity in patient's red blood cells. However, the implication of a change of the k(cat) produced by the mutation cannot be discarded, since we could not determine the k(cat) value of the mutated PGAM.

Adult↗

Glucose phosphate isomerase deficiency: enzymatic and familial characterization of Arg346His mutation.

Homozygous glucose phosphate isomerase (GPI) deficiency is one of the most important genetic disorders responsible for chronic non-spherocytic hemolytic anemia (CNSHA), a red blood cell autosomal recessive genetic disorder which causes severe metabolic alterations. In this work, we studied a patient with CNSHA due to an 82% loss of GPI activity resulting from the homozygous missense replacement in cDNA position 1040G>A, which leads to substitution of the protein residue A346H mutation. The enzyme is present in a dimeric form necessary for normal activity; the A346H mutation causes a loss of GPI capability to dimerize, which renders the enzyme more susceptible to thermolability and produces significant changes in erythrocyte metabolism.

Adolescent↗

Modeling the helicase domain of Brome mosaic virus 1a replicase.

Brome mosaic virus (BMV) is a representative member of positive-strand RNA viruses. The 1a replicase from BMV is a membrane protein of unknown structure with a methyltransferase N-terminal domain and a putative helicase activity in the C-terminal domain. In order to make a functional prediction of the helicase activity of the BMV 1a C-terminal domain, we have built a model of its structure. The use of fold recognition servers hinted at two different superfamilies of helicases [superfamily 1 (SF1) and superfamily 2 (SF2)] as putative templates for the C-terminal fragment of BMV 1a. A structural model of BMV 1a in SF2 was obtained by means of a fold recognition server (3D-PSSM). On the other hand, we used the helicase motifs described in the literature to construct a model of the structure of the BMV 1a C-terminal domain as a member of the SF1. The biological functionality and statistic potentials were used to discriminate between the two models. The results illustrate that the use of sequence profiles and patterns helps modeling. Accordingly, the C-terminal domain of BMV 1a is a potential member of the SF1 of helicases, and it can be modeled with the structure of a member of the UvrD family of helicases. The helicase mechanism was corroborated by the model and this supports the hypothesis that BMV 1a should have helicase activity.

Amino Acid Motifs↗

Positive selection in MAOA gene is human exclusive: determination of the putative amino acid change selected in the human lineage.

Monoamine oxidase A (MAOA) is the X-linked gene responsible for deamination and subsequent degradation of several neurotransmitters and other amines. Among other activities, the gene has been shown to play a role in locomotion, circadian rhythm, and pain sensitivity and to have a critical influence on behavior and cognition. Previous studies have reported a non-neutral evolution of the gene attributable to positive selection in the human lineage. To determine whether this selection was human-exclusive or shared with other species, we performed a population genetic analysis of the pattern of nucleotide variation in non-human species, including bonobo, chimpanzee, gorilla, and orangutan. Footprints of positive selection were absent in all analyzed species, suggesting that positive selection has been recent and unique to humans. To determine which human-unique genetic changes could have been responsible for this differential evolution, the coding region of the gene was compared between human, chimpanzee, and gorilla. Only one human exclusive non-conservative change is present in the gene: Glu151Lys. This human substitution affects protein dimerization according to a three-dimensional structural model that predicts a non-negligible functional shift. This is the only candidate position at present to have been selected to fixation in humans during an episode of positive selection. Divergence analysis among species has shown that, even under positive selection in the human lineage, the MAOA gene did not experience accelerated evolution in any of the analyzed lineages, and that tools such as K(a)/ K(s) would not have detected the selective history of the gene.

Amino Acid Sequence↗

Classification of common functional loops of kinase super-families.

A structural classification of loops has been obtained from a set of 141 protein structures classified as kinases. A total of 1813 loops was classified into 133 subclasses (9 betabeta(links), 15 betabeta(hairpins), 31 alpha-alpha, 46 alpha-beta and 32 beta-alpha). Functional information and specific features relating subclasses and function were included in the classification. Functional loops such as the P-loop (shared by different folds) or the Gly-rich-loop, among others, were classified into structural motifs. As a result, a common mechanism of catalysis and substrate binding was proved for most kinases. Additionally, the multiple-alignment of loop sequences made within each subclass was shown to be useful for comparative modeling of kinase loops. The classification is summarized in a kinase loop database located at http://sbi.imim.es/archki.

Amino Acid Motifs↗

ArchDB: automated protein loop classification as a tool for structural genomics.

The annotation of protein function has become a crucial problem with the advent of sequence and structural genomics initiatives. A large body of evidence suggests that protein structural information is frequently encoded in local sequences, and that folds are mainly made up of a number of simple local units of super-secondary structural motifs, consisting of a few secondary structures and their connecting loops. Moreover, protein loops play an important role in protein function. Here we present ArchDB, a classification database of structural motifs, consisting of one loop plus its bracing secondary structures. ArchDB currently contains 12,665 super-secondary elements classified into 1496 motif subclasses. The database provides an easy way to retrieve functional information from protein structures sharing a common motif, to search motifs found in a given SCOP family, superfamily or fold, or to search by keywords on proteins with classified loops. The ArchDB database of loops is located at http://sbi.imim.es/archdb.

Amino Acid Motifs↗

SNOW: standard nomenclature wizard to help searching for (bio) chemical standardized names.

UNLABELLED: When developing bioinformatical tools dealing with enzymatic activity, metabolism or enzymatic networks, the problem of the lack of a clear nomenclature for biochemical compounds often arises. This problem leads us to develop a small web-based tool (SNOW, Standard NOmenclature Wizard) which may help to find recommended and trivial names or the correct closest spelling for a query compound name, if it exists. AVAILABILITY: Web-based interface available at http://ibb.uab.es/snow/ SUPPLEMENTARY INFORMATION: http://ibb.uab.es/snow/snow_moreinfo.html

Algorithms↗

Molecular dynamics simulation of highly charged proteins: comparison of the particle-particle particle-mesh and reaction field methods for the calculation of electrostatic interactions.

Molecular dynamics (MD) simulations of the activation domain of porcine procarboxypeptidase B (ADBp) were performed to examine the effect of using the particle-particle particle-mesh (P3M) or the reaction field (RF) method for calculating electrostatic interactions in simulations of highly charged proteins. Several structural, thermodynamic, and dynamic observables were derived from the MD trajectories, including estimated entropies and solvation free energies and essential dynamics (ED). The P3M method leads to slightly higher atomic positional fluctuations and deviations from the crystallographic structure, along with somewhat lower values of the total energy and solvation free energy. However, the ED analysis of the system leads to nearly identical results for both simulations. Because of the strong similarity between the results, both methods appear well suited for the simulation of highly charged globular proteins in explicit solvent. However, the lower computational demand of the RF method in the present implementation represents a clear advantage over the P3M method.

Algorithms↗

Human procarboxypeptidase B: three-dimensional structure and implications for thrombin-activatable fibrinolysis inhibitor (TAFI).

Besides their classical role in alimentary protein degradation, zinc-dependant carboxypeptidases also participate in more selective regulatory processes like prohormone and neuropeptide processing or fibrinolysis inhibition in blood plasma. Human pancreatic procarboxypeptidase B (PCPB) is the prototype for those human exopeptidases that cleave off basic C-terminal residues and are secreted as inactive zymogens. One such protein is thrombin-activatable fibrinolysis inhibitor (TAFI), also known as plasma PCPB, which circulates in human plasma as a zymogen bound to plasminogen. The structure of human pancreatic PCPB displays a 95-residue pro-segment consisting of a globular region with an open-sandwich antiparallel-alpha antiparallel-beta topology and a C-terminal alpha-helix, which connects to the enzyme moiety. The latter is a 309-amino acid residue catalytic domain with alpha/beta hydrolase topology and a preformed active site, which is shielded by the globular domain of the pro-segment. The fold of the proenzyme is similar to previously reported procarboxypeptidase structures, also in that the most variable region is the connecting segment that links both globular moieties. However, the empty active site of human procarboxypeptidase B has two alternate conformations in one of the zinc-binding residues, which account for subtle differences in some of the key residues for substrate binding. The reported crystal structure, refined with data to 1.6A resolution, permits in the absence of an experimental structure, accurate homology modelling of TAFI, which may help to explain its properties.

Amino Acid Sequence↗

Human betacellulin structure modeled from other members of EGF family.

We have modeled betacellulin (BTC) to gain insight into the structural elements that can explain its properties. The epidermal growth factor (EGF) signal transduction pathway, a significant mediator of several cell functions, is based on four closely related tyrosine kinase receptors. The ErbB receptors are transmembrane glycoproteins and signal transduction is initiated by ligand binding that induces receptor homo- or heterodimerization to form a complex containing two molecules of ligand and two molecules of receptor. The EGF family of ligands can be divided into three groups based on their ability to bind and activate distinct ErbB receptor homo- and heterodimers. Each member of the group formed by BTC, heparin binding EGF (HB-EGF) and epiregulin (EP) can interact with both the EGF receptor (EGFR) and heregulin receptors (ErbB-3 and ErbB-4), and are hence called "bispecific" ligands. BTC and EP also present the distinctive feature that they activate all possible heterodimeric ErbB receptors. BTC has been modeled with the program MODELLER, using human EGF, human transforming growth factor alpha (hTGFalpha), human HB-EGF and human heregulin one alpha (hHRG-1alpha) as templates. The structure of the model as well as that of the templates were optimized and a simulation of 100 ps was run for all. The main structural properties of the model and the templates were compared and in conclusion the hBTC conformation was closely similar to that of hTGFalpha.

Amino Acid Sequence↗

TranScout: prediction of gene expression regulatory proteins from their sequences.

MOTIVATION: The advent of genomics yields thousands of reading frames in search of function. Identification of conserved functional motifs in protein sequences can be helpful for function prediction. RESULTS: A database and a classification of reported DNA-binding protein motifs has been designed. A program ('TranScout') has been developed for the detection and evaluation of conserved motifs in prokaryotic and eukaryotic sequences of proteins with a gene regulatory function. The efficiency of the program is shown in a benchmark against a database obtained from SWISS-PROT without the protein sequences used to train the program. All motifs were detected with a mean average sensitivity of 0.98 and a mean average specificity of 0.92. AVAILABILITY: The program is freely available for use on the internet at http://luz.uab.es/transcout/. The user can find additional information at this site.

Algorithms↗