Mechanism-Based Redesign of GAP to Activate Oncogenic Ras.

Ras GTPases play a crucial role in cell signaling pathways. Mutations of the Ras gene occur in about one third of cancerous cell lines and are often associated with detrimental clinical prognosis. Hot spot residues Gly12, Gly13, and Gln61 cover 97% of oncogenic mutations, which impair the enzymatic activity in Ras. Using QM/MM free energy calculations, we present a two-step mechanism for the GTP hydrolysis catalyzed by the wild-type Ras.GAP complex. We found that the deprotonation of the catalytic water takes place via the Gln61 as a transient Brønsted base. We also determined the reaction profiles for key oncogenic Ras mutants G12D and G12C using QM/MM minimizations, matching the experimentally observed loss of catalytic activity, thereby validating our reaction mechanism. Using the optimized reaction paths, we devised a fast and accurate procedure to design GAP mutants that activate G12D Ras. We replaced GAP residues near the active site and determined the activation barrier for 190 single mutants. We furthermore built a machine learning for ultrafast screening, by fast prediction of the barrier heights, tested both on the single and double mutations. This work demonstrates that fast and accurate screening can be accomplished via QM/MM reaction path optimizations to design protein sequences with increased catalytic activity. Several GAP mutations are predicted to re-enable catalysis in oncogenic G12D, offering a promising avenue to overcome aberrant Ras-driven signal transduction by activating enzymatic activity instead of inhibition. The outlined computational screening protocol is readily applicable for designing ligands and cofactors analogously.

Dynamic regulation of RAS and RAS signaling.

RAS proteins regulate most aspects of cellular physiology. They are mutated in 30% of human cancers and 4% of developmental disorders termed Rasopathies. They cycle between active GTP-bound and inactive GDP-bound states. When active, they can interact with a wide range of effectors that control fundamental biochemical and biological processes. Emerging evidence suggests that RAS proteins are not simple on/off switches but sophisticated information processing devices that compute cell fate decisions by integrating external and internal cues. A critical component of this compute function is the dynamic regulation of RAS activation and downstream signaling that allows RAS to produce a rich and nuanced spectrum of biological outputs. We discuss recent findings how the dynamics of RAS and its downstream signaling is regulated. Starting from the structural and biochemical properties of wild-type and mutant RAS proteins and their activation cycle, we examine higher molecular assemblies, effector interactions and downstream signaling outputs, all under the aspect of dynamic regulation. We also consider how computational and mathematical modeling approaches contribute to analyze and understand the pleiotropic functions of RAS in health and disease.

Combined Free-Energy Calculation and Machine Learning Methods for Understanding Ligand Unbinding Kinetics.

The determination of drug residence times, which define the time an inhibitor is in complex with its target, is a fundamental part of the drug discovery process. Synthesis and experimental measurements of kinetic rate constants are, however, expensive and time consuming. In this work, we aimed to obtain drug residence times computationally. Furthermore, we propose a novel algorithm to identify molecular design objectives based on ligand unbinding kinetics. We designed an enhanced sampling technique to accurately predict the free-energy profiles of the ligand unbinding process, focusing on the free-energy barrier for unbinding. Our method first identifies unbinding paths determining a corresponding set of internal coordinates (ICs) that form contacts between the protein and the ligand; it then iteratively updates these interactions during a series of biased molecular dynamics (MD) simulations to reveal the ICs that are important for the whole of the unbinding process. Subsequently, we performed finite-temperature string simulations to obtain the free-energy barrier for unbinding using the set of ICs as a complex reaction coordinate. Importantly, we also aimed to enable the further design of drugs focusing on improved residence times. To this end, we developed a supervised machine learning (ML) approach with inputs from unbiased "downhill" trajectories initiated near the transition state (TS) ensemble of the string unbinding path. We demonstrate that our ML method can identify key ligand-protein interactions driving the system through the TS. Some of the most important drugs for cancer treatment are kinase inhibitors. One of these kinase targets is cyclin-dependent kinase 2 (CDK2), an appealing target for anticancer drug development. Here, we tested our method using two different CDK2 inhibitors for the potential further development of these compounds. We compared the free-energy barriers obtained from our calculations with those observed in available experimental data. We highlighted important interactions at the distal ends of the ligands that can be targeted for improved residence times. Our method provides a new tool to determine unbinding rates and to identify key structural features of the inhibitors that can be used as starting points for novel design strategies in drug discovery.

Modelling the active SARS-CoV-2 helicase complex as a basis for structure-based inhibitor design.

