Force field dependence of riboswitch dynamics

Riboswitches are noncoding regulatory elements that control gene expression in response to the presence of metabolites, which bind to the aptamer domain. Metabolite binding appears to occur through a combination of conformational selection and induced fit mechanism. This demands to characterize the structural dynamics of the apo state of aptamer domains. In principle, molecular dynamics (MD) simulations can give insights at the atomistic level into the dynamics of the aptamer domain. However, it is unclear to what extent contemporary force fields can bias such insights. Here, we show that the Amber force field ff99 yields the best agreement with detailed experimental observations on differences in the structural dynamics of wild type and mutant aptamer domains of the guanine-sensing riboswitch (Gsw), including a pronounced influence of Mg2+. In contrast, applying ff99 with parmbsc0 and parmχOL modifications (denoted ff10) results in strongly damped motions and overly stable tertiary loop-loop interactions. These results are based on 58 MD simulations with an aggregate simulation time>11 μs, careful modeling of Mg2+ ions, and thorough statistical testing. Our results suggest that the moderate stabilization of the χ-anti region in ff10 can have an unwanted damping effect on functionally relevant structural dynamics of marginally stable RNA systems. This suggestion is supported by crystal structure analyses of Gsw aptamer domains that reveal χ torsions with high-anti values in the most mobile regions. We expect that future RNA force field development will benefit from considering marginally stable RNA systems and optimization toward good representations of dynamics in addition to structural characteristics.

Two heteronuclear dipolar results at the price of one: quantifying Na/P contacts in phosphosilicate glasses and biomimetic hydroxy-apatite

The analysis of S{I} recoupling experiments applied to amorphous solids yields a heteronuclear second moment M(2)(S-I) that represents the effective through-space dipolar interaction between the detected S spins and the neighboring I-spin species. We show that both M(2)(S-I) and M(2)(I-S) values are readily accessible from a sole S{I} or I{S} experiment, which may involve either S or I detection, and is naturally selected as the most favorable option under the given experimental conditions. For the common case where I has half-integer spin, an I{S} REDOR implementation is preferred to the S{I} REAPDOR counterpart. We verify the procedure by (23)Na{(31)P} REDOR and (31)P{(23)Na} REAPDOR NMR applied to Na(2)O-CaO-SiO(2)-P(2)O(5) glasses and biomimetic hydroxyapatite, where the M(2)(P-Na) values directly determined by REAPDOR agree very well with those derived from the corresponding M(2)(Na-P) results measured by REDOR. Moreover, we show that dipolar second moments are readily extracted from the REAPDOR NMR protocol by a straightforward numerical fitting of the initial dephasing data, in direct analogy with the well-established procedure to determine M(2)(S-I) values from REDOR NMR experiments applied to amorphous materials; this avoids the problems with time-consuming numerically exact simulations whose accuracy is limited for describing the dynamics of a priori unknown multi-spin systems in disordered structures.

Identification of chikungunya virus nsP2 protease inhibitors using structure-base approaches

The nsP2 protease of chikungunya virus (CHIKV) is one of the essential components of viral replication and it plays a crucial role in the cleavage of polyprotein precursors for the viral replication process. Therefore, it is gaining attention as a potential drug design target against CHIKV. Based on the recently determined crystal structure of the nsP2 protease of CHIKV, this study identified potential inhibitors of the virus using structure-based approaches with a combination of molecular docking, virtual screening and molecular dynamics (MD) simulations. The top hit compounds from database searching, using the NCI Diversity Set II, with targeting at five potential binding sites of the nsP2 protease, were identified by blind dockings and focused dockings. These complexes were then subjected to MD simulations to investigate the stability and flexibility of the complexes and to gain a more detailed insight into the interactions between the compounds and the enzyme. The hydrogen bonds and hydrophobic contacts were characterized for the complexes. Through structural alignment, the catalytic residues Cys1013 and His1083 were identified in the N-terminal region of the nsP2 protease. The absolute binding free energies were estimated by the linear interaction energy approach and compared with the binding affinities predicted with docking. The results provide valuable information for the development of inhibitors for CHIKV.

A novel fragment based strategy for membrane active antimicrobials against MRSA

Membrane active antimicrobials are a promising new generation of antibiotics that hold the potential to avert antibiotic resistance. However, poor understanding of the action mechanism and the lack of general design principles have impeded their development. Here we extend the concept of fragment based drug design and propose a pharmacophore model based on first principles for the design of membrane active antimicrobials against Gram positive pathogens. Elaborating on a natural xanthone-based hydrophobic scaffold, two derivatives of the pharmacophore model are proposed, and these demonstrate excellent antimicrobial activity. Rigorous molecular dynamics simulations combined with biophysical experiments suggest a three-step mechanism of action (absorption-translocation-disruption) which allows us to identify key factors for the practical optimization of each fragment of the pharmacophore. Moreover, the model matches the structures of several membrane active antimicrobials which are currently in clinical trials. Our model provides a novel and rational approach for the design of bactericidal molecules that target the bacterial membrane.

A salt-bridge stabilized C-terminal hook is critical for the dimerization of a Bowman Birk inhibitor

