MMPBSA.py: An Efficient Program for End-State Free Energy Calculations

Rapid, Precise, and Reproducible Prediction of Peptide-MHC Binding Affinities from Molecular Dynamics That Correlate Well with Experiment

The presentation of potentially pathogenic peptides by major histocompatibility complex (MHC) molecules is one of the most important processes in adaptive immune defense. Prediction of peptide-MHC (pMHC) binding affinities is therefore a principal objective of theoretical immunology. Machine learning techniques achieve good results if substantial experimental training data are available. Approaches based on structural information become necessary if sufficiently similar training data are unavailable for a specific MHC allele, although they have often been deemed to lack accuracy. In this study, we use a free energy method to rank the binding affinities of 12 diverse peptides bound by a class I MHC molecule HLA-A*02:01. The method is based on enhanced sampling of molecular dynamics calculations in combination with a continuum solvent approximation and includes estimates of the configurational entropy based on either a one or a three trajectory protocol. It produces precise and reproducible free energy estimates which correlate well with experimental measurements. If the results are combined with an amino acid hydrophobicity scale, then an extremely good ranking of peptide binding affinities emerges. Our approach is rapid, robust, and applicable to a wide range of ligand-receptor interactions without further adjustment.

Stability of Iowa mutant and wild type Aβ-peptide aggregates

Recent experiments indicate a connection between the structure of amyloid aggregates and their cytotoxicity as related to neurodegenerative diseases. Of particular interest is the Iowa Mutant, which causes early-onset of Alzheimer's disease. While wild-type Amyloid β-peptides form only parallel beta-sheet aggregates, the mutant also forms meta-stable antiparallel beta sheets. Since these structural variations may cause the difference in the pathological effects of the two Aβ-peptides, we have studied in silico the relative stability of the wild type and Iowa mutant in both parallel and antiparallel forms. We compare regular molecular dynamics simulations with such where the viscosity of the samples is reduced, which, we show, leads to higher sampling efficiency. By analyzing and comparing these four sets of all-atom molecular dynamics simulations, we probe the role of the various factors that could lead to the structural differences. Our analysis indicates that the parallel forms of both wild type and Iowa mutant aggregates are stable, while the antiparallel aggregates are meta-stable for the Iowa mutant and not stable for the wild type. The differences result from the direct alignment of hydrophobic interactions in the in-register parallel oligomers, making them more stable than the antiparallel aggregates. The slightly higher thermodynamic stability of the Iowa mutant fibril-like oligomers in its parallel organization over that in antiparallel form is supported by previous experimental measurements showing slow inter-conversion of antiparallel aggregates into parallel ones. Knowledge of the mechanism that selects between parallel and antiparallel conformations and determines their relative stability may open new avenues for the development of therapies targeting familial forms of early-onset Alzheimer's disease.

A 3D-RISM/RISM study of the oseltamivir binding efficiency with the wild-type and resistance-associated mutant forms of the viral influenza B neuraminidase

The binding affinity of oseltamivir to the influenza B neuraminidase and to its variants with three single substitutions, E119G, R152K, and D198N, is investigated by the MM/3D-RISM method. The binding affinity or the binding free energy of ligand to receptor was found to be determined by a subtle balance of two major contributions that largely cancel out each other: the ligand-receptor interactions and the dehydration free energy. The theoretical results of the binding affinity of the drug to the mutants reproduced the observed trend in the resistivity, measured by IC50 ; the high-level resistance of E119G and R152K, and the low-level resistance of D198N. For E119G and R152K, reduction of the direct drug-target interaction, especially at the mutated residue, is the main source of high-level oseltamivir resistance. This phenomenon, however, is not found in the D198N strain, which is located in the framework of the active-site.

Binding free energy based structural dynamics analysis of HIV-1 RT RNase H-inhibitor complexes

Accurate prediction of binding free energies associated with small molecules binding to a receptor is a major challenge in drug design processes. To achieve this goal many computational methods have been developed ranging from highly efficient empirical based docking schemes to high accuracy methods based on e.g. free energy calculations. In this study, binding affinity predictions for a set of HIV-1 RNase H inhibitors have been performed using MM-PB(GB)/SA methods. The current study describes in detail how the choice of initial ligand structures, e.g. protonation states, impacts the predicted ranking of the compounds. In addition we study the structural dynamics of the RNase H complexes using molecular dynamics. The role of each residue contribution to the overall binding free energy is also explored and used to explain the variations in the inhibition potency. The results reported here can be useful for design of small molecules against RNase H activity in the development of effective drugs for HIV treatment.

An Investigation on the Fundamental Interaction between Abeta Peptides and the AT-Rich DNA

DNA damage is ubiquitous in all mammalian cells with the occurrence of more than 60,000 times per day per cell. In particular, DNA damage in neurons is found to accumulate with age and has been suggested to interfere with the synthesis of functional proteins. Moreover, recent studies have found through transgenic mice that human amyloid precursor protein causes an increase in DNA double-strand breaks (DSBs) with the effect of a prolongation in DNA repair. It is surmised that amyloid β (Aβ) exacerbates the DNA DSBs in neurons, possibly engendering neuronal dysfunction as a result. However, a good understanding on the holistic interaction mechanisms and the manner in which Aβ intertwines with DNA damage is still in its infancy. In our study, we found that DNA with an AT-rich sequence has a very low binding affinity toward Aβ by means of molecular dynamics simulation. While we have pursued a particular sequence of DNA in this study, other DNA sequences are expected to affect the interaction and binding affinity between DNA and Aβ, and will be pursued in our further research. Nonetheless, we have uncovered favorable interaction between the positively charged side chain of Aβ and the two ends of DNA. The latest experiment reveals that many of the double-stranded breaks in neurons can be fixed via DNA repair mechanisms but not in the case that Aβs are present. It is found that the increased numbers of DSBs prevail in active neurons. Here, on the basis of the favorable interaction between Aβ and the two ends of DNA, we propose the possibility that Aβ prevents DNA repair via binding directly to the break ends of the DNA, which further exacerbates DNA damage. Moreover, we have found that the base pair oxygen of the DNA has a greater preference to form hydrogen bonds than the backbone oxygen with Aβ at the two ends. Thus, we postulate that Aβ could serve to prevent the repair of AT-rich DNA, and it is unlikely to cause its breakage or affect its binding toward histone. Another important observation from our study is that AT-rich DNA has very little or no influence on Aβ oligomerization. Finally, even though we do not observe any dramatic DNA conformational change in the presence of Aβ, we do observe an increase in diversity of the DNA structural parameters such as groove width, local base step, and torsional angles in lieu of Aβ interactions.

