Ab initio protein folding simulations using atomic burials as informational intermediates between sequence and structure

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

EPR spectroscopic investigation of membrane protein structure and folding on light harvesting complex LHCIIb

Structure and folding of membrane proteins are important issues in molecular and cell biology. In this work new approaches are developed to characterize the structure of folded, unfolded and partially folded membrane proteins. These approaches combine site-directed spin labeling and pulse EPR techniques. The major plant light harvesting complex LHCIIb was used as a model system. Measurements of longitudinal and transversal relaxation times of electron spins and of hyperfine couplings to neighboring nuclei by electron spin echo envelope modulation(ESEEM) provide complementary information about the local environment of a single spin label. By double electron electron resonance (DEER) distances in the nanometer range between two spin labels can be determined. The results are analyzed in terms of relative water accessibilities of different sites in LHCIIb and its geometry. They reveal conformational changes as a function of micelle composition. This arsenal of methods is used to study protein folding during the LHCIIb self assembly and a spatially and temporally resolved folding model is proposed. The approaches developed here are potentially applicable for studying structure and folding of any protein or other self-assembling structure if site-directed spin labeling is feasible and the time scale of folding is accessible to freeze-quench techniques.

Dimensionamento preliminare di un UAV e progettazione di un sistema di folding wing

L'elaborato descrive la metodologia utilizzata per il dimensionamento preliminare di un velivolo a pilotaggio remoto caratterizzato da un sistema propulsivo costituito da due motori a benzina. Successivamente viene trattata la progettazione concettuale di metodi per la riduzione degli ingombri dovuti all'ala quando il velivolo è al suolo.

2D and 3D numerical modelling of multilayer detachment folding and salt tectonics

Numerical modelling was performed to study the dynamics of multilayer detachment folding and salt tectonics. In the case of multilayer detachment folding, analytically derived diagrams show several folding modes, half of which are applicable to crustal scale folding. 3D numerical simulations are in agreement with 2D predictions, yet fold interactions result in complex fold patterns. Pre-existing salt diapirs change folding patterns as they localize the initial deformation. If diapir spacing is much smaller than the dominant folding wavelength, diapirs appear in fold synclines or limbs.rnNumerical models of 3D down-building diapirism show that sedimentation rate controls whether diapirs will form and influences the overall patterns of diapirism. Numerical codes were used to retrodeform modelled salt diapirs. Reverse modelling can retrieve the initial geometries of a 2D Rayleigh-Taylor instability with non-linear rheologies. Although intermediate geometries of down-built diapirs are retrieved, forward and reverse modelling solutions deviate. rnFinally, the dynamics of fold-and-thrusts belts formed over a tilted viscous detachment is studied and it is demonstrated that mechanical stratigraphy has an impact on the deformation style, switching from thrust- to folding-dominated. The basal angle of the detachment controls the deformation sequence of the fold-and-thrust belt and results are consistent with critical wedge theory.rn

Apoptotic cell death in the lactating mammary gland is enhanced by a folding variant of alpha-lactalbumin

Apoptosis is essential to eliminate secretory epithelial cells during the involution of the mammary gland. The environmental regulation of this process is however, poorly understood. This study tested the effect of HAMLET (human alpha-lactalbumin made lethal to tumor cells) on mammary cells. Plastic pellets containing HAMLET were implanted into the fourth inguinal mammary gland of lactating mice for 3 days. Exposure of mammary tissue to HAMLET resulted in morphological changes typical for apoptosis and in a stimulation of caspase-3 activity in alveolar epithelial cells near the HAMLET pellets but not more distant to the pellet or in contralateral glands. The effect was specific for HAMLET and no effects were observed when mammary glands were exposed to native a-lactalbumin or fatty acid alone. HAMLET also induced cell death in vitro in a mouse mammary epithelial cell line. The results suggest that HAMLET can mediate apoptotic cell death in mammary gland tissue.

Compound heterozygous mutations affect protein folding and function in patients with congenital sucrase-isomaltase deficiency

BACKGROUND & AIMS: Congenital sucrase-isomaltase (SI) deficiency is an autosomal-recessive intestinal disorder characterized by a drastic reduction or absence of sucrase and isomaltase activities. Previous studies have indicated that single mutations underlie individual phenotypes of the disease. We investigated whether compound heterozygous mutations, observed in some patients, have a role in disease pathogenesis. METHODS: We introduced mutations into the SI complementary DNA that resulted in the amino acid substitutions V577G and G1073D (heterozygous mutations found in one group of patients) or C1229Y and F1745C (heterozygous mutations found in another group). The mutant genes were expressed transiently, alone or in combination, in COS cells and the effects were assessed at the protein, structural, and subcellular levels. RESULTS: The mutants SI-V577G, SI-G1073D, and SI-F1745C were misfolded and could not exit the endoplasmic reticulum, whereas SI-C1229Y was transported only to the Golgi apparatus. Co-expression of mutants found on each SI allele in patients did not alter the protein's biosynthetic features or improve its enzymatic activity. Importantly, the mutations C1229Y and F1745C, which lie in the sucrase domains of SI, prevented its targeting to the cell's apical membrane but did not affect protein folding or isomaltase activity. CONCLUSIONS: Compound heterozygosity is a novel pathogenic mechanism of congenital SI deficiency. The effects of mutations in the sucrase domain of SIC1229Y and SIF1745C indicate the importance of a direct interaction between isomaltase and sucrose and the role of sucrose as an intermolecular chaperone in the intracellular transport of SI.

