GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

The dynamics of interleukin-8 and its interaction with human CXC receptor I peptide

Interleukin-8 (CXCL8, IL-8) is a proinflammatory chemokine important for the regulation of inflammatory and immune responses via its interaction with G-protein coupled receptors, including CXC receptor 1 (CXCR1). CXCL8 exists as both a monomer and as a dimer at physiological concentrations, yet the molecular basis of CXCL8 interaction with its receptor as well as the importance of CXCL8 dimer formation remain poorly characterized. Although several biological studies have indicated that both the CXCL8 monomer and dimer are active, biophysical studies have reported conflicting results regarding the binding of CXCL8 to CXCR1. To clarify this problem, we expressed and purified a peptide (hCXCR1pep) corresponding to the N-terminal region of human CXCR1 (hCXCR1) and utilized nuclear magnetic resonance (NMR) spectroscopy to interrogate the binding of wild-type CXCL8 and a previously reported mutant (CXCL8M) that stabilizes the monomeric form. Our data reveal that the CXCL8 monomer engages hCXCR1pep with a slightly higher affinity than the CXCL8 dimer, but that the CXCL8 dimer does not dissociate upon binding hCXCR1pep. These investigations also showed that CXCL8 is dynamic on multiple timescales, which may help explain the versatility in this interleukin for engaging its target receptors.

From a structural average to the conformational ensemble of a DNA bulge

Direct experimental measurements of conformational ensembles are critical for understanding macromolecular function, but traditional biophysical methods do not directly report the solution ensemble of a macromolecule. Small-angle X-ray scattering interferometry has the potential to overcome this limitation by providing the instantaneous distance distribution between pairs of gold-nanocrystal probes conjugated to a macromolecule in solution. Our X-ray interferometry experiments reveal an increasing bend angle of DNA duplexes with bulges of one, three, and five adenosine residues, consistent with previous FRET measurements, and further reveal an increasingly broad conformational ensemble with increasing bulge length. The distance distributions for the AAA bulge duplex (3A-DNA) with six different Au-Au pairs provide strong evidence against a simple elastic model in which fluctuations occur about a single conformational state. Instead, the measured distance distributions suggest a 3A-DNA ensemble with multiple conformational states predominantly across a region of conformational space with bend angles between 24 and 85 degrees and characteristic bend directions and helical twists and displacements. Additional X-ray interferometry experiments revealed perturbations to the ensemble from changes in ionic conditions and the bulge sequence, effects that can be understood in terms of electrostatic and stacking contributions to the ensemble and that demonstrate the sensitivity of X-ray interferometry. Combining X-ray interferometry ensemble data with molecular dynamics simulations gave atomic-level models of representative conformational states and of the molecular interactions that may shape the ensemble, and fluorescence measurements with 2-aminopurine-substituted 3A-DNA provided initial tests of these atomistic models. More generally, X-ray interferometry will provide powerful benchmarks for testing and developing computational methods.

Evolutionarily conserved and conformationally constrained short peptides might serve as DNA recognition elements in intrinsically disordered regions

Despite recent advances, it is yet not clear how intrinsically disordered regions in proteins recognize their targets without any defined structures. Short linear motifs had been proposed to mediate molecular recognition by disordered regions; however, the underlying structural prerequisite remains elusive. Moreover, the role of short linear motifs in DNA recognition has not been studied. We report a repertoire of short evolutionarily Conserved Recognition Elements (CoREs) in long intrinsically disordered regions, which have very distinct amino-acid propensities from those of known motifs, and exhibit a strong tendency to retain their three-dimensional conformations compared to adjacent regions. The majority of CoREs directly interact with the DNA in the available 3D structures, which is further supported by literature evidence, analyses of ΔΔG values of DNA-binding energies and threading-based prediction of DNA binding potential. CoREs were enriched in cancer-associated missense mutations, further strengthening their functional nature. Significant enrichment of glycines in CoREs and the preference of glycyl ϕ-Ψ values within the left-handed bridge range in the l-disallowed region of the Ramachandran plot suggest that Gly-to-nonGly mutations within CoREs might alter the backbone conformation and consequently the function, a hypothesis that we reconciled using available mutation data. We conclude that CoREs might serve as bait for DNA recognition by long disordered regions and that certain mutations in these peptides can disrupt their DNA binding potential and consequently the protein function. We further hypothesize that the preferred conformations of CoREs and of glycyl residues therein might play an important role in DNA binding. The highly ordered nature of CoREs hints at a therapeutic strategy to inhibit malicious molecular interactions using small molecules mimicking CoRE conformations.

How to tackle the issues in free energy simulations of long amphiphiles interacting with lipid membranes: convergence and local membrane deformations

One of the great challenges in membrane biophysics is to find a means to foster the transport of drugs across complex membrane structures. In this spirit, we elucidate methodological challenges associated with free energy computations of complex chainlike molecules across lipid membranes. As an appropriate standard molecule to this end, we consider 7-nitrobenz-2-oxa-1,3-diazol-4-yl-labeled fatty amine, NBD-Cn, which is here dealt with as a homologous series with varying chain lengths. We found the membrane-water interface region to be highly sensitive to details in free energy computations. Despite considerable simulation times, we observed substantial hysteresis, the cause being the small frequency of insertion/desorption events of the amphiphile's alkyl chain in the membrane interface. The hysteresis was most pronounced when the amphiphile was pulled from water to the membrane and compromised the data that were not in line with experiments. The subtleties in umbrella sampling for computing distance along the transition path were also observed to be potential causes of artifacts. With the PGD (pull geometry distance) scheme, in which the distance from the molecule was computed to a reference plane determined by an average over all lipids in the membrane, we found marked deformations in membrane structure when the amphiphile was close to the membrane. The deformations were weaker with the PGC (pull geometry cylinder) method, where the reference plane is chosen based on lipids that are within a cylinder of radius 1.7 nm from the amphiphile. Importantly, the free energy results given by PGC were found to be qualitatively consistent with experimental data, while the PGD results were not. We conclude that with long amphiphiles there is reason for concern with regard to computations of their free energy profiles. The membrane-water interface is the region where the greatest care is warranted.

