Identification and Structure-Activity Relationships of Novel Compounds that Potentiate the Activities of Antibiotics in Escherichia coli

In Gram-negative bacteria, efflux pumps are able to prevent effective cellular concentrations from being achieved for a number of antibiotics. Small molecule adjuvants that act as efflux pump inhibitors (EPIs) have the potential to reinvigorate existing antibiotics that are currently ineffective due to efflux mechanisms. Through a combination of rigorous experimental screening and in silico virtual screening, we recently identified novel classes of EPIs that interact with the membrane fusion protein AcrA, a critical component of the AcrAB-TolC efflux pump in Escherichia coli. Herein, we present initial optimization efforts and structure-activity relationships around one of those previously described hits, NSC 60339 (1). From these efforts we identified two compounds, SLUPP-225 (17h) and SLUPP-417 (17o), which demonstrate favorable properties as potential EPIs in E. coli cells including the ability to penetrate the outer membrane, improved inhibition of efflux relative to 1, and potentiation of the activity of novobiocin and erythromycin.

The tilt-dependent potential of mean force of a pair of DNA oligomers from all-atom molecular dynamics simulations

Electrostatic interactions between DNA molecules have been extensively studied experimentally and theoretically, but several aspects (e.g. its role in determining the pitch of the cholesteric DNA phase) still remain unclear. Here, we performed large-scale all-atom molecular dynamics simulations in explicit water and 150 mM sodium chloride, to reconstruct the potential of mean force (PMF) of two DNA oligomers 24 base pairs long as a function of their interaxial angle and intermolecular distance. We find that the potential of mean force is dominated by total DNA charge, and not by the helical geometry of its charged groups. The theory of homogeneously charged cylinders fits well all our simulation data, and the fit yields the optimal value of the total compensated charge on DNA to  ≈65% of its total fixed charge (arising from the phosphorous atoms), close to the value expected from Manning's theory of ion condensation. The PMF calculated from our simulations does not show a significant dependence on the handedness of the angle between the two DNA molecules, or its size is on the order of [Formula: see text]. Thermal noise for molecules of the studied length seems to mask the effect of detailed helical charge patterns of DNA. The fact that in monovalent salt the effective interaction between two DNA molecules is independent on the handedness of the tilt may suggest that alternative mechanisms are required to understand the cholesteric phase of DNA.

Abstract P4-22-11: Combined targeted therapies for advanced triple negative breast cancer: A phase II trial of nab-paclitaxel and bevacizumab followed by maintenance targeted therapy with bevacizumab and erlotinib

Background: Chemotherapy remains the mainstay of therapy for patients with metastatic triple negative breast cancer (TNBC). We hypothesized that the addition of biologic agents targeting key pathways (bevacizumab targeting angiogenesis and erlotinib directed against EGFR) may prolong progression free survival (PFS) and offer a novel treatment strategy free from chemotherapy for patients with metastatic TNBC.

Methods: Patients with TNBC receiving initial therapy for metastatic disease were eligible for this multicenter phase II trial (NCT00733408) conducted at an academic center and affiliated, community practice sites. Induction therapy included nab-paclitaxel 100 mg/m2 IV Qweek (wk) and Bevacizumab 10 mg/kg IV Q2wks x 24 weeks. Patients free of progression at 24 wks began maintenance therapy with bevacizumab 10 mg/kg IV Q2wks and erlotinib 150 mg po daily until progression with radiographic assessment every 8 wks. Primary objective was PFS with secondary objectives of response rate, overall survival (OS) and safety. All eligible patients were included in the analysis of PFS and OS. Response was evaluated among patients with measurable disease by RECIST 1.1 with central review. Patients with inadequate disease assessments were coded as non-responders. Kaplan-Meier method was used to estimate PFS and OS with patients censored at date of last tumor assessment (PFS) or date of last follow up (OS).

Results: From April 2009 – December 2015, 58 patients (median age 54, range 33-83) were enrolled; 56 (97%) had measurable disease, and all had metastatic TNBC by local assessment. 33 (57%) patients completed induction; 22 (38%) came off study during induction; 3 (5%) continue on maintenance therapy. 4 patients discontinued therapy prior to first assessment. As of June 8, 2016, 53 patients (91%) have progressed. Median follow up for surviving patients is 14.5 months (range 4.1-65.4). Median PFS is 7.7 months (95% CI 5.7, 9.5). Of 56 patients with measurable disease, 38 (66%) had partial response (PR); 10 (17%) with stable disease for clinical benefit rate (CBR) of 86%. Median OS is 18.2 months (95% CI 16.3, 24.5). Most common grade 3-4 toxicities during induction were neutropenia [17 (29%), 1 grade 4], fatigue [13 (22%), all grade 3], leukopenia [7 (12%), all grade 3], and neuropathy [7 (12%), all grade 3]. Rash was most common ≥ grade 3 toxicity during maintenance [4 (7%), grade 3]. One patient experienced clinical CHF during maintenance month 16 requiring bevacizumab discontinuation. Conclusions: Nab-paclitaxel and bevacizumab followed by maintenance targeted therapy with bevacizumab and erlotinib was well tolerated. While the observed PFS did not meet pre-specified criteria of interest, the majority of patients experienced clinical benefit (86%) with 30 (57%) receiving maintenance targeted therapy. Correlative studies are ongoing. Supported by Genentech (OSI4266s), Celegene (AX-CL-BRST-PI-003828) and Janssen.

