Among the variants of SARS-CoV-2, some are more infectious than the Wild-type. Interestingly, these mutations enable the virus to evade the therapeutic efforts. Hence, there is a need for candidate drug molecules that can potently bind with all the variants. We have adopted a strategy combining virtual screening, molecular docking followed by rigorous sampling by metadynamics simulations to find candidate molecules. From our results we found four highly potent drug candidates that can bind to the Spike-RBD of all the variants of the virus. Additionally, we also found that certain signature residues on the RBM region commonly bind to each of these inhibitors. Thus, our study not only gives information on the chemical compounds, but also residues on the proteins which could be targeted for future drug and vaccine development studies.

SARS-CoV-2-encoded accessory protein ORF3a was found to be a conserved coronavirus protein that shows crucial roles in apoptosis in cells as well as in virus release and replications. To complete the knowledge and identify the unknown of this protein, further comprehensive research is needed to clarify the leading role of ORF3a in the functioning of the coronavirus. One of the efficient approaches to determining the functionality of this protein is to investigate the mechanical properties and study its structural dynamics in the presence of physical stimuli. Herein, performing all-atom steered molecular dynamics (SMD) simulations, the mechanical properties of the force-bearing components of the ORF3a channel are calculated in different physiological conditions. As variations occurring in ORF3a may lead to alteration in protein structure and function, the G49V mutation was also simulated to clarify the relationship between the mechanical properties and chemical stability of the protein by comparing the behavior of the wild-type and mutant Orf3a. From a physiological conditions point of view, it was observed that in the solvated system, the presence of water molecules reduces Young's modulus of TM1 by ∼30 %. Our results also show that by substitution of Gly49 with valine, Young's modulus of the whole helix increases from 1.61 ± 0.20 to 2.08 ± 0.15 GPa, which is consistent with the calculated difference in free energy of wild-type and mutant helices. In addition to finding a way to fight against Covid-19 disease, understanding the mechanical behavior of these biological nanochannels can lead to the development of the potential applications of the ORF3a protein channel, such as tunable nanovalves in smart drug delivery systems, nanofilters in the new generation of desalination systems, and promising applications in DNA sequencing.

Inhibiting the biological activity of SARS-CoV-2 Mpro can prevent viral replication. In this context, a hybrid approach using knowledge- and physics-based methods was proposed to characterize potential inhibitors for SARS-CoV-2 Mpro. Initially, supervised machine learning (ML) models were trained to predict a ligand-binding affinity of ca. 2 million compounds with the correlation on a test set of . Atomistic simulations were then used to refine the outcome of the ML model. Using LIE/FEP calculations, nine compounds from the top 100 ML inhibitors were suggested to bind well to the protease with the domination of van der Waals interactions. Furthermore, the binding affinity of these compounds is also higher than that of nirmatrelvir, which was recently approved by the US FDA to treat COVID-19. In addition, the ligands altered the catalytic triad Cys145 - His41 - Asp187, possibly disturbing the biological activity of SARS-CoV-2.

Thiazole orange (TO), an asymmetric cyanine dye, has been widely used in biomolecular detection and imaging of DNA/ RNA in gels, due to its unique fluorogenic behavior: fluorescence of free dye in aqueous solution is very weak but emission can be significantly enhanced in nucleic-acid-bound dye. Herein we describe the ultrafast excited-state dynamics of free TO in aqueous solution by exploiting both a femtosecond upconversion spectrophotofluorometer and a picosecond time-correlated single-photon counting (TCSPC) apparatus. For the first time, the fluorescence lifetime of TO monomer in water was found to be ∼1 ps, mixed with concurrent solvent relaxation (which was confirmed by the experimental results of TO in DMSO). Even at moderate concentration, this lifetime has an amplitude (a measure of molecular fraction) that significantly dominates other lifetimes, and this is the origin of weak steady state fluorescence of free TO in water. We also found a novel slower decay component around 34 ps. Interestingly and in addition, the lifetime component on the 30-40 ps timescale was also found in TO--Cyclodextrin (CD) complexes. The fraction of this component increased with the addition of -CD. Cyclodextrin has been reported to promote the aggregation of TO. Thus, although a very coincidental match of this time constant by one for a torsional process within the cavity can not be ruled out, we ascribe the shared 30-40 ps component to the lifetime of a highly quenched TO dimer experiencing intra-and inter-molecular rearrangement.