The crystal structure of Giardia duodenalis 14-3-3 in the apo form: when protein post-translational modifications make the difference

The 14-3-3s are a family of dimeric evolutionary conserved pSer/pThr binding proteins that play a key role in multiple biological processes by interacting with a plethora of client proteins. Giardia duodenalis is a flagellated protozoan that affects millions of people worldwide causing an acute and chronic diarrheal disease. The single giardial 14-3-3 isoform (g14-3-3), unique in the 14-3-3 family, needs the constitutive phosphorylation of Thr214 and the polyglycylation of its C-terminus to be fully functional in vivo. Alteration of the phosphorylation and polyglycylation status affects the parasite differentiation into the cyst stage. To further investigate the role of these post-translational modifications, the crystal structure of the g14-3-3 was solved in the unmodified apo form. Oligomers of g14-3-3 were observed due to domain swapping events at the protein C-terminus. The formation of filaments was supported by TEM. Mutational analysis, in combination with native PAGE and chemical cross-linking, proved that polyglycylation prevents oligomerization. In silico phosphorylation and molecular dynamics simulations supported a structural role for the phosphorylation of Thr214 in promoting target binding. Our findings highlight unique structural features of g14-3-3 opening novel perspectives on the evolutionary history of this protein family and envisaging the possibility to develop anti-giardial drugs targeting g14-3-3.

Stabilization of the GluCl ligand-gated ion channel in the presence and absence of ivermectin

Improving our understanding of the mechanisms and effects of anesthetics is a critically important part of neuroscience. The currently dominant theory is that anesthetics and similar molecules act by binding to Cys-loop receptors in the postsynaptic terminal of nerve cells and potentiate or inhibit their function. Although structures for some of the most important mammalian channels have still not been determined, a number of important results have been derived from work on homologous cationic channels in bacteria. However, partly due to the lack of a nervous system in bacteria, there are a number of questions about how these results relate to higher organisms. The recent determination of a structure of the eukaryotic chloride channel, GluCl, is an important step toward accurate modeling of mammalian channels, because it is more similar in function to human Cys-loop receptors such as GABAAR or GlyR. One potential issue with using GluCl to model other receptors is the presence of the large ligand ivermectin (IVM) positioned between all five subunits. Here, we have performed a series of microsecond molecular simulations to study how the dynamics and structure of GluCl change in the presence versus absence of IVM. When the ligand is removed, subunits move at least 2 Å closer to each other compared to simulations with IVM bound. In addition, the pore radius shrinks to 1.2 Å, all of which appears to support a model where IVM binding between subunits stabilizes an open state, and that the relaxed nonIVM conformations might be suitable for modeling other channels. Interestingly, the presence of IVM also has an effect on the structure of the important loop C located at the neurotransmitter-binding pocket, which might help shed light on its partial agonist behavior.

ST-analyzer: a web-based user interface for simulation trajectory analysis

Molecular dynamics (MD) simulation has become one of the key tools to obtain deeper insights into biological systems using various levels of descriptions such as all-atom, united-atom, and coarse-grained models. Recent advances in computing resources and MD programs have significantly accelerated the simulation time and thus increased the amount of trajectory data. Although many laboratories routinely perform MD simulations, analyzing MD trajectories is still time consuming and often a difficult task. ST-analyzer, http://im.bioinformatics.ku.edu/st-analyzer, is a standalone graphical user interface (GUI) toolset to perform various trajectory analyses. ST-analyzer has several outstanding features compared to other existing analysis tools: (i) handling various formats of trajectory files from MD programs, such as CHARMM, NAMD, GROMACS, and Amber, (ii) intuitive web-based GUI environment--minimizing administrative load and reducing burdens on the user from adapting new software environments, (iii) platform independent design--working with any existing operating system, (iv) easy integration into job queuing systems--providing options of batch processing either on the cluster or in an interactive mode, and (v) providing independence between foreground GUI and background modules--making it easier to add personal modules or to recycle/integrate pre-existing scripts utilizing other analysis tools. The current ST-analyzer contains nine main analysis modules that together contain 18 options, including density profile, lipid deuterium order parameters, surface area per lipid, and membrane hydrophobic thickness. This article introduces ST-analyzer with its design, implementation, and features, and also illustrates practical analysis of lipid bilayer simulations.

Structural basis of the proinflammatory signaling complex mediated by TSLP

Thymic stromal lymphopoietin (TSLP), a cytokine produced by epithelial cells at barrier surfaces, is pivotal for the development of widespread chronic inflammatory disorders such as asthma and atopic dermatitis. The structure of the mouse TSLP-mediated signaling complex reveals how TSLP establishes extensive interfaces with its cognate receptor (TSLPR) and the shared interleukin 7 receptor α-chain (IL-7Rα) to evoke membrane-proximal receptor-receptor contacts poised for intracellular signaling. Binding of TSLP to TSLPR is a mechanistic prerequisite for recruitment of IL-7Rα to the high-affinity ternary complex, which we propose is coupled to a structural switch in TSLP at the crossroads of the cytokine-receptor interfaces. Functional interrogation of TSLP-receptor interfaces points to putative interaction hotspots that could be exploited for antagonist design. Finally, we derive the structural rationale for the functional duality of IL-7Rα and establish a consensus for the geometry of ternary complexes mediated by interleukin 2 (IL-2)-family cytokines.

Specific non-local interactions are not necessary for recovering native protein dynamics

