Environment-specific amino acid substitution tables: tertiary templates and prediction of protein folds.

The local environment of an amino acid in a folded protein determines the acceptability of mutations at that position. In order to characterize and quantify these structural constraints, we have made a comparative analysis of families of homologous proteins. Residues in each structure are classified according to amino acid type, secondary structure, accessibility of the side chain, and existence of hydrogen bonds from the side chains. Analysis of the pattern of observed substitutions as a function of local environment shows that there are distinct patterns, especially for buried polar residues. The substitution data tables are available on diskette with Protein Science. Given the fold of a protein, one is able to predict sequences compatible with the fold (profiles or templates) and potentially to discriminate between a correctly folded and misfolded protein. Conversely, analysis of residue variation across a family of aligned sequences in terms of substitution profiles can allow prediction of secondary structure or tertiary environment.     

Direct observation by X-ray analysis of the tetrahedral "intermediate" of aspartic proteinases.

We report the X-ray analysis at 2.0 A resolution for crystals of the aspartic proteinase endothiapepsin (EC 3.4.23.6) complexed with a potent difluorostatone-containing tripeptide renin inhibitor (CP-81,282). The scissile bond surrogate, an electrophilic ketone, is hydrated in the complex. The pro-(R) (statine-like) hydroxyl of the tetrahedral carbonyl hydrate is hydrogen-bonded to both active-site aspartates 32 and 215 in the position occupied by a water in the native enzyme. The second hydroxyl oxygen of the hydrate is hydrogen-bonded only to the outer oxygen of Asp 32. These experimental data provide a basis for a model of the tetrahedral intermediate in aspartic proteinase-mediated cleavage of the amide bond. This indicates a mechanism in which Asp 32 is the proton donor and Asp 215 carboxylate polarizes a bound water for nucleophilic attack. The mechanism involves a carboxylate (Asp 32) that is stabilized by extensive hydrogen bonding, rather than an oxyanion derivative of the peptide as in serine proteinase catalysis.     

Molecular and cellular remodeling of HepG2 cells upon treatment with antitubercular drugs

Drug-induced liver injury (DILI) is an adverse outcome of the currently used tuberculosis treatment regimen, which results in patient noncompliance, poor treatment outcomes, and the emergence of drug-resistant tuberculosis. DILI is primarily caused by the toxicity of the drugs and their metabolites, which affect liver cells, biliary epithelial cells, and liver vasculature. However, the precise mechanism behind the cellular damage attributable to first-line antitubercular drugs (ATDs), as well as the effect of toxicity on the cell survival strategies, is yet to be elucidated. In the current study, HepG2 cells upon treatment with a high concentration of ATDs showed increased perforation within the cell, cuboidal shape, and membrane blebbing as compared with control/untreated cells. It was observed that ATD-induced toxicity in HepG2 cells leads to altered mitochondrial membrane permeability, which was depicted by the decreased fluorescence intensity of the MitoRed tracker dye at higher drug concentrations. In addition, high doses of ATDs caused cell damage through an increase in reactive oxygen species production in HepG2 cells and a simultaneous reduction in glutathione levels. Further, high dose of isoniazid (50–200 mM), pyrazinamide (50–200 mM), and rifampicin (20–100 µM) causes cell apoptosis and affects cell survival during toxic conditions by decreasing the expression of potent autophagy markers Atg5, Atg7, and LC3B. Thus, ATD-mediated toxicity contributes to the reduced ability of hepatocytes to tolerate cellular damage caused by altered mitochondrial membrane permeability, increased apoptosis, and decreased autophagy. These findings further emphasize the need to develop adjuvant therapies that can mitigate ATD-induced toxicity for the effective treatment of tuberculosis.     

Integrative structure and function of the yeast exocyst complex

