Molecular dynamic simulations of the metallo-beta-lactamase from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor

Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites

Zn²⁺ is one of the most versatile biologically available metal ions, but accurate modeling of Zn²⁺ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn²⁺ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn²⁺ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-β-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn²⁺ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn²⁺ site geometries that could guide the further development of biomolecular Zn²⁺ models.

Visualizing correlated motion with HDBSCAN clustering

Correlated motion analysis provides a method for understanding communication between and dynamic similarities of biopolymer residues and domains. The typical equal-time correlation matrices-frequently visualized with pseudo-colorings or heat maps-quickly convey large regions of highly correlated motion but hide more subtle similarities of motion. Here we propose a complementary method for visualizing correlations within proteins (or general biopolymers) that quickly conveys intuition about which residues have a similar dynamic behavior. For grouping residues, we use the recently developed non-parametric clustering algorithm HDBSCAN. Although the method we propose here can be used to group residues using correlation as a similarity matrix-the most straightforward and intuitive method-it can also be used to more generally determine groups of residues which have similar dynamic properties. We term these latter groups "Dynamic Domains", as they are based not on spatial closeness but rather closeness in the column space of a correlation matrix. We provide examples of this method across three human proteins of varying size and function-the Nf-Kappa-Beta essential modulator, the clotting promoter Thrombin and the mismatch repair protein (dimer) complex MutS-alpha. Although the examples presented here are from all-atom molecular dynamics simulations, this visualization technique can also be used on correlations matrices built from any ensembles of conformations from experiment or computation.

Molecular dynamics simulations of aptamer-binding reveal generalized allostery in thrombin

Thrombin is an attractive target for antithrombotic therapy due to its central role in thrombosis and hemostasis as well as its role in inducing tumor growth, metastasis, and tumor invasion. The thrombin-binding DNA aptamer (TBA), is under investigation for anticoagulant drugs. Although aptamer binding experiments have been revealed various effects on thrombin's enzymatic activities, the detailed picture of the thrombin's allostery from TBA binding is still unclear. To investigate thrombin's response to the aptamer-binding at the molecular level, we compare the mechanical properties and free energy landscapes of the free and aptamer-bound thrombin using microsecond-scale all-atom GPU-based molecular dynamics simulations. Our calculations on residue fluctuations and coupling illustrate the allosteric effects of aptamer-binding at the atomic level, highlighting the exosite II, 60s, γ and the sodium loops, and the alpha helix region in the light chains involved in the allosteric changes. This level of details clarifies the mechanisms of previous experimentally demonstrated phenomena, and provides a prediction of the reduced autolysis rate after aptamer-binding. The shifts in thrombin's ensemble of conformations and free energy surfaces after aptamer-binding demonstrate that the presence of bound-aptamer restricts the conformational freedom of thrombin suggesting that conformational selection, i.e. generalized allostery, is the dominant mechanism of thrombin-aptamer binding. The profound perturbation on thrombin's mechanical and thermodynamic properties due to the aptamer-binding, which was revealed comprehensively as a generalized allostery in this work, may be exploited in further drug discovery and development.

Spectroscopy analysis and molecular dynamics studies on the binding of penicillin V and sulbactam to beta-lactamase II from Bacillus cereus

The molecular recognition and interaction of beta-lactamase II from Bacillus cereus (Bc II) with penicillin V (PV) and sulbactam (Sul) especially conformational changes of Bc II in the binding process were studied through spectroscopy analysis in combination with molecular dynamics (MD) simulation. The results show that in the binding process, a new coordination bond is observed between the Zn₂ of Bc II and the carboxyl-O of PV or Sul by replacing His204. Electrostatic interaction between Zn₂ and the ligand provide main driving force for the binding affinity. Compared with apo Bc II, there are mainly four loops showing significant conformational changes in ligand-bound Bc II. A weak conformational transformation from β-sheets to random coils is observed in the loop2 of ligand-bound Bc II. The conformational transformation may depend on the functional group and binding pose of the ligand, giving the binding pocket greater flexibility and accordingly allowing for an induced fit of the enzyme-ligand binding site around the newly introduced ligand. The change in the loop2 of ligand-bound Bc II may lead to the opening of the binding pocket of Bc II. Therefore, loop2 can be considered a gate for control of ligand access in Bc II, hence its dynamic response should be considered in new drug design and development.

