A sensitive and scalable fluorescence anisotropy single stranded RNA targeting approach for monitoring riboswitch conformational states.

The capacity of riboswitches to undergo conformational changes in response to binding their native ligands is closely tied to their functional roles and is an attractive target for antimicrobial drug design. Here, we established a probe-based fluorescence anisotropy assay to monitor riboswitch conformational switching with high sensitivity and throughput. Using the Bacillus subtillis yitJ S-Box (SAM-I), Fusobacterium nucleatum impX RFN element of (FMN) and class-I cyclic-di-GMP from Vibrio cholerae riboswitches as model systems, we developed short fluorescent DNA probes that specifically recognize either ligand-free or -bound riboswitch conformational states. We showed that increasing concentrations of native ligands cause measurable and reproducible changes in fluorescence anisotropy that correlate with riboswitch conformational changes observed by native gel analysis. Furthermore, we applied our assay to several ligand analogues and confirmed that it can discriminate between ligands that bind, triggering the native conformational change, from those that bind without causing the conformational change. This new platform opens the possibility of high-throughput screening compound libraries to identify potential new antibiotics that specifically target functional conformational changes in riboswitches.

Docking Ligands into Flexible and Solvated Macromolecules. 8. Forming New Bonds─Challenges and Opportunities.

Over the years, structure-based design programs and specifically docking small molecules to proteins have become prominent in drug discovery. However, many of these computational tools have been developed to primarily dock enzyme inhibitors (and ligands to other protein classes) relying heavily on hydrogen bonds and electrostatic and hydrophobic interactions. In reality, many drug targets either feature metal ions, can be targeted covalently, or are simply not even proteins (e.g., nucleic acids). Herein, we describe several new features that we have implemented into Fitted to broaden its applicability to a wide range of covalent enzyme inhibitors and to metalloenzymes, where metal coordination is essential for drug binding. This updated version of our docking program was tested for its ability to predict the correct binding mode of drug-sized molecules in a large variety of proteins. We also report new datasets that were essential to demonstrate areas of success and those where additional efforts are required. This resource could be used by other program developers to assess their own software.

Augmented base pairing networks encode RNA-small molecule binding preferences.

RNA-small molecule binding is a key regulatory mechanism which can stabilize 3D structures and activate molecular functions. The discovery of RNA-targeting compounds is thus a current topic of interest for novel therapies. Our work is a first attempt at bringing the scalability and generalization abilities of machine learning methods to the problem of RNA drug discovery, as well as a step towards understanding the interactions which drive binding specificity. Our tool, RNAmigos, builds and encodes a network representation of RNA structures to predict likely ligands for novel binding sites. We subject ligand predictions to virtual screening and show that we are able to place the true ligand in the 71st-73rd percentile in two decoy libraries, showing a significant improvement over several baselines, and a state of the art method. Furthermore, we observe that augmenting structural networks with non-canonical base pairing data is the only representation able to uncover a significant signal, suggesting that such interactions are a necessary source of binding specificity. We also find that pre-training with an auxiliary graph representation learning task significantly boosts performance of ligand prediction. This finding can serve as a general principle for RNA structure-function prediction when data is scarce. RNAmigos shows that RNA binding data contains structural patterns with potential for drug discovery, and provides methodological insights for possible applications to other structure-function learning tasks. The source code, data and a Web server are freely available at http://rnamigos.cs.mcgill.ca.

Structural Polymorphism of Guanine Quadruplex-Containing Regions in Human Promoters.

Intramolecular guanine quadruplexes (G4s) are non-canonical nucleic acid structures formed by four guanine (G)-rich tracts that assemble into a core of stacked planar tetrads. G4-forming DNA sequences are enriched in gene promoters and are implicated in the control of gene expression. Most G4-forming DNA contains more G residues than can simultaneously be incorporated into the core resulting in a variety of different possible G4 structures. Although this kind of structural polymorphism is well recognized in the literature, there remain unanswered questions regarding possible connections between G4 polymorphism and biological function. Here we report a detailed bioinformatic survey of G4 polymorphism in human gene promoter regions. Our analysis is based on identifying G4-containing regions (G4CRs), which we define as stretches of DNA in which every residue can form part of a G4. We found that G4CRs with higher degrees of polymorphism are more tightly clustered near transcription sites and tend to contain G4s with shorter loops and bulges. Furthermore, we found that G4CRs with well-characterized biological functions tended to be longer and more polymorphic than genome-wide averages. These results represent new evidence linking G4 polymorphism to biological function and provide new criteria for identifying biologically relevant G4-forming regions from genomic data.

