Bi-allelic SNAPC4 variants dysregulate global alternative splicing and lead to neuroregression and progressive spastic paraparesis.

The vast majority of human genes encode multiple isoforms through alternative splicing, and the temporal and spatial regulation of those isoforms is critical for organismal development and function. The spliceosome, which regulates and executes splicing reactions, is primarily composed of small nuclear ribonucleoproteins (snRNPs) that consist of small nuclear RNAs (snRNAs) and protein subunits. snRNA gene transcription is initiated by the snRNA-activating protein complex (SNAPc). Here, we report ten individuals, from eight families, with bi-allelic, deleterious SNAPC4 variants. SNAPC4 encoded one of the five SNAPc subunits that is critical for DNA binding. Most affected individuals presented with delayed motor development and developmental regression after the first year of life, followed by progressive spasticity that led to gait alterations, paraparesis, and oromotor dysfunction. Most individuals had cerebral, cerebellar, or basal ganglia volume loss by brain MRI. In the available cells from affected individuals, SNAPC4 abundance was decreased compared to unaffected controls, suggesting that the bi-allelic variants affect SNAPC4 accumulation. The depletion of SNAPC4 levels in HeLa cell lines via genomic editing led to decreased snRNA expression and global dysregulation of alternative splicing. Analysis of available fibroblasts from affected individuals showed decreased snRNA expression and global dysregulation of alternative splicing compared to unaffected cells. Altogether, these data suggest that these bi-allelic SNAPC4 variants result in loss of function and underlie the neuroregression and progressive spasticity in these affected individuals.

Newly identified disorder of copper metabolism caused by variants in CTR1, a high-affinity copper transporter.

The high-affinity copper transporter CTR1 is encoded by CTR1 (SLC31A1), a gene locus for which no detailed genotype-phenotype correlations have previously been reported. We describe identical twin male infants homozygous for a novel missense variant NM_001859.4:c.284 G > A (p.Arg95His) in CTR1 with a distinctive autosomal recessive syndrome of infantile seizures and neurodegeneration, consistent with profound central nervous system copper deficiency. We used clinical, biochemical and molecular methods to delineate the first recognized examples of human CTR1 deficiency. These included clinical phenotyping, brain imaging, assays for copper, cytochrome c oxidase (CCO), and mitochondrial respiration, western blotting, cell transfection experiments, confocal and electron microscopy, protein structure modeling and fetal brain and cerebral organoid CTR1 transcriptome analyses. Comparison with two other critical mediators of cellular copper homeostasis, ATP7A and ATP7B, genes associated with Menkes disease and Wilson disease, respectively, revealed that expression of CTR1 was highest. Transcriptome analyses identified excitatory neurons and radial glia as brain cell types particularly enriched for copper transporter transcripts. We also assessed the effects of Copper Histidinate in the patients' cultured cells and in the patients, under a formal clinical protocol. Treatment normalized CCO activity and enhanced mitochondrial respiration in vitro, and was associated with modest clinical improvements. In combination with present and prior studies, these infants' clinical, biochemical and molecular phenotypes establish the impact of this novel variant on copper metabolism and cellular homeostasis and illuminate a crucial role for CTR1 in human brain development. CTR1 deficiency represents a newly defined inherited disorder of brain copper metabolism.

Modulation of free energy landscapes as a strategy for the design of antimicrobial peptides.

Computational design of antimicrobial peptides (AMPs) is a promising area of research for developing novel agents against drug-resistant bacteria. AMPs are present naturally in many organisms, from bacteria to humans, a time-tested mechanism that makes them attractive as effective antibiotics. Depending on the environment, AMPs can exhibit α-helical or ß-sheet conformations, a mix of both, or lack secondary structure; they can be linear or cyclic. Prediction of their structures is challenging but critical for rational design. Promising AMP leads can be developed using essentially two approaches: traditional modeling of the physicochemical mechanisms that determine peptide behavior in aqueous and membrane environments and knowledge-based, e.g., machine learning (ML) techniques, that exploit ever-growing AMP databases. Here, we explore the conformational landscapes of two recently ML-designed AMPs, characterize the dependence of these landscapes on the medium conditions, and identify features in peptide and membrane landscapes that mediate protein-membrane association. For both peptides, we observe greater conformational diversity in an aqueous solvent than in a less polar solvent, and one peptide is seen to alter its conformation more dramatically than the other upon the change of solvent. Our results support the view that structural rearrangement in response to environmental changes is central to the mechanism of membrane-structure disruption by linear peptides. We expect that the design of AMPs by ML will benefit from the incorporation of peptide conformational substates as quantified here with molecular simulations.

