CaV3.1 channel pore pseudo-symmetry revealed by selectivity filter mutations in its domains I/II.

There is growing evidence indicating that the pore structure of voltage-gated ion channels (VGICs) influences gating besides their conductance. Regarding low voltage-activated (LVA) Ca2+ channels, it has been demonstrated that substitutions of the pore aspartate (D) by a glutamate (D-to-E substitution) in domains III and IV alter channel gating properties such as a positive shift in the channel activation voltage dependence. In the present report, we evaluated the effects of E-to-D substitution in domains I and II on the CaV3.1 channel gating properties. Our results indicate that substitutions in these two domains differentially modify the gating properties of CaV3.1 channels. The channel with a single mutation in domain I (DEDD) presented slower activation and faster inactivation kinetics and a slower recovery from inactivation, as compared with the WT channel. In contrast, the single mutant in domain II (EDDD) presented a small but significant negative shift of activation voltage dependence with faster activation and slower inactivation kinetics. Finally, the double mutant channel (DDDD) presented somehow intermediate properties with respect to the two single mutants but with fastest deactivation kinetics. Overall, our results indicate that single amino acid modification of the selectivity filter of LVA Ca2+ channels in distinct domains differentially influence their gating properties, supporting a pore pseudo-symmetry.

Neck linker docking is critical for Kinesin-1 force generation in cells but at a cost to motor speed and processivity.

Kinesin force generation involves ATP-induced docking of the neck linker (NL) along the motor core. However, the roles of the proposed steps of NL docking, cover-neck bundle (CNB) and asparagine latch (N-latch) formation, during force generation are unclear. Furthermore, the necessity of NL docking for transport of membrane-bound cargo in cells has not been tested. We generated kinesin-1 motors impaired in CNB and/or N-latch formation based on molecular dynamics simulations. The mutant motors displayed reduced force output and inability to stall in optical trap assays but exhibited increased speeds, run lengths, and landing rates under unloaded conditions. NL docking thus enhances force production but at a cost to speed and processivity. In cells, teams of mutant motors were hindered in their ability to drive transport of Golgi elements (high-load cargo) but not peroxisomes (low-load cargo). These results demonstrate that the NL serves as a mechanical element for kinesin-1 transport under physiological conditions.

Mitochondrial ROS regulate thermogenic energy expenditure and sulfenylation of UCP1.

Brown and beige adipose tissues can dissipate chemical energy as heat through thermogenic respiration, which requires uncoupling protein 1 (UCP1). Thermogenesis from these adipocytes can combat obesity and diabetes, encouraging investigation of factors that control UCP1-dependent respiration in vivo. Here we show that acutely activated thermogenesis in brown adipose tissue is defined by a substantial increase in levels of mitochondrial reactive oxygen species (ROS). Remarkably, this process supports in vivo thermogenesis, as pharmacological depletion of mitochondrial ROS results in hypothermia upon cold exposure, and inhibits UCP1-dependent increases in whole-body energy expenditure. We further establish that thermogenic ROS alter the redox status of cysteine thiols in brown adipose tissue to drive increased respiration, and that Cys253 of UCP1 is a key target. UCP1 Cys253 is sulfenylated during thermogenesis, while mutation of this site desensitizes the purine-nucleotide-inhibited state of the carrier to adrenergic activation and uncoupling. These studies identify mitochondrial ROS induction in brown adipose tissue as a mechanism that supports UCP1-dependent thermogenesis and whole-body energy expenditure, which opens the way to improved therapeutic strategies for combating metabolic disorders.

Protein:Ligand binding free energies: A stringent test for computational protein design.

A computational protein design method is extended to allow Monte Carlo simulations where two ligands are titrated into a protein binding pocket, yielding binding free energy differences. These provide a stringent test of the physical model, including the energy surface and sidechain rotamer definition. As a test, we consider tyrosyl-tRNA synthetase (TyrRS), which has been extensively redesigned experimentally. We consider its specificity for its substrate l-tyrosine (l-Tyr), compared to the analogs d-Tyr, p-acetyl-, and p-azido-phenylalanine (ac-Phe, az-Phe). We simulate l- and d-Tyr binding to TyrRS and six mutants, and compare the structures and binding free energies to a more rigorous "MD/GBSA" procedure: molecular dynamics with explicit solvent for structures and a Generalized Born + Surface Area model for binding free energies. Next, we consider l-Tyr, ac- and az-Phe binding to six other TyrRS variants. The titration results are sensitive to the precise rotamer definition, which involves a short energy minimization for each sidechain pair to help relax bad contacts induced by the discrete rotamer set. However, when designed mutant structures are rescored with a standard GBSA energy model, results agree well with the more rigorous MD/GBSA. As a third test, we redesign three amino acid positions in the substrate coordination sphere, with either l-Tyr or d-Tyr as the ligand. For two, we obtain good agreement with experiment, recovering the wildtype residue when l-Tyr is the ligand and a d-Tyr specific mutant when d-Tyr is the ligand. For the third, we recover His with either ligand, instead of wildtype Gln.

