HMGB1–Carbenoxolone Interactions: Dynamics Insights from Combined Nuclear Magnetic Resonance and Molecular Dynamics

The interplay of protein dynamics and molecular recognition is of fundamental importance in biological processes. Atomic‐resolution insights into these phenomena may provide new opportunities for drug discovery. Herein, we have combined NMR relaxation experiments and residual dipolar coupling (RDC) measurements with molecular dynamics (MD) simulations to study the effects of the anti‐inflammatory drug carbenoxolone (CBNX) on the conformational properties and on the internal dynamics of a subdomain (box A) of high‐mobility group B protein (HMGB1). ¹⁵N relaxation data show that CBNX binding enhances the fast pico‐ to nanosecond motions of a loop and partially removes the internal motional anisotropy of the first two helices of box A. Dipolar wave analysis of amide RDC data shows that ligand binding induces helical distortions. In parallel, increased mobility of the loop upon ligand binding is highlighted by the essential dynamics analysis (EDA) of MD simulations. Moreover, simulations detect two possible orientations for CBNX, which induces two possible conformations of helix H3, one being similar to the free form and the second one causing a partial helical distortion. Finally, we introduce a new approach for the analysis of the internal coordination of protein residues that is consistent with experimental data and allows us to pinpoint which substructures of box A are dynamically affected by CBNX. The observations reported here may be useful for understanding the role of protein dynamics in binding at atomic resolution.     

Study of the interaction of GFG tripeptide with cesium perfluorooctanoate micelles by means of NMR spectroscopy and MD simulations

The interaction of glycyl-phenylalanyl-glycine (GFG) with bilayers formed by cesium perfluorooctanoate (CsPFO) in water was investigated in the isotropic phase by means of 1H NMR and molecular dynamics (MD) simulations. Details on the preferential location of the different residues of GFG were obtained from selective variations of chemical shift with peptide concentration and of line width in the presence of the paramagnetic ion Mn2+. The analysis of 1H NMR spectra recorded at different concentrations and temperatures allowed the association constant and the enthalpy change upon binding to be evaluated. MD simulations highlighted the hydrogen bonds formed between the different GFG functional groups and the micelle. Both NMR and MD gave indications of high affinity of GFG with the micelle, with the N-terminal residue anchoring on the surface via hydrogen bonds with the micelle COO(-) groups.     

Conformational selection or induced fit for Brinker and DNA recognition

Brinker is the key target protein of the Drosophila Decapentaplegic morphogen signalling pathway. Brinker is widely expressed and can bind with DNA. NMR spectra suggest that apo-Brinker is intrinsically unstructured and undergoes a folding transition upon DNA-binding. However, the coupled mechanism of binding and folding is poorly understood. Here, we performed molecular dynamics (MD) simulations for both bound and apo-Brinker to study the mechanism. Room-temperature MD simulations suggest that Brinker becomes more rigid and stable upon DNA-binding. Kinetic analysis of high-temperature MD simulations shows that both bound and apo-Brinker unfold via a two-state process. The time scale of tertiary unfolding is significantly different between bound and apo-Brinker. The predicted Φ-values suggest that there are more residues with native-like transition state ensembles (TSEs) for bound Brinker than for apo-Brinker. The average RMSD differences between bound and apo-Brinker and Kolmogorov-Smirnov (KS) test analysis illustrate that Brinker folding upon DNA-binding might obey induced-fit mechanism based on MD simulations. These methods can be used for the research of other biomolecular folding upon ligand-binding.     

What induces pocket openings on protein surface patches involved in protein–protein interactions?

