Efficient methods for filtering and ranking fragments for the prediction of structurally variable regions in proteins

The prediction of protein 3D structures close to insertions and deletions or, more generally, loop prediction, is still one of the major challenges in homology modeling projects. In this article, we developed ranking criteria and selection filters to improve knowledge-based loop predictions. These criteria were developed and optimized for a test data set containing 678 insertions and deletions. The examples are, in principle, predictable from the used loop database with an RMSD < 1 A and represent realistic modeling situations. Four noncorrelated criteria for the selection of fragments are evaluated. A fast prefilter compares the distance between the anchor groups in the template protein with the stems of the fragments. The RMSD of the anchor groups is used for fitting and ranking of the selected loop candidates. After fitting, repulsive close contacts of loop candidates with the template protein are used for filtering, and fragments with backbone torsion angles, which are unfavorable according to a knowledge-based potential, are eliminated. By the combined application of these filter criteria to the test set, it was possible to increase the percentage of predictions with a global RMSD < 1 A to over 50% among the first five ranks, with average global RMSD values for the first rank candidate that are between 1.3 and 2.2 A for different loop lengths. Compared to other examples described in the literature, our large numbers of test cases are not self-predictions, where loops are placed in a protein after a peptide loop has been cut out, but are attempts to predict structural changes that occur in evolution when a protein is affected by insertions and deletions.     

Positioning of anchor groups in protein loop prediction: the importance of solvent accessibility and secondary structure elements

The prediction of loop regions in the process of protein structure prediction by homology is still an unsolved problem. In an earlier publication, we could show that the correct placement of the amino acids serving as an anchor group to be connected by a loop fragment with a predicted geometry is a highly important step and an essential requirement within the process (Lessel and Schomburg, Proteins 1999; 37:56-64). In this article, we present an analysis of the quality of possible loop predictions with respect to gap length, fragment length, amino acid type, secondary structure, and solvent accessibility. For 550 insertions and 544 deletions, we test all possible positions for anchor groups with an inserted loop of a length between 3 and 12 amino acids. We could show that approximately 80% of the indel regions could be predicted within 1.5 A RMSD from a knowledge-based loop data base if criteria for the correct localization of anchor groups could be found and the loops can be sorted correctly. From our analysis, several conclusions regarding the optimal placement of anchor groups become obvious: (1) The correct placement of anchor groups is even more important for longer gap lengths, (2) medium length fragments (length 5-8) perform better than short or long ones, (3) the placement of anchor groups at hydrophobic amino acids gives a higher chance to include the best possible loop, (4) anchor groups within secondary structure elements, in particular beta-sheets are suitable, (5) amino acids with lower solvent accessibility are better anchor group. A preliminary test using a combination of the anchor group positioning criteria deduced from our analysis shows very promising results.     

Inflammatory Bowel Disease and Risk of Adverse Pregnancy Outcomes

Background and Objectives

Existing data on pregnancy complications in inflammatory bowel disease (IBD) are inconsistent. To address these inconsistencies, we investigated potential associations between IBD, IBD-related medication use during pregnancy, and pregnancy loss, pre-eclampsia, preterm delivery, Apgar score, and congenital abnormalities.

Methods

We conducted a cohort study in >85,000 Danish National Birth Cohort women who were pregnant in the period 1996-2002 and had information on IBD, IBD-related medication use (systemic or local corticosteroids, 5-aminosalicylates), pregnancy outcomes and potential confounders. We evaluated associations between IBD and adverse pregnancy/birth outcomes using Cox regression and log-linear binomial regression.

Results

IBD was strongly and significantly associated with severe pre-eclampsia, preterm premature rupture of membranes and medically indicated preterm delivery in women using systemic corticosteroids during pregnancy (hazard ratios [HRs] >7). IBD was also associated with premature preterm rupture of membranes in women using local corticosteroid medications (HR 3.30, 95% confidence interval [CI] 1.33-8.20) and with medically indicated preterm delivery (HR 1.91, 95% CI 0.99-3.68) in non-medicated women. Furthermore, IBD was associated with low 5-minute Apgar score in term infants (risk ratio [RR] 2.19, 95% CI 1.03-4.66). Finally, Crohn’s disease (but not ulcerative colitis) was associated with major congenital abnormalities in the offspring (RR 1.85, 95% CI 1.06-3.21). No child with a congenital abnormality born to a woman with IBD was exposed to systemic corticosteroids in utero.

