RESUMO
Individual oncogenic KRAS mutants confer distinct differences in biochemical properties and signaling for reasons that are not well understood. KRAS activity is closely coupled to protein dynamics and is regulated through two interconverting conformations: state 1 (inactive, effector binding deficient) and state 2 (active, effector binding enabled). Here, we use 31P NMR to delineate the differences in state 1 and state 2 populations present in WT and common KRAS oncogenic mutants (G12C, G12D, G12V, G13D, and Q61L) bound to its natural substrate GTP or a commonly used nonhydrolyzable analog GppNHp (guanosine-5'-[(ß,γ)-imido] triphosphate). Our results show that GppNHp-bound proteins exhibit significant state 1 population, whereas GTP-bound KRAS is primarily (90% or more) in state 2 conformation. This observation suggests that the predominance of state 1 shown here and in other studies is related to GppNHp and is most likely nonexistent in cells. We characterize the impact of this differential conformational equilibrium of oncogenic KRAS on RAF1 kinase effector RAS-binding domain and intrinsic hydrolysis. Through a KRAS G12C drug discovery, we have identified a novel small-molecule inhibitor, BBO-8956, which is effective against both GDP- and GTP-bound KRAS G12C. We show that binding of this inhibitor significantly perturbs state 1-state 2 equilibrium and induces an inactive state 1 conformation in GTP-bound KRAS G12C. In the presence of BBO-8956, RAF1-RAS-binding domain is unable to induce a signaling competent state 2 conformation within the ternary complex, demonstrating the mechanism of action for this novel and active-conformation inhibitor.
Assuntos
Proteínas Proto-Oncogênicas p21(ras) , Proteínas ras , Proteínas Proto-Oncogênicas p21(ras)/genética , Proteínas Proto-Oncogênicas p21(ras)/metabolismo , Proteínas ras/metabolismo , Guanosina Trifosfato/metabolismo , Espectroscopia de Ressonância Magnética , Transdução de Sinais , MutaçãoRESUMO
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.
Assuntos
Membrana Celular/enzimologia , Lipídeos/química , Aprendizado de Máquina , Simulação de Dinâmica Molecular , Multimerização Proteica , Proteínas Proto-Oncogênicas p21(ras)/química , Transdução de Sinais , HumanosRESUMO
Protein-ligand interactions are essential to drug discovery and drug development efforts. Desirable on-target or multitarget interactions are the first step in finding an effective therapeutic, while undesirable off-target interactions are the first step in assessing safety. In this work, we introduce a novel ligand-based featurization and mapping of human protein pockets to identify closely related protein targets and to project novel drugs into a hybrid protein-ligand feature space to identify their likely protein interactions. Using structure-based template matches from PDB, protein pockets are featured by the ligands that bind to their best co-complex template matches. The simplicity and interpretability of this approach provide a granular characterization of the human proteome at the protein-pocket level instead of the traditional protein-level characterization by family, function, or pathway. We demonstrate the power of this featurization method by clustering a subset of the human proteome and evaluating the predicted cluster associations of over 7000 compounds.
Assuntos
Proteoma , Humanos , Ligação Proteica , Sítios de Ligação , Conformação Proteica , Ligantes , Análise por ConglomeradosRESUMO
An Online tool for Fragment-based Molecule Parametrization (OFraMP) is described. OFraMP is a web application for assigning atomic interaction parameters to large molecules by matching sub-fragments within the target molecule to equivalent sub-fragments within the Automated Topology Builder (ATB, atb.uq.edu.au) database. OFraMP identifies and compares alternative molecular fragments from the ATB database, which contains over 890,000 pre-parameterized molecules, using a novel hierarchical matching procedure. Atoms are considered within the context of an extended local environment (buffer region) with the degree of similarity between an atom in the target molecule and that in the proposed match controlled by varying the size of the buffer region. Adjacent matching atoms are combined into progressively larger matched sub-structures. The user then selects the most appropriate match. OFraMP also allows users to manually alter interaction parameters and automates the submission of missing substructures to the ATB in order to generate parameters for atoms in environments not represented in the existing database. The utility of OFraMP is illustrated using the anti-cancer agent paclitaxel and a dendrimer used in organic semiconductor devices. OFraMP applied to paclitaxel (ATB ID 35922).
Assuntos
Software , Bases de Dados FactuaisRESUMO
During the activation of mitogen-activated protein kinase (MAPK) signaling, the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF bind to active RAS at the plasma membrane. The orientation of RAS at the membrane may be critical for formation of the RAS-RBDCRD complex and subsequent signaling. To explore how RAS membrane orientation relates to the protein dynamics within the RAS-RBDCRD complex, we perform multiscale coarse-grained and all-atom molecular dynamics (MD) simulations of KRAS4b bound to the RBD and CRD domains of RAF-1, both in solution and anchored to a model plasma membrane. Solution MD simulations describe dynamic KRAS4b-CRD conformations, suggesting that the CRD has sufficient flexibility in this environment to substantially change its binding interface with KRAS4b. In contrast, when the ternary complex is anchored to the membrane, the mobility of the CRD relative to KRAS4b is restricted, resulting in fewer distinct KRAS4b-CRD conformations. These simulations implicate membrane orientations of the ternary complex that are consistent with NMR measurements. While a crystal structure-like conformation is observed in both solution and membrane simulations, a particular intermolecular rearrangement of the ternary complex is observed only when it is anchored to the membrane. This configuration emerges when the CRD hydrophobic loops are inserted into the membrane and helices α3-5 of KRAS4b are solvent exposed. This membrane-specific configuration is stabilized by KRAS4b-CRD contacts that are not observed in the crystal structure. These results suggest modulatory interplay between the CRD and plasma membrane that correlate with RAS/RAF complex structure and dynamics, and potentially influence subsequent steps in the activation of MAPK signaling.
Assuntos
Cisteína , Proteínas Proto-Oncogênicas c-raf , Sítios de Ligação , Membrana Celular/metabolismo , Cisteína/metabolismo , Proteínas Quinases Ativadas por Mitógeno/metabolismo , Ligação Proteica , Proteínas Proto-Oncogênicas c-raf/química , Proteínas Proto-Oncogênicas c-raf/metabolismo , Proteínas Proto-Oncogênicas p21(ras)/metabolismo , Solventes/metabolismoRESUMO
BACKGROUND: The 2020 consensus guidelines for vancomycin therapeutic monitoring recommend using Bayesian estimation targeting the ratio of the area under the curve over 24 hours to minimum inhibitory concentration as an optimal approach to individualize therapy in pediatric patients. To support institutional guideline implementation in children, the objective of this study was to comprehensively assess and compare published population-based pharmacokinetic (PK) vancomycin models and available Bayesian estimation tools, specific to neonatal and pediatric patients. METHODS: PubMed and Embase databases were searched from January 1994 to December 2020 for studies in which a vancomycin population PK model was developed to determine clearance and volume of distribution in neonatal and pediatric populations. Available Bayesian software programs were identified and assessed from published articles, software program websites, and direct communication with the software company. In the present review, 14 neonatal and 20 pediatric models were included. Six programs (Adult and Pediatric Kinetics, BestDose, DoseMeRx, InsightRx, MwPharm++, and PrecisePK) were evaluated. RESULTS: Among neonatal models, Frymoyer et al and Capparelli et al used the largest PK samples to generate their models, which were externally validated. Among the pediatric models, Le et al used the largest sample size, with multiple external validations. Of the Bayesian programs, DoseMeRx, InsightRx, and PrecisePK used clinically validated neonatal and pediatric models. CONCLUSIONS: To optimize vancomycin use in neonatal and pediatric patients, clinicians should focus on selecting a model that best fits their patient population and use Bayesian estimation tools for therapeutic area under the -curve-targeted dosing and monitoring.
Assuntos
Software , Vancomicina , Adulto , Antibacterianos/farmacocinética , Área Sob a Curva , Teorema de Bayes , Criança , Humanos , Recém-Nascido , Cinética , Testes de Sensibilidade Microbiana , Vancomicina/farmacocinéticaRESUMO
Predicting accurate protein-ligand binding affinities is an important task in drug discovery but remains a challenge even with computationally expensive biophysics-based energy scoring methods and state-of-the-art deep learning approaches. Despite the recent advances in the application of deep convolutional and graph neural network-based approaches, it remains unclear what the relative advantages of each approach are and how they compare with physics-based methodologies that have found more mainstream success in virtual screening pipelines. We present fusion models that combine features and inference from complementary representations to improve binding affinity prediction. This, to our knowledge, is the first comprehensive study that uses a common series of evaluations to directly compare the performance of three-dimensional (3D)-convolutional neural networks (3D-CNNs), spatial graph neural networks (SG-CNNs), and their fusion. We use temporal and structure-based splits to assess performance on novel protein targets. To test the practical applicability of our models, we examine their performance in cases that assume that the crystal structure is not available. In these cases, binding free energies are predicted using docking pose coordinates as the inputs to each model. In addition, we compare these deep learning approaches to predictions based on docking scores and molecular mechanic/generalized Born surface area (MM/GBSA) calculations. Our results show that the fusion models make more accurate predictions than their constituent neural network models as well as docking scoring and MM/GBSA rescoring, with the benefit of greater computational efficiency than the MM/GBSA method. Finally, we provide the code to reproduce our results and the parameter files of the trained models used in this work. The software is available as open source at https://github.com/llnl/fast. Model parameter files are available at ftp://gdo-bioinformatics.ucllnl.org/fast/pdbbind2016_model_checkpoints/.
Assuntos
Redes Neurais de Computação , Proteínas , Ligantes , Ligação Proteica , Proteínas/metabolismo , SoftwareRESUMO
We present a new approach to estimate the binding affinity from given three-dimensional poses of protein-ligand complexes. In this scheme, every protein-ligand atom pair makes an additive free-energy contribution. The sum of these pairwise contributions then gives the total binding free energy or the logarithm of the dissociation constant. The pairwise contribution is calculated by a function implemented via a neural network that takes the properties of the two atoms and their distance as input. The pairwise function is trained using a portion of the PDBbind 2018 data set. The model achieves good accuracy for affinity predictions when evaluated with PDBbind 2018 and with the CASF-2016 benchmark, comparing favorably to many scoring functions such as that of AutoDock Vina. The framework here may be extended to incorporate other factors to further improve its accuracy and power.
Assuntos
Desenho de Fármacos , Redes Neurais de Computação , Ligantes , Simulação de Acoplamento Molecular , Ligação ProteicaRESUMO
Accurately predicting small molecule partitioning and hydrophobicity is critical in the drug discovery process. There are many heterogeneous chemical environments within a cell and entire human body. For example, drugs must be able to cross the hydrophobic cellular membrane to reach their intracellular targets, and hydrophobicity is an important driving force for drug-protein binding. Atomistic molecular dynamics (MD) simulations are routinely used to calculate free energies of small molecules binding to proteins, crossing lipid membranes, and solvation but are computationally expensive. Machine learning (ML) and empirical methods are also used throughout drug discovery but rely on experimental data, limiting the domain of applicability. We present atomistic MD simulations calculating 15,000 small molecule free energies of transfer from water to cyclohexane. This large data set is used to train ML models that predict the free energies of transfer. We show that a spatial graph neural network model achieves the highest accuracy, followed closely by a 3D-convolutional neural network, and shallow learning based on the chemical fingerprint is significantly less accurate. A mean absolute error of â¼4 kJ/mol compared to the MD calculations was achieved for our best ML model. We also show that including data from the MD simulation improves the predictions, tests the transferability of each model to a diverse set of molecules, and show multitask learning improves the predictions. This work provides insight into the hydrophobicity of small molecules and ML cheminformatics modeling, and our data set will be useful for designing and testing future ML cheminformatics methods.
Assuntos
Aprendizado Profundo , Simulação de Dinâmica Molecular , Entropia , Humanos , Interações Hidrofóbicas e Hidrofílicas , TermodinâmicaRESUMO
We have implemented the Martini force field within Lawrence Livermore National Laboratory's molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein-lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns. The ddcMD code has been made open source and is available on GitHub at https://github.com/LLNL/ddcMD.
RESUMO
Gain-of-function mutations in PCSK9 (proprotein convertase subtilisin/kexin type 9) lead to reduced uptake of LDL (low density lipoprotein) cholesterol and, therefore, increased plasma LDL levels. However, the mechanism by which these mutants reduce LDL reuptake is not fully understood. Here, we have used molecular dynamics simulations, MM/PBSA (Molecular Mechanics/Poisson-Boltzmann Surface Area) binding affinity calculations, and residue interaction networks, to investigate the protein-protein interaction (PPI) disruptive effects of two of PCSK9's gain-of-function mutations, Ser127Arg and Asp374Tyr on the PCSK9 and LDL receptor complex. In addition to these PPI disruptive mutants, a third, non-interface mutation (Arg496Trp) is included as a positive control. Our results indicate that Ser127Arg and Asp374Tyr confer significantly improved binding affinity, as well as different binding modes, when compared to the wild-type. These PPI disruptive mutations lie between the EGF(A) (epidermal growth factor precursor homology domain A) of the LDL receptor and the catalytic domain of PCSK9 (Asp374Tyr) and between the prodomain of PCSK9 and the ß-propeller of the LDL receptor (Ser127Arg). The interactions involved in these two interfaces result in an LDL receptor that is sterically inhibited from entering its closed conformation. This could potentially implicate the prodomain as a target for small molecule inhibitors.
Assuntos
Pró-Proteína Convertase 9/química , Receptores de LDL/química , Domínio Catalítico , Humanos , Simulação de Dinâmica Molecular , Mutação , Pró-Proteína Convertase 9/genética , Pró-Proteína Convertase 9/metabolismo , Ligação Proteica , Receptores de LDL/genética , Receptores de LDL/metabolismoRESUMO
Membrane protein functions can be altered by subtle changes in the host lipid bilayer physical properties. Gramicidin channels have emerged as a powerful system for elucidating the underlying mechanisms of membrane protein function regulation through changes in bilayer properties, which are reflected in the thermodynamic equilibrium distribution between nonconducting gramicidin monomers and conducting bilayer-spanning dimers. To improve our understanding of how subtle changes in bilayer thickness alter the gramicidin monomer and dimer distributions, we performed extensive atomistic molecular dynamics simulations and fluorescence-quenching experiments on gramicidin A (gA). The free-energy calculations predicted a nonlinear coupling between the bilayer thickness and channel formation. The energetic barrier inhibiting gA channel formation was sharply increased in the thickest bilayer (1,2-dierucoyl-sn-glycero-3-phosphocholine). This prediction was corroborated by experimental results on gramicidin channel activity in bilayers of different thickness. To further explore the mechanism of channel formation, we performed extensive unbiased molecular dynamics simulations, which allowed us to observe spontaneous gA dimer formation in lipid bilayers. The simulations revealed structural rearrangements in the gA subunits and changes in lipid packing, as well as water reorganization, that occur during the dimerization process. Together, the simulations and experiments provide new, to our knowledge, insights into the process and mechanism of gramicidin channel formation, as a prototypical example of the bilayer regulation of membrane protein function.
Assuntos
Dimerização , Gramicidina/química , Bicamadas Lipídicas/química , Fluorescência , Cinética , Simulação de Dinâmica Molecular , Termodinâmica , Água/químicaRESUMO
Collaboration between gastroenterologists and rheumatologists is recommended for the correct management of patients with associated spondyloarthritis (SpA) and inflammatory bowel disease (IBD). We aimed to establish the appropriateness of several red flags for a prompt specialist referral. A systematic review of the literature was performed using the GRADE method to describe the prevalence of co-existing IBD-SpA and the diagnostic accuracy of red flags proposed by a steering committee. Then, a consensus among expert gastroenterologists and rheumatologists (10 in the steering committee and 13 in the expert panel) was obtained using the RAND method to confirm the appropriateness of each red flag as 'major' (one sufficient for patient referral) or 'minor' (at least three needed for patient referral) criteria for specialist referral. The review of the literature confirmed the high prevalence of co-existing IBD-SpA. Positive and negative predictive values of red flags were not calculated, given the lack of available data. A consensus among gastroenterology and rheumatology specialists was used to confirm the appropriateness of each red flag. Major criteria to refer patients with SpA to the gastroenterologist included: rectal bleeding, chronic abdominal pain, perianal fistula or abscess, chronic diarrhoea and nocturnal symptoms. Major criteria to refer patients with IBD to the rheumatologist included: chronic low back pain, dactylitis, enthesitis and pain/swelling of peripheral joints. Several major and minor red flags have been identified for the diagnosis of co-existing IBD-SpA. The use of red flags in routine clinical practice may avoid diagnostic delay and reduce clinic overload.
Assuntos
Gastroenterologistas , Doenças Inflamatórias Intestinais/diagnóstico , Reto/patologia , Reumatologistas , Espondilite Anquilosante/diagnóstico , Dor Abdominal , Consenso , Diarreia , Doença , Prova Pericial , Hemorragia , Humanos , Doenças Inflamatórias Intestinais/epidemiologia , Guias de Prática Clínica como Assunto , Prevalência , Encaminhamento e Consulta , Espondilite Anquilosante/epidemiologiaRESUMO
Protein engineering to alter recognition underlying ligand binding and activity has enormous potential. Here, ligand binding for Escherichia coli phosphoenolpyruvate carboxykinase (PEPCK), which converts oxaloacetate into CO2 and phosphoenolpyruvate as the first committed step in gluconeogenesis, was engineered to accommodate alternative ligands as an exemplary system with structural information. From our identification of bicarbonate binding in the PEPCK active site at the supposed CO2 binding site, we probed binding of nonnative ligands with three oxygen atoms arranged to resemble the bicarbonate geometry. Crystal structures of PEPCK and point mutants with bound nonnative ligands thiosulfate and methanesulfonate along with strained ATP and reoriented oxaloacetate intermediates and unexpected bicarbonate were determined and analyzed. The mutations successfully altered the bound ligand position and orientation and its specificity: mutated PEPCKs bound either thiosulfate or methanesulfonate but never both. Computational calculations predicted a methanesulfonate binding mutant and revealed that release of the active site ordered solvent exerts a strong influence on ligand binding. Besides nonnative ligand binding, one mutant altered the Mn2+ coordination sphere: instead of the canonical octahedral ligand arrangement, the mutant in question had an only five-coordinate arrangement. From this work, critical features of ligand binding, position, and metal ion cofactor geometry required for all downstream events can be engineered with small numbers of mutations to provide insights into fundamental underpinnings of protein-ligand recognition. Through structural and computational knowledge, the combination of designed and random mutations aids in the robust design of predetermined changes to ligand binding and activity to engineer protein function.
Assuntos
Proteínas de Escherichia coli/química , Proteínas de Escherichia coli/metabolismo , Fosfoenolpiruvato Carboxiquinase (ATP)/química , Fosfoenolpiruvato Carboxiquinase (ATP)/metabolismo , Substituição de Aminoácidos , Domínio Catalítico/genética , Cristalografia por Raios X , Escherichia coli/enzimologia , Escherichia coli/genética , Proteínas de Escherichia coli/genética , Ligação de Hidrogênio , Cinética , Ligantes , Modelos Moleculares , Mutagênese Sítio-Dirigida , Fosfoenolpiruvato Carboxiquinase (ATP)/genética , Conformação Proteica , Engenharia de Proteínas , Eletricidade Estática , Especificidade por SubstratoRESUMO
Solithromycin is a novel fluoroketolide antibiotic which was under investigation for the treatment of community-acquired bacterial pneumonia (CABP). A phase 1 study was performed to characterize the pharmacokinetics (PK) and safety of solithromycin in children. Eighty-four subjects (median age, 6 years [age range, 4 days to 17 years]) were administered intravenous (i.v.) or oral (capsules or suspension) solithromycin (i.v., 6 to 8 mg/kg of body weight; capsules/suspension, 14 to 16 mg/kg on days 1 and 7 to 15 mg/kg on days 2 to 5). PK samples were collected after the first and multidose administration. Data from 83 subjects (662 samples) were combined with previously collected adolescent PK data (n = 13; median age, 16 years [age range, 12 to 17 years]) following capsule administration to perform a population PK analysis. A 2-compartment PK model characterized the data well, and postmenstrual age was the only significant covariate after accounting for body size differences. Dosing simulations suggested that 8 mg/kg i.v. daily and oral dosing of 20 mg/kg on day 1 (800-mg adult maximum) followed by 10 mg/kg on days 2 to 5 (400-mg adult maximum) would achieve a pediatric solithromycin exposure consistent with the exposures observed in adults. Seventy-six treatment-emergent adverse events (TEAEs) were reported in 40 subjects. Diarrhea (6 subjects) and infusion site pain or phlebitis (3 subjects) were the most frequently reported adverse events related to treatment. Two subjects experienced TEAEs of increased hepatic enzymes that were deemed not to be related to the study treatment. (The phase 1 pediatric studies discussed in this paper have been registered at ClinicalTrials.gov under identifiers NCT01966055 and NCT02268279.).
Assuntos
Macrolídeos/efeitos adversos , Macrolídeos/farmacocinética , Triazóis/efeitos adversos , Triazóis/farmacocinética , Administração Intravenosa , Administração Oral , Adolescente , Adulto , Idoso , Idoso de 80 Anos ou mais , Criança , Pré-Escolar , Feminino , Humanos , Lactente , Recém-Nascido , Macrolídeos/administração & dosagem , Masculino , Pessoa de Meia-Idade , Triazóis/administração & dosagem , Adulto JovemRESUMO
Hazard identification is the key step in risk assessment and management of manufactured nanomaterials (NM). However, the rapid commercialisation of nano-enabled products continues to out-pace the development of a prudent risk management mechanism that is widely accepted by the scientific community and enforced by regulators. However, a growing body of academic literature is developing promising quantitative methods. Two approaches have gained significant currency. Bayesian networks (BN) are a probabilistic, machine learning approach while the weight of evidence (WoE) statistical framework is based on expert elicitation. This comparative study investigates the efficacy of quantitative WoE and Bayesian methodologies in ranking the potential hazard of metal and metal-oxide NMs-TiO2, Ag, and ZnO. This research finds that hazard ranking is consistent for both risk assessment approaches. The BN and WoE models both utilize physico-chemical, toxicological, and study type data to infer the hazard potential. The BN exhibits more stability when the models are perturbed with new data. The BN has the significant advantage of self-learning with new data; however, this assumes all input data is equally valid. This research finds that a combination of WoE that would rank input data along with the BN is the optimal hazard assessment framework.
Assuntos
Substâncias Perigosas/análise , Substâncias Perigosas/química , Nanoestruturas/química , Medição de Risco/métodos , Algoritmos , Teorema de Bayes , Fenômenos Químicos , Modelos Teóricos , Método de Monte Carlo , Reprodutibilidade dos Testes , Gestão de Riscos/métodosRESUMO
Membrane lipid composition varies greatly within submembrane compartments, different organelle membranes, and also between cells of different cell stage, cell and tissue types, and organisms. Environmental factors (such as diet) also influence membrane composition. The membrane lipid composition is tightly regulated by the cell, maintaining a homeostasis that, if disrupted, can impair cell function and lead to disease. This is especially pronounced in the brain, where defects in lipid regulation are linked to various neurological diseases. The tightly regulated diversity raises questions on how complex changes in composition affect overall bilayer properties, dynamics, and lipid organization of cellular membranes. Here, we utilize recent advances in computational power and molecular dynamics force fields to develop and test a realistically complex human brain plasma membrane (PM) lipid model and extend previous work on an idealized, "average" mammalian PM. The PMs showed both striking similarities, despite significantly different lipid composition, and interesting differences. The main differences in composition (higher cholesterol concentration and increased tail unsaturation in brain PM) appear to have opposite, yet complementary, influences on many bilayer properties. Both mixtures exhibit a range of dynamic lipid lateral inhomogeneities ("domains"). The domains can be small and transient or larger and more persistent and can correlate between the leaflets depending on lipid mixture, Brain or Average, as well as on the extent of bilayer undulations.
Assuntos
Membrana Celular/metabolismo , Lipídeos de Membrana/química , Lipídeos de Membrana/metabolismo , Neurônios/citologia , Humanos , Modelos Moleculares , Conformação MolecularRESUMO
The γ-aminobutyric acid type A receptor (GABAA-R) is a major inhibitory neuroreceptor that is activated by the binding of GABA. The structure of the GABAA-R is well characterized, and many of the binding site residues have been identified. However, most of these residues are obscured behind the C-loop that acts as a cover to the binding site. Thus, the mechanism by which the GABA molecule recognizes the binding site, and the pathway it takes to enter the binding site are both unclear. Through the completion and detailed analysis of 100 short, unbiased, independent molecular dynamics simulations, we have investigated this phenomenon of GABA entering the binding site. In each system, GABA was placed quasi-randomly near the binding site of a GABAA-R homology model, and atomistic simulations were carried out to observe the behavior of the GABA molecules. GABA fully entered the binding site in 19 of the 100 simulations. The pathway taken by these molecules was consistent and non-random; the GABA molecules approach the binding site from below, before passing up behind the C-loop and into the binding site. This binding pathway is driven by long-range electrostatic interactions, whereby the electrostatic field acts as a 'funnel' that sweeps the GABA molecules towards the binding site, at which point more specific atomic interactions take over. These findings define a nuanced mechanism whereby the GABAA-R uses the general zwitterionic features of the GABA molecule to identify a potential ligand some 2 nm away from the binding site.
Assuntos
Receptores de GABA-A/química , Receptores de GABA-A/metabolismo , Ácido gama-Aminobutírico/metabolismo , Animais , Sítios de Ligação , Biologia Computacional , Simulação por Computador , Humanos , Ativação do Canal Iônico , Ligantes , Modelos Moleculares , Simulação de Dinâmica Molecular , Domínios e Motivos de Interação entre Proteínas , Subunidades Proteicas , Eletricidade EstáticaRESUMO
Use of the highly toxic and easily prepared rodenticide tetramethylenedisulfotetramine (TETS) was banned after thousands of accidental or intentional human poisonings, but it is of continued concern as a chemical threat agent. TETS is a noncompetitive blocker of the GABA type A receptor (GABAAR), but its molecular interaction has not been directly established for lack of a suitable radioligand to localize the binding site. We synthesized [(14)C]TETS (14 mCi/mmol, radiochemical purity >99%) by reacting sulfamide with H(14)CHO and s-trioxane then completion of the sequential cyclization with excess HCHO. The outstanding radiocarbon sensitivity of accelerator mass spectrometry (AMS) allowed the use of [(14)C]TETS in neuroreceptor binding studies with rat brain membranes in comparison with the standard GABAAR radioligand 4'-ethynyl-4-n-[(3)H]propylbicycloorthobenzoate ([(3)H]EBOB) (46 Ci/mmol), illustrating the use of AMS for characterizing the binding sites of high-affinity (14)C radioligands. Fourteen noncompetitive antagonists of widely diverse chemotypes assayed at 1 or 10 µM inhibited [(14)C]TETS and [(3)H]EBOB binding to a similar extent (r(2) = 0.71). Molecular dynamics simulations of these 14 toxicants in the pore region of the α1ß2γ2 GABAAR predict unique and significant polar interactions for TETS with α1T1' and γ2S2', which are not observed for EBOB or the GABAergic insecticides. Several GABAAR modulators similarly inhibited [(14)C]TETS and [(3)H]EBOB binding, including midazolam, flurazepam, avermectin Ba1, baclofen, isoguvacine, and propofol, at 1 or 10 µM, providing an in vitro system for recognizing candidate antidotes.
Assuntos
Hidrocarbonetos Aromáticos com Pontes/metabolismo , Antagonistas de Receptores de GABA-A/metabolismo , Receptores de GABA-A/metabolismo , Amidas/química , Animais , Ligação Competitiva/efeitos dos fármacos , Compostos Bicíclicos Heterocíclicos com Pontes/química , Compostos Bicíclicos Heterocíclicos com Pontes/metabolismo , Hidrocarbonetos Aromáticos com Pontes/síntese química , Hidrocarbonetos Aromáticos com Pontes/química , Isótopos de Carbono , Radioisótopos de Carbono , Formaldeído/química , Agonistas GABAérgicos/farmacologia , Antagonistas de Receptores de GABA-A/química , Compostos Heterocíclicos/química , Humanos , Hipnóticos e Sedativos/farmacologia , Inseticidas/química , Inseticidas/metabolismo , Ácidos Isonicotínicos/farmacologia , Modelos Moleculares , Conformação Molecular , Estrutura Molecular , Propofol/farmacologia , Piridoxina/farmacologia , Ensaio Radioligante , Ratos , Enxofre/química , Complexo Vitamínico B/farmacologiaRESUMO
OBJECTIVES: To identify frequent MRI features of parathyroid adenomas (PTAs) in patients with primary hyperparathyroidism (PHPT) using a fast protocol with a 3 T magnet. METHODS: Thirty-eight patients with PHPT underwent a 3 T-MR. All patients had positive US and Tc-99 sestamibi, for a total number of 46 PTAs. T2-weighted IDEAL-FSE and T1 IDEAL-sequences, before and after contrast, were performed. Five features of PTAs were recognised: hyperintensity, homogeneous or "marbled" appearance and elongated morphology on T2-sequences; cleavage plane from thyroid gland on T2-outphase; rapid enhancement in post-contrast T1. Image quality for T2-weighted IDEAL FSE and usefulness for IDEAL post-contrast T1-weighted and T2-outphase sequences were also graded. RESULTS: PTAs were hyperintense in T2-sequences in 44/46 (95.7%), "marbled" in 30/46 (65.2%) and elongated in 38/46 (82.6%) patients. Cleavage plane was observed in 36/46 (78.3%), and rapid enhancement in 20/46 (43.5%) patients. T2-sequences showed both excellent fat suppression and image quality (average scores of 3.2 and 3.1). T2-outphase images demonstrated to be quite useful (score 2.8), whereas, post-contrast T1 images showed a lower degree of utility (score 2.4). CONCLUSIONS: A fast protocol with 3.0-T MRI, recognising most common features of PTAs, may be used as a second-line method in the preoperative detection of PTAs. KEY POINTS: 3 T MRI protocol based on T2-weighted IDEAL FSE sequences was used. T2-hyperintensity and elongated morphology are common features of PTAs. 3 T MRI could be used in the preoperative detection of PTAs.