Computational design of a thermolabile uracil-DNA glycosylase of Escherichia coli.

Polymerase chain reaction (PCR) is a powerful tool to diagnose infectious diseases. Uracil DNA glycosylase (UDG) is broadly used to remove carryover contamination in PCR. However, UDG can contribute to false negative results when not inactivated completely, leading to DNA degradation during the amplification step. In this study, we designed novel thermolabile UDG derivatives by supercomputing molecular dynamic simulations and residual network analysis. Based on enzyme activity analysis, thermolability, thermal stability, and biochemical experiments of Escherichia coli-derived UDG and 22 derivatives, we uncovered that the UDG D43A mutant eliminated the false negative problem, demonstrated high efficiency, and offered great benefit for use in PCR diagnosis. We further obtained structural and thermodynamic insights into the role of the D43A mutation, including perturbed protein structure near D43; weakened pairwise interactions of D43 with K42, N46, and R80; and decreased melting temperature and native fraction of the UDG D43A mutant compared with wild-type UDG.

Development of quantitative and continuous measure for severity degree of Alzheimer's disease evaluated from MRI images of 761 human brains.

BACKGROUND: Alzheimer's disease affects profoundly the quality of human behavior and cognition. The very broad distribution of its severity across various human subjects requires the quantitative diagnose of Alzheimer's disease beyond the conventional tripartite classification of cohorts such as cognitively normal (CN), mild cognitive impairment (MCI), Alzheimer's disease (AD). The unfolding of such broad distributions by the quantitative and continuous degree of AD severity is necessary for the precise diagnose in the cross-sectional study of different stages in AD. RESULTS: We conducted the massive reanalysis on MRI images of 761 human brains based on the accumulated bigdata of Alzheimer's Disease Neuroimaging Initiative. The score matrix of cortical thickness profile at cortex points of subjects was constructed by statistically learning the cortical thickness data of 761 human brains. We also developed a new and simple algebraic predictor which provides the quantitative and continuous degree of AD severity of subjects along the scale from 0 for fully CN to 1 for fully AD state. The mathematical measure of a new predictor for the degree of AD severity is presented based on a covariance correlation matrix of cortical thickness profile between human subjects. One can remove the uncertainty in the determination of different stages in AD by the quantitative degree of AD severity and thus go far beyond the tripartite classification of cohorts. CONCLUSIONS: We unfold the nature of broad distribution of AD severity of subjects even within a given cohort by the scale from 0 for fully CN to 1 for fully AD state. The quantitative and continuous degree of AD severity developed in this study would be a good practical measure for diagnosing the different stages in AD severity.

Structural and molecular insight into the pH-induced low-permeability of the voltage-gated potassium channel Kv1.2 through dewetting of the water cavity.

Understanding the gating mechanism of ion channel proteins is key to understanding the regulation of cell signaling through these channels. Channel opening and closing are regulated by diverse environmental factors that include temperature, electrical voltage across the channel, and proton concentration. Low permeability in voltage-gated potassium ion channels (Kv) is intimately correlated with the prolonged action potential duration observed in many acidosis diseases. The Kv channels consist of voltage-sensing domains (S1-S4 helices) and central pore domains (S5-S6 helices) that include a selectivity filter and water-filled cavity. The voltage-sensing domain is responsible for the voltage-gating of Kv channels. While the low permeability of Kv channels to potassium ion is highly correlated with the cellular proton concentration, it is unclear how an intracellular acidic condition drives their closure, which may indicate an additional pH-dependent gating mechanism of the Kv family. Here, we show that two residues E327 and H418 in the proximity of the water cavity of Kv1.2 play crucial roles as a pH switch. In addition, we present a structural and molecular concept of the pH-dependent gating of Kv1.2 in atomic detail, showing that the protonation of E327 and H418 disrupts the electrostatic balance around the S6 helices, which leads to a straightening transition in the shape of their axes and causes dewetting of the water-filled cavity and closure of the channel. Our work offers a conceptual advancement to the regulation of the pH-dependent gating of various voltage-gated ion channels and their related biological functions.

Uncovering multiloci-ordering by algebraic property of Laplacian matrix and its Fiedler vector.

