Complementarity of hydrophobic properties in ATP-protein binding: a new criterion to rank docking solutions

Recognition of quinolone antibiotics by the multidrug efflux transporter MexB of Pseudomonas aeruginosa

The drug/proton antiporter MexB is the engine of the major efflux pump MexAB-OprM in Pseudomonas aeruginosa. This protein is known to transport a large variety of compounds, including antibiotics, thus conferring a multi-drug resistance phenotype. Due to the difficulty of producing co-crystals, only two X-ray structures of MexB in a complex with ligands are available to date, and mechanistic aspects are largely hypothesized based on the body of data collected for the homologous protein AcrB of Escherichia coli. In particular, a recent study (Ornik-Cha, Wilhelm, Kobylka et al., Nat. Commun., 2021, 12, 6919) reported a co-crystal structure of AcrB in a complex with levofloxacin, an antibiotic belonging to the important class of (fluoro)-quinolones. In this work, we performed a systematic ensemble docking campaign coupled to the cluster analysis and molecular-mechanics optimization of docking poses to study the interaction between 36 quinolone antibiotics and MexB. We additionally investigated surface complementarity between each molecule and the transporter and thoroughly assessed the computational protocol adopted against the known experimental data. Our study reveals different binding preferences of the investigated compounds towards the sub-sites of the large deep binding pocket of MexB, supporting the hypothesis that MexB substrates oscillate between different binding modes with similar affinity. Interestingly, small changes in the molecular structure translate into significant differences in MexB–quinolone interactions. All the predicted binding modes are available for download and visualization at the following link: https://www.dsf.unica.it/dock/mexb/quinolones.

Putative binding modes (BMs) of quinolones to the bacterial efflux transporter MexB were identified. Multiple interaction patterns are possible, supporting the hypothesis that substrates oscillate between different BMs with similar affinity.

Small Amphiphilic Peptides: Activity Against a Broad Range of Drug-Resistant Bacteria and Structural Insight into Membranolytic Properties

We report the synthesis and antibacterial activities of a series of amphiphilic membrane-active peptides composed, in part, of various nongenetically coded hydrophobic amino acids. The lead cyclic peptides, 8C and 9C, showed broad-spectrum activity against drug-resistant Gram-positive (minimum inhibitory concentration (MIC) = 1.5-6.2 μg/mL) and Gram-negative (MIC = 12.5-25 μg/mL) bacteria. The cytotoxicity study showed the predominant lethal action of the peptides against bacteria as compared with mammalian cells. A plasma stability study revealed approximately 2-fold higher stability of lead cyclic peptides as compared to their linear counterparts after 24 h of incubation. A calcein dye leakage experiment revealed the membranolytic effect of the cyclic peptides. Nuclear magnetic resonance spectroscopy and molecular dynamics simulation studies of the interaction of the peptides with the phospholipid bilayer provided a solid structural basis to explain the membranolytic action of the peptides with atomistic details. These results highlight the potential of newly designed amphiphilic peptides as the next generation of peptide-based antibiotics.

A molecular modelling approach to understand the effect of co-evolutionary mutations (V344M, I354L) identified in the PB2 subunit of influenza A 2009 pandemic H1N1 virus on m7GTP ligand binding

The cap binding domain of the polymerase basic 2 (PB2) subunit of influenza polymerases plays a critical role in mediating the 'cap-snatching' mechanism by binding the 5' cap of host pre-mRNAs during viral mRNA transcription. Monitoring variations in the PB2 protein is thus vital for evaluating the pathogenic potential of the virus. Based on selection pressure analysis of PB2 gene sequences of the pandemic H1N1 (pH1N1) viruses of the period 2009-2014, we identified a site, 344V/M, in the vicinity of the cap binding pocket showing evidence of adaptive evolution and another co-evolving residue, 354I/L, in close vicinity. Modelling of the three-dimensional structure of the pH1N1 PB2 cap binding domain, docking of the pre-mRNA cap analogue m7GTP and molecular dynamics simulation studies of the docked complexes performed for four PB2 variants observed showed that the complex possessing V344M with I354L possessed better ligand binding affinity due to additional hydrogen bond contacts between m7GTP and the key residues His432 and Arg355 that was attributed to a displacement of the 424 loop and a flip of the side chain of Arg355, respectively. The co-evolutionary mutations identified (V344M, I354L) were found to be established in the PB2 gene of the pH1N1 viral population over the period 2010-2014. The study demonstrates the molecular basis for the enhanced m7GTP ligand binding affinity with the 344M-354L synergistic combination in PB2. Furthermore, the insight gained into understanding the molecular mechanism of cap binding in pH1N1 viruses may be useful for designing novel drugs targeting the PB2 cap binding domain.

