Improvements in Precision of Relative Binding Free Energy Calculations Afforded by the Alchemical Enhanced Sampling (ACES) Approach.

Accurate predictions of how strongly small molecules bind to proteins, such as those afforded by relative binding free energy (RBFE) calculations, can greatly increase the efficiency of the hit-to-lead and lead optimization stages of the drug discovery process. The success of such calculations, however, relies heavily on their precision. Here, we show that a recently developed alchemical enhanced sampling (ACES) approach can consistently improve the precision of RBFE calculations on a large and diverse set of proteins and small molecule ligands.

ACES: Optimized Alchemically Enhanced Sampling.

We present an alchemical enhanced sampling (ACES) method implemented in the GPU-accelerated AMBER free energy MD engine. The methods hinges on the creation of an "enhanced sampling state" by reducing or eliminating selected potential energy terms and interactions that lead to kinetic traps and conformational barriers while maintaining those terms that curtail the need to otherwise sample large volumes of phase space. For example, the enhanced sampling state might involve transforming regions of a ligand and/or protein side chain into a noninteracting "dummy state" with internal electrostatics and torsion angle terms turned off.

AMBER Free Energy Tools: A New Framework for the Design of Optimized Alchemical Transformation Pathways.

We develop a framework for the design of optimized alchemical transformation pathways in free energy simulations using nonlinear mixing and a new functional form for so-called "softcore" potentials. We describe the implementation and testing of this framework in the GPU-accelerated AMBER software suite. The new optimized alchemical transformation pathways integrate a number of important features, including (1) the use of smoothstep functions to stabilize behavior near the transformation end points, (2) consistent power scaling of Coulomb and Lennard-Jones (LJ) interactions with unitless control parameters to maintain balance of electrostatic attractions and exchange repulsions, (3) pairwise form based on the LJ contact radius for the effective interaction distance with separation-shifted scaling, and (4) rigorous smoothing of the potential at the nonbonded cutoff boundary.

AMBER Drug Discovery Boost Tools: Automated Workflow for Production Free-Energy Simulation Setup and Analysis (ProFESSA).

We report an automated workflow for production free-energy simulation setup and analysis (ProFESSA) using the GPU-accelerated AMBER free-energy engine with enhanced sampling features and analysis tools, part of the AMBER Drug Discovery Boost package that has been integrated into the AMBER22 release. The workflow establishes a flexible, end-to-end pipeline for performing alchemical free-energy simulations that brings to bear technologies, including new enhanced sampling features and analysis tools, to practical drug discovery problems. ProFESSA provides the user with top-level control of large sets of free-energy calculations and offers access to the following key functionalities: (1) automated setup of file infrastructure; (2) enhanced conformational and alchemical sampling with the ACES method; and (3) network-wide free-energy analysis with the optional imposition of cycle closure and experimental constraints.

Room-temperature serial synchrotron crystallography of Drosophila cryptochrome.

Fixed-target serial crystallography allows the high-throughput collection of diffraction data from small crystals at room temperature. This methodology is particularly useful for difficult samples that have sensitivity to radiation damage or intolerance to cryoprotection measures; fixed-target methods also have the added benefit of low sample consumption. Here, this method is applied to the structure determination of the circadian photoreceptor cryptochrome (CRY), previous structures of which have been determined at cryogenic temperature.

Introducing a New Bond-Forming Activity in an Archaeal DNA Polymerase by Structure-Guided Enzyme Redesign.

DNA polymerases have evolved to feature a highly conserved activity across the tree of life: formation of, without exception, internucleotidyl O-P linkages. Can this linkage selectivity be overcome by design to produce xenonucleic acids? Here, we report that the structure-guided redesign of an archaeal DNA polymerase, 9°N, exhibits a new activity undetectable in the wild-type enzyme: catalyzing the formation of internucleotidyl N-P linkages using 3'-NH-ddNTPs. Replacing a metal-binding aspartate in the 9°N active site with asparagine was key to the emergence of this unnatural enzyme activity.

Mechanistic insight into light-dependent recognition of Timeless by Drosophila Cryptochrome.

Cryptochrome (CRY) entrains the fly circadian clock by binding to Timeless (TIM) in light. Undocking of a helical C-terminal tail (CTT) in response to photoreduction of the CRY flavin cofactor gates TIM recognition. We present a generally applicable select western-blot-free tagged-protein interaction (SWFTI) assay that allowed the quantification of CRY binding to TIM in dark and light.

Peripheral Methionine Residues Impact Flavin Photoreduction and Protonation in an Engineered LOV Domain Light Sensor.

