Modeling of loops in protein structures

Comparative protein structure prediction is limited mostly by the errors in alignment and loop modeling. We describe here a new automated modeling technique that significantly improves the accuracy of loop predictions in protein structures. The positions of all nonhydrogen atoms of the loop are optimized in a fixed environment with respect to a pseudo energy function. The energy is a sum of many spatial restraints that include the bond length, bond angle, and improper dihedral angle terms from the CHARMM-22 force field, statistical preferences for the main-chain and side-chain dihedral angles, and statistical preferences for nonbonded atomic contacts that depend on the two atom types, their distance through space, and separation in sequence. The energy function is optimized with the method of conjugate gradients combined with molecular dynamics and simulated annealing. Typically, the predicted loop conformation corresponds to the lowest energy conformation among 500 independent optimizations. Predictions were made for 40 loops of known structure at each length from 1 to 14 residues. The accuracy of loop predictions is evaluated as a function of thoroughness of conformational sampling, loop length, and structural properties of native loops. When accuracy is measured by local superposition of the model on the native loop, 100, 90, and 30% of 4-, 8-, and 12-residue loop predictions, respectively, had <2 A RMSD error for the mainchain N, C(alpha), C, and O atoms; the average accuracies were 0.59 +/- 0.05, 1.16 +/- 0.10, and 2.61 +/- 0.16 A, respectively. To simulate real comparative modeling problems, the method was also evaluated by predicting loops of known structure in only approximately correct environments with errors typical of comparative modeling without misalignment. When the RMSD distortion of the main-chain stem atoms is 2.5 A, the average loop prediction error increased by 180, 25, and 3% for 4-, 8-, and 12-residue loops, respectively. The accuracy of the lowest energy prediction for a given loop can be estimated from the structural variability among a number of low energy predictions. The relative value of the present method is gauged by (1) comparing it with one of the most successful previously described methods, and (2) describing its accuracy in recent blind predictions of protein structure. Finally, it is shown that the average accuracy of prediction is limited primarily by the accuracy of the energy function rather than by the extent of conformational sampling.     

Prediction of the secondary structures of stefins and cystatins, the low-molecular mass protein inhibitors of cysteine proteinases

A procedure for classifying proteins of known sequence into structurally similar groups was developed on the basis of the Argos parametric approach. It is shown that stefins and cystatins constitute two structurally well resolved, but homologous groups of proteins. Furthermore, it is very probable that segments of secondary structures within each family are conserved, although significant differences between stefins and cystatins are indicated at the level of secondary structure. Next, secondary structures of all sequenced stefins and cystatins were predicted and used in the construction of secondary structures of the "typical stefin" and the "typical cystatin". Results were interpreted in the light of evolution and inhibition mechanism: Alignment of the "typical stefin" versus the "typical cystatin" secondary structure segments suggests that the divergence of stefin and cystatin families did not occur by a gene fusion event, but only by a mechanism of substitution, insertion and/or deletion. The central region of low-molecular mass cystatins, which is assumed to interact with cysteine proteinases, is predicted to be in a beta-sheet conformation. This resembles the beta-sheet in the active site of "standard mechanism" serine proteinases inhibitors.     

The collaborative study on the genetics of alcoholism: Genetics

This review describes the genetic approaches and results from the family-based Collaborative Study on the Genetics of Alcoholism (COGA). COGA was designed during the linkage era to identify genes affecting the risk for alcohol use disorder (AUD) and related problems, and was among the first AUD-focused studies to subsequently adopt a genome-wide association (GWAS) approach. COGA's family-based structure, multimodal assessment with gold-standard clinical and neurophysiological data, and the availability of prospective longitudinal phenotyping continues to provide insights into the etiology of AUD and related disorders. These include investigations of genetic risk and trajectories of substance use and use disorders, phenome-wide association studies of loci of interest, and investigations of pleiotropy, social genomics, genetic nurture, and within-family comparisons. COGA is one of the few AUD genetics projects that includes a substantial number of participants of African ancestry. The sharing of data and biospecimens has been a cornerstone of the COGA project, and COGA is a key contributor to large-scale GWAS consortia. COGA's wealth of publicly available genetic and extensive phenotyping data continues to provide a unique and adaptable resource for our understanding of the genetic etiology of AUD and related traits.     

Tertiary structural constraints on protein evolutionary diversity: templates, key residues and structure prediction

