MMC: A Monte Carlo laboratory.

This Note describes features of the program MMC, several of which are unique to MMC, developed over the past five decades. These include sampling in three different ensembles, biased moves, and some non-conventional analysis techniques.

Foldability and chameleon propensity of fold-switching protein sequences.

It has been shown recently by Porter & Looger that a significant number of proteins exist that can form more than one stable fold. This note examines the sequences of these fold-switching proteins by (a) calculating their foldability scores recently introduced by the present author and (b) comparing the propensity of chameleon sequences in fold switchers and in non fold switchers. It has been found that the average foldability score of the fold switchers indicates weaker foldability. As for the propensity of chameleon sequences of length 5 to 7 it was found, somewhat surprisingly, that there is only a very small difference between the fold switchers and the non fold switchers. Furthermore, when looking at amino acid propensities, for several amino acids there was even an opposing trend in the deviation of their propensities from the overall amino acid propensities.

On predicting foldability of a protein from its sequence.

Several properties of amino acid sequences corresponding to proteins that are known to fold are compared to those of randomly generated sequences and to sequences of intrinsically disordered proteins in order to find properties that distinguish folding sequences from the rest. The properties studied included helix and sheet propensities from secondary structure prediction, adjacency correlations, directionality correlations, as well as propensities of all possible triplets and quadruplets. Small differences between known folded and random sequences were observed for the adjacency and directional correlations, and significant differences were seen on the triplet and especially on the quadruplet propensities. Based on the differences in the adjacency, triplet or quadruplet propensities folding scores were defined and used to test the accuracy of foldability prediction based on these statistics. The best predictions were obtained from the quadruplet propensities.

Retro-inverso D-peptides as a novel targeted immunotherapy for Type 1 diabetes.

Over the past four decades, the number of people with Type 1 Diabetes (T1D) has increased by 4% per year, making it an important public health challenge. Currently, no curative therapy exists for T1D and the only available treatment is insulin replacement. HLA-DQ8 has been shown to present antigenic islet peptides driving the activation of CD4+ T-cells in T1D patients. Specifically, the insulin peptide InsB:9-23 activates self-reactive CD4+ T-cells, causing pancreatic beta cell destruction. The aim of the current study was to identify retro-inverso-d-amino acid based peptides (RI-D-peptides) that can suppress T-cell activation by blocking the presentation of InsB:9-23 peptide within HLA-DQ8 pocket. We identified a RI-D-peptide (RI-EXT) that inhibited InsB:9-23 binding to recombinant HLA-DQ8 molecule, as well as its binding to DQ8 expressed on human B-cells. RI-EXT prevented T-cell activation in a cellular antigen presentation assay containing human DQ8 cells loaded with InsB:9-23 peptide and murine T-cells expressing a human T-cell receptor specific for the InsB:9-23-DQ8 complex. Moreover, RI-EXT blocked T-cell activation by InsB:9-23 in a humanized DQ8 mice both ex vivo and in vivo, as shown by decreased production of IL-2 and IFN-γ and reduced lymphocyte proliferation. Interestingly, RI-EXT also blocked lymphocyte activation and proliferation by InsB:9-23 in PBMCs isolated from recent onset DQ8-T1D patients. In summary, we discovered a RI-D-peptide that blocks InsB:9-23 binding to HLA-DQ8 and its presentation to T-cells in T1D. These findings set the stage for using our approach as a novel therapy for patients with T1D and potentially other autoimmune diseases.

Inhibitor of CBP Histone Acetyltransferase Downregulates p53 Activation and Facilitates Methylation at Lysine 27 on Histone H3.

