Mutations in STIL, Encoding a Pericentriolar and Centrosomal Protein, Cause Primary Microcephaly

Primary microcephaly (MCPH) is an autosomal-recessive congenital disorder characterized by smaller-than-normal brain size and mental retardation. MCPH is genetically heterogeneous with six known loci: MCPH1-MCPH6. We report mapping of a novel locus, MCPH7, to chromosome 1p32.3-p33 between markers D1S2797 and D1S417, corresponding to a physical distance of 8.39 Mb. Heterogeneity analysis of 24 families previously excluded from linkage to the six known MCPH loci suggested linkage of five families (20.83%) to the MCPH7 locus. In addition, four families were excluded from linkage to the MCPH7 locus as well as all of the six previously known loci, whereas the remaining 15 families could not be conclusively excluded or included. The combined maximum two-point LOD score for the linked families was 5.96 at marker D1S386 at theta = 0.0. The combined multipoint LOD score was 6.97 between markers D1S2797 and D1S417. Previously, mutations in four genes, MCPH1, CDK5RAP2, ASPM, and CENPJ, that code for centrosomal proteins have been shown to cause this disorder. Three different homozygous mutations in STIL, which codes for a pericentriolar and centrosomal protein, were identified in patients from three of the five families linked to the MCPH7 locus; all are predicted to truncate the STIL protein. Further, another recently ascertained family was homozygous for the same mutation as one of the original families. There was no evidence for a common haplotype. These results suggest that the centrosome and its associated structures are important in the control of neurogenesis in the developing human brain.

Mutations in WDR62, encoding a centrosomal and nuclear protein, in Indian primary microcephaly families with cortical malformations

Primary microcephaly is an autosomal recessive disorder characterized by smaller than normal brain size and mental retardation. It is genetically heterogeneous with seven loci: MCPH1-MCPH7. We have previously reported genetic analysis of 35 families, including the identification of the MCPH7 gene STIL. Of the 35 families, three families showed linkage to the MCPH2 locus. Recent whole-exome sequencing studies have shown that the WDR62 gene, located in the MCPH2 candidate region, is mutated in patients with severe brain malformations. We therefore sequenced the WDR62 gene in our MCPH2 families and identified two novel homozygous protein truncating mutations in two families. Affected individuals in the two families had pachygyria, microlissencephaly, band heterotopias, gyral thickening, and dysplastic cortex. Using immunofluorescence study, we showed that, as with other MCPH proteins, WDR62 localizes to centrosomes in A549, HepG2, and HaCaT cells. In addition, WDR62 was also localized to nucleoli. Bioinformatics analysis predicted two overlapping nuclear localization signals and multiple WD-40 repeats in WDR62. Two other groups have also recently identified WDR62 mutations in MCPH2 families. Our results therefore add further evidence that WDR62 is the MCPH2 gene. The present findings will be helpful in genetic diagnosis of patients linked to the MCPH2 locus.

Mutations in hpyAVIBM, C5 cytosine DNA methyltransferase from Helicobacter pylori result in relaxed specificity

The genome of Helicobacter pylori is rich in restrictionmodification (RM) systems. Approximately 4% of the genome codes for components of RM systems. hpyAVIBM, which codes for a phase-variable C5 cytosine methyltransferase (MTase) from H. pylori, lacks a cognate restriction enzyme. Over-expression of M.HpyAVIB in Escherichia coli enhances the rate of mutations. However, when the catalytically inactive F9N or C82W mutants of M.HpyAVIB were expressed in E. coli, mutations were not observed. The M.HpyAVIB gene itself was mutated to give rise to different variants of the MTase. M.HpyAVIB variants were purified and differences in kinetic properties and specificity were observed. Intriguingly, purified MTase variants showed relaxed substrate specificity. Homologues of hpyAVIBM homologues amplified and sequenced from different clinical isolates showed similar variations in sequence. Thus, hpyAVIBM presents an interesting example of allelic variations in H. pylori where changes in the nucleotide sequence result in proteins with new properties.

