Analysis of conformational motions and related key residue interactions responsible for a specific function of proteins with elastic network model

Protein collective motions play a critical role in many biochemical processes. How to predict the functional motions and the related key residue interactions in proteins is important for our understanding in the mechanism of the biochemical processes. Normal mode analysis (NMA) of the elastic network model (ENM) is one of the effective approaches to investigate the structure-encoded motions in proteins. However, the motion modes revealed by the conventional NMA approach do not necessarily correspond to a specific function of protein. In the present work, a new analysis method was proposed to identify the motion modes responsible for a specific function of proteins and then predict the key residue interactions involved in the functional motions by using a perturbation approach. In our method, an internal coordinate that accounts for the specific function was introduced, and the Cartesian coordinate space was transformed into the internal/Cartesian space by using linear approximation, where the introduced internal coordinate serves as one of the axes of the coordinate space. NMA of ENM in this internal/Cartesian space was performed and the function-relevant motion modes were identified according to their contributions to the specific function of proteins. Then the key residue interactions important for the functional motions of the protein were predicted as the interactions whose perturbation largely influences the fluctuation along the internal coordinate. Using our proposed methods, the maltose transporter (MalFGK₂) from E. Coli was studied. The functional motions and the key residue interactions that are related to the channel-gating function of this protein were successfully identified.     

Orexinergic Neurotransmission in Temperature Responses to Methamphetamine and Stress: Mathematical Modeling as a Data Assimilation Approach

Experimental Data

Orexinergic neurotransmission is involved in mediating temperature responses to methamphetamine (Meth). In experiments in rats, SB-334867 (SB), an antagonist of orexin receptors (OX1R), at a dose of 10 mg/kg decreases late temperature responses (t>60 min) to an intermediate dose of Meth (5 mg/kg). A higher dose of SB (30 mg/kg) attenuates temperature responses to low dose (1 mg/kg) of Meth and to stress. In contrast, it significantly exaggerates early responses (t<60 min) to intermediate and high doses (5 and 10 mg/kg) of Meth. As pretreatment with SB also inhibits temperature response to the stress of injection, traditional statistical analysis of temperature responses is difficult.

Mathematical Modeling

We have developed a mathematical model that explains the complexity of temperature responses to Meth as the interplay between excitatory and inhibitory nodes. We have extended the developed model to include the stress of manipulations and the effects of SB. Stress is synergistic with Meth on the action on excitatory node. Orexin receptors mediate an activation of on both excitatory and inhibitory nodes by low doses of Meth, but not on the node activated by high doses (HD). Exaggeration of early responses to high doses of Meth involves disinhibition: low dose of SB decreases tonic inhibition of HD and lowers the activation threshold, while the higher dose suppresses the inhibitory component. Using a modeling approach to data assimilation appears efficient in separating individual components of complex response with statistical analysis unachievable by traditional data processing methods.     

Identifying a Kinase Network Regulating FGF14:Nav1.6 Complex Assembly Using Split-Luciferase Complementation

Kinases play fundamental roles in the brain. Through complex signaling pathways, kinases regulate the strength of protein:protein interactions (PPI) influencing cell cycle, signal transduction, and electrical activity of neurons. Changes induced by kinases on neuronal excitability, synaptic plasticity and brain connectivity are linked to complex brain disorders, but the molecular mechanisms underlying these cellular events remain for the most part elusive. To further our understanding of brain disease, new methods for rapidly surveying kinase pathways in the cellular context are needed. The bioluminescence-based luciferase complementation assay (LCA) is a powerful, versatile toolkit for the exploration of PPI. LCA relies on the complementation of two firefly luciferase protein fragments that are functionally reconstituted into the full luciferase enzyme by two interacting binding partners. Here, we applied LCA in live cells to assay 12 kinase pathways as regulators of the PPI complex formed by the voltage-gated sodium channel, Nav1.6, a transmembrane ion channel that elicits the action potential in neurons and mediates synaptic transmission, and its multivalent accessory protein, the fibroblast growth factor 14 (FGF14). Through extensive dose-dependent validations of structurally-diverse kinase inhibitors and hierarchical clustering, we identified the PI3K/Akt pathway, the cell-cycle regulator Wee1 kinase, and protein kinase C (PKC) as prospective regulatory nodes of neuronal excitability through modulation of the FGF14:Nav1.6 complex. Ingenuity Pathway Analysis shows convergence of these pathways on glycogen synthase kinase 3 (GSK3) and functional assays demonstrate that inhibition of GSK3 impairs excitability of hippocampal neurons. This combined approach provides a versatile toolkit for rapidly surveying PPI signaling, allowing the discovery of new modular pathways centered on GSK3 that might be the basis for functional alterations between the normal and diseased brain.     