Influence of V54M mutation in giant muscle protein titin: a computational screening and molecular dynamics approach

Recent genetic studies have revealed the impact of mutations in associated genes for cardiac sarcomere components leading to dilated cardiomyopathy (DCM). The cardiac sarcomere is composed of thick and thin filaments and a giant muscle protein known as titin or connectin. Titin interacts with T-cap/telethonin in the Z-line region and plays a vital role in regulating sarcomere assembly. Initially, we screened all the variants associated with giant protein titin and analyzed their impact with the aid of pathogenicity and stability prediction methods. V54M mutation found in the hydrophobic core region of the protein associated with abnormal clinical phenotype leads to DCM was selected for further analysis. To address this issue, we mapped the deleterious mutant V54M, modeled the mutant protein complex, and deciphered the impact of mutation on binding with its partner telethonin in the titin crystal structure of PDB ID: 1YA5 with the aid of docking analysis. Furthermore, two run molecular dynamics simulation was initiated to understand the mechanistic action of V54M mutation in altering the protein structure, dynamics, and stability. According to the results obtained from the repeated 50 ns trajectory files, the overall effect of V54M mutation was destabilizing and transition of bend to coil in the secondary structure was observed. Furthermore, MMPBSA elucidated that V54M found in the Z-line region of titin decreases the binding affinity of titin to Z-line proteins T-cap/telethonin thereby hindering the protein-protein interaction.

The Role of Active Site Flexible Loops in Catalysis and of Zinc in Conformational Stability of Bacillus cereus 569/H/9 β-Lactamase

Metallo-β-lactamases catalyze the hydrolysis of most β-lactam antibiotics and hence represent a major clinical concern. The development of inhibitors for these enzymes is complicated by the diversity and flexibility of their substrate-binding sites, motivating research into their structure and function. In this study, we examined the conformational properties of the Bacillus cereus β-lactamase II in the presence of chemical denaturants using a variety of biochemical and biophysical techniques. The apoenzyme was found to unfold cooperatively, with a Gibbs free energy of stabilization (ΔG(0)) of 32 ± 2 kJ·mol(-1) For holoBcII, a first non-cooperative transition leads to multiple interconverting native-like states, in which both zinc atoms remain bound in an apparently unaltered active site, and the protein displays a well organized compact hydrophobic core with structural changes confined to the enzyme surface, but with no catalytic activity. Two-dimensional NMR data revealed that the loss of activity occurs concomitantly with perturbations in two loops that border the enzyme active site. A second cooperative transition, corresponding to global unfolding, is observed at higher denaturant concentrations, with ΔG(0) value of 65 ± 1.4 kJ·mol(-1) These combined data highlight the importance of the two zinc ions in maintaining structure as well as a relatively well defined conformation for both active site loops to maintain enzymatic activity.

Insights into protein - DNA interactions, stability and allosteric communications: a computational study of mutSα-DNA recognition complexes

DNA mismatch repair proteins (MMR) maintain genetic stability by recognizing and repairing mismatched bases and insertion/deletion loops mistakenly incorporated during DNA replication, and initiate cellular response to certain types of DNA damage. Loss of MMR in mammalian cells has been linked to resistance to certain DNA damaging chemotherapeutic agents, as well as to increase risk of cancer. Mismatch repair pathway is considered to involve the concerted action of at least 20 proteins. The most abundant MMR mismatch-binding factor in eukaryotes, MutSα, recognizes and initiates the repair of base-base mismatches and small insertion/deletion. We performed molecular dynamics simulations on mismatched and damaged MutSα-DNA complexes. A comprehensive DNA binding site analysis of relevant conformations shows that MutSα proteins recognize the mismatched and platinum cross-linked DNA substrates in significantly different modes. Distinctive conformational changes associated with MutSα binding to mismatched and damaged DNA have been identified and they provide insight into the involvement of MMR proteins in DNA-repair and DNA-damage pathways. Stability and allosteric interactions at the heterodimer interface associated with the mismatch and damage recognition step allow for prediction of key residues in MMR cancer-causing mutations. A rigorous hydrogen bonding analysis for ADP molecules at the ATPase binding sites is also presented. Due to extended number of known MMR cancer causing mutations among the residues proved to make specific contacts with ADP molecules, recommendations for further studies on similar mutagenic effects were made.

