Interfacial Water in the SARS Spike Protein: Investigating the Interaction with Human ACE2 Receptor and In Vitro Uptake in A549 Cells.

The severity of global pandemic due to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has engaged the researchers and clinicians to find the key features triggering the viral infection to lung cells. By utilizing such crucial information, researchers and scientists try to combat the spread of the virus. Here, in this work, we performed in silico analysis of the protein-protein interactions between the receptor-binding domain (RBD) of the viral spike protein and the human angiotensin-converting enzyme 2 (hACE2) receptor to highlight the key alteration that happened from SARS-CoV to SARS-CoV-2. We analyzed and compared the molecular differences between spike proteins of the two viruses using various computational approaches such as binding affinity calculations, computational alanine, and molecular dynamics simulations. The binding affinity calculations showed that SARS-CoV-2 binds a little more firmly to the hACE2 receptor than SARS-CoV. The major finding obtained from molecular dynamics simulations was that the RBD-ACE2 interface is populated with water molecules and interacts strongly with both RBD and ACE2 interfacial residues during the simulation periods. The water-mediated hydrogen bond by the bridge water molecules is crucial for stabilizing the RBD and ACE2 domains. Near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS) confirmed the presence of vapor and molecular water phases in the protein-protein interfacial domain, further validating the computationally predicted interfacial water molecules. In addition, we examined the role of interfacial water molecules in virus uptake by lung cell A549 by binding and maintaining the RBD/hACE2 complex at varying temperatures using nanourchins coated with spike proteins as pseudoviruses and fluorescence-activated cell sorting (FACS) as a quantitative approach. The structural and dynamical features presented here may serve as a guide for developing new drug molecules, vaccines, or antibodies to combat the COVID-19 pandemic.

Orientational order and dynamics of interfacial water near a hexagonal boron-nitride sheet: An ab initio molecular dynamics study.

Structural and dynamical properties of interfacial water molecules near a hexagonal boron nitride sheet (h-BN) are investigated by means of Born-Oppenheimer molecular dynamics simulations. Orientational profiles in the interfacial regions reveal two distinct types of water molecules near the BN surface. Depending on the positions of the water molecules, on top of either N or B atoms, one type contains water molecules that are oriented with one OH bond pointing toward the N atoms and the other type contains water molecules that remain parallel to the BN sheet. Distinct hydrogen bonding and stabilization energies of these two types of water molecules are found from our calculations. In order to see the effects of dispersion interactions, simulations are performed with the BLYP (Becke-Lee-Yang-Parr) functional and also BLYP with Grimme's D3 corrections (BLYP-D3). An enhancement of water ordering near the surface is observed with the inclusion of dispersion corrections. Further analysis of the diffusion coefficients, rotational time correlation functions, and hydrogen bond dynamics shows that water molecules near the h-BN sheet move faster compared to bulk water molecules both translationally and rotationally. The water molecules in the first layer are found to show substantial lateral diffusion. The escape dynamics of water from the solvation layer at the BN surface is also looked at in the current study. We have also investigated some of the electronic properties of interfacial water such as the charge density and dipole moment. It is found that the water molecules at the surface of the BN sheet have a lower dipole moment than bulk molecules.

Wetting and dewetting of narrow hydrophobic channels by orthogonal electric fields: Structure, free energy, and dynamics for different water models.

Wetting and dewetting of a (6,6) carbon nanotube in presence of an orthogonal electric field of varying strengths are studied by means of molecular dynamics simulations using seven different models of water. We have looked at filling of the channel, occupancy and structure of water inside it, associated free energy profiles, and also dynamical properties like the time scales of collective dipole flipping and residence dynamics. For the current systems where the entire simulation box is under the electric field, the nanotube is found to undergo electrodrying, i.e., transition from filled to empty states on increase of the electric field. The free energy calculations show that the empty state is the most stable one at higher electric field as it raptures the hydrogen bond environment inside the carbon nanotube by reorienting water molecules to its direction leading to a depletion of water molecules inside the channel. We investigated the collective flipping of water dipoles inside the channel and found that it follows a fast stepwise mechanism. On the dynamical side, the dipole flipping is found to occur at a faster rate with increase of the electric field. Also, the rate of water flow is found to decrease dramatically as the field strength is increased. The residence time of water molecules inside the channel is also found to decrease with increasing electric field. Although the effects of electric field on different water models are found to be qualitatively similar, the quantitative details can be different for different models. In particular, the dynamics of water molecules inside the channel can vary significantly for different water models. However, the general behavior of wetting and dewetting transitions, enhanced dipole flips, and shorter residence times on application of an orthogonal electric field hold true for all water models considered in the current work.