Legume Bowman-Birk inhibitors (BBIs) that inhibit mammalian proteases exist as dimers in solution. The structural basis governing dimerization of HGI-III (horsegram seed BBI) was investigated. An intra-monomer salt bridge (D76-K71) stabilizes an atypical hook-like conformation at the C-terminus. We postulate that this hook, positions D75 to enable an inter-monomer salt-bridge D75(a)-K24(b), which results in dimerization. We verify this by K71A and D76A mutations of HGI-III. The mutants were both monomers, likely due to destabilization of the C-terminal hook. Dimerization was sustained in a double mutant K71D/D76K that was anticipated to form a similar hook critical for dimerization. Conversely, K24(b) that interacts with D75(a) of the loop is the specificity determining residue that interacts with trypsin to inhibit its activity. The inter-monomer salt bridge D75(a)-K24(b) must be disrupted for the inhibition of trypsin, requiring HGI-III to transition into a monomer. Size exclusion studies and a model of HGI-III-trypsin complex support this notion. Interestingly, isoforms of the inhibitor present in germinated seeds (HGGIs) are monomers; and most strikingly, the C-termini of these inhibitors are truncated with the loss the C-terminal hook critical for dimerization. The tendency of HGI-III to self-associate seems to relate to its physiological function of a storage protein.

What can we learn from molecular dynamics simulations for GPCR drug design?

Recent years have seen a tremendous progress in the elucidation of experimental structural information for G-protein coupled receptors (GPCRs). Although for the vast majority of pharmaceutically relevant GPCRs structural information is still accessible only by homology models the steadily increasing amount of structural information fosters the application of structure-based drug design tools for this important class of drug targets. In this article we focus on the application of molecular dynamics (MD) simulations in GPCR drug discovery programs. Typical application scenarios of MD simulations and their scope and limitations will be described on the basis of two selected case studies, namely the binding of small molecule antagonists to the human CC chemokine receptor 3 (CCR3) and a detailed investigation of the interplay between receptor dynamics and solvation for the binding of small molecules to the human muscarinic acetylcholine receptor 3 (hM3R).

Molecular mechanism of the enhanced viral fitness contributed by secondary mutations in the hemagglutinin protein of oseltamivir resistant H1N1 influenza viruses: modeling studies of antibody and receptor binding

The envelope protein hemagglutinin (HA) of influenza viruses is primarily associated with host antibody and receptor interactions. The HA protein is known to maintain a functional balance with neuraminidase (NA), the other major envelope protein. Prior to 2007-2008, human seasonal H1N1 viruses possessing the NA H274Y mutation, which confers oseltamivir resistance, generally had low growth capability. Subsequently, secondary mutations that compensate for the deleterious effect of the NA H274Y mutation have been identified. The molecular mechanism of how the defect could be counteracted by these secondary mutations is not fully understood. We studied here the effect of three such mutations (T86K, K144E and R192K) in the HA protein, which are located at either the HA receptor binding site or in the H1N1 antigenic sites. Molecular docking and dynamics studies showed that, of the three mutations, the R192K mutation could have mediated neutralizing antibody escape and decreased receptor binding affinity, either or both of which may have contributed to increased viral fitness. The study suggests the molecular basis of enhanced viral fitness induced by secondary mutations in the evolution of oseltamivir-resistant influenza strains.

A theoretical view of the C3d:CR2 binding controversy

The C3d:CR2(SCR1-2) interaction plays an important role in bridging innate and adaptive immunity, leading to enhanced antibody production at sites of complement activation. Over the past decade, there has been much debate over the binding mode of this interaction. An initial cocrystal structure (PDB: 1GHQ) was published in 2001, in which the only interactions observed were between the SCR2 domain of CR2 and a side-face of C3d whereas a cocrystal structure (PDB: 3OED) published in 2011 showed both the SCR1 and SCR2 domains of CR2 interacting with an acidic patch on the concave surface of C3d. The initial 1GHQ structure is at odds with the majority of existing biochemical data and the publication of the 3OED structure renewed uncertainty regarding the physiological relevance of 1GHQ, suggesting that crystallization may have been influenced by the presence of zinc acetate in the crystallization process. In our study, we used a variety of computational approaches to gain insight into the binding mode between C3d and CR2 and demonstrate that the binding site at the acidic patch (3OED) is electrostatically more favorable, exhibits better structural and dissociative stability, specifically at the SCR1 domain, and has higher binding affinity than the 1GHQ binding mode. We also observe that nonphysiological zinc ions enhance the formation of the C3d:CR2 complex at the side face of C3d (1GHQ) through increases in electrostatic favorability, intermolecular interactions, dissociative character and overall energetic favorability. These results provide a theoretical basis for the association of C3d:CR2 at the acidic cavity of C3d and provide an explanation for binding of CR2 at the side face of C3d in the presence of nonphysiological zinc ions.

Molecular dynamics simulations: from structure function relationships to drug discovery

Molecular dynamics (MD) simulation is an emerging in silico technique with potential applications in diverse areas of pharmacology. Over the past three decades MD has evolved as an area of importance for understanding the atomic basis of complex phenomena such as molecular recognition, protein folding, and the transport of ions and small molecules across membranes. The application of MD simulations in isolation and in conjunction with experimental approaches have provided an increased understanding of protein structure-function relationships and demonstrated promise in drug discovery.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

How do the protonation states of E296 and D312 in OmpF and D299 and D315 in homologous OmpC affect protein structure and dynamics? Simulation studies

In this study, the structural and dynamic properties of two major porins (OmpF and OmpC) in Escherichia coli are investigated using molecular dynamics (MD) simulations. Both porins have the extracellular loop L3 folded halfway through the pore to form a constriction area. The solute influx and efflux are controlled by the L3 movement. E296 and D312 in OmpF and homologous D299 and D315 in OmpC located on the barrel wall are found to play a key role in L3 gating activity. All possible charged states of both E296(D299) and D312(315) are applied in this study to observe changes in overall structure and especially L3 movement. The results show that different protonation states of both residues cause the large-scale deviations in structure and pore cavity especially in OmpF. Fully charged E296(D299) and D312(315) increase the protein flexibility significantly. Deprotonating at least one of E296(D299) and D312(315) helps to fasten L3 to the barrel wall and maintain pore size. Lacking of interactions with D312(315) can lead to the pore closure in OmpF. Comparing with OmpC, not only is OmpF less stable, but it is also more sensitive to the charge states of both E296(D299) and D312(315).