A point mutation in the nuclease domain of MLH3 eliminates repeat expansions in a mouse stem cell model of the Fragile X-related disorders.

The Fragile X-related disorders (FXDs) are Repeat Expansion Diseases, genetic disorders that result from the expansion of a disease-specific microsatellite. In those Repeat Expansion Disease models where it has been examined, expansion is dependent on functional mismatch repair (MMR) factors, including MutLγ, a heterodimer of MLH1/MLH3, one of the three MutL complexes found in mammals and a minor player in MMR. In contrast, MutLα, a much more abundant MutL complex that is the major contributor to MMR, is either not required for expansion or plays a limited role in expansion in many model systems. How MutLγ acts to generate expansions is unclear given its normal role in protecting against microsatellite instability and while MLH3 does have an associated endonuclease activity, whether that contributes to repeat expansion is uncertain. We show here, using a gene-editing approach, that a point mutation that eliminates the endonuclease activity of MLH3 eliminates expansions in an FXD mouse embryonic stem cell model. This restricts the number of possible models for repeat expansion and supports the idea that MutLγ may be a useful druggable target to reduce somatic expansion in those disorders where it contributes to disease pathology.

Nanosecond-Timescale Dynamics and Conformational Heterogeneity in Human GCK Regulation and Disease.

Human glucokinase (GCK) is the prototypic example of an emerging class of proteins with allosteric-like behavior that originates from intrinsic polypeptide dynamics. High-resolution NMR investigations of GCK have elucidated millisecond-timescale dynamics underlying allostery. In contrast, faster motions have remained underexplored, hindering the development of a comprehensive model of cooperativity. Here, we map nanosecond-timescale dynamics and structural heterogeneity in GCK using a combination of unnatural amino acid incorporation, time-resolved fluorescence, and 19F nuclear magnetic resonance spectroscopy. We find that a probe inserted within the enzyme's intrinsically disordered loop samples multiple conformations in the unliganded state. Glucose binding and disease-associated mutations that suppress cooperativity alter the number and/or relative population of these states. Together, the nanosecond kinetics characterized here and the millisecond motions known to be essential for cooperativity provide a dynamical framework with which we address the origins of cooperativity and the mechanism of activated, hyperinsulinemia-associated, noncooperative variants.

HLA and autoantibodies define scleroderma subtypes and risk in African and European Americans and suggest a role for molecular mimicry.

Systemic sclerosis (SSc) is a clinically heterogeneous autoimmune disease characterized by mutually exclusive autoantibodies directed against distinct nuclear antigens. We examined HLA associations in SSc and its autoantibody subsets in a large, newly recruited African American (AA) cohort and among European Americans (EA). In the AA population, the African ancestry-predominant HLA-DRB1*08:04 and HLA-DRB1*11:02 alleles were associated with overall SSc risk, and the HLA-DRB1*08:04 allele was strongly associated with the severe antifibrillarin (AFA) antibody subset of SSc (odds ratio = 7.4). These African ancestry-predominant alleles may help explain the increased frequency and severity of SSc among the AA population. In the EA population, the HLA-DPB1*13:01 and HLA-DRB1*07:01 alleles were more strongly associated with antitopoisomerase (ATA) and anticentromere antibody-positive subsets of SSc, respectively, than with overall SSc risk, emphasizing the importance of HLA in defining autoantibody subtypes. The association of the HLA-DPB1*13:01 allele with the ATA+ subset of SSc in both AA and EA patients demonstrated a transancestry effect. A direct correlation between SSc prevalence and HLA-DPB1*13:01 allele frequency in multiple populations was observed (r = 0.98, P = 3 × 10-6). Conditional analysis in the autoantibody subsets of SSc revealed several associated amino acid residues, mostly in the peptide-binding groove of the class II HLA molecules. Using HLA α/ß allelic heterodimers, we bioinformatically predicted immunodominant peptides of topoisomerase 1, fibrillarin, and centromere protein A and discovered that they are homologous to viral protein sequences from the Mimiviridae and Phycodnaviridae families. Taken together, these data suggest a possible link between HLA alleles, autoantibodies, and environmental triggers in the pathogenesis of SSc.

