Molecular Dynamics Free Energy Simulations Research Articles

Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations

Riboswitches can regulate gene expression by direct and specific interactions with ligands and have recently attracted interest as potential drug targets for antibacterial. In this work, molecular dynamics (MD) simulations, free energy perturbation (FEP) and molecular mechanics generalized Born surface area (MM-GBSA) methods were integrated to probe the effect of mutations on the binding of ligands to guanine riboswitch (GR). The results not only show that binding free energies predicted by FEP and MM-GBSA obtain an excellent correlation, but also indicate that mutations involved in the current study can strengthen the binding affinity of ligands GR. Residue-based free energy decomposition was applied to compute ligand-nucleotide interactions and the results suggest that mutations highly affect interactions of ligands with key nucleotides U22, U51 and C74. Dynamics analyses based on MD trajectories indicate that mutations not only regulate the structural flexibility but also change the internal motion modes of GR, especially for the structures J12, J23 and J31, which implies that the aptamer domain activity of GR is extremely plastic and thus readily tunable by nucleotide mutations. This study is expected to provide useful molecular basis and dynamics information for the understanding of the function of GR and possibility as potential drug targets for antibacterial.

Synthesis, molecular docking, binding free energy calculation and molecular dynamics simulation studies of benzothiazol-2-ylcarbamodithioates as Staphylococcus aureus MurD inhibitors

A new series of benzothiazol-2-ylcarbamodithioate functional compounds 5a-f has been designed, synthesized and characterized by spectral data. These compounds were screened for their in vitro antibacterial activity against strains of Staphylococcus aureus (NCIM 5021, NCIM 5022 and methicillin-resistant isolate 43300), Bacillus subtilis (NCIM 2545), Escherichia coli (NCIM 2567), Klebsiella pneumoniae (NCIM 2706) and Psudomonas aeruginosa (NCIM 2036). Compounds 5a and 5d exhibited significant activity against all the tested bacterial strains. Specifically, compounds 5a and 5d showed potent activity against K. pneumoniae (NCIM 2706), while compound 5a also displayed potent activity against S. aureus (NCIM 5021). Compound 5d showed minimum IC50 value of 13.37 μM against S. aureus MurD enzyme. Further, the binding interactions of compounds 5a-f in the catalytic pocket have been investigated using the extra-precision molecular docking and binding free energy calculation by MM-GBSA approach. A 30 ns molecular dynamics simulation of 5d/modeled S. aureus MurD enzyme was performed to determine the stability of the predicted binding conformation.

Phosphorene: A promising metal free cathode material for proton exchange membrane fuel cell

The design of new cathode materials is vital to develop the high performance proton exchange membrane fuel cell (PEMFC). The adsorption properties of the species involved in the oxygen reduction reaction (ORR) and several ORR mechanisms on phosphorene were presented with using the first-principle calculations. Results show that the H adsorption is more preferable than O2, and the OOH spontaneously form from the adsorbed H and O2. The HOOH from the hydrogenation of the OOH species could be further hydrogenated into OH and H2O, and the OH would be finally hydrogenated into the (second) H2O. The reaction barrier of rate-determining step is 1.02 eV. Moreover, the presented kinetic processes are confirmed by the thermodynamic simulation from Gibbs free energy calculations and molecular dynamics simulations. These results indicate that phosphorene as one of new 2D materials could be a promising metal free cathode material for PEMFC.

Thermodynamics of the NMDA Receptor Amino-Terminal Domain

The N-methyl-D-aspartate receptor (NMDAR) is a ligand-gated ion channel found in post-synaptic neurons that is critically important in learning and memory. This tetrameric membrane protein has two extracellular domains, a ligand-binding domain (LBD) and an amino-terminal domain (ATD), that control the movement of ions through the channel. Specifically, conformational changes in the ATD layer have been linked to the allosteric regulation of channel gating. Previous structural studies of the GluN1/GluN2B NMDAR have identified a ∼20° opening of the GluN2B ATD and a ∼15° rotation of the GluN1 and GluN2B ATD subunits with respect to the dimer interface in the presence of allosteric modulators. In this investigation, we constructed models of both GluN1 and GluN2B ATD monomers and the GluN1/GluN2B dimer with and without allosteric modulators. We then selected order parameters that capture the range of conformational motion adopted by these domains. By performing free energy molecular dynamics simulations, we gain insight into the thermodynamics that drive these observed structural changes.

