The tilting motion of the central core reveals the transport mechanism of the sarco/endoplasmic reticulum Ca2+-ATPase.

The sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) transports two Ca2+ ions per ATP hydrolyzed from the cytoplasm to the lumen. However, how the ATP hydrolysis remotely drives the Ca2+ transport is unclear. In the SERCA1a crystal structures, the ATP hydrolysis is accompanied by the notably increasing tilting angle of the central core (CC) and the Ca2+ transport, and the CC tilting angle dramatically decreases in the E2 to E1 transition. We demonstrated that the significantly increasing tilting motion of the CC drove the Ca2+ release in the molecular dynamics simulation of the R836A variant, and the dramatic spontaneous decrease in the CC tilting angle of the E2 state triggers the restart of the SERCA1a's transport cycle. The repulsion between the phosphorylated D351 and the phosphate groups in ADP triggers the release of ADP from the SERCA1a headpiece. We proposed a novel SERCA transport mechanism in which ATP hydrolysis drives a significant tilting motion of the CC, which drives Ca2+ transport and the A domain rotational motion in the E1P-ADP-2Ca2+ to E2P transition. The dramatic spontaneous decrease in the CC tilting angle of the E2 state drives the restart of the transport cycle.

Pharmacophore Oriented MP2 Characterization of Charge Distribution for Anti-SARS-CoV-2 Inhibitor Nirmatrelvir.

Quantum mechanical second order Møller-Plesset (MP2) perturbation theory and density functional theory (DFT) Becke, 3-parameter, Lee-Yang-Parr (B3LYP) and Minnesota 2006 local functional (M06L) calculations were performed to optimize structure of nirmatrelvir and compute the Merz-Kollman electrostatic potential (MK ESP), natural population analysis (NPA), Hirshfeld, charge model 5 (CM5), and mulliken partial charges. The mulliken partial charge distribution of nirmatrelvir exhibits a poor correlation with the MK ESP charges in MP2, B3LYP, and M06L calculations respectively. The NPA, Hirshfeld, and CM5 partial charge scheme of nirmatrelvir indicate a reasonable correlation with MK ESP charge assignments in B3LYP and M06L calculations. The above correlations were not improved by the inclusion of implicit solvation model. The MK ESP and CM5 partial charges show a strong correlation between the results of MP2 and two DFT methods. The three optimized structures present a certain degree of differences from the crystal bioactive conformation of nirmatrelvir, suggesting the nirmatrelvir-enzyme complex is formed in the induced-fit model. The Reactivity of warhead electrophilic nitrile is justified by the relatively weaker strength of π bonds in the MP2 calculations. The nirmatrelvir hydrogen bond acceptors consistently show strong delocalization of lone pair electrons in three calculations, whereas hydrogen bond donors are found to have high polarization on the heavy nitrogen atoms in MP2 computations. This work helps to parametrize the force field of nirmatrelvir and improve accuracy of molecular docking and rational inhibitor design.

QM/MM study of N501 involved intermolecular interaction between SARS-CoV-2 receptor binding domain and antibody of human origin.

Intermolecular interaction between key residue N501 of the epitope on SARS-CoV-2 RBD and screening antibody B38 was studied using the QM/MM and QM approach. The QM/MM optimized geometry shows that angle X-H---Y is 165° for O-H---O between mAb light chain S30 and RBD N501. High level MP2 calculations indicated the interaction between RBD N501 and S30 of B38 Fab light chain provide a relatively strong attractive force of - 3.32 kcal/mol, whereas the hydrogen bond between RBD Q498 and S30 was quantified as 0.10 kcal/mol. The decrease in ESP partial charge on hydrogen atom of hydroxyl group on S30 drops from 0.38 a.u. to 0.31 a.u., exhibiting the sharing of 0.07 a.u. from the lone pair electron oxygen of N501 due to hydrogen bond formation. The NBO occupancy of hydrogen atom also decreases from 25.79 % to 22.93 % in the hydroxyl H-O NBO bond of S30. However, the minor change of NBO hybridization of hydroxyl oxygen of S30 from sp3.00 to sp3.05 implies the rigidity of hydrogen bond tetrahedral geometry in the relative dynamic protein complex. The O-H---O angle is 165° which is close but not exactly linear. The structural requirement for sp3 hybridization of oxygen for hydroxyl group on S30 and dimension of protein likely prevent O-H---O from adopting linear geometry. The hydrogen bond strengths were also calculated using a variety of DFT methods, and the result of - 3.33 kcal/mol from the M06L method is the closest to that of the MP2 calculation. Results of this work may aid in the COVID-19 vaccine and drug screening.