Mutations affecting beta-tubulin folding and degradation.

Revertants of a colcemid-resistant Chinese hamster ovary cell line with an altered (D45Y) beta-tubulin have allowed the identification of four cis-acting mutations (L187R, Y398C, a 12-amino acid in-frame deletion, and a C-terminal truncation) that act by destabilizing the mutant tubulin and preventing it from incorporating into microtubules. These unstable beta-tubulins fail to form heterodimers and are predominantly found in association with the chaperonin CCT, suggesting that they cannot undergo productive folding. In agreement with these in vivo observations, we show that the defective beta-tubulins do not stably interact with cofactors involved in the tubulin folding pathway and, hence, fail to exchange with beta-tubulin in purified alphabeta heterodimers. Treatment of cells with MG132 causes an accumulation of the aberrant tubulins, indicating that improperly folded beta-tubulin is degraded by the proteasome. Rapid degradation of the mutant tubulin does not elicit compensatory changes in wild-type tubulin synthesis or assembly. Instead, loss of beta-tubulin from the mutant allele causes a 30-40% decrease in cellular tubulin content with no obvious effect on cell growth or survival.

STABILITY AND FOLDING OF MUTANT FORMS OF ISO-1 CYTOCHROME-C FROM SACCHAROMYCES CEREVISIAE

Partially functional forms of iso-1-cytochrome c from Saccharomyces cerevisiae were obtained by replacements of the evolutionarily conserved proline 71 with valine, isoleucine and threonine (Ernst et.al.,1985). Pro-71 lies at the juncture of two short helical regions and is believed to be important for proper local polypeptide chain folding within the iso-1-cytochrome c structure.^ To study folding in the absence of intermolecular disulfide dimer formation the free sulfhydryl group of Cys-102 was modified in both wild type and mutant proteins with an alkylating reagent, methyl methanethiosulfonate. Spectral analysis of the wild type and mutant proteins shows that the native-like functional (or partially functional) folded structure of cytochrome c is retained in the chemically modified derivatives. The replacement of Pro-71 with valine, isoleucine or threonine reduces the intensity of the 696 nm absorbance band which is an indicator of the Met-80 ligation to the heme. Thermal stability and guanidine hydrochloride unfolding studies of the mutant proteins shows a destabilization of the protein as a result of mutation. The degree of destabilization depends on the chemical nature of the substituent amino acid in the mutant protiens.^ Kinetics of folding/unfolding reactions of the proteins were monitored by fluorescence changes using stopped flow mixing to obtain guanidine hydrochloride concentration jumps ending below, within, and above the transition zone. The replacement of Pro-71 alters the rate on one of the fastest phases, $\tau\sb3$, while the two other phases, $\tau\sb1$ & $\tau\sb2$, remain the same.^ Slow refolding kinetic studies indicate that replacement of Pro-71 does not completely eliminate the absorbance or fluorescence detected slow phases leading to the conclusion that Pro-71 is not involved in the generation of the slow phases in the folding kinetics of iso-1-cytochrome c.^ The alkaline conformational change involving the disappearance of the 696 nm absorbance band occurs with increasing pH in the alkaline pH region (Davis et al., 1974). The apparent pK of this conformational change in mutant proteins is shifted as much as two pH units compared to wild type. The equilibrium and kinetic data of alkaline transition for the wild type follows a simple mechanism proposed by Davis et al., (1974) for horse heart cytochrome c. A more complex mechanism is proposed for the behavior of the mutant proteins. ^

IGHD II: A Novel GH-1 Gene Mutation (GH-L76P) Severely Affects GH Folding, Stability, and Secretion

CONTEXT The autosomal dominant form of GH deficiency (IGHD II) is characterized by markedly reduced GH secretion combined with low concentrations of IGF-1 leading to short stature. OBJECTIVE Structure-function analysis of a missense mutation in the GH-1 gene converting codon 76 from leucine (L) to proline (P) yielding a mutant GH-L76P peptide. DESIGN, SETTINGS, AND PATIENTS Heterozygosity for GH-L76P/wt-GH was identified in a nonconsanguineous Spanish family. The index patients, two siblings, a boy and a girl, were referred for assessment of their short stature (-3.2 and -3.8 SD). Their grandmother, father, and aunt were also carrying the same mutation and showed severe short stature; therefore, IGHD II was diagnosed. INTERVENTIONS AND RESULTS AtT-20 cells coexpressing both wt-GH and GH-L76P showed a reduced GH secretion (P < .001) after forskolin stimulation when compared with the cells expressing only wt-GH. In silico mutagenesis and molecular dynamics simulations presented alterations of correct folding and mutant stability compared with wt-GH. Therefore, further structural analysis of the GH-L76P mutant was performed using expressed and purified proteins in Escherichia coli by thermofluor assay and fast degradation proteolysis assay. Both assays revealed that the GH-L76P mutant is unstable and misfolded compared to wt-GH confirming the bioinformatic model prediction. CONCLUSIONS This is the first report of a family suffering from short stature caused by IGHD II, which severely affects intracellular GH folding and stability as well as secretion, highlighting the necessity of functional analysis of any GH variant for defining new mechanisms as a cause for IGHD II.