Mechanistic Insights into S-Depalmitolyse Activity of Cln5 Protein Linked to Neurodegeneration and Batten Disease: A QM/MM Study.

Ceroid lipofuscinosis neuronal protein 5 (Cln5) is encoded by the CLN5 gene. The genetic variants of this gene are associated with the CLN5 form of Batten disease. Recently, the first crystal structure of Cln5 was reported. Cln5 shows cysteine palmitoyl thioesterase S-depalmitoylation activity, which was explored via fluorescent emission spectroscopy utilizing the fluorescent probe DDP-5. In this work, the mechanism of the reaction between Cln5 and DDP-5 was studied computationally by applying a QM/MM methodology at the ωB97X-D/6-31G(d,p):AMBER level. The results of our study clearly demonstrate the critical role of the catalytic triad Cys280-His166-Glu183 in S-depalmitoylation activity. This is evidenced through a comparison of the pathways catalyzed by the Cys280-His166-Glu183 triad and those with only Cys280 involved. The computed reaction barriers are in agreement with the catalytic efficiency. The calculated Gibb's free-energy profile suggests that S-depalmitoylation is a rate-limiting step compared to the preceding S-palmitoylation, with barriers of 26.1 and 25.3 kcal/mol, respectively. The energetics were complemented by monitoring the fluctuations in the electron density distribution through NBO charges and bond strength alterations via local mode stretching force constants during the catalytic pathways. This comprehensive protocol led to a more holistic picture of the reaction mechanism at the atomic level. It forms the foundation for future studies on the effects of gene mutations on both the S-palmitoylation and S-depalmitoylation steps, providing valuable data for the further development of enzyme replacement therapy, which is currently the only FDA-approved therapy for childhood neurodegenerative diseases, including Batten disease.

Unraveling the effect of aromaticity for the dynamics of excited states of single benzene fluorophores.

The photophysical properties of a series of recently synthesized single benzene fluorophores were investigated using ensemble density functional theory calculations. The energetic stability of the ground and excited state species were counterposed against the aromaticity index derived from local vibrational modes. It was found that the large Stokes shift of the fluorophores (up to ca. 5800 cm - 1 ) originates from the effect of electron donating and electron withdrawing substituents rather than π -delocalization and related (anti-)aromaticity. On the basis of nonadiabatic molecular dynamics simulations, the absence of fluorescence from one of the regioisomers was explained by the occurrence of easily accessible S 1 /S 0 conical intersections below the vertical excitation energy level. It is demonstrated in the manuscript that the analysis of local mode force constants and the related aromaticity index represent a useful tool for the characterization of π -delocalization effects in π -conjugated compounds.

Iron-histidine bonding in bishistidyl hemoproteins-A local vibrational mode study.

We investigated the intrinsic strength of distal and proximal FeN bonds for both ferric and ferrous oxidation states of bishistidyl hemoproteins from bacteria, animals, human, and plants, including two cytoglobins, ten hemoglobins, two myoglobins, six neuroglobins, and six phytoglobins. As a qualified measure of bond strength, we used local vibrational force constants k a (FeN) based on local mode theory developed in our group. All calculations were performed with a hybrid QM/MM ansatz. Starting geometries were taken from available x-ray structures. k a (FeN) values were correlated with FeN bond lengths and covalent bond character. We also investigated the stiffness of the axial NFeN bond angle. Our results highlight that protein effects are sensitively reflected in k a (FeN), allowing one to compare trends in diverse protein groups. Moreover, k a (NFeN) is a perfect tool to monitor changes in the axial heme framework caused by different protein environments as well as different Fe oxidation states.

Chemical bonding in Uranium-based materials: A local vibrational mode case study of Cs 2 UO 2 Cl 4 and UCl 4 crystals.