Heterotropic roles of divalent cations in the establishment of allostery and affinity maturation of integrin αXß2.

Allosteric activation and silencing of leukocyte ß2-integrins transpire through cation-dependent structural changes, which mediate integrin biosynthesis and recycling, and are essential to designing leukocyte-specific drugs. Stepwise addition of Mg2+ reveals two mutually coupled events for the αXß2 ligand-binding domain-the αX I-domain-corresponding to allostery establishment and affinity maturation. Electrostatic alterations in the Mg2+-binding site establish long-range couplings, leading to both pH- and Mg2+-occupancy-dependent biphasic stability change in the αX I-domain fold. The ligand-binding sensorgrams show composite affinity events for the αX I-domain accounting for the multiplicity of the αX I-domain conformational states existing in the solution. On cell surfaces, increasing Mg2+ concentration enhanced adhesiveness of αXß2. This work highlights how intrinsically flexible pH- and cation-sensitive architecture endows a unique dynamic continuum to the αI-domain structure on the intact integrin, thereby revealing the importance of allostery establishment and affinity maturation in both extracellular and intracellular integrin events.

Antibody mix-and-read assays based on fluorescence intensity probes.

Antibodies and Fc fusion proteins are a rapidly growing class of pharmaceuticals. Cell culture and purification process development and operation require frequent measurement of product concentrations, commonly by complex enzyme-linked immunosorbent assay and high-performance liquid chromatography methods. Here we report a fast (<30 s), and simple antibody Fc assay based on mix-and-read reporting by fluorescence emission. A soluble fluorescein-labeled Fc-affinity reporter produced by standard peptide synthesis is mixed with an Fc-containing sample to produce an immediate shift in both fluorescence polarization and intensity, compatible with on- and at-line measurements and microbioreactor monitoring. We observed significant shifts in fluorescence intensity in Chinese hamster ovary cell culture fluid spiked with IgG and detected an adalimumab biosimilar down to 100 ng/mL (10-4 g/L), despite the interferents in the complex sample matrix. Neither the fluorescence polarization nor the fluorescence intensity assay is significantly affected by the addition of clarified lysate of 2 million CHO-k1 cells/mL, suggesting applicability even to cultures of low viability. Biochemical and molecular docking approaches suggest that the fluorescence intensity enhancement is caused by changes in the fluorophore's local microenvironment upon binding to IgG Fc, especially by interactions with Fc His433.Abbreviations: CCF: Cell Culture Fluid; CHO: Chinese Hamster Ovary cells; ELISA: Enzyme Linked Immunosorbent Assay; Fc: Fragment Crystallizable of antibody; HPLC: High-Performance Liquid Chromatography; HPßCD: hydroxypropyl-ß-cyclodextrin; IgG: ImmunoglobulinG; mAb: Monoclonal Antibody; PBS: Phosphate-Buffered Saline; PDB: Protein Data Bank; SpA: Staphylococcal protein A; SpG: Staphylococcal protein G.

Screening of Focused Compound Library Targeting Liver X Receptors in Pancreatic Cancer Identified Ligands with Inverse Agonist and Degrader Activity.