The RNA helicase (non-structural protein 13, NSP13) of SARS-CoV-2 is essential for viral replication, and it is highly conserved among the coronaviridae family, thus a prominent drug target to treat COVID-19. We present here structural models and dynamics of the helicase in complex with its native substrates based on thorough analysis of homologous sequences and existing experimental structures. We performed and analysed microseconds of molecular dynamics (MD) simulations, and our model provides valuable insights to the binding of the ATP and ssRNA at the atomic level. We identify the principal motions characterising the enzyme and highlight the effect of the natural substrates on this dynamics. Furthermore, allosteric binding sites are suggested by our pocket analysis. Our obtained structural and dynamical insights are important for subsequent studies of the catalytic function and for the development of specific inhibitors at our characterised binding pockets for this promising COVID-19 drug target.

Resolving sub-angstrom ambient motion through reconstruction from vibrational spectra.

Metal/organic-molecule interactions underpin many key chemistries but occur on sub-nm scales where nanoscale visualisation techniques tend to average over heterogeneous distributions. Single molecule imaging techniques at the atomic scale have found it challenging to track chemical behaviour under ambient conditions. Surface-enhanced Raman spectroscopy can optically monitor the vibrations of single molecules but understanding is limited by the complexity of spectra and mismatch between theory and experiment. We demonstrate that spectra from an optically generated metallic adatom near a molecule of interest can be inverted into dynamic sub-Å metal-molecule interactions using a comprehensive model, revealing anomalous diffusion of a single atom. Transient metal-organic coordination bonds chemically perturb molecular functional groups > 10 bonds away. With continuous improvements in computational methods for modelling large and complex molecular systems, this technique will become increasingly applicable to accurately tracking more complex chemistries.

Toward Understanding CB[7]-Based Supramolecular Diels-Alder Catalysis.

Cucurbiturils (CBs) are robust and versatile macrocyclic compounds, often used as molecular hosts in complex supramolecular systems. In previous work, remarkable catalytic activity has been observed for asymmetric cycloadditions under very mild conditions. Herein, we investigate the nature of supramolecular catalysis using DFT calculations and QM/MM techniques. We discuss induced conformational changes, electrostatic shielding effects from the highly polar aqueous environment and cooperativity in hydrogen bonding of the substrates in explicit water using QM/MM simulation techniques. Our results show little specificity for the chosen molecules, suggesting an excellent opportunity to expand the scope for catalytic use of these supramolecular macrocyclic containers.

Optical probes of molecules as nano-mechanical switches.

Molecular electronics promises a new generation of ultralow-energy information technologies, based around functional molecular junctions. Here, we report optical probing that exploits a gold nanoparticle in a plasmonic nanocavity geometry used as one terminal of a well-defined molecular junction, deposited as a self-assembled molecular monolayer on flat gold. A conductive transparent cantilever electrically contacts individual nanoparticles while maintaining optical access to the molecular junction. Optical readout of molecular structure in the junction reveals ultralow-energy switching of ∼50 zJ, from a nano-electromechanical torsion spring at the single molecule level. Real-time Raman measurements show these electronic device characteristics are directly affected by this molecular torsion, which can be explained using a simple circuit model based on junction capacitances, confirmed by density functional theory calculations. This nanomechanical degree of freedom is normally invisible and ignored in electrical transport measurements but is vital to the design and exploitation of molecules as quantum-coherent electronic nanodevices.

Metal-Free C-H Borylation of N-Heteroarenes by Boron Trifluoride.

Organoboron compounds are essential reagents in modern C-C coupling reactions. Their synthesis via catalytic C-H borylation by main group elements is emerging as a powerful tool alternative to transition metal based catalysis. Herein, a straightforward metal-free synthesis of aryldifluoroboranes from BF3 and heteroarenes is reported. The reaction is assisted by sterically hindered amines and catalytic amounts of thioureas. According to computational studies the reaction proceeds via frustrated Lewis pair (FLP) mechanism. The obtained aryldifluoroboranes are further stabilized against destructive protodeborylation by converting them to the corresponding air stable tetramethylammonium organotrifluoroborates.

Nuclear Quantum Effects from the Analysis of Smoothed Trajectories: Pilot Study for Water.