Mutations in adenine-binding pockets enhance catalytic properties of NAD(P)H-dependent enzymes

NAD(P)H-dependent enzymes are ubiquitous in metabolism and cellular processes and are also of great interest for pharmaceutical and industrial applications. Here, we present a structure-guided enzyme engineering strategy for improving catalytic properties of NAD(P)H-dependent enzymes toward native or native-like reactions using mutations to the enzyme's adenine-binding pocket, distal to the site of catalysis. Screening single-site saturation mutagenesis libraries identified mutations that increased catalytic efficiency up to 10-fold in 7 out of 10 enzymes. The enzymes improved in this study represent three different cofactor-binding folds (Rossmann, DHQS-like, and FAD/NAD binding) and utilize both NADH and NADPH. Structural and biochemical analyses show that the improved activities are accompanied by minimal changes in other properties (cooperativity, thermostability, pH optimum, uncoupling), and initial tests on two enzymes (ScADH6 and EcFucO) show improved functionality in Escherichia coli.

Structural Insight into Specificity of Interactions between Nonconventional Three-finger Weak Toxin from Naja kaouthia (WTX) and Muscarinic Acetylcholine Receptors

Weak toxin from Naja kaouthia (WTX) belongs to the group of nonconventional "three-finger" snake neurotoxins. It irreversibly inhibits nicotinic acetylcholine receptors and allosterically interacts with muscarinic acetylcholine receptors (mAChRs). Using site-directed mutagenesis, NMR spectroscopy, and computer modeling, we investigated the recombinant mutant WTX analogue (rWTX) which, compared with the native toxin, has an additional N-terminal methionine residue. In comparison with the wild-type toxin, rWTX demonstrated an altered pharmacological profile, decreased binding of orthosteric antagonist N-methylscopolamine to human M1- and M2-mAChRs, and increased antagonist binding to M3-mAChR. Positively charged arginine residues located in the flexible loop II were found to be crucial for rWTX interactions with all types of mAChR. Computer modeling suggested that the rWTX loop II protrudes to the M1-mAChR allosteric ligand-binding site blocking the entrance to the orthosteric site. In contrast, toxin interacts with M3-mAChR by loop II without penetration into the allosteric site. Data obtained provide new structural insight into the target-specific allosteric regulation of mAChRs by "three-finger" snake neurotoxins.

Mg²⁺ coordinating dynamics in Mg:ATP fueled motor proteins

The coordination of Mg(2+) with the triphosphate group of adenosine triphosphate (ATP) in motor proteins is investigated using data mining and molecular dynamics. The possible coordination structures available from crystal data for actin, myosin, RNA polymerase, DNA polymerase, DNA helicase, and F1-ATPase are verified and investigated further by molecular dynamics. Coordination states are evaluated using structural analysis and quantified by radial distribution functions, coordination numbers, and pair interaction energy calculations. The results reveal a diverse range of both transitory and stable coordination arrangements between Mg(2+) and ATP. The two most stable coordinating states occur when Mg(2+) coordinates two or three oxygens from the triphosphate group of ATP. Evidence for five-site coordination is also reported involving water in addition to the triphosphate group. The stable states correspond to a pair interaction energy of either ∼-2750 kJ/mol or -3500 kJ/mol. The role of water molecules in the hydration shell surrounding Mg(2+) is also reported.