Pancreatic ductal adenocarcinoma (PDAC) is the predominant form of pancreatic cancer. PDACs harbor oncogenic mutations in the KRAS gene, and ongoing efforts to directly target its mutant protein product to inhibit tumor growth are a priority not only in pancreatic cancer but in other malignancies such as lung and colorectal cancers where KRAS is also commonly mutated. An alternative strategy to directly targeting KRAS is to identify and target druggable receptors involved in dysregulated cancer hallmarks downstream of KRAS dysregulation. Liver X receptors (LXRs) are members of the nuclear receptor family of ligand-modulated transcription factors and are involved in the regulation of genes which function in key cancer-related processes, including cholesterol transport, lipid and glucose metabolism, and inflammatory and immune responses. Modulation of LXRs via small molecule ligands has emerged as a promising approach for directly targeting tumor cells or the stromal and immune cells within the tumor microenvironment. We have previously shown that only one of the two LXR subtypes (LXRß) is expressed in pancreatic cancer cells, and targeting LXR with available synthetic ligands blocked the proliferation of PDAC cells and tumor formation. In a screen of a focused library of drug-like small molecules predicted to dock in the ligand-binding pocket of LXRß, we identified two novel LXR ligands with more potent antitumor activity than current LXR agonists used in our published studies. Characterization of the two lead compounds (GAC0001E5 and GAC0003A4) indicates that they function as LXR inverse agonists which inhibit their transcriptional activity. Prolonged treatments with novel ligands further revealed their function as LXR "degraders" which significantly reduced LXR protein levels in all three PDAC cell lines tested. These findings support the utility of these novel inhibitors in basic research on ligand design, allosteric mechanisms, and LXR functions and their potential application as treatments for advanced pancreatic cancer and other recalcitrant malignancies.

Impact of lymphoma-linked Asn11Tyr point mutation on the interaction between Bcl-2 and a BH3 mimetic: Insights from molecular dynamics simulation.

Frequent mutations in the Bcl-2 anti-apoptotic protein are often implicated in diffuse large B-cell lymphoma (DLBCL), a disease profoundly resistant to drugs. Bcl-2-competitive inhibitors, "BH3 mimetics," activate apoptosis by interfering with the interactions between pro-apoptotic BH3 domains and the hydrophobic groove of Bcl-2. The aim of our research is to determine the potential of DLBCL-linked N11Y mutation to facilitate resistance against a "BH3 mimetic" using molecular dynamics simulation. Binding free energy calculations suggest a significant decrease in the binding affinity in the mutant model. In-depth analysis of the models using residue interaction network, dynamic cross-correlation, and free energy landscape approaches reveal that the mutation modifies the conformations of key residues, thereby altering the shape of the hydrophobic groove. This subsequently changes the ligand orientation and counteracts the phenomenon of LB region unwinding, a crucial event observed in the wild-type model. Lowest frequency motions captured by principal component analysis reflect the stretching of the groove for efficient ligand accommodation in the wild-type complex but not in the mutant model. This is the first in silico study that unravels the mechanism of drug resistance induced by a Bcl-2 mutation, which could be of great relevance while designing and tailoring therapeutics.

RNAdemocracy: an ensemble method for RNA secondary structure prediction using consensus scoring.

The identification of RNA secondary structure has been an important tool for the characterization of nucleic acids. Computational structure prediction has been an effective approach toward this end, but improvement of established methods is often slow and reliant on redundant methodology. Here we present a novel consensus scoring approach, created to incorporate inputs from an array of established methods with the goal of producing outputs that contain mutual structures from these programs. This method is implemented in RNAdemocracy, a python program capable of competing with existing methods. This ensemble approach was limited by commonalities in established methods like parameter sourcing, which may lead to agreement error, an unavoidable outcome due to the limit of available RNA structure datasets. The modular construction of RNAdemocracy allows for its easy upgrading and customization to suit user's needs. RNAdemocracy, while capable of accurate predictions, is best suited to guide users to regions of the sequence space that exhibit agreement instead of a totally reliant predictor of structure. It is also capable of grading predictions for potential accuracy by providing a percentage of consensus between contributing methods in the final structure.

Insights into the substrate binding specificity of quorum-quenching acylase PvdQ.