Membrane association of the PTEN tumor suppressor: neutron scattering and MD simulations reveal the structure of protein-membrane complexes

Neutron reflection (NR) from planar interfaces is an emerging technology that provides unique and otherwise inaccessible structural information on disordered molecular systems such as membrane proteins associated with fluid bilayers, thus addressing one of the remaining challenges of structural biology. Although intrinsically a low-resolution technique, using structural information from crystallography or NMR allows the construction of NR models that describe the architecture of protein-membrane complexes at high resolution. In addition, a combination of these methods with molecular dynamics (MD) simulations has the potential to reveal the dynamics of protein interactions with the bilayer in atomistic detail. We review recent advances in this area by discussing the application of these techniques to the complex formed by the PTEN phosphatase with the plasma membrane. These studies provide insights in the cellular regulation of PTEN, its interaction with PI(4,5)P2 in the inner plasma membrane and the pathway by which its substrate, PI(3,4,5)P3, accesses the PTEN catalytic site.

Extension of the free energy workflow FEW towards implicit solvent/implicit membrane MM-PBSA calculations

BACKGROUND

The number of high-resolution structures of pharmacologically relevant membrane proteins has been strongly increasing. This makes computing relative affinities of chemically similar compounds binding to a membrane protein possible in order to guide decision making in drug design. However, the preparation step of such calculations is time-consuming and complex.

METHODS

We extended the free energy workflow tool FEW, available in AMBER, towards facilitating the setup of molecular dynamics simulations with explicit membrane, and the setup and execution of effective binding energy calculations according to a 1-trajectory implicit solvent/implicit membrane MM-PBSA approach for multiple ligands binding to the same membrane protein.

RESULTS

We validated the implemented protocol initially on two model systems, a sodium ion in the presence of an implicit membrane slab and a proton traversing the M2 proton-channel of the influenza A virus. For the latter, we found a good agreement for several important events along the proton pathway with those obtained in a recent computational study. Finally, we performed a case study on effective binding energy calculations for a set of inhibitors binding to the M2 proton-channel.

CONCLUSIONS

From the case study, we estimate a considerable speed up in the setup and analysis times for implicit solvent/implicit membrane MM-PBSA calculations by the extended version of FEW compared to a manual preparation.

GENERAL SIGNIFICANCE

Together with the overall runtime and the analysis results, this suggests that such type of calculations can be valuable in later stages of drug design projects on membrane proteins. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Two homologous neutrophil serine proteases bind to POPC vesicles with different affinities: When aromatic amino acids matter

Neutrophil serine proteases Proteinase 3 (PR3) and human neutrophil elastase (HNE) are homologous antibiotic serine proteases of the polymorphonuclear neutrophils. Despite sharing a 56% sequence identity they have been shown to have different functions and localizations in the neutrophils. In particular, and in contrast to HNE, PR3 has been detected at the outer leaflet of the plasma membrane and its membrane expression is a risk factor in a number of chronic inflammatory diseases. Although a plethora of studies performed in various cell-based assays have been reported, the mechanism by which PR3, and possibly HNE bind to simple membrane models remains unclear. We used surface plasmon resonance (SPR) experiments to measure and compare the affinity of PR3 and HNE for large unilamellar vesicles composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). We also conducted 500-nanosecond long molecular dynamics simulations of each enzyme at the surface of a POPC bilayer to map the interactions between proteins and lipids and rationalize the difference in affinity observed in the SPR experiment. We find that PR3 binds strongly to POPC large unilamellar vesicles (Kd=9.2×10(-7)M) thanks to the insertion of three phenylalanines, one tryptophan and one leucine beyond the phosphate groups of the POPC lipids. HNE binds in a significantly weaker manner (Kd>10(-5)M) making mostly electrostatic interactions via lysines and arginines and inserting only one leucine between the hydrophobic lipid tails. Our results support the early reports that PR3, unlike HNE, is able to directly and strongly anchor directly to the neutrophil membrane.

Grape anthocyanin oligomerization: a putative mechanism for red color stabilization?

The equilibrium forms of malvidin-3-O-glucoside trimer present in grape skins were studied in aqueous solution at different pH values through UV-Visible spectroscopy. It was observed that the reactivity of this compound is strongly dominated by acid-base chemistry (pKa1=3.61±0.03; pKa2=6.83±0.06), with the reaction sequence hydration-tautomerization-isomerization accounting less than 10% of the overall reactivity. The decrease of hydration of this flavylium cation derivative when compared to the original anthocyanins results from the formation of a cluster of the pigment with a high-energy of solvation that inhibits the access of water molecules to the flavylium cation core preventing by this way the hydration reactions. Overall of these results raise the hypothesis that polymerization may be a natural stabilization mechanism for the red color of anthocyanins.

Discovery of Potential, Non-Toxic Influenza Virus Inhibitor by Computational Techniques

Influenza infection continues to be a major problem in many parts of the world. Rimantadine is a first-line drug used to treat the influenza infection by targeting M2 proton channel. However, S31N mutation in M2 proton channel diminishes the efficiency of rimantadine and creates resistance. To address this issue, the present study was aimed to screen the effective lead candidate against drug resistance strain of influenza from DrugBank database. Initially, the lead molecules were filtered using Lipinski rule of five and the drug likeliness property. Subsequently, the data reduction was carried out by employing molecular docking study. Finally, molecular dynamics simulations techniques were performed to validate the lead compound. Most importantly, the -p LD50 of the screened lead molecule was calculated using CORAL software to estimate the Rat oral toxicity. Accordingly, memantine may possibly become a promising lead compound of rimantadine-resistant influenza virus strain.