Phase transitions and structure analysis in wild-type, A30P, E46K, and A53T mutants of α-synuclein.

α-Synuclein has been implicated in the development of neural plaques in Parkinson's Disease and Lewy-Body Dementia. This paper reports on the structural phase change behavior exhibited over a relevant range of temperatures in canonical protein Monte Carlo simulations for wild-type α-synuclein and three of its familial variants. We performed and analyzed these simulations to determine residue occupancy variations above and below this phase transition. From this analysis, we found regions above the phase transition temperature that consistently exhibited increased propensity for formation of long-chain beta-sheets, suggesting a possible role in α-synuclein aggregation.

Assessment of Mutational Effects on Peptide Stability through Confinement Simulations.

The evaluation of free energy differences between specific states of a system is of fundamental interest in the study of (bio)chemical systems. Herein, we examine the use of the recently introduced confinement method (CM) to evaluate relative free energy changes upon protein/peptide mutations. CM is a path-independent technique that involves the transformation of a configurational state of the system into an ideal crystal permitting the direct computation of free energy differences. We illustrate the method by evaluating the differential stabilities between native and mutant sequences of a model peptide that has been extensively characterized by experimental approaches, the GB1 hairpin. We show a good correlation between calculated and experimental relative stabilities and discuss other possible applications of this method in the context of complex molecular conversions.

Computational study of the W260A activating mutant of Src tyrosine kinase.

Tyrosine kinases are enzymes playing a critical role in cellular signaling. Mutations causing increased in kinase activity are often associated with cancer and various pathologies. One example in Src tyrosine kinases is offered by the substitution of the highly conserved tryptophan 260 by an alanine (W260A), which has been shown to cause an increase in activity. Here, molecular dynamics simulations based on atomic models are carried out to characterize the conformational changes in the linker region and the catalytic (kinase) domain of Src kinase to elucidate the impact of the W260A mutation. Umbrella sampling calculations show that the conformation of the linker observed in the assembled down-regulated state of the kinase is most favored when the kinase domain is in the inactive state, whereas the conformation of the linker observed in the re-assembled up-regulated state of the kinase is favored when the kinase domain is in the unphosphorylated active-like state. The calculations further indicate that there are only small differences between the WT and W260A mutant. In both cases, the intermediates states are very similar and the down-regulated inactive conformation is the most stable state. However, the calculations also show that the free energy cost to reach the unphosphorylated active-like conformation is slightly smaller for the W260A mutant compared with WT. A simple kinetic model is developed and submitted to a Bayesian Monte Carlo analysis to illustrate how such small differences can contribute to accelerate the trans-autophosphorylation reaction and yield a large increase in the activity of the mutant as observed experimentally.

Crystal structure of the dynamin tetramer.

The mechanochemical protein dynamin is the prototype of the dynamin superfamily of large GTPases, which shape and remodel membranes in diverse cellular processes. Dynamin forms predominantly tetramers in the cytosol, which oligomerize at the neck of clathrin-coated vesicles to mediate constriction and subsequent scission of the membrane. Previous studies have described the architecture of dynamin dimers, but the molecular determinants for dynamin assembly and its regulation have remained unclear. Here we present the crystal structure of the human dynamin tetramer in the nucleotide-free state. Combining structural data with mutational studies, oligomerization measurements and Markov state models of molecular dynamics simulations, we suggest a mechanism by which oligomerization of dynamin is linked to the release of intramolecular autoinhibitory interactions. We elucidate how mutations that interfere with tetramer formation and autoinhibition can lead to the congenital muscle disorders Charcot-Marie-Tooth neuropathy and centronuclear myopathy, respectively. Notably, the bent shape of the tetramer explains how dynamin assembles into a right-handed helical oligomer of defined diameter, which has direct implications for its function in membrane constriction.

Cytosine unstacking and strand slippage at an insertion-deletion mutation sequence in an overhang-containing DNA duplex.