Quorum sensing is a cell to cell signaling mechanism that enables them to coordinate their behaviors in a density-dependent manner mediated by small diffusible signaling molecules, which can control the virulence and biofilm gene expression in many Gram-negative and positive bacteria. N-acyl homoserine lactone acylase PvdQ from human opportunistic pathogen Pseudomonas aeruginosa is a quorum-quenching enzyme that can hydrolyze the amide bond of the quorum signaling N-acyl homoserine lactones (AHLs) thereby degrading the signaling molecules, turning off the biofilm phenotype and resulting in a reduction of bacterial virulence. Previous studies demonstrated that PvdQ has different preferences for N-acyl substrates with different acyl chain lengths and substituents. However, the substrate binding specificity determinants of the quorum-quenching enzyme PvdQ with the different bacterial ligands are unknown and unintuitive. Further, elucidation of these determinants can lead to mutants with efficiency and broader substrate promiscuity. To investigate this question, a computational study was carried out combining multiple molecular docking methods, molecular dynamics simulations, residue interaction network analysis, and binding free energy calculations. The main findings are: firstly, the results from pKa predictions support that the pKa of the N-terminus of Serß1 was depressed due to the surrounding residues. Multiple molecular docking studies provide useful information about the detailed binding modes and binding affinities. Secondly, 300 ns molecular dynamics simulations were carried out to analyze the overall molecular motions of substrate-bound and substrate-free PvdQ. The specific interactions between the active site of PvdQ and different ligands revealed the determinants for the preference among the ligands. A systematic comparison and analysis of the protein dynamic fingerprint of each complex demonstrated that binding of the most favorable ligand, C12-homoserine lactone (C12-HSL), reduced the global motions of the complex and maintained the correct arrangement of the catalytic site. Further, the residue interaction network analysis of each system illustrated that there are more communication contacts and pathways between the residues in the C12-HSL complex as compared to complexes with the other ligands. The binding of the C12-HSL ligand facilitates structural communication between the two knobs and the active site. While the binding of the other ligands tend to impair specific communication pathways between the two knobs and the active site, and lead to a catalytically inefficient state. Finally, simulation results from free energy landscape and binding free energy analysis revealed that the C12-HSL ligand has the lowest binding free energy and greater stability than the less favored ligands. Each of the following residues: Serß1, Hisß23, Pheß24, Metß30, Pheß32, Leuß50, Asnß57, Thrß69, Valß70, Trpß162, Trpß186, Asnß269, Argß297 and Leuα146, play different roles in substrate binding specificity. This is the first computational study that provides molecular information for structure-dynamic-function relationships of PvdQ with different ligands and demonstrates determinants of bacterial substrate binding specificity.

Immobilization and unbinding investigation of the antigen-antibody complex using theoretical and experimental techniques.

Experimental results for the antibody known as immunoglobulin G - IgG interacting with phenobarbital were obtained via atomic force microscopy (AFM) and thereafter investigated using computer simulation modeling tools. Using molecular dynamics simulation and docking calculations, the energetically stable configurations of an immobilized antibody over a silicon surface were searched. Six stable configurations of the immobilized antibody over the silicon nitride surface covered by linker molecules were found. Although, only three of them (P1, P2, P5) maintained the Fragment antigen binding available for antigen interaction. Therefore, these configurations were equilibrated after reaching 100 ns molecular dynamics trajectory. The average interaction energy between the surface and the immunoglobulin G - IgG antibody in the P1, P2 and P5 configurations were -62.4 ±â€¯2.4 kcal/mol; -54.3 ±â€¯5.7 kcal/mol, and -360.9 ±â€¯4.2 kcal/mol respectively. Phenobarbital was docked within the Fab domain of P1, P2, and P5 immobilized configurations and equilibrated with molecular dynamics for binding energy estimation. Then, steered molecular dynamics was performed to evaluate unbinding energy pathway between phenobarbital and IgG in each of the three-oriented IgG configurations. No significant differences were observed in the rupture force values (EP1 = 591 ±â€¯13 pN, EP2 = 605 ±â€¯18 pN, and EP5 = 610 ±â€¯45 pN). In comparison, the average AFM experimental results were (641.6 ±â€¯363.3 pN). Therefore, it is worth noting that P5 is the configuration with highest protein-surface interaction. Therefore, the force value calculated for the P5 orientation is statistically more favorable and it is the one to be compared to the experimental data. The agreement between experimental and theoretical results indicates a favorable presented for this study opening new perspectives for antigen-antibody evaluation.