Conclusion

Women with IBD are at increased risk of severe pre-eclampsia, medically indicated preterm delivery, preterm premature rupture of membranes, and delivering infants with low Apgar score and major congenital malformations. These associations are only partly explained by severe disease as reflected by systemic corticosteroid use.     

The MHC2TA -168A>G gene polymorphism is not associated with rheumatoid arthritis in Austrian patients

An association between susceptibility to rheumatoid arthritis (RA) and a common -168A>G polymorphism in the MHC2TA gene with differential major histocompatibility complex (MHC) II molecule expression was recently reported in a Swedish population. The objective of the present study was to replicate this finding by examining the -168A>G polymorphism in an Austrian case-control study. Three hundred and sixty-two unrelated RA cases and 351 sex-matched and age-matched controls as well as 1,709 Austrian healthy individuals were genotyped. All participants were from the same ethnic background. Genotyping was performed using 5' allelic discrimination assays. The association between susceptibility to RA and the -168A>G single nucleotide polymorphism was examined by chi-square test. Comparison was made assuming a dominant effect (AG + GG genotypes versus AA genotype). In contrast to the primary report, the frequency of MHC2TA -168G allele carriers was not significantly different between patients and controls in the Austrian cohort. The homozygous MHC2TA -168 GG genotype was more frequent in matched controls than in Austrian RA patients. There was no association between the presence of RA-specific autoantibodies and the MHC2TA -168 GG genotype. In this cohort of Austrian patients, no association between the MHC2TA polymorphism and RA was found.     

The Sec1/Munc18 Protein Groove Plays a Conserved Role in Interaction with Sec9p/SNAP‐25

The Sec1/Munc18 (SM) proteins constitute a conserved family with essential functions in SNARE‐mediated membrane fusion. Recently, a new protein–protein interaction site in Sec1p, designated the groove, was proposed. Here, we show that a sec1 groove mutant yeast strain, sec1(w24), displays temperature‐sensitive growth and secretion defects. The yeast Sec1p and mammalian Munc18‐1 grooves were shown to play an important role in the interaction with the SNAREs Sec9p and SNAP‐25b, respectively. Incubation of SNAP‐25b with the Munc18‐1 groove mutant resulted in a lag in the kinetics of SNARE complex assembly in vitro when compared with wild‐type Munc18‐1. The SNARE regulator SRO7 was identified as a multicopy suppressor of sec1(w24) groove mutant and an intact Sec1p groove was required for the plasma membrane targeting of Sro7p–SNARE complexes. Simultaneous inactivation of Sec1p groove and SRO7 resulted in reduced levels of exocytic SNARE complexes. Our results identify the groove as a conserved interaction surface in SM proteins. The results indicate that this structural element is important for interactions with Sec9p/SNAP‐25 and participates, in concert with Sro7p, in the initial steps of SNARE complex assembly.      

Field‐based comparison of ligand and coactivator binding sites of nuclear receptors

A structure‐based comparison of the ligand‐binding domains of 35 nuclear receptors from five different subfamilies is presented. Their ligand and coactivator binding sites are characterized using knowledge‐based contact preference fields for hydrophobic and hydrophilic interactions implemented in the MOE modeling environment. Additionally, for polar knowledge‐based field points the preference for negative or positive electrostatic interactions is estimated using the Poisson‐Boltzmann equation. These molecular‐interaction fields are used to cluster the nuclear receptor family based on similarities of their binding sites. By analyzing the similarities and differences of hydrophobic and polar fields in binding pockets of related receptors it is possible to identify conserved interactions in ligand and coactivator binding pockets, which support e.g. design of specific ligands during lead optimization or virtual screening as docking filter. Examples of remarkable similarities between ligand binding sites of members from phylogenetically different nuclear receptor families (RXR, RAR, HNF4, NR5) and differences between closely related subtypes (LXR, RAR, TR) are discussed in more detail. Significant similarities and differences of coactivator binding sites are shown for NR3Cs, LXRs and PPARs. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 884–894, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Population based study of rates of multiple pregnancies in Denmark, 1980-94.