The Local Vibrational Mode Analysis, initially applied to diverse molecular systems, was extended to periodic systems in 2019. This work introduces an enhanced version of the LModeA software, specifically designed for the comprehensive analysis of two and three-dimensional periodic structures. Notably, a novel interface with the Crystal package was established, enabling a seamless transition from molecules to periodic systems using a unified methodology. Two distinct sets of uranium-based systems were investigated: (i) the evolution of the Uranyl ion (UO 2 2 + ) traced from its molecular configurations to the solid state, exemplified by Cs 2 UO 2 Cl 4 and (ii) Uranium tetrachloride (UCl 4 ) in both its molecular and crystalline forms. The primary focus was on exploring the impact of crystal packing on key properties, including IR and Raman spectra, structural parameters, and an in-depth assessment of bond strength utilizing local mode perspectives. This work not only demonstrates the adaptability and versatility of LModeA for periodic systems but also highlights its potential for gaining insights into complex materials and aiding in the design of new materials through fine-tuning.

Metal-ring interactions in group 2 ansa-metallocenes: assessed with the local vibrational mode theory.

Ansa-metallocenes, a vital class of organometallic compounds, have attracted significant attention due to their diverse structural motifs and their pivotal roles in catalysis and materials science. We investigated 37 distinct group 2 ansa-metallocenes at the B3LYP-D3/def2-TZVP level of theory. Utilizing local mode force constants derived from our local vibrational mode theory, including a special force constant directly targeting the metal-ring interaction, we could unveil latent structural differences between solvated and non-solvated metallocenophanes and the influence of the solvent on complex stability and structure. We could quantify the intrinsic strength of the metal-cyclopentadienyl (M-Cp) bonds and the influence of the bridging motifs on the stiffness of the Cp-M-Cp angles, another determinant of complex stability. LMA was complemented by the analysis of electronic density, utilizing the quantum theory of atoms in molecules (QTAIM), which confirmed both the impact of solvent coordination on the strength of the M-Cp bond(s) and the influence of the bridging motif on the Cp-M-Cp angles. The specific effect of the ansa-motif on the M-Cp interaction was further elucidated by a comparison with linear/bent metallocene structures. In summary, our results identify the local mode analysis as an efficient tool for unraveling the intricate molecular properties of ansa-metallocenes and their unique structural features.

Ligand Characterization and DNA Intercalation of Ru(II) Polypyridyl Complexes: A Local Vibrational Mode Study.

We investigated in this work ruthenium-ligand bonding across the RuN framework in 12 Ru(II) polypyridyl complexes in the gas phase and solution for both singlet and triplet states, in addition to their affinity for DNA binding through π-π stacking interactions with DNA nucleobases. As a tool to assess the intrinsic strength of the ruthenium-ligand bonds, we determined local vibrational force constants via our local vibrational mode analysis software. We introduced a novel local force constant that directly accounts for the intrinsic strength of the π-π stacking interaction between DNA and the intercalated Ru(II) complex. According to our findings, [Ru(phen)2(dppz)]2+ and [Ru(phen)2(11-CN-dppz)]2+ provide an intriguing trade-off between photoinduced complex excitation and the strength of the subsequent π-π stacking interaction with DNA. [Ru(phen)2(dppz)]2+ displays a small singlet-triplet splitting and a strong π-π stacking interaction in its singlet state, suggesting a favorable photoexcitation but potentially weaker interaction with DNA in the excited state. Conversely, [Ru(phen)2(11-CN-dppz)]2+ exhibits a larger singlet-triplet splitting and a stronger π-π stacking interaction with DNA in its triplet state, indicating a less favorable photoinduced transition but a stronger interaction with DNA postexcitation. We hope our study will inspire future experimental and computational work aimed at the design of novel Ru-polypyridyl drug candidates and that our new quantitative measure of π-π stacking interactions in DNA will find a general application in the field.

Noncovalent π Interactions in Mutated Aquomet-Myoglobin Proteins: A QM/MM and Local Vibrational Mode Study.