Tumor suppressor p53-directed apoptosis triggers loss of normal cells, which contributes to the side-effects from anticancer therapies. Thus, small molecules with potential to downregulate the activation of p53 could minimize pathology emerging from anticancer therapies. Acetylation of p53 by the histone acetyltransferase (HAT) domain is the hallmark of coactivator CREB-binding protein (CBP) epigenetic function. During genotoxic stress, CBP HAT-mediated acetylation is essential for the activation of p53 to transcriptionally govern target genes, which control cellular responses. Here, we present a small molecule, NiCur, which blocks CBP HAT activity and downregulates p53 activation upon genotoxic stress. Computational modeling reveals that NiCur docks into the active site of CBP HAT. On CDKN1A promoter, the recruitment of p53 as well as RNA Polymerase II and levels of acetylation on histone H3 were diminished by NiCur. Specifically, NiCur reduces the levels of acetylation at lysine 27 on histone H3, which concomitantly increases the levels of trimethylation at lysine 27. Finally, NiCur attenuates p53-directed apoptosis by inhibiting the Caspase 3 activity and cleavage of Poly (ADP-ribose) polymerase (PARP) in normal gastrointestinal epithelial cells. Collectively, NiCur demonstrates the potential to reprogram the chromatin landscape and modulate biological outcomes of CBP-mediated acetylation under normal and disease conditions.

Identifying a Small Molecule Blocking Antigen Presentation in Autoimmune Thyroiditis.

We previously showed that an HLA-DR variant containing arginine at position 74 of the DRß1 chain (DRß1-Arg74) is the specific HLA class II variant conferring risk for autoimmune thyroid diseases (AITD). We also identified 5 thyroglobulin (Tg) peptides that bound to DRß1-Arg74. We hypothesized that blocking the binding of these peptides to DRß1-Arg74 could block the continuous T-cell activation in thyroiditis needed to maintain the autoimmune response to the thyroid. The aim of the current study was to identify small molecules that can block T-cell activation by Tg peptides presented within DRß1-Arg74 pockets. We screened a large and diverse library of compounds and identified one compound, cepharanthine that was able to block peptide binding to DRß1-Arg74. We then showed that Tg.2098 is the dominant peptide when inducing experimental autoimmune thyroiditis (EAT) in NOD mice expressing human DRß1-Arg74. Furthermore, cepharanthine blocked T-cell activation by thyroglobulin peptides, in particular Tg.2098 in mice that were induced with EAT. For the first time we identified a small molecule that can block Tg peptide binding and presentation to T-cells in autoimmune thyroiditis. If confirmed cepharanthine could potentially have a role in treating human AITD.

Rescore protein-protein docked ensembles with an interface contact statistics.

The recently developed statistical measure for the type of residue-residue contact at protein complex interfaces, based on a parameter-free definition of contact, has been used to define a contact score that is correlated with the likelihood of correctness of a proposed complex structure. Comparing the proposed contact scores on the native structure and on a set of model structures the proposed measure was shown to generally favor the native structure but in itself was not able to reliably score the native structure to be the best. Adjusting the scores of redocking experiments with the contact score showed that the adjusted score was able to move up the ranking of the native-like structure among the proposed complexes when the native-like was not ranked the best by the respective program. Tests on docking of unbound proteins compared the contact scores of the complexes with the contact score of the crystal structure again showing the tendency of the contact score to favor native-like conformations. The possibility of using the contact score to improve the determination of biological dimers in a crystal structure was also explored. Proteins 2017; 85:235-241. © 2016 Wiley Periodicals, Inc.

Comparison of a quantum random number generator with pseudorandom number generators for their use in molecular Monte Carlo simulations.

Four pseudorandom number generators were compared with a physical, quantum-based random number generator using the NIST suite of statistical tests, which only the quantum-based random number generator could successfully pass. We then measured the effect of the five random number generators on various calculated properties in different Markov-chain Monte Carlo simulations. Two types of systems were tested: conformational sampling of a small molecule in aqueous solution and liquid methanol under constant temperature and pressure. The results show that poor quality pseudorandom number generators produce results that deviate significantly from those obtained with the quantum-based random number generator, particularly in the case of the small molecule in aqueous solution setup. In contrast, the widely used Mersenne Twister pseudorandom generator and a 64-bit Linear Congruential Generator with a scrambler produce results that are statistically indistinguishable from those obtained with the quantum-based random number generator. © 2017 Wiley Periodicals, Inc.

In Silico Modeling of Novel Drug Ligands for Treatment of Concussion Associated Tauopathy.