The SARS-CoV-2 coronavirus (COVID-19) that is causing the massive global pandemic exhibits similar human cell invasion mechanism as the coronavirus SARS-CoV, which had significantly lower fatalities. The cell membrane protein Angiotensin-converting enzyme 2 (ACE2) is the initiation point for both the coronavirus infections in humans. Here, we model the molecular interactions and mechanical properties of ACE2 with both SARS-CoV and COVID-19 spike protein receptor-binding domains (RBD). We report that the COVID-19 spike RBD interacts with ACE2 more strongly and at only two protein residues, as compared to multi-residue interaction of the SARS-CoV. Although both coronaviruses stiffen the ACE2, the impact of COVID-19 is six times larger, which points towards differences in the severity of the reported respiratory distress. The recognition of specific residues of ACE2 attachments to coronaviruses is important as the residues suggest potential sites of intervention to inhibit attachment and subsequent entry of the COVID-19 into human host cells.

In the present investigation, molecular docking studies have been performed using AutoDock Vina to investigate the role of ligand-binding affinity at the hydrophobic pocket of COVID-19. The knowledge of the binding of protein receptors with ligand molecules is essential in drug discovery processes. Hydroxamic acids with reported biological activity, have been investigated for docking to an important target, SARS-CoV-2, in order to predict their therapeutic efficacy. The spike protein of the coronavirus is responsible for the attachment to host cells and a positive-sense single-strand RNA, (+)ssRNA, is a genetic material that can be translated into protein in the host cell. We modeled the structure of SARS-CoV-2 with the ligands, hydroxamic acids. They show binding capability with both, Spike protein and (+)ssRNA. The twain exhibit negative binding energies which signify that reactions are spontaneous, strong, and fast. The present research proposed hydroxamic acids as molecules which can be used for the development of anti-virals therapeutics against SARS-CoV-2.

We present the first exact solution to the time-independent Schrödinger equation of a model Hamiltonian consisting of a vibrational mode coupled to three electronic states. This Hamiltonian serves as a generic model for photo-induced electronic transfer reactions. The solution is non-perturbative and can be applied to ET reactions with weak and strong electronic and vibrational coupling strengths. This work suggests a new direction towards understanding the vibronic effects in ET dynamics beyond the non-adiabatic limit and Condon approximation.

Coronavirus disease 2019 (COVID-19) has caused more than 840,000 deaths as of 31 August 2020 in the whole world. The COVID-19 main protease (M) has been validated as an attractive target for drug design. In this work, the binding mechanisms of five protease inhibitors (e.g., danoprevir, darunavir, ASC09, lopinavir and ritonavir) to COVID-19 M were investigated. Based on the docking score, five protease inhibitors structures were selected for further evaluation. It is found that most of the selected drug molecules bind stably to the COVID-19 M from the molecular dynamics simulation. Moreover, the MM/PBSA free energy calculations suggest that lopinavir with positive charge might be most active against COVID-19 M.

We calculate the nearest-neighbour enthalpies and entropies at 5 salt concentrations of 18 RNA sequences, each for at least 9 different species concentrations, totalling 757 melting temperatures, using a melting temperature optimization method. These new parameters do not need to be salt-corrected and are shown to provide overall improved melting temperature predictions. They show a marked quadratic dependence with salt concentrations which are compensated to form linear Gibbs free energies. Two different parameter schemes were tested, with fixed or variable initial parameters. We have found that using variable initial parameters provides better predictive results than using salt correction factors and that the prediction uncertainty is considerably reduced for a validation set of independent sequences. An interpolation scheme is introduced to generate model parameters for arbitrary salt concentrations which performs better against a validation set than predictions using salt corrections.

Graphical representation methods constitute a class of alignment-free techniques for comparative study of biomolecular sequences. In this brief commentary, we study how well some of these methods can discriminate among closely related genes.