Protein dynamics and function is strongly connected to the energy flow taking place. Myoglobin (Mb) and its mutations are ideal systems to study the process of vibrational energy transfer (VET) at the molecular level. Anti-Stokes ultraviolet resonance Raman studies using a tryptophan (Trp) probe, introduced at different Mb positions by amino acid replacement, have suggested that the amount of VET depends on the position of the Trp probe relative to the heme group. Inspired by this experimental work, we explored the strength of noncovalent π interactions, as well as covalent interactions for both the axial and distal ligands bound to iron in aquomet-Mb with the local vibrational mode analysis (LMA), originally developed by Konkoli and Cremer. Two sets of noncovalent interactions were investigated: (1) the interaction between the water ligand and Trp rings and (2) the interaction between the Trp and the porphyrin rings of the heme group. We assessed the strength of these noncovalent interactions via a special local mode force constant. Various Trp-modified water-bound ferric Mb proteins in the ground state were studied (6 in total) using gas-phase and QM/MM calculations followed by LMA. Our results disclose that VET is indeed dependent on the position of the Trp probe relative to the heme group but also on the tautomeric nature of distal histidine. They provide new guidelines on how to assess noncovalent π interactions in proteins utilizing LMA and how to use these data to explore VET, and more generally protein dynamics and function.

Linear versus Bent Uranium(II) Metallocenes─A Local Vibrational Mode Study.

Uranium metallocenes have recently attracted attention driven by their use as catalysts in organometallic synthesis. In addition to bent U(IV) and U(III), an U(II) metallocene [(η5-C5i Pr5)2U] was synthesized with an unusual linear Cp-U-Cp angle. In this work, we investigated 22 U(II) metallocenes, (i) assessing the intrinsic strength of the U-ring interactions in these complexes with a novel bond strength measure based on our local vibrational mode analysis and (ii) systematically exploring what makes these U(II) metallocenes bent. We included relativistic effects through the NESCau Hamiltonian and complemented the local mode analysis with natural bonding orbital (NBO) and quantum theory of atoms in molecules (QTAIM) data. Our study led to the following results: (i) reduction of bulky U-ring ligand substituents does not lead to bent complexes for alkyl substituents (iPr and iBu) in contrast to SiMe3 ring substituents, which are all bent. (ii) The most bent complexes are [(η5-C5H4SiMe3)2U] (130°) and [η5-P5H5)2U] (143°). (iii) Linear complexes showed one hybridized NBO with s/d character, while bent structures were characterized by s/d/f mixing. (iv) We did not observe a correlation between the strength of the U-ring interaction and the amount of the ring-U-ring bend; the strongest interaction was found for [η5-Cp)2U] and the weakest for [η5-P5H5)2U]. In conclusion, our results provide a foundation for the design of U(II) metallocenes with specific physicochemical properties and increased reactivity.

Hydrogen Sulfide Ligation in Hemoglobin I of Lucina pectinata─A QM/MM and Local Mode Study.

In this study, we investigated the interaction between the H2S ligand and the heme pocket of hemoglobin I (HbI) of Lucina pectinata for the wild-type protein; three known mutations where distal glutamine is replaced by hydrophobic valine (Gln64Val) and hydrophilic histidine in both protonation forms (Gln64Hisϵ and Gln64Hisδ); five known mutations of the so-called phenyl cage, replacing the hydrophobic phenylalanines Phe29 and Phe43 with tyrosine (Tyr), valine (Val), or leucine (Leu); and two additional mutations, Phe68Tyr and Phe68Val, in order to complement previous studies with new insights about the binding mechanism at the molecular level. A particular focus was on the intrinsic strengths of the chemical bonds involved, utilizing local vibrational force constants based on combined quantum mechanical-molecular mechanical calculations. Wild-type protein and mutations clustered into two distinct groups: Group 1 protein systems with a proton acceptor in the distal protein pocket, close to one of the H2S bonds, and Group 2 protein systems without a hydrogen acceptor close by in the active site of the protein. According to our results, the interactions between H2S and HbI of Lucina pectinata involve two important elements, namely, binding of H2S to Fe of the heme group, followed by the proton transfer from the HS bond to the distal residue. The distal residue is additionally stabilized by a second proton transfer from the distal residue to COO- of the propionate group in heme. We could identify the FeS bond as a key player and discovered that the strength of this bond depends on two mutual factors, namely, the strength of the HS bond involved in the proton transfer and the electrostatic field of the protein pocket qualifying the FeS bond as a sensitive probe for monitoring changes in H2S ligation upon protein mutations. We hope our study will inspire and guide future experimental studies, targeting new promising mutations such as Phe68Tyr, Phe68Val, or Phe43Tyr/Phe68Val.