Targeted reduction of the EGFR protein, but not inhibition of its kinase activity, induces mitophagy and death of cancer cells through activation of mTORC2 and Akt.

The oncogenic epidermal growth factor receptor (EGFR) is commonly overexpressed in solid cancers. The tyrosine kinase activity of EGFR has been a major therapeutic target for cancer; however, the efficacy of EGFR tyrosine kinase inhibitors to treat cancers has been challenged by innate and acquired resistance at the clinic. Accumulating evidence suggests that EGFR possesses kinase-independent pro-survival functions, and that cancer cells are more vulnerable to reduction of EGFR protein than to inhibition of its kinase activity. The molecular mechanism underlying loss-of-EGFR-induced cell death remains largely unknown. In this study, we show that, unlike inhibiting EGFR kinase activity that is known to induce pro-survival non-selective autophagy, downregulating EGFR protein, either by siRNA, or by a synthetic EGFR-downregulating peptide (Herdegradin), kills prostate and ovarian cancer cells via selective mitophagy by activating the mTORC2/Akt axis. Furthermore, Herdegradin induced mitophagy and inhibited the growth of orthotopic ovarian cancers in mice. This study identifies anti-mitophagy as a kinase-independent function of EGFR, reveals a novel function of mTORC2/Akt axis in promoting mitophagy in cancer cells, and offers a novel approach for pharmacological downregulation of EGFR protein as a potential treatment for EGFR-positive cancers.

Inhibition of Cholera Toxin and Other AB Toxins by Polyphenolic Compounds.

Cholera toxin (CT) is an AB-type protein toxin that contains a catalytic A1 subunit, an A2 linker, and a cell-binding B homopentamer. The CT holotoxin is released into the extracellular environment, but CTA1 attacks a target within the cytosol of a host cell. We recently reported that grape extract confers substantial resistance to CT. Here, we used a cell culture system to identify twelve individual phenolic compounds from grape extract that inhibit CT. Additional studies determined the mechanism of inhibition for a subset of the compounds: two inhibited CT binding to the cell surface and even stripped CT from the plasma membrane of a target cell; two inhibited the enzymatic activity of CTA1; and four blocked cytosolic toxin activity without directly affecting the enzymatic function of CTA1. Individual polyphenolic compounds from grape extract could also generate cellular resistance to diphtheria toxin, exotoxin A, and ricin. We have thus identified individual toxin inhibitors from grape extract and some of their mechanisms of inhibition against CT.

Structural mutation analysis of PTEN and its genotype-phenotype correlations in endometriosis and cancer.

The phosphatase and tensin homolog deleted on chromosome ten (PTEN) gene encodes a tumor suppressor phosphatase that has recently been found to be frequently mutated in patients with endometriosis, endometrial cancer and ovarian cancer. Here, we present the first computational analysis of 13 somatic missense PTEN mutations associated with these phenotypes. We found that a majority of the mutations are associated in conserved positions within the active site and are clustered within the signature motif, which contain residues that play a crucial role in loop conformation and are essential for catalysis. In silico analyses were utilized to identify the putative effects of these mutations. In addition, coarse-grained models of both wild-type (WT) PTEN and mutants were constructed using elastic network models to explore the interplay of the structural and global dynamic effects that the mutations have on the relationship between genotype and phenotype. The effects of the mutations reveal that the local structure and interactions affect polarity, protein structure stability, electrostatic surface potential, and global dynamics of the protein. Our results offer new insight into the role in which PTEN missense mutations contribute to the molecular mechanism and genotypic-phenotypic correlation of endometriosis, endometrial cancer, and ovarian cancer. Proteins 2016; 84:1625-1643. © 2016 Wiley Periodicals, Inc.

Functional Implications of the spectrum of BCL2 mutations in Lymphoma.