An analysis of substitution, deletion and insertion mutations in cancer genes

Cancer-associated mutations in cancer genes constitute a diverse set of mutations associated with the disease. To gain insight into features of the set, substitution, deletion and insertion mutations were analysed at the nucleotide level, from the COSMIC database. The most frequent substitutions were c -> t, g -> a, g -> t, and the most frequent codon changes were to termination codons. Deletions more than insertions, FS (frameshift) indels more than I-F (in-frame) ones, and single-nucleotide indels, were frequent. FS indels cause loss of significant fractions of proteins. The 5'-cut in FS deletions, and 5'-ligation in FS insertions, often occur between pairs of identical bases. Interestingly, the cut-site and 3'-ligation in insertions, and 3'-cut and join-pair in deletions, were each found to be the same significantly often (p < 0.001). It is suggested that these features aid the incorporation of indel mutations. Tumor suppressors undergo larger numbers of mutations, especially disruptive ones, over the entire protein length, to inactivate two alleles. Proto-oncogenes undergo fewer, less-disruptive mutations, in selected protein regions, to activate a single allele. Finally, catalogues, in ranked order, of genes mutated in each cancer, and cancers in which each gene is mutated, were created. The study highlights the nucleotide level preferences and disruptive nature of cancer mutations.

Rationalization and prediction of drug resistant mutations in targets for clinical anti-tubercular drugs

Resistance to therapy limits the effectiveness of drug treatment in many diseases. Drug resistance can be considered as a successful outcome of the bacterial struggle to survive in the hostile environment of a drug-exposed cell. An important mechanism by which bacteria acquire drug resistance is through mutations in the drug target. Drug resistant strains (multi-drug resistant and extensively drug resistant) of Mycobacterium tuberculosis are being identified at alarming rates, increasing the global burden of tuberculosis. An understanding of the nature of mutations in different drug targets and how they achieve resistance is therefore important. An objective of this study is to first decipher sequence as well as structural bases for the observed resistance in known drug resistant mutants and then to predict positions in each target that are more prone to acquiring drug resistant mutations. A curated database containing hundreds of mutations in the 38 drug targets of nine major clinical drugs, associated with resistance is studied here. Mutations have been classified into those that occur in the binding site itself, those that occur in residues interacting with the binding site and those that occur in outer zones. Structural models of the wild type and mutant forms of the target proteins have been analysed to seek explanations for reduction in drug binding. Stability analysis of an entire array of 19 mutations at each of the residues for each target has been computed using structural models. Conservation indices of individual residues, binding sites and whole proteins are computed based on sequence conservation analysis of the target proteins. The analyses lead to insights about which positions in the polypeptide chain have a higher propensity to acquire drug resistant mutations. Thus critical insights can be obtained about the effect of mutations on drug binding, in terms of which amino acid positions and therefore which interactions should not be heavily relied upon, which in turn can be translated into guidelines for modifying the existing drugs as well as for designing new drugs. The methodology can serve as a general framework to study drug resistant mutants in other micro-organisms as well.

CSSI-PRO: a method for secondary structure type editing, assignment and estimation in proteins using linear combination of backbone chemical shifts

Estimation of secondary structure in polypeptides is important for studying their structure, folding and dynamics. In NMR spectroscopy, such information is generally obtained after sequence specific resonance assignments are completed. We present here a new methodology for assignment of secondary structure type to spin systems in proteins directly from NMR spectra, without prior knowledge of resonance assignments. The methodology, named Combination of Shifts for Secondary Structure Identification in Proteins (CSSI-PRO), involves detection of specific linear combination of backbone H-1(alpha) and C-13' chemical shifts in a two-dimensional (2D) NMR experiment based on G-matrix Fourier transform (GFT) NMR spectroscopy. Such linear combinations of shifts facilitate editing of residues belonging to alpha-helical/beta-strand regions into distinct spectral regions nearly independent of the amino acid type, thereby allowing the estimation of overall secondary structure content of the protein. Comparison of the predicted secondary structure content with those estimated based on their respective 3D structures and/or the method of Chemical Shift Index for 237 proteins gives a correlation of more than 90% and an overall rmsd of 7.0%, which is comparable to other biophysical techniques used for structural characterization of proteins. Taken together, this methodology has a wide range of applications in NMR spectroscopy such as rapid protein structure determination, monitoring conformational changes in protein-folding/ligand-binding studies and automated resonance assignment.