Mutations in Dynamic Structural Elements Alter the Kinetics and Fidelity of the Multifunctional Class II Lanthipeptide Synthetase, HalM2.

Class II lanthipeptide synthetases (LanM enzymes) catalyze the multistep post-translational modification of genetically encoded precursor peptides into macrocyclic (often antimicrobial) lanthipeptides. The reaction sequence involves dehydration of serine/threonine residues, followed by intramolecular addition of cysteine thiols onto the nascent dehydration sites to construct thioether bridges. LanMs utilize two separate active sites in an iterative yet highly coordinated manner to maintain a remarkable level of regio- and stereochemical control over the multistep maturation. The mechanisms underlying this biosynthetic fidelity remain enigmatic. We recently demonstrated that proper function of the haloduracin ß synthetase (HalM2) requires dynamic structural elements scattered across the surface of the enzyme. Here, we perform kinetic simulations, structural analysis of reaction intermediates, hydrogen-deuterium exchange mass spectrometry studies, and molecular dynamics simulations to investigate the contributions of these dynamic HalM2 structural elements to biosynthetic efficiency and fidelity. Our studies demonstrate that a large, conserved loop (HalM2 residues P349-P405) plays essential roles in defining the precursor peptide binding site, facilitating efficient peptide dehydration, and guiding the order of thioether ring formation. Moreover, mutations near the interface of the HalM2 dehydratase and cyclase domains perturb cyclization fidelity and result in aberrant thioether topologies that cannot be corrected by the wild type enzyme, suggesting an element of kinetic control in the normal cyclization sequence. Overall, this work provides the most comprehensive correlation of the structural and functional properties of a LanM enzyme reported to date and should inform mechanistic studies of the biosynthesis of other ribosomally synthesized and post-translationally modified peptide natural products.

Automated, customizable and efficient identification of 3D base pair modules with BayesPairing.

RNA structures possess multiple levels of structural organization. A secondary structure, made of Watson-Crick helices connected by loops, forms a scaffold for the tertiary structure. The 3D structures adopted by these loops are therefore critical determinants shaping the global 3D architecture. Earlier studies showed that these local 3D structures can be described as conserved sets of ordered non-Watson-Crick base pairs called RNA structural modules. Unfortunately, the computational efficiency and scope of the current 3D module identification methods are too limited yet to benefit from all the knowledge accumulated in the module databases. We present BayesPairing, an automated, efficient and customizable tool for (i) building Bayesian networks representing RNA 3D modules and (ii) rapid identification of 3D modules in sequences. BayesPairing uses a flexible definition of RNA 3D modules that allows us to consider complex architectures such as multi-branched loops and features multiple algorithmic improvements. We benchmarked our methods using cross-validation techniques on 3409 RNA chains and show that BayesPairing achieves up to ∼70% identification accuracy on module positions and base pair interactions. BayesPairing can handle a broader range of motifs (versatility) and offers considerable running time improvements (efficiency), opening the door to a broad range of large-scale applications.

Use of Extended-Hückel Descriptors for Rapid and Accurate Predictions of Conjugated Torsional Energy Barriers.

Over the past few decades, virtual high-throughput screening (vHTS) and molecular dynamics simulations have become effective and widely used tools in the initial stages of drug discovery efforts. These methods allow a great number of druglike molecules to be screened quickly and inexpensively. Unfortunately, however, the accuracies of both these methods rely on the quality of the underlying molecular mechanics force fields (FFs), which are often poor. This major weakness originates from the reliance of FFs on a finite list of specific parameters, called atom types, which have low transferability between molecules. In particular, the torsional energy barriers of druglike molecules are notoriously difficult to predict. Continuing our endeavor to understand factors affecting the torsional energy barriers of small molecules and quantify them, we showed that descriptors calculated using the extended-Hückel method could be used to rapidly assign accurate torsion parameters for conjugated molecules. This method, called H-TEQ 4.5, was developed using a set of 684 conjugated molecules. It was subsequently validated on a test set of 200 diverse molecules and produced an average root-mean-square error (rmse) of 1.01 kcal·mol-1, with respect to the reference quantum mechanic torsional profiles. For comparison, GAFF2, MMFF94, and MAB produced average rmse's of 3.49, 1.50, and 1.77 kcal·mol-1, respectively. H-TEQ 4.5 is also computationally inexpensive, running just under 0.25 ms for a biphenyl molecule on a home computer, allowing it to be used for vHTS of large libraries of compounds. Overall, H-TEQ 4.5 solved the problems associated with the transferability of torsion parameters for conjugated molecules. This method was incorporated into the Molecular Operating Environment and will be available for a wide variety of applications.

