Water network chemistry to exploit the nature of catalytic water molecules in Mtb DNA gyrase: a computational study to understand the binding mechanism of fluoroquinolones

The dynamics of DNA gyrase and mutants of DNA gyrA such as G88A, A90V, S91P, D94A, D94G, D94N, D94Y; and double-point mutant (S91P-D94G), are meticulously investigated using computational approaches. Molecular dynamics (MD) and hydration thermodynamics have shed light on the fundamental, mechanistic basis of mutations on the conformational stability of Quinolone Binding Pocket (QBP) of DNA gyrase. Analysis of MD results revealed the displacement of a single crystal water molecule (HOH201) from the catalytic site of wild-type (WT) and mutants of DNA gyrA. This prompted our research group to probe the five crystal water molecules present in the QBP of the enzyme using water thermodynamics. Hydration thermodynamics analysis revealed the displacement of HOH201 due to unstable thermodynamic signatures. Further, the analysis highlighted significant changes in thermodynamic signatures and locations of five crystal water hydration sites upon mutation. Integrated MD simulations and water thermodynamics provided promising insights into the conformational changes and inaccessibility of the catalytic water molecule that can influence the design of DNA gyrase inhibitors.Communicated by Ramaswamy H. Sarma.

Water Networks Repopulate Protein-Ligand Interfaces with Temperature

High-resolution crystal structures highlight the importance of water networks in protein-ligand interactions. However, as these are typically determined at cryogenic temperature, resulting insights may be structurally precise but not biologically accurate. By collecting 10 matched room-temperature and cryogenic datasets of the biomedical target Hsp90α, we identified changes in water networks that impact protein conformations at the ligand binding interface. Water repositioning with temperature repopulates protein ensembles and ligand interactions. We introduce Flipper conformational barcodes to identify temperature-sensitive regions in electron density maps. This revealed that temperature-responsive states coincide with ligand-responsive regions and capture unique binding signatures that disappear upon cryo-cooling. Our results have implications for discovering Hsp90 selective ligands, and, more generally, for the utility of hidden protein and water conformations in drug discovery.

Structural dynamics support electrostatic interactions in the active site of Adenylate Kinase

Electrostatic preorganization as well as structural and dynamic heterogeneity are often used to rationalize the remarkable catalytic efficiency of enzymes. However, they are often presented as incompatible because the generation of permanent electrostatic effects implies that the protein structure remains rigid. Here, we use a metric, electric fields, that can treat electrostatic contributions and dynamics effects on equal footing, for a unique perspective on enzymatic catalysis. We find that the residues that contribute the most to electrostatic interactions with the substrate in the active site of Adenylate Kinase (our working example) are also the most flexible residues. Further, entropy-tuning mutations raise flexibility at the picosecond timescale where more conformations can be visited on short time periods, thereby softening the sharp heterogeneity normally visible at the microsecond timescale.

Engineering the ligand specificity of the human galectin-1 by incorporation of tryptophan analogs

Galectin-1 is a β-galactoside-binding lectin with manifold biological functions. A single tryptophan residue (W68) in its carbohydrate binding site plays a major role in ligand binding and is highly conserved among galectins. To fine tune galectin-1 specificity, we introduced several non-canonical tryptophan analogs at this position of human galectin-1 and analyzed the resulting variants using glycan microarrays. Two variants containing 7-azatryptophan and 7-fluorotryptophan showed a reduced affinity for 3'-sulfated oligosaccharides. Their interaction with different ligands was further analyzed by fluorescence polarization competition assay. Using molecular modeling we provide structural clues that the change in affinities comes from modulated interactions and solvation patterns. Thus, we show that the introduction of subtle atomic mutations in the ligand binding site of galectin-1 is an attractive approach for fine-tuning its interactions with different ligands.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.