Recapturing the Correlated Motions of Protein Using Coarse- Grained Models

Long-range interactions and allostery are important for many biological processes. Increasing numbers of studies, both experimental and computational, show that internal dynamics may play an important role in such behaviors. Investigating the dynamical effects of proteins, how- ever, is a challenging problem using all-atom molecular dynamics because of the length-scales and timescales involved. As a result, coarse-grained models are often implemented. Herein, we use three well-defined coarse-grained models: Go, Martini and Cafemol, and a small model protein Eglin C, which is readily studied via all-atom molecular dynamics, to examine if these coarse grained models can explore the dynamics of Eglin C accurately as well as to see how these models respond to mutations. We found that all three models can recapture the dynamics of Eglin C to a significant extent - where we focus on root-mean square fluctuations and correlated motions as dynamical measures - but that the Cafemol and Go models are superior. The best agreement with all-atom simulations is for structured regions of Eglin C.

A variety of roles for versatile zinc in metallo-β-lactamases

β-Lactamases inactivate the important β-lactam antibiotics by catalysing the hydrolysis of the β-lactam ring, thus. One class of these enzymes, the metallo-β-lactamases, bind two zinc ions at the active site and these play important roles in the catalytic mechanism.

Probing the effect of the non-active-site mutation Y229W in New Delhi metallo-β-lactamase-1 by site-directed mutagenesis, kinetic studies, and molecular dynamics simulations

New Delhi metallo-β-lactmase-1 (NDM-1) has attracted extensive attention for its high catalytic activities of hydrolyzing almost all β-lactam antibiotics. NDM-1 shows relatively higher similarity to subclass B1 metallo-β-lactmases (MβLs), but its residue at position 229 is identical to that of B2/B3 MβLs, which is a Tyr instead of a B1-MβL-conserved Trp. To elucidate the possible role of Y229 in the bioactivity of NDM-1, we performed mutagenesis study and molecular dynamics (MD) simulations. Although residue Y229 is spatially distant from the active site and not contacting directly with the substrate or zinc ions, the Y229W mutant was found to have higher k_cat and K_m values than those of wild-type NDM-1, resulting in 1∼7 fold increases in k_cat/K_m values against tested antibiotics. In addition, our MD simulations illustrated the enhanced flexibility of Loop 2 upon Y229W mutation, which could increase the kinetics of both substrate entrance (kon) and product egress (koff). The enhanced flexibility of Loop 2 might allow the enzyme to adjust the geometry of its active site to accommodate substrates with different structures, broadening its substrate spectrum. This study indicated the possible role of the residue at position 229 in the evolution of NDM-1.

Destabilization of the MutSα's protein-protein interface due to binding to the DNA adduct induced by anticancer agent carboplatin via molecular dynamics simulations

DNA mismatch repair (MMR) proteins maintain genetic integrity in all organisms by recognizing and repairing DNA errors. Such alteration of hereditary information can lead to various diseases, including cancer. Besides their role in DNA repair, MMR proteins detect and initiate cellular responses to certain type of DNA damage. Its response to the damaged DNA has made the human MMR pathway a useful target for anticancer agents such as carboplatin. This study indicates that strong, specific interactions at the interface of MutSα in response to the mismatched DNA recognition are replaced by weak, non-specific interactions in response to the damaged DNA recognition. Data suggest a severe impairment of the dimerization of MutSα in response to the damaged DNA recognition. While the core of MutSα is preserved in response to the damaged DNA recognition, the loss of contact surface and the rearrangement of contacts at the protein interface suggest a different packing in response to the damaged DNA recognition. Coupled in response to the mismatched DNA recognition, interaction energies, hydrogen bonds, salt bridges, and solvent accessible surface areas at the interface of MutSα and within the subunits are uncoupled or asynchronously coupled in response to the damaged DNA recognition. These pieces of evidence suggest that the loss of a synchronous mode of response in the MutSα's surveillance for DNA errors would possibly be one of the mechanism(s) of signaling the MMR-dependent programed cell death much wanted in anticancer therapies. The analysis was drawn from dynamics simulations.

Non-specificity and synergy at the binding site of the carboplatin-induced DNA adduct via molecular dynamics simulations of the MutSα-DNA recognition complex