We previously showed for the proteins BCL-X_L, IL-2, and MDM2 that transient pockets at their protein–protein binding interfaces can be identified by applying the PASS algorithm to molecular dynamics (MD) snapshots. We now investigated which aspects of the natural conformational dynamics of proteins induce the formation of such pockets. The pocket detection protocol was applied to three different conformational ensembles for the same proteins that were extracted from MD simulations of the inhibitor bound crystal conformation in water and the free crystal/NMR structure in water and in methanol. Additional MD simulations studied the impact of backbone mobility. The more efficient CONCOORD or normal mode analysis (NMA) techniques gave significantly smaller pockets than MD simulations, whereas tCONCOORD generated pockets comparable to those observed in MD simulations for two of the three systems. Our findings emphasize the influence of solvent polarity and backbone rearrangements on the formation of pockets on protein surfaces and should be helpful in future generation of transient pockets as putative ligand binding sites at protein–protein interfaces. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Implications of fast-time scale dynamics of human DNA/RNA cytosine methyltransferases (DNMTs) for protein function

The role of protein dynamics in the control of substrate recognition, catalysis, and protein–protein interactions is often underestimated. Recently, a number of studies have examined the contribution of protein dynamics to the thermodynamics of ligand binding in detail, mostly using NMR relaxation measurements and molecular dynamics (MD) simulations. The results unequivocally demonstrate that conformational dynamics play a pivotal role in the properties and functions of proteins, and ignoring this contribution is likely to lead to substantial errors when explaining the biological function of proteins and in predictions of the binding affinities of their cognate ligands. However, the details of the interplay between structure and dynamics and the way it affects the biological function of the target protein remain poorly understood. In this study, the changes in fast (picosecond-to-nanosecond time scale) dynamics of catalytic domains of four human cytosine DNA methyltransferases (DNMTs) were studied using molecular dynamics (MD) simulations. The results provide insight into the protein dynamics changes that occur upon binding of the cofactor, S-adenosylmethionine (SAM). Contrary to expectations, increased amplitude of motions of backbone amide (N–H) and terminal heavy atom (C–C) bond vectors was observed in all studied DNMTs upon binding of SAM. These results imply that the cofactor binding causes a global increase in the extent of protein dynamics in the short time scale. This global dynamic change constitutes a favourable entropic contribution to the free energy of SAM binding. These results suggest that cytosine DNA methyltransferases may exploit changes in their fast scale dynamics to reduce the entropic cost of the substrate binding.     

Conformational Selection in Glycomimetics: Human Galectin‐1 Only Recognizes syn‐Ψ‐Type Conformations of β‐1,3‐Linked Lactose and Its C‐Glycosyl Derivative

The human lectin galectin‐1 (hGal‐1) translates sugar signals, that is, β‐galactosides, into effects on the level of cells, for example, growth regulation, and has become a model for studying binding of biopharmaceutically relevant derivatives. Bound‐state conformations of Galβ‐C‐(1→3)‐Glcβ‐OMe ( 1 ) and its βGal‐(1→3)‐βGlc‐OMe disaccharide parent compound were studied by using NMR spectroscopy (transferred (TR)‐NOESY data), assisted by docking experiments and molecular dynamics (MD) simulations. The molecular recognition process involves a conformational selection event. Although free C‐glycoside access four distinct conformers in solution, hGal‐1 recognizes shape of a local minimum of compound 1 , the syn‐Φ/syn‐Ψ conformer, not the structure at global minimum. MD simulations were run to explain, in structural terms, the observed geometry of the complex.     

Conformational Dynamics of apo‐GlnBP Revealed by Experimental and Computational Analysis

The glutamine binding protein (GlnBP) binds l ‐glutamine and cooperates with its cognate transporters during glutamine uptake. Crystal structure analysis has revealed an open and a closed conformation for apo‐ and holo‐GlnBP, respectively. However, the detailed conformational dynamics have remained unclear. Herein, we combined NMR spectroscopy, MD simulations, and single‐molecule FRET techniques to decipher the conformational dynamics of apo‐GlnBP. The NMR residual dipolar couplings of apo‐GlnBP were in good agreement with a MD‐derived structure ensemble consisting of four metastable states. The open and closed conformations are the two major states. This four‐state model was further validated by smFRET experiments and suggests the conformational selection mechanism in ligand recognition of GlnBP.     

Localized and Collective Motions in HET‐s(218‐289) Fibrils from Combined NMR Relaxation and MD Simulation