Detailed analysis and follow-up studies of a high-throughput screening for indoleamine 2,3-dioxygenase 1 (IDO1) inhibitors

Indoleamine 2,3-dioxygenase 1 (IDO1) is a key regulator of immune responses and therefore an important therapeutic target for the treatment of diseases that involve pathological immune escape, such as cancer. Here, we describe a robust and sensitive high-throughput screen (HTS) for IDO1 inhibitors using the Prestwick Chemical Library of 1200 FDA-approved drugs and the Maybridge HitFinder Collection of 14,000 small molecules. Of the 60 hits selected for follow-up studies, 14 displayed IC50 values below 20 μM under the secondary assay conditions, and 4 showed an activity in cellular tests. In view of the high attrition rate we used both experimental and computational techniques to identify and to characterize compounds inhibiting IDO1 through unspecific inhibition mechanisms such as chemical reactivity, redox cycling, or aggregation. One specific IDO1 inhibitor scaffold, the imidazole antifungal agents, was chosen for rational structure-based lead optimization, which led to more soluble and smaller compounds with micromolar activity.

Molecular dynamics study on conformational differences between dGMP and 8-oxo-dGMP: Effects of metal ions

The modified nucleotide base 7,8-dihydro-8-oxo-guanine (8-oxo-G) is one of the major sources of spontaneous mutagenesis. Nucleotide-sanitizing enzymes, such as the MutT homolog-1 (MTH1) and nudix-type motif 5 (NUDT5), selectively remove 8-oxo-G from the cellular pool of nucleotides. Previous studies showed that, although the syn conformation generally predominates in purine nucleotides with a bulky substituent at the 8-position, 8-oxo-dGMP binds to both MTH1 and NUDT5 in the anti conformation. This study was initiated to investigate the possibility that 8-oxo-dGMP itself may adopt the anti conformation. Molecular dynamics simulations of mononucleotides (dGMP, 8-oxo-dGMP) in aqueous solution were performed. 8-oxo-dGMP adopted the anti conformation as well as the syn conformation, and the proportion of adopting the anti conformation increased in the presence of metal ions. When 8-oxo-dGMP was in the anti conformation, a metal ion was located between the oxygen atom of phosphate and the oxygen atom at the 8-position of 8-oxo-G. The types of stable anti conformations of 8-oxo-dGMP differed, depending on the ionic radii and charges of coexisting ions. These data suggested a role for metal ions, other than as cofactors for the hydrolysis of the di- and tri-phosphate forms of mononucleotides; that the metal ions help retain the anti conformation of the N-glycosidic torsion angle of 8-oxo-dGMP to promote the binding between the 8-oxo-G deoxynucleotide and the nucleotide-sanitizing enzymes.

Molecular mechanisms in the selective basal activation of pyrabactin receptor 1: Comparative analysis of mutants

Pyrabactin receptors (PYR) play a central role in abscisic acid (ABA) signal transduction; they are ABA receptors that inhibit type 2C protein phosphatases (PP2C). Molecular aspects contributing to increased basal activity of PYR against PP2C are studied by molecular dynamics (MD) simulations. An extensive series of MD simulations of the apo-form of mutagenized PYR1 as a homodimer and in complex with homology to ABA-insensitive 1 (HAB1) phosphatase are reported. In order to investigate the detailed molecular mechanisms mediating PYR1 activity, the MD data was analyzed by essential collective dynamics (ECD), a novel approach that allows the identification, with atomic resolution, of persistent dynamic correlations based on relatively short MD trajectories. Employing the ECD method, the effects of select mutations on the structure and dynamics of the PYR1 complexes were investigated and considered in the context of experimentally determined constitutive activities against HAB1. Approaches to rationally design constitutively active PYR1 constructs to increase PP2C inhibition are discussed.

Structural and energetic analysis to provide insight residues of CYP2C9, 2C11 and 2E1 involved in valproic acid dehydrogenation selectivity

Docking and molecular dynamics (MD) simulation have been two computational techniques used to gain insight about the substrate orientation within protein active sites, allowing to identify potential residues involved in the binding and catalytic mechanisms. In this study, both methods were combined to predict the regioselectivity in the binding mode of valproic acid (VPA) on three cytochrome P-450 (CYP) isoforms CYP2C9, CYP2C11, and CYP2E1, which are involved in the biotransformation of VPA yielding reactive hepatotoxic intermediate 2-n-propyl-4-pentenoic acid (4nVPA). There are experimental data about hydrogen atom abstraction of the C4-position of VPA to yield 4nVPA, however, there are not structural evidence about the binding mode of VPA and 4nVPA on CYPs. Therefore, the complexes between these CYP isoforms and VPA or 4nVPA were studied to explore their differences in binding and energetic stabilization. Docking results showed that VPA and 4nVPA are coupled into CYPs binding site in a similar conformation, but it does not explain the VPA hydrogen atom abstraction. On the other hand, MD simulations showed a set of energetic states that reorient VPA at the first ns, then making it susceptible to a dehydrogenation reaction. For 4nVPA, multiple binding modes were observed in which the different states could favor either undergo other reaction mechanism or ligand expulsion from the binding site. Otherwise, the energetic and entropic contribution point out a similar behavior for the three CYP complexes, showing as expected a more energetically favorable binding free energy for the complexes between CYPs and VPA than with 4nVPA.

Interactions of amikacin with the RNA model of the ribosomal A-site: computational, spectroscopic and calorimetric studies

