Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background: Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue. Results: Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues. Conclusion: We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Germline BRCA1 mutation analysis in Indian breast/ovarian cancer families

Most of the predisposition to hereditary breast and ovarian cancer has been attributed to inherited defects in two tumor suppressor genes BRCA1 and BRCA2. To explore the contribution of BRCA1 mutations to hereditary breast cancer among Indian women, we examined the coding sequence of the BRCA1 gene in 14 breast cancer patients with a positive family history of breast and/or ovarian cancer. Mutation analysis was carried out using conformation sensitive gel electrophoresis (CSGE) followed by sequencing. Three mutations (21%) in the BRCA1 gene were identified. Two of them are novel mutations of which one is a missense mutation in exon 7 near the RING finger domain, while the other is a one base pair deletion in exon 11 which results in protein truncation. The third mutation, 185delAG, has been previously described in Ashkenazi Jewish families. To our knowledge this is the first report of a study of germline BRCA1 mutation analysis in familial breast cancer in India. Our data from 14 different families suggests a lower prevalence but definite involvement of germline mutations in the BRCA1 gene among Indian women with breast cancer and a family history of breast cancer.

Clinical phenotype and the lack of mutations in the CHRNG, CHRND, and CHRNA1 genes in two Indian families with Escobar syndrome

The objective of this study was to report the clinical phenotype and genetic analysis of two Indian families with Escobar syndrome (ES). The diagnosis of ES in both families was made on the basis of published clinical features. Blood samples were collected from members of both families and used in genomic DNA isolation. The entire coding regions and intron-exon junctions of the ES gene CHRNG (cholinergic receptor, nicotinic, gamma), and two other related genes, CHRND and CHRNA1, were amplified and sequenced to search for mutations in both families. Both families show a typical form of ES. Sequencing of the entire coding regions including the intron-exon junctions of the three genes did not yield any mutations in these families. In conclusion, it is possible that the mutations in these genes are located in the promoter or deep intronic regions that we failed to identify or the ES in these families is caused by mutations in a different gene. The lack of mutations in CHRNG has also been reported in several families, suggesting the possibility of at least one more gene for this syndrome. Clin Dysmorphol 22:54-58 (C) 2013 Wolters Kluwer Health vertical bar Lippincott Williams & Wilkins.

Investigations of Ramachandran disallowed conformations in protein domain families

In peptide and protein structures, occurrence of (phi,psi.) angles in the disallowed region of the Ramachandran map almost always suggests local regions of error or poor accuracy. However, very rarely genuine disallowed conformations occur as noted in the current study in proteins of known structure available at ultra-high resolution (<= 1.2 (A) over circle). In the current work, extent of conservation of genuine disallowed conformations in homologous proteins of known structures has been analyzed. From a dataset of 124 protein domain families, with structure of at least one constituent member in each family available at a resolution of 1.2 (A) over circle or better, we have analyzed the conservation of 221 disallowed conformations. It is observed that the disallowed conformation is only moderately conservedin protein domain families. In the gross dataset no particular residue type adopting disallowed conformation elicit high conservation of residue type though there are alignment positions in the dataset with complete conservation of both the residue type and the disallowed conformation. Conserved disallowed conformation in protein domain families play biologically significant role in roughly 50% of the cases. The residues with the disallowed conformation or its flanking residues are often located within or around the functional site of the protein. (C) 2013 Elsevier B.V. All rights reserved.

Minimizing total weighted tardiness on heterogeneous batch processors with incompatible job families

In this paper, we address a scheduling problem for minimizing total weighted tardiness. The background for the paper is derived from the automobile gear manufacturing process. We consider the bottleneck operation of heat treatment stage of gear manufacturing. Real-life scenarios like unequal release times, incompatible job families, nonidentical job sizes, heterogeneous batch processors, and allowance for job splitting have been considered. We have developed a mathematical model which takes into account dynamic starting conditions. The problem considered in this study is NP-hard in nature, and hence heuristic algorithms have been proposed to address it. For real-life large-size problems, the performance of the proposed heuristic algorithms is evaluated using the method of estimated optimal solution available in literature. Extensive computational analyses reveal that the proposed heuristic algorithms are capable of consistently obtaining near-optimal statistically estimated solutions in very reasonable computational time.

A Global Constraint On Relative Rotation To Avoid Lumped Compliant Mechanisms In Topology Optimization

Some of the well known formulations for topology optimization of compliant mechanisms could lead to lumped compliant mechanisms. In lumped compliance, most of the elastic deformation in a mechanism occurs at few points, while rest of the mechanism remains more or less rigid. Such points are referred to as point-flexures. It has been noted in literature that high relative rotation is associated with point-flexures. In literature we also find a formulation of local constraint on relative rotations to avoid lumped compliance. However it is well known that a global constraint is easier to handle than a local constraint, by a numerical optimization algorithm. The current work presents a way of putting global constraint on relative rotations. This constraint is also simpler to implement since it uses linearized rotation at the center of finite-elements, to compute relative rotations. I show the results obtained by using this constraint oil the following benchmark problems - displacement inverter and gripper.

Wave propagation in multi-walled carbon nanotube

In this paper, wave propagation in multi-walled carbon nanotubes (MWNTs) are studied by modeling them as continuum multiple shell coupled through van der Waals force of interaction. The displacements, namely, axial, radial and circumferential displacements vary along the circumferential direction. The wave propagation are simulated using the wavelet based spectral finite element (WSFE) method. This technique involves Daubechies scaling function approximation in time and spectral element approach. The WSFE Method allows the study of wave properties in both time and frequency domains. This is in contrast to the conventional Fourier transform based analysis which are restricted to frequency domain analysis. Here, first, the wavenumbers and wave speeds of carbon nanotubes (CNTs) are Studied to obtain the characteristics of the waves. These group speeds have been compared with those reported in literature. Next, the natural frequencies of a single-walled carbon nanotube (SWNT) are studied for different values of the radius. The frequencies of the first five modes vary linearly with the radius of the SWNT. Finally, the time domain responses are simulated for SWNT and three-walled carbon nanotubes.