The objective of this study was to develop an in silico screening model for characterization of potential novel ligands from commercial drug libraries able to functionally activate certain olfactory receptors (ORs), which are members of the class A rhodopsin-like family of G protein couple receptors (GPCRs), in the brain of murine models of concussion. We previously found that concussions may significantly influence expression of certain ORs, for example, OR4M1 in subjects with a history of concussion/traumatic brain injury (TBI). In this study, we built a 3-D OR4M1 model and used it in in silico screening of potential novel ligands from commercial drug libraries. We report that in vitro activation of OR4M1 with the commercially available ZINC library compound 10915775 led to a significant attenuation of abnormal tau phosphorylation in embryonic cortico-hippocampal neuronal cultures derived from NSE-OR4M1 transgenic mice, possibly through modulation of the JNK signaling pathway. The attenuation of abnormal tau phosphorylation was rather selective since ZINC10915775 significantly decreased tau phosphorylation on tau Ser202/T205 (AT8 epitope) and tau Thr212/Ser214 (AT100 epitope), but not on tau Ser396/404 (PHF-1 epitope). Moreover, no response of ZINC10915775 was found in control hippocampal neuronal cultures derived from wild type littermates. Our in silico model provides novel means to pharmacologically modulate select ubiquitously expressed ORs in the brain through high affinity ligand activation to prevent and eventually to treat concussion induced down regulation of ORs and subsequent cascade of tau pathology. J. Cell. Biochem. 117: 2241-2248, 2016. © 2016 Wiley Periodicals, Inc.

The role of protein "Stability patches" in molecular recognition: A case study of the human growth hormone-receptor complex.

Dynamic characteristics of protein surfaces are among the factors determining their functional properties, including their potential participation in protein-protein interactions. The presence of clusters of static residues-"stability patches" (SPs)-is a characteristic of protein surfaces involved in intermolecular recognition. The mechanism, by with SPs facilitate molecular recognition, however, remains unclear. Analyzing the surface dynamic properties of the growth hormone and of its high-affinity variant we demonstrated that reshaping of the SPs landscape may be among the factors accountable for the improved affinity of this variant to the receptor. We hypothesized that SPs facilitate molecular recognition by moderating the conformational entropy of the unbound state, diminishing enthalpy-entropy compensation upon binding, and by augmenting the favorable entropy of desolvation. SPs mapping emerges as a valuable tool for investigating the structural basis of the stability of protein complexes and for rationalizing experimental approaches, such as affinity maturation, aimed at improving it.

Characterization of the Binding Site of Aspartame in the Human Sweet Taste Receptor.

The sweet taste receptor, a heterodimeric G protein-coupled receptor comprised of T1R2 and T1R3, binds sugars, small molecule sweeteners, and sweet proteins to multiple binding sites. The dipeptide sweetener, aspartame binds in the Venus Flytrap Module (VFTM) of T1R2. We developed homology models of the open and closed forms of human T1R2 and human T1R3 VFTMs and their dimers and then docked aspartame into the closed form of T1R2's VFTM. To test and refine the predictions of our model, we mutated various T1R2 VFTM residues, assayed activity of the mutants and identified 11 critical residues (S40, Y103, D142, S144, S165, S168, Y215, D278, E302, D307, and R383) in and proximal to the binding pocket of the sweet taste receptor that are important for ligand recognition and activity of aspartame. Furthermore, we propose that binding is dependent on 2 water molecules situated in the ligand pocket that bridge 2 carbonyl groups of aspartame to residues D142 and L279. These results shed light on the activation mechanism and how signal transmission arising from the extracellular domain of the T1R2 monomer of the sweet receptor leads to the perception of sweet taste.

Use of circular variance to quantify the deviation of a macromolecule from the spherical shape.

It is shown that the extent of deviation of a molecular shape from spherical can be characterized by comparing the distribution of the circular variances, a measure originally proposed to quantify angular spread, of the vectors from each atom to the rest of the molecule to the circular variance of a collection of atoms filling the unit sphere. Different measures for quantifying the difference between distribution are proposed and compared.

The full-length TSH receptor is stabilized by TSH ligand.