The elastic network model (ENM) is a widely used method to study native protein dynamics by normal mode analysis (NMA). In ENM we need information about all pairwise distances, and the distance between contacting atoms is restrained to the native value. Therefore ENM requires O(N²) information to realize its dynamics for a protein consisting of N amino acid residues. To see if (or to what extent) such a large amount of specific structural information is required to realize native protein dynamics, here we introduce a novel model based on only O(N) restraints. This model, named the ‘contact number diffusion’ model (CND), includes specific distance restraints for only local (along the amino acid sequence) atom pairs, and semi-specific non-local restraints imposed on each atom, rather than atom pairs. The semi-specific non-local restraints are defined in terms of the non-local contact numbers of atoms. The CND model exhibits the dynamic characteristics comparable to ENM and more correlated with the explicit-solvent molecular dynamics simulation than ENM. Moreover, unrealistic surface fluctuations often observed in ENM were suppressed in CND. On the other hand, in some ligand-bound structures CND showed larger fluctuations of buried protein atoms interacting with the ligand compared to ENM. In addition, fluctuations from CND and ENM show comparable correlations with the experimental B-factor. Although there are some indications of the importance of some specific non-local interactions, the semi-specific non-local interactions are mostly sufficient for reproducing the native protein dynamics.

Scalable Evaluation of Polarization Energy and Associated Forces in Polarizable Molecular Dynamics: II.Towards Massively Parallel Computations using Smooth Particle Mesh Ewald

In this paper, we present a scalable and efficient implementation of point dipole-based polarizable force fields for molecular dynamics (MD) simulations with periodic boundary conditions (PBC). The Smooth Particle-Mesh Ewald technique is combined with two optimal iterative strategies, namely, a preconditioned conjugate gradient solver and a Jacobi solver in conjunction with the Direct Inversion in the Iterative Subspace for convergence acceleration, to solve the polarization equations. We show that both solvers exhibit very good parallel performances and overall very competitive timings in an energy-force computation needed to perform a MD step. Various tests on large systems are provided in the context of the polarizable AMOEBA force field as implemented in the newly developed Tinker-HP package which is the first implementation for a polarizable model making large scale experiments for massively parallel PBC point dipole models possible. We show that using a large number of cores offers a significant acceleration of the overall process involving the iterative methods within the context of spme and a noticeable improvement of the memory management giving access to very large systems (hundreds of thousands of atoms) as the algorithm naturally distributes the data on different cores. Coupled with advanced MD techniques, gains ranging from 2 to 3 orders of magnitude in time are now possible compared to non-optimized, sequential implementations giving new directions for polarizable molecular dynamics in periodic boundary conditions using massively parallel implementations.

The effect of structural parameters and positive charge distance on the interaction free energy of antimicrobial peptides with membrane surface

Many attempts have been made to find hints explaining the relationship between physicochemical and structural properties of antimicrobial peptides (AMPs) which are relevant to their antimicrobial activities. We here found that there is a difference in the percentages of hydrophobic, hydrophilic, and charged residues between AMPs killing both bacteria and fungi (Group A) and AMPs that only kill bacteria (Group B). The percentage of charged residues in Group A AMPs is highly elevated, while in Group B the percentage of hydrophobic residues is increased. This result suggests a sequence-based mechanism of selectivity for AMPs. Moreover, we examined how the distance between basic residues affects the interaction free energy of AMPs with the membrane surface, since most of the known AMPs act by membrane perturbation. We measured the average distance between basic residues throughout the 3D structure of AMPs by defining Dpr parameter and calculated the interaction free energy for 10 AMPs that interacted with the DPPC membrane using molecular dynamics simulation. We found that the changes of the interaction free energy correlates with the change of Dpr by a linear regression coefficient of r(2 )= .47 and a cubic regression coefficient of r(2 )= .70.

CHARMM-GUI PACE CG Builder for solution, micelle, and bilayer coarse-grained simulations

Coarse-grained (CG) and multiscale simulations are widely used to study large biological systems. However, preparing the simulation system is time-consuming when the system has multiple components, because each component must be arranged carefully as in protein/micelle or protein/bilayer systems. We have developed CHARMM-GUI PACE CG Builder for building solution, micelle, and bilayer systems using the PACE force field, a united-atom (UA) model for proteins, and the Martini CG force field for water, ions, and lipids. The robustness of PACE CG Builder is validated by simulations of various systems in solution (α3D, fibronectin, and lysozyme), micelles (Pf1, DAP12-NKG2C, OmpA, and DHPC-only micelle), and bilayers (GpA, OmpA, VDAC, MscL, OmpF, and lipid-only bilayers for six lipids). The micelle’s radius of gyration, the bilayer thickness, and the per-lipid area in bilayers are comparable to the values from previous all-atom and CG simulations. Most tested proteins have root-mean squared deviations of less than 3 Å. We expect PACE CG Builder to be a useful tool for modeling/refining large, complex biological systems at the mixed UA/CG level.

Extensive conformational heterogeneity within protein cores

Basic principles of statistical mechanics require that proteins sample an ensemble of conformations at any nonzero temperature. However, it is still common to treat the crystallographic structure of a protein as the structure of its native state, largely because high-resolution structural characterization of protein flexibility remains a profound challenge. To assess the typical degree of conformational heterogeneity within folded proteins, we construct Markov state models describing the thermodynamics and kinetics of proteins ranging from 72 to 263 residues in length. Each of these models is built from hundreds of microseconds of atomically detailed molecular dynamics simulations. Examination of the side-chain degrees of freedom reveals that almost every residue visits at least two rotameric states over this time frame, with rotamer transition rates spanning a wide range of time scales (from nanoseconds to tens of microseconds). We also report substantial backbone dynamics on time scales longer than are typically addressed by experimental measures of protein flexibility, such as NMR order parameters. Finally, we demonstrate that these extensive rearrangements are consistent with NMR and crystallographic data, which supports the validity of our models. Altogether, these results depict the interior of proteins not as well-ordered solids, as is often imagined, but instead as dense fluids, which undergo substantial structural fluctuations despite their high packing fraction.