Coarse-grained free-energy simulations of conformational state transitions in an adenosine 5′-triphosphate-binding cassette exporter

ATP-binding cassette exporters transport many substrates out of cellular membranes via alternating between inward-facing and outward-facing conformations. Despite extensive research efforts over the past decades, understanding of the molecular mechanism remains elusive. As these large-scale conformational movements are global and collective, we have previously performed extensive coarse-grained molecular dynamics simulations of the potential of mean force along the conformational transition pathway [J. Phys. Chem. B 119, 1295 (2015)]. However, the occluded conformational state, in which both the internal and external gate are closed, was not determined in the calculated free energy profile. In this work, we extend the above methods to the calculation of the free energy profile along the reaction coordinate, d1–d2, which are the COM distances between the two sides of the internal (d1) and the external gate (d2). The potential of mean force is thus obtained to identify the transition pathway, along which several outward-facing, inward-facing, and occluded state structures are predicted in good agreement with structural experiments. Our coarse-grained molecular dynamics free-energy simulations demonstrate that the internal gate is closed before the external gate is open during the inward-facing to outward-facing transition and vice versa during the inward-facing to outward-facing transition. Our results capture the unidirectional feature of substrate translocation via the exporter, which is functionally important in biology. This finding is different from the previous result, in which both the internal and external gates are open reported in an X-ray experiment [Proc. Natl. Acad. Sci. USA 104, 19005 (2007)]. Our study sheds light on the molecular mechanism of the state transitions in an ATP-binding cassette exporter.

How methyl-sugar interactions determine DNA structure and flexibility.

The sequence dependent structure and flexibility of the DNA double helix is of key importance for gene expression and DNA packing and it can be modulated by DNA modifications. The presence of a C5′-methyl group in thymine or the frequent C5′-methylated-cytosine affects the DNA fine structure, however, the underlying mechanism and steric origins have remained largely unexplained. Employing Molecular Dynamics free energy simulations that allow switching on or off interactions with the methyl groups in several DNA sequences, we systematically identified the physical origin of the coupling between methyl groups and DNA backbone fine structure. Whereas methyl-solvent and methyl–nucleobase interactions were found to be of minor importance, the methyl group interaction with the 5′ neighboring sugar was identified as main cause for influencing the population of backbone substates. The sterical methyl sugar clash prevents the formation of unconventional stabilizing hydrogen bonds between nucleobase and backbone. The technique was also used to study the contribution of methyl groups to DNA flexibility and served to explain why the presence of methyl sugar clashes in thymine and methyl-cytosine can result in an overall local increase of DNA flexibility.

Neutron and X-ray crystal structures of Lactobacillus brevis alcohol dehydrogenase reveal new insights into hydrogen-bonding pathways.

Lactobacillus brevis alcohol dehydrogenase (LbADH) is a well studied homotetrameric enzyme which catalyzes the enantioselective reduction of prochiral ketones to the corresponding secondary alcohols. LbADH is stable and enzymatically active at elevated temperatures and accepts a broad range of substrates, making it a valuable tool in industrial biocatalysis. Here, the expression, purification and crystallization of LbADH to generate large, single crystals with a volume of up to 1 mm3 suitable for neutron diffraction studies are described. Neutron diffraction data were collected from an H/D-exchanged LbADH crystal using the BIODIFF instrument at the Heinz Maier-Leibnitz Zentrum (MLZ), Garching, Germany to a resolution dmin of 2.15 Å in 16 days. This allowed the first neutron crystal structure of LbADH to be determined. The neutron structure revealed new details of the hydrogen-bonding network originating from the ion-binding site of LbADH and provided new insights into the reasons why divalent magnesium (Mg2+) or manganese (Mn2+) ions are necessary for its activity. X-ray diffraction data were obtained from the same crystal at the European Synchrotron Radiation Facility (ESRF), Grenoble, France to a resolution dmin of 1.48 Å. The high-resolution X-ray structure suggested partial occupancy of Mn2+ and Mg2+ at the ion-binding site. This is supported by the different binding affinity of Mn2+ and Mg2+ to the tetrameric structure calculated via free-energy molecular-dynamics simulations.