Inhibitory effects and oxidation of 6-methylcoumarin, 7-methylcoumarin and 7-formylcoumarin via human CYP2A6 and its mouse and pig orthologous enzymes

1. Information about the metabolism of compounds is essential in drug discovery and development, risk assessment of chemicals and further development of predictive methods. 2. In vitro and in silico methods were applied to evaluate the metabolic and inhibitory properties of 6-methylcoumarin, 7-methylcoumarin and 7-formylcoumarin with human CYP2A6, mouse CYP2A5 and pig CYP2A19. 3. 6-Methylcoumarin was oxidized to fluorescent 7-hydroxy-6-methylcoumarin by CYP2A6 (Km: 0.64-0.91 µM; Vmax: 0.81-0.89 min(-1)) and by CYP2A5 and CYP2A19. The reaction was almost completely inhibited at 10 µM 7-methylcoumarin in liver microsomes of human and mouse, but in pig only 40% inhibition was obtained with the anti-CYP2A5 antibody or with methoxsalen and pilocarpine. 7-Methylcoumarin was a mechanism-based inhibitor for CYP2A6, but not for the mouse and pig enzymes. 7-Formylcoumarin was a mechanism-based inhibitor for CYP2As of all species. 4. Docking and molecular dynamics simulations of 6-methylcoumarin and 7-methylcoumarin in the active sites of CYP2A6 and CYP2A5 demonstrated a favorable orientation of the 7-position of 6-methylcoumarin towards the heme moiety. Several orientations of 7-methylcoumarin were possible in CYP2A6 and CYP2A5. 5. These results indicate that the active site of CYP2A6 has unique interaction properties for ligands and differs in this respect from CYP2A5 and CYP2A19.

Mechanistic studies for tri-targeted inhibition of enzymes involved in cholesterol biosynthesis by green tea polyphenols

In the present study, we found that three enzymes, MVK, MDD and FPPS, in the mevalonate pathway (MVP) of cholesterol biosynthesis, can be simultaneously inhibited by two green tea polyphenols ((-)-epicatechin-3-gallate, ECG; (-)-epigallocatechin-3-gallate, EGCG). Molecular dynamics simulations and pharmacophore studies were carried out to elucidate the tri-targeted inhibition mechanisms. Our results indicate that similar triangular binding pockets exist in all three enzymes, which is essential for their binding with polyphenols. Two distinct binding poses for ECG and EGCG were observed in our MD simulations. These results shed light on the potential for further selective and multi-targeted inhibitor design for the treatment of hyperlipidemia.

Elucidation of Ligand-Dependent Modulation of Disorder-Order Transitions in the Oncoprotein MDM2

Numerous biomolecular interactions involve unstructured protein regions, but how to exploit such interactions to enhance the affinity of a lead molecule in the context of rational drug design remains uncertain. Here clarification was sought for cases where interactions of different ligands with the same disordered protein region yield qualitatively different results. Specifically, conformational ensembles for the disordered lid region of the N-terminal domain of the oncoprotein MDM2 in the presence of different ligands were computed by means of a novel combination of accelerated molecular dynamics, umbrella sampling, and variational free energy profile methodologies. The resulting conformational ensembles for MDM2, free and bound to p53 TAD (17-29) peptide identify lid states compatible with previous NMR measurements. Remarkably, the MDM2 lid region is shown to adopt distinct conformational states in the presence of different small-molecule ligands. Detailed analyses of small-molecule bound ensembles reveal that the ca. 25-fold affinity improvement of the piperidinone family of inhibitors for MDM2 constructs that include the full lid correlates with interactions between ligand hydrophobic groups and the C-terminal lid region that is already partially ordered in apo MDM2. By contrast, Nutlin or benzodiazepinedione inhibitors, that bind with similar affinity to full lid and lid-truncated MDM2 constructs, interact additionally through their solubilizing groups with N-terminal lid residues that are more disordered in apo MDM2.

Life as we know it depends on interactions between proteins. There is substantial evidence that many interactions between proteins involve very flexible protein regions. These disordered regions may undergo disorder/order transitions upon forming an interaction with another protein. Many successful approaches to medicinal chemistry are based on mimicking the interactions of biological molecules with man-made small molecules. However how drug-like small-molecules may modulate protein disorder is currently poorly understood, largely because it is difficult to measure in details this type of interaction with experimental methods. Here we have used computer simulations to resolve with great details the process by which different small-molecules modulate the flexibility of a disordered region of the protein MDM2. This protein is overexpressed in many cancers and small-molecules that recognize MDM2 have been developed over the last decade as possible novel anti-cancer agents. We show that the flexible MDM2 “lid” region adopts different conformational states in the presence of different small-molecules. Our results suggest why some classes of small-molecules form favorable interactions with the lid region, whereas others do not. These findings may prove crucial to develop new and more effective MDM2 inhibitors, and more generally to help drug designers target disordered proteins regions with small-molecules.

Detailed Per-residue Energetic Analysis Explains the Driving Force for Microtubule Disassembly