Base unstacking in template strands, when accompanied by strand slippage, can result in deletion mutations during strand extension by nucleic acid polymerases. In a GCCC mutation hot-spot sequence, which was previously identified to have a 50% probability of causing such mutations during DNA replication by a Y-family polymerase, a single-base deletion mutation could result from such unstacking of any one of its three template cytosines. In this study, the intrinsic energetic differences in unstacking among these three cytosines in a solvated DNA duplex overhang model were examined using umbrella sampling molecular dynamics simulations. The free energy profiles obtained show that cytosine unstacking grows progressively more unfavorable as one moves inside the duplex from the 5'-end of the overhang template strand. Spontaneous strand slippage occurs in response to such base unstacking in the direction of both the major and minor grooves for all three cytosines. Unrestrained simulations run from three distinct strand-slipped states and one non-strand-slipped state suggest that a more duplexlike environment can help stabilize strand slippage. The possible underlying reasons and biological implications of these observations are discussed in the context of nucleic acid replication active site dynamics.

Controlling dissociation channels of gas-phase protein complexes using charge manipulation.

Coarse-grained simulations with charge hopping were performed for a positively charged tetrameric transthyretin (TTR) protein complex with a total charge of +20. Charges were allowed to move among basic amino acid sites as well as N-termini. Charge distributions and radii of gyration were calculated for complexes simulated at two temperatures, 300 and 600 K, under different scenarios. One scenario treated the complex in its normal state allowing charge to move to any basic site. Another scenario blocked protonation of all the N-termini except one. A final scenario used the complex in its normal state but added a basic-site containing tether (charge tag) near the N-terminus of one chain. The differences in monomer unfolding and charging were monitored in all three scenarios and compared. The simulation results show the importance of the N-terminus in leading the unfolding of the monomer units; a process that follows a zipper-like mechanism. Overall, experimentally modifying the complex by adding a tether or blocking the protonation of N-termini may give the potential for controlling the unraveling and subsequent dissociation of protein complexes.

Genetic assessment of the role of AcrB ß-hairpins in the assembly of the TolC-AcrAB multidrug efflux pump of Escherichia coli.

The tripartite AcrAB-TolC multidrug efflux pump of Escherichia coli is the central conduit for cell-toxic compounds and contributes to antibiotic resistance. While high-resolution structures of all three proteins have been solved, much remains to be learned as to how the individual components come together to form a functional complex. In this study, we investigated the importance of the AcrB ß-hairpins belonging to the DN and DC subdomains, which are presumed to dock with TolC, in complex stability and activity of the complete pump. Our data show that the DN subdomain ß-hairpin residues play a more critical role in complex stability and activity than the DC subdomain hairpin residues. The failure of the AcrB DN ß-hairpin deletion mutant to engage with TolC leads to the drug hypersensitivity phenotype, which is reversed by compensatory alterations in the lipoyl and ß-barrel domains of AcrA. Moreover, AcrA and TolC mutants that induce TolC opening also reverse the drug hypersensitivity phenotype of the AcrB ß-hairpin mutants, indicating a failure by the AcrB mutant to interact and thus induce TolC opening on its own. Together, these data suggest that both AcrB ß-hairpins and AcrA act to stabilize the tripartite complex and induce TolC opening for drug expulsion.

Effect of proline mutations on the monomer conformations of amylin.

The formation of human islet amyloid polypeptide (hIAPP) is implicated in the loss of pancreatic ß-cells in type II diabetes. Rat amylin, which differs from human amylin at six residues, does not lead to formation of amyloid fibrils. Pramlintide is a synthetic analog of human amylin that shares three proline substitutions with rat amylin. Pramlintide has a much smaller propensity to form amyloid aggregates and has been widely prescribed in amylin replacement treatment. It is known that the three prolines attenuate ß-sheet formation. However, the detailed effects of these proline substitutions on full-length hIAPP remain poorly understood. In this work, we use molecular simulations and bias-exchange metadynamics to investigate the effect of proline substitutions on the conformation of the hIAPP monomer. Our results demonstrate that hIAPP can adopt various ß-sheet conformations, some of which have been reported in experiments. The proline substitutions perturb the formation of long ß-sheets and reduce their stability. More importantly, we find that all three proline substitutions of pramlintide are required to inhibit ß conformations and stabilize the α-helical conformation. Fewer substitutions do not have a significant inhibiting effect.

Aggregation formation in the polyglutamine diseases: protection at a cost?

Mutant protein aggregation is a hallmark of many neurodegenerative diseases, including the polyglutamine disorders. Although the correlation between aggregation formation and disease pathology originally suggested that the visible inclusions seen in patient tissue might directly contribute to pathology, additional studies failed to confirm this hypothesis. Current opinion in the field of polyglutamine disease research now favors a model in which large inclusions are cytoprotective and smaller oligomers or misfolded monomers underlie pathogenesis. Nonetheless, therapies aimed at reducing or preventing aggregation show promise. This review outlines the debate about the role of aggregation in the polyglutamine diseases as it has unfolded in the literature and concludes with a brief discussion on the manipulation of aggregation formation and clearance mechanisms as a means of therapeutic intervention.