Citation Format: Specht JM, Gadi VK, Gralow JR, Korde LA, Linden HM, Salazar LG, Rodler ET, Cundy A, Buening BJ, Baker KK, Redman MW, Kurland BF, Garrison MA, Smith JC, vanHaelst C, Anderson JE. Combined targeted therapies for advanced triple negative breast cancer: A phase II trial of nab-paclitaxel and bevacizumab followed by maintenance targeted therapy with bevacizumab and erlotinib [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P4-22-11.

Reviving Antibiotics: Efflux Pump Inhibitors That Interact with AcrA, a Membrane Fusion Protein of the AcrAB-TolC Multidrug Efflux Pump

Antibiotic resistance is a major threat to human welfare. Inhibitors of multidrug efflux pumps (EPIs) are promising alternative therapeutics that could revive activities of antibiotics and reduce bacterial virulence. Identification of new druggable sites for inhibition is critical for the development of effective EPIs, especially in light of constantly emerging resistance. Here, we describe EPIs that interact with periplasmic membrane fusion proteins, critical components of efflux pumps that are responsible for the activation of the transporter and the recruitment of the outer-membrane channel. The discovered EPIs bind to AcrA, a component of the prototypical AcrAB-TolC pump, change its structure in vivo, inhibit efflux of fluorescent probes, and potentiate the activities of antibiotics in Escherichia coli and other Gram-negative bacteria. Our findings expand the chemical and mechanistic diversity of EPIs, suggest the mechanism for regulation of the efflux pump assembly and activity, and provide a promising path for reviving the activities of antibiotics in resistant bacteria.

Structure-based design of broadly protective group a streptococcal M protein-based vaccines

BACKGROUND

A major obstacle to the development of broadly protective M protein-based group A streptococcal (GAS) vaccines is the variability within the N-terminal epitopes that evoke potent bactericidal antibodies. The concept of M type-specific protective immune responses has recently been challenged based on the observation that multivalent M protein vaccines elicited cross-reactive bactericidal antibodies against a number of non-vaccine M types of GAS. Additionally, a new "cluster-based" typing system of 175M proteins identified a limited number of clusters containing closely related M proteins. In the current study, we used the emm cluster typing system, in combination with computational structure-based peptide modeling, as a novel approach to the design of potentially broadly protective M protein-based vaccines.

METHODS

M protein sequences (AA 16-50) from the E4 cluster containing 17 emm types of GAS were analyzed using de novo 3-D structure prediction tools and the resulting structures subjected to chemical diversity analysis to identify sequences that were the most representative of the 3-D physicochemical properties of the M peptides in the cluster. Five peptides that spanned the range of physicochemical attributes of all 17 peptides were used to formulate synthetic and recombinant vaccines. Rabbit antisera were assayed for antibodies that cross-reacted with E4 peptides and whole bacteria by ELISA and for bactericidal activity against all E4GAS.

RESULTS

The synthetic vaccine rabbit antisera reacted with all 17 E4M peptides and demonstrated bactericidal activity against 15/17 E4GAS. A recombinant hybrid vaccine containing the same E4 peptides also elicited antibodies that cross-reacted with all E4M peptides.

CONCLUSIONS

Comprehensive studies using structure-based design may result in a broadly protective M peptide vaccine that will elicit cluster-specific and emm type-specific antibody responses against the majority of clinically relevant emm types of GAS.

A computationally identified compound antagonizes excess FGF-23 signaling in renal tubules and a mouse model of hypophosphatemia

Fibroblast growth factor-23 (FGF-23) interacts with a binary receptor complex composed of α-Klotho (α-KL) and FGF receptors (FGFRs) to regulate phosphate and vitamin D metabolism in the kidney. Excess FGF-23 production, which causes hypophosphatemia, is genetically inherited or occurs with chronic kidney disease. Among other symptoms, hypophosphatemia causes vitamin D deficiency and the bone-softening disorder rickets. Current therapeutics that target the receptor complex have limited utility clinically. Using a computationally driven, structure-based, ensemble docking and virtual high-throughput screening approach, we identified four novel compounds predicted to selectively inhibit FGF-23-induced activation of the FGFR/α-KL complex. Additional modeling and functional analysis found that Zinc13407541 bound to FGF-23 and disrupted its interaction with the FGFR1/α-KL complex; experiments in a heterologous cell expression system showed that Zinc13407541 selectivity inhibited α-KL-dependent FGF-23 signaling. Zinc13407541 also inhibited FGF-23 signaling in isolated renal tubules ex vivo and partially reversed the hypophosphatemic effects of excess FGF-23 in a mouse model. These chemical probes provide a platform to develop lead compounds to treat disorders caused by excess FGF-23.

A Distal Disulfide Bridge in OXA-1 β-Lactamase Stabilizes the Catalytic Center and Alters the Dynamics of the Specificity Determining Ω Loop

Widespread antibiotic resistance, particularly when mediated by broad-spectrum β-lactamases, has major implications for public health. Substitutions in the active site often allow broad-spectrum enzymes to accommodate diverse types of β-lactams. Substitutions observed outside the active site are thought to compensate for the loss of thermal stability. The OXA-1 clade of class D β-lactamases contains a pair of conserved cysteines located outside the active site that forms a disulfide bond in the periplasm. Here, the effect of the distal disulfide bond on the structure and dynamics of OXA-1 was investigated via 4 μs molecular dynamics simulations. The results reveal that the disulfide promotes the preorganized orientation of the catalytic residues and affects the conformation of the functionally important Ω loop. Furthermore, principal component analysis reveals differences in the global dynamics between the oxidized and reduced forms, especially in the motions involving the Ω loop. A dynamical network analysis indicates that, in the oxidized form, in addition to its role in ligand binding, the KTG family motif is a central hub of the global dynamics. As activity of OXA-1 has been measured only in the reduced form, we suggest that accurate assessment of its functional profile would require oxidative conditions mimicking periplasm.

Determination of functional collective motions in a protein at atomic resolution using coherent neutron scattering

Protein function often depends on global, collective internal motions. However, the simultaneous quantitative experimental determination of the forms, amplitudes, and time scales of these motions has remained elusive. We demonstrate that a complete description of these large-scale dynamic modes can be obtained using coherent neutron-scattering experiments on perdeuterated samples. With this approach, a microscopic relationship between the structure, dynamics, and function in a protein, cytochrome P450cam, is established. The approach developed here should be of general applicability to protein systems.

Neutron structure of human carbonic anhydrase II in complex with methazolamide: mapping the solvent and hydrogen-bonding patterns of an effective clinical drug

Carbonic anhydrases (CAs; EC 4.2.1.1) catalyze the interconversion of CO2 and HCO3-, and their inhibitors have long been used as diuretics and as a therapeutic treatment for many disorders such as glaucoma and epilepsy. Acetazolamide (AZM) and methazolamide (MZM, a methyl derivative of AZM) are two of the classical CA inhibitory drugs that have been used clinically for decades. The jointly refined X-ray/neutron structure of MZM in complex with human CA isoform II (hCA II) has been determined to a resolution of 2.2 Å with an Rcryst of ∼16.0%. Presented in this article, along with only the second neutron structure of a clinical drug-bound hCA, is an in-depth structural comparison and analyses of differences in hydrogen-bonding network, water-molecule orientation and solvent displacement that take place upon the binding of AZM and MZM in the active site of hCA II. Even though MZM is slightly more hydrophobic and displaces more waters than AZM, the overall binding affinity (Ki) for both of the drugs against hCA II is similar (∼10 nM). The plausible reasons behind this finding have also been discussed using molecular dynamics and X-ray crystal structures of hCA II-MZM determined at cryotemperature and room temperature. This study not only allows a direct comparison of the hydrogen bonding, protonation states and solvent orientation/displacement of AZM and MZM, but also shows the significant effect that the methyl derivative has on the solvent organization in the hCA II active site.

The unconscious, myth, and the rule of law: Reflections on the persistence of gender inequality

Social order, to remain stable, needs the voluntary compliance of the majority of the population. Such consent requires normative justification. The rational foundation of the rule of law and the democratic state rests on the presumption of the equality of every citizen. Male domination of females nevertheless remains universal even in the most advanced democratic nation states because it is legitimized by the shared assumption that patriarchy reflects the will of God or is dictated by nature. Freud's diagnosis of patriarchy as a collective neurosis of the group mind negates every possible normative justification that can be made for gender hierarchy. Freud made extensive references to myth in developing his analysis of the neurotic foundations of social order. An analysis of the structure of myth suggests that ideological seduction rather than God, nature or biology determines male dominance.

Stress Symptoms One Year after 9/11/01: A Follow-up

This study investigated whether stress symptoms related to attention deficit previously found for a sample ( n = 149) of Chicago-area college students 1 to 5 weeks after the September 11th attacks were present in a comparable sample ( n = 129) 12 to 14 months later. As hypothesized, the later participants had a significantly lower mean score on the scale, Attention Deficit of the Smith Stress Symptoms Inventory, than a different group of participants tested immediately after the attacks. Mean Attention Deficit scores for groups after one year and immediately before September 11th ( n = 320) did not differ statistically. Findings are consistent with the interpretation that elevations in scores on Attention Deficit immediately after September 11th were indeed associated with the catastrophic attacks and not artifacts of time of examination. Consistent with previous research with this inventory, disaster-related traumatic stress scores are lower over time, although here no control group was included.

Ensemble-based docking: From hit discovery to metabolism and toxicity predictions

This paper describes and illustrates the use of ensemble-based docking, i.e., using a collection of protein structures in docking calculations for hit discovery, the exploration of biochemical pathways and toxicity prediction of drug candidates. We describe the computational engineering work necessary to enable large ensemble docking campaigns on supercomputers. We show examples where ensemble-based docking has significantly increased the number and the diversity of validated drug candidates. Finally, we illustrate how ensemble-based docking can be extended beyond hit discovery and toward providing a structural basis for the prediction of metabolism and off-target binding relevant to pre-clinical and clinical trials.

Toward Quantitatively Accurate Calculation of the Redox-Associated Acid-Base and Ligand Binding Equilibria of Aquacobalamin

Redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. A major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co-ligand binding equilibrium constants (Kon/off), pKas, and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co(III), Co(II), and Co(I) species, respectively, and the second model features saturation of each vacant axial coordination site on Co(II) and Co(I) species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pKas and 2.3 log units for two log Kon/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co-axial ligand binding, leading to substantial errors in predicted pKas and Kon/off values. These findings demonstrate the effectiveness of the present approach for computing electrochemical and thermodynamic properties of a complex transition metal-containing cofactor.

Quasielastic neutron scattering in biology: Theory and applications

Neutrons scatter quasielastically from stochastic, diffusive processes, such as overdamped vibrations, localized diffusion and transitions between energy minima. In biological systems, such as proteins and membranes, these relaxation processes are of considerable physical interest. We review here recent methodological advances and applications of quasielastic neutron scattering (QENS) in biology, concentrating on the role of molecular dynamics simulation in generating data with which neutron profiles can be unambiguously interpreted. We examine the use of massively-parallel computers in calculating scattering functions, and the application of Markov state modeling. The decomposition of MD-derived neutron dynamic susceptibilities is described, and the use of this in combination with NMR spectroscopy. We discuss dynamics at very long times, including approximations to the infinite time mean-square displacement and nonequilibrium aspects of single-protein dynamics. Finally, we examine how neutron scattering and MD can be combined to provide information on lipid nanodomains. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

The role of histone tails in the nucleosome: a computational study

Histone tails play an important role in gene transcription and expression. We present here a systematic computational study of the role of histone tails in the nucleosome, using replica exchange molecular dynamics simulations with an implicit solvent model and different well-established force fields. We performed simulations for all four histone tails, H4, H3, H2A, and H2B, isolated and with inclusion of the nucleosome. The results confirm predictions of previous theoretical studies for the secondary structure of the isolated tails but show a strong dependence on the force field used. In the presence of the entire nucleosome for all force fields, the secondary structure of the histone tails is destabilized. Specific contacts are found between charged lysine and arginine residues and DNA phosphate groups and other binding sites in the minor and major DNA grooves. Using cluster analysis, we found a single dominant configuration of binding to DNA for the H4 and H2A histone tails, whereas H3 and H2B show multiple binding configurations with an equal probability. The leading stabilizing contribution for those binding configurations is the attractive interaction between the positively charged lysine and arginine residues and the negatively charged phosphate groups, and thus the resulting charge neutralization. Finally, we present results of molecular dynamics simulations in explicit solvent to confirm our conclusions. Results from both implicit and explicit solvent models show that large portions of the histone tails are not bound to DNA, supporting the complex role of these tails in gene transcription and expression and making them possible candidates for binding sites of transcription factors, enzymes, and other proteins.

Evidence for Osteocalcin Binding and Activation of GPRC6A in β-Cells

The possibility that G protein-coupled receptor family C member A (GPRC6A) is the osteocalcin (Ocn)-sensing G protein-coupled receptor that directly regulates pancreatic β-cell functions is controversial. In the current study, we found that Ocn and an Ocn-derived C-terminal hexapeptide directly activate GPRC6A-dependent ERK signaling in vitro. Computational models probe the structural basis of Ocn binding to GPRC6A and predict that the C-terminal hexapeptide docks to the extracellular side of the transmembrane domain of GPRC6A. Consistent with the modeling, mutations in the computationally identified binding pocket of GPRC6A reduced Ocn and C-terminal hexapeptide activation of this receptor. In addition, selective deletion of Gprc6a in β-cells (Gprc6a(β)(-cell-cko)) by crossing Gprc6a(flox/flox) mice with Ins2-Cre mice resulted in reduced pancreatic weight, islet number, insulin protein content, and insulin message expression. Both islet size and β-cell proliferation were reduced in Gprc6a(β)(-cell-cko) compared with control mice. Gprc6a(β)(-cell-cko) exhibited abnormal glucose tolerance, but normal insulin sensitivity. Islets isolated from Gprc6a(β)(-cell-cko) mice showed reduced insulin simulation index in response to Ocn. These data establish the structural basis for Ocn direct activation of GPRC6A and confirm a role for GPRC6A in regulating β-cell proliferation and insulin secretion.

Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.

Conformations of Low-Molecular-Weight Lignin Polymers in Water

Low-molecular-weight lignin binds to cellulose during the thermochemical pretreatment of biomass for biofuel production, which prevents the efficient hydrolysis of the cellulose to sugars. The binding properties of lignin are influenced strongly by the conformations it adopts. Here, we use molecular dynamics simulations in aqueous solution to investigate the dependence of the shape of lignin polymers on chain length and temperature. Lignin is found to adopt collapsed conformations in water at 300 and 500 K. However, at 300 K, a discontinuous transition is found in the shape of the polymer as a function of the chain length. Below a critical degree of polymerization, Nc =15, the polymer adopts less spherical conformations than above Nc. The transition disappears at high temperatures (500 K) at which only spherical shapes are adopted. An implication relevant to cellulosic biofuel production is that lignin will self-aggregate even at high pretreatment temperatures.

General trends of dihedral conformational transitions in a globular protein

Dihedral conformational transitions are analyzed systematically in a model globular protein, cytochrome P450cam, to examine their structural and chemical dependences through combined conventional molecular dynamics (cMD), accelerated molecular dynamics (aMD) and adaptive biasing force (ABF) simulations. The aMD simulations are performed at two acceleration levels, using dihedral and dual boost, respectively. In comparison with cMD, aMD samples protein dihedral transitions approximately two times faster on average using dihedral boost, and ∼ 3.5 times faster using dual boost. In the protein backbone, significantly higher dihedral transition rates are observed in the bend, coil, and turn flexible regions, followed by the β bridge and β sheet, and then the helices. Moreover, protein side chains of greater length exhibit higher transition rates on average in the aMD-enhanced sampling. Side chains of the same length (particularly Nχ = 2) exhibit decreasing transition rates with residues when going from hydrophobic to polar, then charged and aromatic chemical types. The reduction of dihedral transition rates is found to be correlated with increasing energy barriers as identified through ABF free energy calculations. These general trends of dihedral conformational transitions provide important insights into the hierarchical dynamics and complex free energy landscapes of functional proteins.

Enhanced sampling simulation analysis of the structure of lignin in the THF–water miscibility gap

Using temperature replica-exchange molecular dynamics, we characterize a globule-to-coil transition for a softwood-like lignin biopolymer in a tetrahydrofuran (THF)–water cosolvent system at temperatures at which the cosolvent undergoes a de-mixing transition.

A Structural Study of CESA1 Catalytic Domain of Arabidopsis Cellulose Synthesis Complex: Evidence for CESA Trimers

A cellulose synthesis complex with a "rosette" shape is responsible for synthesis of cellulose chains and their assembly into microfibrils within the cell walls of land plants and their charophyte algal progenitors. The number of cellulose synthase proteins in this large multisubunit transmembrane protein complex and the number of cellulose chains in a microfibril have been debated for many years. This work reports a low resolution structure of the catalytic domain of CESA1 from Arabidopsis (Arabidopsis thaliana; AtCESA1CatD) determined by small-angle scattering techniques and provides the first experimental evidence for the self-assembly of CESA into a stable trimer in solution. The catalytic domain was overexpressed in Escherichia coli, and using a two-step procedure, it was possible to isolate monomeric and trimeric forms of AtCESA1CatD. The conformation of monomeric and trimeric AtCESA1CatD proteins were studied using small-angle neutron scattering and small-angle x-ray scattering. A series of AtCESA1CatD trimer computational models were compared with the small-angle x-ray scattering trimer profile to explore the possible arrangement of the monomers in the trimers. Several candidate trimers were identified with monomers oriented such that the newly synthesized cellulose chains project toward the cell membrane. In these models, the class-specific region is found at the periphery of the complex, and the plant-conserved region forms the base of the trimer. This study strongly supports the "hexamer of trimers" model for the rosette cellulose synthesis complex that synthesizes an 18-chain cellulose microfibril as its fundamental product.

AM1/d Parameters for Magnesium in Metalloenzymes

AM1/d parameters are derived for magnesium, optimized for modeling reactions in metalloenzymes. The parameters are optimized with a Monte Carlo procedure so as to reproduce the geometries and energies of a training set calculated with density functional theory. The training set consists of compounds with magnesium coordinated to the oxygen atom of typical biological ligands. Optimization of AM1 parameters without extension to d functions leaves serious errors. The new AM1/d parameters provide a clear improvement in accuracy compared to the standard semiempirical methods AM1 and MNDO/d and will be particularly useful for modeling reactions in large biological systems at low computational cost.

Scaling of Multimillion-Atom Biological Molecular Dynamics Simulation on a Petascale Supercomputer

A strategy is described for a fast all-atom molecular dynamics simulation of multimillion-atom biological systems on massively parallel supercomputers. The strategy is developed using benchmark systems of particular interest to bioenergy research, comprising models of cellulose and lignocellulosic biomass in an aqueous solution. The approach involves using the reaction field (RF) method for the computation of long-range electrostatic interactions, which permits efficient scaling on many thousands of cores. Although the range of applicability of the RF method for biomolecular systems remains to be demonstrated, for the benchmark systems the use of the RF produces molecular dipole moments, Kirkwood G factors, other structural properties, and mean-square fluctuations in excellent agreement with those obtained with the commonly used Particle Mesh Ewald method. With RF, three million- and five million-atom biological systems scale well up to ∼30k cores, producing ∼30 ns/day. Atomistic simulations of very large systems for time scales approaching the microsecond would, therefore, appear now to be within reach.

Transition Networks for the Comprehensive Characterization of Complex Conformational Change in Proteins

Functionally relevant transitions between native conformations of a protein can be complex, involving, for example, the reorganization of parts of the backbone fold, and may occur via a multitude of pathways. Such transitions can be characterized by a transition network (TN), in which the experimentally determined end state structures are connected by a dense network of subtransitions via low-energy intermediates. We show here how the computation of a TN can be achieved for a complex protein transition. First, an efficient hierarchical procedure is used to uniformly sample the conformational subspace relevant to the transition. Then, the best path which connects the end states is determined as well as the rate-limiting ridge on the energy surface which separates them. Graph-theoretical algorithms permit this to be achived by computing the barriers of only a small number out of the many subtransitions in the TN. These barriers are computed using the Conjugate Peak Refinement method. The approach is illustrated on the conformational switch of Ras p21. The best and the 12 next-best transition pathways, having rate-limiting barriers within a range of 10 kcal/mol, were identified. Two main energy ridges, which respectively involve rearrangements of the switch I and switch II loops, show that switch I must rearrange by threading Tyr32 underneath the protein backbone before the rate-limiting switch II rearrangement can occur, while the details of the switch II rearrangement differ significantly among the low-energy pathways.

The Hubbard dimer: a density functional case study of a many-body problem

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn-Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe ansatz local density approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H2.

Structural and Functional Evidence for Testosterone Activation of GPRC6A in Peripheral Tissues

G protein-coupled receptor (GPCR) family C group 6 member A (GPRC6A) is a multiligand GPCR that is activated by cations, L-amino acids, and osteocalcin. GPRC6A plays an important role in the regulation of testosterone (T) production and energy metabolism in mice. T has rapid, transcription-independent (nongenomic) effects that are mediated by a putative GPCR. We previously found that T can activate GPRC6A in vitro, but the possibility that T is a ligand for GPRC6A remains controversial. Here, we demonstrate direct T binding to GPRC6A and construct computational structural models of GPRC6A that are used to identify potential binding poses of T. Mutations of the predicted binding site residues were experimentally found to block T activation of GPRC6A, in agreement with the modeling. Using Gpr6ca(-/-) mice, we confirmed that loss of GPRC6A resulted in loss of T rapid signaling responses and elucidated several biological functions regulated by GPRC6A-dependent T rapid signaling, including T stimulation of insulin secretion in pancreatic islets and enzyme expression involved in the biosynthesis of T in Leydig cells. Finally, we identified a stereo-specific effect of an R-isomer of a selective androgen receptor modulator that is predicted to bind to and shown to activate GPRC6A but not androgen receptor. Together, our data show that GPRC6A directly mediates the rapid signaling response to T and uncovers previously unrecognized endocrine networks.

Structure and dynamics of a compact state of a multidomain protein, the mercuric ion reductase

The functional efficacy of colocalized, linked protein domains is dependent on linker flexibility and system compaction. However, the detailed characterization of these properties in aqueous solution presents an enduring challenge. Here, we employ a novel, to our knowledge, combination of complementary techniques, including small-angle neutron scattering, neutron spin-echo spectroscopy, and all-atom molecular dynamics and coarse-grained simulation, to identify and characterize in detail the structure and dynamics of a compact form of mercuric ion reductase (MerA), an enzyme central to bacterial mercury resistance. MerA possesses metallochaperone-like N-terminal domains (NmerA) tethered to its catalytic core domain by linkers. The NmerA domains are found to interact principally through electrostatic interactions with the core, leashed by the linkers so as to subdiffuse on the surface over an area close to the core C-terminal Hg(II)-binding cysteines. How this compact, dynamical arrangement may facilitate delivery of Hg(II) from NmerA to the core domain is discussed.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

Motional displacements in proteins: The origin of wave-vector-dependent values

The average mean-square displacement, 〈r(2)〉, of H atoms in a protein is frequently determined using incoherent neutron-scattering experiments. 〈r(2)〉 is obtained from the observed elastic incoherent dynamic structure factor, S(i)(Q,ω=0), assuming the form S(i)(Q,ω=0) =exp(-Q(2)〈r(2)〉/3). This is often referred to as the Gaussian approximation (GA) to S(i)(Q,ω=0). 〈r(2)〉 obtained in this way depends on the value of the wave vector, Q considered. Equivalently, the observed S(i)(Q,ω=0) deviates from the GA. We investigate the origin of the Q dependence of 〈r(2)〉 by evaluating the scattering functions in different approximations using molecular dynamics (MD) simulation of the protein lysozyme. We find that keeping only the Gaussian term in a cumulant expansion of S(Q,ω) is an accurate approximation and is not the origin of the Q dependence of 〈r(2)〉. This is demonstrated by showing that the term beyond the Gaussian is negligible and that the GA is valid for an individual atom in the protein. Rather, the Q dependence (deviation from the GA) arises from the dynamical heterogeneity of the H in the protein. Specifically it arises from representing, in the analysis of data, this diverse dynamics by a single average scattering center that has a single, average 〈r(2)〉. The observed Q dependence of 〈r(2)〉 can be used to provide information on the dynamical heterogeneity in proteins.

HackaMol: An Object-Oriented Modern Perl Library for Molecular Hacking on Multiple Scales

HackaMol is an open source, object-oriented toolkit written in Modern Perl that organizes atoms within molecules and provides chemically intuitive attributes and methods. The library consists of two components: HackaMol, the core that contains classes for storing and manipulating molecular information, and HackaMol::X, the extensions that use the core. The core is well-tested, well-documented, and easy to install across computational platforms. The goal of the extensions is to provide a more flexible space for researchers to develop and share new methods. In this application note, we provide a description of the core classes and two extensions: HackaMol::X::Calculator, an abstract calculator that uses code references to generalize interfaces with external programs, and HackaMol::X::Vina, a structured class that provides an interface with the AutoDock Vina docking program.

Why genetic modification of lignin leads to low-recalcitrance biomass

Molecular dynamics simulations show genetically modified lignin to associate less with hemicellulose than does wild type.

Examination of the perception of sweet- and bitter-like taste qualities in sucralose preferring and avoiding rats

Sucralose avoiding rats detect a bitter-like taste quality in concentrations of sucralose that are strongly preferred over water by sucralose preferring rats. Here, we investigated whether sucralose preferrers (SP) also detect a bitter-like quality in sucralose that may be masked by an increased perception of sucralose's sweet-like quality. A microstructural analysis of sucralose intake revealed that, at concentrations they avoided in preference tests, sucralose avoiders (SA) consumed smaller and fewer bouts of sucralose than SP. Interestingly, the concentration-dependent increase in sucralose preference in SP was not associated with larger bouts or increased lick rate, two measures that are expected to increase with increasing perceived sweetness. This suggests that SP can detect an aversive quality in sucralose, but this perception of a presumably bitter-like quality may be masked by increased salience of a sweet-like quality that sustains high levels of intake in SP. Further evidence for increased sweet-taste perception in SP, relative to SA, was obtained in a second study in which SP consumed more of a palatable sweet-milk diet than SA. These are the first data to suggest that SP are not blind to the bitter-like quality in sucralose, and that there may be differences in sweet-taste perception between SP and SA.

Tri-peptide reference structures for the calculation of relative solvent accessible surface area in protein amino acid residues

Relative amino acid residue solvent accessibility values allow the quantitative comparison of atomic solvent-accessible surface areas in different residue types and physical environments in proteins and in protein structural alignments. Geometry-optimised tri-peptide structures in extended solvent-exposed reference conformations have been obtained for 43 amino acid residue types at a high level of quantum chemical theory. Significant increases in side-chain solvent accessibility, offset by reductions in main-chain atom solvent exposure, were observed for standard residue types in partially geometry-optimised structures when compared to non-minimised models built from identical sets of proper dihedral angles abstracted from the literature. Optimisation of proper dihedral angles led most notably to marked increases of up to 54% in proline main-chain atom solvent accessibility compared to literature values. Similar effects were observed for fully-optimised tri-peptides in implicit solvent. The relief of internal strain energy was associated with systematic variation in N, C(α) and C(β) atom solvent accessibility across all standard residue types. The results underline the importance of optimisation of 'hard' degrees of freedom (bond lengths and valence bond angles) and improper dihedral angle values from force field or other context-independent reference values, and impact on the use of standardised fixed internal co-ordinate geometry in sampling approaches to the determination of absolute values of protein amino acid residue solvent accessibility. Quantum chemical methods provide a useful and accurate alternative to molecular mechanics methods to perform energy minimisation of peptides containing non-standard (chemically modified) amino acid residues frequently present in experimental protein structure data sets, for which force field parameters may not be available. Reference tri-peptide atomic co-ordinate sets including hydrogen atoms are made freely available.

X-ray structure of a Hg2+ complex of mercuric reductase (MerA) and quantum mechanical/molecular mechanical study of Hg2+ transfer between the C-terminal and buried catalytic site cysteine pairs

Mercuric reductase, MerA, is a key enzyme in bacterial mercury resistance. This homodimeric enzyme captures and reduces toxic Hg²⁺ to Hg⁰, which is relatively unreactive and can exit the cell passively. Prior to reduction, the Hg²⁺ is transferred from a pair of cysteines (C558′ and C559′ using Tn501 numbering) at the C-terminus of one monomer to another pair of cysteines (C136 and C141) in the catalytic site of the other monomer. Here, we present the X-ray structure of the C-terminal Hg²⁺ complex of the C136A/C141A double mutant of the Tn501 MerA catalytic core and explore the molecular mechanism of this Hg transfer with quantum mechanical/molecular mechanical (QM/MM) calculations. The transfer is found to be nearly thermoneutral and to pass through a stable tricoordinated intermediate that is marginally less stable than the two end states. For the overall process, Hg²⁺ is always paired with at least two thiolates and thus is present at both the C-terminal and catalytic binding sites as a neutral complex. Prior to Hg²⁺ transfer, C141 is negatively charged. As Hg²⁺ is transferred into the catalytic site, a proton is transferred from C136 to C559′ while C558′ becomes negatively charged, resulting in the net transfer of a negative charge over a distance of ∼7.5 Å. Thus, the transport of this soft divalent cation is made energetically feasible by pairing a competition between multiple Cys thiols and/or thiolates for Hg²⁺ with a competition between the Hg²⁺ and protons for the thiolates.

Multi-conformer ensemble docking to difficult protein targets

Large-scale ensemble docking is investigated using five proteins from the Directory of Useful Decoys (DUD, dud.docking.org ) for which docking to crystal structures has proven difficult. Molecular dynamics trajectories are produced for each protein and an ensemble of representative conformational structures extracted from the trajectories. Docking calculations are performed on these selected simulation structures and ensemble-based enrichment factors compared with those obtained using docking in crystal structures of the same protein targets or random selection of compounds. Simulation-derived snapshots are found with improved enrichment factors that increase the chemical diversity of docking hits for four of the five selected proteins. A combination of all the docking results obtained from molecular dynamics simulation followed by selection of top-ranking compounds appears to be an effective strategy for increasing the number and diversity of hits when using docking to screen large libraries of chemicals against difficult protein targets.

Eppur si muove! The 2013 Nobel Prize in Chemistry

The 2013 Nobel Prize in Chemistry has been awarded to Martin Karplus, Michael Levitt, and Arieh Warshel for their work on developing computational methods to study complex chemical systems. Their work has led to mechanistic critical insights into chemical systems both large and small and has enabled progress in a number of different fields, including structural biology.

Replica-exchange molecular dynamics simulations of cellulose solvated in water and in the ionic liquid 1-butyl-3-methylimidazolium chloride

Ionic liquids have become a popular solvent for cellulose pretreatment in biorefineries due to their efficiency in dissolution and their reusability. Understanding the interactions between cations, anions, and cellulose is key to the development of better solvents and the improvement of pretreatment conditions. While previous studies described the interactions between ionic liquids and cellulose fibers, shedding light on the initial stages of the cellulose dissolution process, we study the end state of that process by exploring the structure and dynamics of a single cellulose decamer solvated in 1-butyl-3-methyl-imidazolium chloride (BmimCl) and in water using replica-exchange molecular dynamics. In both solvents, global structural features of the cellulose chain are similar. However, analyses of local structural properties show that cellulose explores greater conformational variability in the ionic liquid than in water. For instance, in BmimCl the cellulose intramolecular hydrogen bond O3H'···O5 is disrupted more often resulting in greater flexibility of the solute. Our results indicate that the cellulose chain is more dynamic in BmimCl than in water, which may play a role in the favorable dissolution of cellulose in the ionic liquid. Calculation of the configurational entropy of the cellulose decamer confirms its higher conformational flexibility in BmimCl than in water at elevated temperatures.