Microtubules are long filamentous hollow cylinders whose surfaces form lattice structures of αβ-tubulin heterodimers. They perform multiple physiological roles in eukaryotic cells and are targets for therapeutic interventions. In our study, we carried out all-atom molecular dynamics simulations for arbitrarily long microtubules that have either GDP or GTP molecules in the E-site of β-tubulin. A detailed energy balance of the MM/GBSA inter-dimer interaction energy per residue contributing to the overall lateral and longitudinal structural stability was performed. The obtained results identified the key residues and tubulin domains according to their energetic contributions. They also identified the molecular forces that drive microtubule disassembly. At the tip of the plus end of the microtubule, the uneven distribution of longitudinal interaction energies within a protofilament generates a torque that bends tubulin outwardly with respect to the cylinder's axis causing disassembly. In the presence of GTP, this torque is opposed by lateral interactions that prevent outward curling, thus stabilizing the whole microtubule. Once GTP hydrolysis reaches the tip of the microtubule (lateral cap), lateral interactions become much weaker, allowing tubulin dimers to bend outwards, causing disassembly. The role of magnesium in the process of outward curling has also been demonstrated. This study also showed that the microtubule seam is the most energetically labile inter-dimer interface and could serve as a trigger point for disassembly. Based on a detailed balance of the energetic contributions per amino acid residue in the microtubule, numerous other analyses could be performed to give additional insights into the properties of microtubule dynamic instability.

Insights into the Interactions of Fasciola hepatica Cathepsin L3 with a Substrate and Potential Novel Inhibitors through In Silico Approaches

Background

Fasciola hepatica is the causative agent of fascioliasis, a disease affecting grazing animals, causing economic losses in global agriculture and currently being an important human zoonosis. Overuse of chemotherapeutics against fascioliasis has increased the populations of drug resistant parasites. F. hepatica cathepsin L3 is a protease that plays important roles during the life cycle of fluke. Due to its particular collagenolytic activity it is considered an attractive target against the infective phase of F. hepatica.

Methodology/Principal Findings

Starting with a three dimensional model of FhCL3 we performed a structure-based design of novel inhibitors through a computational study that combined virtual screening, molecular dynamics simulations, and binding free energy (ΔG_bind) calculations. Virtual screening was carried out by docking inhibitors obtained from the MYBRIDGE-HitFinder database inside FhCL3 and human cathepsin L substrate-binding sites. On the basis of dock-scores, five compounds were predicted as selective inhibitors of FhCL3. Molecular dynamic simulations were performed and, subsequently, an end-point method was employed to predict ΔG_bind values. Two compounds with the best ΔG_bind values (-10.68 kcal/mol and -7.16 kcal/mol), comparable to that of the positive control (-10.55 kcal/mol), were identified. A similar approach was followed to structurally and energetically characterize the interface of FhCL3 in complex with a peptidic substrate. Finally, through pair-wise and per-residue free energy decomposition we identified residues that are critical for the substrate/ligand binding and for the enzyme specificity.

Conclusions/Significance

The present study is the first computer-aided drug design approach against F. hepatica cathepsins. Here we predict the principal determinants of binding of FhCL3 in complex with a natural substrate by detailed energetic characterization of protease interaction surface. We also propose novel compounds as FhCL3 inhibitors. Overall, these results will foster the future rational design of new inhibitors against FhCL3, as well as other F. hepatica cathepsins.

Fascioliosis is considered an emerging disease in humans, causing important losses in global agriculture through the infection of livestock animals. The outcome of resistant parasites has increased the search for new drugs which may contribute to disease control. In recent decades, Fasciola cathepsins (FhCs) have been defined as the principal virulence factors of this parasite. Despite being in the same protein family, they have different specificities and, thus, distinct roles throughout the fluke life cycle. Differences in specificity have been attributed to a few variations in the sequence of key FhCs subsites. Currently, the structure-based drug design of inhibitors against Fasciola cathepsin Ls (FhCLs) with unknown structures is possible due to the availability of the three-dimensional structure of FhCL1. Our detailed structural analysis of the major infective juvenile enzyme (FhCL3) identifies the molecular determinants for protein binding. Also, novel potential inhibitors against FhCL3 are proposed, which might reduce host invasion and penetration processes. These compounds are predicted to interact with the binding site of the enzyme, therefore they could prevent substrate processing by competitive inhibition. The structure-based drug design strategy described here will be useful for the development of new potent and selective inhibitors against other FhCs.

Intercalation processes of copper complexes in DNA

The family of anticancer complexes that include the transition metal copper known as Casiopeínas® shows promising results. Two of these complexes are currently in clinical trials. The interaction of these compounds with DNA has been observed experimentally and several hypotheses regarding the mechanism of action have been developed, and these include the generation of reactive oxygen species, phosphate hydrolysis and/or base-pair intercalation. To advance in the understanding on how these ligands interact with DNA, we present a molecular dynamics study of 21 Casiopeínas with a DNA dodecamer using 10 μs of simulation time for each compound. All the complexes were manually inserted into the minor groove as the starting point of the simulations. The binding energy of each complex and the observed representative type of interaction between the ligand and the DNA is reported. With this extended sampling time, we found that four of the compounds spontaneously flipped open a base pair and moved inside the resulting cavity and four compounds formed stacking interactions with the terminal base pairs. The complexes that formed the intercalation pocket led to more stable interactions.

Structure Based Modeling of Small Molecules Binding to the TLR7 by Atomistic Level Simulations

Toll-Like Receptors (TLR) are a large family of proteins involved in the immune system response. Both the activation and the inhibition of these receptors can have positive effects on several diseases, including viral pathologies and cancer, therefore prompting the development of new compounds. In order to provide new indications for the design of Toll-Like Receptor 7 (TLR7)-targeting drugs, the mechanism of interaction between the TLR7 and two important classes of agonists (imidazoquinoline and adenine derivatives) was investigated through docking and Molecular Dynamics simulations. To perform the computational analysis, a new model for the dimeric form of the receptors was necessary and therefore created. Qualitative and quantitative differences between agonists and inactive compounds were determined. The in silico results were compared with previous experimental observations and employed to define the ligand binding mechanism of TLR7.

Investigate the binding of catechins to trypsin using docking and molecular dynamics simulation

