Selection among site-dependent structurally constrained substitution models of protein evolution by approximate Bayesian computation.

MOTIVATION: The selection among substitution models of molecular evolution is fundamental for obtaining accurate phylogenetic inferences. At the protein level, evolutionary analyses are traditionally based on empirical substitution models but these models make unrealistic assumptions and are being surpassed by structurally constrained substitution (SCS) models. The SCS models often consider site-dependent evolution, a process that provides realism but complicates their implementation into likelihood functions that are commonly used for substitution model selection. RESULTS: We present a method to perform selection among site-dependent SCS models, also among empirical and site-dependent SCS models, based on the approximate Bayesian computation (ABC) approach and its implementation into the computational framework ProteinModelerABC. The framework implements ABC with and without regression adjustments and includes diverse empirical and site-dependent SCS models of protein evolution. Using extensive simulated data, we found that it provides selection among SCS and empirical models with acceptable accuracy. As illustrative examples, we applied the framework to analyze a variety of protein families observing that SCS models fit them better than the corresponding best-fitting empirical substitution models. AVAILABILITY AND IMPLEMENTATION: ProteinModelerABC is freely available from https://github.com/DavidFerreiro/ProteinModelerABC, can run in parallel and includes a graphical user interface. The framework is distributed with detailed documentation and ready-to-use examples.

Substitution Models of Protein Evolution with Selection on Enzymatic Activity.

Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme-substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.

The evolution of the HIV-1 protease folding stability.

The evolution of structural proteins is generally constrained by the folding stability. However, little is known about the particular capacity of viral proteins to accommodate mutations that can potentially affect the protein stability and, in general, the evolution of the protein stability over time. As an illustrative model case, here, we investigated the evolution of the stability of the human immunodeficiency virus (HIV-1) protease (PR), which is a common HIV-1 drug target, under diverse evolutionary scenarios that include (1) intra-host virus evolution in a cohort of seventy-five patients sampled over time, (2) intra-host virus evolution sampled before and after specific PR-based treatments, and (3) inter-host evolution considering extant and ancestral (reconstructed) PR sequences from diverse HIV-1 subtypes. We also investigated the specific influence of currently known HIV-1 PR resistance mutations on the PR folding stability. We found that the HIV-1 PR stability fluctuated over time within a constant and wide range in any studied evolutionary scenario, accommodating multiple mutations that partially affected the stability while maintaining activity. We did not identify relationships between change of PR stability and diverse clinical parameters such as viral load, CD4+ T-cell counts, and a surrogate of time from infection. Counterintuitively, we predicted that nearly half of the studied HIV-1 PR resistance mutations do not significantly decrease stability, which, together with compensatory mutations, would allow the protein to adapt without requiring dramatic stability changes. We conclude that the HIV-1 PR presents a wide structural plasticity to acquire molecular adaptations without affecting the overall evolution of stability.

Evaluating Causes of Current Genetic Gradients of Modern Humans of the Iberian Peninsula.

The history of modern humans in the Iberian Peninsula includes a variety of population arrivals sometimes presenting admixture with resident populations. Genetic data from current Iberian populations revealed an overall east-west genetic gradient that some authors interpreted as a direct consequence of the Reconquista, where Catholic Kingdoms expanded their territories toward the south while displacing Muslims. However, this interpretation has not been formally evaluated. Here, we present a qualitative analysis of the causes of the current genetic gradient observed in the Iberian Peninsula using extensive spatially explicit computer simulations based on a variety of evolutionary scenarios. Our results indicate that the Neolithic range expansion clearly produces the orientation of the observed genetic gradient. Concerning the Reconquista (including political borders among Catholic Kingdoms and regions with different languages), if modeled upon a previous Neolithic expansion, it effectively favored the orientation of the observed genetic gradient and shows local isolation of certain regions (i.e., Basques and Galicia). Despite additional evolutionary scenarios could be evaluated to more accurately decipher the causes of the Iberian genetic gradient, here we show that this gradient has a more complex explanation than that previously hypothesized.

Multiple Drug Treatments That Increase cAMP Signaling Restore Long-Term Memory and Aberrant Signaling in Fragile X Syndrome Models.

Fragile X is the most common monogenic disorder associated with intellectual disability (ID) and autism spectrum disorders (ASD). Additionally, many patients are afflicted with executive dysfunction, ADHD, seizure disorder and sleep disturbances. Fragile X is caused by loss of FMRP expression, which is encoded by the FMR1 gene. Both the fly and mouse models of fragile X are also based on having no functional protein expression of their respective FMR1 homologs. The fly model displays well defined cognitive impairments and structural brain defects and the mouse model, although having subtle behavioral defects, has robust electrophysiological phenotypes and provides a tool to do extensive biochemical analysis of select brain regions. Decreased cAMP signaling has been observed in samples from the fly and mouse models of fragile X as well as in samples derived from human patients. Indeed, we have previously demonstrated that strategies that increase cAMP signaling can rescue short term memory in the fly model and restore DHPG induced mGluR mediated long term depression (LTD) in the hippocampus to proper levels in the mouse model (McBride et al., 2005; Choi et al., 2011, 2015). Here, we demonstrate that the same three strategies used previously with the potential to be used clinically, lithium treatment, PDE-4 inhibitor treatment or mGluR antagonist treatment can rescue long term memory in the fly model and alter the cAMP signaling pathway in the hippocampus of the mouse model.