The pattern of residue substitution in divergently evolving families of globular proteins is highly variable. At each position in a fold there are constraints on the identities of amino acids from both the three-dimensional structure and the function of the protein. To characterize and quantify the structural constraints, we have made a comparative analysis of families of homologous globular proteins. Residues are classified according to amino acid type, secondary structure, accessibility of the sidechain, and existence of hydrogen bonds from sidechain to other sidechains or peptide carbonyl or amide functions. There are distinct patterns of substitution especially where residues are both solvent inaccessible and hydrogen bonded through their sidechains. The patterns of residue substitution can be used to construct templates or to identify 'key' residues if one or more structures are known. Conversely, analysis of conversation and substitution across a large family of aligned sequences in terms of substitution profiles can allow prediction of tertiary environment or indicate a functional role. Similar analyses can be used to test the validity of putative structures if several homologous sequences are available.     

The effect of a single injection of irinotecan on the development of enamel in the Wistar rats

Cancer is the second most frequent cause of death in children. Because the prognosis for childhood malignancies has improved, attention has now focused on long‐term consequences of cancer treatment. The immediate effects of chemotherapy on soft tissues have been well described; however, there is less information about long‐term effects of chemotherapy on the development of dental tissues. To test the association between the effect of chemotherapy on enamel development, we examined two groups of rats: one that had received an intraperitoneal dose of 200 mg/kg of irinotecan, whereas the other (control) group had received vehicle only. Rats were killed at 6, 48 and 96 hr post‐injection; the mandibles dissected out, fixed for histological evaluation and scanned for mineralization defects by Micro‐CT. Our results showed structural changes in the ameloblast layer along with a significant reduction in mineralization and thickness of enamel at 96 hr after chemotherapy. These data demonstrate that irinotecan induces structural changes in forming enamel that become apparent after anticancer chemotherapy treatment.     

Vanilloid receptor-related osmotically activated channel (VR-OAC), a candidate vertebrate osmoreceptor

The detection of osmotic stimuli is essential for all organisms, yet few osmoreceptive proteins are known, none of them in vertebrates. By employing a candidate-gene approach based on genes encoding members of the TRP superfamily of ion channels, we cloned cDNAs encoding the vanilloid receptor-related osmotically activated channel (VR-OAC) from the rat, mouse, human, and chicken. This novel cation-selective channel is gated by exposure to hypotonicity within the physiological range. In the central nervous system, the channel is expressed in neurons of the circumventricular organs, neurosensory cells responsive to systemic osmotic pressure. The channel also occurs in other neurosensory cells, including inner-ear hair cells, sensory neurons, and Merkel cells.     

ModBase: a database of comparative protein structure models

ModBase is a database of evaluated and annotated comparative protein structure models. The database also includes fold assignments and alignments on which the models were based.     

The structure of chagasin in complex with a cysteine protease clarifies the binding mode and evolution of an inhibitor family

Protein inhibitors of proteolytic enzymes regulate proteolysis and prevent the pathological effects of excess endogenous or exogenous proteases. Cysteine proteases are a large family of enzymes found throughout the plant and animal kingdoms. Disturbance of the equilibrium between cysteine proteases and natural inhibitors is a key event in the pathogenesis of cancer, rheumatoid arthritis, osteoporosis, and emphysema. A family (I42) of cysteine protease inhibitors (http://merops.sanger.ac.uk) was discovered in protozoan parasites and recently found widely distributed in prokaryotes and eukaryotes. We report the 2.2 A crystal structure of the signature member of the I42 family, chagasin, in complex with a cysteine protease. Chagasin has a unique variant of the immunoglobulin fold with homology to human CD8alpha. Interactions of chagasin with a target protease are reminiscent of the cystatin family inhibitors. Protein inhibitors of cysteine proteases may have evolved more than once on nonhomologous scaffolds.     

Immunoglobulin motif DNA recognition and heterodimerization of the PEBP2/CBF Runt domain.

The polyomavirus enhancer binding protein 2 (PEBP2) or core binding factor (CBF) is a heterodimeric enhancer binding protein that is associated with genetic regulation of hematopoiesis and osteogenesis. Aberrant forms of PEBP2/CBF are implicated in the cause of the acute human leukemias and in a disorder of bone development known as cleidocranial dysplasia. The common denominator in the natural and mutant forms of this protein is a highly conserved domain of PEBP2/CBF alpha, termed the Runt domain (RD), which is responsible for both DNA binding and heterodimerization with the beta subunit of PEBP2/CBF. The three-dimensional structure of the RD bound to DNA has been determined to be an S-type immunoglobulin fold, establishing a structural relationship between the RD and the core DNA binding domains of NF-kappaB, NFAT1, p53 and the STAT proteins. NMR spectroscopy of a 43.6 kD RD-beta-DNA ternary complex identified the surface of the RD in contact with the beta subunit, suggesting a mechanism for the enhancement of RD DNA binding by beta. Analysis of leukemogenic mutants within the RD provides molecular insights into the role of this factor in leukemogenesis and cleidocranial dysplasia.