Amikacin is a 2-deoxystreptamine aminoglycoside antibiotic possessing a unique l-HABA (l-(-)-γ-amino-α-hydroxybutyric acid) group and applied in the treatment of hospital-acquired infections. Amikacin influences bacterial translation by binding to the decoding region of the small ribosomal subunit that overlaps with the binding site of aminoacylated-tRNA (A-site). Here, we have characterized thermodynamics of interactions of amikacin with a 27-mer RNA oligonucleotide mimicking the aminoglycoside binding site in the bacterial ribosome. We applied isothermal titration and differential scanning calorimetries, circular dichroism and thermal denaturation experiments, as well as computer simulations. Thermal denaturation studies have shown that amikacin affects only slightly the melting temperatures of the A-site mimicking RNA model suggesting a moderate stabilization of RNA by amikacin. Isothermal titration calorimetry gives the equilibrium dissociation constants for the binding reaction between amikacin and the A-site oligonucleotide in the micromolar range with a favorable enthalpic contribution. However, for amikacin we observe a positive entropic contribution to binding, contrary to other aminoglycosides, paromomycin and ribostamycin. Circular dichroism spectra suggest that the observed increase in entropy is not caused by structural changes of RNA because amikacin binding does not destabilize the helicity of the RNA model. To investigate the origins of this positive entropy change we performed all-atom molecular dynamics simulations in explicit solvent for the 27-mer RNA oligonucleotide mimicking one A-site and the crystal structure of an RNA duplex containing two A-sites. We observed that the diversity of the conformational states of the l-HABA group sampled in the simulations of the complex was larger than for the free amikacin in explicit water. Therefore, the larger flexibility of the l-HABA group in the bound form may contribute to an increase of entropy upon binding.

Conformational changes in dopamine transporter intracellular regions upon cocaine binding and dopamine translocation

The dopamine transporter (DAT), a member of the neurotransmitter:sodium symporter family, mediates the reuptake of dopamine at the synaptic cleft. DAT is the primary target for psychostimulants such as cocaine and amphetamine. We previously demonstrated that cocaine binding and dopamine transport alter the accessibility of Cys342 in the third intracellular loop (IL3). To study the conformational changes associated with the functional mechanism of the transporter, we made cysteine substitution mutants, one at a time, from Phe332 to Ser351 in IL3 of the background DAT construct, X7C, in which 7 endogenous cysteines were mutated. The accessibility of the 20 engineered cysteines to polar charged sulfhydryl reagents was studied in the absence and presence of cocaine or dopamine. Of the 11 positions that reacted with methanethiosulfonate ethyl ammonium, as evidenced by inhibition of ligand binding, 5 were protected against this inhibition by cocaine and dopamine (S333C, S334C, N336C, M342C and T349C), indicating that reagent accessibility is affected by conformational changes associated with inhibitor and substrate binding. In some of the cysteine mutants, transport activity is disrupted, but can be rescued by the presence of zinc, most likely because the distribution between inward- and outward-facing conformations is restored by zinc binding. The experimental data were interpreted in the context of molecular models of DAT in both the inward- and outward-facing conformations. Differences in the solvent accessible surface area for individual IL3 residues calculated for these states correlate well with the experimental accessibility data, and suggest that protection by ligand binding results from the stabilization of the outward-facing configuration. Changes in the residue interaction networks observed from the molecular dynamics simulations also revealed the critical roles of several positions during the conformational transitions. We conclude that the IL3 region of DAT undergoes significant conformational changes in transitions necessary for both cocaine binding and substrate transport.

Implicit membrane treatment of buried charged groups: application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Electrophoretic mobilities of neutral analytes and electroosmotic flow markers in aqueous solutions of Hofmeister salts

Small neutral organic compounds have traditionally the role of EOF markers in electrophoresis, as they are expected to have zero electrophoretic mobility in external electric fields. The BGE contains, however, ions that have unequal affinities to the neutral molecules, which in turn results in their mobilization. In this study we focused on two EOF markers-thiourea and DMSO, as well as on N-methyl acetamide (NMA) as a model of the peptide bond. By means of CE and all atom molecular dynamics simulations we explored mobilization of these neutral compounds in large set of Hofmeister salts. Employing a statistical mechanics approach, we were able to reproduce by simulations the experimental electrophoretic mobility coefficients. We also established the role of the chemical composition of marker and the BGE on the measured electrophoretic mobility coefficient. For NMA, we interpreted the results in terms of the relative affinities of cations versus anions to the peptide bond.

Current state-of-the-art molecular dynamics methods and applications

Molecular dynamics simulations are used to describe the patterns, strength, and properties of protein behavior, drug-receptor interactions, the solvation of molecules, the conformational changes that a protein or molecule may undergo under various conditions, and other events that require the systematic evaluation of molecular properties in dynamic molecular systems. Only few years ago proteins were considered to be rigid body structures with very limited conformational flexibility. However, it is now clear that proteins are highly dynamic structures, the internal organization of which is the key to their 3D spatial arrangement and hence biological function. The study of protein dynamics in the lab is a very complicated, expensive, and time-consuming process. Therefore, a lot of effort and hope lies with the computers and the in silico study of protein structure and molecular dynamics. Herein, an effort has been made to describe the ever-evolving field of molecular dynamics, the different algorithms, and force fields that are being used as well as to provide some insight on what the near future holds for this auspicious field of computational structural biology.

P-I class metalloproteinase from Bothrops moojeni venom is a post-proline cleaving peptidase with kininogenase activity: insights into substrate selectivity and kinetic behavior