MutSα is the most abundant mismatch-binding factor of human DNA mismatch repair (MMR) proteins. MMR maintains genetic stability by recognizing and repairing DNA defects. Failure to accomplish their function may lead to cancer. In addition, MutSα recognizes at least some types of DNA damage making it a target for anticancer agents. Here, complementing scarce experimental data, we report unique hydrogen-bonding motifs associated with the recognition of the carboplatin induced DNA damage by MutSα. These data predict that carboplatin and cisplatin induced damaging DNA adducts are recognized by MutSα in a similar manner. Our simulations also indicate that loss of base pairing at the damage site results in (1) non-specific binding and (2) changes in the atomic flexibility at the lesion site and beyond. To further quantify alterations at MutSα-DNA interface in response to damage recognition, non-bonding interactions and salt bridges were investigated. These data indicate (1) possible different packing and (2) disruption of the salt bridges at the MutSα-DNA interface in the damaged complex. These findings (1) underscore the general observation of disruptions at the MutSα-DNA interface and (2) highlight the nature of the anticancer effect of the carboplatin agent. The analysis was carried out from atomistic simulations.

Predicting the impact of deleterious mutations in the protein kinase domain of FGFR2 in the context of function, structure, and pathogenesis--a bioinformatics approach

Fibroblast growth factor receptor 2 (FGFR2) controls a wide range of biological functions by regulating the cellular proliferation, survival, migration and differentiation. A growing body of preclinical data demonstrated that deregulation of the FGFR signalling through genetic modification was observed in various types of cancers. However, the extent to which genetic modifications interfere with gene regulation and their involvement in cancer susceptibility remains largely unknown. In this work, we performed in silico profiling of harmful non-synonymous single nucleotide polymorphisms (SNPs) in the protein kinase domain of FGFR2. Tolerance index, position-specific independent count score, change in free energy score (ΔΔG), Eris and FoldX indicated that seven mutations were found to be deleterious and may alter the protein function and structure. Furthermore, based on physico-chemical properties, two mutations K659N and R747H were found to be most deleterious in protein kinase domain and taken for further structural analysis. Docking study showed a complete loss of binding affinity followed by interference in hydrogen bonding and surrounding residues due to K659N and R747H mutations. In order to elucidate the mechanism behind the impact of mutation that can generate a ripple effect throughout the protein structure and ultimately affect the function, in-depth molecular dynamics simulation and principal component analysis were performed. The obtained results indicate that K659N and R747H mutations have a distinct effect on the dynamic behaviour of FGFR2 protein. Our strategy may be helpful for understanding SNP effects on proteins with function and their role in human genetic diseases and for the development of novel pharmacological strategies.

Allosteric inhibition of VIM metallo-β-lactamases by a camelid nanobody

MβL (metallo-β-lactamase) enzymes are usually produced by multi-resistant Gram-negative bacterial strains and have spread worldwide. An approach on the basis of phage display was used to select single-domain antibody fragments (VHHs, also called nanobodies) that would inhibit the clinically relevant VIM (Verona integron-encoded MβL)-4 MβL. Out of more than 50 selected nanobodies, only the NbVIM_38 nanobody inhibited VIM-4. The paratope, inhibition mechanism and epitope of the NbVIM_38 nanobody were then characterized. An alanine scan of the NbVIM_38 paratope showed that its binding was driven by hydrophobic amino acids. The inhibitory concentration was in the micromolar range for all β-lactams tested. In addition, the inhibition was found to follow a mixed hyperbolic profile with a predominantly uncompetitive component. Moreover, substrate inhibition was recorded only after nanobody binding. These kinetic data are indicative of a binding site that is distant from the active site. This finding was confirmed by epitope mapping analysis that was performed using peptides, and which identified two stretches of amino acids in the L6 loop and at the end of the α2 helix. Because this binding site is distant from the active site and alters both the substrate binding and catalytic properties of VIM-4, this nanobody can be considered as an allosteric inhibitor.

The molecular origin of the MMR-dependent apoptosis pathway from dynamics analysis of MutSα-DNA complexes