Deciphering fine molecular details of proteins' structure and function with a Protein Surface Topography (PST) method

Molecular surfaces are the key players in biomolecular recognition and interactions. Nowadays, it is trivial to visualize a molecular surface and surface-distributed properties in three-dimensional space. However, such a representation trends to be biased and ambiguous in case of thorough analysis. We present a new method to create 2D spherical projection maps of entire protein surfaces and manipulate with them--protein surface topography (PST). It permits visualization and thoughtful analysis of surface properties. PST helps to easily portray conformational transitions, analyze proteins' properties and their dynamic behavior, improve docking performance, and reveal common patterns and dissimilarities in molecular surfaces of related bioactive peptides. This paper describes basic usage of PST with an example of small G-proteins conformational transitions, mapping of caspase-1 intersubunit interface, and intrinsic "complementarity" in the conotoxin-acetylcholine binding protein complex. We suggest that PST is a beneficial approach for structure-function studies of bioactive peptides and small proteins.

Conformational dynamics of ATP/Mg:ATP in motor proteins via data mining and molecular simulation

The conformational diversity of ATP/Mg:ATP in motor proteins was investigated using molecular dynamics and data mining. Adenosine triphosphate (ATP) conformations were found to be constrained mostly by inter cavity motifs in the motor proteins. It is demonstrated that ATP favors extended conformations in the tight pockets of motor proteins such as F(1)-ATPase and actin whereas compact structures are favored in motor proteins such as RNA polymerase and DNA helicase. The incorporation of Mg(2+) leads to increased flexibility of ATP molecules. The differences in the conformational dynamics of ATP/Mg:ATP in various motor proteins was quantified by the radius of gyration. The relationship between the simulation results and those obtained by data mining of motor proteins available in the protein data bank is analyzed. The data mining analysis of motor proteins supports the conformational diversity of the phosphate group of ATP obtained computationally.

ATP changes the fluorescence lifetime of cyan fluorescent protein via an interaction with His148

Recently, we described that ATP induces changes in YFP/CFP fluorescence intensities of Fluorescence Resonance Energy Transfer (FRET) sensors based on CFP-YFP. To get insight into this phenomenon, we employed fluorescence lifetime spectroscopy to analyze the influence of ATP on these fluorescent proteins in more detail. Using different donor and acceptor pairs we found that ATP only affected the CFP-YFP based versions. Subsequent analysis of purified monomers of the used proteins showed that ATP has a direct effect on the fluorescence lifetime properties of CFP. Since the fluorescence lifetime analysis of CFP is rather complicated by the existence of different lifetimes, we tested a variant of CFP, i.e. Cerulean, as a monomer and in our FRET constructs. Surprisingly, this CFP variant shows no ATP concentration dependent changes in the fluorescence lifetime. The most important difference between CFP and Cerulean is a histidine residue at position 148. Indeed, changing this histidine in CFP into an aspartic acid results in identical fluorescence properties as observed for the Cerulean fluorescent based FRET sensor. We therefore conclude that the changes in fluorescence lifetime of CFP are affected specifically by possible electrostatic interactions of the negative charge of ATP with the positively charged histidine at position 148. Clearly, further physicochemical characterization is needed to explain the sensitivity of CFP fluorescence properties to changes in environmental (i.e. ATP concentrations) conditions.

[Molecular docking: role of intermolecular contacts in formation of complexes of proteins with nucleotides and peptides]

Knowledge of 3D-structure of protein-ligand complex is a major prerequisite for understanding the functioning mechanism of cellular proteins and membrane receptors. This is also of a great help in rational drug design projects. In the present paper we briefly review the molecular docking approaches used to predict possible orientation of a ligand in the protein binding site. The recent trends to improve the accuracy and efficiency of docking algorithms are demonstrated with the results obtained in Laboratory of Biomolecular Modeling. Particular attention is paid to protein-ligand hydrophobic and stacking interactions responsible for molecular recognition of ligand fragments. Such type of interactions are not always adequately represented in scoring criteria of docking applications that leads to mismatch in 3D-structure complexes predictions. That is why further inquiry of methods to account for these interactions is now the area of active research.