Predicting Positions of Bridging Water Molecules in Nucleic Acid-Ligand Complexes.

Over the past two decades, interests in DNA and RNA as drug targets have been growing rapidly. Following the trends observed with protein drug targets, computational approaches for drug design have been developed for this new class of molecules. Our efforts toward the development of a universal docking program, Fitted, led us to focus on nucleic acids. Throughout the development of this docking program, efforts were directed toward displaceable water molecules which must be accurately located for optimal docking-based drug discovery. However, although there is a plethora of methods to place water molecules in and around protein structures, there is, to the best of our knowledge, no such fully automated method for nucleic acids, which are significantly more polar and solvated than proteins. We report herein a new method, Splash'Em (Solvation Potential Laid around Statistical Hydration on Entire Macromolecules) developed to place water molecules within the binding cavity of nucleic acids. This fast method was shown to have high agreement with water positions in crystal structures and will therefore provide essential information to medicinal chemists.

Atom Type Independent Modeling of the Conformational Energy of Benzylic, Allylic, and Other Bonds Adjacent to Conjugated Systems.

Applications of computational methods to predict binding affinities for protein/drug complexes are routinely used in structure-based drug discovery. Applications of these methods often rely on empirical force fields (FFs) and their associated parameter sets and atom types. However, it is widely accepted that FFs cannot accurately cover the entire chemical space of drug-like molecules, due to the restrictive cost of parametrization and the poor transferability of existing parameters. To address these limitations, initiatives have been carried out to develop more transferable methods, in order to allow for more rigorous descriptions of any drug-like molecule. We have previously reported H-TEQ, a method which does not rely on atom types and incorporates well established chemical principles to assign parameters to organic molecules. The previous implementation of H-TEQ (a torsional barrier prediction method) only covered saturated and lone pair containing molecules; here, we report our efforts to incorporate conjugated systems into our model. The next step was the evaluation of the introduction of unsaturations. The developed model (H-TEQ3.0) has been validated on a wide variety of molecules containing heteroaromatic groups, alkyls, and fused ring systems. Our method performs on par with one of the most commonly used FFs (GAFF2), without relying on atom types or any prior parametrization.

Torsional Energy Barriers of Biaryls Could Be Predicted by Electron Richness/Deficiency of Aromatic Rings; Advancement of Molecular Mechanics toward Atom-Type Independence.

Biaryl molecules are ubiquitous pharmacophores found in natural products and pharmaceuticals. In spite of this, existing molecular mechanics force fields are unable to accurately reproduce their torsional energy profiles, except for a few well-parametrized cases. This effectively limits the ability of structure-based drug design methods to correctly identify hits involving biaryls with confidence (e.g., during virtual screening, employing docking and/or molecular dynamics simulations). Continuing in our endeavor to quantify organic chemistry principles, we showed that the torsional energy profile of biaryl compounds could be computed on-the-fly based on the electron richness/deficiency of the aromatic rings. This method, called H-TEQ 4.0, was developed using a set of 131 biaryls. It was subsequently validated on a separate set of 100 diverse biaryls, including multisubstituted, bicyclic and tricyclic druglike molecules, and produced an average root-mean-square error (RMSE) of 0.95 kcal·mol-1. For comparison, GAFF2 produced an RMSE of 3.88 kcal·mol-1, owing to problems associated with the transferability of torsion parameters. The success of H-TEQ 4.0 provided further evidence that force fields could transition to become atom-type independent, providing that the correct chemical principles are used. Overall, this method solved the problem of transferability of biaryl torsion parameters, while simultaneously improving the overall accuracy of the force field.