Dynamical Crossover of Interfacial Water upon Melting of a Lipid Bilayer: Influence of Different Parts of the Headgroups.

All-atom molecular dynamics simulations of a 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer in contact with liquid water were performed at different temperatures ranging from 285 to 320 K. We have investigated the heterogeneity and dynamical transitions in interfacial water as the lipid bilayer undergoes a melting transition. Results are obtained for water at the outer surface of the bilayer and for those buried more deeply in the lipid chains of the bilayer. It is found that lipid bilayer melting influences both the structure and dynamics of interfacial water. The number of interfacial water molecules shows a jump in the melting of the bilayer. The temperature dependence of the diffusivity and orientational relaxation of interfacial water molecules exhibits a dynamical crossover upon melting of the bilayer. The extent of dynamical crossover is found to be rather strong with significant changes in activation barriers for interfacial water around the carbonyl groups, which are deeply buried toward the lipid chains of the bilayer. The dynamical crossover gradually decreases as one moves further away from the outer surface, and it essentially vanishes for water in the region of 5-10 Å from the outer surface. It is found that the lipid melting-induced dynamical crossover of interfacial water is significant only for water that is in close proximity to the bilayer surface or deeply buried into it. The current results reveal that water molecules in different parts of the interface respond differently on melting of the bilayer. The current study also shows that the carbonyl-bound water molecules can play an important role in the phase transition of the bilayer as the temperature is raised through its melting point.

In silico screening of phytochemical compounds and FDA drugs as potential inhibitors for NSP16/10 5' methyl transferase activity.

The recent global pandemic associated with the highly contagious novel coronavirus (SARS-CoV-2) has led to an unpredictable loss of life and economy worldwide, and the discovery of antiviral drugs is an urgent necessity. For the discovery of new drug leads and for the treatment of various diseases, natural products and purified photochemical from medicinal plants are used. The RNA cap was methylated by two S-adenosyl-L-methionine (SAM)-dependent methyltransferases of SARS coronavirus (SARS-CoV-2), catalyzed by NSP16 2'-O-Mtase. Natural substrate SAM, 128 Phytocompounds retrieved from the Phytocompounds database, and 11 standard FDA-approved HIV drugs reclaimed from the PubChem database are subjected to docking analysis. The docking study was done using AutoDock Vina. Further, admetSAR and DruLiTO servers are used to analyze the drug-likeness properties. The NSP16/10 structure and natural substrate SAM, Phytocompounds Withanolide (WTL), and HIV standard drug Dolutegravir (DLT) as hit compounds were identified by molecular dynamics using the Gromacs GPU-enabled package. To examine the effectiveness of the identified drugs versus COVID-19, further in vitro and in vivo studies are required. Communicated by Ramaswamy H. Sarma.

In silico studies disclose the underlying link between binding affinity and redox potential in laccase isoforms.

Laccases are copper-containing enzymes belonging to the family of multicopper oxidases (MCOs). All MCOs use molecular oxygen to oxidize a wide range of organic compounds by radical catalysis. One of the key fundamental properties of laccases is having high or low redox potentials depending on the active site organization. Several experimental studies have been done to rationalize the high and low redox potential laccases (LRPL), however, molecular understanding is still lacking. In this work, we explored the proteomic profile of laccases produced in the fungal cultures, specifically induced with lignocellulosic biomass such as rice straw. This study was undertaken to explain the differences in the high redox and low redox potential values of different laccases using in-silico approaches. Proteomic profiling and structural and sequence analysis revealed a low level of similarity among them. Docking analyses and molecular dynamics simulation analysis revealed that high redox potential laccases (HRPL) are having good binding affinity compared to low or medium redox potential laccases (MRPL). The stability of these complexes was further analyzed based on reactive distances, active site volume comparison and a number of tunnel formations that were observed to be significantly higher for HRPL. Our results indicate that the number of tunnel formations calculated from the simulation's trajectories and available water molecules at the T3 site directly correlates with the laccases' redox potentials. This study will be helpful and provide valuable inputs for the designing of new laccases to improve lignin degradation.Communicated by Ramaswamy H. Sarma.