A molecular dynamics description of the conformational flexibility of the L-iduronate ring in glycosaminoglycans

For a synthetic hexasaccharide model it is shown that the conformational flexibility of the L-iduronate ring in glycosaminoglycans can be adequately described by using the PME methodology together with simulation protocols suitable for highly charged systems.     

The future of molecular dynamics simulations in drug discovery

Molecular dynamics simulations can now track rapid processes—those occurring in less than about a millisecond—at atomic resolution for many biologically relevant systems. These simulations appear poised to exert a significant impact on how new drugs are found, perhaps even transforming the very process of drug discovery. We predict here future results we can expect from, and enhancements we need to make in, molecular dynamics simulations over the coming 25 years, and in so doing set out several Grand Challenges for the field. In the context of the problems now facing the pharmaceutical industry, we ask how we can best address drug discovery needs of the next quarter century using molecular dynamics simulations, and we suggest some possible approaches.     

The conformational flexibility of nucleic acid bases paired in gas phase: a Car-Parrinello molecular dynamics study

The flexibilities of pyrimidine and imidazole rings in the paired nucleobases are investigated using Car-Parrinello molecular dynamics simulation in gas phase. The pairing influence on the stiffness of rings is analyzed based on the molecular structure of the nucleobases and constraints caused by pairing. We prove that the flexibilities of pyrimidine rings in isolated state have subtle correlation with the degree of aromaticity of the rings. The pairings in nucleic base pairs cause the rings to be more rigid for G, T, and U but more flexible for A and the same for C.     

Calculation of ligand binding free energies from molecular dynamics simulations

A recently developed method for predicting binding affinities in ligand–receptor complexes, based on interaction energy averaging and conformational sampling by molecular dynamics simulation, is presented. Polar and nonpolar contributions to the binding free energy are approximated by a linear scaling of the corresponding terms in the average intermolecular interaction energy for the bound and free states of the ligand. While the method originally assumed the validity of electrostatic linear response, we show that incorporation of systematic deviations from linear response derived from free energy perturbation calculations enhances the accuracy of the approach. The method is applied to complexes of wild-type and mutant human dihydrofolate reductases with 2,4-diaminopteridine and 2,4-diaminoquinazoline inhibitors. It is shown that a binding energy accuracy of about 1 kcal/mol is attainable even for multiply ionized compounds, such as methotrexate, for which electrostatic interactions energies are very large. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 77–88, 1998     

Solvent effect on the syn/anti conformational stability: A comparison between conformational bias Monte Carlo and molecular dynamics methods

The solvent effect on the syn/anti population ratio of the mesityl oxide (MOx) was investigated using a new implementation of conformational bias Monte Carlo (CBMC) and molecular dynamics (MD) methods. It was observed by a previous theoretical work (Theor. Chem. Acc. (2012) 131:1214) that in gas-phase the MOx exists dominantly in syn-form and in aqueous solution in anti-form. The syn/anti free energy difference in the gas phase was used in the intramolecular parametrization and a rotational barrier of approximately 10 kcal mol⁻¹ was found. Molecular systems with barriers of this order of magnitude have been studied by experimental techniques. However, they have not been discussed yet comparing CBMC and MD simulations. In this work, we show that the intramolecular geometrical information such as bond lengths, angles and torsional angles sampled with CBMC and MD methods are equivalent. Nonetheless, only the CBMC simulations sample appropriately the syn/anti population ratio. With the CBMC configurations in gas phase, it was obtained 95% in syn-form and 5% in anti-form regardless the initial conformation. An inversion of the population was found in water, 25% in syn-form and 75% in anti-form. Comparing the gas phase and in-water CBMC sampling, it was observed that the MOx spends typically approximately 110 successive MC cycles in anti-form and approximately 2300 in syn-form in gas phase. While it was much larger with explicit water, approximately 400 times more for anti-form and approximately 6 times more for syn-form. We argue that this strong stabilization of the anti-form in aqueous solution, does not come from the MOx-water hydrogen bonds interactions, because they are the same for both conformations. Instead, the stabilization comes from the dipole-dipole interaction caused by a larger dipole moment of the MOx in the anti-form, 7.2 D, than in the syn-form, 5.2 D. With the MD sampled configurations in both conditions, we observe that the syn/anti conformational change is a very rare event due to the rotational barrier, which is approximately 17 times larger than the thermal energy. Therefore, the MD sampling of the MOx is not appropriated because it is strongly dependent on the initial conformation even for large simulations with 150 ns up to 400 ns for the isolated solute and for solute–solvent systems.     

Identifying ligand binding sites and poses using GPU-accelerated Hamiltonian replica exchange molecular dynamics