Mutations in the translocated BCL2 gene are often detected in diffuse large B-cell lymphomas (DLBCLs), indicating both their significance and pervasiveness. Large series genome sequencing of more than 200 DLBCLs has identified frequent BCL2 mutations clustered in the exons coding for the BH4 domain and the folded loop domain (FLD) of the protein. However, BCL2 mutations are mostly contemplated to represent bystander events with negligible functional impact on the pathogenesis of DLBCL. BCL2 arbitrates apoptosis through a classic interaction between its hydrophobic groove forming BH1-3 domains and the BH3 domain of pro-apoptotic members of the BCL2 family. The effects of mutations are mainly determined by the ability of the mutated BCL2 to mediate apoptosis by this inter-member protein binding. Nevertheless, BCL2 regulates diverse non-canonical pathways that are unlikely to be explained by canonical interactions. In this review, first, we identify recurrent missense mutations in the BH4 domain and the FLD reported in independent lymphoma sequencing studies. Second, we discuss the probable consequences of mutations on the binding ability of BCL2 to non-BCL2 family member proteins crucial for 1) maintaining mitochondrial energetics and calcium hemostasis such as VDAC, IP3R, and RyR and 2) oncogenic pathways implicated in the acquisition of the 'hallmarks of cancer' such as SOD, Raf-1, NFAT, p53, HIF-1α, and gelsolin. The study also highlights the likely ramifications of mutations on binding of BCL2 antagonists and BH3 profiling. Based on our analysis, we believe that an in-depth focus on BCL2 interactions mediated by these domains is warranted to elucidate the functional significance of missense mutations in DLBCL. In summary, we provide an extensive overview of the pleiotropic functions of BCL2 mediated by its physical binding interaction with other proteins and the various ways BCL2 mutations would affect the normal function of the cell leading to the development of DLBCL.

Biophysical characteristics of cholera toxin and Escherichia coli heat-labile enterotoxin structure and chemistry lead to differential toxicity.

The biophysical chemistry of macromolecular complexes confer their functional characteristics. We investigate the mechanisms that make the AB5 holotoxin of Vibrio cholerae (CT) a significantly more pathogenic molecule than the enterotoxin of Escherichia coli (LT) with which it shares 88% similarity and whose structure is homologous with a backbone RMSD of 0.84 Å and imposes its deleterious effects though the same process to constitutively ADP-ribosylate adenylate cyclase. We present computational data that characterizes the impact of amino acid variations in the A2 tail, which helps to explain experimental data that demonstrate CT's higher toxicity. A hydrophobic patch on the B pentamer interface and its interactions with the A subdomain are partially disrupted by the substitution of an aspartic acid (LT) for glycine in CT. CT's holotoxin has less solvent accessible surface area (94 Å(2) vs 54 Å(2)) and higher contact area (280 Å(2) vs 241 Å(2)) with S228, which is a gatekeeper, partially controlling the diffusion of water into the pore. CT excludes water from the top of the central pore whereas LT allows much more water to interact. These biophysical properties of the toxins lead to their differential toxicity and resulting impact to human health.

Associations between the Rho kinase-1 catalytic and PH domain regulatory unit.

Rho-associated kinase, or ROCK, is an important mediator of ventricular remodeling in cardiac hypertrophy. It has a kinase catalytic domain, a coiled-coil domain and a Pleckstrin-Homology domain (PH domain) with a C1 domain insert. The C-terminal region including the PH domain and C1 domain insert is involved in an autoregulatory role for ROCK. We sought to evaluate whether a self association complex could form using computational docking approaches. We found that both the PH domain and the C1 domain could dock with the catalytic domain and we further found that they could dock in poses that are complementary to each other forming a three domain complex. We also confirmed a binding response using a surface plasmon resonance experimental approach. Information about the regulation of ROCK might lead to new strategies to develop lead inhibitor compounds to modulate cardiac remodeling.

Molecular modeling of enzyme attachment on AFM probes.