Snake venom metalloproteinases (SVMPs) belonging to P-I class are able to hydrolyze extracellular matrix proteins and coagulation factors triggering local and systemic reactions by multiple molecular mechanisms that are not fully understood. BmooMPα-I, a P-I class SMVP from Bothrops moojeni venom, was active upon neuro- and vaso-active peptides including angiotensin I, bradykinin, neurotensin, oxytocin and substance P. Interestingly, BmooMPα-I showed a strong bias towards hydrolysis after proline residues, which is unusual for most of characterized peptidases. Moreover, the enzyme showed kininogenase activity similar to that observed in plasma and cells by kallikrein. FRET peptide assays indicated a relative promiscuity at its S2-S'2 subsites, with proline determining the scissile bond. This unusual post-proline cleaving activity was confirmed by the efficient hydrolysis of the synthetic combinatorial library MCA-GXXPXXQ-EDDnp, described as resistant for canonical peptidases, only after Pro residues. Structural analysis of the tripeptide LPL complexed with BmooMPα-I, generated by molecular dynamics simulations, assisted in defining the subsites and provided the structural basis for subsite preferences such as the restriction of basic residues at the S2 subsite due to repulsive electrostatic effects and the steric impediment for large aliphatic or aromatic side chains at the S1 subsite. These new functional and structural findings provided a further understanding of the molecular mechanisms governing the physiological effects of this important class of enzymes in envenomation process.

Prediction of the interaction between spermidine and the G-G mismatch containing acceptor stem in tRNA(Ile): molecular modeling, density functional theory, and molecular dynamics study

Polyamines, putrescine, spermidine (SPD), and spermine are closely linked to cell growth, and highly regulate the levels of transcription, translation and protein turnover. We propose that SPD stimulates the formation of Ile-tRNA(Ile) by inducing a selective structural change of the G-G mismatch containing acceptor stem in tRNA(Ile). Here, we provide insight into how SPD recognizes and stabilizes the G-G mismatch containing acceptor stem in tRNA(Ile) with molecular modeling (MM), density functional theory (DFT) calculations, and molecular dynamics (MD) simulations. The results of the MM and DFT calculations indicate that the negatively charged region of the G-G mismatch containing acceptor stem in tRNA(Ile) is preferentially recognized by positively charged SPD. In addition, MD simulations indicate that all of the positively charged amino groups of SPD under physiological conditions (N1(NH3(+)), N5(NH2(+)), and N10(NH3(+)) could form hydrogen bonds with tRNA(Ile) and trigger the SPD-induced stabilization and structural change of the G-G mismatch containing acceptor stem in tRNA(Ile). Thus, this approach should be useful for determining the preferential binding site and appropriate binding mode of polyamines on tRNA(Ile).

Qualitative computational bioanalytics: assembly of viral channel-forming peptides around mono and divalent ions

A fine-grained docking protocol was used to generate a bundle-like structure of the bitopic membrane protein Vpu from HIV-1. Vpu is a type I membrane protein with 81 amino acids. It is proposed that Vpu forms ion- and substrate-conducting bundles, which are located at the plasma membrane in the infected cell. The Vpu1-32 peptide that includes the transmembrane domain (TMD) is assembled into homo-pentameric bundles around prepositioned Na, K, Ca or Cl ions. For bundles with the lowest energy, the TMDs generate a hydrophobic pore. Bundles in which Ser-24 faces the pore have higher energy. The tilt of the helices in the lowest energy bundles is larger than bundles with serines facing the pore. Left-handed bundles are lowest in energy where the ions are located at the serines.

Testing synthetic amyloid-β aggregation inhibitor using single molecule atomic force spectroscopy

Alzheimer's disease is a neurodegenerative disease with no known cure and few effective treatment options. The principal neurotoxic agent is an oligomeric form of the amyloid-β peptide and one of the treatment options currently being studied is the inhibition of amyloid aggregation. In this work, we test a novel pseudopeptidic aggregation inhibitor designated as SG1. SG1 has been designed to bind at the amyloid-β self-recognition site and prevent amyloid-β from misfolding into β sheet. We used atomic force spectroscopy, a nanoscale measurement technique, to quantify the binding forces between two single amyloid peptide molecules. For the first time, we demonstrate that single molecule atomic force spectroscopy can be used to assess the effectiveness of amyloid aggregation inhibitors by measuring the experimental yield of binding and can potentially be used as a screening technique for quick testing of efficacy of inhibitor drugs for amyloid aggregation.

Homology Models and Molecular Dynamics Simulations of Main Proteinase from Coronavirus Associated with Severe Acute Respiratory Syndrome (SARS)

In this study, two structural models (denoted as M^proST and M^proSH) of the main proteinase (M^pro) from the novel coronavirus associated with severe acute respiratory syndrome (SARS‐CoV) were constructed based on the crystallographic structures of M^pro from transmissible gastroenteritis coronavirus (TGEV) (M^proT) and human coronavirus HcoV‐229E (M^proH), respectively. Various 200 ps molecular dynamics simulations were subsequently performed to investigate the dynamics behaviors of several structural features. Both M^proST and M^proSH exhibit similar folds as their respective template proteins. These structural models reveal three distinct functional domains as well as an intervening loop connecting domains II and III as found in both template proteins. In addition, domain III of these structures exhibits the least secondary structural conservation. A catalytic cleft containing the substrate binding subsites S1 and the S2 between domains I and II are also observed in these structural models. Although these structures share many common features, the most significant difference occurs at the S2 subsite, where the amino acid residues lining up this subsite are least conserved. It may be a critical challenge for designing anti‐SARS drugs by simply screening the known database of proteinase inhibitors.

A thermodynamic approach to alamethicin pore formation

The structure and energetics of alamethicin Rf30 monomer to nonamer in cylindrical pores of 5 to 11Å radius are investigated using molecular dynamics simulations in an implicit membrane model that includes the free energy cost of acyl chain hydrophobic area exposure. Stable, low energy pores are obtained for certain combinations of radius and oligomeric number. The trimer and the tetramer formed 6Å pores that appear closed while the larger oligomers formed open pores at their optimal radius. The hexamer in an 8Å pore and the octamer in an 11Å pore give the lowest effective energy per monomer. However, all oligomers beyond the pentamer have comparable energies, consistent with the observation of multiple conductance levels. The results are consistent with the widely accepted "barrel-stave" model. The N terminal portion of the molecule exhibits smaller tilt with respect to the membrane normal than the C terminal portion, resulting in a pore shape that is a hybrid between a funnel and an hourglass. Transmembrane voltage has little effect on the structure of the oligomers but enhances or decreases their stability depending on its orientation. Antiparallel bundles are lower in energy than the commonly accepted parallel ones and could be present under certain experimental conditions. Dry aggregates (without an aqueous pore) have lower average effective energy than the corresponding aggregates in a pore, suggesting that alamethicin pores may be excited states that are stabilized in part by voltage and in part by the ion flow itself.

Correlation of biological activity with computationally derived structural features from transmembrane hetero-dimers of HIV-1 Vpu with host factors

Vpu is an 81 amino acid type I integral membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). It is identified to support viral release by potentially forming ion and substrate conducting channels and by modulating the function of host factors. The focus is on the interaction of the transmembrane domains of Vpu with those of host factors using a combination of molecular dynamics simulations and docking approach. Binding poses and adopted tilt angles of the dimers are analyzed and correlated with experimentally derived activity data from literature. Vpu activity is driven by dimerization with the host protein via its alanine rim Ala-8/11/15/19. Tight binding is shown by an almost parallel alignment of the helices in the dimers. Less parallel alignment is proposed to correlate with lower activity. This article is part of a Special Issue entitled: Viral Membrane Proteins - Channels for Cellular Networking.