Quantum Mechanical Assessment of Protein-Ligand Hydrogen Bond Strength Patterns: Insights from Semiempirical Tight-Binding and Local Vibrational Mode Theory.

Hydrogen bonds (HB)s are the most abundant motifs in biological systems. They play a key role in determining protein-ligand binding affinity and selectivity. We designed two pharmaceutically beneficial HB databases, database A including ca. 12,000 protein-ligand complexes with ca. 22,000 HBs and their geometries, and database B including ca. 400 protein-ligand complexes with ca. 2200 HBs, their geometries, and bond strengths determined via our local vibrational mode analysis. We identified seven major HB patterns, which can be utilized as a de novo QSAR model to predict the binding affinity for a specific protein-ligand complex. Glycine was reported as the most abundant amino acid residue in both donor and acceptor profiles, and N-H⋯O was the most frequent HB type found in database A. HBs were preferred to be in the linear range, and linear HBs were identified as the strongest. HBs with HB angles in the range of 100-110°, typically forming intramolecular five-membered ring structures, showed good hydrophobic properties and membrane permeability. Utilizing database B, we found a generalized Badger's relationship for more than 2200 protein-ligand HBs. In addition, the strength and occurrence maps between each amino acid residue and ligand functional groups open an attractive possibility for a novel drug-design approach and for determining drug selectivity and affinity, and they can also serve as an important tool for the hit-to-lead process.

Unusual Intramolecular Motion of ReH92- in K2ReH9 Crystal: Circle Dance and Three-Arm Turnstile Mechanisms Revealed by Computational Studies.

The nonahydridorhenate dianion ReH92- is a unique rhenium polyhydride complex due to its remarkably high coordination number; however, its detailed polytopal rearrangement process in either solution or crystal is so far unclear. In this work, our quantum chemical calculations have identified two previously unreported fluxional mechanisms for the ReH92- dianion in the K2ReH9 crystal: three-arm turnstile rotation and circle dance mechanism. These two polytopal rearrangements in the crystal offer an alternative interpretation to the pulse and wide-line NMR spectra (Farrar et al. J. Chem. Phys. 1969, 51, 3595). The previously postulated hindered rotation of the whole ReH92- dianion in K2ReH9 (White et al. J. Chem. Soc., Faraday Trans. 2 1972, 68, 1414) turns out to be a combination of the above-mentioned two elementary fluxional processes. In addition, our calculations have confirmed the Muetterties' D3h⇌C4v rearrangement as the intramolecular motion for the ReH92- dianion in solution.

Dynamics of the Energy Transfer Process in Eu(III) Complexes Containing Polydentate Ligands Based on Pyridine, Quinoline, and Isoquinoline as Chromophoric Antennae.