The immobilization of enzymes on atomic force microscope tip (AFM tip) surface is a crucial step in the development of nanobiosensors to be used in detection process. In this work, an atomistic modeling of the attachment of the acetyl coenzyme A carboxylase (ACC enzyme) on a functionalized AFM tip surface is proposed. Using electrostatic considerations, suitable enzyme-surface orientations with the active sites of the ACC enzyme available for interactions with bulk molecules were found. A 50 ns molecular dynamics trajectory in aqueous solution was obtained and surface contact area, hydrogen bonding and protein stability were analyzed. The enzyme-surface model proposed here with minor adjustment can be applied to study antigen-antibody interactions as well as enzyme immobilization on silica for chromatography applications.

New insights into the binding and catalytic mechanisms of Bacillus thuringiensis lactonase: insights into B. thuringiensis AiiA mechanism.

The lactonase enzyme (AiiA) produced by Bacillus thuringiensis serves to degrade autoinducer-1 (AI-1) signaling molecules in what is an evolved mechanism by which to compete with other bacteria. Bioassays have been previously performed to determine whether the AI-1 aliphatic tail lengths have any effect on AiiA's bioactivity, however, data to date are conflicting. Additionally, specific residue contributions to the catalytic activity of AiiA provide for some interesting questions. For example, it has been proposed that Y194 serves to provide an oxyanion hole to AI-1 which is curious given the fact the substrate spans two Zn(2+) ions. These ions might conceivably provide enough charge to promote both ligand stability and the carbonyl activation necessary to drive a nucleophilic attack. To investigate these questions, multiple molecular dynamics simulations were performed across a family of seven acylated homoserine lactones (AHL) along with their associated intermediate and product states. Distance analyses and interaction energy analyses were performed to investigate current bioassay data. Our simulations are consistent with experimental studies showing that AiiA degrades AHLs in a tail length independent manner. However, the presence of the tail is required for activity. Also, the putative oxyanion hole function of Y194 toward the substrate is not observed in any of the reactant or product state simulation trajectories, but does seem to show efficacy in stabilizing the intermediate state. Last, we argue through ionization state analyses, that the proton shuttling necessary for catalytic activity might be mediated by both water and substrate-based intra-molecular proton transfer. Based on this argument, an alternate catalytic mechanism is proposed.

Interactions of acetylcholine binding site residues contributing to nicotinic acetylcholine receptor gating: role of residues Y93, Y190, K145 and D200.

The nicotinic acetylcholine receptor exhibits multiple conformational states, resting (channel closed), active (channel open) and desensitized (channel closed). The resting state may be distinguished from the active and desensitized states by the orientation of loop C in the extracellular ligand binding domain (LBD). Homology modeling was used to generate structures of the Torpedo californica α2ßδγ nAChR that initially represent the resting state (loop C open) and the desensitized state (loop C closed). Molecular dynamics (MD) simulations were performed on the extracellular LBD on each nAChR conformational state, with and without the agonist anabaseine present in each binding site (the αγ and the αδ sites). Three MD simulations of 10ns each were performed for each of the four conditions. Comparison of dynamics revealed that in the presence of agonist, loop C was drawn inward and attains a more stable conformation. Examination of side-chain interactions revealed that residue αY190 exhibited hydrogen-bonding interactions either with residue αY93 in the ligand binding site or with residue αK145 proximal to the binding site. αK145 also exhibited side chain (salt bridge) interactions with αD200 and main chain interactions with αY93. Residues αW149, αY198, γY116/δT119, γL118/δL121 and γL108/δL111 appear to play the role of stabilizing ligand in the binding site. In MD simulations for the desensitized state, the effect of ligand upon the interactions among αK145, αY190, and αY93 as well as ligand-hydrogen-bonding to αW149 were more pronounced at the αγ interface than at the αδ interface. Differences in affinity for the desensitized state were determined experimentally to be 10-fold. The changes in side chain interactions observed for the two conformations and induced by ligand support a model wherein hydrogen bond interactions between αD200 and αY93 are broken and rearrange to form a salt-bridge between αK145 and αD200 and hydrogen bond interactions between αY93 and αY190 and between αK145 and αY190.