Jeff Wereszczynski

Statistical mechanics and molecular dynamics in evaluating thermodynamic properties of biomolecular recognition.

Molecular recognition plays a central role in biochemical processes. Although well studied, understanding the mechanisms of recognition is inherently difficult due to the range of potential interactions, the molecular rearrangement associated with binding, and the time and length scales involved. Computational methods have the potential for not only complementing experiments that have been performed, but also in guiding future ones through their predictive abilities. In this review, we discuss how molecular dynamics (MD) simulations may be used in advancing our understanding of the thermodynamics that drive biomolecular recognition. We begin with a brief review of the statistical mechanics that form a basis for these methods. This is followed by a description of some of the most commonly used methods: thermodynamic pathways employing alchemical transformations and potential of mean force calculations, along with end-point calculations for free energy differences, and harmonic and quasi-harmonic analysis for entropic calculations. Finally, a few of the fundamental findings that have resulted from these methods are discussed, such as the role of configurational entropy and solvent in intermolecular interactions, along with selected results of the model system T4 lysozyme to illustrate potential and current limitations of these methods.

Nucleotide-dependent mechanism of Get3 as elucidated from free energy calculations

The unique topology of tail-anchored (TA) proteins precludes them from utilizing the well-studied cotranslational translocation mechanism of most transmembrane proteins, forcing them into a distinct, posttranslational pathway. In yeast, this process is the guided entry of TA-proteins (GET) pathway, which utilizes a combination of cytosolic and transmembrane proteins to identify a TA protein, transfer it, and insert it into the endoplasmic reticulum membrane. At the center of this mechanism is the Get3 homodimer, which transfers a TA protein between the two GET phases by leveraging energy gained in ATP binding and hydrolysis to undergo significant structural changes from “open” to “closed” conformations. We present all-atom molecular dynamics simulations of Get3 in multiple nucleotide states, and through rigorous potential of mean force calculations, compute the free energy landscape of the Get3 opening/closing pathway. Results agree well with experiments on the nucleotide bias of Get3 open and closed structures in the crystallographically observed no-nucleotide, two ATP, and two ADP states, and also reveal their populations in the asymmetric one ATP and one ADP cases. Structures also compare well with the recently observed “semiopen” conformation and suggest that Get3 may sample this state free in solution and not just when bound to Get1, as observed in experiments. Finally, we present evidence for a unique, “wide-open” conformation of Get3. These calculations describe the nucleotide-dependent thermodynamics of Get3 in solution, and improve our understanding of its mechanism in each phase of the GET cycle.

Impact of calcium on N1 influenza neuraminidase dynamics and binding free energy.

The highly pathogenic influenza strains H5N1 and H1N1 are currently treated with inhibitors of the viral surface protein neuraminidase (N1). Crystal structures of N1 indicate a conserved, high affinity calcium binding site located near the active site. The specific role of this calcium in the enzyme mechanism is unknown, though it has been shown to be important for enzymatic activity and thermostability. We report molecular dynamics (MD) simulations of calcium‐bound and calcium‐free N1 complexes with the inhibitor oseltamivir (marketed as the drug Tamiflu), independently using both the AMBER FF99SB and GROMOS96 force fields, to give structural insight into calcium stabilization of key framework residues. Y347, which demonstrates similar sampling patterns in the simulations of both force fields, is implicated as an important N1 residue that can “clamp” the ligand into a favorable binding pose. Free energy perturbation and thermodynamic integration calculations, using two different force fields, support the importance of Y347 and indicate a +3 to +5 kcal/mol change in the binding free energy of oseltamivir in the absence of calcium. With the important role of structure‐based drug design for neuraminidase inhibitors and the growing literature on emerging strains and subtypes, inclusion of this calcium for active site stability is particularly crucial for computational efforts such as homology modeling, virtual screening, and free energy methods. Proteins 2010.

On structural transitions, thermodynamic equilibrium, and the phase diagram of DNA and RNA duplexes under torque and tension

A precise understanding of the flexibility of double stranded nucleic acids and the nature of their deformed conformations induced by external forces is important for a wide range of biological processes including transcriptional regulation, supercoil and catenane removal, and site-specific recombination. We present, at atomic resolution, a simulation of the dynamics involved in the transitions from B-DNA and A-RNA to Pauling (P) forms and to denatured states driven by application of external torque and tension. We then calculate the free energy profile along a B- to P-transition coordinate and from it, compute a reversible pathway, i.e., an isotherm of tension and torque pairs required to maintain P-DNA in equilibrium. The reversible isotherm maps correctly onto a phase diagram derived from single molecule experiments, and yields values of elongation, twist, and twist-stretch coupling in agreement with measured values. We also show that configurational entropy compensates significantly for the large electrostatic energy increase due to closer-packed P backbones. A similar set of simulations applied to RNA are used to predict a novel structure, P-RNA, with its associated free energy, equilibrium tension, torque and structural parameters, and to assign the location, on the phase-diagram, of a putative force–torque-dependent RNA “triple point.”

Discovery of Staphylococcus aureus sortase A inhibitors using virtual screening and the relaxed complex scheme.