In this work, we investigated from a theoretical point of view the dynamics of the energy transfer process from the ligand to Eu(III) ion for 12 isomeric species originating from six different complexes differing by nature of the ligand and the total charge. The cationic complexes present the general formula [Eu(L)(H2O)2]+ (where L = bpcd2- = N,N'-bis(2-pyridylmethyl)-trans-1,2-diaminocyclohexane N,N'-diacetate; bQcd2- = N,N'-bis(2-quinolinmethyl)-trans-1,2-diaminocyclohexane N,N'-diacetate; and bisoQcd2- = N,N'-bis(2-isoquinolinmethyl)-trans-1,2-diaminocyclohexane N,N'-diacetate), while the neutral complexes present the Eu(L)(H2O)2 formula (where L = PyC3A3- = N-picolyl-N,N',N'-trans-1,2-cyclohexylenediaminetriacetate; QC3A3- = N-quinolyl-N,N',N'-trans-1,2-cyclohexylenediaminetriacetate; and isoQC3A3- = N-isoquinolyl-N,N',N'-trans-1,2-cyclohexylenediaminetriacetate). Time-dependent density functional theory (TD-DFT) calculations provided the energy of the ligand excited donor states, distances between donor and acceptor orbitals involved in the energy transfer mechanism (RL), spin-orbit coupling matrix elements, and excited-state reorganization energies. The intramolecular energy transfer (IET) rates for both singlet-triplet intersystem crossing and ligand-to-metal (and vice versa) involving a multitude of ligand and Eu(III) levels and the theoretical overall quantum yields (ϕovl) were calculated (the latter for the first time without the introduction of experimental parameters). This was achieved using a blend of DFT, Judd-Ofelt theory, IET theory, and rate equation modeling. Thanks to this study, for each isomeric species, the most efficient IET process feeding the Eu(III) excited state, its related physical mechanism (exchange interaction), and the reasons for a better or worse overall energy transfer efficiency (ηsens) in the different complexes were determined. The spectroscopically measured ϕovl values are in good agreement with the ones obtained theoretically in this work.

Allosteric control of ACE2 peptidase domain dynamics.

The Angiotensin Converting Enzyme 2 (ACE2) assists the regulation of blood pressure and is the main target of the coronaviruses responsible for SARS and COVID19. The catalytic function of ACE2 relies on the opening and closing motion of its peptidase domain (PD). In this study, we investigated the possibility of allosterically controlling the ACE2 PD functional dynamics. After confirming that ACE2 PD binding site opening-closing motion is dominant in characterizing its conformational landscape, we observed that few mutations in the viral receptor binding domain fragments were able to impart different effects on the binding site opening of ACE2 PD. This showed that binding to the solvent exposed area of ACE2 PD can effectively alter the conformational profile of the protein, and thus likely its catalytic function. Using a targeted machine learning model and relative entropy-based statistical analysis, we proposed the mechanism for the allosteric perturbation that regulates the ACE2 PD binding site dynamics at atomistic level. The key residues and the source of the allosteric regulation of ACE PD dynamics are also presented.

Structural study of 1- and 2-naphthol: new insights into the non-covalent H-H interaction in cis-1-naphthol.

Previous microwave studies of naphthol monomers were supplemented by measuring spectra of all 13C mono-substituted isotopologues of the cis- and trans-conformers of 1-naphthol and 2-naphthol in their natural abundances. The resulting data were utilized to determine substitution structures and so-called semi-experimental effective structures. Results from electronic structure calculations show that the OH group of cis-1-naphthol points ≈6° out of plane, which is consistent with the inertial defect data of cis- and trans-1-naphthol. The non-planarity of cis-1-naphthol is a result of a close-contact H-atom-H-atom interaction. This type of H-H interaction has been the subject of much controversy in the past and we provide here an in-depth theoretical analysis of it. The naphthol system is particularly well-suited for such analysis as it provides internal standards with its four different isomers. The methods used include quantum theory of atoms in molecules, non-covalent interactions, independent gradient model, local vibrational mode, charge model 5, and natural bond orbital analyses. We demonstrate that the close-contact H-H interaction is neither a purely attractive nor repulsive interaction, but rather a mixture of the two.

Automatic Generation of Local Vibrational Mode Parameters: From Small to Large Molecules and QM/MM Systems.

