RESUMEN
RAS is a signaling protein associated with the cell membrane that is mutated in up to 30% of human cancers. RAS signaling has been proposed to be regulated by dynamic heterogeneity of the cell membrane. Investigating such a mechanism requires near-atomistic detail at macroscopic temporal and spatial scales, which is not possible with conventional computational or experimental techniques. We demonstrate here a multiscale simulation infrastructure that uses machine learning to create a scale-bridging ensemble of over 100,000 simulations of active wild-type KRAS on a complex, asymmetric membrane. Initialized and validated with experimental data (including a new structure of active wild-type KRAS), these simulations represent a substantial advance in the ability to characterize RAS-membrane biology. We report distinctive patterns of local lipid composition that correlate with interfacially promiscuous RAS multimerization. These lipid fingerprints are coupled to RAS dynamics, predicted to influence effector binding, and therefore may be a mechanism for regulating cell signaling cascades.
Asunto(s)
Membrana Celular/enzimología , Lípidos/química , Aprendizaje Automático , Simulación de Dinámica Molecular , Multimerización de Proteína , Proteínas Proto-Oncogénicas p21(ras)/química , Transducción de Señal , HumanosRESUMEN
The small GTPase KRAS is localized at the plasma membrane where it functions as a molecular switch, coupling extracellular growth factor stimulation to intracellular signaling networks. In this process, KRAS recruits effectors, such as RAF kinase, to the plasma membrane where they are activated by a series of complex molecular steps. Defining the membrane-bound state of KRAS is fundamental to understanding the activation of RAF kinase and in evaluating novel therapeutic opportunities for the inhibition of oncogenic KRAS-mediated signaling. We combined multiple biophysical measurements and computational methodologies to generate a consensus model for authentically processed, membrane-anchored KRAS. In contrast to the two membrane-proximal conformations previously reported, we identify a third significantly populated state using a combination of neutron reflectivity, fast photochemical oxidation of proteins (FPOP), and NMR. In this highly populated state, which we refer to as "membrane-distal" and estimate to comprise â¼90% of the ensemble, the G-domain does not directly contact the membrane but is tethered via its C-terminal hypervariable region and carboxymethylated farnesyl moiety, as shown by FPOP. Subsequent interaction of the RAF1 RAS binding domain with KRAS does not significantly change G-domain configurations on the membrane but affects their relative populations. Overall, our results are consistent with a directional fly-casting mechanism for KRAS, in which the membrane-distal state of the G-domain can effectively recruit RAF kinase from the cytoplasm for activation at the membrane.
Asunto(s)
Proteínas Proto-Oncogénicas p21(ras)/metabolismo , Quinasas raf/metabolismo , Membrana Celular/metabolismo , Simulación de Dinámica MolecularRESUMEN
Cyclooxygenase-2 [(COX-2) or prostaglandin endoperoxide H2 synthase-2 (PTGS-2)] induces the production of prostaglandins as part of the host-immune response to infections. Although a number of studies have demonstrated the effects of COX-2 promoter variants on autoimmune and inflammatory diseases, their role in malaria remains undefined. As such, we investigated the relationship between four COX-2 promoter variants (COX-2 -512 C > T, -608 T > C, -765 G > C, and -1195 A > G) and susceptibility to malaria and severe malarial anemia (SMA) upon enrollment and longitudinally over a 36-month follow-up period. All-cause mortality was also explored. The investigation was carried out in children (n = 1081, age; 2-70 months) residing in a holoendemic Plasmodium falciparum transmission region of western Kenya. At enrollment, genotypes/haplotypes (controlling for anemia-promoting covariates) did not reveal any strong effects on susceptibility to either malaria or SMA. Longitudinal analyses showed decreased malaria episodes in children who inherited the -608 CC mutant allele (RR = 0.746, P = 1.811 × 10-4) and -512C/-608T/-765G/-1195G (CTGG) haplotype (RR = 0.856, P = 0.011), and increased risk in TTCA haplotype carriers (RR = 1.115, P = 0.026). Over the follow-up period, inheritance of the rare TTCG haplotype was associated with enhanced susceptibility to both malaria (RR = 1.608, P = 0.016) and SMA (RR = 5.714, P = 0.004), while carriage of the rare TTGG haplotype increased the risk of malaria (RR = 1.755, P = 0.007), SMA (RR = 8.706, P = 3.97 × 10-4), and all-cause mortality (HR = 110.000, P = 0.001). Collectively, these results show that SNP variations in the COX-2 promoter, and their inherited combinations, are associated with the longitudinal risk of malaria, SMA, and all-cause mortality among children living in a high transmission area for P. falciparum.
Asunto(s)
Anemia/genética , Ciclooxigenasa 2/genética , Malaria Falciparum/genética , Malaria/genética , Polimorfismo de Nucleótido Simple/genética , Regiones Promotoras Genéticas/genética , Anemia/mortalidad , Niño , Preescolar , Femenino , Predisposición Genética a la Enfermedad , Genotipo , Haplotipos , Humanos , Lactante , Kenia , Estudios Longitudinales , Malaria/inmunología , Malaria/mortalidad , Malaria/transmisión , Malaria Falciparum/inmunología , Malaria Falciparum/mortalidad , Malaria Falciparum/transmisión , Masculino , RiesgoRESUMEN
Drug discovery faces a crisis. The industry has used up the "obvious" space in which to find novel drugs for biomedical applications, and productivity is declining. One strategy to combat this is rational approaches to expand the search space without relying on chemical intuition, to avoid rediscovery of similar spaces. In this work, we present proof of concept of an approach to rationally identify a "chemical vocabulary" related to a specific drug activity of interest without employing known rules. We focus on the pressing concern of multidrug resistance in Pseudomonas aeruginosa by searching for submolecules that promote compound entry into this bacterium. By synergizing theory, computation, and experiment, we validate our approach, explain the molecular mechanism behind identified fragments promoting compound entry, and select candidate compounds from an external library that display good permeation ability.
Asunto(s)
Antibacterianos , Vocabulario , Algoritmos , Antibacterianos/farmacología , Bacterias Gramnegativas , Aprendizaje Automático , Pseudomonas aeruginosaRESUMEN
It is a challenge to obtain an accurate model of the state-to-state dynamics of a complex biological system from molecular dynamics (MD) simulations. In recent years, Markov state models have gained immense popularity for computing state-to-state dynamics from a pool of short MD simulations. However, the assumption that the underlying dynamics on the reduced space is Markovian induces a systematic bias in the model, especially in biomolecular systems with complicated energy landscapes. To address this problem, we have devised a new approach we call quasistationary distribution kinetic Monte Carlo (QSD-KMC) that gives accurate long time state-to-state evolution while retaining the entire time resolution even when the dynamics is highly non-Markovian. The proposed method is a kinetic Monte Carlo approach that takes advantage of two concepts: (i) the quasistationary distribution, the distribution that results when a trajectory remains in one state for a long time (the dephasing time), such that the next escape is Markovian, and (ii) dynamical corrections theory, which properly accounts for the correlated events that occur as a trajectory passes from state to state before it settles again. In practice, this is achieved by specifying, for each escape, the intermediate states and the final state that has resulted from the escape. Implementation of QSD-KMC imposes stricter requirements on the lengths of the trajectories than in a Markov state model approach as the trajectories must be long enough to dephase. However, the QSD-KMC model produces state-to-state trajectories that are statistically indistinguishable from an MD trajectory mapped onto the discrete set of states for an arbitrary choice of state decomposition. Furthermore, the aforementioned concepts can be used to construct a Monte Carlo approach to optimize the state boundaries regardless of the initial choice of states. We demonstrate the QSD-KMC method on two one-dimensional model systems, one of which is a driven nonequilibrium system, and on two well-characterized biomolecular systems.
Asunto(s)
Simulación de Dinámica Molecular , Método de Montecarlo , CinéticaRESUMEN
OBJECTIVES: Thyroid abnormalities in patients with bipolar disorder (BD) have been linked to lithium treatment for decades, yet other drugs have been less well studied. Our objective was to compare hypothyroidism risk for lithium versus the anticonvulsants and second-generation antipsychotics commonly prescribed for BD. METHODS: Administrative claims data on 24,574 patients with BD were analyzed with competing risk survival analysis. Inclusion criteria were (i) one year of no prior hypothyroid diagnosis nor BD drug treatment, (ii) followed by at least one thyroid test during BD monotherapy on lithium carbonate, mood-stabilizing anticonvulsants (lamotrigine, valproate, oxcarbazepine, or carbamazepine) or antipsychotics (aripiprazole, olanzapine, risperidone, or quetiapine). The outcome was cumulative incidence of hypothyroidism per drug, in the presence of the competing risk of ending monotherapy, adjusted for age, sex, physician visits, and thyroid tests. RESULTS: Adjusting for covariates, the four-year cumulative risk of hypothyroidism for lithium (8.8%) was 1.39-fold that of the lowest risk therapy, oxcarbazepine (6.3%). Lithium was non-statistically significantly different from quetiapine. While lithium conferred a higher risk when compared to all other treatments combined as a group, hypothyroidism risk error bars overlapped for all drugs. Treatment (p = 3.86e-3), age (p = 6.91e-10), sex (p = 3.93e-7), and thyroid testing (p = 2.79e-87) affected risk. Patients taking lithium were tested for hypothyroidism 2.26-3.05 times more frequently than those on other treatments. CONCLUSIONS: Thyroid abnormalities occur frequently in patients with BD regardless of treatment. Therefore, patients should be regularly tested for clinical or subclinical thyroid abnormalities on all therapies and treated as indicated to prevent adverse effects of hormone imbalances on mood.
Asunto(s)
Antimaníacos/efectos adversos , Antimaníacos/uso terapéutico , Trastorno Bipolar/tratamiento farmacológico , Hipotiroidismo/inducido químicamente , Adulto , Anticonvulsivantes/efectos adversos , Anticonvulsivantes/uso terapéutico , Antipsicóticos/efectos adversos , Antipsicóticos/uso terapéutico , Estudios de Cohortes , Femenino , Humanos , Compuestos de Litio/efectos adversos , Compuestos de Litio/uso terapéutico , Masculino , Riesgo , Análisis de SupervivenciaRESUMEN
Bacteremia and malaria coinfection is a common and life-threatening condition in children residing in sub-Saharan Africa. We previously showed that coinfection with Gram negative (G[-]) enteric Bacilli and Plasmodium falciparum (Pf[+]) was associated with reduced high-density parasitemia (HDP, >10,000 parasites/µL), enhanced respiratory distress, and severe anemia. Since inflammatory mediators are largely unexplored in such coinfections, circulating cytokines were determined in four groups of children (n = 206, aged <3 yrs): healthy; Pf[+] alone; G[-] coinfected; and G[+] coinfected. Staphylococcus aureus and non-Typhi Salmonella were the most frequently isolated G[+] and G[-] organisms, respectively. Coinfected children, particularly those with G[-] pathogens, had lower parasite burden (peripheral and geometric mean parasitemia and HDP). In addition, both coinfected groups had increased IL-4, IL-5, IL-7, IL-12, IL-15, IL-17, IFN-γ, and IFN-α and decreased TNF-α relative to malaria alone. Children with G[-] coinfection had higher IL-1ß and IL-1Ra and lower IL-10 than the Pf[+] group and higher IFN-γ than the G[+] group. To determine how the immune response to malaria regulates parasitemia, cytokine production was investigated with a multiple mediation model. Cytokines with the greatest mediational impact on parasitemia were IL-4, IL-10, IL-12, and IFN-γ. Results here suggest that enhanced immune activation, especially in G[-] coinfected children, acts to reduce malaria parasite burden.
Asunto(s)
Bacteriemia/microbiología , Bacteriemia/parasitología , Coinfección/sangre , Coinfección/microbiología , Coinfección/parasitología , Malaria Falciparum/microbiología , Malaria Falciparum/parasitología , Bacteriemia/sangre , Preescolar , Femenino , Humanos , Lactante , Recién Nacido , Inflamación/sangre , Inflamación/microbiología , Inflamación/parasitología , Interferón-alfa/sangre , Interferón gamma/sangre , Interleucina-10/sangre , Interleucina-12/sangre , Interleucina-15/sangre , Interleucina-17/sangre , Interleucina-4/sangre , Interleucina-5/sangre , Interleucina-7/sangre , Malaria Falciparum/sangre , Masculino , Salmonella/patogenicidad , Staphylococcus aureus/patogenicidad , Factor de Necrosis Tumoral alfa/sangreRESUMEN
The ability Gram-negative pathogens have at adapting and protecting themselves against antibiotics has increasingly become a public health threat. Data-driven models identifying molecular properties that correlate with outer membrane (OM) permeation and growth inhibition while avoiding efflux could guide the discovery of novel classes of antibiotics. Here we evaluate 174 molecular descriptors in 1260 antimicrobial compounds and study their correlations with antibacterial activity in Gram-negative Pseudomonas aeruginosa. The descriptors are derived from traditional approaches quantifying the compounds' intrinsic physicochemical properties, together with, bacterium-specific from ensemble docking of compounds targeting specific MexB binding pockets, and all-atom molecular dynamics simulations in different subregions of the OM model. Using these descriptors and the measured inhibitory concentrations, we design a statistical protocol to identify predictors of OM permeation/inhibition. We find consistent rules across most of our data highlighting the role of the interaction between the compounds and the OM. An implementation of the rules uncovered in our study is shown, and it demonstrates the accuracy of our approach in a set of previously unseen compounds. Our analysis sheds new light on the key properties drug candidates need to effectively permeate/inhibit P. aeruginosa, and opens the gate to similar data-driven studies in other Gram-negative pathogens.
RESUMEN
We present an ensemble transfer learning method to predict suicide from Veterans Affairs (VA) electronic medical records (EMR). A diverse set of base models was trained to predict a binary outcome constructed from reported suicide, suicide attempt, and overdose diagnoses with varying choices of study design and prediction methodology. Each model used twenty cross-sectional and 190 longitudinal variables observed in eight time intervals covering 7.5 years prior to the time of prediction. Ensembles of seven base models were created and fine-tuned with ten variables expected to change with study design and outcome definition in order to predict suicide and combined outcome in a prospective cohort. The ensemble models achieved c-statistics of 0.73 on 2-year suicide risk and 0.83 on the combined outcome when predicting on a prospective cohort of [Formula: see text] 4.2 M veterans. The ensembles rely on nonlinear base models trained using a matched retrospective nested case-control (Rcc) study cohort and show good calibration across a diversity of subgroups, including risk strata, age, sex, race, and level of healthcare utilization. In addition, a linear Rcc base model provided a rich set of biological predictors, including indicators of suicide, substance use disorder, mental health diagnoses and treatments, hypoxia and vascular damage, and demographics.
Asunto(s)
Carcinoma de Células Renales , Neoplasias Renales , Veteranos , Humanos , Veteranos/psicología , Estudios Retrospectivos , Estudios Transversales , Estudios Prospectivos , Intento de Suicidio , Aprendizaje AutomáticoRESUMEN
Interdependence across time and length scales is common in biology, where atomic interactions can impact larger-scale phenomenon. Such dependence is especially true for a well-known cancer signaling pathway, where the membrane-bound RAS protein binds an effector protein called RAF. To capture the driving forces that bring RAS and RAF (represented as two domains, RBD and CRD) together on the plasma membrane, simulations with the ability to calculate atomic detail while having long time and large length- scales are needed. The Multiscale Machine-Learned Modeling Infrastructure (MuMMI) is able to resolve RAS/RAF protein-membrane interactions that identify specific lipid-protein fingerprints that enhance protein orientations viable for effector binding. MuMMI is a fully automated, ensemble-based multiscale approach connecting three resolution scales: (1) the coarsest scale is a continuum model able to simulate milliseconds of time for a 1 µm2 membrane, (2) the middle scale is a coarse-grained (CG) Martini bead model to explore protein-lipid interactions, and (3) the finest scale is an all-atom (AA) model capturing specific interactions between lipids and proteins. MuMMI dynamically couples adjacent scales in a pairwise manner using machine learning (ML). The dynamic coupling allows for better sampling of the refined scale from the adjacent coarse scale (forward) and on-the-fly feedback to improve the fidelity of the coarser scale from the adjacent refined scale (backward). MuMMI operates efficiently at any scale, from a few compute nodes to the largest supercomputers in the world, and is generalizable to simulate different systems. As computing resources continue to increase and multiscale methods continue to advance, fully automated multiscale simulations (like MuMMI) will be commonly used to address complex science questions.
Asunto(s)
Proteínas de la Membrana , Simulación de Dinámica Molecular , Proteínas de la Membrana/química , Membrana Celular/metabolismo , Aprendizaje Automático , LípidosRESUMEN
Fitting Susceptible-Infected-Recovered (SIR) models to incidence data is problematic when not all infected individuals are reported. Assuming an underlying SIR model with general but known distribution for the time to recovery, this paper derives the implied differential-integral equations for observed incidence data when a fixed fraction of newly infected individuals are not observed. The parameters of the resulting system of differential equations are identifiable. Using these differential equations, we develop a stochastic model for the conditional distribution of current disease incidence given the entire past history of reported cases. We estimate the model parameters using Bayesian Markov Chain Monte-Carlo sampling of the posterior distribution. We use our model to estimate the transmission rate and fraction of asymptomatic individuals for the current Coronavirus 2019 outbreak in eight American Countries: the United States of America, Brazil, Mexico, Argentina, Chile, Colombia, Peru, and Panama, from January 2020 to May 2021. Our analysis reveals that the fraction of reported cases varies across all countries. For example, the reported incidence fraction for the United States of America varies from 0.3 to 0.6, while for Brazil it varies from 0.2 to 0.4.
Asunto(s)
COVID-19/epidemiología , Argentina/epidemiología , Teorema de Bayes , Brasil/epidemiología , Chile/epidemiología , Colombia/epidemiología , Humanos , Incidencia , Cadenas de Markov , México/epidemiología , Panamá/epidemiología , Perú/epidemiología , Procesos Estocásticos , Estados Unidos/epidemiologíaRESUMEN
Discovery of reliable signatures for the empirical diagnosis of neurological diseases-both infectious and non-infectious-remains unrealized. One of the primary challenges encountered in such studies is the lack of a comprehensive database representative of a signature background that exists in healthy individuals, and against which an aberrant event can be assessed. For neurological insults and injuries, it is important to understand the normal profile in the neuronal (cerebrospinal fluid) and systemic fluids (e.g., blood). Here, we present the first comparative multi-omic human database of signatures derived from a population of 30 individuals (15 males, 15 females, 23-74 years) of serum and cerebrospinal fluid. In addition to empirical signatures, we also assigned common pathways between serum and CSF. Together, our findings provide a cohort against which aberrant signature profiles in individuals with neurological injuries/disease can be assessed-providing a pathway for comprehensive diagnostics and therapeutics discovery.
Asunto(s)
Enfermedades del Sistema Nervioso , Proteómica , Líquido Cefalorraquídeo , Estudios de Cohortes , Femenino , Humanos , Masculino , Metabolómica , NeuronasRESUMEN
Severe malarial anemia (SMA) is a leading cause of childhood morbidity and mortality in holoendemic Plasmodium falciparum transmission regions. To gain enhanced understanding of predisposing factors for SMA, we explored the relationship between complement component 3 (C3) missense mutations [rs2230199 (2307C>G, Arg>Gly102) and rs11569534 (34420G>A, Gly>Asp1224)], malaria, and SMA in a cohort of children (n = 1617 children) over 36 months of follow-up. Variants were selected based on their ability to impart amino acid substitutions that can alter the structure and function of C3. The 2307C>G mutation results in a basic to a polar residue change (Arg to Gly) at position 102 (ß-chain) in the macroglobulin-1 (MG1) domain, while 34420G>A elicits a polar to acidic residue change (Gly to Asp) at position 1224 (α-chain) in the thioester-containing domain. After adjusting for multiple comparisons, longitudinal analyses revealed that inheritance of the homozygous mutant (GG) at 2307 enhanced the risk of SMA (RR = 2.142, 95%CI: 1.229-3.735, P = 0.007). The haplotype containing both wild-type alleles (CG) decreased the incident risk ratio of both malaria (RR = 0.897, 95%CI: 0.828-0.972, P = 0.008) and SMA (RR = 0.617, 95%CI: 0.448-0.848, P = 0.003). Malaria incident risk ratio was also reduced in carriers of the GG (Gly102Gly1224) haplotype (RR = 0.941, 95%CI: 0.888-0.997, P = 0.040). Collectively, inheritance of the missense mutations in MG1 and thioester-containing domain influence the longitudinal risk of malaria and SMA in children exposed to intense Plasmodium falciparum transmission.
Asunto(s)
Anemia , Complemento C3 , Malaria Falciparum , Anemia/genética , Anemia/parasitología , Niño , Complemento C3/genética , Predisposición Genética a la Enfermedad , Humanos , Malaria Falciparum/complicaciones , Malaria Falciparum/genética , Mutación , Plasmodium falciparumRESUMEN
The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (â¼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with â¼1.5 million atoms is discussed in the context of MuMMI.
Asunto(s)
Membrana Dobles de Lípidos , Simulación de Dinámica Molecular , Humanos , Membrana Dobles de Lípidos/química , Estructura Secundaria de Proteína , Proteínas/químicaRESUMEN
Seasonal flu is an acute respiratory disease that exacts a massive toll on human populations, healthcare systems and economies. The disease is caused by an enveloped Influenza virus containing eight ribonucleoprotein (RNP) complexes. Each RNP incorporates multiple copies of nucleoprotein (NP), a fragment of the viral genome (vRNA), and a viral RNA-dependent RNA polymerase (POL), and is responsible for packaging the viral genome and performing critical functions including replication and transcription. A complete model of an Influenza RNP in atomic detail can elucidate the structural basis for viral genome functions, and identify potential targets for viral therapeutics. In this work we construct a model of a complete Influenza A RNP complex in atomic detail using multiple sources of structural and sequence information and a series of homology-modeling techniques, including a motif-matching fragment assembly method. Our final model provides a rationale for experimentally-observed changes to viral polymerase activity in numerous mutational assays. Further, our model reveals specific interactions between the three primary structural components of the RNP, including potential targets for blocking POL-binding to the NP-vRNA complex. The methods developed in this work open the possibility of elucidating other functionally-relevant atomic-scale interactions in additional RNP structures and other biomolecular complexes.
Asunto(s)
Virus de la Influenza A/metabolismo , Modelos Biológicos , Nucleoproteínas/metabolismo , ARN Viral/metabolismo , ARN Polimerasas Dirigidas por ADN/metabolismo , Modelos Moleculares , Conformación de Ácido Nucleico , Unión Proteica , Multimerización de Proteína , ARN Viral/química , Relación Estructura-ActividadRESUMEN
Antibiotic-resistant bacteria rapidly spread in clinical and natural environments and challenge our modern lifestyle. A major component of defense against antibiotics in Gram-negative bacteria is a drug permeation barrier created by active efflux across the outer membrane. We identified molecular determinants defining the propensity of small peptidomimetic molecules to avoid and inhibit efflux pumps in Pseudomonas aeruginosa, a human pathogen notorious for its antibiotic resistance. Combining experimental and computational protocols, we mapped the fate of the compounds from structure-activity relationships through their dynamic behavior in solution, permeation across both the inner and outer membranes, and interaction with MexB, the major efflux transporter of P. aeruginosa We identified predictors of efflux avoidance and inhibition and demonstrated their power by using a library of traditional antibiotics and compound series and by generating new inhibitors of MexB. The identified predictors will enable the discovery and optimization of antibacterial agents suitable for treatment of P. aeruginosa infections.IMPORTANCE Efflux pump avoidance and inhibition are desired properties for the optimization of antibacterial activities against Gram-negative bacteria. However, molecular and physicochemical interactions defining the interface between compounds and efflux pumps remain poorly understood. We identified properties that correlate with efflux avoidance and inhibition, are predictive of similar features in structurally diverse compounds, and allow researchers to distinguish between efflux substrates, inhibitors, and avoiders in P. aeruginosa The developed predictive models are based on the descriptors representative of different clusters comprising a physically intuitive combination of properties. Molecular shape (represented by acylindricity), amphiphilicity (anisotropic polarizability), aromaticity (number of aromatic rings), and the partition coefficient (LogD) are physicochemical predictors of efflux inhibitors, whereas interactions with Pro668 and Leu674 residues of MexB distinguish between inhibitors/substrates and efflux avoiders. The predictive models and efflux rules are applicable to compounds with unrelated chemical scaffolds and pave the way for development of compounds with the desired efflux interface properties.
Asunto(s)
Antibacterianos/farmacología , Proteínas de la Membrana Bacteriana Externa/química , Farmacorresistencia Bacteriana Múltiple/efectos de los fármacos , Proteínas de Transporte de Membrana/química , Modelos Biológicos , Peptidomiméticos/farmacología , Pseudomonas aeruginosa/efectos de los fármacos , Antibacterianos/síntesis química , Antibacterianos/metabolismo , Proteínas de la Membrana Bacteriana Externa/antagonistas & inhibidores , Proteínas de la Membrana Bacteriana Externa/genética , Proteínas de la Membrana Bacteriana Externa/metabolismo , Sitios de Unión , Transporte Biológico/efectos de los fármacos , Expresión Génica , Cinética , Proteínas de Transporte de Membrana/genética , Proteínas de Transporte de Membrana/metabolismo , Pruebas de Sensibilidad Microbiana , Modelos Moleculares , Peptidomiméticos/síntesis química , Peptidomiméticos/metabolismo , Análisis de Componente Principal , Unión Proteica , Conformación Proteica en Hélice alfa , Conformación Proteica en Lámina beta , Dominios y Motivos de Interacción de Proteínas , Pseudomonas aeruginosa/genética , Pseudomonas aeruginosa/metabolismo , Relación Estructura-Actividad , TermodinámicaRESUMEN
Dense surface glycosylation on the HIV-1 envelope (Env) protein acts as a shield from the adaptive immune system. However, the molecular complexity and flexibility of glycans make experimental studies a challenge. Here we have integrated high-throughput atomistic modeling of fully glycosylated HIV-1 Env with graph theory to capture immunologically important features of the shield topology. This is the first complete all-atom model of HIV-1 Env SOSIP glycan shield that includes both oligomannose and complex glycans, providing physiologically relevant insights of the glycan shield. This integrated approach including quantitative comparison with cryo-electron microscopy data provides hitherto unexplored details of the native shield architecture and its difference from the high-mannose glycoform. We have also derived a measure to quantify the shielding effect over the antigenic protein surface that defines regions of relative vulnerability and resilience of the shield and can be harnessed for rational immunogen design.
RESUMEN
[This corrects the article DOI: 10.1371/journal.pone.0225858.].
RESUMEN
Around the world, scavenging birds such as vultures and condors have been experiencing drastic population declines. Scavenging birds have a distinct digestive process to deal with higher amounts of bacteria in their primary diet of carcasses in varying levels of decay. These observations motivate us to present an analysis of captive and healthy California condor (Gymnogyps californianus) microbiomes to characterize a population raised together under similar conditions. Shotgun metagenomic DNA sequences were analyzed from fecal and cloacal samples of captive birds. Classification of shotgun DNA sequence data with peptide signatures using the Sequedex package provided both phylogenetic and functional profiles, as well as individually annotated reads for targeted confirmatory analysis. We observed bacterial species previously associated with birds and gut microbiomes, including both virulent and opportunistic pathogens such as Clostridium perfringens, Propionibacterium acnes, Shigella flexneri, and Fusobacterium mortiferum, common flora such as Lactobacillus johnsonii, Lactobacillus ruminus, and Bacteroides vulgatus, and mucosal microbes such as Delftia acidovorans, Stenotrophomonas maltophilia, and Corynebacterium falsnii. Classification using shotgun metagenomic reads from phylogenetic marker genes was consistent with, and more specific than, analysis based on 16S rDNA data. Classification of samples based on either phylogenetic or functional profiles of genomic fragments differentiated three types of samples: fecal, mature cloacal and immature cloacal, with immature birds having approximately 40% higher diversity of microbes.
Asunto(s)
Bacterias , Aves/microbiología , Metagenoma , Microbiota/fisiología , Animales , Bacterias/clasificación , Bacterias/genética , Bacterias/crecimiento & desarrolloRESUMEN
BACKGROUND: Severe malarial anaemia (SMA) is a leading cause of childhood mortality in holoendemic Plasmodium falciparum regions. METHODS: To gain an improved understanding of SMA pathogenesis, whole genome and transcriptome profiling was performed in Kenyan children (n=144, 3-36months) with discrete non-SMA and SMA phenotypes. Leukocyte associated immunoglobulin like receptor 1 (LAIR1) emerged as a predictor of susceptibility to SMA (P<1×10-2, OR: 0.44-1.37), and was suppressed in severe disease (-1.69-fold, P=0.004). To extend these findings, the relationship between LAIR1 polymorphisms [rs6509867 (16231C>A); rs2287827 (18835G>A)] and clinical outcomes were investigated in individuals (n=1512, <5years) at enrolment and during a 36-month longitudinal follow-up. FINDINGS: Inheritance of the 16,231 recessive genotype (AA) increased susceptibility to SMA at enrolment (OR=1.903, 95%CI: 1.252-2.891, P=0.003), and longitudinally (RR=1.527, 95%CI: 1.119-2.083, P=0.008). Carriage of the 18,835 GA genotype protected against SMA cross-sectionally (OR=0.672, 95%CI: 0.480-0.9439, P=0.020). Haplotype carriage (C16231A/G18835A) also altered cross-sectional susceptibility to SMA: CG (OR=0.717, 95%CI: 0.527-0.9675, P=0.034), CA (OR=0.745, 95%CI: 0.536-1.036, P=0.080), and AG (OR=1.641, 95%CI: 1.160-2.321, P=0.005). Longitudinally, CA carriage was protective against SMA (RR=0.715, 95%CI: 0.554-0.923, P=0.010), while AG carriage had an additive effect on enhanced SMA risk (RR=1.283, 95%CI: 1.057-1.557, P=0.011). Variants that protected against SMA had elevated LAIR1 transcripts, while those with enhanced risk had lower expression (P<0.05). Inheritance of 18,835 GA reduced all-cause mortality by 44.8% (HR=0.552, 95%CI: 0.329-0.925, P=0.024), while AG haplotype carriage increased susceptibility by 68% (HR=1.680, 95%CI: 1.020-2.770, P=0.040). INTERPRETATION: These findings suggest LAIR1 is important for modulating susceptibility to SMA and all-cause childhood mortality.