Biallelic HEPHL1 variants impair ferroxidase activity and cause an abnormal hair phenotype.

Maintenance of the correct redox status of iron is functionally important for critical biological processes. Multicopper ferroxidases play an important role in oxidizing ferrous iron, released from the cells, into ferric iron, which is subsequently distributed by transferrin. Two well-characterized ferroxidases, ceruloplasmin (CP) and hephaestin (HEPH) facilitate this reaction in different tissues. Recently, a novel ferroxidase, Hephaestin like 1 (HEPHL1), also known as zyklopen, was identified. Here we report a child with compound heterozygous mutations in HEPHL1 (NM_001098672) who presented with abnormal hair (pili torti and trichorrhexis nodosa) and cognitive dysfunction. The maternal missense mutation affected mRNA splicing, leading to skipping of exon 5 and causing an in-frame deletion of 85 amino acids (c.809_1063del; p.Leu271_ala355del). The paternal mutation (c.3176T>C; p.Met1059Thr) changed a highly conserved methionine that is part of a typical type I copper binding site in HEPHL1. We demonstrated that HEPHL1 has ferroxidase activity and that the patient's two mutations exhibited loss of this ferroxidase activity. Consistent with these findings, the patient's fibroblasts accumulated intracellular iron and exhibited reduced activity of the copper-dependent enzyme, lysyl oxidase. These results suggest that the patient's biallelic variants are loss-of-function mutations. Hence, we generated a Hephl1 knockout mouse model that was viable and had curly whiskers, consistent with the hair phenotype in our patient. These results enhance our understanding of the function of HEPHL1 and implicate altered ferroxidase activity in hair growth and hair disorders.

Static Kinks or Flexible Hinges: Multiple Conformations of Bent DNA Bound to Integration Host Factor Revealed by Fluorescence Lifetime Measurements.