Nuclear magnetic resonance (NMR) relaxation data and molecular dynamics (MD) simulations are combined to characterize the dynamics of the fungal prion HET‐s(218‐289) in its amyloid form. NMR data is analyzed with the dynamics detector method, which yields timescale‐specific information. An analogous analysis is performed on MD trajectories. Because specific MD predictions can be verified as agreeing with the NMR data, MD was used for further interpretation of NMR results: for the different timescales, cross‐correlation coefficients were derived to quantify the correlation of the motion between different residues. Short timescales are the result of very local motions, while longer timescales are found for longer‐range correlated motion. Similar trends on ns‐ and μs‐timescales suggest that μs motion in fibrils is the result of motion correlated over many fibril layers.     

A comparison of different initialization protocols to obtain statistically independent molecular dynamics simulations

We study how the results of molecular dynamics (MD) simulations are affected by various choices during the setup, e.g., the starting velocities, the solvation, the location of protons, the conformation of His, Asn, and Gln residues, the protonation and titration of His residues, and the treatment of alternative conformations. We estimate the binding affinity of ligands to four proteins calculated with the MM/GBSA method (molecular mechanics combined with a generalized Born and surface area solvation energy). For avidin and T4 lysozyme, all variations gave similar results within 2 kJ/mol. For factor Xa, differences in the solvation or in the selection of alternative conformations gave results that are significantly different from those of the other approaches by 4-6 kJ/mol, whereas for galectin-3, changes in the conformations, rotations, and protonation gave results that differed by 10 kJ/mol, but only if residues close to the binding site were modified. This shows that the results of MM/GBSA calculations are reasonably reproducible even if the MD simulations are set up with different software. Moreover, we show that the sampling of phase space can be enhanced by solvating the systems with different equilibrated water boxes, in addition to the common use of different starting velocities. If different conformations are available in the crystal structure, they should also be employed to enhance the sampling. Protonation, ionization, and conformations of Asn, Gln, and His may also be used to enhance sampling, but great effort should be spent to obtain as reliable predictions as possible close to the active site.     

Backbone motions of free and pheromone-bound major urinary protein I studied by molecular dynamics simulation

Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.     

The Plasticity of the Carbohydrate Recognition Domain Dictates the Exquisite Mechanism of Binding of Human Macrophage Galactose-Type Lectin

The human macrophage galactose-type lectin (MGL), expressed on macrophages and dendritic cells (DCs), modulates distinct immune cell responses by recognizing N-acetylgalactosamine (GalNAc) containing structures present on pathogens, self-glycoproteins, and tumor cells. Herein, NMR spectroscopy and molecular dynamics (MD) simulations were used to investigate the structural preferences of MGL against different GalNAc-containing structures derived from the blood group A antigen, the Forssman antigen, and the GM2 glycolipid. NMR spectroscopic analysis of the MGL carbohydrate recognition domain (MGL-CRD, C181-H316) in the absence and presence of methyl α-GalNAc (α-MeGalNAc), a simple monosaccharide, shows that the MGL-CRD is highly dynamic and its structure is strongly altered upon ligand binding. This plasticity of the MGL-CRD structure explains the ability of MGL to accommodate different GalNAc-containing molecules. However, key differences are observed in the recognition process depending on whether the GalNAc is part of the blood group A antigen, the Forssman antigen, or GM2-derived structures. These results are in accordance with molecular dynamics simulations that suggest the existence of a distinct MGL binding mechanism depending on the context of GalNAc moiety presentation. These results afford new perspectives for the rational design of GalNAc modifications that fine tune MGL immune responses in distinct biological contexts, especially in malignancy.     

Comprehensive Characterization of the Self-Folding Cavitand Dynamics