Elucidating Aggregation Pathways in the Donor-Acceptor Type Molecules p-DTS(FBTTh2)2 and p-SIDT(FBTTh2)2.

We investigated the aggregation behavior of the donor-acceptor molecules p-DTS(FBTTh2)2 (T1) and p-SIDT(FBTTh2)2 (H1) in MTHF solutions. Using optical spectroscopy, we found that T1 forms aggregates in solution while H1 aggregates only when processed as a thin film, but not in solution. Free energy molecular dynamics (MD) simulations based on force fields derived from quantum-mechanical density functional theory fully reproduce this difference. Our simulations reveal that this difference is not due to the lengthy carbon side chains. Rather, the molecular symmetry of T1 allows for an aggregated state in which the central donor units are spatially well-separated while a similar configuration is sterically impossible for H1. As a consequence, any aggregation of H1 necessarily involves aggregation of the central donors which requires, as a first step, stripping the central donor of its protective MTHF solvation shell. This unfavorable process leads to a significant kinetic hindrance for aggregation and explains the strongly differing aggregation behavior of T1/H1 in MTHF despite their otherwise similar structures.

Identification of PfENR inhibitors: A hybrid structure-based approach in conjunction with molecular dynamics simulations.

In the current study, we have constructed receptor-based pharmacophore models by exploiting the Plasmodium falciparum enoyl-acyl carrier protein reductase (PfENR) structural proteome. The derived models were subjected to a series of validation procedures to list the representative hypotheses that can be used for the screening of the Drug-like Diverse Database. A set of 739 molecules was retrieved and analyzed for the adsorption, distribution, metabolism, excretion and toxicity (ADMET) and drug-likeness attributes. The filtered drug-like molecules (64) were then subjected to molecular docking and HYDE assessment studies. The hybrid structure-based approach yielded 4 molecules, UKR1308259, ENA1096786, UKR403454, and ASI51224, as PfENR inhibitors. The stability of these inhibitors was assessed using molecular mechanics-generalized born surface area approach-based free binding energy calculations and molecular dynamics simulations. Molecular mechanics-generalized born surface area calculations and molecular dynamics simulations showed that UKR1308259, ENA1096786, and ASI51224 were more potent PfENR inhibitors. The rationale behind the current work was to identify orally available inhibitor molecules with diverse scaffolds that could serve as initial leads for the drug design against PfENR.

Spectroscopic and Molecular Dynamics Simulation Study of Lysozyme in the Aqueous Mixture of Ethanol: Insights into the Nonmonotonic Change of the Structure of Lysozyme.

Structural change of lysozyme in the aqueous mixture of ethanol (ETH) has been investigated using emission, circular dichroism spectroscopy, free-energy molecular dynamics (MD) simulation, and contact map analysis methods. The emission and circular dichroism data of protein depict the nonmonotonic change suggesting that the structure as well as local environment near the Trp of lysozyme modifies differently for different compositions of the ETH-water mixture. The free-energy MD simulation shows that the change in the average size of lysozyme in the aqueous mixture of ETH also shows a nonmonotonic nature. The free-energy profile of lysozyme in the low concentration of ETH suggests that the distance between helices increases (χETH ≈ 0.07) first and decreases again (χETH ≈ 0.2) becoming almost similar to the native structure. Around χETH ≈ 0.5, the size of lysozyme increases significantly probably leading to the unfolding of the protein. With further increase of ETH (χETH ≈ 0.7), size of lysozyme decreases suggesting the refolding of the protein. Contact map as well as solvent organization analysis depicts that ETH gets preferentially solvated by the hydrophobic core of lysozyme which weaken the hydrophobic interactions of protein, resulting in the increase in its size. The aggregation of ETH dominated at the higher concentration of ETH (χETH ≈ 0.7); hence, the, weakening of hydrophobic core of protein by ETH is reduced, which probably lead to the refolding of lysozyme.