LModeAGen, a new protocol for the automatic determination of a nonredundant, complete set of local vibrational modes is reported, which is based on chemical graph concepts. Whereas local mode properties can be calculated for a selection of parameters targeting specific local modes of interest, a complete set of nonredundant local mode parameters is requested for the adiabatic connection scheme (ACS), relating each local vibrational mode with a normal mode counterpart, and for the decomposition of normal modes (CNM) in terms of local mode contributions, a unique way to analyze vibrational spectra. So far, nonredundant parameter sets have been generated manually following chemical intuition or from a set of redundant parameters in a trial-and-error fashion, which has hampered the study of larger systems with hundreds of parameters. LModeAGen was successfully applied for a test set of 11 systems, ranging from small molecules to the large QM (>100 atoms) subsystem of carbomonoxy-neuroglobin protein, described with a hybrid QM/MM method. The ωB97X-D/aug-cc-pVDZ, M06L/def2-TZVP, and QM/MM ωB97X-D/6-31G(d,p)/AMBER model chemistries were adopted for the description of the molecules in the test set. Our new protocol is an important step forward for a routine ACS and CNM analysis of the vibrational spectra of complex and large systems with hundreds of atoms, providing new access to important encoded electronic structure information.

The Local Vibrational Mode Theory and Its Place in the Vibrational Spectroscopy Arena.

This Feature Article starts highlighting some recent experimental and theoretical advances in the field of IR and Raman spectroscopy, giving a taste of the breadth and dynamics of this striving field. The local mode theory is then reviewed, showing how local vibrational modes are derived from fundamental normal modes. New features are introduced that add to current theoretical efforts: (i) a unique measure of bond strength based on local mode force constants ranging from bonding in single molecules in different environments to bonding in periodic systems and crystals and (ii) a new way to interpret vibrational spectra by pinpointing and probing interactions between particular bond stretching contributions to the normal modes. All of this represents a means to work around the very nature of normal modes, namely that the vibrational motions in polyatomic molecules are delocalized. Three current focus points of the local mode analysis are reported, demonstrating how the local mode analysis extracts important information hidden in vibrational spectroscopy data supporting current experiments: (i) metal-ligand bonding in heme proteins, such as myoglobin and neuroglobin; (ii) disentanglement of DNA normal modes; and (iii) hydrogen bonding in water clusters and ice. Finally, the use of the local mode analysis by other research groups is summarized. Our vision is that in the future local mode analysis will be routinely applied by the community and that this Feature Article serves as an incubator for future collaborations between experiment and theory.

Bonding in nitrile photo-dissociating ruthenium drug candidates-A local vibrational mode study.

In this work, we investigated bonding features of 15 ruthenium(II) nitrile complexes of the type [Ru(tpy)(L)-(CH3CN)]n+, containing the tridentate tpy ligand (tpy = 2,2':6',2″-terpyridine) and various bidentate ancillary ligands L; 12 compounds originally synthesized by Loftus et al. [J. Phys. Chem. C 123, 10291-10299 (2019)] and three new complexes. We utilized local vibrational force constants derived from the local mode theory as a quantitative measure of bond strength complemented with the topological analysis of the electron density and the natural bond orbital analysis. Loftus et al. suggested that nitrile dissociation occurs after light induced singlet-triplet transition of the original complexes and they used as a measure of nitrile release efficiency quantum yields for ligand exchange in water. They observed larger quantum yields for complexes with smaller singlet-triplet energy gaps. The major goal of this work was to assess how the Ru-NC and Ru-L bond strengths in these 15 compounds relate to and explain the experimental data of Loftus et al., particularly focusing on the question whether there is a direct correlation between Ru-NC bond strength and measured quantum yield. Our study provides the interesting result that the compounds with the highest quantum yields also have the strongest Ru-NC bonds suggesting that breaking the Ru-NC bond is not the driving force for the delivery process rather than the change of the metal framework as revealed by first results of a unified reaction valley approach investigation of the mechanism. Compounds with the highest quantum yield show larger electronic structure changes upon singlet-triplet excitation, i.e., larger changes in bond strength, covalency, and difference between the singlet and triplet HOMOs, with exception of the compound 12. In summary, this work provides new insights into the interplay of local properties and experimental quantum yields forming in synergy a useful tool for fine tuning of existing and future design of new nitrile releasing ruthenium compounds. We hope that this work will bring theoretical and experimental studies closer together and serves as an incubator for future collaboration between computational chemists and their experimental colleagues.