Structural interrogation of phosphoproteome identified by mass spectrometry reveals allowed and disallowed regions of phosphoconformation

Background

High-throughput mass spectrometric (HT-MS) study is the method of choice for monitoring global changes in proteome. Data derived from these studies are meant for further validation and experimentation to discover novel biological insights. Here we evaluate use of relative solvent accessible surface area (rSASA) and DEPTH as indices to assess experimentally determined phosphorylation events deposited in PhosphoSitePlus.

Results

Based on accessibility, we map these identifications on allowed (accessible) or disallowed (inaccessible) regions of phosphoconformation. Surprisingly a striking number of HT-MS/MS derived events (1461/5947 sites or 24.6%) are present in the disallowed region of conformation. By considering protein dynamics, autophosphorylation events and/or the sequence specificity of kinases, 13.8% of these phosphosites can be moved to the allowed region of conformation. We also demonstrate that rSASA values can be used to increase the confidence of identification of phosphorylation sites within an ambiguous MS dataset.

Conclusion

While MS is a stand-alone technique for the identification of vast majority of phosphorylation events, identifications within disallowed region of conformation will benefit from techniques that independently probe for phosphorylation and protein dynamics. Our studies also imply that trapping alternate protein conformations may be a viable alternative to the design of inhibitors against mutation prone drug resistance kinases.

Modified Amber Force Field Correctly Models the Conformational Preference for Tandem GA pairs in RNA

Molecular mechanics with all-atom models was used to understand the conformational preference of tandem guanine-adenine (GA) noncanonical pairs in RNA. These tandem GA pairs play important roles in determining stability, flexibility, and structural dynamics of RNA tertiary structures. Previous solution structures showed that these tandem GA pairs adopt either imino (cis Watson-Crick/Watson-Crick A-G) or sheared (trans Hoogsteen/sugar edge A-G) conformations depending on the sequence and orientation of the adjacent closing base pairs. The solution structures (GCGGACGC)2 [Biochemistry, 1996, 35, 9677-9689] and (GCGGAUGC)2 [Biochemistry, 2007, 46, 1511-1522] demonstrate imino and sheared conformations for the two central GA pairs, respectively. These systems were studied using molecular dynamics and free energy change calculations for conformational changes, using umbrella sampling. For the structures to maintain their native conformations during molecular dynamics simulations, a modification to the standard Amber ff10 force field was required, which allowed the amino group of guanine to leave the plane of the base [J. Chem. Theory Comput., 2009, 5, 2088-2100] and form out-of-plane hydrogen bonds with a cross-strand cytosine or uracil. The requirement for this modification suggests the importance of out-of-plane hydrogen bonds in stabilizing the native structures. Free energy change calculations for each sequence demonstrated the correct conformational preference when the force field modification was used, but the extent of the preference is underestimated.

Converging free energies of binding in cucurbit[7]uril and octa-acid host-guest systems from SAMPL4 using expanded ensemble simulations

Molecular containers such as cucurbit[7]uril (CB7) and the octa-acid (OA) host are ideal simplified model test systems for optimizing and analyzing methods for computing free energies of binding intended for use with biologically relevant protein-ligand complexes. To this end, we have performed initially blind free energy calculations to determine the free energies of binding for ligands of both the CB7 and OA hosts. A subset of the selected guest molecules were those included in the SAMPL4 prediction challenge. Using expanded ensemble simulations in the dimension of coupling host-guest intermolecular interactions, we are able to show that our estimates in most cases can be demonstrated to fully converge and that the errors in our estimates are due almost entirely to the assigned force field parameters and the choice of environmental conditions used to model experiment. We confirm the convergence through the use of alternative simulation methodologies and thermodynamic pathways, analyzing sampled conformations, and directly observing changes of the free energy with respect to simulation time. Our results demonstrate the benefits of enhanced sampling of multiple local free energy minima made possible by the use of expanded ensemble molecular dynamics and may indicate the presence of significant problems with current transferable force fields for organic molecules when used for calculating binding affinities, especially in non-protein chemistries.

A new insight into mushroom tyrosinase inhibitors: docking, pharmacophore-based virtual screening, and molecular modeling studies

Tyrosinase, a widely spread enzyme in micro-organisms, animals, and plants, participates in two rate-limiting steps in melanin formation pathway which is responsible for skin protection against UV lights' harm whose functional deficiency result in serious dermatological diseases. This enzyme seems to be responsible for neuromelanin formation in human brain as well. In plants, the enzyme leads the browning pathway which is commonly observed in injured tissues that is economically very unfavorable. Among different types of tyrosinase, mushroom tyrosinase has the highest homology with the mammalian tyrosinase and the only commercial tyrosinase available. In this study, ligand-based pharmacophore drug discovery method was applied to rapidly identify mushroom tyrosinase enzyme inhibitors using virtual screening. The model pharmacophore of essential interactions was developed and refined studying already experimentally discovered potent inhibitors employing Docking analysis methodology. After pharmacophore virtual screening and binding modes prediction, 14 compounds from ZINC database were identified as potent inhibitors of mushroom tyrosinase which were classified into five groups according to their chemical structures. The inhibition behavior of the discovered compounds was further studied through Classical Molecular Dynamic Simulations and the conformational changes induced by the presence of the studied ligands were discussed and compared to those of the substrate, tyrosine. According to the obtained results, five novel leads are introduced to be further optimized or directly used as potent inhibitors of mushroom tyrosinase.

The bitter barricading of prostaglandin biosynthesis pathway: understanding the molecular mechanism of selective cyclooxygenase-2 inhibition by amarogentin, a secoiridoid glycoside from Swertia chirayita