Classical Drude Polarizable Force Field Model for Methyl Phosphate and Its Interactions with Mg2.

Phosphate groups are essential components of nucleic acids and proteins, whose interactions with solvent, metal ions, and ionic side chains help control folding and binding. Methyl phosphate (MP) represents a simple analog of phosphate moieties that are post-translation modifications in proteins and present at the termini of nucleic acids, among other environments. In the present study, we optimized parameters for use in polarizable molecular dynamics simulations of MP in its mono- and dianionic forms, MP- ≡ CH3HPO4- and MP2- ≡ CH3PO42-, along with P i2- ≡ HPO42-, in the context of the classical Drude oscillator model. Parameter optimization was done in a manner consistent with the remainder of the Drude molecular mechanics force field, choosing atomic charges and polarizabilities to reproduce molecular properties from quantum mechanics as well as experimental hydration free energies. Optimized parameters were similar to existing dimethyl phosphate parameters, with a few significant differences. The developed parameters were then used to compute magnesium binding affinities in aqueous solution, using alchemical molecular dynamics free energy simulations. Good agreement with experiment was obtained, and outer sphere binding was shown to be predominant for MP- and MP2-.

The potential impact of carboxylic-functionalized multi-walled carbon nanotubes on trypsin: A Comprehensive spectroscopic and molecular dynamics simulation study.

In this study, we report a detailed experimental, binding free energy calculation and molecular dynamics (MD) simulation investigation of the interactions of carboxylic-functionalized multi-walled carbon nanotubes (COOH-f-MWCNTs) with porcine trypsin (pTry). The enzyme exhibits decreased thermostability at 330K in the presence of COOH-f-MWCNTs. Furthermore, the activity of pTry also decreases in the presence of COOH-f-MWCNTs. The restricted diffusion of the substrate to the active site of the enzyme was observed in the experiment. The MD simulation analysis suggested that this could be because of the blocking of the S1 pocket of pTry, which plays a vital role in the substrate selectivity. The intrinsic fluorescence of pTry is quenched with increase in the COOH-f-MWCNTs concentration. Circular dichroism (CD) and UV–visible absorption spectroscopies indicate the ability of COOH-f-MWCNTs to experience conformational change in the native structure of the enzyme. The binding free energy calculations also show that electrostatics, π-cation, and π-π stacking interactions play important roles in the binding of the carboxylated CNTs with pTry. The MD simulation results demonstrated that the carboxylated CNTs adsorb to the enzyme stronger than the CNT without the–COOH groups. Our observations can provide an example of the nanoscale toxicity of COOH-f-MWCNTs for proteins, which is a critical issue for in vivo application of COOH-f-MWCNTs.

Principle of K+/Na+ selectivity in the active site of group II intron at various stages of self-splicing pathway

Group II introns are ribozymes which can catalyze its own splicing and relegate itself. They share common structural features and are evolutionarily related to eukaryotic spliceosome. Hence, group II introns are an excellent model system for understanding the mechanism of RNA-splicing in gene expression. Recent advancement of structural studies has provided x-ray structures of group II intron at different stages (Pre-hydrolytic, post-hydrolytic and free intron) of splicing pathway and revealed heteronuclear metal ion cluster (two potassium; K1 and K2, two magnesium; M1 and M2) in the active site as common structural feature of group II introns. It is believed that these four metal ions are very crucial for catalysis. K+ is very important for the function and buffers containing only Na+ can destroy the function of group II introns. Despite the available 3D structures of different stages of self-splicing, the energetic origin of K+ selectivity over Na+ is not known. However, these structures provide sufficiently good models for directly computing the energetics of K+ vs Na+ discrimination that is required for understanding self-splicing. Here, we report extensive (∼550 ns) of structure-based molecular dynamics free energy simulations that quantitatively estimated K1 vs Na1 discrimination. We show that the strength of discrimination varies along the self-splicing pathway. The free and pre-hydrolytic states are highly selective for K1 and the selectivity is lost in the post-hydrolytic state. The simulations further reveal that Na1 in the free and pre-hydrolytic state is trapped with an unsatisfied first coordination shell in the active site, which is responsible for large discrimination. The Na1 in the active site of post-hydrolytic state allow water entry and satisfy the Na1 bonding requirement, resulting in low discrimination. The results not only give insights into the experimentally unresolved Na1 bound complexes but also provide the link between structures and computed energetics.