Modelling Infectious Hematopoietic Necrosis Virus Dispersion from Marine Salmon Farms in the Discovery Islands,British Columbia,Canada

Finite volume ocean circulation and particle tracking models are used to simulate water-borne transmission of infectious hematopoietic necrosis virus (IHNV) among Atlantic salmon (Salmo salar) farms in the Discovery Islands region of British Columbia, Canada. Historical simulations for April and July 2010 are carried out to demonstrate the seasonal impact of river discharge, wind, ultra-violet (UV) radiation, and heat flux conditions on near-surface currents, viral dispersion and survival. Numerical particles released from infected farm fish in accordance with IHNV shedding rates estimated through laboratory experiments are dispersed by model oceanic flows. Viral particles are inactivated by ambient UV radiation levels and by the natural microbial community at rates derived through laboratory studies. Viral concentration maps showing temporal and spatial changes are produced and combined with lab-determined minimum infectious dosages to estimate the infective connectivity among farms. Results demonstrate that neighbouring naïve farms can become exposed to IHNV via water-borne transport from an IHNV diseased farm, with a higher risk in April than July, and that many events in the sequence of farm outbreaks in 2001-2002 are consistent with higher risks in our farm connectivity matrix. Applications to other diseases, transfers between farmed and wild fish, and the effect of vaccinations are also discussed.     

Retinol kinetics in unsupplemented and vitamin A-retinoic acid supplemented neonatal rats: a preliminary model

Vitamin A (VA) metabolism in neonates is virtually uncharacterized. Our objective was to develop a compartmental model of VA metabolism in unsupplemented and VA-supplemented neonatal rats. On postnatal day 4, pups (n = 3/time) received 11,12-[³H]retinol orally, in either oil (control) or VA combined with retinoic acid (VARA) [VA (∼6 mg/kg body weight) + 10% retinoic acid]. Plasma and tissues were collected at 14 time points up to 14 days after dose administration. VARA supplementation rapidly, but transiently, increased total retinol mass in plasma, liver, and lung. It decreased the peak fraction of the dose in plasma. A multi-compartmental model developed to fit plasma [³H]retinol data predicted more extensive recycling of retinol between plasma and tissues in neonates compared with that reported in adults (144 vs. 12–13 times). In VARA pups, the recycling number for retinol between plasma and tissues (100 times) and the time that retinol spent in plasma were both lower compared with controls; VARA also stimulated the uptake of plasma VA into extravascular tissues. A VARA perturbation model indicated that the effect of VARA in stimulating VA uptake into tissues in neonates is both dramatic and transient.     

Survival rates of mutant genes under artificial selection using individual and family information

We have used diffusion and branching process methods to investigate fixation rates, probabilities of survival per generation, and times to fixation of mutant genes under different selection methods incorporating individual and family information. Diffusion approximations fit well to simulated results even for large selection coefficients. Methods that give much weight to family information, such as BLUP evaluation which is widely used in animal breeding, reduce fixation rates of mutant genes because of the reduced effective population sizes. In general, it is observed that even mutants with relatively small heterozygous effects (say 0.1 phenotypic standard deviation) are practically ‘safe’ (i.e. their probability of loss from one generation to the next is smaller than, say, 10%) after just a few generations, typically less than 10. For methods of selection with larger effective size, such as within-family selection, the mutant is ‘safe’ in the population somewhat earlier but eventual fixation takes a longer time. Finally we evaluate the amount by which the use of marker assisted selection reduces the fixation probability of newly arisen mutants.