Swertia chirayita, a medicinal herb inhabiting the challenging terrains and high altitudes of the Himalayas, is a rich source of essential phytochemical isolates. Amarogentin, a bitter secoiridoid glycoside from S. chirayita, shows varied activity in several patho-physiological conditions, predominantly in leishmaniasis and carcinogenesis. Experimental analysis has revealed that amarogentin downregulates the cyclooxygenase-2 (COX-2) activity and helps to curtail skin carcinogenesis in mouse models; however, there exists no account on selective inhibition of the inducible cyclooxygenase (COX) isoform by amarogentin. Hence the computer-aided drug discovery methods were used to unravel the COX-2 inhibitory mechanism of amarogentin and to check its selectivity for the inducible isoform over the constitutive one. The generated theoretical models of both isoforms were subjected to molecular docking analysis with amarogentin and twenty-one other Food and Drug Authority (FDA) approved lead molecules. The post-docking binding energy profile of amarogentin was comparable to the binding energy profiles of the FDA approved selective COX-2 inhibitors. Subsequent molecular dynamics simulation analysis delineated the difference in the stability of both complexes, with amarogentin-COX-2 complex being more stable after 40ns simulation. The total binding free energy calculated by MMGBSA for the amarogentin-COX-2 complex was −52.35 KCal/mol against a binding free energy of −8.57 KCal/mol for amarogentin-COX-1 complex, suggesting a possible selective inhibition of the COX-2 protein by the natural inhibitor. Amarogentin achieves this potential selectivity by small, yet significant, structural differences inherent to the binding cavities of the two isoforms. Hypothetically, it might block the entry of the natural substrates in the hydrophobic binding channel of the COX-2, inhibiting the cyclooxygenation step. To sum up briefly, this work highlights the mechanism of the possible selective COX-2 inhibition by amarogentin and endorses the possibility of obtaining efficient, futuristic and targeted therapeutic agents for relieving inflammation and malignancy from this phytochemical source.

Role of berberine in anti-bacterial as a high-affinity LPS antagonist binding to TLR4/MD-2 receptor

BACKGROUND

Berberine is an isoquinoline alkaloid mainly extracted from Rhizoma Coptidis and has been shown to possess a potent inhibitory activity against bacterial. However, the role of berberine in anti-bacterial action has not been extensively studied.

METHODS

The animal model was established to investigate the effects of berberine on bacterial and LPS infection. Docking analysis, Molecular dynamics simulations and Real-time RT-PCR analysis was adopted to investigate the molecular mechanism.

RESULTS

Treatment with 40 mg/kg berberine significantly increased the survival rate of mice challenged with Salmonella typhimurium (LT2), but berberine show no effects in bacteriostasis. Further study indicated that treatment with 0.20 g/kg berberine markedly increased the survival rate of mice challenged with 2 EU/ml bacterial endotoxin (LPS) and postpone the death time of the dead mice. Moreover, pretreatment with 0.05 g/kg berberine significantly lower the increasing temperature of rabbits challenged with LPS. The studies of molecular mechanism demonstrated that Berberine was able to bind to the TLR4/MD-2 receptor, and presented higher affinity in comparison with LPS. Furthermore, berberine could significantly suppressed the increasing expression of NF-κB, IL-6, TNFα, and IFNβ in the RAW264.7 challenged with LPS.

CONCLUSION

Berberine can act as a LPS antagonist and block the LPS/TLR4 signaling from the sourse, resulting in the anti-bacterial action.

Sequence and 3D structure based analysis of TNT degrading proteins in Arabidopsis thaliana

TNT, accidentally released at several manufacturing sites, contaminates ground water and soil. It has a toxic effect to algae and invertebrate, and chronic exposure to TNT also causes harmful effects to human. On the other hand, many plants including Arabidopsis thaliana have the ability to metabolize TNT either completely or at least to a reduced less toxic form. In A. thaliana, the enzyme UDP glucosyltransferase (UDPGT) can further conjugate the reduced forms 2-HADNT and 4-HADNT (2-hydroxylamino-4, 6- dinitrotoluene and 4-hydroxylamino-2, 6- dinitrotoluene) of TNT. Based on the experimental analysis, existing literature and phylogenetic analysis, it is evident that among 107 UDPGT proteins only six are involved in the TNT degrading process. A total of 13 UDPGT proteins including five of these TNT degrading proteins fall within the same group of phylogeny. Thus, these 13 UDPGT proteins have been classified into two groups, TNT-degrading and TNT-non-degrading proteins. To understand the differences in TNT-degrading capacities; using homology modeling we first predicted two structures, taking one representative sequence from both the groups. Next, we performed molecular docking of the modeled structure and TNT reduced form 2-hydroxylamino-4, 6- dinitrotoluene (2-HADNT). We observed that while the Trp residue located within the active site region of the TNT- degrading protein showed π-Cation interaction; such type of interaction was absent in TNT-non-degrading protein, as the respective Trp residue lay outside of the pocket in this case. We observed the conservation of this π-Cation interaction during MD simulation of TNT-degrading protein. Thus, the position and the orientation of the active site residue Trp could explain the presence and absence of TNT-degrading capacity of the UDPGT proteins.

Role of κ→λ light-chain constant-domain switch in the structure and functionality of A17 reactibody

The engineering of catalytic function in antibodies requires precise information on their structure. Here, results are presented that show how the antibody domain structure affects its functionality. The previously designed organophosphate-metabolizing reactibody A17 has been re-engineered by replacing its constant κ light chain by the λ chain (A17λ), and the X-ray structure of A17λ has been determined at 1.95 Å resolution. It was found that compared with A17κ the active centre of A17λ is displaced, stabilized and made more rigid owing to interdomain interactions involving the CDR loops from the VL and VH domains. These VL/VH domains also have lower mobility, as deduced from the atomic displacement parameters of the crystal structure. The antibody elbow angle is decreased to 126° compared with 138° in A17κ. These structural differences account for the subtle changes in catalytic efficiency and thermodynamic parameters determined with two organophosphate ligands, as well as in the affinity for peptide substrates selected from a combinatorial cyclic peptide library, between the A17κ and A17λ variants. The data presented will be of interest and relevance to researchers dealing with the design of antibodies with tailor-made functions.