Complete Kinetic Characterization of Enzyme Inhibition in a Single Isothermal Titration Calorimetric Experiment.

Techniques for rapidly measuring both the strength and mode of enzyme inhibitors are crucial to lead generation and optimization in drug development. Isothermal titration calorimetry (ITC) is emerging as a powerful tool for measuring enzyme kinetics with distinct advantages over traditional techniques. ITC measures heat flow, a feature of nearly all chemical reactions, and gives an instantaneous readout of enzyme velocity, eliminating the need for artificial substrates or postreaction processing. In principle, ITC is an ideal method for characterizing enzyme inhibition. However, existing ITC experiments are not well-suited to rapid throughput and few studies to date have employed this approach. We have developed a new ITC experiment, in which substrate and inhibitor are premixed in the injection syringe, that yields complete kinetic characterization of an enzyme inhibitor in an hour or less. This corresponds to savings in time and material of 5-fold or greater compared to previous ITC methods. We validated the approach using the trypsin inhibitor benzamidine as a model system, recapitulating both its competitive inhibition mode and binding constant. Our approach combines the rapid throughput of optimized spectroscopic assays with the universality and precision of ITC-based methods, providing substantially improved inhibitor characterization for biochemistry and drug development applications.

Exploring Atypical Fluorine-Hydrogen Bonds and Their Effects on Nucleoside Conformations.

The ability of fluorine to serve as a hydrogen-bond acceptor has been debated for many years. Short fluorine-hydrogen contacts are thought to play a key role in stabilizing some complex supramolecular systems. To directly probe the existence of fluorine-hydrogen bonds, we have performed NMR spectroscopy and computational modeling on a series of C2'-fluorinated nucleosides. Specifically, quantum mechanics/molecular mechanics (QM/MM) analysis and [19 F,1 H] HMBC NMR experiments provided direct evidence for a C-H⋅⋅⋅F hydrogen bond in a 2'-F,4'-C-α-alkyl-ribonucleoside analogue. This interaction was also supported by QTAIM and NBO analyses, which confirmed a bond critical point for the C-H⋅⋅⋅F interaction (0.74 kcal mol-1 ). In contrast, although conformational analysis and NMR experiments of 2'-deoxy-2'-fluoro-arabinonucleosides indicated a close proximity between the 2'-fluorine and the H6/8 protons of the nucleobase, molecular simulations did not provide evidence for a C-H⋅⋅⋅F hydrogen bond.

Adjusting the Structure of 2'-Modified Nucleosides and Oligonucleotides via C4'-α-F or C4'-α-OMe Substitution: Synthesis and Conformational Analysis.

We report the first syntheses of three nucleoside analogues, namely, 2',4'-diOMe-rU, 2'-OMe,4'-F-rU, and 2'-F,4'-OMe-araU, via stereoselective introduction of fluorine or methoxy functionalities at the C4'-α-position of a 4',5'-olefinic intermediate. Conformational analyses of these nucleosides and comparison to other previously reported 2',4'-disubstituted nucleoside analogues make it possible to evaluate the effect of fluorine and methoxy substitution on the sugar pucker, as assessed by NMR, X-ray diffraction, and computational methods. We found that C4'-α-F/OMe substituents reinforce the C3'-endo ( north) conformation of 2'-OMe-rU. Furthermore, the predominant C2'-endo ( south/ east) conformation of 2'-F-araU switches to C3'-endo upon introduction of these substituents at C4'. The nucleoside analogues were incorporated into DNA and RNA oligonucleotides via standard phosphoramidite chemistry, and their effects on the thermal stability of homo- and heteroduplexes were assessed via UV thermal melting experiments. We found that 4'-substituents can modulate the binding affinity of the parent 2'-modified oligomers, inducing a mildly destabilizing or stabilizing effect depending on the duplex type. This study expands the spectrum of oligonucleotide modifications available for rational design of oligonucleotide therapeutics.

Atom Types Independent Molecular Mechanics Method for Predicting the Conformational Energy of Small Molecules.

