Molecular dynamics study of a calmodulin-like protein with an IQ peptide: spontaneous refolding of the protein around the peptide

The Calmodulin (CaM) is a small (16.7 kDa), highly acidic protein that is crucial to all eukaryotes by serving as a prototypical calcium sensor. In the present study, we investigated, through molecular dynamics simulations, the dynamics of a complex between the Mlc1p protein, which is a CaM-like protein, and the IQ4 peptide. This protein-peptide interaction is of high importance because IQ motifs are widely distributed among different kinds of CaM-binding proteins. The Mlc1p-IQ4 complex, which had been resolved by crystallography to 2.1 A, confers to a Ca(+2)-independent stable structure. During the simulations, the complex undergoes a complicated modulation process, which involves bending of the angles between the alpha-helices of the protein, breaking of the alpha-helical structure of the IQ4 peptide into two sections, and formation of new contact points between the protein and the peptide. The dynamics of the process consist of fast sub picosecond events and much slower ones that take a few nanoseconds to completion. Our study expands the information embedded in the crystal structure of the Mlc1p-IQ4 complex by describing its dynamic behavior as it evolves from the crystal structure to a form stable in solution. The article shows that careful application of molecular dynamics simulations can be used for extending the structural information presented by the crystal structure, thereby revealing the dynamic configuration of the protein in its physiological environment.     

Early Closure of a Long Loop in the Refolding of Adenylate Kinase: A Possible Key Role of Non-Local Interactions in the Initial Folding Steps

Most globular protein chains, when transferred from high to low denaturant concentrations, collapse instantly before they refold to their native state. The initial compaction of the protein molecule is assumed to have a key effect on the folding pathway, but it is not known whether the earliest structures formed during or instantly after collapse are defined by local or by non-local interactions—that is, by secondary structural elements or by loop closure of long segments of the protein chain. Stable closure of one or several long loops can reduce the chain entropy at a very early stage and can prevent the protein from following non-productive pathways whose number grows exponentially with the length of the protein chain. In Escherichia coli adenylate kinase (AK), about seven long loops define the topology of the native structure. We selected four loop-forming sections of the chain and probed the time course of loop formation during refolding of AK. We labeled the termini of the loop segments with tryptophan and cysteine-5-amidosalicylic acid. This donor-acceptor pair of probes used with fluorescence resonance excitation energy transfer spectroscopy (FRET) is suitable for detecting very short distances and thus is able to distinguish between random and specific compactions. Refolding of AK was initiated by stopped-flow mixing, followed simultaneously by donor and acceptor fluorescence, and analyzed in terms of energy transfer efficiency and distance. In the collapsed state of AK, observed after the 5-ms dead time of the instrument, one of the selected segments shows a native-like separation of its termini; it forms a loop already in the collapsed state. A second segment that includes the first but is longer by 15 residues shows an almost native-like separation of its termini. In contrast, a segment that is shorter but part of the second segment shows a distance separation of its termini as high as a segment that spans almost the whole protein chain. We conclude that a specific network of non-local interactions, the closure of one or several loops, can play an important role in determining the protein folding pathway at its early phases.     

Structural coupling between the Rho-insert domain of Cdc42 and the geranylgeranyl binding site of RhoGDI

The small GTPase proteins are components of the intracellular signaling system, alternating between active (membrane-bound and GTP-loaded) and inactive (GDP-loaded and cytosolic) states. In the inactive state, the proteins are soluble in the cytoplasm. To compensate for the energetic penalty of extraction of the hydrophobic moiety from the membrane phase, the inactive state is stabilized via formation of a complex with the RhoGDI proteins that provide a hydrophobic pocket for the binding of the hydrophobic moieties. The signals delivered by the Rho subfamily involve a specific, short, highly exposed α-helix (Rho-insert), located close to the GDP binding site. Upon simulating the complex in solution, we observed that the Rho-insert domain of Cdc42 can assume two basic orientations. One is the canonical one, as detected in both crystals and NMR spectra of concentrated protein solutions. The second orientation appears only in the RhoGDI-Cdc42 complex where the GER moiety of Cdc42 is properly inserted into the specific binding site of RhoGDI. Any impairment of the GER-RhoGDI interactions, such as insertion of specific mutations in the hydrophobic binding site, abolished the coupling between the proteins and the Rho-insert domain, preserving its canonical orientation as in the crystalline structure. The noncanonical conformation of the Rho-insert domain is not a simulation artifact, as it appears in crystals of plant Rho proteins (ROP4, ROP5, and ROP7). In accord with the notion that the Rho-insert domain participates in downstream signaling, we propose that the deformation of the Rho-insert is part of the signal transmissions.     