Crystallographic analysis of Neisseria meningitidis PorB extracellular loops potentially implicated in TLR2 recognition

Among all Neisseriae species, Neisseria meningitidis and Neisseria gonorrhoeae are the only human pathogens, causative agents of bacterial meningitis and gonorrhoea, respectively. PorB, a pan-Neisseriae trimeric porin that mediates diffusive transport of essential molecules across the bacterial outer membrane, is also known to activate host innate immunity via Toll-like receptor 2 (TLR2)-mediated signaling. The molecular mechanism of PorB binding to TLR2 is not known, but it has been hypothesized that electrostatic interactions contribute to ligand/receptor binding. Strain-specific sequence variability in the surface-exposed loops of PorB which are potentially implicated in TLR2 binding, may explain the difference in TLR2-mediated cell activation in vitro by PorB homologs from the commensal Neisseriae lactamica and the pathogen N. meningitidis. Here, we report a comparative structural analysis of PorB from N. meningitidis serogroup B strain 8765 (63% sequence homology with PorB from N. meningitidis serogroup W135) and a mutant in which amino acid substitutions in the extracellular loop 7 lead to significantly reduced TLR2-dependent activity in vitro. We observe that this mutation both alters the loop conformation and causes dramatic changes of electrostatic surface charge, both of which may affect TLR2 recognition and signaling.

Homology modeling and virtual screening for antagonists of protease from yellow head virus

Yellow head virus (YHV) is one of the causative agents of shrimp viral disease. The prevention of YHV infection in shrimp has been developed by various methods, but it is still insufficient to protect the mass mortality in shrimp. New approaches for the antiviral drug development for viral infection have been focused on the inhibition of several potent viral enzymes, and thus the YHV protease is one of the interesting targets for developing antiviral drugs according to the pivotal roles of the enzyme in an early stage of viral propagation. In this study, a theoretical modeling of the YHV protease was constructed based on the folds of several known crystal structures of other viral proteases, and was subsequently used as a target for virtual screening—molecular docking against approximately 1364 NCI structurally diversity compounds. A complex between the protease and the hit compounds was investigated for intermolecular interactions by molecular dynamics simulations. Five best predicted compounds (NSC122819, NSC345647, NSC319990, NSC50650, and NSC5069) were tested against bacterial expressed YHV. The NSC122819 showed the best inhibitory characteristic among the candidates, while others showed more than 50 % of inhibition in the assay condition. These compounds could potentially be inhibitors for curing YHV infection.

Two polymorphisms facilitate differences in plasticity between two chicken major histocompatibility complex class I proteins

Major histocompatibility complex class I molecules (MHC I) present peptides to cytotoxic T-cells at the surface of almost all nucleated cells. The function of MHC I molecules is to select high affinity peptides from a large intracellular pool and they are assisted in this process by co-factor molecules, notably tapasin. In contrast to mammals, MHC homozygous chickens express a single MHC I gene locus, termed BF2, which is hypothesised to have co-evolved with the highly polymorphic tapasin within stable haplotypes. The BF2 molecules of the B15 and B19 haplotypes have recently been shown to differ in their interactions with tapasin and in their peptide selection properties. This study investigated whether these observations might be explained by differences in the protein plasticity that is encoded into the MHC I structure by primary sequence polymorphisms. Furthermore, we aimed to demonstrate the utility of a complimentary modelling approach to the understanding of complex experimental data. Combining mechanistic molecular dynamics simulations and the primary sequence based technique of statistical coupling analysis, we show how two of the eight polymorphisms between BF2*15∶01 and BF2*19∶01 facilitate differences in plasticity. We show that BF2*15∶01 is intrinsically more plastic than BF2*19∶01, exploring more conformations in the absence of peptide. We identify a protein sector of contiguous residues connecting the membrane bound α3 domain and the heavy chain peptide binding site. This sector contains two of the eight polymorphic residues. One is residue 22 in the peptide binding domain and the other 220 is in the α3 domain, a putative tapasin binding site. These observations are in correspondence with the experimentally observed functional differences of these molecules and suggest a mechanism for how modulation of MHC I plasticity by tapasin catalyses peptide selection allosterically.

Atomistic mechanisms of huntingtin N-terminal fragment insertion on a phospholipid bilayer revealed by molecular dynamics simulations

The huntingtin protein is characterized by a segment of consecutive glutamines (Q(N)) that is responsible for its fibrillation. As with other amyloid proteins, misfolding of huntingtin is related to Huntington's disease through pathways that can involve interactions with phospholipid membranes. Experimental results suggest that the N-terminal 17-amino-acid sequence (htt(NT)) positioned just before the Q(N) region is important for the binding of huntingtin to membranes. Through all-atom explicit solvent molecular dynamics simulations, we unveil the structure and dynamics of the htt(NT)Q(N) fragment on a phospholipid membrane at the atomic level. We observe that the insertion dynamics of this peptide can be described by four main steps-approach, reorganization, anchoring, and insertion-that are very diverse at the atomic level. On the membrane, the htt(NT) peptide forms a stable α-helix essentially parallel to the membrane with its nonpolar side-chains-mainly Leu-4, Leu-7, Phe-11 and Leu-14-positioned in the hydrophobic core of the membrane. Salt-bridges involving Glu-5, Glu-12, Lys-6, and Lys-15, as well as hydrogen bonds involving Thr-3 and Ser-13 with the phospholipids also stabilize the structure and orientation of the htt(NT) peptide. These observations do not significantly change upon adding the Q(N) region whose role is rather to provide, through its hydrogen bonds with the phospholipids' head group, a stable scaffold facilitating the partitioning of the htt(NT) region in the membrane. Moreover, by staying accessible to the solvent, the amyloidogenic Q(N) region could also play a key role for the oligomerization of htt(NT)Q(N) on phospholipid membranes.