To explore the inhibitory mechanism of catechins for digestive enzymes, we investigated the binding mode of catechins to a typical digestive enzyme-trypsin and analyzed the structure-activity relationship of catechins, using an integration of molecular docking, molecular dynamics simulation and binding free energy calculation. We found that catechins with different structures bound to a conservative pocket S1 of trypsin, which is comprised of residues 189–195, 214–220 and 225–228. In the trypsin-catechin complexes, Asp189 by forming strong hydrogen bonding, and Gln192, Trp215 and Gly216 through hydrophobic interactions, all significantly contribute to the binding of catechins. The number and the position of hydroxyl and aromatic groups, the structure of stereoisomers, and the orientation of catechins in the binding pocket S1 of trypsin all affect the binding affinity. The binding affinity is in the order of Epigallocatechin gallate (EGCG) > Epicatechin gallate (ECG) > Epicatechin (EC) > Epigallocatechin (EGC), and 2R-3R EGCG shows the strongest binding affinity out of other stereoisomers. Meanwhile, the synergic conformational changes of residues and catechins were also analyzed. These findings will be helpful in understanding the knowledge of interactions between catechins and trypsin and referable for the design of novel polyphenol based functional food and nutriceutical formulas.

Screening of novel inhibitors targeting lactate dehydrogenase A via four molecular docking strategies and dynamics simulations

Lactate dehydrogenase A (LDHA) is a metabolic enzyme which catalyzes the interconversion of lactate and pyruvate in the glycolysis pathway, thus playing key roles in aerobic glycolysis. The inhibition of LDHA by small molecules has become an attractive strategy for anticancer therapy in recent years. However, very few LDHA inhibitors have been reported, even though a great deal of effort has directed into identifying LDHA inhibitors using structure-based approaches. Therefore, high-throughput and high-accuracy screening approaches are still urgently needed in order to target LDHA effectively. In the present work, after establishing that our docking strategies performed well using test datasets, we screened 32791 Specs products for their docking scores with the substrate-binding pocket and, separately, the cofactor-binding pocket of LDHA. We subsequently identified 76 hits (i.e., ligands that show low docking scores) for the cofactor-binding pocket and 27 hits for the substrate-binding pocket. Two representative compounds, ZINC20036549 and ZINC19369718, were then chosen for further MD simulation analysis, and we found that these compounds maintained their inhibitory activity during the MD simulations. Meanwhile, we found that ZINC19369718 interacts with a novel binding site close to the active site, and that this interaction may inhibit the catalytic activity of LDHA. Together, these results offer not only a new paradigm for identifying Specs drug-like products for novel therapeutic use but they also provide further opportunity to adopt LDHA inhibition as a strategy for cancer therapy.

The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities

INTRODUCTION

The molecular mechanics energies combined with the Poisson-Boltzmann or generalized Born and surface area continuum solvation (MM/PBSA and MM/GBSA) methods are popular approaches to estimate the free energy of the binding of small ligands to biological macromolecules. They are typically based on molecular dynamics simulations of the receptor-ligand complex and are therefore intermediate in both accuracy and computational effort between empirical scoring and strict alchemical perturbation methods. They have been applied to a large number of systems with varying success.

AREAS COVERED

The authors review the use of MM/PBSA and MM/GBSA methods to calculate ligand-binding affinities, with an emphasis on calibration, testing and validation, as well as attempts to improve the methods, rather than on specific applications.

EXPERT OPINION

MM/PBSA and MM/GBSA are attractive approaches owing to their modular nature and that they do not require calculations on a training set. They have been used successfully to reproduce and rationalize experimental findings and to improve the results of virtual screening and docking. However, they contain several crude and questionable approximations, for example, the lack of conformational entropy and information about the number and free energy of water molecules in the binding site. Moreover, there are many variants of the method and their performance varies strongly with the tested system. Likewise, most attempts to ameliorate the methods with more accurate approaches, for example, quantum-mechanical calculations, polarizable force fields or improved solvation have deteriorated the results.

The Interaction of Heparin Tetrasaccharides with Chemokine CCL5 Is Modulated by Sulfation Pattern and pH

Interactions between chemokines such as CCL5 and glycosaminoglycans (GAGs) are essential for creating haptotactic gradients to guide the migration of leukocytes into inflammatory sites, and the GAGs that interact with CCL5 with the highest affinity are heparan sulfates/heparin. The interaction between CCL5 and its receptor on monocytes, CCR1, is mediated through residues Arg-17 and -47 in CCL5, which overlap with the GAG-binding (44)RKNR(47) "BBXB" motifs. Here we report that heparin and tetrasaccharide fragments of heparin are able to inhibit CCL5-CCR1 binding, with IC50 values showing strong dependence on the pattern and extent of sulfation. Modeling of the CCL5-tetrasaccharide complexes suggested that interactions between specific sulfate and carboxylate groups of heparin and residues Arg-17 and -47 of the protein are essential for strong inhibition; tetrasaccharides lacking the specific sulfation pattern were found to preferentially bind CCL5 in positions less favorable for inhibition of the interaction with CCR1. Simulations of a 12-mer heparin fragment bound to CCL5 indicated that the oligosaccharide preferred to interact simultaneously with both (44)RKNR(47) motifs in the CCL5 homodimer and engaged residues Arg-47 and -17 from both chains. Direct engagement of these residues by the longer heparin oligosaccharide provides a rationalization for its effectiveness as an inhibitor of CCL5-CCR1 interaction. In this mode, histidine (His-23) may contribute to CCL5-GAG interactions when the pH drops just below neutral, as occurs during inflammation. Additionally, an examination of the contribution of pH to modulating CCL5-heparin interactions suggested a need for careful interpretation of experimental results when experiments are performed under non-physiological conditions.

Molecular mechanism of HIV-1 Tat interacting with human dopamine transporter

Nearly 70% of HIV-1-infected individuals suffer from HIV-associated neurocognitive disorders (HAND). HIV-1 transactivator of transcription (Tat) protein is known to synergize with abused drugs and exacerbate the progression of central nervous system (CNS) pathology. Cumulative evidence suggest that the HIV-1 Tat protein exerts the neurotoxicity through interaction with human dopamine transporter (hDAT) in the CNS. Through computational modeling and molecular dynamics (MD) simulations, we develop a three-dimensional (3D) structural model for HIV-1 Tat binding with hDAT. The model provides novel mechanistic insights concerning how HIV-1 Tat interacts with hDAT and inhibits dopamine uptake by hDAT. In particular, according to the computational modeling, Tat binds most favorably with the outward-open state of hDAT. Residues Y88, K92, and Y470 of hDAT are predicted to be key residues involved in the interaction between hDAT and Tat. The roles of these hDAT residues in the interaction with Tat are validated by experimental tests through site-directed mutagensis and dopamine uptake assays. The agreement between the computational and experimental data suggests that the computationally predicted hDAT-Tat binding mode and mechanistic insights are reasonable and provide a new starting point to design further pharmacological studies on the molecular mechanism of HIV-1-associated neurocognitive disorders.