One of the most important electrostatic interactions between molecules is most definitely the hydrogen bond. Understanding the basis of this interaction may offer us the insight needed to understand its effect on the macroscopic scale. Hydrogen bonding is for example the reason for anomalous properties in compounds like water and naturally life as we know it. The strength of the bond depends on numerous factors, among them the electronegativity of participating atoms. In this work we calculated the strength of hydrogen bonds between hydrides of the upper-right part of the periodic table (C, N, O, F, P, S, Cl, As, Se, Br) using quantum-chemical methods. The aim was to determine what influences the strength of strong and weak hydrogen bonds in simple hydrides. Various relationships were checked. A relation between the strength of the bond and the electronegativity of the participating atoms was found. We also observed a correlation between the strength of hydrogen bonds and the inter-atomic distances, along with the dependence on the charge transfer on the atom of the donor. We also report characteristic geometries of different dimers.

Durability of organo-lead halide perovskite are important issue for its practical application in a solar cells. In this study, using density functional theory (DFT) and molecular dynamics, we theoretically investigated a crystal structure, electronic structure, and ionic diffusivity of the partially substituted cubic MAXPbI (MA = CHNH, X = NH or (NH)CH or Cs). Our calculation results indicate that a partial substitution of MA induces a lattice distortion, resulting in preventing MA or X from the diffusion between A sites in the perovskite. DFT calculations show that electronic structures of the investigated partially substituted perovskites were similar with that of MAPbI, while their bandgaps slightly decrease compared to that of MAPbI Our results mean that partial substitution in halide perovskite is effective technique to suppress diffusion of intrinsic ions and tune the band gap.

Molecular dynamics simulations and NMR spectroscopy were used to compare the binding of two β-blocker drugs to the chiral molecular micelle poly-(sodium undecyl-(L)-leucine-valine). The molecular micelle is used as a chiral selector in capillary electrophoresis. This study is part of a larger effort to understand the mechanism of chiral recognition in capillary electrophoresis by characterizing the molecular micelle binding of chiral compounds with different geometries and charges. Propranolol and atenolol were chosen because their structures are similar, but their chiral interactions with the molecular micelle are different. Molecular dynamics simulations showed both propranolol enantiomers inserted their aromatic rings into the molecular micelle core and that (S)-propranolol associated more strongly with the molecular micelle than (R)-propranolol. This difference was attributed to stronger molecular micelle hydrogen bonding interactions experienced by (S)-propranolol. Atenolol enantiomers were found to bind near the molecular micelle surface and to have similar molecular micelle binding free energies.

Infrared and Raman spectra of 2,3,5,6-tetrafluoropyridine (TFPy) were recorded and vibrational frequencies were assigned for its S electronic ground states. and density functional theory (DFT) calculations were used to complement the experimental work. The lowest electronic excited state of this molecule was investigated with ultraviolet absorption spectroscopy and theoretical CASSCF calculations. The band origin was found to be at 35,704.6 cm in the ultraviolet absorption spectrum. A slightly puckered structure with a barrier to planarity of 30 cm was predicted by CASSCF calculations for the S(π,π*) state. Lower frequencies for the out-of-plane ring bending vibrations for the electronic excited state result from the weaker bonding within the pyridine ring.

Enzymes are the most efficient chemical catalysts known, but the exact nature of chemical barrier crossing in enzymes is not fully understood. Application of transition state theory to enzymatic reactions indicates that the rates of all possible reaction paths, weighted by their relative probabilities, must be considered in order to achieve an accurate calculation of the overall rate. Previous studies in our group have shown a single mechanism for enzymatic barrier passage in human heart lactate dehydrogenase (LDH). To ensure that this result was not due to our methodology insufficiently sampling reactive phase space, we implement high-perturbation transition path sampling in both microcanonical and canonical regimes for the reaction catalyzed by human heart LDH. We find that, although multiple, distinct paths through reactive phase space are possible for this enzymatic reaction, one specific reaction path is dominant. Since the frequency of these paths in a canonical ensemble is inversely proportional to the free energy barriers separating them from other regions of phase space, we conclude that the rarer reaction paths are likely to have a negligible contribution. Furthermore, the non-dominate reaction paths correspond to altered reactive conformations and only occur after multiple steps of high perturbation, suggesting that these paths may be the result of non-biologically significant changes to the structure of the enzymatic active site.

Here we report sub-100-ps structural dynamics of horse heart myoglobin revealed by time-resolved X-ray solution scattering. By applying the time-slicing scheme to the measurement and subsequent deconvolution, we investigate the protein structural dynamics that occur faster than the X-ray temporal pulse width of synchrotrons (~100 ps). The singular value decomposition analysis of the experimental data suggests that two structurally distinguishable intermediates are formed within 100 ps. In particular, the global structural change occurring on the time scale of 70 ps is identified.