The receptor for thyroid stimulating hormone (TSHR), a GPCR, is the primary antigen in autoimmune hyperthyroidism (Graves' disease) caused by stimulating TSHR antibodies. While we have previously published a full length model of the TSHR, including its leucine rich domain (LRD), linker region (LR) and transmembrane domain (TMD), to date, only a partial LRD (aa 21-261) stabilized with TSHR autoantibodies has been crystallized. Recently, however, cryo-EM structures of the full-length TSHR have been published but they include only an incomplete LR. We have now utilized the cryo-EM models, added disulfide bonds to the LR and performed longer (3000 ns) molecular dynamic (MD) simulations to update our previous model of the entire full-length TSHR, with and without the presence of TSH ligand. As in our earlier work, the new model was embedded in a lipid membrane and was solvated with water and counterions. We found that the 3000 ns Molecular Dynamic simulations showed that the structure of the LRD and TMD were remarkably constant while the LR, known more commonly as the "hinge region", again showed significant flexibility, forming several transient secondary structural elements. Analysis of the new simulations permitted a detailed examination of the effect of TSH binding on the structure of the TSHR. We found a structure-stabilizing effect of TSH, including increased stability of the LR, which was clearly demonstrated by analyzing several intrinsic receptor properties including hydrogen bonding, fluctuation of the LRD orientation, and radius of gyration. In conclusion, we were able to quantify the flexibility of the TSHR and show its increased stability after TSH binding. These data indicated the important role of ligands in directing the signaling structure of a receptor.

Severe anemia in the Nan mutant mouse caused by sequence-selective disruption of erythroid Kruppel-like factor.

Studies of mouse models of anemia have long provided fundamental insights into red blood cell formation and function. Here we show that the semidominant mouse mutation Nan ("neonatal anemia") carries a single amino acid change (E339D) within the second zinc finger of the erythroid Krüppel-like factor (EKLF), a critical erythroid regulatory transcription factor. The mutation alters the DNA-binding specificity of EKLF so that it no longer binds promoters of a subset of its DNA targets. Remarkably, even when mutant Nan and wild-type EKLF alleles are expressed at equivalent levels, the mutant form selectively interferes with expression of EKLF target genes whose promoter elements it no longer binds. This interference yields a distorted genetic output and selective protein deficiencies that differ from those seen in EKLF-heterozygous and EKLF-null red blood cells and presents a unique and unexpected mechanism of inherited disease.

Functional Water Channels Within the TSH Receptor: A New Paradigm for TSH Action With Disease Implications.

The thyroid-stimulating hormone receptor (TSHR) transmembrane domain (TMD) is found in the plasma membrane and consists of lipids and water molecules. To understand the role of TSHR-associated water molecules, we used molecular dynamic simulations of the TMD and identified a network of putative receptor-associated transmembrane water channels. This result was confirmed with extended simulations of the full-length TSHR with and without TSH ligand binding. While the transport time observed in the simulations via the TSHR protein was slower than via the lipid bilayer itself, we found that significantly more water traversed via the TSHR than via the lipid bilayer, which more than doubled with the binding of TSH. Using rat thyroid cells (FRTL-5) and a calcein fluorescence technique, we measured cell volumes after blockade of aquaporins 1 and 4, the major thyroid cell water transporters. TSH showed a dose-dependent ability to influence water transport, and similar effects were observed with stimulating TSHR autoantibodies. Small molecule TSHR agonists, which are allosteric activators of the TMD, also enhanced water transport, illustrating the role of the TMD in this phenomenon. Furthermore, the water channel pathway was also mapped across 2 activation motifs within the TSHR TMD, suggesting how water movement may influence activation of the receptor. In pathophysiological conditions such as hypothyroidism and hyperthyroidism where TSH concentrations are highly variable, this action of TSH may greatly influence water movement in thyroid cells and many other extrathyroidal sites where the TSHR is expressed, thus affecting normal cellular function.

Substrate binding and inhibition of the anion exchanger 1 transporter.

Anion exchanger 1 (AE1), a member of the solute carrier (SLC) family, is the primary bicarbonate transporter in erythrocytes, regulating pH levels and CO2 transport between lungs and tissues. Previous studies characterized its role in erythrocyte structure and provided insight into transport regulation. However, key questions remain regarding substrate binding and transport, mechanisms of drug inhibition and modulation by membrane components. Here we present seven cryo-EM structures in apo, bicarbonate-bound and inhibitor-bound states. These, combined with uptake and computational studies, reveal important molecular features of substrate recognition and transport, and illuminate sterol binding sites, to elucidate distinct inhibitory mechanisms of research chemicals and prescription drugs. We further probe the substrate binding site via structure-based ligand screening, identifying an AE1 inhibitor. Together, our findings provide insight into mechanisms of solute carrier transport and inhibition.

Predicting protein interactions by Brownian dynamics simulations.

We present a newly adapted Brownian-Dynamics (BD)-based protein docking method for predicting native protein complexes. The approach includes global BD conformational sampling, compact complex selection, and local energy minimization. In order to reduce the computational costs for energy evaluations, a shell-based grid force field was developed to represent the receptor protein and solvation effects. The performance of this BD protein docking approach has been evaluated on a test set of 24 crystal protein complexes. Reproduction of experimental structures in the test set indicates the adequate conformational sampling and accurate scoring of this BD protein docking approach. Furthermore, we have developed an approach to account for the flexibility of proteins, which has been successfully applied to reproduce the experimental complex structure from the structure of two unbounded proteins. These results indicate that this adapted BD protein docking approach can be useful for the prediction of protein-protein interactions.

Tools for Characterizing Proteins: Circular Variance, Mutual Proximity, Chameleon Sequences, and Subsequence Propensities.

For the characterization of various aspects of protein structures, four useful concepts are discussed: chameleon sequences, circular variance, mutual proximity, and a subsequence-based foldability score. These concepts were used in estimating foldability of globular, intrinsically disordered and fold-switching proteins, properties of protein-protein interfaces, quantifying sphericity, helping to improve protein-protein docking scores, and estimating the effect of mutations on stability. A conjecture about the Achilles' heel of proteins is presented as well.

Modeling TSH Receptor Dimerization at the Transmembrane Domain.

Biophysical studies have established that the thyrotropin (TSH) receptor (TSHR) undergoes posttranslational modifications including dimerization. Following our earlier simulation of a TSHR-transmembrane domain (TMD) monomer (called TSHR-TMD-TRIO) we have now proceeded with a molecular dynamics simulation (MD) of TSHR-TMD dimerization using this improved membrane-embedded model. The starting structure was the TMD protein with all extracellular and intracellular loops and internal waters, which was placed in the relative orientation of the model originally generated with Brownian dynamics. Furthermore, this model was embedded in a DPPC lipid bilayer further solvated with water and added salt. Data from the MD simulation studies showed that the dimeric subunits stayed in the same relative orientation and distance during the 1000 ns of study. Comparison of representative conformations of the individual monomers when dimerized with the conformations from the monomer simulation showed subtle differences as represented by the backbone root mean square deviations. Differences in the conformations of the ligand-binding sites, suggesting variable affinities for these "hot spots," were also revealed by comparing the docking scores of 46 small-molecule ligands that included known TSHR agonists and antagonists as well as their derivatives. These data add further insight into the tendency of the TSHR-TMD to form dimeric and oligomeric structures and show that the differing conformations influence small-molecule binding sites within the TMD.

Computational model of the full-length TSH receptor.

(GPCR)The receptor for TSH receptor (TSHR), a G protein coupled receptor (GPCR), is of particular interest as the primary antigen in autoimmune hyperthyroidism (Graves' disease) caused by stimulating TSHR antibodies. To date, only one domain of the extracellular region of the TSHR has been crystallized. We have run a 1000 ns molecular dynamic simulation on a model of the entire TSHR generated by merging the extracellular region of the receptor, obtained using artificial intelligence, with our recent homology model of the transmembrane domain, embedded it in a lipid membrane and solvated it with water and counterions. The simulations showed that the structure of the transmembrane and leucine-rich domains were remarkably constant while the linker region (LR), known more commonly as the 'hinge region,' showed significant flexibility, forming several transient secondary structural elements. Furthermore, the relative orientation of the leucine-rich domain with the rest of the receptor was also seen to be variable. These data suggest that this LR is an intrinsically disordered protein. Furthermore, preliminary data simulating the full TSHR model complexed with its ligand (TSH) showed that (a) there is a strong affinity between the LR and TSH ligand and (b) the association of the LR and the TSH ligand reduces the structural fluctuations in the LR. This full-length model illustrates the importance of the LR in responding to ligand binding and lays the foundation for studies of pathologic TSHR autoantibodies complexed with the TSHR to give further insight into their interaction with the flexible LR.