Structural and functional roles of glycosylation in fungal laccase from Lentinus sp

Laccases are multi-copper oxidases that catalyze the oxidation of various organic and inorganic compounds by reducing O2 to water. Here we report the crystal structure at 1.8 Å resolution of a native laccase (designated nLcc4) isolated from a white-rot fungus Lentinus sp. nLcc4 is composed of three cupredoxin-like domains D1-D3 each folded into a Greek key β-barrel topology. T1 and T2/T3 copper binding sites and three N-glycosylated sites at Asn75, Asn238, and Asn458 were elucidated. Initial rate kinetic analysis revealed that the kcat, Km, and kcat/Km of nLcc4 with substrate ABTS were 3,382 s-1, 65.0 ± 6.5 μM, and 52 s-1μM-1, respectively; and the values with lignosulfonic acid determined using isothermal titration calorimetry were 0.234 s-1, 56.7 ± 3.2 μM, and 0.004 s-1μM-1, respectively. Endo H-deglycosylated nLcc4 (dLcc4), with only one GlcNAc residue remaining at each of the three N-glycosylation sites in the enzyme, exhibited similar kinetic efficiency and thermal stability to that of nLcc4. The isolated Lcc4 gene contains an open reading frame of 1563 bp with a deduced polypeptide of 521 amino acid residues including a predicted signaling peptide of 21 residues at the N-terminus. Recombinant wild-type Lcc4 and mutant enzymes N75D, N238D and N458D were expressed in Pichia pastoris cells to evaluate the effect on enzyme activity by single glycosylation site deficiency. The mutant enzymes secreted in the cultural media of P. pastoris cells were observed to maintain only 4-50% of the activity of the wild-type laccase. Molecular dynamics simulations analyses of various states of (de-)glycosylation in nLcc support the kinetic results and suggest that the local H-bond networks between the domain connecting loop D2-D3 and the glycan moieties play a crucial role in the laccase activity. This study provides new insights into the role of glycosylation in the structure and function of a Basidiomycete fungal laccase.

Improvements, trends, and new ideas in molecular docking: 2012-2013 in review

Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher-level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods.

Structural determinants of host specificity of complement Factor H recruitment by Streptococcus pneumoniae

Many human pathogens have strict host specificity, which affects not only their epidemiology but also the development of animal models and vaccines. Complement Factor H (FH) is recruited to pneumococcal cell surface in a human-specific manner via the N-terminal domain of the pneumococcal protein virulence factor choline-binding protein A (CbpAN). FH recruitment enables Streptococcus pneumoniae to evade surveillance by human complement system and contributes to pneumococcal host specificity. The molecular determinants of host specificity of complement evasion are unknown. In the present study, we show that a single human FH (hFH) domain is sufficient for tight binding of CbpAN, present the crystal structure of the complex and identify the critical structural determinants for host-specific FH recruitment. The results offer new approaches to the development of better animal models for pneumococcal infection and redesign of the virulence factor for pneumococcal vaccine development and reveal how FH recruitment can serve as a mechanism for both pneumococcal complement evasion and adherence.

Molecular dynamics study of the inhibitory effects of ChEMBL474807 on the enzymes GSK-3β and CDK-2

Indirubin derivatives and analogs comprise a significant group of ATP-competitive inhibitors. The inhibitory effects of ChEMBL474807 (1-(4-amino-1,2,5-oxadiazol-3-yl)-5-(piperidin-1-ylmethyl)-N′-(pyridin-4-ylmethylene)-1H-1,2,3-triazole-4-carbohydrazide) on two enzymes, namely glycogen synthase kinase-3β (GSK-3β) and cyclin-dependent kinase-2 (CDK-2), were analyzed. The close resemblance of the amino acid sequences of these two enzymes (with 25 % identity and 41 % similarity) explains why indirubin derivatives are inhibitors of both of the enzymes studied. The docking and molecular dynamics investigation performed here led to the identification of the interactions responsible for stabilizing the ligand ChEMBL474807 at the active sites of the enzymes considered. The structural and energetic data collected during our investigations clearly indicate that there are important differences in the behavior of the ligand at the two active sites investigated here.

Structure based approach for understanding organism specific recognition of protein-RNA complexes

Background

Protein-RNA interactions perform diverse functions within the cell. Understanding the recognition mechanism of protein-RNA complexes has been a challenging task in molecular and computational biology. In earlier works, the recognition mechanisms have been studied for a specific complex or using a set of non–redundant complexes. In this work, we have constructed 18 sets of same protein-RNA complexes belonging to different organisms from Protein Data Bank (PDB). The similarities and differences in each set of complexes have been revealed in terms of various sequence and structure based features such as root mean square deviation, sequence homology, propensity of binding site residues, variance, conservation at binding sites, binding segments, binding motifs of amino acid residues and nucleotides, preferred amino acid-nucleotide pairs and influence of neighboring residues for binding.

Results

We found that the proteins of mesophilic organisms have more number of binding sites than thermophiles and the binding propensities of amino acid residues are distinct in E. coli, H. sapiens, S. cerevisiae, thermophiles and archaea. Proteins prefer to bind with RNA using a single residue segment in all the organisms while RNA prefers to use a stretch of up to six nucleotides for binding with proteins. We have developed amino acid residue-nucleotide pair potentials for different organisms, which could be used for predicting the binding specificity. Further, molecular dynamics simulation studies on aspartyl tRNA synthetase complexed with aspartyl tRNA showed specific modes of recognition in E. coli, T. thermophilus and S. cerevisiae.

Conclusion