Absolute proton hydration free energy, surface potential of water, and redox potential of the hydrogen electrode from first principles: QM/MM MD free-energy simulations of sodium and potassium hydration.

The absolute intrinsic hydration free energy GH+,wat◦ of the proton, the surface electric potential jump χwat◦ upon entering bulk water, and the absolute redox potential VH+,wat◦ of the reference hydrogen electrode are cornerstone quantities for formulating single-ion thermodynamics on absolute scales. They can be easily calculated from each other but remain fundamentally elusive, i.e., they cannot be determined experimentally without invoking some extra-thermodynamic assumption (ETA). The Born model provides a natural framework to formulate such an assumption (Born ETA), as it automatically factors out the contribution of crossing the water surface from the hydration free energy. However, this model describes the short-range solvation inaccurately and relies on the choice of arbitrary ion-size parameters. In the present study, both shortcomings are alleviated by performing first-principle calculations of the hydration free energies of the sodium (Na+) and potassium (K+) ions. The calculations rely on thermodynamic integration based on quantum-mechanical molecular-mechanical (QM/MM) molecular dynamics (MD) simulations involving the ion and 2000 water molecules. The ion and its first hydration shell are described using a correlated ab initio method, namely resolution-of-identity second-order Møller-Plesset perturbation (RIMP2). The next hydration shells are described using the extended simple point charge water model (SPC/E). The hydration free energy is first calculated at the MM level and subsequently increased by a quantization term accounting for the transformation to a QM/MM description. It is also corrected for finite-size, approximate-electrostatics, and potential-summation errors, as well as standard-state definition. These computationally intensive simulations provide accurate first-principle estimates for GH+,wat◦, χwat◦, and VH+,wat◦, reported with statistical errors based on a confidence interval of 99%. The values obtained from the independent Na+ and K+ simulations are in excellent agreement. In particular, the difference between the two hydration free energies, which is not an elusive quantity, is 73.9 ± 5.4 kJ mol-1 (K+ minus Na+), to be compared with the experimental value of 71.7 ± 2.8 kJ mol-1. The calculated values of GH+,wat◦, χwat◦, and VH+,wat◦ (-1096.7 ± 6.1 kJ mol-1, 0.10 ± 0.10 V, and 4.32 ± 0.06 V, respectively, averaging over the two ions) are also in remarkable agreement with the values recommended by Reif and Hünenberger based on a thorough analysis of the experimental literature (-1100 ± 5 kJ mol-1, 0.13 ± 0.10 V, and 4.28 ± 0.13 V, respectively). The QM/MM MD simulations are also shown to provide an accurate description of the hydration structure, dynamics, and energetics.

In silico quest of selective naphthyl-based CREBBP bromodomain inhibitor.