Conformation of the LL and LD hairpin bends with internal hydrogen bonds in proteins and peptides

The conformation of three linked peptide units having an internal 4 → 1 type of hydrogen bond has been studied in detail, and the low energy conformations are listed. These conformations all lead to the reversal of the chain direction, and may therefore be called as “hairpin bends” or “U-bends”. Since this bend can occur at the end of two chains hydrogen-bonded in the antiparallel β-conformation, it is also known as the “β-bend”. Two types of conformation are possible when the residues at the second and third Cα atoms are both of type L (the LL bend), while only one type is possible for the LD and the DL bend. The LL bend can also accommodate the sequences LG, GL, GG (G = glycine), while the LD bend can accommodate the sequences LG, GD and GG. The conformations for the sequences DD and DL are exact inverses (or mirror images) of those for the sequences LL and LD, respectively, and have dihedral angles (phi2, ψ2), (phi3, ψ3) of the same magnitudes, but of opposite signs as those for the former types, which are listed, along with the characteristics (length, angle and energy) of the hydrogen bonds. A comparison of the theoretical predictions with experimental data (from X-ray diffraction and NMR studies) on proteins and peptides, show reasonably good agreement. However, a systematic trend is observable in the experimental data, slightly deviating from theory, which indicates that some deformations occur in the shapes of the peptide units forming the bend, differing from that of the standard planar peptide unit.

An explanation for the rare occurrence of cis peptide units in proteins and polypeptides

The rarity of occurrence of cis peptide units is only partially explained by the higher intrinsic energy of the cis over the trans form, which provides a probability of 0·01 for cis peptide units to occur. An additional factor is the conformational restriction imposed by the occurrence of a cis peptide unit in a chain of trans units. Taking a section of three peptide units having the sequences trans-trans-trans (ttt) and trans-cis-trans (tct), conformational energy calculations indicate that the latter can occur only to an extent of 0·1%, unless there occurs the sequence X-Pro, in which case it is of the order of 30%. This explains the extreme rarity of cis peptide units, in general; however, it follows that even with non-prolyl residues, cis peptide units are not forbidden, but can occur in some rare examples and should be looked for.

Orthogonal beta-beta motifs in proteins

A super-secondary structural motif comprising two orthogonally oriented beta-strands connected by short linking segments of <5 residues has been identified from a data set of 65 independent protein crystal structures. Of the 42 examples from 14 proteins, a vast majority have only a single residue as the linking element. Analysis of the conformational angles at the junction reveals that the recently described type VIII beta-turn occurs frequently at the connecting hinge, while the type II beta-turn is also fairly common.

Pi-turns in proteins and peptides: Classification, conformation, occurrence, hydration and sequence.