Protonation of trimethylamine N-oxide (TMAO) is required for stabilization of RNA tertiary structure

The osmolyte trimethylamine N-oxide (TMAO) stabilizes the tertiary but not the secondary structures of RNA. However, molecular dynamics simulations performed on the PreQ1 riboswitch showed that TMAO destabilizes the tertiary riboswitch structure, leading us to hypothesize that the presence of RNA could result in enhanced population of the protonated form, TMAOP. Constant pH replica exchange simulations showed that a percentage of TMAO is indeed protonated, thus contributing to the stability of the tertiary but not the secondary structure of PreQ1. TMAOP results in an unfavorable dehydration of phosphodiester backbone, which is compensated by electrostatic attraction between TMAOP and the phosphate groups. In addition, TMAOP interacts with specific sites in the tertiary RNA structure, mimicking the behavior of positively charged ions and of the PreQ1 ligand in stabilizing RNA. Finally, we predict that TMAO-induced stabilization of RNA tertiary structures should be strongly pH dependent.

Computational modeling of the p7 monomer from HCV and its interaction with small molecule drugs

Hepatitis C virus p7 protein is a 63 amino acid polytopic protein with two transmembrane domains (TMDs) and one of the prime targets for anti HCV drug development. A bio-inspired modeling pathway is used to generate plausible computational models of the two TMDs forming the monomeric protein model. A flexible region between Leu-13 and Gly-15 is identified for TMD1_1-32 and a region around Gly-46 to Trp-48 for TMD2_36-58. Mutations of the tyrosine residues in TMD2_36-58 into phenylalanine and serine are simulated to identify their role in shaping TMD2. Lowest energy structures of the two TMDs connected with the loop residues are used for a posing study in which small molecule drugs BIT225, amantadine, rimantadine and NN-DNJ, are identified to bind to the loop region. BIT225 is identified to interact with the backbone of the functionally important residues Arg-35 and Trp-36.

Insights from the docking and molecular dynamics simulation of the Phosphopantetheinyl transferase (PptT) structural model from Mycobacterium tuberculosis

A great challenge is posed to the treatment of tuberculosis due to the evolution of multidrug-resistant (MDR) and extensively drugresistant (XDR) strains of Mycobacterium tuberculosis in recent times. The complex cell envelope of the bacterium contains unusual structures of lipids which protects the bacterium from host enzymes and escape immune response. To overcome the drug resistance, targeting "drug targets" which have a critical role in growth and virulence factor is a novel approach for better tuberculosis treatment. The enzyme Phosphopantetheinyl transferase (PptT) is an attractive drug target as it is primarily involved in post translational modification of various types-I polyketide synthases and assembly of mycobactin, which is required for lipid virulence factors. Our in silico studies reported that the structural model of M.tuberculosis PptT characterizes the structure-function activity. The refinement of the model was carried out with molecular dynamics simulations and was analyzed with root mean square deviation (RMSD), and radius of gyration (Rg). This confirmed the structural behavior of PptT in dynamic system. Molecular docking with substrate coenzyme A (CoA) identified the binding pocket and key residues His93, Asp114 and Arg169 involved in PptT-CoA binding. In conclusion, our results show that the M.tuberculosis PptT model and critical CoA binding pocket initiate the inhibitor design of PptT towards tuberculosis treatment.

Modeling, docking and dynamics simulations of a non-specific lipid transfer protein from Peganum harmala L

Non-specific lipid transfer proteins (ns-LTPs), ubiquitously found in various types of plants, have been well-known to transfer amphiphilic lipids and promote the lipid exchange between mitochondria and microbody. In this study, an in silico analysis was proposed to study ns-LTP in Peganum harmala L., which may belong to ns-LTP1 family, aiming at constructing its three-dimensional structure. Moreover, we adopted MEGA to analyze ns-LTPs and other species phylogenetically, which brought out an initial sequence alignment of ns-LTPs. In addition, we used molecular docking and molecular dynamics simulations to further investigate the affinities and stabilities of ns-LTP with several ligands complexes. Taken together, our results about ns-LTPs and their ligand-binding activities can provide a better understanding of the lipid-protein interactions, indicating some future applications of ns-LTP-mediated transport.

Molecular dynamics simulations of isoleucine-release pathway in GAF domain of N-CodY from Bacillus Subtilis