OBJECTIVE: To study trends in multiple pregnancies not explained by changes in maternal age and parity patterns. DESIGN: Trends in population based figures for multiple pregnancies in Denmark studied from complete national records on parity history and vital status. POPULATION: 497,979 Danish women and 803,019 pregnancies, 1980-94. MAIN OUTCOME MEASURES: National rates of multiple pregnancies, infant mortality, and stillbirths controlled for maternal age and parity. Special emphasis on primiparous women > or = 30 years of age, who are most likely to undergo fertility treatment. RESULTS: The national incidence of multiple pregnancies increased 1.7-fold during 1980-94, the increase primarily in 1989-94 and almost exclusively in primiparous women aged > or = 30 years, for whom the adjusted population based twinning rate increased 2.7-fold and the triplet rate 9.1-fold. During 1989-94, the adjusted yearly increase in multiple pregnancies for these women was 19% (95% confidence interval 16% to 21%) and in dizygotic twin pregnancies 25% (21% to 28%). The proportion of multiple births among infant deaths in primiparous women > or = 30 years increased from 11.5% to 26.9% during the study period. The total infant mortality, however, did not increase for these women because of a simultaneous significant decrease in infant mortality among singletons. CONCLUSIONS: A relatively small group of women has drastically changed the overall national rates of multiple pregnancies. The introduction of new treatments to enhance fertility has probably caused these changes and has also affected the otherwise decreasing trend in infant mortality. Consequently, the resources, both economical and otherwise, associated with these treatments go well beyond those invested in specific fertility enhancing treatments.     

Environmental risk factors differ between rheumatoid arthritis with and without auto-antibodies against cyclic citrullinated peptides

The aim of this study was to evaluate new and previously hypothesised non-genetic risk factors for serologic subtypes of rheumatoid arthritis (RA) defined by the presence or absence of auto-antibodies to cyclic citrullinated peptides (CCP). In a national case-control study, we included 515 patients recently diagnosed with RA according to the American College of Rheumatology 1987 classification criteria and 769 gender- and age-matched population controls. Telephone interviews provided information about non-genetic exposures, and serum samples for patients were tested for anti-CCP-antibodies. Associations between exposure variables and risk of anti-CCP-positive and anti-CCP-negative RA were evaluated using logistic regression. A series of RA subtype-specific risk factors were identified. Tobacco smoking (odds ratio [OR] = 1.65; 95% confidence interval: 1.03–2.64, for >20 versus 0 pack-years) was selectively associated with risk of anti-CCP-positive RA, whereas alcohol consumption exhibited an inverse dose-response association with this RA subtype (OR = 1.98, 1.22–3.19, for 0 versus >0–5 drinks per week). Furthermore, coffee consumption (OR = 2.18; 1.07–4.42, for >10 versus 0 cups per day), ever use of oral contraceptives (OR = 1.65; 1.06–2.57) and having a first-degree relative with schizophrenia (OR = 4.18; 1.54–11.3) appeared more strongly associated with risk of anti-CCP-positive RA. Obesity was selectively associated with risk of anti-CCP-negative RA, with obese individuals being at more than 3-fold increased risk of this subtype compared with normal-weight individuals (OR = 3.45; 1.73–6.87). Age at menarche was the only examined factor that was significantly associated with both serologic subtypes of RA (p-trends = 0.01); women with menarche at age ≥ 15 years had about twice the risk of either RA subtype compared with women with menarche at age ≤ 12 years. Major differences in risk factor profiles suggest distinct etiologies for anti-CCP-positive and anti-CCP-negative RA.