Saccharomyces cerevisiae phosphoenolpyruvate carboxykinase: the relevance of Glu299 and Leu460 for nucleotide binding

A homology model of Saccharomyces cerevisiae phosphoenolpyruvate (PEP) carboxykinase (ATP + oxaloacetate right arrow over left arrow ADP + PEP + CO(2)) in complex with its substrates shows that the isobutyl group of Leu460 is in close proximity to the adenine ring of the nucleotide, while the carboxyl group of Glu299 is within hydrogen-bonding distance of the ribose 2'OH. The Leu460Ala mutation caused three-fold and seven-fold increases in the K (m) for ADPMn(-) and ATPMn(2-), respectively, while the Glu299Ala mutation had no effect. Binding studies showed losses of approximately 2 kcal mol(-1) in the nucleotide binding affinity due to the Leu460Ala mutation and no effect for the Glu299Ala mutation. PEP carboxykinase utilized 2'deoxyADP and 2'deoxyATP as substrates with kinetic and equilibrium dissociation constants very similar to those of ADP and ATP, respectively. These results show that the hydrophobic interaction between Leu460 and the adenine ring of the nucleotide significantly contributed to the nucleotide affinity of the enzyme. The 2'deoxy nucleotide studies and the lack of an effect of the Glu299Ala mutation in nucleotide binding suggest that the possible hydrogen bond contributed by Glu299 and the ribose 2'OH group may not be relevant for nucleotide binding.

Three-dimensional structures reveal multiple ADP/ATP binding modes for a synthetic class of artificial proteins

The creation of synthetic enzymes with predefined functions represents a major challenge in future synthetic biology applications. Here, we describe six structures of de novo proteins that have been determined using protein crystallography to address how simple enzymes perform catalysis. Three structures are of a protein, DX, selected for its stability and ability to tightly bind ATP. Despite the addition of ATP to the crystallization conditions, the presence of a bound but distorted ATP was found only under excess ATP conditions, with ADP being present under equimolar conditions or when crystallized for a prolonged period of time. A bound ADP cofactor was evident when Asp was substituted for Val at residue 65, but ATP in a linear configuration is present when Phe was substituted for Tyr at residue 43. These new structures complement previously determined structures of DX and the protein with the Phe 43 to Tyr substitution [Simmons, C. R., et al. (2009) ACS Chem. Biol. 4, 649-658] and together demonstrate the multiple ADP/ATP binding modes from which a model emerges in which the DX protein binds ATP in a configuration that represents a transitional state for the catalysis of ATP to ADP through a slow, metal-free reaction capable of multiple turnovers. This unusual observation suggests that design-free methods can be used to generate novel protein scaffolds that are tailor-made for catalysis.

Relevance of Arg457 for the nucleotide affinity of Saccharomyces cerevisiae phosphoenolpyruvate carboxykinase

Phosphoenolpyruvate carboxykinases catalyze one of the first steps in the biosynthesis of glucose and depending on the enzyme origin, preferentially use adenine or guanine nucleotides as substrates. The Saccharomyces cerevisiae enzyme has a marked preference for ADP (or ATP) over other nucleotides. Homology models of the enzyme in complex with ADP or ATP show that the guanidinium group of Arg457 is close to the adenine base, suggesting that this group might be involved in the stabilization of the nucleotide substrate. To evaluate this we have performed the mutation Arg457Met, replacing the positively charged guanidinium group by a neutral residue. The mutated enzyme retained the structural characteristics of the wild-type protein. Fluorescence titration experiments showed that mutation causes a loss of 1.7 kcal mol(-1) in the binding affinity of the enzyme for ADPMn. Similarly, kinetic analyses of the mutated enzyme showed 50-fold increase in K(m) for ADPMn, with minor alterations in the other kinetic parameters. These results show that Arg457 is an important factor for nucleotide binding by S. cerevisiae PEP carboxykinase.