The conformational equilibria and guest exchange process of a resorcin[4]arene derived self-folding cavitand receptor have been characterized in detail by molecular dynamics simulations (MD) and ¹H EXSY NMR experiments. A multi-timescale strategy for exploring the fluxional behaviour of this system has been constructed, exploiting conventional MD and accelerated MD (aMD) techniques. The use of aMD allows the reconstruction of the folding/unfolding process of the receptor by sampling high-energy barrier processes unattainable by conventional MD simulations. We obtained MD trajectories sampling events occurring at different timescales from ns to s: 1) rearrangement of the directional hydrogen bond seam stabilizing the receptor, 2) folding/unfolding of the structure transiting partially open intermediates, and 3) guest departure from different folding stages. Most remarkably, reweighing of the biased aMD simulations provided kinetic barriers that are in very good agreement with those determined experimentally by ¹H NMR. These results constitute the first comprehensive characterization of the complex dynamic features of cavitand receptors. Our approach emerges as a valuable rational design tool for synthetic host-guest systems     

Interaction of a tripeptide with cesium perfluorooctanoate micelles

The interaction of alanyl-phenylalanyl-alanine (Ala-Phe-Ala) with the micelles formed by cesium perfluorooctanoate (CsPFO) in water was studied in the isotropic phase by means of 1H NMR and by molecular dynamics (MD) simulations. Information on the location of the peptide was experimentally obtained from selective variations in Ala-Phe-Ala chemical shifts and from differential line broadening in the presence of the paramagnetic ion Mn2+. The peptide-micelle association constant was estimated analyzing the chemical shift variations of the most sensitive Ala-Phe-Ala resonances with the peptide concentration. MD simulations of Ala-Phe-Ala in the micellar environment confirmed the experimental observations, identifying the hydrogen bonding interactions of the different peptide moieties with the micelle, yielding a binding constant close to the experimental one. NOESY experiments suggest that the peptide in the micellar environment does not adopt a preferred conformation but is mainly unstructured. Details on the conformational behavior of the peptide in the micellar solution observed through MD were consistent with a different conformational equilibrium in the proximity of the micelle. Information on Ala-Phe-Ala dynamics was obtained from 1H T1 data and compared to MD simulation results on the overall tumbling motion.     

Synthesis and antimicrobial bioassays of 1,3,4-thiadiazole sulfone derivatives containing amide moiety: A study based on molecular dynamics (MD) simulations,MM/GBSA,and molecular docking

The backbone structure (1,3,4-thiadiazole sulfone derivatives containing amide moiety) of target compounds was determined by modification and optimization of the theoretical design based on commercial chemical carboxin, including molecular docking, scaffold hopping, ligand expansion, etc.In this paper, 23 target compounds were synthesized by the combination of theoretical design and chemical synthesis, and characterized by ¹H NMR, ¹³C NMR and HR MS. Addtionally, the antibacterical bioassay showed that most target compounds performed excellent inhibition on Xanthomonas axonopodis pv. citri (Xac) and Xanthomonas oryzae pv. oryzae (Xoo) in vitro. Meanwhile, molecular docking, molecular dynamics (MD) simulations, and studies on ligand/protein (carboxin/2FBW and 4n/2FBW) complex systems were displayed, and the interaction patterns of ligand/protein complex system were predicted by molecular docking. Besides, the ligand/protein complex system was subject to MD simulation. The analysis of molecular dynamics such as RMSD values suggested that compound/2FBW complexes were stable. MM/GBSA (Molecular mechanics generalized born surface area) dynamic binding affinity results revealed that the active residues (TYR58, HIS26, ARG43, SER39, etc.) played an essential part in the binding of the compound(s) to form a stable low-energy ligand/protein complex, while the MD trajectories demonstrated that the interactions of drugs with 2FBW affected the tertiary structure and increased the stability of protein. Besides, compound 4n also showed control efficacies (curative and protective) on Xoo in vivo, where the curative efficacy was 35.91% and the protective efficacy was 18.97%. In a word, this study showed that 1,3,4-thiadiazole sulfone derivatives containing amide moiety designed based on the structure of carboxin were promising agricultural antibacterial agents, featuring certain stability of binding affinity to proteins and carboxin.     

Improvement of Antibody Activity by Controlling Its Dynamics Using the Glycan–Lectin Interaction