We previously implemented a well-known qualitative chemical principle into an accurate quantitative model computing relative potential energies of conformers. According to this principle, hyperconjugation strength correlates with electronegativity of donors and acceptors. While this earlier version of our model applies to σ bonds, lone pairs, disregarded in this earlier version, also have a major impact on the conformational preferences of molecules. Among the well-established principles used by organic chemists to rationalize some organic chemical behaviors are the anomeric effect, the alpha effect, basicity, and nucleophilicity. These effects are directly related to the presence of lone pairs. We report herein our effort to incorporate lone pairs into our model to extend its applicability domain to any saturated small molecules. The developed model H-TEQ 2 has been validated on a wide variety of molecules from polyaromatic molecules to carbohydrates and molecules with high heteroatoms/carbon ratios. Interestingly, this method, in contrast to common force field-based methods, does not rely on atom types and is virtually applicable to any organic molecules.

Accurately Modeling the Conformational Preferences of Nucleosides.

Sugar puckering of nucleosides impacts nucleic acid structures; hence their biological function. Similarly, nucleoside-based therapeutics may adopt different conformations affecting their binding affinity, DNA incorporation, and excision rates. As a result, significant efforts have been made to develop nucleoside analogues adopting specific conformations to improve bioactivity and pharmacokinetic profiles of the corresponding nucleoside-containing drugs. Understanding and ultimately predicting these conformational preferences would significantly help in the design of more effective structures. We report herein a computational study based on hybrid QM/MM umbrella sampling simulations that allow the accurate prediction of the sugar conformational preferences of chemically modified nucleosides in solution. Moreover, we pair these simulations with natural bond orbital (NBO) analysis to gain key insights into the role of substituents in the conformational preferences of these nucleosides.

4'-C-Methoxy-2'-deoxy-2'-fluoro Modified Ribonucleotides Improve Metabolic Stability and Elicit Efficient RNAi-Mediated Gene Silencing.

We designed novel 4'-modified 2'-deoxy-2'-fluorouridine (2'-F U) analogues with the aim to improve nuclease resistance and potency of therapeutic siRNAs by introducing 4'-C-methoxy (4'-OMe) as the alpha (C4'α) or beta (C4'ß) epimers. The C4'α epimer was synthesized by a stereoselective route in six steps; however, both α and ß epimers could be obtained by a nonstereoselective approach starting from 2'-F U. 1H NMR analysis and computational investigation of the α-epimer revealed that the 4'-OMe imparts a conformational bias toward the North-East sugar pucker, due to intramolecular hydrogen bonding and hyperconjugation effects. The α-epimer generally conceded similar thermal stability as unmodified nucleotides, whereas the ß-epimer led to significant destabilization. Both 4'-OMe epimers conferred increased nuclease resistance, which can be explained by the close proximity between 4'-OMe substituent and the vicinal 5'- and 3'-phosphate group, as seen in the X-ray crystal structure of modified RNA. siRNAs containing several C4'α-epimer monomers in the sense or antisense strands triggered RNAi-mediated gene silencing with efficiencies comparable to that of 2'-F U.

Measuring Rapid Time-Scale Reaction Kinetics Using Isothermal Titration Calorimetry.

Isothermal titration calorimetry (ITC) is a powerful tool for acquiring both thermodynamic and kinetic data for biological interactions including molecular recognition and enzymatic catalysis. ITC-based kinetics measurements typically focus on reactions taking place over long time scales (tens of minutes or hours) in order to avoid complications due to the finite length of time needed detect heat flow in the calorimeter cell. While progress has been made toward analyzing more rapid reaction kinetics by ITC, the capabilities and limitations of this approach have not been thoroughly tested to date. Here, we report that the time resolution of commercial instruments is on the order of 0.2 s or less. We successfully performed rapid ITC kinetics assays with durations of just tens of seconds using the enzyme trypsin. This is substantially shorter than previous ITC enzyme measurements. However, we noticed that for short reaction durations, standard assumptions regarding the ITC instrument response led to significant deviations between calculated and measured ITC peak shapes. To address this issue, we developed an ITC empirical response model (ITC-ERM) that quantitatively reproduces ITC peak shapes for all reaction durations. Applying the ITC-ERM approach to another enzyme (prolyl oligopeptidase), we unexpectedly discovered non-Michaelis-Menten kinetics in short time-scale measurements that are absent in more typical long time-scale experiments and are obscured in short time-scale experiments when standard assumptions regarding the instrument response are made. This highlights the potential of ITC measurements of rapid time scale kinetics in conjunction with the ITC-ERM approach to shed new light on biological dynamics.