MOTIVATION: The loci-ordering, based on two-point recombination fractions for a pair of loci, is the most important step in constructing a reliable and fine genetic map. RESULTS: Using the concept from complex graph theory, here we propose a Laplacian ordering approach which uncovers the loci-ordering of multiloci simultaneously. The algebraic property for a Fiedler vector of a Laplacian matrix, constructed from the recombination fraction of the loci-ordering for 26 loci of barley chromosome IV, 846 loci of Arabidopsis thaliana and 1903 loci of Malus domestica, together with the variable threshold uncovers their loci-orders. It offers an alternative yet robust approach for ordering multiloci. AVAILABILITY AND IMPLEMENTATION: Source code program with data set is available as supplementary data and also in a software category of the website (http://biophysics.dgist.ac.kr) CONTACT: crkim@pusan.ac.kr or iksoochang@dgist.ac.kr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structural Conversion of Aß17-42 Peptides from Disordered Oligomers to U-Shape Protofilaments via Multiple Kinetic Pathways.

Discovering the mechanisms by which proteins aggregate into fibrils is an essential first step in understanding the molecular level processes underlying neurodegenerative diseases such as Alzheimer's and Parkinson's. The goal of this work is to provide insights into the structural changes that characterize the kinetic pathways by which amyloid-ß peptides convert from monomers to oligomers to fibrils. By applying discontinuous molecular dynamics simulations to PRIME20, a force field designed to capture the chemical and physical aspects of protein aggregation, we have been able to trace out the entire aggregation process for a system containing 8 Aß17-42 peptides. We uncovered two fibrillization mechanisms that govern the structural conversion of Aß17-42 peptides from disordered oligomers into protofilaments. The first mechanism is monomeric conversion templated by a U-shape oligomeric nucleus into U-shape protofilament. The second mechanism involves a long-lived and on-pathway metastable oligomer with S-shape chains, having a C-terminal turn, en route to the final U-shape protofilament. Oligomers with this C-terminal turn have been regarded in recent experiments as a major contributing element to cell toxicity in Alzheimer's disease. The internal structures of the U-shape protofilaments from our PRIME20/DMD simulation agree well with those from solid state NMR experiments. The approach presented here offers a simple molecular-level framework to describe protein aggregation in general and to visualize the kinetic evolution of a putative toxic element in Alzheimer's disease in particular.

Molecular Basis of the Membrane Interaction of the ß2e Subunit of Voltage-Gated Ca(2+) Channels.

The auxiliary ß subunit plays an important role in the regulation of voltage-gated calcium (CaV) channels. Recently, it was revealed that ß2e associates with the plasma membrane through an electrostatic interaction between N-terminal basic residues and anionic phospholipids. However, a molecular-level understanding of ß-subunit membrane recruitment in structural detail has remained elusive. In this study, using a combination of site-directed mutagenesis, liposome-binding assays, and multiscale molecular-dynamics (MD) simulation, we developed a physical model of how the ß2e subunit is recruited electrostatically to the plasma membrane. In a fluorescence resonance energy transfer assay with liposomes, binding of the N-terminal peptide (23 residues) to liposome was significantly increased in the presence of phosphatidylserine (PS) and phosphatidylinositol 4,5-bisphosphate (PIP2). A mutagenesis analysis suggested that two basic residues proximal to Met-1, Lys-2 (K2) and Trp-5 (W5), are more important for membrane binding of the ß2e subunit than distal residues from the N-terminus. Our MD simulations revealed that a stretched binding mode of the N-terminus to PS is required for stable membrane attachment through polar and nonpolar interactions. This mode obtained from MD simulations is consistent with experimental results showing that K2A, W5A, and K2A/W5A mutants failed to be targeted to the plasma membrane. We also investigated the effects of a mutated ß2e subunit on inactivation kinetics and regulation of CaV channels by PIP2. In experiments with voltage-sensing phosphatase (VSP), a double mutation in the N-terminus of ß2e (K2A/W5A) increased the PIP2 sensitivity of CaV2.2 and CaV1.3 channels by ∼3-fold compared with wild-type ß2e subunit. Together, our results suggest that membrane targeting of the ß2e subunit is initiated from the nonspecific electrostatic insertion of N-terminal K2 and W5 residues into the membrane. The PS-ß2e interaction observed here provides a molecular insight into general principles for protein binding to the plasma membrane, as well as the regulatory roles of phospholipids in transporters and ion channels.

Impact of sequence on the molecular assembly of short amyloid peptides.

The goal of this work is to understand how the sequence of a protein affects the likelihood that it will form an amyloid fibril and the kinetics along the fibrillization pathway. The focus is on very short fragments of amyloid proteins since these play a role in the fibrillization of the parent protein and can form fibrils themselves. Discontinuous molecular dynamics simulations using the PRIME20 force field were performed of the aggregation of 48-peptide systems containing SNQNNF (PrP (170-175)), SSTSAA (RNaseA(15-20)), MVGGVV (Aß(35-40)), GGVVIA (Aß(37-42)), and MVGGVVIA (Aß(35-42)). In our simulations SNQQNF, SSTTSAA, and MVGGVV form large numbers of fibrillar structures spontaneously (as in experiment). GGVVIA forms ß-sheets that do not stack into fibrils (unlike experiment). The combination sequence MVGGVVIA forms less fibrils than MVGGVV, hindered by the presence of the hydrophobic residues at the C-terminal. Analysis of the simulation kinetics and energetics reveals why MVGGVV forms fibrils and GGVVIA does not, and why adding I and A to MVGGVVIA reduces fibrillization and enhances amorphous aggregation into oligomeric structures. The latter helps explain why Aß(1-42) assembles into more complex oligomers than Aß(1-40), a consequence of which is that it is more strongly associated with Alzheimer's disease.

Development of equation of motion deciphering locomotion including omega turns of Caenorhabditis elegans.

Locomotion is a fundamental behavior of Caenorhabditis elegans (C. elegans). Previous works on kinetic simulations of animals helped researchers understand the physical mechanisms of locomotion and the muscle-controlling principles of neuronal circuits as an actuator part. It has yet to be understood how C. elegans utilizes the frictional forces caused by the tension of its muscles to perform sequenced locomotive behaviors. Here, we present a two-dimensional rigid body chain model for the locomotion of C. elegans by developing Newtonian equations of motion for each body segment of C. elegans. Having accounted for friction-coefficients of the surrounding environment, elastic constants of C. elegans, and its kymogram from experiments, our kinetic model (ElegansBot) reproduced various locomotion of C. elegans such as, but not limited to, forward-backward-(omega turn)-forward locomotion constituting escaping behavior and delta-turn navigation. Additionally, ElegansBot precisely quantified the forces acting on each body segment of C. elegans to allow investigation of the force distribution. This model will facilitate our understanding of the detailed mechanism of various locomotive behaviors at any given friction-coefficients of the surrounding environment. Furthermore, as the model ensures the performance of realistic behavior, it can be used to research actuator-controller interaction between muscles and neuronal circuits.

Understanding the molecular mechanism of pathogenic variants of BIR2 domain in XIAP-deficient inflammatory bowel disease.

X-linked inhibitor of apoptosis protein (XIAP) deficiency causes refractory inflammatory bowel disease. The XIAP protein plays a pivotal role in the pro-inflammatory response through the nucleotide-binding oligomerization domain-containing signaling pathway that is important in mucosal homeostasis. We analyzed the molecular mechanism of non-synonymous pathogenic variants (PVs) of XIAP BIR2 domain. We generated N-terminally green fluorescent protein-tagged XIAP constructs of representative non-synonymous PVs. Co-immunoprecipitation and fluorescence cross-correlation spectroscopy showed that wild-type XIAP and RIP2 preferentially interacted in live cells, whereas all non-synonymous PV XIAPs failed to interact properly with RIP2. Structural analysis showed that various structural changes by mutations, such as hydrophobic core collapse, Zn-finger loss, and spatial rearrangement, destabilized the two loop structures (174-182 and 205-215) that critically interact with RIP2. Subsequently, it caused a failure of RIP2 ubiquitination and loss of protein deficiency by the auto-ubiquitination of all XIAP mutants. These findings could enhance our understanding of the role of XIAP mutations in XIAP-deficient inflammatory bowel disease and may benefit future therapeutic strategies.

Functional implications of an intermeshing cogwheel-like interaction between TolC and MacA in the action of macrolide-specific efflux pump MacAB-TolC.

Macrolide-specific efflux pump MacAB-TolC has been identified in diverse gram-negative bacteria including Escherichia coli. The inner membrane transporter MacB requires the outer membrane factor TolC and the periplasmic adaptor protein MacA to form a functional tripartite complex. In this study, we used a chimeric protein containing the tip region of the TolC α-barrel to investigate the role of the TolC α-barrel tip region with regard to its interaction with MacA. The chimeric protein formed a stable complex with MacA, and the complex formation was abolished by substitution at the functionally essential residues located at the MacA α-helical tip region. Electron microscopic study delineated that this complex was made by tip-to-tip interaction between the tip regions of the α-barrels of TolC and MacA, which correlated well with the TolC and MacA complex calculated by molecular dynamics. Taken together, our results demonstrate that the MacA hexamer interacts with TolC in a tip-to-tip manner, and implies the manner by which MacA induces opening of the TolC channel.

Crystal structure of Fushi tarazu factor 1 ligand binding domain/Fushi tarazu peptide complex identifies new class of nuclear receptors.

The interaction between the orphan nuclear receptor FTZ-F1 (Fushi tarazu factor 1) and the segmentation gene protein FTZ is critical for specifying alternate parasegments in the Drosophila embryo. Here, we have determined the structure of the FTZ-F1 ligand-binding domain (LBD)·FTZ peptide complex using x-ray crystallography. Strikingly, the ligand-binding pocket of the FTZ-F1 LBD is completely occupied by helix 6 (H6) of the receptor, whereas the cofactor FTZ binds the co-activator cleft site of the FTZ-F1 LBD. Our findings suggest that H6 is essential for transcriptional activity of FTZ-F1; this is further supported by data from mutagenesis and activity assays. These data suggest that FTZ-F1 might belong to a novel class of ligand-independent nuclear receptors. Our findings are intriguing given that the highly homologous human steroidogenic factor-1 and liver receptor homolog-1 LBDs exhibit sizable ligand-binding pockets occupied by putative ligand molecules.

PACSY, a relational database management system for protein structure and chemical shift analysis.

PACSY (Protein structure And Chemical Shift NMR spectroscopY) is a relational database management system that integrates information from the Protein Data Bank, the Biological Magnetic Resonance Data Bank, and the Structural Classification of Proteins database. PACSY provides three-dimensional coordinates and chemical shifts of atoms along with derived information such as torsion angles, solvent accessible surface areas, and hydrophobicity scales. PACSY consists of six relational table types linked to one another for coherence by key identification numbers. Database queries are enabled by advanced search functions supported by an RDBMS server such as MySQL or PostgreSQL. PACSY enables users to search for combinations of information from different database sources in support of their research. Two software packages, PACSY Maker for database creation and PACSY Analyzer for database analysis, are available from http://pacsy.nmrfam.wisc.edu.

Spontaneous formation of twisted Aß(16-22) fibrils in large-scale molecular-dynamics simulations.

Protein aggregation is associated with fatal neurodegenerative diseases, including Alzheimer's and Parkinson's. Mapping out kinetics along the aggregation pathway could provide valuable insights into the mechanisms that drive oligomerization and fibrillization, but that is beyond the current scope of computational research. Here we trace out the full kinetics of the spontaneous formation of fibrils by 48 Aß(16-22) peptides, following the trajectories in molecular detail from an initial random configuration to a final configuration of twisted protofilaments with cross-ß-structure. We accomplish this by performing large-scale molecular-dynamics simulations based on an implicit-solvent, intermediate-resolution protein model, PRIME20. Structural details such as the intersheet distance, perfectly antiparallel ß-strands, and interdigitating side chains analogous to a steric zipper interface are explained by and in agreement with experiment. Two characteristic fibrillization mechanisms - nucleation/templated growth and oligomeric merging/structural rearrangement - emerge depending on the temperature.

Computer simulation study of amyloid fibril formation by palindromic sequences in prion peptides.

We simulate the aggregation of large systems containing palindromic peptides from the Syrian hamster prion protein SHaPrP 113-120 (AGAAAAGA) and the mouse prion protein MoPrP 111-120 (VAGAAAAGAV) and eight sequence variations: GAAAAAAG, (AG)(4) , A8, GAAAGAAA, A10, V10, GAVAAAAVAG, and VAVAAAAVAV The first two peptides are thought to act as the Velcro that holds the parent prion proteins together in amyloid structures and can form fibrils themselves. Kinetic events along the fibrillization pathway influence the types of structures that occur and variations in the sequence affect aggregation kinetics and fibrillar structure. Discontinuous molecular dynamics simulations using the PRIME20 force field are performed on systems containing 48 peptides starting from a random coil configuration. Depending on the sequence, fibrillar structures form spontaneously over a range of temperatures, below which amorphous aggregates form and above which no aggregation occurs. AGAAAAGA forms well organized fibrillar structures whereas VAGAAAAGAV forms less well organized structures that are partially fibrillar and partially amorphous. The degree of order in the fibrillar structure stems in part from the types of kinetic events leading up to its formation, with AGAAAAGA forming less amorphous structures early in the simulation than VAGAAAAGAV. The ability to form fibrils increases as the chain length and the length of the stretch of hydrophobic residues increase. However as the hydrophobicity of the sequence increases, the ability to form well-ordered structures decreases. Thus, longer hydrophobic sequences form slightly disordered aggregates that are partially fibrillar and partially amorphous. Subtle changes in sequence result in slightly different fibril structures.

Uncovering symmetry-breaking vector and reliability order for assigning secondary structures of proteins from atomic NMR chemical shifts in amino acids.

Unravelling the complex correlation between chemical shifts of (13) C (α), (13) C (ß), (13) C', (1) H (α), (15) N, (1) H ( N ) atoms in amino acids of proteins from NMR experiment and local structural environments of amino acids facilitates the assignment of secondary structures of proteins. This is an important impetus for both determining the three-dimensional structure and understanding the biological function of proteins. The previous empirical correlation scores which relate chemical shifts of (13) C (α), (13) C (ß), (13) C', (1) H (α), (15) N, (1) H ( N ) atoms to secondary structures resulted in progresses toward assigning secondary structures of proteins. However, the physical-mathematical framework for these was elusive partly due to both the limited and orthogonal exploration of higher-dimensional chemical shifts of hetero-nucleus and the lack of physical-mathematical understanding underlying those correlation scores. Here we present a simple multi-dimensional hetero-nuclear chemical shift score function (MDHN-CSSF) which captures systematically the salient feature of such complex correlations without any references to a random coil state of proteins. We uncover the symmetry-breaking vector and its reliability order not only for distinguishing different secondary structures of proteins but also for capturing the delicate sensitivity interplayed among chemical shifts of (13) C (α), (13) C (ß), (13) C', (1) H (α), (15) N, (1) H ( N ) atoms simultaneously, which then provides a straightforward framework toward assigning secondary structures of proteins. MDHN-CSSF could correctly assign secondary structures of training (validating) proteins with the favourable (comparable) Q3 scores in comparison with those from the previous correlation scores. MDHN-CSSF provides a simple and robust strategy for the systematic assignment of secondary structures of proteins and would facilitate the de novo determination of three-dimensional structures of proteins.

A simple and exact Laplacian clustering of complex networking phenomena: application to gene expression profiles.

Unraveling of the unified networking characteristics of complex networking phenomena is of great interest yet a formidable task. There is currently no simple strategy with a rigorous framework. Using an analogy to the exact algebraic property for a transition matrix of a master equation in statistical physics, we propose a method based on a Laplacian matrix for the discovery and prediction of new classes in the unsupervised complex networking phenomena where the class of each sample is completely unknown. Using this proposed Laplacian approach, we can simultaneously discover different classes and determine the identity of each class. Through an illustrative test of the Laplacian approach applied to real datasets of gene expression profiles, leukemia data [Golub TR, et al. (1999) Science 286:531-537], and lymphoma data [Alizadeh AA, et al. (2000) Nature 403:503-511], we demonstrate that this approach is accurate and robust with a mathematical and physical realization. It offers a general framework for characterizing any kind of complex networking phenomenon in broad areas irrespective of whether they are supervised or unsupervised.

Cooperative folding kinetics of BBL protein and peripheral subunit-binding domain homologues.

Recent experiments claiming that Naf-BBL protein follows a global downhill folding raised an important controversy as to the folding mechanism of fast-folding proteins. Under the global downhill folding scenario, not only do proteins undergo a gradual folding, but folding events along the continuous folding pathway also could be mapped out from the equilibrium denaturation experiment. Based on the exact calculation using a free energy landscape, relaxation eigenmodes from a master equation, and Monte Carlo simulation of an extended Muñoz-Eaton model that incorporates multiscale-heterogeneous pairwise interactions between amino acids, here we show that the very nature of a two-state cooperative transition such as a bimodal distribution from an exact free energy landscape and biphasic relaxation kinetics manifest in the thermodynamics and folding-unfolding kinetics of BBL and peripheral subunit-binding domain homologues. Our results provide an unequivocal resolution to the fundamental controversy related to the global downhill folding scheme, whose applicability to other proteins should be critically reexamined.

Impaired formation of high-order gephyrin oligomers underlies gephyrin dysfunction-associated pathologies.

Gephyrin is critical for the structure, function, and plasticity of inhibitory synapses. Gephyrin mutations have been linked to various neurological disorders; however, systematic analyses of the functional consequences of these mutations are lacking. Here, we performed molecular dynamics simulations of gephyrin to predict how six reported point mutations might change the structural stability and/or function of gephyrin. Additional in silico analyses revealed that the A91T and G375D mutations reduce the binding free energy of gephyrin oligomer formation. Gephyrin A91T and G375D displayed altered clustering patterns in COS-7 cells and nullified the inhibitory synapse-promoting effect of gephyrin in cultured neurons. However, only the G375D mutation reduced gephyrin interaction with GABAA receptors and neuroligin-2 in mouse brain; it also failed to normalize deficits in GABAergic synapse maintenance and neuronal hyperactivity observed in hippocampal dentate gyrus-specific gephyrin-deficient mice. Our results provide insights into biochemical, cell-biological, and network-activity effects of the pathogenic G375D mutation.

Interpretation of XIAP Variants of Uncertain Significance in Paediatric Patients with Refractory Crohn's Disease.

BACKGROUND AND AIMS: Mutations in XIAP can lead to the development of treatment-refractory severe paediatric Crohn's disease [CD], for which haematopoietic stem cell transplantation is the primary therapeutic option. The interpretation of variants of uncertain significance [VUSs] in XIAP needs to be scrutinized. METHODS: Targeted next-generation sequencing was performed for 33 male paediatric patients with refractory CD admitted at a tertiary referral hospital. To obtain functional data, biomolecular cell assays and supercomputing molecular dynamics simulations were performed. RESULTS: Nine unrelated male patients harboured hemizygous XIAP variants. Four known pathogenic variants and one novel pathogenic variant [p.Lys168Serfs*12] were identified in five patients, and two novel VUSs [p.Gly205del and p.Pro260Ser] and one known VUS [p.Glu350del] were identified in the remaining four. Among children with VUSs, only the subject with p.Gly205del exhibited defective NOD2 signalling. Using molecular dynamics simulation, we determined that the altered backbone torsional energy of C203 in XIAP of p.G205del was ~2 kcal/mol, suggesting loss of zinc binding in the mutant XIAP protein and poor coordination between the mutant XIAP and RIP2 proteins. Elevated auto-ubiquitination of zinc-depleted p.G205del XIAP protein resulted in XIAP protein deficiency. CONCLUSION: A high prevalence of XIAP deficiency was noted among children with refractory CD. Advanced functional studies decreased the subjectivity in the case-level interpretation of XIAP VUSs and directed consideration of haematopoietic stem cell transplantation.

Extending the PRIME model for protein aggregation to all 20 amino acids.

We extend PRIME, an intermediate-resolution protein model previously used in simulations of the aggregation of polyalanine and polyglutamine, to the description of the geometry and energetics of peptides containing all 20 amino acid residues. The 20 amino acid side chains are classified into 14 groups according to their hydrophobicity, polarity, size, charge, and potential for side chain hydrogen bonding. The parameters for extended PRIME, called PRIME 20, include hydrogen-bonding energies, side chain interaction range and energy, and excluded volume. The parameters are obtained by applying a perceptron-learning algorithm and a modified stochastic learning algorithm that optimizes the energy gap between 711 known native states from the PDB and decoy structures generated by gapless threading. The number of independent pair interaction parameters is chosen to be small enough to be physically meaningful yet large enough to give reasonably accurate results in discriminating decoys from native structures. The most physically meaningful results are obtained with 19 energy parameters.