Proton-coupled electron transfer reactions play critical roles in many aspects of sensory phototransduction. In the case of flavoprotein light sensors, reductive quenching of flavin excited states initiates chemical and conformational changes that ultimately transmit light signals to downstream targets. These reactions generally require neighboring aromatic residues and proton-donating side chains for rapid and coordinated electron and proton transfer to flavin.

Tuning flavin environment to detect and control light-induced conformational switching in Drosophila cryptochrome.

Light-induction of an anionic semiquinone (SQ) flavin radical in Drosophila cryptochrome (dCRY) alters the dCRY conformation to promote binding and degradation of the circadian clock protein Timeless (TIM). Specific peptide ligation with sortase A attaches a nitroxide spin-probe to the dCRY C-terminal tail (CTT) while avoiding deleterious side reactions. Pulse dipolar electron-spin resonance spectroscopy from the CTT nitroxide to the SQ shows that flavin photoreduction shifts the CTT ~1 nm and increases its motion, without causing full displacement from the protein.

Confluence of theory and experiment reveals the catalytic mechanism of the Varkud satellite ribozyme.

The Varkud satellite ribozyme catalyses site-specific RNA cleavage and ligation, and serves as an important model system to understand RNA catalysis. Here, we combine stereospecific phosphorothioate substitution, precision nucleobase mutation and linear free-energy relationship measurements with molecular dynamics, molecular solvation theory and ab initio quantum mechanical/molecular mechanical free-energy simulations to gain insight into the catalysis. Through this confluence of theory and experiment, we unify the existing body of structural and functional data to unveil the catalytic mechanism in unprecedented detail, including the degree of proton transfer in the transition state.

Evidence for a Catalytic Strategy to Promote Nucleophile Activation in Metal-Dependent RNA-Cleaving Ribozymes and 8-17 DNAzyme.