Based on structural analysis and molecular dynamics simulations we suggest that the mode of recognition depends on the type of the organism in a protein-RNA complex.

Reviewers

This article was reviewed by Sandor Pongor, Gajendra Raghava and Narayanaswamy Srinivasan.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0039-8) contains supplementary material, which is available to authorized users.

Tunable assembly of amyloid-forming peptides into nanosheets as a retrovirus carrier

Using and engineering amyloid as nanomaterials are blossoming trends in bionanotechnology. Here, we show our discovery of an amyloid structure, termed "amyloid-like nanosheet," formed by a key amyloid-forming segment of Alzheimer's Aβ. Combining multiple biophysical and computational approaches, we proposed a structural model for the nanosheet that is formed by stacking the amyloid fibril spines perpendicular to the fibril axis. We further used the nanosheet for laboratorial retroviral transduction enhancement and directly visualized the presence of virus on the nanosheet surface by electron microscopy. Furthermore, based on our structural model, we designed nanosheet-forming peptides with different functionalities, elucidating the potential of rational design for amyloid-based materials with novel architecture and function.

Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Acα(2-3)Gal: a study by in silico mutations and molecular dynamics simulations

Wheat germ agglutinin (WGA) is a plant lectin, which specifically recognizes the sugars NeuNAc and GlcNAc. Mutated WGA with enhanced binding specificity can be used as biomarkers for cancer. In silico mutations are performed at the active site of WGA to enhance the binding specificity towards sialylglycans, and molecular dynamics simulations of 20 ns are carried out for wild type and mutated WGAs (WGA1, WGA2, and WGA3) in complex with sialylgalactose to examine the change in binding specificity. MD simulations reveal the change in binding specificity of wild type and mutated WGAs towards sialylgalactose and bound conformational flexibility of sialylgalactose. The mutated polar amino acid residues Asn114 (S114N), Lys118 (G118K), and Arg118 (G118R) make direct and water mediated hydrogen bonds and hydrophobic interactions with sialylgalactose. An analysis of possible hydrogen bonds, hydrophobic interactions, total pair wise interaction energy between active site residues and sialylgalactose and MM-PBSA free energy calculation reveals the plausible binding modes and the role of water in stabilizing different binding modes. An interesting observation is that the binding specificity of mutated WGAs (cyborg lectin) towards sialylgalactose is found to be higher in double point mutation (WGA3). One of the substituted residues Arg118 plays a crucial role in sugar binding. Based on the interactions and energy calculations, it is concluded that the order of binding specificity of WGAs towards sialylgalactose is WGA3 > WGA1 > WGA2 > WGA. On comparing with the wild type, double point mutated WGA (WGA3) exhibits increased specificity towards sialylgalactose, and thus, it can be effectively used in targeted drug delivery and as biological cell marker in cancer therapeutics.

Pentacycloundecane lactam vs lactone norstatine type protease HIV inhibitors: binding energy calculations and DFT study

Background

Novel pentacycloundecane (PCU)-lactone-CO-EAIS peptide inhibitors were designed, synthesized, and evaluated against wild-type C-South African (C-SA) HIV-1 protease. Three compounds are reported herein, two of which displayed IC₅₀ values of less than 1.00 μM. A comparative MM-PB(GB)SA binding free energy of solvation values of PCU-lactam and lactone models and their enantiomers as well as the PCU-lactam-NH-EAIS and lactone-CO-EAIS peptide inhibitors and their corresponding diastereomers complexed with South African HIV protease (C-SA) was performed. This will enable us to rationalize the considerable difference between inhibitory concentration (IC₅₀) of PCU-lactam-NH-EAIS and PCU-lactone-CO-EAIS peptides.

Results

The PCU-lactam model exhibited more negative calculated binding free energies of solvation than the PCU-lactone model. The same trend was observed for the PCU-peptide inhibitors, which correspond to the experimental activities for the PCU-lactam-NH-EAIS peptide (IC₅₀ = 0.076 μM) and the PCU-lactone-CO-EAIS peptide inhibitors (IC₅₀ = 0.850 μM). Furthermore, a density functional theory (DFT) study on the natural atomic charges of the nitrogen and oxygen atoms of the three PCU-lactam, PCU-lactim and PCU-lactone models were performed using natural bond orbital (NBO) analysis. Electrostatic potential maps were also used to visualize the electron density around electron-rich regions. The asymmetry parameter (η) and quadrupole coupling constant (χ) values of the nitrogen and oxygen nuclei of the model compounds were calculated at the same level of theory. Electronic molecular properties including polarizability and electric dipole moments were also calculated and compared. The Gibbs theoretical free solvation energies of solvation (∆G_solv) were also considered.

Conclusions

A general trend is observed that the lactam species appears to have a larger negative charge distribution around the heteroatoms, larger quadrupole constant, dipole moment and better solvation energy, in comparison to the PCU-lactone model. It can be argued that these characteristics will ensure better eletronic interaction between the lactam and the receptor, corresponding to the observed HIV protease activities in terms of experimental IC₅₀ data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-015-0115-5) contains supplementary material, which is available to authorized users.

A marine analgesic peptide, Contulakin-G, and neurotensin are distinct agonists for neurotensin receptors: uncovering structural determinants of desensitization properties

Neurotensin receptors have been studied as molecular targets for the treatment of pain, schizophrenia, addiction, or cancer. Neurotensin (NT) and Contulakin-G, a glycopeptide isolated from a predatory cone snail Conus geographus, share a sequence similarity at the C-terminus, which is critical for activation of neurotensin receptors. Both peptides are potent analgesics, although affinity and agonist potency of Contulakin-G toward neurotensin receptors are significantly lower, as compared to those for NT. In this work, we show that the weaker agonist properties of Contulakin-G result in inducing significantly less desensitization of neurotensin receptors and preserving their cell-surface density. Structure-activity relationship (SAR) studies suggested that both glycosylation and charged amino acid residues in Contulakin-G or NT played important roles in desensitizing neurotensin receptors. Computational modeling studies of human neurotensin receptor NTS1 and Contulakin-G confirmed the role of glycosylation in weakening interactions with the receptors. Based on available SAR data, we designed, synthesized, and characterized an analog of Contulakin-G in which the glycosylated amino acid residue, Gal-GalNAc-Thr10, was replaced by memantine-Glu10 residue. This analog exhibited comparable agonist potency and weaker desensitization properties as compared to that of Contulakin-G, while producing analgesia in the animal model of acute pain following systemic administration. We discuss our study in the context of feasibility and safety of developing NT therapeutic agents with improved penetration across the blood-brain barrier. Our work supports engineering peptide-based agonists with diverse abilities to desensitize G-protein coupled receptors and further emphasizes opportunities for conotoxins as novel pharmacological tools and drug candidates.