Uranium: The Nuclear Fuel Cycle and Beyond.

This review summarizes the recent developments regarding the use of uranium as nuclear fuel, including recycling and health aspects, elucidated from a chemical point of view, i.e., emphasizing the rich uranium coordination chemistry, which has also raised interest in using uranium compounds in synthesis and catalysis. A number of novel uranium coordination features are addressed, such the emerging number of U(II) complexes and uranium nitride complexes as a promising class of materials for more efficient and safer nuclear fuels. The current discussion about uranium triple bonds is addressed by quantum chemical investigations using local vibrational mode force constants as quantitative bond strength descriptors based on vibrational spectroscopy. The local mode analysis of selected uranium nitrides, N≡U≡N, U≡N, N≡U=NH and N≡U=O, could confirm and quantify, for the first time, that these molecules exhibit a UN triple bond as hypothesized in the literature. We hope that this review will inspire the community interested in uranium chemistry and will serve as an incubator for fruitful collaborations between theory and experimentation in exploring the wealth of uranium chemistry.

Papaya Leaf Extracts as Potential Dengue Treatment: An In-Silico Study.

Dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) cause serious public health problems, with nearly 390 million people affected and 20,000 deaths per year in tropical and subtropical countries. Despite numerous attempts, no antiviral drug or vaccine is currently available to combat the manifestation. The challenge of discovering an efficient vaccine is enhanced by the surplus presence of efficient vectors and drug resistance from the virus. For centuries, papaya (Carica papaya) extracts have been traditionally used to treat DF, DHF, and DSS. In the present study, we systematically investigated seven compounds isolated from papaya leaf extract with regard to their potential as inhibitors for non-structural (NS) proteins, NS3 and NS5, which play a crucial role in viral RNA replication. The computational tools applied stretched across classical molecular docking, molecular dynamics (MD) simulations and SwissADME used to calculate binding affinities; binding free energies; Absorption, Distribution, Metabolism, and Excretion (ADME); and drug-likeness properties, thus, identifying Kaempferol, Chlorogenic acid, and Quercetin as potential candidates, with Kaempferol and Quercetin scoring best. Therefore, for the Kaempferol and Quercetin complexes, hybrid quantum mechanical/molecular mechanical (QM/MM) geometry and frequency calculations were performed, followed by the local mode analysis developed in our group to quantify Kaempferol-NS and Quercetin-NS hydrogen bonding. Given the non-toxic nature and the wide availability of the Kaempferol and Quercetin papaya extract in almost all of the susceptible regions, and our results showing high NS3 and NS5 binding affinities and energies, strong hydrogen bonding with both NS3 and NS5, and excellent ADME properties, we suggest Kaempferol and Quercetin as a strong NS3 and NS5 inhibitor to be further investigated in vitro.

Dihydrogen Bonding-Seen through the Eyes of Vibrational Spectroscopy.

In this work, we analyzed five groups of different dihydrogen bonding interactions and hydrogen clusters with an H3+ kernel utilizing the local vibrational mode theory, developed by our group, complemented with the Quantum Theory of Atoms-in-Molecules analysis to assess the strength and nature of the dihydrogen bonds in these systems. We could show that the intrinsic strength of the dihydrogen bonds investigated is primarily related to the protonic bond as opposed to the hydridic bond; thus, this should be the region of focus when designing dihydrogen bonded complexes with a particular strength. We could also show that the popular discussion of the blue/red shifts of dihydrogen bonding based on the normal mode frequencies is hampered from mode-mode coupling and that a blue/red shift discussion based on local mode frequencies is more meaningful. Based on the bond analysis of the H3+(H2)n systems, we conclude that the bond strength in these crystal-like structures makes them interesting for potential hydrogen storage applications.