The reader proteins like bromodomains have recently gained increased attentions in the area of epigenetic drug discovery, as they are the potent regulators in gene transcription process. Among the other bromodomains, cAMP response element-binding protein (CREB) binding protein or CREBBP bomodomain involved in various cancer progressions and therefore, efforts to develop specific inhibitors of CREBBP bomodomain are of clinical value. In this study, we tried to identify selective CREBBP bromodomain inhibitor, which was accomplished by using molecular docking, free energy calculation and molecular dynamics (MD) simulation studies, considering a series of naphthyl based compounds. The docking procedure was validated by comparing root mean square deviations (RMSDs) of crystallographic complex to docked complex. Favorable electrostatic interactions with the Arg1173 side chain were considered to attain selectivity for CREBBP bromodomain over other human bromodomain subfamilies. We found that naphthyl-based compounds have greater binding affinities towards the CREBBP bromodomain, and formed non-bonded interactions with various side chain residues that are important for bromodomain inhibition. From detailed investigation by induced fit docking, compound 31 was found to have favorable electrostatic interactions with the Arg1173 side chain by forming conventional hydrogen bonds. This result was further confirmed by analyzing hydrogen bond occupancy and bonding distance during the molecular dynamics simulation. We believe that these findings offer useful insight for the designing of target specific new bromodomain inhibitor and also promote further structure guided synthesis of analogues for identification of potent CREBBP bromodomain inhibitors as well as detailed in vitro and in vivo analyses.

Effect of phosphorylation and single nucleotide polymorphisms on caspase substrates processing.

Posttranslational modifications that involve either reversible covalent modification of proteins or irreversible proteolysis are central to the regulation of key cellular mechanisms, including apoptosis, cell-cycle regulation and signal transduction. There is mounting evidence suggesting cross-talk between proteases and kinases. For instance: caspases, a class of proteases involved in programmed cell death-apoptosis, cleave a large set of various types of proteins. Simultaneously, kinases restrict caspase activity by phosphorylating their protein substrates in the vicinity of cleavage site. In addition, the caspase cleavage pattern in target proteins may be modified as a result of single nucleotide polymorphisms (SNPs) in the coding gene. This may either create a novel cleavage site, or increase/decrease the cleavage efficiency of a substrate. Such point mutations are often associated with the onset of disease. In this study, we predicted how phosphorylation and SNPs affect known human caspase proteolytic events collected in the CASBAH and Degrabase databases by applying Random Forest caspases' substrates prediction method, as implemented in the CaspDB, and the molecular dynamics free energy simulations approach. Our analysis confirms several experimental observations. Phosphorylation could have both positive or negative regulatory effects depending on its position with respect to the caspase cleavage site. For instance, we demonstrate that phosphorylation at P1' is the most detrimental for proteolytic efficiency of caspases. Phosphorylation at the P2 and P2' positions also negatively affect the cleavage events. In addition, we uncovered SNPs in 11 caspase substrates capable of completely abolishing the cleavage site due to polymorphism at the P1 position. The findings presented here may be useful for determining the link between aberrant proteolysis and disease.

An MM and QM Study of Biomimetic Catalysis of Diels-Alder Reactions Using Cyclodextrins.

We performed a computational investigation of the mechanism by which cyclodextrins (CDs) catalyze Diels-Alder reactions between 9-anthracenemethanol and N-cyclohexylmaleimide. Hydrogen bonds (Hbonds) between N-cyclohexylmaleimide and the hydroxyl groups of cyclodextrins were suggested to play an important role in this catalytic process. However, our free energy calculations and molecular dynamics simulations showed that these Hbonds are not stable, and quantum mechanical calculations suggested that the reaction is not promoted by these Hbonds. The binding of 9-anthracenemethanol and N-cyclohexylmaleimide to cyclodextrins was the key to the catalytic process. Cyclodextrins act as a container to hold the two reactants in the cavity, pre-organize them for the reactions, and thus reduce the entropy penalty to the activation free energy. Dimethyl-β-CD was a better catalyst for this specific reaction than β-CD because of its stronger van der Waals interaction with the pre-organized reactants and its better performance in reducing the activation energy. This computational work sheds light on the mechanism of the catalytic reaction by cyclodextrins and introduces new perspectives of supramolecular catalysis.

In-silico interaction studies suggest RND efflux pump mediates polymyxin resistance in Acinetobacter baumannii