Nuclear quantum effects have significant contributions to thermodynamic quantities and structural properties; furthermore, very expensive methods are necessary for their accurate computation. In most calculations, these effects, for instance, zero-point energies, are simply neglected or only taken into account within the quantum harmonic oscillator approximation. Herein, we present a new method, Generalized Smoothed Trajectory Analysis, to determine nuclear quantum effects from molecular dynamics simulations. The broad applicability is demonstrated with the examples of a harmonic oscillator and different states of water. Ab initio molecular dynamics simulations have been performed for ideal gas up to the temperature of 5000 K. Classical molecular dynamics have been carried out for hexagonal ice, liquid water, and vapor at atmospheric pressure. With respect to the experimental heat capacity, our method outperforms previous calculations in the literature in a wide temperature range at lower computational cost than other alternatives. Dynamic and structural nuclear quantum effects of water are also discussed.

Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H+ and Mg2+ ions in enzyme reactions.

Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg2+ ions, and recent advances in their transient involvement in the catalytic mechanisms.

Structural basis of second-generation HIV integrase inhibitor action and viral resistance.

Although second-generation HIV integrase strand-transfer inhibitors (INSTIs) are prescribed throughout the world, the mechanistic basis for the superiority of these drugs is poorly understood. We used single-particle cryo-electron microscopy to visualize the mode of action of the advanced INSTIs dolutegravir and bictegravir at near-atomic resolution. Glutamine-148→histidine (Q148H) and glycine-140→serine (G140S) amino acid substitutions in integrase that result in clinical INSTI failure perturb optimal magnesium ion coordination in the enzyme active site. The expanded chemical scaffolds of second-generation compounds mediate interactions with the protein backbone that are critical for antagonizing viruses containing the Q148H and G140S mutations. Our results reveal that binding to magnesium ions underpins a fundamental weakness of the INSTI pharmacophore that is exploited by the virus to engender resistance and provide a structural framework for the development of this class of anti-HIV/AIDS therapeutics.

Organocatalytic Access to a cis-Cyclopentyl-γ-amino Acid: An Intriguing Model of Selectivity and Formation of a Stable 10/12-Helix from the Corresponding γ/α-Peptide.

In this study, we have developed a highly enantioselective organocatalytic route to the (1S,2R)-2-(aminomethyl)cyclopentane-1-carboxylic acid monomer precursor, which has a cis-configuration between the C- and N-termini around the cyclopentane core. Kinetic measurements show that the product distribution changes over time due to epimerization of the C1 center. Computations suggest the cis-selectivity is a result of selective C-C bond formation, while subsequent steps appear to influence the selectivity at higher temperature. The resulting γ-amino acid residue was incorporated into a novel γ/α-peptide, which forms a well-ordered 10/12-helix with alternate H-bond directionality in spite of the smallest value of the ζ-angle yet observed for a helix of this type. This highly defined structure is also a result of the narrow range of potential ζ-angles in our monomer. In contrast, the larger range of potential ζ-values observed for the corresponding trans-system can be fulfilled by several competing helical structures.

Development of a True Transition State Force Field from Quantum Mechanical Calculations.

Transition state force fields (TSFF) treated the TS structure as an artificial minimum on the potential energy surface in the past decades. The necessary parameters were developed either manually or by the Quantum-to-molecular mechanics method (Q2MM). In contrast with these approaches, here we propose to model the TS structures as genuine saddle points at the molecular mechanics level. Different methods were tested on small model systems of general chemical reactions such as protonation, nucleophilic attack, and substitution, and the new procedure led to more accurate models than the Q2MM-type parametrization. To demonstrate the practicality of our approach, transferrable parameters have been developed for Mo-catalyzed olefin metathesis using quantum mechanical properties as reference data. Based on the proposed strategy, any force field can be extended with true transition state force field (TTSFF) parameters, and they can be readily applied in several molecular mechanics programs as well.

A Catalyst Designed for the Enantioselective Construction of Methyl- and Alkyl-Substituted Tertiary Stereocenters.

Tertiary methyl-substituted stereocenters are present in numerous biologically active natural products. Reported herein is a catalytic enantioselective method for accessing these chiral building blocks using the Mukaiyama-Michael reaction between silyl ketene thioacetals and acrolein. To enable remote enantioface control on the nucleophile, a new iminium catalyst, optimized by three-parameter tuning and by identifying substituent effects on enantioselectivity, was designed. The catalytic process allows rapid access to chiral thioesters, amides, aldehydes, and ketones bearing an α-methyl stereocenter with excellent enantioselectivities, and allowed rapid access to the C4-C13 segment of (-)-bistramide A. DFT calculations rationalized the observed sense and level of enantioselectivity.