Exocyst is an evolutionarily conserved hetero‐octameric tethering complex that plays a variety of roles in membrane trafficking, including exocytosis, endocytosis, autophagy, cell polarization, cytokinesis, pathogen invasion, and metastasis. Exocyst serves as a platform for interactions between the Rab, Rho, and Ral small GTPases, SNARE proteins, and Sec1/Munc18 regulators that coordinate spatial and temporal fidelity of membrane fusion. However, its mechanism is poorly described at the molecular level. Here, we determine the molecular architecture of the yeast exocyst complex by an integrative approach, based on a 3D density map from negative‐stain electron microscopy (EM) at ~16 Å resolution, 434 disuccinimidyl suberate and 1‐ethyl‐3‐(3‐dimethylaminopropyl)carbodiimide hydrochloride cross‐links from chemical‐crosslinking mass spectrometry, and partial atomic models of the eight subunits. The integrative structure is validated by a previously determined cryo‐EM structure, cross‐links, and distances from in vivo fluorescence microscopy. Our subunit configuration is consistent with the cryo‐EM structure, except for Sec5. While not observed in the cryo‐EM map, the integrative model localizes the N‐terminal half of Sec3 near the Sec6 subunit. Limited proteolysis experiments suggest that the conformation of Exo70 is dynamic, which may have functional implications for SNARE and membrane interactions. This study illustrates how integrative modeling based on varied low‐resolution structural data can inform biologically relevant hypotheses, even in the absence of high‐resolution data.     

Using Integrative Modeling Platform to compute,validate, and archive a model of a protein complex structure

Biology is advanced by producing structural models of biological systems, such as protein complexes. Some systems are recalcitrant to traditional structure determination methods. In such cases, it may still be possible to produce useful models by integrative structure determination that depends on simultaneous use of multiple types of data. An ensemble of models that are sufficiently consistent with the data is produced by a structural sampling method guided by a data‐dependent scoring function. The variation in the ensemble of models quantified the uncertainty of the structure, generally resulting from the uncertainty in the input information and actual structural heterogeneity in the samples used to produce the data. Here, we describe how to generate, assess, and interpret ensembles of integrative structural models using our open source Integrative Modeling Platform program ( https://integrativemodeling.org ).     

Genome‐wide association study identifies loci associated with liability to alcohol and drug dependence that is associated with variability in reward‐related ventral striatum activity in African‐ and European‐Americans

Genetic influences on alcohol and drug dependence partially overlap, however, specific loci underlying this overlap remain unclear. We conducted a genome‐wide association study (GWAS) of a phenotype representing alcohol or illicit drug dependence (ANYDEP) among 7291 European‐Americans (EA; 2927 cases) and 3132 African‐Americans (AA: 1315 cases) participating in the family‐based Collaborative Study on the Genetics of Alcoholism. ANYDEP was heritable (h ² in EA = 0.60, AA = 0.37). The AA GWAS identified three regions with genome‐wide significant (GWS; P < 5E‐08) single nucleotide polymorphisms (SNPs) on chromosomes 3 (rs34066662, rs58801820) and 13 (rs75168521, rs78886294), and an insertion‐deletion on chromosome 5 (chr5:141988181). No polymorphisms reached GWS in the EA. One GWS region (chromosome 1: rs1890881) emerged from a trans‐ancestral meta‐analysis (EA + AA) of ANYDEP, and was attributable to alcohol dependence in both samples. Four genes (AA: CRKL, DZIP3, SBK3; EA: P2RX6) and four sets of genes were significantly enriched within biological pathways for hemostasis and signal transduction. GWS signals did not replicate in two independent samples but there was weak evidence for association between rs1890881 and alcohol intake in the UK Biobank. Among 118 AA and 481 EA individuals from the Duke Neurogenetics Study, rs75168521 and rs1890881 genotypes were associated with variability in reward‐related ventral striatum activation. This study identified novel loci for substance dependence and provides preliminary evidence that these variants are also associated with individual differences in neural reward reactivity. Gene discovery efforts in non‐European samples with distinct patterns of substance use may lead to the identification of novel ancestry‐specific genetic markers of risk.     

LigBase: a database of families of aligned ligand binding sites in known protein sequences and structures

A database comprising all ligand-binding sites of known structure aligned with all related protein sequences and structures is described. Currently, the database contains approximately 50000 ligand-binding sites for small molecules found in the Protein Data Bank (PDB). The structure-structure alignments are obtained by the Combinatorial Extension (CE) program (Shindyalov and Bourne, Protein Eng., 11, 739-747, 1998) and sequence-structure alignments are extracted from the ModBase database of comparative protein structure models for all known protein sequences (Sanchez et al., Nucleic Acids Res., 28, 250-253, 2000). It is possible to search for binding sites in LigBase by a variety of criteria. LigBase reports summarize ligand data including relevant structural information from the PDB file, such as ligand type and size, and contain links to all related protein sequences in the TrEMBL database. Residues in the binding sites are graphically depicted for comparison with other structurally defined family members. LigBase provides a resource for the analysis of families of related binding sites.