Bacterial efflux pumps have emerged as antibiotic resistance determinants and confers multi-drug resistance to a broad range of antimicrobials as well as non-antibiotic substances. A study about translocation of antibiotic molecules through the efflux transporter, will contribute in determining substrate specificity. In the present study, we have explored RND family efflux pump extensively found in Acinetobacter baumannii i.e. AdeABC. Besides, another well studied RND efflux pump, AcrAB-TolC together with a non-RND efflux pump, NorM was investigated for comparative analysis. We employed a series of computational techniques ranging from molecular docking to binding free energy estimation and molecular dynamics simulations to determine the binding affinity for different classes of drugs, namely aminoglycosides, polymyxins, β-lactams, tetracyclines, glycylcyclines, quinolones and metronidazole with AdeB, AcrB, and NorM efflux proteins. Our results revealed that class polymyxins has the highest binding affinity with the RND efflux pumps i.e. AcrAB-TolC and AdeABC as well as non-RND efflux pump, NorM. The experimental validation study demonstrated bigger zone of inhibition in presence of efflux pump inhibitor than polymyxin alone thus unveiling its specificity toward efflux pump. The reported experimental data comprising of minimum inhibitory concentration of antibiotics toward these efflux pumps also support our finding based on in silico approach. To recapitulate the outcome, polymyxins shows maximum specificity toward RND as well as non-RND efflux pump and may unlatch the way to rationally develop new potential antibacterial agents as well as efflux pump inhibitors in order to combat resistance.

Unwinding Induced Melting of Double-Stranded DNA Studied by Free Energy Simulations.

DNA unwinding plays a major role in many biological processes, such as replication, transcription, and repair. It can lead to local melting and strand separation and can serve as a key mechanism to promote access to the separate strands of a double-stranded DNA. While DNA unwinding has been investigated extensively by DNA cyclization and single-molecule studies on a length-scale of kilo base pairs, it is neither fully understood at the base pair level nor at the level of molecular interactions. By employing a torque acting on the termini of DNA oligonucleotides during molecular dynamics free energy simulations, we locally unwind the central part of a DNA beyond an elastic (harmonic) regime. The simulations reproduce experimental results on the twist elasticity in the harmonic regime (characterized by a mostly quadratic free energy change with respect to changes in twist) and a deformation up to 7° was found as a limit of the harmonic response. Beyond this limit the free energy increase per twist change dropped dramatically coupled to local base pair disruptions and significant deformation of the nucleic acid backbone structure. Restriction of the DNA bending flexibility resulted in a stiffer harmonic response and an earlier onset of the anharmonic response. Whereas local melting with a complete disruption of base pairing and flipping of nucleotides was observed in case of an AT rich central segment strong backbone changes and changes in the stacking arrangements were observed in case of a GC rich segment. Unrestrained MD simulations starting from locally melted DNA reformed regular B-DNA after 50-300 ns simulation time. The simulations may have important implications for understanding DNA recognition processes coupled with significant structural alterations.

Induced fit docking, free energy calculation and molecular dynamics studies on Mycobacterium tuberculosis alanine racemase inhibitor

Alar, a Pyridoxal 5′-phosphate (PLP)-dependent bacterial enzyme is responsible for the racemisation of L-alanine into D-alanine which is essential for the peptidoglycan biosynthesis in both Gram-positive and Gram-negative bacteria. In the present study, we performed induced fit docking, binding free energy calculation and molecular dynamics simulation to elucidate the Alar inhibition potential of 1,2,4-thiadiazolidine-3,5-dione-based inhibitor 1. The inhibitor binds to the hydrophobic groove of Alar and the binding was found to be stable throughout 20-ns MD simulation. Induced fit docking result showed that Lys42, Tyr46, Tyr175 and Tyr364 residues are primarily responsible for the stabilisation of inhibitor–protein complex. Further, high negative van der Waals binding free energy value of –38.88 kcal/mol, indicated it as the main driving force for the inhibitor binding. Based on the information obtained from this study, we designed few molecules as potent Alar inhibitor. In order to gain structural insight and to validate the stability of complex, we performed 20-ns MD simulation of the designed molecule D1. Results obtained from this study can be used for the design of M. tuberculosis Alar potent inhibitors lacking affinity for the co-factor PLP.