Molecular dynamics (MD) simulations were used to investigate the binding of 1,1'-binaphthyl-2,2'-diyl hydrogenphosphate (BNP) enantiomers to the molecular micelle poly-(sodium undecyl-(L,L)-leucine-valine) (poly(SULV)). Poly(SULV) is used as a chiral selector in capillary electrophoresis separations. Four poly(SULV) binding pockets were identified and either (R)-BNP or (S)-BNP were docked into each pocket. MD simulations were then used to identify the preferred BNP binding site. Within the preferred site, both enantiomers formed hydrogen bonds with poly(SULV) and penetrated into the poly(SULV) core. Comparisons of BNP enantiomer binding to the preferred poly(SULV) pocket showed that (S)-BNP formed stronger hydrogen bonds, moved deeper into the binding site, and had a lower poly(SULV) binding free energy than the (R) enantiomer. Finally, MD simulation results were in agreement with capillary electrophoresis and NMR experiments. Each technique showed (S)-BNP interacted more strongly with poly(SULV) than (R)-BNP and that the site of chiral recognition was near the poly(SULV) leucine chiral center.

We have developed the method of picosecond Laue crystallography and used this capability to probe ligand dynamics in tetrameric R-state hemoglobin (Hb). Time-resolved, 2 Å-resolution electron density maps of photolyzed HbCO reveal the time-dependent population of CO in the binding (A) and primary docking (B) sites of both α and β subunits from 100 ps to 10 μs. The proximity of the B site in the β subunit is about 0.25 Å closer to its A binding site, and its rebinding rate (~300 μs) is six times faster, suggesting distal control of the rebinding dynamics. Geminate rebinding in the β subunit exhibits both prompt and delayed geminate phases. We developed a microscopic model to quantitatively explain the observed kinetics, with three states for the α subunit and four states for the β subunit. This model provides a consistent framework for interpreting rebinding kinetics reported in prior studies of both HbCO and HbO.

Though there have been many studies on photosensitizers coupled to model complexes of the [FeFe]-hydrogenases, few have looked at how the models react upon exposure to light. To extract photoreaction information, ultrafast time-resolved UV/visible pump, IR probe spectroscopy was performed on Fe(μ-SCH)(CO)(PMe) () dissolved in heptane and acetonitrile and the photochemical dynamics were determined. Excitation with 532 and 355 nm light produces bleaches and new absorptions that decay to half their original intensity with time constants of 300 ± 120 ps and 380 ± 210 ps in heptane and acetonitrile, respectively. These features persist to the microsecond timescale. The dynamics of are assigned to formation of an initial set of photoproducts, which were a mixture of excited-state tricarbonyl isomers. These isomers decay into another set of long-lived photoproducts in which approximately half the excited-state tricarbonyl isomers recombine with CO to form another complex mixture of tricarbonyl and tetracarbonyl isomers.

One of the most commonly used nonsteroidal anti-inflammatory active pharmaceutical ingredient called Meloxicam has been characterized spectroscopically both by Terahertz (THz) time domain spectroscopy (THz-TDS) and by Fourier Transform Infrared (FTIR) spectroscopy in far-IR regions of electromagnetic spectrum; 0.2 THz to 20 THz. While many relatively sharp features are observed in the far-IR range between 2 THz to 20 THz as expected for being an organic substance, very distinct and relatively strong absorption bands are also observed at 1.00, 1.66, 2.07 and 2.57 THz in the THz range. These well separated, defined, and fairly strong spectral features can be used for discrimination and quantification of Meloxicam in drug analysis. Frequency dependent refractive index of the drug was determined in a range of 0.2 THz and 2.7 THz, where an almost constant index was observed with an average index of 1.75. Powder XRD, and solid-state Density Functional Theory (SS-DFT) calculations were utilized to determine the crystalline form of the Meloxicam sample in its enolic crystalline form. Single molecule DFT calculations were also performed in all four possible structures of Meloxicam. In addition, the capability of THz waves transmission through common packaging materials is demonstrated for possibility of future on-site analysis. The results suggest that drug analysis will be possible to perform not only at every stage of manufacturing without destruction but also directly at the shelf of a market after development of portable THz technologies.