We present a method to identify small molecule ligand binding sites and poses within a given protein crystal structure using GPU-accelerated Hamiltonian replica exchange molecular dynamics simulations. The Hamiltonians used vary from the physical end state of protein interacting with the ligand to an unphysical end state where the ligand does not interact with the protein. As replicas explore the space of Hamiltonians interpolating between these states, the ligand can rapidly escape local minima and explore potential binding sites. Geometric restraints keep the ligands from leaving the vicinity of the protein and an alchemical pathway designed to increase phase space overlap between intermediates ensures good mixing. Because of the rigorous statistical mechanical nature of the Hamiltonian exchange framework, we can also extract binding free energy estimates for all putative binding sites. We present results of this methodology applied to the T4 lysozyme L99A model system for three known ligands and one non-binder as a control, using an implicit solvent. We find that our methodology identifies known crystallographic binding sites consistently and accurately for the small number of ligands considered here and gives free energies consistent with experiment. We are also able to analyze the contribution of individual binding sites to the overall binding affinity. Our methodology points to near term potential applications in early-stage structure-guided drug discovery.     

Coarse‐grained molecular dynamics simulations of protein–ligand binding

Coarse‐grained molecular dynamics (CGMD) simulations with the MARTINI force field were performed to reproduce the protein–ligand binding processes. We chose two protein–ligand systems, the levansucrase–sugar (glucose or sucrose), and LinB–1,2‐dichloroethane systems, as target systems that differ in terms of the size and shape of the ligand‐binding pocket and the physicochemical properties of the pocket and the ligand. Spatial distributions of the Coarse‐grained (CG) ligand molecules revealed potential ligand‐binding sites on the protein surfaces other than the real ligand‐binding sites. The ligands bound most strongly to the real ligand‐binding sites. The binding and unbinding rate constants obtained from the CGMD simulation of the levansucrase–sucrose system were approximately 10 times greater than the experimental values; this is mainly due to faster diffusion of the CG ligand in the CG water model. We could obtain dissociation constants close to the experimental values for both systems. Analysis of the ligand fluxes demonstrated that the CG ligand molecules entered the ligand‐binding pockets through specific pathways. The ligands tended to move through grooves on the protein surface. Thus, the CGMD simulations produced reasonable results for the two different systems overall and are useful for studying the protein–ligand binding processes. © 2014 Wiley Periodicals, Inc.     

An investigation of the conformational dynamics of ion-pair intermediates in the syn-elimination of 8a-methyldecahydronaphthalen-4a-ols

The (4aβ,8aβ)-8a-methyl(5β D)decahydronaphthalen-4a-ol has been prepared and its reaction with H₂SO₄-Ac₂O-HOAc shown to occur with a k_h/k_D of 2.2 (±0.2) by loss of a proton (deuteron) syn to the departing oxy-anion. The results are used as a probe of carbocation-anion conformational mobility.     

Subtype selectivity and flexibility of ionotropic glutamate receptors upon antagonist ligand binding

The binding modes of a set of known ionotropic glutamate receptor antagonist-ligands have been studied using homology modeling, molecular docking, molecular dynamics (MD) simulations and ab initio quantum mechanical calculations. The core structure of the studied ligands is the decahydroisoquinoline ring, which has a carboxylic acid group at position three and different negatively-charged substituents (R) at position six. The binding affinities of these molecules have been reported earlier. From the current study, the carboxylate group of the decahydroisoquinoline ring hydrogen bonds with Arg485, the amino group with Pro478 and Thr480, and the negatively charged substituent R interacts with the positively charged N-terminus of helix-F. The subtype selectivity of these ligands seems to be strongly dependent on the amino acid at position 650 (GluR2: leucine, GluR5: valine), which affects the conformation of the ligand and ligand-receptor interactions, but depends considerably on the size of the R-group of the ligand. In addition, the MD simulations also revealed that the relative positions of the S1 and S2 domains can alter significantly showing different "closure" and "rotational movements" depending on the antagonist-ligand that is bound. Accordingly, molecular docking of antagonist ligands into static crystal structures cannot sufficiently explain ligand binding and subtype selectivity.     

Protein flexibility and species specificity in structure-based drug discovery: dihydrofolate reductase as a test system

In structure-based drug discovery, researchers would like to identify all possible scaffolds for a given target. However, techniques that push the boundaries of chemical space could lead to many false positives or inhibitors that lack specificity for the target. Is it possible to broadly identify the appropriate chemical space for the inhibitors and yet maintain target specificity? To address this question, we have turned to dihydrofolate reductase (DHFR), a well-studied metabolic enzyme of pharmacological relevance. We have extended our multiple protein structure (MPS) method for receptor-based pharmacophore models to use multiple X-ray crystallographic structures. Models were created for DHFR from human and Pneumocystis carinii. These models incorporate a fair degree of protein flexibility and are highly selective for known DHFR inhibitors over drug-like non-inhibitors. Despite sharing a highly conserved active site, the pharmacophore models reflect subtle differences between the human and P. carinii forms, which identify species-specific, high-affinity inhibitors. We also use structures of DHFR from Candida albicans as a counter example. The available crystal structures show little flexibility, and the resulting models give poorer performance in identifying species-specific inhibitors. Therapeutic success for this system may depend on achieving species specificity between the related human host and these key fungal targets. The MPS technique is a promising advance for structure-based drug discovery for DHFR and other proteins of biomedical interest.     