MD and NMR analyses of choline and TMA binding to duplex DNA: on the origins of aberrant sequence-dependent stability by alkyl cations in aqueous and water-free solvents

It has been known for decades that alkylammonium ions, such as tetramethyl ammonium (TMA), alter the usual correlation between DNA GC-content and duplex stability. In some cases it is even possible for an AT-rich duplex to be more stable than a GC-rich duplex of the same length. There has been much speculation regarding the origin of this aberration in sequence-dependent DNA duplex stability, but no clear resolution. Using a combination of molecular dynamics simulations and NMR spectroscopy we demonstrate that choline (2-hydroxy-N,N,N-trimethylethanaminium) and TMA are preferentially localized in the minor groove of DNA duplexes at A·T base pairs and these same ions show less pronounced localization in the major groove compared to what has been demonstrated for alkali and alkali earth metal ions. Furthermore, free energy calculations show that single-stranded GC-rich sequences exhibit more favorable solvation by choline than single-stranded AT-rich sequences. The sequence-specific nature of choline and TMA binding provides a rationale for the enhanced stability of AT-rich sequences when alkyl-ammonium ions are used as the counterions of DNA. Our combined theoretical and experimental study provides one of the most detailed pictures to date of cations localized along DNA in the solution state, and provides insights that go beyond understanding alkyl-ammonium ion binding to DNA. In particular, because choline and TMA bind to DNA in a manner that is found to be distinct from that previously reported for Na(+), K(+), Mg(2+), and Ca(2+), our results reveal the important but underappreciated role that most other cations play in sequence-specific duplex stability.

Structure of RecX protein complex with the presynaptic RecA filament: Molecular dynamics simulations and small angle neutron scattering

Using molecular modeling techniques we have built the full atomic structure and performed molecular dynamics simulations for the complexes formed by Escherichia coli RecX protein with a single-stranded oligonucleotide and with RecA presynaptic filament. Based on the modeling and SANS experimental data a sandwich-like filament structure formed two chains of RecX monomers bound to the opposite sides of the single stranded DNA is proposed for RecX::ssDNA complex. The model for RecX::RecA::ssDNA include RecX binding into the grove of RecA::ssDNA filament that occurs mainly via Coulomb interactions between RecX and ssDNA. Formation of RecX::RecA::ssDNA filaments in solution was confirmed by SANS measurements which were in agreement with the spectra computed from the molecular dynamics simulations.

Affinity comparison of different THCA synthase to CBGA using modeling computational approaches

The Δ^9-Tetrahydrocannabinol (THCA) is the primary psychoactive compound of Cannabis Sativa. It is produced by Δ^1- Tetrahydrocannabinolic acid synthase (THCA) which catalyzes the oxidative cyclization of cannabigerolic acid (CBGA) the precursor of the THCA. In this study, we were interested by the three dimensional structure of THCA synthase protein. Generation of models were done by MODELLER v9.11 and homology modeling with Δ1-tetrahydrocannabinolic acid (THCA) synthase X ray structure (PDB code 3VTE) on the basis of sequences retrieved from GenBank. Procheck, Errat, and Verify 3D tools were used to verify the reliability of the six 3D models obtained, the overall quality factor and the Prosa Z-score were also used to check the quality of the six modeled proteins. The RMSDs for C-alpha atoms, main-chain atoms, side-chain atoms and all atoms between the modeled structures and the corresponding template ranged between 0.290 Å-1.252 Å, reflecting the good quality of the obtained models. Our study of the CBGA-THCA synthase docking demonstrated that the active site pocket was successfully recognized using computational approach. The interaction energy of CBGA computed in ‘fiber types’ proteins ranged between -4.1 95 kcal/mol and -5.95 kcal/mol whereas in the ‘drug type’ was about -7.02 kcal/mol to -7.16 kcal/mol, which maybe indicate the important role played by the interaction energy of CBGA in the determination of the THCA level in Cannabis Sativa L. varieties. Finally, we have proposed an experimental design in order to explore the binding energy source of ligand-enzyme in Cannabis Sativa and the production level of the THCA in the absence of any information regarding the correlation between the enzyme affinity and THCA level production. This report opens the doors to more studies predicting the binding site pocket with accuracy from the perspective of the protein affinity and THCA level produced in Cannabis Sativa.

Virtual screening of integrase inhibitors by large scale binding free energy calculations: the SAMPL4 challenge

As part of the SAMPL4 blind challenge, filtered AutoDock Vina ligand docking predictions and large scale binding energy distribution analysis method binding free energy calculations have been applied to the virtual screening of a focused library of candidate binders to the LEDGF site of the HIV integrase protein. The computational protocol leveraged docking and high level atomistic models to improve enrichment. The enrichment factor of our blind predictions ranked best among all of the computational submissions, and second best overall. This work represents to our knowledge the first example of the application of an all-atom physics-based binding free energy model to large scale virtual screening. A total of 285 parallel Hamiltonian replica exchange molecular dynamics absolute protein-ligand binding free energy simulations were conducted starting from docked poses. The setup of the simulations was fully automated, calculations were distributed on multiple computing resources and were completed in a 6-weeks period. The accuracy of the docked poses and the inclusion of intramolecular strain and entropic losses in the binding free energy estimates were the major factors behind the success of the method. Lack of sufficient time and computing resources to investigate additional protonation states of the ligands was a major cause of mispredictions. The experiment demonstrated the applicability of binding free energy modeling to improve hit rates in challenging virtual screening of focused ligand libraries during lead optimization.