The Drosophila DmGluRA is required for social interaction and memory.

Metabotropic glutamate receptors (mGluRs) have well-established roles in cognition and social behavior in mammals. Whether or not these roles have been conserved throughout evolution from invertebrate species is less clear. Mammals have eight mGluRs whereas Drosophila has a single DmGluRA, which has both Gi and Gq coupled signaling activity. We have utilized Drosophila to examine the role of DmGluRA in social behavior and various phases of memory. We have found that flies that are homozygous or heterozygous for loss of function mutations of DmGluRA have impaired social behavior in male Drosophila. Futhermore, flies that are heterozygous for loss of function mutations of DmGluRA have impaired learning during training, immediate-recall memory, short-term memory, and long-term memory as young adults. This work demonstrates a role for mGluR activity in both social behavior and memory in Drosophila.

Pharmacological and genetic reversal of age-dependent cognitive deficits attributable to decreased presenilin function.

Alzheimer's disease (AD) is the leading cause of cognitive loss and neurodegeneration in the developed world. Although its genetic and environmental causes are not generally known, familial forms of the disease (FAD) are attributable to mutations in a single copy of the Presenilin (PS) and amyloid precursor protein genes. The dominant inheritance pattern of FAD indicates that it may be attributable to gain or change of function mutations. Studies of FAD-linked forms of presenilin (psn) in model organisms, however, indicate that they are loss of function, leading to the possibility that a reduction in PS activity might contribute to FAD and that proper psn levels are important for maintaining normal cognition throughout life. To explore this issue further, we have tested the effect of reducing psn activity during aging in Drosophila melanogaster males. We have found that flies in which the dosage of psn function is reduced by 50% display age-onset impairments in learning and memory. Treatment with metabotropic glutamate receptor (mGluR) antagonists or lithium during the aging process prevented the onset of these deficits, and treatment of aged flies reversed the age-dependent deficits. Genetic reduction of Drosophila metabotropic glutamate receptor (DmGluRA), the inositol trisphosphate receptor (InsP(3)R), or inositol polyphosphate 1-phosphatase also prevented these age-onset cognitive deficits. These findings suggest that reduced psn activity may contribute to the age-onset cognitive loss observed with FAD. They also indicate that enhanced mGluR signaling and calcium release regulated by InsP(3)R as underlying causes of the age-dependent cognitive phenotypes observed when psn activity is reduced.

Age-dependent cognitive impairment in a Drosophila fragile X model and its pharmacological rescue.

Fragile X syndrome afflicts 1 in 2,500 individuals and is the leading heritable cause of mental retardation worldwide. The overriding clinical manifestation of this disease is mild to severe cognitive impairment. Age-dependent cognitive decline has been identified in Fragile X patients, although it has not been fully characterized nor examined in animal models. A Drosophila model of this disease has been shown to display phenotypes bearing similarity to Fragile X symptoms. Most notably, we previously identified naive courtship and memory deficits in young adults with this model that appear to be due to enhanced metabotropic glutamate receptor (mGluR) signaling. Herein we have examined age-related cognitive decline in the Drosophila Fragile X model and found an age-dependent loss of learning during training. We demonstrate that treatment with mGluR antagonists or lithium can prevent this age-dependent cognitive impairment. We also show that treatment with mGluR antagonists or lithium during development alone displays differential efficacy in its ability to rescue naive courtship, learning during training and memory in aged flies. Furthermore, we show that continuous treatment during aging effectively rescues all of these phenotypes. These results indicate that the Drosophila model recapitulates the age-dependent cognitive decline observed in humans. This places Fragile X in a category with several other diseases that result in age-dependent cognitive decline. This demonstrates a role for the Drosophila Fragile X Mental Retardation Protein (dFMR1) in neuronal physiology with regard to cognition during the aging process. Our results indicate that misregulation of mGluR activity may be causative of this age onset decline and strengthens the possibility that mGluR antagonists and lithium may be potential pharmacologic compounds for counteracting several Fragile X symptoms.

Pharmacological rescue of synaptic plasticity, courtship behavior, and mushroom body defects in a Drosophila model of fragile X syndrome.

Fragile X syndrome is a leading heritable cause of mental retardation that results from the loss of FMR1 gene function. A Drosophila model for Fragile X syndrome, based on the loss of dfmr1 activity, exhibits phenotypes that bear similarity to Fragile X-related symptoms. Herein, we demonstrate that treatment with metabotropic glutamate receptor (mGluR) antagonists or lithium can rescue courtship and mushroom body defects observed in these flies. Furthermore, we demonstrate that dfmr1 mutants display cognitive deficits in experience-dependent modification of courtship behavior, and treatment with mGluR antagonists or lithium restores these memory defects. These findings implicate enhanced mGluR signaling as the underlying cause of the cognitive, as well as some of the behavioral and neuronal, phenotypes observed in the Drosophila Fragile X model. They also raise the possibility that compounds having similar effects on metabotropic glutamate receptors may ameliorate cognitive and behavioral defects observed in Fragile X patients.