Mixed Monte Carlo/Molecular Dynamics simulations of the prion protein.

In this paper we present the results of mixed Monte Carlo/Molecular Dynamics (MC/MD) simulations of the D178N mutant of the human prion protein. We have used the MC moves for polypeptide sampling known as Concerted Rotations with Angles (CRA) to selectively sample the region of the prion protein comprising the ß-sheet and one of the α-helices. The results indicate that the MC/MD simulations sample the phase space substantially faster than regular Molecular Dynamics simulations starting with the same initial conditions. This work further indicates the MC/MD technique as a potentially powerful simulation tool, allowing the selective sampling of a region of a physical system that is deemed important.

Comprehensive functional assessment of MLH1 variants of unknown significance.

Lynch syndrome is associated with germline mutations in DNA mismatch repair (MMR) genes. Up to 30% of DNA changes found are variants of unknown significance (VUS). Our aim was to assess the pathogenicity of eight MLH1 VUS identified in patients suspected of Lynch syndrome. All of them are novel or not previously characterized. For their classification, we followed a strategy that integrates family history, tumor pathology, and control frequency data with a variety of in silico and in vitro analyses at RNA and protein level, such as MMR assay, MLH1 and PMS2 expression, and subcellular localization. Five MLH1 VUS were classified as pathogenic: c.[248G>T(;)306G>C], c.[780C>G;788A>C], and c.791-7T>A affected mRNA processing, whereas c.218T>C (p.L73P) and c.244A>G [corrected] (p.T82A) impaired MMR activity. Two other VUS were considered likely neutral: the silent c.702G>A variant did not affect mRNA processing or stability, and c.974G>A (p.R325Q) did not influence MMR function. In contrast, variant c.25C>T (p.R9W) could not be classified, as it associated with intermediate levels of MMR activity. Comprehensive functional assessment of MLH1 variants was useful in their classification and became relevant in the diagnosis and genetic counseling of carrier families.

Expression of recombinant Ara h 6 in Pichia pastoris but not in Escherichia coli preserves allergic effector function and allows assessment of specific mutations.

SCOPE: Ara h 6 has recently been recognized as an important peanut allergen. Recombinant allergens have been used for analysis of IgE binding, but have not been used to analyze the allergic effector activity that is more relevant to allergic reactions. METHODS AND RESULTS: Ara h 6 was expressed as a recombinant protein in both Escherichia coli and Pichia pastoris (rAra h 6-E. coli and rAra h 6-Pichia, respectively). Effector activity was assayed by measuring degranulation of RBL SX-38 cells sensitized with IgE from patients with severe peanut allergy. Compared to native Ara h 6 (nAra h 6), rAra h 6-Pichia had intact effector function whereas rAra h 6-E. coli had significantly reduced function. The lower effector activity in rAra h 6-E. coli compared to nAra h 6 and rAra h 6-Pichia did not appear to be due to differences in posttranslational modifications (analyzed by mass spectrometry and staining for carbohydrates) and may be due to subtle alteration(s) of folding seen on CD analysis and on nonreduced gels. Finally, we introduced point mutations in four important IgE-binding linear epitopes of Ara h 6 and found dramatically reduced allergic effector activity. CONCLUSION: Our studies demonstrate the utility of fully functional rAra h 6-Pichia as a starting point for analysis of specific mutations that adversely affect allergic effector function.

Mechanistic assessment of DNA ligase as an antibacterial target in Staphylococcus aureus.