Discovery and Statistical Genotyping of Copy-Number Variation from Whole-Exome Sequencing Depth

Sequencing of gene-coding regions (the exome) is increasingly used for studying human disease, for which copy-number variants (CNVs) are a critical genetic component. However, detecting copy number from exome sequencing is challenging because of the noncontiguous nature of the captured exons. This is compounded by the complex relationship between read depth and copy number; this results from biases in targeted genomic hybridization, sequence factors such as GC content, and batching of samples during collection and sequencing. We present a statistical tool (exome hidden Markov model [XHMM]) that uses principal-component analysis (PCA) to normalize exome read depth and a hidden Markov model (HMM) to discover exon-resolution CNV and genotype variation across samples. We evaluate performance on 90 schizophrenia trios and 1,017 case-control samples. XHMM detects a median of two rare (<1%) CNVs per individual (one deletion and one duplication) and has 79% sensitivity to similarly rare CNVs overlapping three or more exons discovered with microarrays. With sensitivity similar to state-of-the-art methods, XHMM achieves higher specificity by assigning quality metrics to the CNV calls to filter out bad ones, as well as to statistically genotype the discovered CNV in all individuals, yielding a trio call set with Mendelian-inheritance properties highly consistent with expectation. We also show that XHMM breakpoint quality scores enable researchers to explicitly search for novel classes of structural variation. For example, we apply XHMM to extract those CNVs that are highly likely to disrupt (delete or duplicate) only a portion of a gene.     

Phenotypic changes in satellite glial cells in cultured trigeminal ganglia

Satellite glial cells (SGCs) are specialized cells that form a tight sheath around neurons in sensory ganglia. In recent years, there is increasing interest in SGCs and they have been studied in both intact ganglia and in tissue culture. Here we studied phenotypic changes in SGCs in cultured trigeminal ganglia from adult mice, containing both neurons and SGCs, using phase optics, immunohistochemistry and time-lapse photography. Cultures were followed for up to 14 days. After isolation virtually every sensory neuron is ensheathed by SGCs, as in the intact ganglia. After one day in culture, SGCs begin to migrate away from their parent neurons, but in most cases the neurons still retain an intact glial cover. At later times in culture, there is a massive migration of SGCs away from the neurons and they undergo clear morphological changes, and at 7 days they become spindle-shaped. At one day in culture SGCs express the glial marker glutamine synthetase, and also the purinergic receptor P2X7. From day 2 in culture the glutamine synthetase expression is greatly diminished, whereas that of P2X7 is largely unchanged. We conclude that SGCs retain most of their characteristics for about 24?h after culturing, but undergo major phenotypic changes at later times.     

Effect of smoking on circulating angiogenic factors in high risk pregnancies

Objective

Changes in maternal concentrations of the anti-angiogenic factors, soluble fms-like tyrosine kinase 1 (sFlt1) and soluble endoglin (sEng), and the pro-angiogenic placental growth factor (PlGF) precede the development of preeclampsia in healthy women. The risk of preeclampsia is reduced in women who smoke during pregnancy. The objective of this study was to investigate whether smoking affects concentrations of angiogenic factors (sFlt1, PlGF, and sEng) in women at high risk for developing preeclampsia.

Study Design

We performed a secondary analysis of serum samples from 993 high-risk women (chronic hypertension, diabetes, multifetal gestation, and previous preeclampsia) in a preeclampsia prevention trial. sFlt1, sEng and PlGF were measured in serum samples obtained at study entry, which was prior to initiation of aspirin (median 19.0 weeks'' [interquartile range of 16.0–22.6 weeks'']). Smoking status was determined by self-report.

Results

sFlt1 was not significantly different in smokers from any high-risk groups compared to their nonsmoking counterparts. PlGF was higher among smokers compared to nonsmokers among diabetic women (142.7 [77.4–337.3] vs 95.9 [48.5–180.7] pg/ml, p = 0.005) and women with a history of preeclampsia (252.2 [137.1–486.0] vs 152.2 [73.6–253.7] pg/ml, p = 0.001). sEng was lower in smokers with multifetal gestations (5.8 [4.6–6.5] vs 6.8 [5.5–8.7] ng/ml, p = 0.002) and trended lower among smokers with diabetes (4.9 [3.8–5.6] vs 5.3 [4.3–6.3] ng/ml, p = 0.05). Smoking was not associated with a lower incidence of preeclampsia in any of these groups.

Conclusions

In certain high-risk groups, smoking is associated with changes in the concentrations of these factors towards a pro-angiogenic direction during early pregnancy; however, there was no apparent association between smoking and the development of preeclampsia in our cohort.     

Insight into the interaction sites between fatty acid binding proteins and their ligands

