RESUMEN
Measurement errors occur very commonly in practice. After fitting the model, influence diagnostics is an important step in statistical data analysis. The most frequently used diagnostic method for measurement error models is the local influence. However, this methodology may fail to detect masked influential observations. To overcome this limitation, we propose the use of the conformal normal curvature with the forward search algorithm. The results are presented through easy to interpret plots considering different perturbation schemes. The proposed methodology is illustrated with three real data sets and one simulated data set, two of which have been previously analyzed in the literature. The third data set deals with the stability of the hygroscopic solid dosage in pharmaceutical processes to ensure the maintenance of product safety quality. In this application, the analytical mass balance is subject to measurement errors, which require attention in the modeling process and diagnostic analysis.
RESUMEN
Proficiency testing (PT) determines the performance of individual laboratories for specific tests or measurements and it is used to monitor the reliability of laboratories measurements. PT plays a highly valuable role as it provides objective evidence of the competence of the participant laboratories. In this paper, we propose a multivariate calibration model to assess equivalence among laboratories measurements in PT. Our method allows to deal with multivariate data, where the item under test is measured at different levels. Although intuitive, the proposed model is nonergodic, which means that the asymptotic Fisher information matrix is random. As a consequence, a detailed asymptotic analysis was carried out to establish the strategy for comparing the results of the participating laboratories. To illustrate, we apply our method to analyze the data from the Brazilian engine test group, PT program, where the power of an engine was measured by eight laboratories at several levels of rotation.
RESUMEN
Understanding the functional effect of Single Amino acid Substitutions (SAS), derived from the occurrence of single nucleotide variants (SNVs), and their relation to disease development is a major issue in clinical genomics. Despite the existence of several bioinformatic algorithms and servers that predict if a SAS is pathogenic or not, they give little or no information at all on the reasons for pathogenicity prediction and on the actual predicted effect of the SAS on the protein function. Moreover, few actual methods take into account structural information when available for automated analysis. Moreover, many of these algorithms are able to predict an effect that no necessarily translates directly into pathogenicity. VarQ is a bioinformatic pipeline that incorporates structural information for the detailed analysis and prediction of SAS effect on protein function. It is an online tool which uses UniProt id and automatically analyzes known and user provided SAS for their effect on protein activity, folding, aggregation and protein interactions, among others. We show that structural information, when available, can improve the SAS pathogenicity diagnosis and more important explain its causes. We show that VarQ is able to correctly reproduce previous analysis of RASopathies related mutations, saving extensive and time consuming manual curation. VarQ assessment was performed over a set of previously manually curated RASopathies (diseases that affects the RAS/MAPK signaling pathway) related variants, showing its ability to correctly predict the phenotypic outcome and its underlying cause. This resource is available online at http://varq.qb.fcen.uba.ar/. Supporting Information & Tutorials may be found in the webpage of the tool.
RESUMEN
BACKGROUND: Congenital Hypopituitarism is caused by genetic and environmental factors. Over 30 genes have been implicated in isolated and/or combined pituitary hormone deficiency. The etiology remains unknown for up to 80% of the patients, but most cases have been analyzed by limited candidate gene screening. Mutations in the PROP1 gene are the most common known cause, and the frequency of mutations in this gene varies greatly by ethnicity. We designed a custom array to assess the frequency of mutations in known hypopituitarism genes and new candidates, using single molecule molecular inversion probes sequencing (smMIPS). METHODS: We used this panel for the first systematic screening for causes of hypopituitarism in children. Molecular inversion probes were designed to capture 693 coding exons of 30 known genes and 37 candidate genes. We captured genomic DNA from 51 pediatric patients with CPHD (n = 43) or isolated GH deficiency (IGHD) (n = 8) and their parents and conducted next generation sequencing. RESULTS: We obtained deep coverage over targeted regions and demonstrated accurate variant detection by comparison to whole-genome sequencing in a control individual. We found a dominant mutation GH1, p.R209H, in a three-generation pedigree with IGHD. CONCLUSIONS: smMIPS is an efficient and inexpensive method to detect mutations in patients with hypopituitarism, drastically limiting the need for screening individual genes by Sanger sequencing.
RESUMEN
Predicting function from sequence is an important goal in current biological research, and although, broad functional assignment is possible when a protein is assigned to a family, predicting functional specificity with accuracy is not straightforward. If function is provided by key structural properties and the relevant properties can be computed using the sequence as the starting point, it should in principle be possible to predict function in detail. The truncated hemoglobin family presents an interesting benchmark study due to their ubiquity, sequence diversity in the context of a conserved fold and the number of characterized members. Their functions are tightly related to O2 affinity and reactivity, as determined by the association and dissociation rate constants, both of which can be predicted and analyzed using in-silico based tools. In the present work we have applied a strategy, which combines homology modeling with molecular based energy calculations, to predict and analyze function of all known truncated hemoglobins in an evolutionary context. Our results show that truncated hemoglobins present conserved family features, but that its structure is flexible enough to allow the switch from high to low affinity in a few evolutionary steps. Most proteins display moderate to high oxygen affinities and multiple ligand migration paths, which, besides some minor trends, show heterogeneous distributions throughout the phylogenetic tree, again suggesting fast functional adaptation. Our data not only deepens our comprehension of the structural basis governing ligand affinity, but they also highlight some interesting functional evolutionary trends.
Asunto(s)
Hemoglobinas Truncadas , Secuencia de Aminoácidos , Biología Computacional , Evolución Molecular , Modelos Lineales , Modelos Moleculares , Datos de Secuencia Molecular , Oxígeno/metabolismo , Filogenia , Alineación de Secuencia , Hemoglobinas Truncadas/química , Hemoglobinas Truncadas/genética , Hemoglobinas Truncadas/fisiologíaRESUMEN
BACKGROUND: Understanding the molecular mechanism through which proteins are functional at extreme high and low temperatures is one of the key issues in structural biology. To investigate this phenomenon, we have focused on two instructive truncated hemoglobins from Thermobifida fusca (Tf-trHbO) and Mycobacterium tuberculosis (Mt-trHbO); although the two proteins are structurally nearly identical, only the former is stable at high temperatures. METHODS: We used molecular dynamics simulations at different temperatures as well as thermal melting profile measurements of both wild type proteins and two mutants designed to interchange the amino acid residue, either Pro or Gly, at E3 position. RESULTS: The results show that the presence of a Pro at the E3 position is able to increase (by 8°) or decrease (by 4°) the melting temperature of Mt-trHbO and Tf-trHbO, respectively. We observed that the ProE3 alters the structure of the CD loop, making it more flexible. CONCLUSIONS: This gain in flexibility allows the protein to concentrate its fluctuations in this single loop and avoid unfolding. The alternate conformations of the CD loop also favor the formation of more salt-bridge interactions, together augmenting the protein's thermostability. GENERAL SIGNIFICANCE: These results indicate a clear structural and dynamical role of a key residue for thermal stability in truncated hemoglobins.
Asunto(s)
Modelos Moleculares , Mycobacterium tuberculosis/metabolismo , Estabilidad Proteica , Hemoglobinas Truncadas/química , Actinomycetales/química , Actinomycetales/metabolismo , Calor , Humanos , Simulación de Dinámica Molecular , Mycobacterium tuberculosis/química , Hemoglobinas Truncadas/aislamiento & purificación , Hemoglobinas Truncadas/metabolismoRESUMEN
Internal water molecules play an active role in ligand uptake regulation, since displacement of retained water molecules from protein surfaces or cavities by incoming ligands can promote favorable or disfavorable effects over the global binding process. Detection of these water molecules by X-ray crystallography is difficult given their positional disorder and low occupancy. In this work, we employ a combination of molecular dynamics simulations and ligand rebinding over a broad time range to shed light into the role of water molecules in ligand migration and binding. Computational studies on the unliganded structure of the thermostable truncated hemoglobin from Thermobifida fusca (Tf-trHbO) show that a water molecule is in the vicinity of the iron heme, stabilized by WG8 with the assistance of YCD1, exerting a steric hindrance for binding of an exogenous ligand. Mutation of WG8 to F results in a significantly lower stabilization of this water molecule and in subtle dynamical structural changes that favor ligand binding, as observed experimentally. Water is absent from the fully hydrophobic distal cavity of the triple mutant YB10F-YCD1F-WG8F (3F), due to the lack of residues capable of stabilizing it nearby the heme. In agreement with these effects on the barriers for ligand rebinding, over 97% of the photodissociated ligands are rebound within a few nanoseconds in the 3F mutant case. Our results demonstrate the specific involvement of water molecules in shaping the energetic barriers for ligand migration and binding.