The i + 5-->i hydrogen bonded turn conformation (pi-turn) with the fifth residue adopting alpha L conformation is frequently found at the C-terminus of helices in proteins and hence is speculated to be a "helix termination signal." An analysis of the occurrence of i + 5-->i hydrogen bonded turn conformation at any general position in proteins (not specifically at the helix C-terminus), using coordinates of 228 protein crystal structures determined by X-ray crystallography to better than 2.5 A resolution is reported in this paper. Of 486 detected pi-turn conformations, 367 have the (i + 4)th residue in alpha L conformation, generally occurring at the C-terminus of alpha-helices, consistent with previous observations. However, a significant number (111) of pi-turn conformations occur with (i + 4)th residue in alpha R conformation also, generally occurring in alpha-helices as distortions either at the terminii or at the middle, a novel finding. These two sets of pi-turn conformations are referred to by the names pi alpha L and pi alpha R-turns, respectively, depending upon whether the (i + 4)th residue adopts alpha L or alpha R conformations. Four pi-turns, named pi alpha L'-turns, were noticed to be mirror images of pi alpha L-turns, and four more pi-turns, which have the (i + 4)th residue in beta conformation and denoted as pi beta-turns, occur as a part of hairpin bend connecting twisted beta-strands. Consecutive pi-turns occur, but only with pi alpha R-turns. The preference for amino acid residues is different in pi alpha L and pi alpha R-turns. However, both show a preference for Pro after the C-termini. Hydrophilic residues are preferred at positions i + 1, i + 2, and i + 3 of pi alpha L-turns, whereas positions i and i + 5 prefer hydrophobic residues. Residue i + 4 in pi alpha L-turns is mainly Gly and less often Asn. Although pi alpha R-turns generally occur as distortions in helices, their amino acid preference is different from that of helices. Poor helix formers, such as His, Tyr, and Asn, also were found to be preferred for pi alpha R-turns, whereas good helix former Ala is not preferred. pi-Turns in peptides provide a picture of the pi-turn at atomic resolution. Only nine peptide-based pi-turns are reported so far, and all of them belong to pi alpha L-turn type with an achiral residue in position i + 4. The results are of importance for structure prediction, modeling, and de novo design of proteins.

Amino acid interaction preferences in proteins

Understanding the key factors that influence the interaction preferences of amino acids in the folding of proteins have remained a challenge. Here we present a knowledge-based approach for determining the effective interactions between amino acids based on amino acid type, their secondary structure, and the contact based environment that they find themselves in the native state structure as measured by their number of neighbors. We find that the optimal information is approximately encoded in a 60 x 60 matrix describing the 20 types of amino acids in three distinct secondary structures (helix, beta strand, and loop). We carry out a clustering scheme to understand the similarity between these interactions and to elucidate a nonredundant set. We demonstrate that the inferred energy parameters can be used for assessing the fit of a given sequence into a putative native state structure.

Conformations of disulfide bridges in proteins

The conformational characteristics of disulfide bridges in proteins have been analyzed using a dataset of 22 protein structures, available at a resolution of 2.0 Å, containing a total of 72 disulfide crosslinks. The parameters used in the analysis include (φ, Ψ) values at Cys residues, bridge dihedral angles χss, χ1i, χ1j, χ2i and χ2j the distances Cαi-Cαj and Cβi-Cβj between the Cα and Cβ atoms of Cys(i) and Cys(j). Eight families of bridge conformations with three or more occurrences have been identified on the basis of these stereochemical parameters. The most populated family corresponds to the "left handed spiral" identified earlier by Richardson ((1981) Adv. Protein Chem. 34, 167–330). Disulfide bridging across antiparallel extended strands is observed in α-lytic protease, crambin, and β-trypsin and this structure is shown to be very similar to those obtained in small cystine peptides. Solvent accessible surface area calculations show that the overwhelming majority of disulfide bridges are inaccessible to solvent.

Peptide mimics for structural features in proteins. Crystal structures of three heptapeptide helices with a C-terminal 6-->1 hydrogen bond.