Antibody dynamics on membranes, such as endocytosis and clustering, are vital in determining antibody functions. In this study, we demonstrated that glycan conjugation can modulate antibody dynamics through the glycan–lectin interaction to regulate its potency. The anti-HER2 antibody, an anti-breast-cancer antibody, was conjugated with galactose-containing N-glycan, and its internalization was suppressed by interaction with galectin-3, leading to enhanced complement-dependent cytotoxic (CDC) activity. This glycan–antibody conjugate is proposed as a new approach to modulate antibody activity and may provide an alternative strategy for redeveloping antibody drugs that do not exhibit sufficient activity.     

Binding structures of tri‐N‐acetyl‐β‐glucosamine in hen egg white lysozyme using molecular dynamics with a polarizable force field

Lysozyme is a well‐studied enzyme that hydrolyzes the β‐(1,4)‐glycosidic linkage of N‐acetyl‐β‐glucosamine (NAG)_n oligomers. The active site of hen egg‐white lysozyme (HEWL) is believed to consist of six subsites, A‐F that can accommodate six sugar residues. We present studies exploring the use of polarizable force fields in conjunction with all‐atom molecular dynamics (MD) simulations to analyze binding structures of complexes of lysozyme and NAG trisaccharide, (NAG)₃. MD trajectories are applied to analyze structures and conformation of the complex as well as protein–ligand interactions, including the hydrogen‐bonding network in the binding pocket. Two binding modes (ABC and BCD) of (NAG)₃ are investigated independently based on a fixed‐charge model and a polarizable model. We also apply molecular mechanics with generalized born and surface area (MM‐GBSA) methods based on MD using both nonpolarizable and polarizable force fields to compute binding free energies. We also study the correlation between root‐mean‐squared deviation and binding free energies of the wildtype and W62Y mutant; we find that for this prototypical system, approaches using the MD trajectories coupled with implicit solvent models are equivalent for polarizable and fixed‐charge models. © 2012 Wiley Periodicals, Inc.     

Validation of the GROMOS 54A7 Force Field Regarding Mixed α/β‐Peptide Molecules

A molecular‐dynamics (MD) simulation study of two heptapeptides containing α‐ and β‐amino acid residues is presented. According to NMR experiments, the two peptides differ in dominant fold when solvated in MeOH: peptide 3 adopts predominantly β‐hairpin‐like conformations, while peptide 8 adopts a 14/15‐helical fold. The MD simulations largely reproduce the experimental data. Application of NOE atom? atom distance restraining improves the agreement with experimental data, but reduces the conformational sampling. Peptide 3 shows a variety of conformations, while still agreeing with the NOE and ³J‐coupling data, whereas the conformational ensemble of peptide 8 is dominated by one helical conformation. The results confirm the suitability of the GROMOS 54A7 force field for simulation or structure refinement of mixed α/β‐peptides in MeOH.     

<Emphasis Type="Italic">In-silico</Emphasis> Screening using Flexible Ligand Binding Pockets: A Molecular Dynamics-based Approach

Summary In-silico screening of flexible ligands against flexible ligand binding pockets (LBP) is an emerging approach in structure-based drug discovery. Here, we describe a molecular dynamics (MD) based docking approach to investigate the influence on the high-throughput in-silico screening of small molecules against flexible ligand binding pockets. In our approach, an ensemble of 51 energetically favorable structures of the LBP of human estrogen receptor α (hERα) were collected from 3 ns MD simulations. In-silico screening of 3500 endocrine disrupting compounds against these flexible ligand binding pockets resulted in thousands of ER–ligand complexes of which 582 compounds were unique. Detailed analysis of MD generated structures showed that only 17 of the LBP residues significantly contribute to the overall binding pocket flexibility. Using the flexible LBP conformations generated, we have identified 32 compounds that bind better to the flexible ligand-binding pockets compared to the crystal structure. These compounds, though chemically divergent, are structurally similar to the natural hormone. Our MD-based approach in conjunction with grid–based distributed computing could be applied routinely for in-silico screening of large databases against any given target.