Numerical studies on columnar-to-equiaxed transition in directional solidification of binary alloys

A numerical study on columnar-to-equiaxed transition (CET) during directional solidification of binary alloys is presented using a macroscopic solidification model. The position of CET is predicted numerically using a critical cooling rate criterion reported in literature. The macroscopic solidification model takes into account movement of solid phase due to buoyancy, and drag effect on the moving solid phase because of fluid motion. The model is applied to simulate the solidification process for binary alloys (Sn-Pb) and to estimate solidification parameters such as position of the liquidus, velocity of the liquidus isotherm, temperature gradient ahead of the liquidus, and cooling rate at the liquidus. Solidification phenomena under two cooling configurations are studied: one without melt convection and the other involvin thermosolutal convection. The numerically predicted positions of CET compare well with those of experiments reported in literature. Melt convection results in higher cooling rate, higher liquidus isotherm velocities, and stimulation of occurrence of CET in comparison to the nonconvecting case. The movement of solid phase aids further the process of CET. With a fixed solid phase, the occurrence of CET based on the same critical cooling rate is delayed and it occurs at a greater distance from the chill.

Free convection heat transfer to mercury in vertical annuli

Data on free convection heat transfer to water and mercury are collected using a test rig in vertical annuli of three radii ratios, the walls of which are maintained at uniform temperatures. A theoretical analysis of the boundary layer equations has been attempted using local similarity transformation and double boundary layer approach. Correlations derived from the present theoretical analysis are compared with the analysis and the experimental data available in literature for non-metallic fluids and also with the present experimental data on water and mercury. Generalised correlations are set up for expressing the ratio of heat transferred by convection to the heat transferred by pure conduction and Nusselt's number, in terms of Grashof, Rayleigh and Prandtl numbers, based on the theoretical analysis and the present data on mercury and water. The present generalised correlations agree with the reported and present data for non-metallic fluids and liquid metals with an average deviation of 9% and maximum deviation of ± 13.7%.

Breakage of viscous and non-Newtonian drops in stirred dispersions

A model of breakage of drops in a stirred vessel has been proposed to account for the effect of rheology of the dispersed phase. The deformation of the drop is represented by a Voigt element. A realistic description of the role of interfacial tension is incorporated by treating it as a restoring force which passes through a maximum as the drop deforms and eventually reaching a zero value at the break point. It is considered that the drop will break when the strain of the drop has reached a value equal to its diameter. An expression for maximum stable drop diameter, dmax, is derived from the model and found to be applicable over a wide range of variables, as well as to data already existing in literature. The model could be naturally extended to predict observed values of dmax when the dispersed phase is a power law fluid or a Bingham plastic.

Analysis of human chorionic gonadotropin-monoclonal antibody interaction in BIAcore

Kinetic studies of macromolecular ligand-ligate interaction have generated ample interest since the advent of plasmon resonance based instruments like BIAcore. Most of the studies reported in literature assume a simple 1 : 1 Langmuir binding and complete reversibility of the system. However we observed that in a high affinity antigen-antibody system [human chorionic gonadotropin-monoclonal antibody (hCG-mAb)] dissociation is insignificant and the sensogram data cannot be used to measure the equilibrium and kinetic parameters. At low concentrations of mAb the complete sensogram could be fitted to a single exponential. Interestingly we found that at higher mAb concentrations, the binding data did not conform to a simple bimolecular model. Instead, the data fitted a two-step model, which may be because of surface heterogeneity of affinity sites. In this paper, we report on the global fit of the sensograms. We have developed a method by which a single two-minute sensogram can be used in high affinity systems to measure the association rate constant of the reaction and the functional capacity of the ligand (hCG) immobilized on the chip. We provide a rational explanation for the discrepancies generally observed in most of the BIAcore sensograms

Behaviour of Bi-material interface cracks in the presence of material nonlinearity of adherends

The well known features of crack face interpenetration/contact at the tip of an interface crack is re-examined using finite element analysis and assuming material nonlinear properties for the adherends. It was assumed in literature that the crack tips are fully open at all load levels in the presence of material nonlinearity of the adherends. Analysis for the case of remote tension shows that even in the presence of material nonlinearity, crack tip closes at small load levels and opens above a certain load level. Mixed-mode fracture parameters are evaluated for the situation when the crack tips are fully open. Due to the presence of nonlinearity, the mixed-mode fracture parameters are measured with the symmetric and anti-symmetric components of J-integral. The present analysis explains the sequence of events at the interface crack tip with progressively increasing remote tension load for the case of adherends with material nonlinear behaviour.

Stability of an unsupported vertical circular excavation in clays under undrained condition

By using the axisymmetric finite elements static limit analysis formulation, proposed recently by the authors, the stability numbers (gamma H/c(o)) for an unsupported vertical circular excavation in clays, whose cohesion increases with depth, have been determined under undrained condition; gamma = unit weight, H., height of the excavation and c(o) = cohesion along ground surface. The results are obtained for various values of H/b and m; where b = the radius of the excavation and m = a non-dimensional parameter which accounts for the rate of the increase of cohesion with depth. The values of the stability numbers increase continuously both with increases in H/b and m. The results obtained in this study compare well with those available in literature.(C) 2009 Elsevier Ltd. All rights reserved.