We report the use of a known pyridochromanone inhibitor with antibacterial activity to assess the validity of NAD(+)-dependent DNA ligase (LigA) as an antibacterial target in Staphylococcus aureus. Potent inhibition of purified LigA was demonstrated in a DNA ligation assay (inhibition constant [K(i)] = 4.0 nM) and in a DNA-independent enzyme adenylation assay using full-length LigA (50% inhibitory concentration [IC(50)] = 28 nM) or its isolated adenylation domain (IC(50) = 36 nM). Antistaphylococcal activity was confirmed against methicillin-susceptible and -resistant S. aureus (MSSA and MRSA) strains (MIC = 1.0 µg/ml). Analysis of spontaneous resistance potential revealed a high frequency of emergence (4 × 10(-7)) of high-level resistant mutants (MIC > 64) with associated ligA lesions. There were no observable effects on growth rate in these mutants. Of 22 sequenced clones, 3 encoded point substitutions within the catalytic adenylation domain and 19 in the downstream oligonucleotide-binding (OB) fold and helix-hairpin-helix (HhH) domains. In vitro characterization of the enzymatic properties of four selected mutants revealed distinct signatures underlying their resistance to inhibition. The infrequent adenylation domain mutations altered the kinetics of adenylation and probably elicited resistance directly. In contrast, the highly represented OB fold domain mutations demonstrated a generalized resistance mechanism in which covalent LigA activation proceeds normally and yet the parameters of downstream ligation steps are altered. A resulting decrease in substrate K(m) and a consequent increase in substrate occupancy render LigA resistant to competitive inhibition. We conclude that the observed tolerance of staphylococcal cells to such hypomorphic mutations probably invalidates LigA as a viable target for antistaphylococcal chemotherapy.

Quantitative analysis of ligand migration from transition networks.

In this work we use transition network analysis for the first time to investigate ligand migration in truncated hemoglobin (trHbN) and obtain kinetic information about the docking-site dynamics in the protein. A comparison with explicit water molecular dynamics simulations (100 ns in total) shows that the rate constants derived from the network analysis are realistic. The transition network analysis provides 1) The time-resolved connectivity network in the protein; 2) The half-lives of the docking sites; 3) The transition timescales between two given docking sites; and 4) The extent of population transfer among different docking sites of the protein as a function of lag time. We investigate the role of the Tyr33 and Gln58 residues in ligand migration by studying ligand migration in four mutants of trHbN. The mutation study suggests that residues Tyr33 and Gln58 stabilize the NO ligand in the Xe2 docking site of trHbN, thus facilitating the efficiency of the NO detoxification reaction.

The prediction of protein-protein interacting sites in genome-wide protein interaction networks: the test case of the human cell cycle.

In this paper we aim at investigating possible correlations between the number of putative interaction patches of a given protein, as inferred by an algorithm that we have developed, and its degree (number of edges of the protein node in a protein interaction network). We focus on the human cell cycle that, as compared with other biological processes, comprises the largest number of proteins whose structure is known at atomic resolution both as monomers and as interacting complexes. For predicting interaction patches we specifically develop a HM-SVM based method reaching 71% overall accuracy with a correlation coefficient value equal to 0.43 on a non redundant set of protein complexes. To test the biological meaning of our predictions, we also explore whether interacting patches contain energetically important residues and/or disease related mutations and find that predicted patches are endowed with both features. Based on this, we propose that mapping the protein with all the predicted interaction patches bridges the molecule to the interactome at the cell level. To test our hypothesis we downloaded interaction data from interaction data bases and find that the number of predicted interaction patches significantly correlates (Pearson correlation value >0.3) with the number of the known interactions (edges) per protein in the human interactome, as contained in MINT and IntAct. We also show that the correlation increases (Pearson correlation value >0.5) when the subcellular co-localization and the co-expression levels of the interacting partners are taken into account.

Disease risk of missense mutations using structural inference from predicted function.

Advancements in sequencing techniques place personalized genomic medicine upon the horizon, bringing along the responsibility of clinicians to understand the likelihood for a mutation to cause disease, and of scientists to separate etiology from nonpathologic variability. Pathogenicity is discernable from patterns of interactions between a missense mutation, the surrounding protein structure, and intermolecular interactions. Physicochemical stability calculations are not accessible without structures, as is the case for the vast majority of human proteins, so diagnostic accuracy remains in infancy. To model the effects of missense mutations on functional stability without structure, we combine novel protein sequence analysis algorithms to discern spatial distributions of sequence, evolutionary, and physicochemical conservation, through a new approach to optimize component selection. Novel components include a combinatory substitution matrix and two heuristic algorithms that detect positions which confer structural support to interaction interfaces. The method reaches 0.91 AUC in ten-fold cross-validation to predict alteration of function for 6,392 in vitro mutations. For clinical utility we trained the method on 7,022 disease associated missense mutations within the Online Mendelian inheritance in man amongst a larger randomized set. In a blinded prospective test to delineate mutations unique to 186 patients with craniosynostosis from those in the 95 highly variant Coriell controls and 1000 age matched controls, we achieved roughly 1/3 sensitivity and perfect specificity. The component algorithms retained during machine learning constitute novel protein sequence analysis techniques to describe environments supporting neutrality or pathology of mutations. This approach to pathogenetics enables new insight into the mechanistic relationship of missense mutations to disease phenotypes in our patients.