An unique catalytic strategy was recently reported for the S ribozyme [Bingaman et al., Nat. Chem.

An Ontology for Facilitating Discussion of Catalytic Strategies of RNA-Cleaving Enzymes.

A predictive understanding of the mechanisms of RNA cleavage is important for the design of emerging technology built from biological and synthetic molecules that have promise for new biochemical and medicinal applications. Over the past 15 years, RNA cleavage reactions involving 2'-O-transphosphorylation have been discussed using a simplified framework introduced by Breaker that consists of four fundamental catalytic strategies (designated α, β, γ, and δ) that contribute to rate enhancement. As more detailed mechanistic data emerge, there is need for the framework to evolve and keep pace.

Elucidation of the Catalytic Mechanism of a Miniature Zinc Finger Hydrolase.

To improve our mechanistic understanding of zinc metalloenzymes, we report a joint computational and experimental study of a minimal carbonic anhydrase (CA) mimic, a 22-residue Zn-finger hydrolase. We combine classical molecular dynamics (MD) simulations, quantum mechanics/molecular mechanics (QM/MM) geometry optimizations, and QM/MM free energy simulations with ambient and high-pressure kinetic measurements to investigate the mechanism of the hydrolysis of the substrate p-nitrophenylacetate (pNPA). The zinc center of the hydrolase prefers a pentacoordinated geometry, as found in most naturally occurring CAs and CA-like enzymes.

Importance of MM Polarization in QM/MM Studies of Enzymatic Reactions: Assessment of the QM/MM Drude Oscillator Model.

Unlabelled: For accurate quantum mechanics/molecular mechanics (QM/MM) studies of enzymatic reactions, it is desirable to include MM polarization, for example by using the Drude oscillator (DO) model. For a long time, such studies were hampered by the lack of well-tested polarizable force fields for proteins. Following up on a recent preliminary QM/MM-DO assessment (J.

Glutamine Amide Flip Elicits Long Distance Allosteric Responses in the LOV Protein Vivid.

Light-oxygen-voltage (LOV) domains sense blue light through the photochemical formation of a cysteinyl-flavin covalent adduct. Concurrent protonation at the flavin N5 position alters the hydrogen bonding interactions of an invariant Gln residue that has been proposed to flip its amide side chain as a critical step in the propagation of conformational change. Traditional molecular dynamics (MD) and replica-exchange MD (REMD) simulations of the well-characterized LOV protein Vivid (VVD) demonstrate that the Gln182 amide indeed reorients by ∼180° in response to either adduct formation or reduction of the isoalloxazine ring to the neutral semiquinone, both of which involve N5 protonation.

Changes in active site histidine hydrogen bonding trigger cryptochrome activation.

Cryptochrome (CRY) is the principal light sensor of the insect circadian clock. Photoreduction of the Drosophila CRY (dCRY) flavin cofactor to the anionic semiquinone (ASQ) restructures a C-terminal tail helix (CTT) that otherwise inhibits interactions with targets that include the clock protein Timeless (TIM). All-atom molecular dynamics (MD) simulations indicate that flavin reduction destabilizes the CTT, which undergoes large-scale conformational changes (the CTT release) on short (25 ns) timescales.

Inverse thio effects in the hepatitis delta virus ribozyme reveal that the reaction pathway is controlled by metal ion charge density.

The hepatitis delta virus (HDV) ribozyme self-cleaves in the presence of a wide range of monovalent and divalent ions. Prior theoretical studies provided evidence that self-cleavage proceeds via a concerted or stepwise pathway, with the outcome dictated by the valency of the metal ion. In the present study, we measure stereospecific thio effects at the nonbridging oxygens of the scissile phosphate under a wide range of experimental conditions, including varying concentrations of diverse monovalent and divalent ions, and combine these with quantum mechanical/molecular mechanical (QM/MM) free energy simulations on the stereospecific thio substrates.

Role of the active site guanine in the glmS ribozyme self-cleavage mechanism: quantum mechanical/molecular mechanical free energy simulations.

The glmS ribozyme catalyzes a self-cleavage reaction at the phosphodiester bond between residues A-1 and G1. This reaction is thought to occur by an acid-base mechanism involving the glucosamine-6-phosphate cofactor and G40 residue. Herein quantum mechanical/molecular mechanical free energy simulations and pKa calculations, as well as experimental measurements of the rate constant for self-cleavage, are utilized to elucidate the mechanism, particularly the role of G40.

Quantum mechanical/molecular mechanical free energy simulations of the self-cleavage reaction in the hepatitis delta virus ribozyme.

The hepatitis delta virus (HDV) ribozyme catalyzes a self-cleavage reaction using a combination of nucleobase and metal ion catalysis. Both divalent and monovalent ions can catalyze this reaction, although the rate is slower with monovalent ions alone. Herein, we use quantum mechanical/molecular mechanical (QM/MM) free energy simulations to investigate the mechanism of this ribozyme and to elucidate the roles of the catalytic metal ion.

Thio effects and an unconventional metal ion rescue in the genomic hepatitis delta virus ribozyme.

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg(2+) ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies of the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn(2+).

Identification of the catalytic Mg²⁺ ion in the hepatitis delta virus ribozyme.

The hepatitis delta virus ribozyme catalyzes an RNA cleavage reaction using a catalytic nucleobase and a divalent metal ion. The catalytic base, C75, serves as a general acid and has a pK(a) shifted toward neutrality. Less is known about the role of metal ions in the mechanism.

Quantum Mechanical/Molecular Mechanical Study of the HDV Ribozyme: Impact of the Catalytic Metal Ion on the Mechanism.

A recent crystal structure of the precleaved HDV ribozyme along with biochemical data support a mechanism for phosphodiester bond self-cleavage in which C75 acts as a general acid and bound Mg(2+) ion acts as a Lewis acid. Herein this precleaved crystal structure is used as the basis for quantum mechanical/molecular mechanical calculations. These calculations indicate that the self-cleavage reaction is concerted with a phosphorane-like transition state when a divalent ion, Mg(2+) or Ca(2+), is bound at the catalytic site but is sequential with a phosphorane intermediate when a monovalent ion, such as Na(+), is at this site.

Mechanistic strategies in the HDV ribozyme: chelated and diffuse metal ion interactions and active site protonation.

The crystal structure of the precleaved form of the hepatitis delta virus (HDV) ribozyme reveals two G•U wobbles near the active site: a rare reverse G•U wobble involving a syn G base, and a standard G•U wobble at the cleavage site. The catalytic mechanism for this ribozyme has been proposed to involve a Mg(2+) ion bound to the reverse G•U wobble, as well as a protonated C75 base. We carried out molecular dynamics simulations to analyze metal ion interaction with the reverse and standard G•U wobbles and to investigate the impact of C75 protonation on the structure and motions of the ribozyme.

Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions.

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis.

Excess entropy and structural transitions in a two-dimensional square-shoulder fluid.

Metropolis Monte Carlo simulations on the square-shoulder fluid of Malescio and Pellicane are used to trace the temperature dependent excess entropy, the heat capacity, and configurational energy along several isochores, including those for which mechanically stable zero-temperature structures have been predicted. Thermodynamic signatures of structural phase transitions are identified along several isochores, in addition to the low-density triangular solid and stripe phase transitions identified earlier. The finite temperature phases illustrate the competition between cluster formation and stripe formation as competing mechanisms for generating minimum free energy configurations as a function of density, consistent with earlier results at zero temperature.