Dynamical insights of Mnk2 kinase activation by phosphorylation to facilitate inhibitor discovery

Aim: Mitogen-activated protein kinase-interacting kinases (Mnks) are emerging anticancer targets. Mnks feature unique structural features, enhancing their importance for selective inhibitor discovery. Nonetheless, the lack of structural details obstruct the development of selective Mnk inhibitors. Results: We disclose the first complete structure model of the activated state of Mnk2. Using all-atom accelerated molecular dynamics, we also demonstrate that its activation by phosphorylation grants access to distinct activation loop conformations, steering the inactive-to-active conformational transformation. Then we propose the binding mode of CGP57380 to active Mnk2, and evaluate key interactions that could be critical for future Mnk-targeted inhibitors. Conclusion: Critical insights of the Mnk2 activation process are gained, while providing a platform for designing Mnk-targeted anticancer agents.

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.

Mutations at tyrosine 88, lysine 92 and tyrosine 470 of human dopamine transporter result in an attenuation of HIV-1 Tat-induced inhibition of dopamine transport

HIV-1 transactivator of transcription (Tat) protein disrupts the dopamine (DA) neurotransmission by inhibiting DA transporter (DAT) function, leading to increased neurocognitive impairment in HIV-1 infected individuals. Through integrated computational modeling and pharmacological studies, we have demonstrated that mutation of tyrosine470 (Y470H) of human DAT (hDAT) attenuates Tat-induced inhibition of DA uptake by changing the transporter conformational transitions. The present study examined the functional influences of other substitutions at tyrosine470 (Y470F and Y470A) and tyrosine88 (Y88F) and lysine92 (K92M), two other relevant residues for Tat binding to hDAT, in Tat-induced inhibitory effects on DA transport. Y88F, K92M and Y470A attenuated Tat-induced inhibition of DA transport, implicating the functional relevance of these residues for Tat binding to hDAT. Compared to wild type hDAT, Y470A and K92M but not Y88F reduced the maximal velocity of [(3)H]DA uptake without changes in the Km. Y88F and K92M enhanced IC50 values for DA inhibition of [(3)H]DA uptake and [(3)H]WIN35,428 binding but decreased IC50 for cocaine and GBR12909 inhibition of [(3)H]DA uptake, suggesting that these residues are critical for substrate and these inhibitors. Y470F, Y470A, Y88F and K92M attenuated zinc-induced increase of [(3)H]WIN35,428 binding. Moreover, only Y470A and K92M enhanced DA efflux relative to wild type hDAT, suggesting mutations of these residues differentially modulate transporter conformational transitions. These results demonstrate Tyr88 and Lys92 along with Tyr470 as functional recognition residues in hDAT for Tat-induced inhibition of DA transport and provide mechanistic insights into identifying target residues on the DAT for Tat binding.

Structural dynamics of native and V260E mutant C-terminal domain of HIV-1 integrase

The C-terminal domain (CTD) of HIV-1 integrase is a five stranded β-barrel resembling an SH3 fold. Mutational studies on isolated CTD and full-length IN have reported V260E mutant as either homo-dimerization defective or affecting the stability and folding of CTD. In this study, molecular dynamics simulation techniques were used to unveil the effect of V260E mutation on isolated CTD monomer and dimer. Both monomeric and dimeric forms of wild type and V260E mutant are highly stable during the simulated period. However, the stabilizing π-stacking interaction between Trp243 and Trp243' at the dimer interface is highly disturbed in CTD-V260E (>6 Å apart). The loss in entropy for dimerization is -30 and -25 kcal/mol for CTD-wt and CTD-V260E respectively signifying a weak hydrophobic interaction and its perturbation in CTD-V260E. The mutant Glu260 exhibits strong attraction/repulsion with all the basic/acidic residues of CTD. In addition to this, the dynamics of CTD-wild type and V260E monomers at 498 K was analyzed to elucidate the effect of V260E mutation on CTD folding. Increase in SASA and reduction in the number of contacts in CTD-V260E during simulation highlights the instability caused by the mutation. In general, V260E mutation affects both multimerization and protein folding with a pronounced effect on protein folding rather than multimerization. This study emphasizes the importance of the hydrophobic nature and SH3 fold of CTD in proper functioning of HIV integrase and perturbing this nature would be a rational approach toward designing more selective and potent allosteric anti-HIV inhibitors.

Exploring resistance mechanisms of HCV NS3/4A protease mutations to MK5172: insight from molecular dynamics simulations and free energy calculations

Mutations at a number of key positions (Ala156, Asp168 and Arg155) of the HCV NS3/4A protease can induce medium to high resistance to MK5172.

Key role of the linker in pyrene-linker-carboxylate surfactants for the efficient aqueous dispersion of multiwalled carbon nanotubes

This study establishes that the structure of the linker group between the hydrophobic and hydrophilic units in the new surfactants 1–13 plays an important role in the dispersibility of multiwalled carbon nanotubes (MWNTs) in water.

The impact of active site protonation on substrate ring conformation in Melanocarpus albomyces cellobiohydrolase Cel7B

Understanding how the protonation state of active site residues impacts the enzyme's structure and substrate conformation is important for improving the efficiency and economic viability of the degradation of cellulosic materials as feedstock for liquid fuel and value-added chemicals.

Structural study of a small molecule receptor bound to dimethyllysine in lysozyme