Active Site Crowding of Cytochrome P450 3A4 as a Strategy To Alter Its Selectivity.

Substrate-promiscuous enzymes are a promising starting point for the development of versatile biocatalysts. In this study, human cytochrome P450 3A4, known for its ability to metabolise hundreds of drugs, was engineered to alter its regio- and stereoselectivity. Rational mutagenesis was used to introduce steric hindrance in a specific manner in the large active site of P450 3A4 and to favour oxidation at a more sterically accessible position on the substrate. Hydroxylation of a synthetic precursor of (R)-lisofylline, a compound under investigation for its anti-inflammatory properties, was chosen as a first proof-of-principle application of our protein engineering strategy. In a second example, increasing active site crowding led to an incremental shift in the selectivity of oxidation from an internal double bond to a terminal phenyl group in a derivative of theobromine. The same correlation between crowding and selectivity was found in a final case focused on the hydroxylation of the steroid sex hormone progesterone.

Medicinal Chemistry Projects Requiring Imaginative Structure-Based Drug Design Methods.

Computational methods for docking small molecules to proteins are prominent in drug discovery. There are hundreds, if not thousands, of documented examples-and several pertinent cases within our research program. Fifteen years ago, our first docking-guided drug design project yielded nanomolar metalloproteinase inhibitors and illustrated the potential of structure-based drug design. Subsequent applications of docking programs to the design of integrin antagonists, BACE-1 inhibitors, and aminoglycosides binding to bacterial RNA demonstrated that available docking programs needed significant improvement. At that time, docking programs primarily considered flexible ligands and rigid proteins. We demonstrated that accounting for protein flexibility, employing displaceable water molecules, and using ligand-based pharmacophores improved the docking accuracy of existing methods-enabling the design of bioactive molecules. The success prompted the development of our own program, Fitted, implementing all of these aspects. The primary motivation has always been to respond to the needs of drug design studies; the majority of the concepts behind the evolution of Fitted are rooted in medicinal chemistry projects and collaborations. Several examples follow: (1) Searching for HDAC inhibitors led us to develop methods considering drug-zinc coordination and its effect on the pKa of surrounding residues. (2) Targeting covalent prolyl oligopeptidase (POP) inhibitors prompted an update to Fitted to identify reactive groups and form bonds with a given residue (e.g., a catalytic residue) when the geometry allows it. Fitted-the first fully automated covalent docking program-was successfully applied to the discovery of four new classes of covalent POP inhibitors. As a result, efficient stereoselective syntheses of a few screening hits were prioritized rather than synthesizing large chemical libraries-yielding nanomolar inhibitors. (3) In order to study the metabolism of POP inhibitors by cytochrome P450 enzymes (CYPs)-for toxicology studies-the program Impacts was derived from Fitted and helped us to reveal a complex metabolism with unforeseen stereocenter isomerizations. These efforts, combined with those of other docking software developers, have strengthened our understanding of the complex drug-protein binding process while providing the medicinal chemistry community with useful tools that have led to drug discoveries. In this Account, we describe our contributions over the past 15 years-within their historical context-to the design of drug candidates, including BACE-1 inhibitors, POP covalent inhibitors, G-quadruplex binders, and aminoglycosides binding to nucleic acids. We also remark the necessary developments of docking programs, specifically Fitted, that enabled structure-based design to flourish and yielded multiple fruitful, rational medicinal chemistry campaigns.

Fluoride-Mediated Desulfonylative Intramolecular Cyclization to Fused and Bridged Bicyclic Compounds: A Complex Mechanism.

We previously reported the synthesis of polysubstituted chiral oxazepanes in three steps from commercially available starting materials. The unexpected reaction of one of these 1,4-oxazepanes in the presence of TBAF provided a 4-oxa-1-azabicyclo[4.1.0]heptane core. This unusual process significantly increased the complexity of the molecular scaffold by introducing a bicyclic core. Surprisingly, the generated bicyclic structure featuring three stereocenters was a mixture of enantiomers with no other diastereomers observed. These striking experimental observations deserved further investigations. A combination of experimental and computational investigations unveiled a complex diastereoselective mechanism. Mechanistic rationale is presented for this observed rearrangement.