Staphylococcus aureus is the leading cause of hospital‐acquired infections in the United States. The emergence of multidrug‐resistant strains of S. aureus has created an urgent need for new antibiotics. Staphylococcus aureus uses the sortase A enzyme to display surface virulence factors suggesting that compounds that inhibit its activity will function as potent anti‐infective agents. Here, we report the identification of several inhibitors of sortase A using virtual screening methods that employ the relaxed complex scheme, an advanced computer‐docking methodology that accounts for protein receptor flexibility. Experimental testing validates that several compounds identified in the screen inhibit the activity of sortase A. A lead compound based on the 2‐phenyl‐2,3‐dihydro‐1H‐perimidine scaffold is particularly promising, and its binding mechanism was further investigated using molecular dynamics simulations and conducting preliminary structure–activity relationship studies.

Simulations of Biased Agonists in the β2 Adrenergic Receptor with Accelerated Molecular Dynamics

The biased agonism of the G protein-coupled receptors (GPCRs), where in addition to a traditional G protein-signaling pathway a GPCR promotes intracellular signals though β-arrestin, is a novel paradigm in pharmacology. Biochemical and biophysical studies have suggested that a GPCR forms a distinct ensemble of conformations signaling through the G protein and β-arrestin. Here we report on the dynamics of the β2 adrenergic receptor bound to the β-arrestin and G protein-biased agonists and the empty receptor to further characterize the receptor conformational changes caused by biased agonists. We use conventional and accelerated molecular dynamics (aMD) simulations to explore the conformational transitions of the GPCR from the active state to the inactive state. We found that aMD simulations enable monitoring of the transition within the nanosecond time scale while capturing the known microscopic characteristics of the inactive states, such as the ionic lock, the inward position of F6.44, and water clusters. Distinct conformational states are shown to be stabilized by each biased agonist. In particular, in simulations of the receptor with the β-arrestin-biased agonist N-cyclopentylbutanepherine, we observe a different pattern of motions in helix 7 when compared to simulations with the G protein-biased agonist salbutamol that involves perturbations of the network of interactions within the NPxxY motif. Understanding the network of interactions induced by biased ligands and the subsequent receptor conformational shifts will lead to development of more efficient drugs.

Structural and Computational Studies of the Staphylococcus aureus Sortase B-Substrate Complex Reveal a Substrate-stabilized Oxyanion Hole.

Background: Sortase enzymes catalyze a transpeptidation reaction that displays bacterial surface proteins. Results: Structural and computational studies reveal how the sortase B enzyme recognizes its sorting signal substrate. Conclusion: Sortase enzymes catalyze transpeptidation using a substrate-stabilized oxyanion hole. Significance: The results of this work could facilitate the rational design of sortase inhibitors. Sortase cysteine transpeptidases covalently attach proteins to the bacterial cell wall or assemble fiber-like pili that promote bacterial adhesion. Members of this enzyme superfamily are widely distributed in Gram-positive bacteria that frequently utilize multiple sortases to elaborate their peptidoglycan. Sortases catalyze transpeptidation using a conserved active site His-Cys-Arg triad that joins a sorting signal located at the C terminus of their protein substrate to an amino nucleophile located on the cell surface. However, despite extensive study, the catalytic mechanism and molecular basis of substrate recognition remains poorly understood. Here we report the crystal structure of the Staphylococcus aureus sortase B enzyme in a covalent complex with an analog of its NPQTN sorting signal substrate, revealing the structural basis through which it displays the IsdC protein involved in heme-iron scavenging from human hemoglobin. The results of computational modeling, molecular dynamics simulations, and targeted amino acid mutagenesis indicate that the backbone amide of Glu224 and the side chain of Arg233 form an oxyanion hole in sortase B that stabilizes high energy tetrahedral catalytic intermediates. Surprisingly, a highly conserved threonine residue within the bound sorting signal substrate facilitates construction of the oxyanion hole by stabilizing the position of the active site arginine residue via hydrogen bonding. Molecular dynamics simulations and primary sequence conservation suggest that the sorting signal-stabilized oxyanion hole is a universal feature of enzymes within the sortase superfamily.

The binding mechanism, multiple binding modes, and allosteric regulation of Staphylococcus aureus Sortase A probed by molecular dynamics simulations

Sortase enzymes are vitally important for the virulence of gram‐positive bacteria as they play a key role in the attachment of surface proteins to the cell wall. These enzymes recognize a specific sorting sequence in proteins destined to be displayed on the surface of the bacteria and catalyze the transpeptidation reaction that links it to a cell wall precursor molecule. Because of their role in establishing pathogenicity, and in light of the recent rise of antibiotic‐resistant bacterial strains, sortase enzymes are novel drug targets. Here, we present a study of the prototypical sortase protein Staphylococcus aureus Sortase A (SrtA). Both conventional and accelerated molecular dynamics simulations of S. aureus SrtA in its apo state and when bound to an LPATG sorting signal (SS) were performed. Results support a binding mechanism that may be characterized as conformational selection followed by induced fit. Additionally, the SS was found to adopt multiple metastable states, thus resolving discrepancies between binding conformations in previously reported experimental structures. Finally, correlation analysis reveals that the SS actively affects allosteric pathways throughout the protein that connect the first and the second substrate binding sites, which are proposed to be located on opposing faces of the protein. Overall, these calculations shed new light on the role of dynamics in the binding mechanism and function of sortase enzymes.

Thermodynamic integration to predict host-guest binding affinities.

An alchemical free energy method with explicit solvent molecular dynamics simulations was applied as part of the blind prediction contest SAMPL3 to calculate binding free energies for seven guests to an acyclic cucurbit-[n]uril host. The predictions included determination of protonation states for both host and guests, docking pose generation, and binding free energy calculations using thermodynamic integration. We found a root mean square error (RMSE) of