Fatty acid binding proteins (FABPs), are evolutionarily conserved small cytoplasmic proteins that occur in many tissue-specific types. One of their primary functions is to facilitate the clearance of the cytoplasmic matrix from free fatty acids and of other detergent-like compounds. Crystallographic studies of FABP proteins have revealed a well defined binding site located deep inside their β-clam structure that is hardly exposed to the bulk solution. However, NMR measurements revealed that, when the protein is equilibrated with its ligands, residues that are clearly located on the outer surface of the protein do interact with the ligand. To clarify this apparent contradiction we applied molecular dynamics simulations to follow the initial steps associated with the FABP–fatty acid interaction using, as a model, the interaction of toad liver basic FABP, or chicken liver bile acid binding protein, with a physiological concentration of palmitate ions. The simulations (~200 ns of accumulated time) show that fatty acid molecules interact, unevenly, with various loci on the protein surface, with the favored regions being the portal and the anti-portal domains. Random encounters with palmitate at these regions led to lasting adsorption to the surface, while encounters at the outer surface of the β-clam were transient. Therefore, we suggest that the protein surface is capable of sequestering free fatty acids from solution, where brief encounters evolve into adsorbed states, which later mature by migration of the ligand into a more specific binding site.     

Structural basis of binding by cyclic nonphosphorylated peptide antagonists of Grb7 implicated in breast cancer progression

Growth-receptor-bound protein (Grb)7 is an adapter protein aberrantly overexpressed, along with the erbB-2 receptor in breast cancer and in other cancers. Normally recruited to focal adhesions with a role in cell migration, it is associated with erbB-2 in cancer cells and is found to exacerbate cancer progression via stimulation of cell migration and proliferation. The G7-18NATE peptide (sequence: WFEGYDNTFPC cyclized via a thioether bond) is a nonphosphorylated peptide that was developed for the specific inhibition of Grb7 by blocking its SH2 domain. Cell-permeable versions of G7-18NATE are effective in the reduction of migration and proliferation in Grb7-overexpressing cells. It thus represents a promising starting point for the development of a therapeutic against Grb7. Here, we report the crystal structure of the G7-18NATE peptide in complex with the Grb7-SH2 domain, revealing the structural basis for its interaction. We also report further rounds of phage display that have identified G7-18NATE analogues with micromolar affinity for Grb7-SH2. These peptides retained amino acids F2, G4, and F9, as well as the YDN motif that the structural biology study showed to be the main residues in contact with the Grb7-SH2 domain. Isothermal titration calorimetry measurements reveal similar and better binding affinity of these peptides compared with G7-18NATE. Together, this study facilitates the optimization of second-generation inhibitors of Grb7.     

Cognitive effects of cellular phones: a possible role of non-radiofrequency radiation factors

Some studies found that cognitive functions of human beings may be altered while exposed to radiofrequency radiation (RFR) emitted by cellular phones. In two recent studies, we have found that experiment duration and exposure side (i.e., phone's location--right or left) may have a major influence on the detection of such effects. In this brief follow-up experiment, 29 right-handed male subjects were divided into two groups. Each subject had two standard cellular phones attached to both sides of his head. The subjects performed a spatial working memory task that required either a left-hand or a right-hand response under one of the two exposure conditions: left side of the head or right side. Contrary to our previous studies, in this work external antennas located far away from the subjects were connected to the cellular phones. This setup prevents any emission of RFR from the internal antenna, thus drastically reducing RFR exposure. Despite that, the results remain similar to those obtained in our previous work. These results indicate that some of the effects previously attributed to RFR can be the result of some confounders.     

Cognitive effects of radiation emitted by cellular phones: The influence of exposure side and time

This study examined the time dependence effects of exposure to radiofrequency radiation (RFR) emitted by standard GSM cellular phones on the cognitive functions of humans. A total of 48 healthy right‐handed male subjects performed a spatial working memory task (that required either a left‐hand or a right‐hand response) while being exposed to one of two GSM phones placed at both sides of the head. The subjects were randomly divided into three groups. Each group was exposed to one of three exposure conditions: left‐side of the head, right‐side, or sham‐exposure. The experiment consisted of 12 blocks of trials. Response times (RTs) and accuracy of the responses were recorded. It was found that the average RT of the right‐hand responses under left‐side exposure condition was significantly longer than those of the right‐side and sham‐exposure groups averaged together during the first two time blocks. These results confirmed the existence of an effect of exposure on RT, as well as the fact that exposure duration (together with the responding hand and the side of exposure) may play an important role in producing detectable RFR effects on performance. Differences in these parameters might be the reason for the failure of certain studies to detect or replicate RFR effects. Bioelectromagnetics 30:198–204, 2009. © 2008 Wiley‐Liss, Inc.