Decoding the structural events in substrate-gating mechanism of eukaryotic prolyl oligopeptidase using normal mode analysis and molecular dynamics simulations

Prolyl oligopeptidase (POP) is a serine protease, unique for its ability to cleave various small oligopeptides shorter than 30 amino acids. POP is an important drug target since it is implicated in various neurological disorders. Although there is structural evidence that bacterial POPs undergo huge interdomain movements and acquire an "open" state in the substrate-unbound form, hitherto, no crystal structure is available in the substrate-unbound domain-open form of eukaryotic POPs. Indeed, there is no difference between the substrate-unbound/bound states of eukaryotic POPs. This raises unanswered questions about whether difference in the substrate access pathway exists between bacterial and eukaryotic POPs. Here, we have used normal mode analysis and molecular dynamics to unravel the mechanism of substrate entry in mammalian POPs, which has been debated until now. Motions observed using normal modes of porcine and bacterial POPs were analyzed and compared, augmented by molecular dynamics of these proteins. Identical to bacterial POPs, interdomain opening was found to be the possible pathway for the substrate-gating in mammals as well. On the basis of our analyses and evidences, a mechanistic model of substrate entry in POPs has been proposed. Up-down movement of N-terminal hydrolase domain resulted in twisting motion of two domains, followed by the conformational changes of interdomain loop regions, which facilitate interdomain opening. Similar to bacterial POPs, an open form of porcine POP is also proposed with domain-closing motion. This work has direct implications for the development of novel inhibitors of mammalian POPs to understand the etiology of various neurological diseases.

Free energy calculations for molecular solids using GROMACS

In this work, we describe a procedure to evaluate the free energy of molecular solids with the GROMACS molecular dynamics package. The free energy is calculated using the Einstein molecule method that can be regarded as a small modification of the Einstein crystal method. Here, the position and orientation of the molecules is fixed by using an Einstein field that binds with harmonic springs at least three non-collinear atoms (or points of the molecule) to their reference positions. The validity of the Einstein field is tested by performing free-energy calculations of methanol, water (ice), and patchy colloids molecular solids. The free energies calculated with GROMACS show a very good agreement with those obtained using Monte Carlo and with previously published results.

Perspectives for hybrid ab initio/molecular mechanical simulations of solutions: from complex chemistry to proton-transfer reactions and interfaces

As a consequence of the ongoing development of enhanced computational resources, theoretical chemistry has become an increasingly valuable field for the investigation of a variety of chemical systems. Simulations employing a hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) technique have been shown to be a particularly promising approach, whenever ultrafast (i.e., picosecond) dynamical properties are to be studied, which are in many cases difficult to access via experimental techniques. Details of the quantum mechanical charge field (QMCF) ansatz, an advanced QM/MM protocol, are discussed and simulation results for various systems ranging from simple ionic hydrates to solvated organic molecules and coordination complexes in solution are presented. A particularly challenging application is the description of proton-transfer reactions in chemical simulations, which is a prerequisite to study acidified and basic systems. The methodical requirements for a combination of the QMCF methodology with a dissociative potential model for the description of the solvent are discussed. Furthermore, the possible extension of QM/MM approaches to solid/liquid interfaces is outlined.

Activation-triggered subunit exchange between CaMKII holoenzymes facilitates the spread of kinase activity

The activation of the dodecameric Ca(2+)/calmodulin dependent kinase II (CaMKII) holoenzyme is critical for memory formation. We now report that CaMKII has a remarkable property, which is that activation of the holoenzyme triggers the exchange of subunits between holoenzymes, including unactivated ones, enabling the calcium-independent phosphorylation of new subunits. We show, using a single-molecule TIRF microscopy technique, that the exchange process is triggered by the activation of CaMKII, and that exchange is modulated by phosphorylation of two residues in the calmodulin-binding segment, Thr 305 and Thr 306. Based on these results, and on the analysis of molecular dynamics simulations, we suggest that the phosphorylated regulatory segment of CaMKII interacts with the central hub of the holoenzyme and weakens its integrity, thereby promoting exchange. Our results have implications for an earlier idea that subunit exchange in CaMKII may have relevance for information storage resulting from brief coincident stimuli during neuronal signaling. DOI: http://dx.doi.org/10.7554/eLife.01610.001.

Silencing motifs in the Clr2 protein from fission yeast, Schizosaccharomyces pombe

The fission yeast, Schizosaccharomyces pombe, is a well-established model for heterochromatin formation, but the exact sequence of events for initiation remains to be elucidated. The essential factors involved include RNA transcribed from repeated sequences together with the methyltransferase Clr4. In addition, histone deacetylases, like Clr3, found in the SHREC complex are also necessary for transcriptional silencing. Clr2 is another crucial factor required for heterochromatin formation found in the SHREC complex. The function of Clr2 has been difficult to establish due to the lack of conserved domains or homology to proteins of known molecular function. Using a bioinformatics approach, three conserved motifs in Clr2 were identified, which contained amino acids important for transcriptional repression. Analysis of clr2 mutant strains revealed a major role for Clr2 in mating-type and rDNA silencing, and weaker effects on centromeric silencing. The effect on mating-type silencing showed variegation in several of the strains with mutated versions of Clr2 indicating an establishment or maintenance defect. Moreover, the critical amino acids in Clr2 were also necessary for transcriptional repression in a minimal system, by the tethering of Clr4 upstream of a reporter gene, inserted into the euchromatic part of the genome. Finally, in silico modeling suggested that the mutations in Clr2 cause disruption of secondary structures in the Clr2 protein. Identification of these critical amino acids in the protein provides a useful tool to explore the molecular mechanism behind the role of Clr2 in heterochromatin formation.