The crystal structure determination of three heptapeptides containing alpha-aminoisobutyryl (Aib) residues as a means of helix stabilization provides a high-resolution characterization of 6-->1 hydrogen-bonded conformations, reminiscent of helix-terminating structural features in proteins. The crystal parameters for the three peptides, Boc-Val-Aib-X-Aib-Ala-Aib-Y-OMe, where X and Y are Phe, Leu (I), Leu, Phe (II) and Leu, Leu (III) are: (I) space group P1, Z = 1, a = 9.903 A, b = 10.709 A, c = 11.969 A, alpha = 102.94 degrees, beta = 103.41 degrees, gamma = 92.72 degrees, R = 4.55%; (II) space group P21, Z = 2, a = 10.052 A, b = 17.653 A, c = 13.510 A, beta = 108.45 degrees, R = 4.49%; (III) space group P1, Z = 2 (two independent molecules IIIa and IIIb in the asymmetric unit), a = 10.833 A, b = 13.850 A, c = 16.928 A, alpha = 99.77 degrees, beta = 105.90 degrees, gamma = 90.64 degrees, R = 8.54%. In all cases the helices form 3(10)/alpha-helical (or 3(10)helical) structures, with helical columns formed by head-to-tail hydrogen bonding. The helices assemble in an all-parallel motif in crystals I and III and in an antiparallel motif in II. In the four crystallographically characterized molecules, I, II, IIIa and IIIb, Aib(6) adopts a left-handed helical (hL) conformation with positive phi, psi values, resulting in 6-->1 hydrogen-bond formation between Aib(2) CO and Leu(7)/Phe(7) NH groups. In addition a 4-->1 hydrogen bond is seen between Aib(3) CO and Aib(6) NH groups. This pattern of hydrogen bonding is often observed at the C-terminus of helices proteins, with the terminal pi-type turn being formed by four residues adopting the hRhRhRhL conformation.

Conformational characteristics of asparaginyl residues in proteins

Backbone conformations at 1064 asparaginyl residues in 123 non-homologous, high-resolution X-ray structures of proteins were analysed. Asn adopts conformations in left-handed x-helical region and other partially allowed regions in the Ramachandran map more readily than any other non-glycyl residue. Asn conformational clusters in the (phi,psi) regions of left-handed alpha-helix, right-handed alpha-helix and extended (beta) strands were investigated in detail for their occurrence in various secondary structures, especially in beta-turn regions. Preferences were observed for Asn conformations in different positions in various beta-turn types, including the first and fourth positions of the turn. Asparaginyl residues with extended conformations are found to occur frequently in irregular regions, although they are expected to occur predominantly in extended strands or in the third position of type II beta-turns. Asn conformations at the N-cap positions of helices strongly prefer extended conformation than alpha(L), which seems to be characteristic of non-glycyl residues at that position. In the linkers connecting two extended strands and those connecting an alpha-helix and an extended strand, Asn with alpha(L) or alpha(R) conformation is more favoured than Asn with the beta-conformation. Analysis of Asn-Asn doublets and Asn-X-Asn triplets permitted identification of conformational families in such sequences. Results of this investigation provide useful hints in modelling Asn-rich regions in proteins such as malaria parasite coat protein. (C) Munksgaard 1994.

Thermodynamic and structural studies of cavity formation in proteins suggest that loss of packing interactions rather than the hydrophobic effect dominates the observed energetics

The hydrophobic effect is widely believed to be an important determinant of protein stability. However, it is difficult to obtain unambiguous experimental estimates of the contribution of the hydrophobic driving force to the overall free energy of folding. Thermodynamic and structural studies of large to small substitutions in proteins are the most direct method of measuring this contribution. We have substituted the buried residue Phe8 in RNase S with alanine, methionine, and norleucine, Binding thermodynamics and structures were characterized by titration calorimetry and crystallography, respectively. The crystal structures of the RNase S F8A, F8M, and F8Nle mutants indicate that the protein tolerates the changes without any main chain adjustments, The correlation of structural and thermodynamic parameters associated with large to small substitutions was analyzed for nine mutants of RNase S as well as 32 additional cavity-containing mutants of T4 lysozyme, human lysozyme, and barnase. Such substitutions were typically found to result in negligible changes in Delta C-p and positive values of both Delta Delta H degrees and aas of folding. Enthalpic effects were dominant, and the sign of Delta Delta S is the opposite of that expected from the hydrophobic effect. Values of Delta Delta G degrees and Delta Delta H degrees correlated better with changes in packing parameters such as residue depth or occluded surface than with the change in accessible surface area upon folding. These results suggest that the loss of packing interactions rather than the hydrophobic effect is a dominant contributor to the observed energetics for large to small substitutions. Hence, estimates of the magnitude of the hydrophobic driving force derived from earlier mutational studies are likely to be significantly in excess of the actual value.