Lysine is a ubiquitous residue on protein surfaces. Post translational modifications of lysine, including methylation to the mono-, di- or trimethylated amine result in chemical and structural alterations that have major consequences for protein interactions and signalling pathways. Small molecules that bind to methylated lysines are potential tools to modify such pathways. To make progress in this direction, detailed structural data of ligands in complex with methylated lysine is required. Here, we report a crystal structure of p-sulfonatocalix[4]arene (sclx₄) bound to methylated lysozyme in which the lysine residues were chemically modified from Lys-NH₃⁺ to Lys-NH(Me₂)⁺. Of the six possible dimethyllysine sites, sclx₄ selected Lys116-Me₂ and the dimethylamino substituent was deeply buried in the calixarene cavity. This complex confirms the tendency for Lys-Me₂ residues to form cation-π interactions, which have been shown to be important in protein recognition of histone tails bearing methylated lysines. Supporting data from NMR spectroscopy and MD simulations confirm the selectivity for Lys116-Me₂ in solution. The structure presented here may serve as a stepping stone to the development of new biochemical reagents that target methylated lysines.

A combined molecular dynamic and quantum mechanic study of the solvent and guest molecule effect on the stability and length of heterocyclic peptide nanotubes

Guest molecules and solvents affect the stability and length of the heterocyclic peptide nanotube through the electrostatic interactions.

Activation helix orientation of the estrogen receptor is mediated by receptor dimerization: evidence from molecular dynamics simulations

ERα dimer formation reshapes the helix 12 conformational landscape and is a leading factor for the activation helix conformation.

Irreversible covalent modification of type I dehydroquinase with a stable Schiff base

Structural and computational studies carried out with two epoxides provide insight into the irreversible inhibition of type I dehydroquinase.

Mathematical and computational modeling in biology at multiple scales

A variety of topics are reviewed in the area of mathematical and computational modeling in biology, covering the range of scales from populations of organisms to electrons in atoms. The use of maximum entropy as an inference tool in the fields of biology and drug discovery is discussed. Mathematical and computational methods and models in the areas of epidemiology, cell physiology and cancer are surveyed. The technique of molecular dynamics is covered, with special attention to force fields for protein simulations and methods for the calculation of solvation free energies. The utility of quantum mechanical methods in biophysical and biochemical modeling is explored. The field of computational enzymology is examined.

Forcefield_NCAA: ab initio charge parameters to aid in the discovery and design of therapeutic proteins and peptides with unnatural amino acids and their application to complement inhibitors of the compstatin family

We describe the development and testing of ab initio derived, AMBER ff03 compatible charge parameters for a large library of 147 noncanonical amino acids including β- and N-methylated amino acids for use in applications such as protein structure prediction and de novo protein design. The charge parameter derivation was performed using the RESP fitting approach. Studies were performed assessing the suitability of the derived charge parameters in discriminating the activity/inactivity between 63 analogs of the complement inhibitor Compstatin on the basis of previously published experimental IC50 data and a screening procedure involving short simulations and binding free energy calculations. We found that both the approximate binding affinity (K*) and the binding free energy calculated through MM-GBSA are capable of discriminating between active and inactive Compstatin analogs, with MM-GBSA performing significantly better. Key interactions between the most potent Compstatin analog that contains a noncanonical amino acid are presented and compared to the most potent analog containing only natural amino acids and native Compstatin. We make the derived parameters and an associated web interface that is capable of performing modifications on proteins using Forcefield_NCAA and outputting AMBER-ready topology and parameter files freely available for academic use at http://selene.princeton.edu/FFNCAA . The forcefield allows one to incorporate these customized amino acids into design applications with control over size, van der Waals, and electrostatic interactions.

Exploring alternate states and oligomerization preferences of coiled-coils by de novo structure modeling

Homomeric coiled-coils can self-assemble into a wide range of structural states with different helix topologies and oligomeric states. In this study, we have combined de novo structure modeling with stability calculations to simultaneously predict structure and oligomeric states of homomeric coiled-coils. For dimers an asymmetric modeling protocol was developed. Modeling without symmetry constraints showed that backbone asymmetry is important for the formation of parallel dimeric coiled-coils. Collectively, our results demonstrate that high-resolution structure of coiled-coils, as well as parallel and antiparallel orientations of dimers and tetramers, can be accurately predicted from sequence. De novo modeling was also used to generate models of competing oligomeric states, which were used to compare stabilities and thus predict the native stoichiometry from sequence. In a benchmark set of 33 coiled-coil sequences, forming dimers to pentamers, up to 70% of the oligomeric states could be correctly predicted. The calculations demonstrated that the free energy of helix folding could be an important factor for determining stability and oligomeric state of homomeric coiled-coils. The computational methods developed here should be broadly applicable to studies of sequence-structure relationships in coiled-coils and the design of higher order assemblies with improved oligomerization specificity.

Energetic contributions of residues to the formation of early amyloid-β oligomers

Low-weight amyloid-β (Aβ) oligomers formed at early stages of oligomerization rather than fibril assemblies seem to be the toxic components that drive neurodegeneration in Alzheimer's disease. Unfortunately, detailed knowledge of the structure of these early oligomers at the residue level is not yet available. In this study, we performed all-atom explicit solvent molecular dynamics simulations to examine the oligomerization process of Aβ10-35 monomers when forming dimers, trimers, tetramers and octamers, with four independent simulations of a total simulated time of 3 μs for each oligomer system. The decomposition of the stability free energy by MM-GBSA methodology allowed us to unravel the network of energetic interactions that stabilize such oligomers. The contribution of the intermonomeric van der Waals term is the most significant energy feature of the oligomerization process, consistent with the so-called hydrophobic effect. Furthermore, the decomposition of the stability free energy into residues and residue-pairwise terms revealed that it is mainly apolar interactions between the three specific hydrophobic fragments 31-35 (C-terminal region), 17-20 (central hydrophobic core) and 12-14 (N-terminal region) that are responsible for such a favourable effect. The conformation in which the hydrophobic cthr-chc interaction is oriented perpendicularly is particularly important. We propose three other model substructures that favour the oligomerization process and can thus be considered as molecular targets for future inhibitors. Understanding Aβ oligomerization at the residue level could lead to more efficient design of inhibitors of this process.