Role of aspartic acid in collagen structure and stability: A molecular dynamics investigation

A molecular dynamics (MD) simulation study has been carried out to understand the stability of the triple helical collagen models. The calculations show that the presence of the aspartic acid residue in different positions leads to the local variation in the structure. Analyses of root-mean-square deviation (RMSD), radial distribution function (RDF), puckering effect, dihedral angle variation, hydrogen bond (H-bond), and conformational changes during molecular dynamics simulation reveal that the local perturbation in the sequences, increase in chain flexibility due to removal of five membered rings in the collagen by aspartic acid, change of intermolecular H-bonding pattern, and differences in the association of water are mainly influencing the nature of stabilization of collagen by aspartic acid.     

Toward functional carboxylate-bridged diiron protein mimics: achieving structural stability and conformational flexibility using a macrocylic ligand framework

A dinucleating macrocycle, H(2)PIM, containing phenoxylimine metal-binding units has been prepared. Reaction of H(2)PIM with [Fe(2)(Mes)(4)] (Mes = 2,4,6-trimethylphenyl) and sterically hindered carboxylic acids, Ph(3)CCO(2)H or Ar(Tol)CO(2)H (2,6-bis(p-tolyl)benzoic acid), afforded complexes [Fe(2)(PIM)(Ph(3)CCO(2))(2)] (1) and [Fe(2)(PIM)(Ar(Tol)CO(2))(2)] (2), respectively. X-ray diffraction studies revealed that these diiron(II) complexes closely mimic the active site structures of the hydroxylase components of bacterial multicomponent monooxygenases (BMMs), particularly the syn disposition of the nitrogen donor atoms and the bridging μ-η(1)η(2) and μ-η(1)η(1) modes of the carboxylate ligands at the diiron(II) centers. Cyclic voltammograms of 1 and 2 displayed quasi-reversible redox couples at +16 and +108 mV vs ferrocene/ferrocenium, respectively. Treatment of 2 with silver perchlorate afforded a silver(I)/iron(III) heterodimetallic complex, [Fe(2)(μ-OH)(2)(ClO(4))(2)(PIM)(Ar(Tol)CO(2))Ag] (3), which was structurally and spectroscopically characterized. Complexes 1 and 2 both react rapidly with dioxygen. Oxygenation of 1 afforded a (μ-hydroxo)diiron(III) complex [Fe(2)(μ-OH)(PIM)(Ph(3)CCO(2))(3)] (4), a hexa(μ-hydroxo)tetrairon(III) complex [Fe(4)(μ-OH)(6)(PIM)(2)(Ph(3)CCO(2))(2)] (5), and an unidentified iron(III) species. Oxygenation of 2 exclusively formed di(carboxylato)diiron(III) compounds, a testimony to the role of the macrocylic ligand in preserving the dinuclear iron center under oxidizing conditions. X-ray crystallographic and (57)Fe M?ssbauer spectroscopic investigations indicated that 2 reacts with dioxygen to give a mixture of (μ-oxo)diiron(III) [Fe(2)(μ-O)(PIM)(Ar(Tol)CO(2))(2)] (6) and di(μ-hydroxo)diiron(III) [Fe(2)(μ-OH)(2)(PIM)(Ar(Tol)CO(2))(2)] (7) units in the same crystal lattice. Compounds 6 and 7 spontaneously convert to a tetrairon(III) complex, [Fe(4)(μ-OH)(6)(PIM)(2)(Ar(Tol)CO(2))(2)] (8), when treated with excess H(2)O.     

Effects of zinc binding on the conformational distribution of the amyloid-beta peptide based on molecular dynamics simulations

A number of experiments suggested that metal binding can promote the aggregation of the amyloid-beta peptide. In this work, the effects of the zinc binding on the conformational distributions of the full length amyloid-beta peptide are investigated on the basis of extensive molecular dynamics simulations. By comparing the conformational distributions of the apo-peptide and the holo-peptide, we show that the zinc binding can affect the conformational distribution of the amyloid-beta monomer dramatically. Compared with the apo-peptide, the holo-peptide samples more beta-strand conformation for the central hydrophobic cluster 17-21. Meanwhile, the formation probabilities of the salt bridge Asp23-Lys28 and the turn comprising 23-28 are also increased significantly at room temperature. Since these local structures are essential for the amyloid-beta aggregation, the observed effects of the zinc binding indicate that the metal induced conformational change of the monomer is one of the possible mechanisms for the metal promoted aggregation of the amyloid-beta peptide.