The cellular response to DNA damage signaling by mismatch-repair (MMR) proteins is incompletely understood. It is generally accepted that MMR-dependent apoptosis pathway in response to DNA damage detection is independent of MMR's DNA repair function. In this study, we investigate correlated motions in response to the binding of mismatched and platinum cross-linked DNA fragments by MutSα, as derived from 50 ns molecular dynamics simulations. The protein dynamics in response to the mismatched and damaged DNA recognition suggests that MutSα signals their recognition through independent pathways providing evidence for the molecular origin of the MMR-dependent apoptosis. MSH2 subunit is indicated to play a key role in signaling both mismatched and damaged DNA recognition; localized and collective motions within the protein allow identifying sites on the MSH2 surface possible involved in recruiting proteins responsible for downstream events. Unlike in the mismatch complex, predicted key communication sites specific for the damage recognition are on the list of known cancer-causing mutations or deletions. This confirms MSH2's role in signaling DNA damage-induced apoptosis and suggests that defects in MMR alone is sufficient to trigger tumorigenesis, supporting the experimental evidence that MMR-damage response function could protect from the early occurrence of tumors. Identifying these particular communication sites may have implications for the treatment of cancers that are not defective for MMR, but are unable to function optimally for MMR-dependent responses following DNA damage such as the case of resistance to cisplatin.

The stability of a model substrate for topoisomerase 1-mediated DNA religation depends on the presence of mismatched base pairs

Topoisomerase 1 (Top1) enzymes regulate DNA superhelicity by forming covalent cleavage complexes that undergo controlled rotation. Substitution of nucleoside analogs at the +1 position of the DNA duplex relative to the Top1 cleavage site inhibits DNA religation. The reduced efficiency for Top1-mediated religation contributes to the anticancer activity of widely used anticancer drugs including fluoropyrimidines and gemcitabine. In the present study, we report that mismatched base pairs at the +1 position destabilize the duplex DNA components for a model Top1 cleavage complex formation even though one duplex component does not directly include a mismatched base pair. Molecular dynamics simulations reveal G-dU and G-FdU mismatched base pairs, but not a G-T mismatched base pair, increase flexibility at the Top1 cleavage site, and affect coupling between the regions required for the religation reaction to occur. These results demonstrate that substitution of dT analogs into the +1 position of the non-scissile strand alters the stability and flexibility of DNA contributing to the reduced efficiency for Top1-mediated DNA religation. These effects are inherent in the DNA duplex and do not require formation of the Top1:DNA complex. These results provide a biophysical rationale for the inhibition of Top1-mediated DNA religation by nucleotide analog substitution.

Molecular dynamics simulations of protein dynamics and their relevance to drug discovery

Molecular dynamics simulations have become increasingly useful in studying biological systems of biomedical interest, and not just in the study of model or toy systems. In this article, the methods and principles of all-atom molecular dynamics will be elucidated with several examples provided of their utility to investigators interested on drug discovery.

Distribution, detection of enterotoxigenic strains and antimicrobial drug susceptibility patterns of bacteroides fragilis group in diarrheic and non-diarrheic feces from brazilian infants

Despite the importance of gastrointestinal diseases and their global distribution, affecting millions of individuals around the world, the role and antimicrobial susceptibility patterns of anaerobic bacteria such as those in the Bacteroides fragilis group (BFG) are still unclear in young children. This study investigated the occurrence and distribution of species in the BFG and enterotoxigenic strains in the fecal microbiota of children and their antimicrobial susceptibility patterns. Diarrheic (n=110) and non-diarrheic (n=65) fecal samples from children aged 0-5 years old were evaluated. BFG strains were isolated and identified by conventional biochemical, physiological and molecular approaches. Alternatively, bacteria and enterotoxigenic strains were detected directly from feces by molecular biology. Antimicrobial drug susceptibility patterns were determined by the agar dilution method according to the guidelines for isolated bacteria. BFG was detected in 64.3% of the fecal samples (55% diarrheic and 80.4% non-diarrheic), and 4.6% were enterotoxigenic. Antimicrobial resistance was observed against ampicillin, ampicillin/sulbactam, piperacillin/tazobactam, meropenem, ceftriaxone, clindamycin and chloramphenicol. The data show that these bacteria are prevalent in fecal microbiota at higher levels in healthy children. The molecular methodology was more effective in identifying the B. fragilis group when compared to the biochemical and physiological techniques. The observation of high resistance levels stimulates thoughts about the indiscriminate use of antimicrobial drugs in early infancy. Further quantitative studies are needed to gain a better understanding of the role of these bacteria in acute diarrhea in children.