Gene regulation depends on proteins that bind to specific DNA sites. Such specific recognition often involves severe DNA deformations, including sharp kinks. It has been unclear how rigid or flexible these protein-induced kinks are. Here, we investigated the dynamic nature of DNA in complex with integration host factor (IHF), a nucleoid-associated architectural protein known to bend one of its cognate sites (35 base pair H') into a U-turn by kinking DNA at two sites. We utilized fluorescence-lifetime-based FRET spectroscopy to assess the distribution of bent conformations in various IHF-DNA complexes. Our results reveal a surprisingly dynamic specific complex: while 78% of the IHF-H' population exhibited FRET efficiency consistent with the crystal structure, 22% exhibited FRET efficiency indicative of unbent or partially bent DNA. This conformational flexibility is modulated by sequence variations in the cognate site. In another site (H1) that lacks the A-tract of H' found on one side of the binding site, the extent of bending in the fully U-bent conformation decreased, and the population in that state decreased to 32%. A similar decrease in the U-bent population was observed with a single base mutation in H' in a consensus region on the other side. Taken together, these results provide important insights into the finely tuned interactions between IHF and its cognate sites that keep the DNA bent (or not) and yield quantitative data on the dynamic equilibrium between different DNA conformations (kinked or not kinked) that depend sensitively on DNA sequence and deformability. Notably, the difference in dynamics between IHF-H' and IHF-H1 reflects the different roles of these complexes in their natural context, in the phage lambda "intasome" (the complex that integrates phage lambda into the E. coli chromosome).

Peptide and Protein Structure Prediction with a Simplified Continuum Solvent Model.

A continuum solvent model based on screened Coulomb potentials has been simplified and parametrized to sample native-like structures in replica-exchange simulations of each of six different peptides and miniproteins. Low-energy, native, and non-native structures were used to iteratively refine 11 parameter values. The centroid of the largest cluster of structures sampled in simulations initiated from an extended conformation represents the predicted structure. The main-chain rms deviation of this prediction from the experimental structure was 0.47 Å for the 12-residue Trp-zip2, 0.86 Å for the 14-residue MBH12, 2.53 Å for the 17-residue U(1-17)T9D, 2.03 Å for the 20-residue BS1, 1.08 Å for the 20-residue Trp-cage, and 3.64 Å for the 35-residue villin headpiece subdomain HP35. The centroid of the sixth largest cluster sampled for HP35 deviated by 0.91 Å. The CHARMM22/CMAP force field was used, with an additional ψ torsion term for residues other than glycine and proline. Six parameters govern the dielectric response of the continuum solvent, and four values of surface tension approximate nonpolar effects. An atom's self-energy and interaction energies are screened independently, each depending on whether the atom is part of a charged group, a neutral hydrogen-bonding main-chain group, or any other neutral group. The parameters inferred result in strong main-chain hydrogen bonds, consistent with the view that protein folding is dominated by the formation of these bonds. (1,2) Conformations of MBH12 and BS1 were excluded from the energy-function refinement, suggesting the parameters, referred to as SCP18, are transferable. An efficient estimate of solvent-accessible surface area is also described.

Loss of function mutations in VARS encoding cytoplasmic valyl-tRNA synthetase cause microcephaly, seizures, and progressive cerebral atrophy.

Progressive microcephaly and neurodegeneration are genetically heterogenous conditions, largely associated with genes that are essential for the survival of neurons. In this study, we interrogate the genetic etiology of two siblings from a non-consanguineous family with severe early onset of neurological manifestations. Whole exome sequencing identified novel compound heterozygous mutations in VARS that segregated with the proband: a missense (c.3192G>A; p.Met1064Ile) and a splice site mutation (c.1577-2A>G). The VARS gene encodes cytoplasmic valyl-tRNA synthetase (ValRS), an enzyme that is essential during eukaryotic translation. cDNA analysis on patient derived fibroblasts revealed that the splice site acceptor variant allele led to nonsense mediated decay, thus resulting in a null allele. Three-dimensional modeling of ValRS predicts that the missense mutation lies in a highly conserved region and could alter side chain packing, thus affecting tRNA binding or destabilizing the interface between the catalytic and tRNA binding domains. Further quantitation of the expression of VARS showed remarkably reduced levels of mRNA and protein in skin derived fibroblasts. Aminoacylation experiments on patient derived cells showed markedly reduced enzyme activity of ValRS suggesting the mutations to be loss of function. Bi-allelic mutations in cytoplasmic amino acyl tRNA synthetases are well-known for their role in neurodegenerative disorders, yet human disorders associated with VARS mutations have not yet been clinically well characterized. Our study describes the phenotype associated with recessive VARS mutations and further functional delineation of the pathogenicity of novel variants identified, which widens the clinical and genetic spectrum of patients with progressive microcephaly.

Enhanced spontaneous DNA twisting/bending fluctuations unveiled by fluorescence lifetime distributions promote mismatch recognition by the Rad4 nucleotide excision repair complex.

Rad4/XPC recognizes diverse DNA lesions including ultraviolet-photolesions and carcinogen-DNA adducts, initiating nucleotide excision repair. Studies have suggested that Rad4/XPC senses lesion-induced helix-destabilization to flip out nucleotides from damaged DNA sites. However, characterizing how DNA deformability and/or distortions impact recognition has been challenging. Here, using fluorescence lifetime measurements empowered by a maximum entropy algorithm, we mapped the conformational heterogeneities of artificially destabilized mismatched DNA substrates of varying Rad4-binding specificities. The conformational distributions, as probed by FRET between a cytosine-analog pair exquisitely sensitive to DNA twisting/bending, reveal a direct connection between intrinsic DNA deformability and Rad4 recognition. High-specificity CCC/CCC mismatch, free in solution, sampled a strikingly broad range of conformations from B-DNA-like to highly distorted conformations that resembled those observed with Rad4 bound; the extent of these distortions increased with bound Rad4 and with temperature. Conversely, the non-specific TAT/TAT mismatch had a homogeneous, B-DNA-like conformation. Molecular dynamics simulations also revealed a wide distribution of conformations for CCC/CCC, complementing experimental findings. We propose that intrinsic deformability promotes Rad4 damage recognition, perhaps by stalling a diffusing protein and/or facilitating 'conformational capture' of pre-distorted damaged sites. Surprisingly, even mismatched DNA specifically bound to Rad4 remains highly dynamic, a feature that may reflect the versatility of Rad4/XPC to recognize many structurally dissimilar lesions.

CELSR2, encoding a planar cell polarity protein, is a putative gene in Joubert syndrome with cortical heterotopia, microophthalmia, and growth hormone deficiency.

Joubert syndrome is a ciliopathy characterized by a specific constellation of central nervous system malformations that result in the pathognomonic "molar tooth sign" on imaging. More than 27 genes are associated with Joubert syndrome, but some patients do not have mutations in any of these genes. Celsr1, Celsr2, and Celsr3 are the mammalian orthologues of the drosophila planar cell polarity protein, flamingo; they play important roles in neural development, including axon guidance, neuronal migration, and cilium polarity. Here, we report bi-allelic mutations in CELSR2 in a Joubert patient with cortical heterotopia, microophthalmia, and growth hormone deficiency. © 2017 Wiley Periodicals, Inc.

Cytoskeletal abnormalities and neutrophil dysfunction in WDR1 deficiency.

Cell motility, division, and structural integrity depend on dynamic remodeling of the cellular cytoskeleton, which is regulated in part by actin polymerization and depolymerization. In 3 families, we identified 4 children with recurrent infections and varying clinical manifestations including mild neutropenia, impaired wound healing, severe stomatitis with oral stenosis, and death. All patients studied had similar distinctive neutrophil herniation of the nuclear lobes and agranular regions within the cytosol. Chemotaxis and chemokinesis were markedly impaired, but staphylococcal killing was normal, and neutrophil oxidative burst was increased both basally and on stimulation. Neutrophil spreading on glass and cell polarization were also impaired. Neutrophil F-actin was elevated fourfold, suggesting an abnormality in F-actin regulation. Two-dimensional differential in-gel electrophoresis identified abnormal actin-interacting protein 1 (Aip1), encoded by WDR1, in patient samples. Biallelic mutations in WDR1 affecting distinct antiparallel ß-strands of Aip1 were identified in all patients. It has been previously reported that Aip1 regulates cofilin-mediated actin depolymerization, which is required for normal neutrophil function. Heterozygous mutations in clinically normal relatives confirmed that WDR1 deficiency is autosomal recessive. Allogeneic stem cell transplantation corrected the immunologic defect in 1 patient. Mutations in WDR1 affect neutrophil morphology, motility, and function, causing a novel primary immunodeficiency.

Structural Basis of Species-Dependent Differential Affinity of 6-Alkoxy-5-Aryl-3-Pyridinecarboxamide Cannabinoid-1 Receptor Antagonists.

6-Alkoxy-5-aryl-3-pyridincarboxamides, including the brain-penetrant compound 14G: [5-(4-chlorophenyl)-6-(cyclopropylmethoxy)-N-[(1R,2R)-2-hydroxy-cyclohexyl]-3-pyridinecarboxamide] and its peripherally restricted analog 14H: [5-(4-chlorophenyl)-N-[(1R,2R)-2-hydroxycyclohexyl]-6-(2-methoxyethoxy)-3-pyridinecarboxamide], have been recently introduced as selective, high-affinity antagonists of the human cannabinoid-1 receptor (hCB1R). Binding analyses revealed two orders of magnitude lower affinity of these compounds for mouse and rat versus human CB1R, whereas the affinity of rimonabant is comparable for all three CB1Rs. Modeling of ligand binding to CB1R and binding assays with native and mutant (Ile105Met) hCB1Rs indicate that the Ile105 to Met mutation in rodent CB1Rs accounts for the species-dependent affinity of 14G: and 14H: . Our work identifies Ile105 as a new pharmacophore component for developing better hCB1R antagonists and invalidates rodent models for assessing the antiobesity efficacy of 14G: and 14H: .

Drastic neofunctionalization associated with evolution of the timezyme AANAT 500 Mya.

Melatonin (N-acetyl-5-methoxytrypamine) is the vertebrate hormone of the night: circulating levels at night are markedly higher than day levels. This increase is driven by precisely regulated increases in acetylation of serotonin in the pineal gland by arylalkylamine N-acetyltransferase (AANAT), the penultimate enzyme in the synthesis of melatonin. This unique essential role of AANAT in vertebrate timekeeping is recognized by the moniker the timezyme. AANAT is also found in the retina, where melatonin is thought to play a paracrine role. Here, we focused on the evolution of AANAT in early vertebrates. AANATs from Agnathans (lamprey) and Chondrichthyes (catshark and elephant shark) were cloned, and it was found that pineal glands and retinas from these groups express a form of AANAT that is compositionally, biochemically, and kinetically similar to AANATs found in bony vertebrates (VT-AANAT). Examination of the available genomes indicates that VT-AANAT is absent from other forms of life, including the Cephalochordate amphioxus. Phylogenetic analysis and evolutionary rate estimation indicate that VT-AANAT evolved from the nonvertebrate form of AANAT after the Cephalochordate-Vertebrate split over one-half billion years ago. The emergence of VT-AANAT apparently involved a dramatic acceleration of evolution that accompanied neofunctionalization after a duplication of the nonvertebrate AANAT gene. This scenario is consistent with the hypotheses that the advent of VT-AANAT contributed to the evolution of the pineal gland and lateral eyes from a common ancestral photodetector and that it was not a posthoc recruitment.

Identifying protein kinase target preferences using mass spectrometry.

A general question in molecular physiology is how to identify candidate protein kinases corresponding to a known or hypothetical phosphorylation site in a protein of interest. It is generally recognized that the amino acid sequence surrounding the phosphorylation site provides information that is relevant to identification of the cognate protein kinase. Here, we present a mass spectrometry-based method for profiling the target specificity of a given protein kinase as well as a computational tool for the calculation and visualization of the target preferences. The mass spectrometry-based method identifies sites phosphorylated in response to in vitro incubation of protein mixtures with active recombinant protein kinases followed by standard phosphoproteomic methodologies. The computational tool, called "PhosphoLogo," uses an information-theoretic algorithm to calculate position-specific amino acid preferences and anti-preferences from the mass-spectrometry data (http://helixweb.nih.gov/PhosphoLogo/). The method was tested using protein kinase A (catalytic subunit α), revealing the well-known preference for basic amino acids in positions -2 and -3 relative to the phosphorylated amino acid. It also provides evidence for a preference for amino acids with a branched aliphatic side chain in position +1, a finding compatible with known crystal structures of protein kinase A. The method was also employed to profile target preferences and anti-preferences for 15 additional protein kinases with potential roles in regulation of epithelial transport: CK2, p38, AKT1, SGK1, PKCδ, CaMK2δ, DAPK1, MAPKAPK2, PKD3, PIM1, OSR1, STK39/SPAK, GSK3ß, Wnk1, and Wnk4.

Molecular and biochemical characterization of a unique mutation in CCS, the human copper chaperone to superoxide dismutase.

Copper (Cu) is a trace metal that readily gains and donates electrons, a property that renders it desirable as an enzyme cofactor but dangerous as a source of free radicals. To regulate cellular Cu metabolism, an elaborate system of chaperones and transporters has evolved, although no human Cu chaperone mutations have been described to date. We describe a child from a consanguineous family who inherited homozygous mutations in the SLC33A1, encoding an acetyl CoA transporter, and in CCS, encoding the Cu chaperone for superoxide dismutase. The CCS mutation, p.Arg163Trp, predicts substitution of a highly conserved arginine residue at position 163, with tryptophan in domain II of CCS, which interacts directly with superoxide dismutase 1 (SOD1). Biochemical analyses of the patient's fibroblasts, mammalian cell transfections, immunoprecipitation assays, and Lys7Δ (CCS homolog) yeast complementation support the pathogenicity of the mutation. Expression of CCS was reduced and binding of CCS to SOD1 impaired. As a result, this mutation causes reduced SOD1 activity and may impair other mechanisms important for normal Cu homeostasis. CCS-Arg163Trp represents the primary example of a human mutation in a gene coding for a Cu chaperone.

Filtering artifacts from lifetime distributions when maximizing entropy using a bootstrapped model.

The maximum entropy method (MEM) has been used in many studies to reliably recover effective lifetimes from kinetics, whether measured experimentally or simulated computationally. Here, recent claims made by Mulligan et al. regarding MEM analyses of kinetics (Anal. Biochem. 421 (2012) 181-190) are shown to be unfounded. Their assertion that their software allows "analysis of datasets too noisy to process by existing iterative search algorithms" is refuted with a MEM analysis of their triexponential test case with increased noise. In addition, it is shown that lifetime distributions recovered from noisy kinetics data with the MEM can be improved by using a simple filter when bootstrapping the prior model. When deriving the bootstrapped model from the lifetime distribution obtained using a uniform model, only the slower processes are represented as Gaussians in the bootstrapped model. Using this new approach, results are clearly superior to those of Mulligan et al. despite the presence of increased noise. In a second example, ambiguity in the interpretation of Poisson kinetics in the presence of scattered excitation light is resolved by filtering the prior model.

Water-exclusion and liquid-structure forces in implicit solvation.

A continuum model of solvation is proposed to describe (i) long-range electrostatic effects of water exclusion resulting from incomplete and anisotropic hydration in crowded environments and (ii) short-range effects of liquid-structure forces on the hydrogen-bond interactions at solute/water interfaces. The model is an extension of the phenomenological screened coulomb potential-based implicit model of solvation. The developments reported here allow a more realistic representation of highly crowded and spatially heterogeneous environments, such as those in the interior of a living cell. Only the solvent is treated as a continuum medium. It is shown that the electrostatic effects of long-range water-exclusion can strongly affect protein-protein binding energies and are then related to the thermodynamics of complex formation. Hydrogen-bond interactions modulated by the liquid structure at interfaces are calibrated based on systematic calculations of potentials of mean force in explicit water. The electrostatic component of the model is parametrized for monovalent, divalent and trivalent ions. The conceptual and practical aspects of the model are discussed based on simulations of protein complexation and peptide folding. The current implementation is ~1.5 times slower than the gas-phase force field and exhibits good parallel performance.

The Tensin-3 protein, including its SH2 domain, is phosphorylated by Src and contributes to tumorigenesis and metastasis.

In cell lines from advanced lung cancer, breast cancer, and melanoma, endogenous tensin-3 contributes to cell migration, anchorage-independent growth, and tumorigenesis. Although SH2 domains have not been reported previously to be phosphorylated, the tensin-3 SH2 domain is a physiologic substrate for Src. Tyrosines in the SH2 domain contribute to the biological activity of tensin-3, and phosphorylation of these tyrosines can regulate ligand binding. In a mouse breast cancer model, tensin-3 tyrosines are phosphorylated in a Src-associated manner in primary tumors, and experimental metastases induced by tumor-derived cell lines depend on endogenous tensin-3. Thus, tensin-3 is implicated as an oncoprotein regulated by Src and possessing an SH2 domain with a previously undescribed mechanism for the regulation of ligand binding.