The GAF domain located in the N-terminal motifs of CodY (N-CodY) is responsible for increasing the affinity of CodY to its target sites on DNA by its interaction with the branched chain amino acids (BCAAs) involving isoleucine, leucine and valine. The study of the interaction of GAF domain with isoleucine gains much attention in recent years, but the mechanism of isoleucine release still remains unclear. In this paper, a conventional molecular dynamics (MD) and force probe molecular dynamics (FPMD) simulations have been performed with the aim to understand how the isoleucine ligand escapes from the GAF domain of N-CodY from Bacillus subtilis. The MD results reveal that the ligand release is a gradual process, which is accompanied by the movement of the loop between β3 and β4. During the periods of ligand escaping from the bottom to the top of binding pocket, isoleucine forms hydrogen bonds one after another with series of residues, such as ARG61, THR96, PHE98, VAL100, GLU101 and ASN102, under the mediation of hydrophobic contacts. The FPMD results show that the easiest way to pull ligand out of the cavity is along x direction (i.e. the direction is opposite to MET62).

Interpretation of 1H and 2H spin-lattice relaxation dispersions: insights from molecular dynamics simulations of polymer melts

We demonstrate that molecular dynamics simulations are a versatile tool to ascertain the interpretation of spin-lattice relaxation data. For (1)H, our simulation approach allows us to separate and to compare intra- and inter-molecular contributions to spin-lattice relaxation dispersions. Dealing with the important example of polymer melts, we show that the intramolecular parts of (1)H spectral densities and correlation functions are governed by rotational motion, while their inter-molecular counterparts provide access to translational motion, in particular, to mean-squared displacements and self-diffusion coefficients. Exploiting that the full microscopic information is available from molecular dynamics simulations, we determine the range of validity of experimental approaches, which often assume Gaussian dynamics, and we provide guidelines for the determination of free parameters required in experimental analyses. For (2)H, we examine the traditional methodology to extract correlation times of complex dynamics from relaxation data. Furthermore, based on knowledge from our computational study, it is shown that measurement of (2)H spin-lattice relaxation dispersions allows one to disentangle the intra- and inter-molecular contributions to the corresponding (1)H data in experimental work. Altogether, our simulation results yield a solid basis for future (1)H and (2)H spin-lattice relaxation analysis.

Network analysis to uncover the structural communication in GPCRs

Protein structure network (PSN) analysis is one of the graph theory-based approaches currently used to investigate the structural communication in biomolecular systems. Information on system dynamics can be provided by atomistic molecular dynamics simulations or coarse-grained Elastic Network Models paired with Normal Mode Analysis (ENM-NMA). This chapter describes the application of PSN analysis to uncover the structural communication in G protein-coupled receptors (GPCRs). Strategies to highlight changes in structural communication upon misfolding mutations, dimerization, and activation are described. Focus is put on the ENM-NMA-based strategy applied to the crystallographic structures of rhodopsin in its inactive (dark) and signaling active (meta II (MII)) states, highlighting clear changes in the PSN and the centrality of the retinal chromophore in differentiating the inactive and active states of the receptor.

Conformational flexibility of N-glycans in solution studied by REMD simulations

Protein-glycan recognition regulates a wide range of biological and pathogenic processes. Conformational diversity of glycans in solution is apparently incompatible with specific binding to their receptor proteins. One possibility is that among the different conformational states of a glycan, only one conformer is utilized for specific binding to a protein. However, the labile nature of glycans makes characterizing their conformational states a challenging issue. All-atom molecular dynamics (MD) simulations provide the atomic details of glycan structures in solution, but fairly extensive sampling is required for simulating the transitions between rotameric states. This difficulty limits application of conventional MD simulations to small fragments like di- and tri-saccharides. Replica-exchange molecular dynamics (REMD) simulation, with extensive sampling of structures in solution, provides a valuable way to identify a family of glycan conformers. This article reviews recent REMD simulations of glycans carried out by us or other research groups and provides new insights into the conformational equilibria of N-glycans and their alteration by chemical modification. We also emphasize the importance of statistical averaging over the multiple conformers of glycans for comparing simulation results with experimental observables. The results support the concept of "conformer selection" in protein-glycan recognition.

Restrained electrostatic potential atomic partial charges for condensed-phase simulations of carbohydrates

Charges derived from fitting a classical Coulomb model to quantum mechanical molecular electrostatic potentials (so called ESP-charges) are frequently used in simulations of macromolecules. Simulational methods that use ESP-charges generally reproduce the geometries of hydrogen bonded complexes, despite the fact that these charges are known to overestimate the strengths of these interactions. Through the use of a restraint function during the fitting of the partial charges to the electrostatic potentials the magnitudes of the charges may be attenuated (so called RESP-charges). For the AMBER force field RESP-charges have been proposed for proteins and nucleic acids. Here we examine a novel approach for determining the RESP-charges for carbohydrates based on molecular dynamics (MD) simulations of crystal structures. During a simulation, the crystallographic unit cell geometry is sensitive to both inter-molecular non-bonded forces and internal torsional rotations. However, for polar molecules, and specifically carbohydrates, the crystal geometries are particularly sensitive to the set of partial atomic charges employed in the simulation. Thus, given a force field in which the van der Waals and torsion terms are well parameterized, it is possible to assess the suitability of a set of partial charges by monitoring the properties of the crystal during an MD simulation. We have examined several charge sets for use with the GLYCAM parameters for carbohydrate and glycoprotein simulations and found that a restraint weight of 0.01 gives the best agreement with the neutron diffraction structure of α-d-glucopyranose. Unrestrained ESP-charges performed poorly as did the charges obtained from Mulliken and distributed multipole analyses of the quantum mechanical HF/6-31G* wavefunctions.